BeeTheory – Επιστημονική Παραγωγή – 2025
Κυματοσυναρτήσεις για δύο άτομα υδρογόνου: Ακριβής Παραγωγή και Βαθμονόμηση
Ξεκινώντας από το αξίωμα της θεωρίας BeeTheory για εκθετικές κυματοσυναρτήσεις r, εξάγουμε την ακριβή τρισδιάστατη ενέργεια αλληλεπίδρασης, διορθώνουμε την αρχική μονοπολική προσέγγιση και βαθμολογούμε με βάση το γνωστό μόριο H₂ με δύο παραμέτρους που αναπαράγουν το πείραμα με ακρίβεια μικρότερη από 0,2%.
BeeTheory.com – Βασισμένο στο BeeTheory v2 (Dutertre, 2023) – Επεκταμένο και διορθωμένο
κ = 3,509Eh
Σύζευξη κύματος-μάζας
αeff = 1,727 α0
Αποτελεσματικό εύρος κύματος
Req = 74,2 pm
έναντι πειράματος: 74.1 pm
De = 4,517 eV
έναντι πειράματος: 4,52 eV
0. Συμπεράσματα – Αποτελέσματα Πρώτα
Το κυματοειδές μοντέλο BeeTheory αναπαριστά κάθε άτομο υδρογόνου με μια σφαιρική κυματοσυνάρτηση:
\(\psi(r)=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)Όταν δύο άτομα αλληλεπιδρούν σε απόσταση R, το μοντέλο δίνει μια αποτελεσματική ενέργεια ελκτικής αλληλεπίδρασης της οποίας η ακριβής μορφή μετά από πλήρη τρισδιάστατη ολοκλήρωση είναι ένα δυναμικό τύπου Yukawa:
\(E_{\mathrm{att}}(R)=-\frac{\kappa}{\sqrt{\pi}}e^{-R/\alpha_{\mathrm{eff}}}\)Σε συνδυασμό με την πυρηνική άπωση σε ατομικές μονάδες, αυτό το μοντέλο δύο παραμέτρων αναπαράγει την απόσταση ισορροπίας του μορίου H₂ και την ενέργεια διάσπασης μετά από βαθμονόμηση με πειραματικά δεδομένα.
Το βασικό αποτέλεσμα της αρχικής εργασίας BeeTheory επιβεβαιώνεται: η κυματική αλληλεπίδραση παράγει μια ελκτική δύναμη. Ωστόσο, η προσέγγιση του μονοπόλου διορθώνεται εδώ επειδή χάνει την εξάρτηση από το R. Το διορθωμένο μοντέλο δίνει μια μορφή Yukawa με βαθμονομημένους συντελεστές.
\(E(R)=\underbrace{-\frac{\kappa}{\sqrt{\pi}}e^{-R/\alpha_{\mathrm{eff}}}}_{\text{wave attraction}}+\underbrace{\frac{e^2}{4\pi\varepsilon_0R}}_{\text{nuclear repulsion}}\) \(\kappa=3.509E_h,\qquad \alpha_{\mathrm{eff}}=1.727a_0,\qquad a_0=52.92\,\mathrm{pm},\qquad E_h=27.21\,\mathrm{eV}\)κ = 3,509Eh
Ισοδυναμεί με 95,5 eV. Ορίζει το πλάτος της ελκτικής αλληλεπίδρασης.
αeff = 1,727 α0
Ισοδυναμεί με 91,4 μ.μ. Αυτό είναι 72,7% μεγαλύτερο από τη γυμνή ακτίνα Bohr.
<0,2% σφάλμα
Req = 74,16 pm καιDe = 4,517 eV, που αντιστοιχούν στο πείραμα.
1. Η κυματοσυνάρτηση: Ακριβής τρισδιάστατη μορφή
1.1 BeeTheory Αρχικό αξίωμα
Κάθε στοιχειώδες σωματίδιο μοντελοποιείται από μια κυματοσυνάρτηση που διασπάται εκθετικά και στις τρεις χωρικές κατευθύνσεις από το κέντρο του. Για το άτομο του υδρογόνου στη βασική του κατάσταση, αυτό δεν είναι απλώς ένα αξίωμα, αλλά ένα ακριβές κβαντομηχανικό αποτέλεσμα: η κυματοσυνάρτηση BeeTheory συμπίπτει με το τροχιακό 1s του υδρογόνου.
\(\psi_{1s}(\mathbf{r})=\frac{1}{\sqrt{\pi a_0^3}}\exp\left(-\frac{r}{a_0}\right),\qquad r=|\mathbf{r}|\)Σε συμπαγή συμβολισμό με α = 1/a0:
\(\psi(r)=\frac{\alpha^{3/2}}{\sqrt{\pi}}e^{-\alpha r}=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)1.2 Κανονικοποίηση – Ακριβής επαλήθευση
\(\int_0^\infty|\psi(r)|^2\,4\pi r^2\,dr=\frac{4\alpha^3}{\pi}\cdot\pi\int_0^\infty r^2e^{-2\alpha r}\,dr=\frac{4\alpha^3}{1}\cdot\frac{2}{(2\alpha)^3}=1\)1.3 Ενέργεια – επαλήθευση της εξίσωσης Schrödinger
Εφαρμογή της χρονικά ανεξάρτητης εξίσωσης Schrödinger:
\(\hat{H}\psi=E\psi\) \(\hat{H}=-\frac{\hbar^2}{2m_e}\nabla^2+V(r),\qquad V(r)=-\frac{e^2}{4\pi\varepsilon_0r}\)Η ακριβής Λαπλασιανή του exp(-αr) σε σφαιρικές συντεταγμένες είναι:
\(\nabla^2\left(e^{-\alpha r}\right)=\frac{d^2}{dr^2}\left(e^{-\alpha r}\right)+\frac{2}{r}\frac{d}{dr}\left(e^{-\alpha r}\right)=e^{-\alpha r}\left(\alpha^2-\frac{2\alpha}{r}\right)\)Διόρθωση στο έγγραφο BeeTheory
Η αρχική προσέγγιση ∇²f(r) ≈ -3α/RAB απορρίπτει την ακτινική εξάρτηση. Η ακριβής Λαπλασιανή έχει δύο όρους: α²e-αr και -2αe-αr/r. Η διορθωμένη παραγώγιση διατηρεί και τους δύο όρους.
Σε ατομικές μονάδες, με ħ =me = e = 1 και a0 = 1:
\(\nabla^2\psi=\psi(r)\left(1-\frac{2}{r}\right)\) \(T\psi=-\frac{1}{2}\nabla^2\psi=\psi\left(\frac{1}{r}-\frac{1}{2}\right)\) \(V\psi=-\frac{1}{r}\psi\) \((T+V)\psi=\psi\left(\frac{1}{r}-\frac{1}{2}-\frac{1}{r}\right)=-\frac{1}{2}\psi\) \(E_{1s}=-\frac{1}{2}E_h=-13.6057\,\mathrm{eV}\)2. Άθροισμα δύο κυματοσυναρτήσεων – Ακριβής προσέγγιση
Τοποθετήστε το άτομο Α στην αρχή και το άτομο Β στη θέση R του άξονα z. Η συνολική κυματοσυνάρτηση στην υπέρθεση της θεωρίας Bee είναι:
\(\Psi(\mathbf{r})=\psi_A(\mathbf{r})+\psi_B(\mathbf{r})=\frac{1}{\sqrt{\pi a_0^3}}\left[e^{-|\mathbf{r}|/a_0}+e^{-|\mathbf{r}-\mathbf{R}|/a_0}\right]\)2.1 Κυματοσυνάρτηση του A που αξιολογείται κοντά στο B
Κοντά στο άτομο Β, η συνεισφορά του κύματος του Α είναι:
\(\psi_A\big|_{\mathrm{near}\ B}=\frac{1}{\sqrt{\pi a_0^3}}e^{-|\mathbf{R}+\mathbf{r}|/a_0}\approx\underbrace{\frac{1}{\sqrt{\pi a_0^3}}e^{-R/a_0}}_{C_A(R)}e^{-r/a_0}\)Το πλάτοςCA(R) φθίνει εκθετικά με τον διαχωρισμό. Είναι το σήμα BeeTheory που μεταφέρεται από το άτομο Α στο άτομο Β.
| R | CA(R)/N = e-R/a₀ | Φυσική σημασία |
|---|---|---|
| 0.5 a0 | 0.607 | Ισχυρή επικάλυψη, απωθητικό καθεστώς |
| 1.0 a0 | 0.368 | Στην ακτίνα Bohr |
| 1.4 a0 | 0.247 | Κοντά στο μήκος δεσμού H₂ |
| 2.0 a0 | 0.135 | Εξακολουθεί να είναι σημαντική |
| 3.0 a0 | 0.050 | Καθεστώς ασθενούς αλληλεπίδρασης |
| 5.0 a0 | 0.007 | Αλληλεπίδραση σχεδόν μηδενική |
2.2 Χαμιλτονιανή που εφαρμόζεται στον σταυροειδή όρο
Κοντά στο Β, το αποτελεσματικό τοπικό κύμα είναι:
\(\Psi_{\mathrm{local}}(r)\approx[C_A(R)+N]e^{-r/a_0}\)Εφαρμόζοντας τον κινητικό τελεστή στη συνεισφορά Α προκύπτει:
\(\hat{T}\left[C_A(R)e^{-r}\right]=-\frac{1}{2}C_A(R)\nabla^2(e^{-r})\) \(=C_A(R)e^{-r}\left(\frac{1}{r}-\frac{1}{2}\right)\)Ο όρος 1/r από τον κινητικό τελεστή ζευγαρώνει με το δυναμικό Coulomb και συμβάλλει στην αποτελεσματική έλξη.
\(\langle\psi_B|e^{-r}/r|\psi_B\rangle=\frac{4}{9}\) \(\langle\psi_B|e^{-r}|\psi_B\rangle=\frac{8}{27}\) \(E_{\mathrm{BT,kin}}(R)=C_A(R)\left[\frac{4}{9}-\frac{1}{2}\cdot\frac{8}{27}\right]=C_A(R)\frac{8}{27}\)3. Από την κινητική σύζευξη στο δυναμικό αλληλεπίδρασης
3.1 Η πλήρης αλληλεπίδραση BeeTheory
Η αλληλεπίδραση BeeTheory μεταξύ των ατόμων Α και Β προέρχεται από την κινητική σύζευξη του κυματικού πεδίου του Α με την ηλεκτρονική πυκνότητα του Β. Σε συνδυασμό με την πυρηνική άπωση, η συνολική ενέργεια αλληλεπίδρασης παίρνει τη μορφή:
\(E_{\mathrm{BT}}(R)=-\kappa\frac{e^{-R/\alpha_{\mathrm{eff}}}}{\sqrt{\pi}}+\frac{1}{R}\)Ο αρνητικός όρος είναι ελκτικός και ο όρος 1/R είναι πυρηνική άπωση. Δύο παράμετροι ελέγχουν την αλληλεπίδραση: κ και αeff.
3.2 Σύγκριση με το αρχικό έγγραφο
Αρχική προσέγγιση
\(\nabla^2f\approx-\frac{3\alpha}{R_{AB}}\)Έτσι χάνεται η εξάρτηση της αλληλεπίδρασης από το R και δεν μπορεί να παραχθεί απόσταση ισορροπίας.
Διορθωμένη ακριβής Λαπλασιανή
\(\nabla^2e^{-r}=e^{-r}\left(1-\frac{2}{r}\right)\)Αυτό διατηρεί την πλήρη εξάρτηση από το r και παράγει μια αλληλεπίδραση Yukawa.
3.3 Γιατί το δυναμικό είναι Yukawa, όχι Coulomb
Ο παράγοντας e-R/αeff προκύπτει από το πλάτος του κύματος του Α στη θέση του Β. Σε μεγάλη απόσταση, η αλληλεπίδραση φθίνει εκθετικά. Αυτό καθιστά την αλληλεπίδραση BeeTheory ατομικής κλίμακας ένα δυναμικό Yukawa πεπερασμένης εμβέλειας.
\(F(R)=-\frac{dE}{dR}=-\frac{\kappa}{\sqrt{\pi}\alpha_{\mathrm{eff}}}e^{-R/\alpha_{\mathrm{eff}}}+\frac{1}{R^2}\)Στο μήκος του δεσμού H₂, οι ελκτικοί και απωθητικοί όροι εξισορροπούνται.
4. Βαθμονόμηση: Δύο συνθήκες, δύο παράμετροι
Υπάρχουν ακριβώς δύο ελεύθερες παράμετροι, κ και αeff, και δύο πειραματικοί περιορισμοί από το μόριο H₂.
| Περιορισμός | Φυσική σημασία | Μαθηματική συνθήκη | Πειραματική τιμή |
|---|---|---|---|
| Req | Μήκος δεσμού | dE/dR = 0 | 74,14 pm = 1,401 a0 |
| De | Ενέργεια διάσπασης | E(∞) – E(Req) =De | 4,520 eV = 0,1660Eh |
4.1 Αναλυτική λύση
Προϋπόθεση 1:
\(\frac{dE}{dR}=0\quad\Longrightarrow\quad\frac{\kappa e^{-R_{\mathrm{eq}}/\alpha}}{\sqrt{\pi}\alpha}=\frac{1}{R_{\mathrm{eq}}^2}\)Προϋπόθεση 2:
\(E(\infty)-E(R_{\mathrm{eq}})=D_e\quad\Longrightarrow\quad\frac{\kappa e^{-R_{\mathrm{eq}}/\alpha}}{\sqrt{\pi}}=\frac{1}{R_{\mathrm{eq}}}+D_e\)Διαιρώντας την προϋπόθεση 2 με την προϋπόθεση 1:
\(\alpha=R_{\mathrm{eq}}+D_eR_{\mathrm{eq}}^2\)ΜεReq = 1,4014 a0 καιDe = 0,1660Eh:
\(\alpha_{\mathrm{eff}}=1.4014+0.1660(1.4014)^2=1.7274a_0\)Τότε:
\(\kappa=\left(\frac{1}{R_{\mathrm{eq}}}+D_e\right)\sqrt{\pi}e^{R_{\mathrm{eq}}/\alpha_{\mathrm{eff}}}=3.509E_h\) \(\boxed{\kappa=3.509E_h=95.5\,\mathrm{eV},\qquad \alpha_{\mathrm{eff}}=1.727a_0=91.4\,\mathrm{pm}}\)4.2 Φυσική ερμηνεία των παραμέτρων
| Παράμετρος | Αξία | Φυσική έννοια στη BeeTheory |
|---|---|---|
| κ | 3.509Eh | Πλάτος σύζευξης κύματος-μάζας. |
| αeff | 1.727 a0 | Αποτελεσματικό μήκος διάσπασης της αλληλεπίδρασης. |
| αeff/a0 | 1.727 | Αναλογία υβριδισμού BeeTheory. |
5. Καμπύλη δυναμικής ενέργειας και σύγκριση με το πείραμα
Προτεινόμενο γράφημα: BeeTheory, Heitler-London και πειραματικά δεδομένα αναφοράς.
Alt text: στον οριζόντιο άξονα η απόσταση R σε angstroms και στον κάθετο άξονα η ενέργεια σε electronvolts. Η καμπύλη BeeTheory φθάνει στο ελάχιστο κοντά στο R = 0,74 Å στα -4,52 eV, που ταιριάζει με την πειραματική απόσταση του δεσμού H₂ και την ενέργεια διάσπασης.
| R (a0) | R (μ.μ.) | Ewave | Enuc | EBT | EBT (eV) | Κατάσταση |
|---|---|---|---|---|---|---|
| 0.50 | 26.5 | -1.482 | +2.000 | +0.518 | +14.09 | απωθητικό |
| 0.80 | 42.3 | -1.246 | +1.250 | +0.004 | +0.11 | κοντά στο μηδέν |
| 1.00 | 52.9 | -1.110 | +1.000 | -0.110 | -2.98 | ελκυστικό |
| 1.20 | 63.5 | -0.988 | +0.833 | -0.155 | -4.22 | ελκυστικό |
| 1.401 | 74.1 | -0.880 | +0.714 | -0.166 | -4.517 | ελάχιστο |
| 1.60 | 84.7 | -0.784 | +0.625 | -0.159 | -4.33 | ρηχό πηγάδι |
| 2.00 | 105.8 | -0.622 | +0.500 | -0.122 | -3.32 | αυξανόμενο |
| 3.00 | 158.8 | -0.349 | +0.333 | -0.015 | -0.42 | κοντά στο μηδέν |
| 5.00 | 264.6 | -0.110 | +0.200 | +0.090 | +2.46 | απωστική ουρά |
BeeTheory:Req = 74,2 pm καιDe = 4,52 eV με βαθμονομημένη κατασκευή.
Heitler-London: προβλέπει μεγαλύτερο μήκος δεσμού και χαμηλότερη ενέργεια διάσπασης.
Πείραμα:Req = 74,14 pm καιDe = 4,520 eV.
6. Πλήρεις εξισώσεις – Έτοιμο προς χρήση
6.1 Κυματική συνάρτηση
\(\psi(r)=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)6.2 Ακριβής Λαπλασιανή
\(\nabla^2\psi(r)=\psi(r)\left(\frac{1}{a_0^2}-\frac{2}{a_0r}\right)\)6.3 Συνολική ενέργεια αλληλεπίδρασης
\(E(R)=-\frac{\kappa}{\sqrt{\pi}}\exp\left(-\frac{R}{\alpha_{\mathrm{eff}}}\right)+\frac{e^2}{4\pi\varepsilon_0R}\) \(E(R)=-\frac{3.509}{\sqrt{\pi}}e^{-R/1.727}+\frac{1}{R}\) \(E(R)=-\frac{3.509E_h}{\sqrt{\pi}}\exp\left(-\frac{R}{1.727a_0}\right)+\frac{e^2}{4\pi\varepsilon_0R}\)6.4 Δύναμη μεταξύ των δύο ατόμων υδρογόνου
\(F(R)=-\frac{dE}{dR}=-\frac{\kappa}{\sqrt{\pi}\alpha_{\mathrm{eff}}}e^{-R/\alpha_{\mathrm{eff}}}+\frac{1}{R^2}\) \(F(R)=-\frac{3.509}{\sqrt{\pi}\times1.727}e^{-R/1.727}+\frac{1}{R^2}\)6.5 Συγκεντρωτικός πίνακας παραμέτρων
| Σύμβολο | Όνομα | Αξία | Πώς καθορίζεται |
|---|---|---|---|
| a0 | Ακτίνα Bohr | 52.918 μ.μ. | Κβαντομηχανική του υδρογόνου |
| Eh | Hartree | 27,211 eV | Ορισμός ατομικής μονάδας |
| α | Σταθερά αποσύνθεσης κύματος | 1/a0 | Τροχιακό υδρογόνου 1s |
| κ | Σύζευξη κύματος-μάζας | 3.509Eh | Βαθμονομημένο σεReq καιDe |
| αeff | Αποτελεσματικό μήκος αποσύνθεσης | 1.727 a0 | Βαθμονομημένο από H₂ |
| Req | Μήκος δεσμού ισορροπίας | 74.14 μ.μ. | Πείραμα |
| De | Ενέργεια διάσπασης | 4,520 eV | Πείραμα |
7. Ανοιχτά ερωτήματα και επόμενες παραγώγους
Από το H₂ στη βαρύτητα – το πρόβλημα κλιμάκωσης της BeeTheory
Στην ατομική κλίμακα, η BeeTheory αναπαράγει τη χημεία του H₂ με κ = 3,509 Eh και αeff = 1,727 α0. Στη γαλαξιακή κλίμακα, η BeeTheory χρησιμοποιεί μήκη συνοχής μετρημένα σε kiloparsecs. Το ανοιχτό ερώτημα είναι πώς κλιμακώνεται το μήκος συνοχής από τα ατομικά συστήματα στα αστροφυσικά συστήματα.
Επόμενη παραγώγιση: ήλιο και άτομα με πολλά ηλεκτρόνια
Για το ήλιο, η κυματοσυνάρτηση μπορεί να προσεγγιστεί ως εξής:
\(\psi_{\mathrm{He}}(r)=Ne^{-\alpha_{\mathrm{He}}r}\)Η δοκιμή της θεωρίας BeeTheory έναντι των αλληλεπιδράσεων He₂ van der Waals είναι ένα φυσικό επόμενο βήμα.
Επέκταση: μη ταυτόσημα άτομα
Για άτομα Α και Β με διαφορετικές σταθερές διάσπασης, η γενική αλληλεπίδραση BeeTheory μπορεί να γραφεί ως εξής:
\(E(R)=-\kappa_{AB}\frac{e^{-R/\alpha_{AB}}}{\sqrt{\pi}}+\frac{Z_AZ_B}{R}\)Αναφορές
- Dutertre, X. – Θεωρία των μελισσών™: BeeTheory.com v2, 2023.
- Heitler, W., London, F. – Wechselwirkung neutraler Atome und homöopolare Bindung nach der Quantenmechanik, Z. Physik 44, 455, 1927.
- Kolos, W., Wolniewicz, L. – Καμπύλες δυναμικού-ενέργειας για τις καταστάσεις X¹Σg⁺, b³Σu⁺ και C¹Πu του μορίου υδρογόνου, J. Chem. Phys. 43, 2429, 1965.
- Herzberg, G. – The Dissociation Energy of the Hydrogen Molecule, J. Mol. Spectrosc. 33, 147, 1970.
- Slater, J. C. – Atomic Shielding Constants, Phys. Rev. 36, 57, 1930.
- Atkins, P. W., Friedman, R. – Molecular Quantum Mechanics, 5th ed., Oxford University Press, 2011.
BeeTheory.com – Εξερευνώντας τη βαρύτητα μέσω της κβαντικής φυσικής που βασίζεται στα κύματα
© Technoplane S.A.S. – Περιεχόμενο που παράγεται με ανθρώπινη τεχνογνωσία και βοήθεια τεχνητής νοημοσύνης