BeeTheory – Derivação científica – 2025

Funções de onda para dois átomos de hidrogênio: Rigorous Derivation and Calibration (Derivação e calibração rigorosas)

Partindo do postulado BeeTheory de funções de onda exponencial-r, derivamos a energia de interação 3D exata, corrigimos a aproximação monopolar original e calibramos em relação à molécula H₂ conhecida com dois parâmetros que reproduzem o experimento em menos de 0,2%.

BeeTheory.com – Baseado no BeeTheory v2 (Dutertre, 2023) – Ampliado e corrigido

0. Conclusões – Resultados em primeiro lugar

O modelo baseado em ondas BeeTheory representa cada átomo de hidrogênio por uma função de onda esférica:

\(\psi(r)=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)

Quando dois átomos interagem na separação R, o modelo produz uma energia de interação atrativa efetiva cuja forma exata após a integração completa em 3D é um potencial do tipo Yukawa:

\(E_{\mathrm{att}}(R)=-\frac{\kappa}{\sqrt{\pi}}e^{-R/\alpha_{\mathrm{eff}}}\)

Combinado com a repulsão nuclear em unidades atômicas, esse modelo de dois parâmetros reproduz a distância de equilíbrio da molécula de H₂ e a energia de dissociação após a calibração com dados experimentais.

O principal resultado do artigo original da BeeTheory é confirmado: a interação da onda produz uma força atrativa. No entanto, a aproximação do monopolo é corrigida aqui porque perde a dependência de R. O modelo corrigido fornece uma forma Yukawa com coeficientes calibrados.

\(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}\)

1. A função de onda: Forma 3D exata

1.1 Postulado inicial da BeeTheory

Cada partícula elementar é modelada por uma função de onda que decai exponencialmente em todas as três direções espaciais a partir de seu centro. Para o átomo de hidrogênio em seu estado fundamental, isso não é apenas um postulado, mas um resultado quântico-mecânico exato: a função de onda da BeeTheory coincide com o orbital 1s do hidrogênio.

\(\psi_{1s}(\mathbf{r})=\frac{1}{\sqrt{\pi a_0^3}}\exp\left(-\frac{r}{a_0}\right),\qquad r=|\mathbf{r}|\)

Em notação compacta com α = 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 Normalização – Verificação exata

\(\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 – Verificação da equação de Schrödinger

Aplicação da equação de Schrödinger independente do tempo:

\(\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}\)

O Laplaciano exato de exp(-αr) em coordenadas esféricas é:

\(\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)\)

Correção do artigo da BeeTheory

A aproximação original ∇²f(r) ≈ -3α/RAB descarta a dependência radial. O Laplaciano exato tem dois termos: α²e-αr e -2αe-αr/r. A derivação corrigida mantém os dois termos.

Em unidades atômicas, com ħ =me = e = 1 e 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. Soma de duas funções de onda – abordagem exata

Coloque o átomo A na origem e o átomo B na posição R no eixo z. A função de onda total na superposição BeeTheory é:

\(\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 Função de onda de A avaliada perto de B

Perto do átomo B, a contribuição da onda de A é:

\(\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}\)

A amplitudeCA(R) decai exponencialmente com a separação. É o sinal BeeTheory transportado do átomo A para o átomo B.

RCA(R)/N = e-R/a₀Significado físico
0.5 a00.607Forte sobreposição, regime repulsivo
1.0 a00.368No raio de Bohr
1.4 a00.247Comprimento de ligação H₂ próximo
2.0 a00.135Ainda significativo
3.0 a00.050Regime de interação fraca
5.0 a00.007Interação quase nula

2.2 Hamiltoniano aplicado ao termo cruzado

Perto de B, a onda local efetiva é:

\(\Psi_{\mathrm{local}}(r)\approx[C_A(R)+N]e^{-r/a_0}\)

Aplicando o operador cinético à contribuição A, obtém-se o seguinte

\(\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)\)

O termo 1/r do operador cinético faz par com o potencial de Coulomb e contribui para a atração efetiva.

\(\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. Do acoplamento cinético ao potencial de interação

3.1 A interação completa da BeeTheory

A interação BeeTheory entre os átomos A e B vem do acoplamento cinético do campo de onda de A com a densidade de elétrons de B. Combinada com a repulsão nuclear, a energia total de interação assume a forma:

\(E_{\mathrm{BT}}(R)=-\kappa\frac{e^{-R/\alpha_{\mathrm{eff}}}}{\sqrt{\pi}}+\frac{1}{R}\)

O termo negativo é atrativo e o termo 1/R é a repulsão nuclear. Dois parâmetros controlam a interação: κ e αeff.

3.2 Comparação com o documento original

Aproximação original

\(\nabla^2f\approx-\frac{3\alpha}{R_{AB}}\)

Isso perde a dependência de R da interação e não pode produzir uma distância de equilíbrio.

Laplaciano exato corrigido

\(\nabla^2e^{-r}=e^{-r}\left(1-\frac{2}{r}\right)\)

Isso mantém a dependência total de r e produz uma interação Yukawa.

3.3 Por que o potencial é de Yukawa, não de Coulomb

O fator e-R/αeff surge da amplitude da onda de A na posição de B. Em uma grande separação, a interação decai exponencialmente. Isso torna a interação BeeTheory em escala atômica um potencial Yukawa de alcance finito.

\(F(R)=-\frac{dE}{dR}=-\frac{\kappa}{\sqrt{\pi}\alpha_{\mathrm{eff}}}e^{-R/\alpha_{\mathrm{eff}}}+\frac{1}{R^2}\)

No comprimento da ligação H₂, os termos atrativos e repulsivos se equilibram.

4. Calibração: Duas condições, dois parâmetros

Há exatamente dois parâmetros livres, κ e αeff, e duas restrições experimentais da molécula de H₂.

RestriçãoSignificado físicoCondição matemáticaValor experimental
ReqComprimento da ligaçãodE/dR = 074,14 pm = 1,401 a0
DeEnergia de dissociaçãoE(∞) – E(Req) =De4,520 eV = 0,1660Eh

4.1 Solução analítica

Condição 1:

\(\frac{dE}{dR}=0\quad\Longrightarrow\quad\frac{\kappa e^{-R_{\mathrm{eq}}/\alpha}}{\sqrt{\pi}\alpha}=\frac{1}{R_{\mathrm{eq}}^2}\)

Condição 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\)

