Teoria pszczół – pochodna naukowa – 2025 r.
Funkcje falowe dla dwóch atomów wodoru: Rygorystyczne wyprowadzenie i kalibracja
Wychodząc od postulatu BeeTheory o wykładniczych funkcjach falowych r, wyprowadzamy dokładną energię interakcji 3D, korygujemy oryginalne przybliżenie monopolu i kalibrujemy względem znanej cząsteczki H₂ z dwoma parametrami, które odtwarzają eksperyment z dokładnością mniejszą niż 0,2%.
BeeTheory.com – Na podstawie BeeTheory v2 (Dutertre, 2023) – Rozszerzone i poprawione
κ = 3,509 Eh
Sprzężenie fala-masa
αeff = 1,727 a0
Efektywny zakres fal
Req = 74,2 pm
vs eksperyment: 74,1 pm
De = 4,517 eV
vs eksperyment: 4,52 eV
0. Wnioski – najpierw wyniki
Model falowy BeeTheory reprezentuje każdy atom wodoru za pomocą sferycznej funkcji falowej:
\(\psi(r)=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)Gdy dwa atomy oddziałują ze sobą w separacji R, model daje efektywną energię oddziaływania przyciągającego, której dokładna postać po pełnym całkowaniu 3D jest potencjałem typu Yukawy:
\(E_{\mathrm{att}}(R)=-\frac{\kappa}{\sqrt{\pi}}e^{-R/\alpha_{\mathrm{eff}}}\)W połączeniu z odpychaniem jądrowym w jednostkach atomowych, ten dwuparametrowy model odtwarza odległość równowagi cząsteczki H₂ i energię dysocjacji po kalibracji do danych eksperymentalnych.
Kluczowy wynik oryginalnego artykułu BeeTheory został potwierdzony: oddziaływanie falowe wytwarza siłę przyciągającą. Jednak przybliżenie monopolu zostało tutaj poprawione, ponieważ traci zależność od R. Poprawiony model daje postać Yukawy ze skalibrowanymi współczynnikami.
\(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,509 Eh
Odpowiednik 95,5 eV. Ustawia amplitudę oddziaływania przyciągającego.
αeff = 1,727 a0
Odpowiada to 91,4 pm. Jest to o 72,7% więcej niż goły promień Bohra.
Błąd <0,2%
Req = 74,16 pm iDe = 4,517 eV, pasujące do eksperymentu.
1. Funkcja falowa: Dokładna forma 3D
1.1 Postulat wyjściowy teorii pszczół
Każda cząstka elementarna jest modelowana przez funkcję falową, która rozpada się wykładniczo we wszystkich trzech kierunkach przestrzennych od jej środka. W przypadku atomu wodoru w stanie podstawowym nie jest to jedynie postulat, ale dokładny wynik kwantowo-mechaniczny: funkcja falowa BeeTheory pokrywa się z orbitalem 1s wodoru.
\(\psi_{1s}(\mathbf{r})=\frac{1}{\sqrt{\pi a_0^3}}\exp\left(-\frac{r}{a_0}\right),\qquad r=|\mathbf{r}|\)W notacji kompaktowej z α = 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 Normalizacja – dokładna weryfikacja
\(\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 Energia – weryfikacja równania Schrödingera
Zastosowanie niezależnego od czasu równania Schrödingera:
\(\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}\)Dokładny Laplacian exp(-αr) we współrzędnych sferycznych wynosi:
\(\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)\)Korekta do artykułu BeeTheory
Oryginalne przybliżenie ∇²f(r) ≈ -3α/RAB odrzuca zależność radialną. Dokładny Laplacian zawiera dwa wyrazy: α²e-αr i -2αe-αr/r. Poprawiona pochodna zachowuje oba wyrazy.
W jednostkach atomowych, z ħ =me = e = 1 i 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. Suma dwóch funkcji falowych – podejście dokładne
Proszę umieścić atom A w punkcie początkowym, a atom B w pozycji R na osi z. Całkowita funkcja falowa w superpozycji BeeTheory wynosi:
\(\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 Funkcja falowa A oszacowana w pobliżu B
W pobliżu atomu B wkład fali A wynosi:
\(\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}\)AmplitudaCA(R) maleje wykładniczo wraz z separacją. Jest to sygnał BeeTheory przenoszony z atomu A do atomu B.
| R | CA(R)/N = e-R/a₀ | Znaczenie fizyczne |
|---|---|---|
| 0.5 a0 | 0.607 | Silne nakładanie się, reżim odpychający |
| 1.0 a0 | 0.368 | Przy promieniu Bohra |
| 1.4 a0 | 0.247 | Długość wiązania blisko H₂ |
| 2.0 a0 | 0.135 | Wciąż znaczące |
| 3.0 a0 | 0.050 | Reżim słabej interakcji |
| 5.0 a0 | 0.007 | Interakcja bliska zeru |
2.2 Hamiltonian zastosowany do członu krzyżowego
W pobliżu B efektywna fala lokalna wynosi:
\(\Psi_{\mathrm{local}}(r)\approx[C_A(R)+N]e^{-r/a_0}\)Zastosowanie operatora kinetycznego do wkładu A daje wynik:
\(\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)\)Wyrażenie 1/r z operatora kinetycznego paruje się z potencjałem Coulomba i przyczynia się do efektywnego przyciągania.
\(\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. Od sprzężenia kinetycznego do potencjału interakcji
3.1 Pełna interakcja BeeTheory
Interakcja BeeTheory między atomami A i B pochodzi z kinetycznego sprzężenia pola falowego A z gęstością elektronową B. W połączeniu z odpychaniem jądrowym, całkowita energia oddziaływania przyjmuje postać:
\(E_{\mathrm{BT}}(R)=-\kappa\frac{e^{-R/\alpha_{\mathrm{eff}}}}{\sqrt{\pi}}+\frac{1}{R}\)Ujemny człon to przyciąganie, a człon 1/R to odpychanie jądrowe. Dwa parametry kontrolują oddziaływanie: κ i αeff.
3.2 Porównanie z oryginalnym artykułem
Oryginalne przybliżenie
\(\nabla^2f\approx-\frac{3\alpha}{R_{AB}}\)Traci to zależność interakcji od R i nie może wytworzyć odległości równowagi.
Poprawiony dokładny Laplacian
\(\nabla^2e^{-r}=e^{-r}\left(1-\frac{2}{r}\right)\)Zachowuje to pełną zależność od r i wytwarza oddziaływanie Yukawy.
3.3 Dlaczego potencjał jest Yukawy, a nie Coulomba?
Współczynnik e-R/αeff wynika z amplitudy fali A w pozycji B. Przy dużej separacji oddziaływanie zanika wykładniczo. Sprawia to, że oddziaływanie BeeTheory w skali atomowej jest potencjałem Yukawy o skończonym zasięgu.
\(F(R)=-\frac{dE}{dR}=-\frac{\kappa}{\sqrt{\pi}\alpha_{\mathrm{eff}}}e^{-R/\alpha_{\mathrm{eff}}}+\frac{1}{R^2}\)Przy długości wiązania H₂ warunki przyciągania i odpychania równoważą się.
4. Kalibracja: Dwa warunki, dwa parametry
Istnieją dokładnie dwa wolne parametry, κ i αeff, oraz dwa ograniczenia eksperymentalne z cząsteczki H₂.
| Ograniczenie | Znaczenie fizyczne | Warunek matematyczny | Wartość eksperymentalna |
|---|---|---|---|
| Proszę | Długość wiązania | dE/dR = 0 | 74,14 pm = 1,401 a0 |
| De | Energia dysocjacji | E(∞) – E(Req) =De | 4,520 eV = 0,1660 Eh |
4.1 Rozwiązanie analityczne
Warunek 1:
\(\frac{dE}{dR}=0\quad\Longrightarrow\quad\frac{\kappa e^{-R_{\mathrm{eq}}/\alpha}}{\sqrt{\pi}\alpha}=\frac{1}{R_{\mathrm{eq}}^2}\)Warunek 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\)Dzielenie warunku 2 przez warunek 1:
\(\alpha=R_{\mathrm{eq}}+D_eR_{\mathrm{eq}}^2\)PrzyReq = 1,4014 a0 iDe = 0,1660 Eh:
\(\alpha_{\mathrm{eff}}=1.4014+0.1660(1.4014)^2=1.7274a_0\)Następnie:
\(\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 Fizyczna interpretacja parametrów
| Parametr | Wartość | Znaczenie fizyczne w BeeTheory |
|---|---|---|
| κ | 3.509 Eh | Amplituda sprzężenia fala-masa. |
| αeff | 1.727 a0 | Efektywna długość rozpadu interakcji. |
| αeff/a0 | 1.727 | Współczynnik hybrydyzacji BeeTheory. |
5. Krzywa energii potencjalnej i porównanie z eksperymentem
Sugerowany wykres: Krzywa energii potencjalnej H₂ porównująca BeeTheory, Heitler-London i eksperymentalne dane referencyjne.
Alt text: Krzywa energii potencjalnej H₂ z odległością R w angstremach na osi poziomej i energią w elektronowoltach na osi pionowej. Krzywa BeeTheory osiąga minimum w pobliżu R = 0,74 Å przy -4,52 eV, dopasowując się do eksperymentalnej odległości wiązania H₂ i energii dysocjacji.
| R (a0) | R (pm) | Ewave | Enuc | EBT | EBT (eV) | Status |
|---|---|---|---|---|---|---|
| 0.50 | 26.5 | -1.482 | +2.000 | +0.518 | +14.09 | odpychający |
| 0.80 | 42.3 | -1.246 | +1.250 | +0.004 | +0.11 | blisko zera |
| 1.00 | 52.9 | -1.110 | +1.000 | -0.110 | -2.98 | atrakcyjny |
| 1.20 | 63.5 | -0.988 | +0.833 | -0.155 | -4.22 | atrakcyjny |
| 1.401 | 74.1 | -0.880 | +0.714 | -0.166 | -4.517 | minimum |
| 1.60 | 84.7 | -0.784 | +0.625 | -0.159 | -4.33 | płytka studnia |
| 2.00 | 105.8 | -0.622 | +0.500 | -0.122 | -3.32 | rosnący |
| 3.00 | 158.8 | -0.349 | +0.333 | -0.015 | -0.42 | blisko zera |
| 5.00 | 264.6 | -0.110 | +0.200 | +0.090 | +2.46 | odpychający ogon |
BeeTheory:Req = 74,2 pm iDe = 4,52 eV według skalibrowanej konstrukcji.
Heitler-London: przewiduje większą długość wiązania i niższą energię dysocjacji.
Eksperyment:Req = 74,14 pm iDe = 4,520 eV.
6. Kompletne równania – gotowe do użycia
6.1 Funkcja falowa
\(\psi(r)=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)6.2 Dokładny Laplacian
\(\nabla^2\psi(r)=\psi(r)\left(\frac{1}{a_0^2}-\frac{2}{a_0r}\right)\)6.3 Całkowita energia interakcji
\(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 Siła pomiędzy dwoma atomami wodoru
\(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 Tabela podsumowująca parametry
| Symbol | Nazwa | Wartość | Jak ustalono |
|---|---|---|---|
| a0 | Promień Bohra | 52.918 pm | Mechanika kwantowa wodoru |
| Eh | Hartree | 27,211 eV | Definicja jednostki atomowej |
| α | Stała zaniku fali | 1/a0 | Orbital wodoru 1s |
| κ | Sprzężenie fala-masa | 3.509 Eh | Skalibrowany doReq iDe |
| αeff | Efektywna długość rozpadu | 1.727 a0 | Skalibrowany na podstawie H₂ |
| Proszę | Długość wiązania równowagi | 74.14 pm | Eksperyment |
| De | Energia dysocjacji | 4,520 eV | Eksperyment |
7. Otwarte pytania i kolejne pochodne
Od H₂ do grawitacji – problem skalowania teorii pszczół
W skali atomowej BeeTheory odtwarza chemię H₂ z κ = 3,509 Eh i αeff = 1,727 a0. W skali galaktycznej BeeTheory wykorzystuje długości koherencji mierzone w kiloparsekach. Otwartym pytaniem jest, jak długość koherencji skaluje się od układów atomowych do układów astrofizycznych.
Następna pochodna: hel i atomy wieloelektronowe
Dla helu funkcja falowa może być przybliżona jako:
\(\psi_{\mathrm{He}}(r)=Ne^{-\alpha_{\mathrm{He}}r}\)Naturalnym kolejnym krokiem jest przetestowanie BeeTheory pod kątem oddziaływań van der Waalsa z He₂.
Rozszerzenie: nieidentyczne atomy
Dla atomów A i B o różnych stałych rozpadu, ogólne oddziaływanie BeeTheory można zapisać jako:
\(E(R)=-\kappa_{AB}\frac{e^{-R/\alpha_{AB}}}{\sqrt{\pi}}+\frac{Z_AZ_B}{R}\)Referencje
- Dutertre, X. – Bee Theory™: Wave-Based Modeling of Gravity, 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. – Potential-Energy Curves for the X¹Σg⁺, b³Σu⁺, and C¹Πu States of the Hydrogen Molecule, 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 – Badanie grawitacji poprzez fizykę kwantową opartą na falach
© Technoplane S.A.S. – Treści tworzone z wykorzystaniem ludzkiej wiedzy i pomocy sztucznej inteligencji