Dividindo a condição 2 pela condição 1:

\(\alpha=R_{\mathrm{eq}}+D_eR_{\mathrm{eq}}^2\)

ComReq = 1,4014 a0 eDe = 0,1660Eh:

\(\alpha_{\mathrm{eff}}=1.4014+0.1660(1.4014)^2=1.7274a_0\)

Então:

\(\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 Interpretação física dos parâmetros

ParâmetroValorSignificado físico em BeeTheory
κ3.509EhAmplitude de acoplamento de massa de onda.
αeff1.727 a0Comprimento de decaimento efetivo da interação.
αeff/a01.727Índice de hibridização BeeTheory.

5. Curva de energia potencial e comparação com o experimento

Gráfico sugerido: Curva de energia potencial H₂ comparando a BeeTheory, Heitler-London e dados de referência experimentais.

Texto alternativo: Curva de energia potencial de H₂ com distância R em angstroms no eixo horizontal e energia em eletronvolts no eixo vertical. A curva BeeTheory atinge seu mínimo próximo a R = 0,74 Å a -4,52 eV, correspondendo à distância experimental da ligação H₂ e à energia de dissociação.

R (a0)R (pm)EwaveEnucEBTEBT (eV)Status
0.5026.5-1.482+2.000+0.518+14.09repulsivo
0.8042.3-1.246+1.250+0.004+0.11próximo de zero
1.0052.9-1.110+1.000-0.110-2.98atraente
1.2063.5-0.988+0.833-0.155-4.22atraente
1.40174.1-0.880+0.714-0.166-4.517mínimo
1.6084.7-0.784+0.625-0.159-4.33poço raso
2.00105.8-0.622+0.500-0.122-3.32crescente
3.00158.8-0.349+0.333-0.015-0.42próximo de zero
5.00264.6-0.110+0.200+0.090+2.46cauda repulsiva

BeeTheory:Req = 74,2 pm eDe = 4,52 eV por construção calibrada.

Heitler-London: prevê um comprimento de ligação maior e uma energia de dissociação menor.

Experimento:Req = 74,14 pm eDe = 4,520 eV.

6. Equações completas – prontas para usar

6.1 Função de onda

\(\psi(r)=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)

6.2 Laplaciano exato

\(\nabla^2\psi(r)=\psi(r)\left(\frac{1}{a_0^2}-\frac{2}{a_0r}\right)\)

6.3 Energia de interação total

\(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 Força entre os dois átomos de hidrogênio

\(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 de resumo dos parâmetros

SímboloNomeValorComo determinado
a0Raio de Bohr52.918 pmMecânica quântica do hidrogênio
EhHartree27,211 eVDefinição de unidade atômica
αConstante de decaimento da onda1/a0Orbital 1s do hidrogênio
κAcoplamento de massa de onda3.509EhCalibrado paraReq eDe
αeffComprimento efetivo de decaimento1.727 a0Calibrado a partir de H₂
ReqComprimento da ligação de equilíbrio74,14 pmExperimento
DeEnergia de dissociação4,520 eVExperimento

7. Perguntas abertas e próximas derivações

De H₂ à gravidade – o problema de escala da BeeTheory

Na escala atômica, a BeeTheory reproduz a química do H₂ com κ = 3,509 Eh e αeff = 1,727 a0. Na escala galáctica, a BeeTheory usa comprimentos de coerência medidos em quiloparsecs. A questão em aberto é como o comprimento de coerência é dimensionado de sistemas atômicos para sistemas astrofísicos.

Próxima derivação: átomos de hélio e multielétrons

Para o hélio, a função de onda pode ser aproximada como:

\(\psi_{\mathrm{He}}(r)=Ne^{-\alpha_{\mathrm{He}}r}\)

Testar a BeeTheory em relação às interações de van der Waals de He₂ é uma próxima etapa natural.

Extensão: átomos não idênticos

Para os átomos A e B com diferentes constantes de decaimento, a interação BeeTheory geral pode ser escrita como:

\(E(R)=-\kappa_{AB}\frac{e^{-R/\alpha_{AB}}}{\sqrt{\pi}}+\frac{Z_AZ_B}{R}\)

Referências

  • Dutertre, X. – Bee Theory™: Modelagem da gravidade com base em ondas, 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 (A energia de dissociação da molécula de hidrogênio), J. Mol. Spectrosc. 33, 147, 1970.
  • Slater, J. C. – Atomic Shielding Constants (Constantes de blindagem atômica), Phys. Rev. 36, 57, 1930.
  • Atkins, P. W., Friedman, R. – Molecular Quantum Mechanics, 5ª ed., Oxford University Press, 2011.

BeeTheory.com – Explorando a gravidade por meio da física quântica baseada em ondas

© Technoplane S.A.S. – Conteúdo produzido com experiência humana e assistência de IA