Módulo 0 — História e Evolução da Química Quântica: Dos Quanta de Planck e Equação de Schrödinger ao Hartree–Fock e DFT Module 0 — History and Evolution of Quantum Chemistry: From Planck's Quanta and Schrödinger's Equation to Hartree–Fock and DFT

Porque é que uma luz vermelha intensíssima não arranca um único eletrão de uma placa de potássio, enquanto uma luz ultravioleta muito fraca arranca eletrões instantaneamente? Porque o efeito fotoelétrico não é uma 'onda contínua a aquecer a superfície', mas sim uma colisão individual $1:1$ entre um fotão de energia $h\nu$ e um eletrão preso pela função trabalho $\Phi$ (energia de ionização): Why does intense red light fail to eject electrons from potassium metal while faint UV light ejects them immediately? Because photoemission is a $1:1$ collision between a single photon of energy $h\nu$ and an electron bound by work function $\Phi$:

💡 Como ler esta equação com intuição física:How to read this equation with physical intuition: Repare como $\lambda = h/p$ explica imediatamente a quantização na caixa de tamanho $L$: para a onda do eletrão não se anular nas paredes, têm de caber meias-ondas inteiras ($L = n\,\lambda/2$), o que força a energia a saltar em degraus discretos $E_n \propto n^2/L^2$! Notice how $\lambda = h/p$ explains quantization in a box of length $L$: standing waves require $L = n\lambda/2$, forcing energy into discrete steps $E_n \propto n^2/L^2$!

A fórmula da partícula na caixa ($\Delta E \propto 1/L^2$) está nos seus olhos neste preciso momento! No $\beta$-caroteno da cenoura, os $22$ eletrões $\pi$ conjugados movem-se ao longo de uma 'caixa molecular 1D' de $L \approx 20\text{ Å}$, cujo salto quântico HOMO$\rightarrow$LUMO absorve luz azul-esverdeada ($\sim 450\text{ nm}$), refletindo o laranja complementar. E nos televisores QLED modernos (Pontos Quânticos, Prémio Nobel da Química 2023 para Bawendi, Brus e Ekimov), nanocristais do mesmo semicondutor ($\text{CdSe}$ ou $\text{InP}$) emitem luz vermelha quando têm $6\text{ nm}$ ($L$ grande $\Rightarrow \Delta E$ pequeno) e luz azul quando têm $2\text{ nm}$ ($L$ pequeno $\Rightarrow \Delta E$ grande)! The particle-in-a-box law ($\Delta E \propto 1/L^2$) explains why $\beta$-carotene's conjugated $\pi$-chain ($L\approx 20\text{ Å}$) absorbs blue light to make carrots look orange, and why semiconductor Quantum Dots (Nobel Prize 2023) emit red light at $6\text{ nm}$ diameter and blue light at $2\text{ nm}$!

Na mecânica clássica especificamos a posição $x(t)$ e o momento $p(t)$ de uma partícula. Na mecânica quântica, toda a informação física do sistema está codificada na função de onda $\psi(\vec{r})$, obtida resolvendo a equação aos valores próprios do operador Hamiltoniano $\hat{H} = \hat{T} + \hat{V}$: In quantum mechanics, all physical information of a stationary state is encoded in the wavefunction $\psi(\vec{r})$, obtained by solving the eigenvalue equation of the Hamiltonian operator $\hat{H} = \hat{T} + \hat{V}$:

💡 Como ler esta equação com intuição física:How to read this equation with physical intuition: Na aproximação de Born–Oppenheimer, resolvemos primeiro a equação eletrónica para cada geometria nuclear fixa $\vec{R}$. A energia resultante $U(\vec{R}) = E_{\text{el}}(\vec{R}) + V_{NN}(\vec{R})$ forma a paisagem topográfica (PES) onde os mínimos são isómeros estáveis e os pontos de sela são estados de transição! Under the Born–Oppenheimer approximation, solving the electronic problem at fixed nuclear positions $\vec{R}$ yields the topographic landscape $U(\vec{R})$ whose valleys are stable isomers and whose mountain passes are transition states!

Num átomo pesado como o Ouro ($\text{Au}, Z=79$) ou o Mercúrio ($\text{Hg}, Z=80$), a carga nuclear é tão grande que os eletrões $1s$ movem-se a cerca de $58\%$ da velocidade da luz ($v \approx 0.58c$)! Pela equação relativista de Dirac, a massa relativista do eletrão aumenta e as orbitais $s$ contraem-se fortemente, blindando melhor o núcleo e fazendo as orbitais $5d$ expandirem-se e subirem em energia. No Ouro, isto reduz o hiato $5d \rightarrow 6s$ do ultravioleta (como na Prata, $Z=47$, que reflete todas as cores visíveis) para o azul ($\sim 2.4\text{ eV}$), dando ao ouro o seu brilho amarelo! E no Mercúrio ($6s^2$), a orbital $6s$ fica tão contraída e inerte que os átomos de $\text{Hg}$ quase não formam ligação metálica, permanecendo líquidos a $25\,^\circ\text{C}$! In Gold ($Z=79$) and Mercury ($Z=80$), $1s$ electrons travel at $\sim 58\%$ of the speed of light! Relativistic contraction of the $6s$ orbital and expansion of the $5d$ shell narrow Gold's $5d\rightarrow 6s$ gap so it absorbs blue light (appearing golden yellow), while making Mercury's $6s^2$ shell so inert that $\text{Hg}$ is a liquid at room temperature!

Como sabemos se uma função de onda aproximada $\tilde{\Phi}$ é boa se não conhecemos a solução exata? Graças ao Teorema Variacional: o valor esperado da energia de qualquer função de onda tentativa normalizada é sempre superior ou igual à energia exata do estado fundamental ($E[\tilde{\Phi}] \ge E_0$). Ao minimizar a energia de um único Determinante de Slater $\Phi_{\text{SD}}$ expandido numa base de orbitais atómicas (Roothaan–Hall, 1951), obtemos o problema matricial autoconsistente: How do we optimize an approximate wavefunction $\tilde{\Phi}$? By the Variational Theorem, the expectation energy of any trial wavefunction is an upper bound to the true ground-state energy ($E[\tilde{\Phi}] \ge E_0$). Minimizing the energy of a single Slater Determinant yields the Roothaan–Hall matrix equations:

💡 Como ler esta equação com intuição física:How to read this equation with physical intuition: O método de Hartree–Fock recupera $\sim 99\%$ da energia eletrónica total, mas como trata os eletrões como se se movessem num campo médio suave, ignora a correlação instantânea entre eletrões de spins opostos. Os restantes $1\%$ — a energia de correlação eletrónica ($E_{\text{corr}}$) — têm a mesma ordem de grandeza das energias de ligação química! Hartree–Fock recovers $\sim 99\%$ of total electronic energy, but its mean-field approximation misses instantaneous correlation between opposite-spin electrons. That remaining $1\%$—the correlation energy ($E_{\text{corr}}$)—is comparable in magnitude to chemical bond energies!

Se desenhar a estrutura de Lewis do oxigénio atmosférico ($\text{O}=\text{O}$), todos os 12 eletrões de valência aparecem emparelhados — o que faria prever que o $\text{O}_2$ fosse diamagnético. Porém, já Michael Faraday tinha observado que o $\text{O}_2$ é paramagnético (quando se despeja $\text{O}_2$ líquido azul-pálido entre os polos de um íman forte, o líquido fica suspenso no ar!). Em 1928, Lennard-Jones e Mulliken mostraram que na teoria das Orbitais Moleculares os dois últimos eletrões do $\text{O}_2$ ocupam o par de orbitais antiligantes degeneradas $(\pi^*_{2p_x})^1(\pi^*_{2p_y})^1$ com spins paralelos ($S=1$, estado fundamental tripleto $^3\Sigma_g^-$) pela Regra de Hund! A simple Lewis structure ($\text{O}=\text{O}$) pairs all 12 valence electrons, falsely predicting $\text{O}_2$ to be diamagnetic. Yet pale-blue liquid $\text{O}_2$ hangs suspended between the poles of a magnet! In 1928, MO theory triumphantly explained why: the highest two electrons occupy degenerate $(\pi^*_{2p_x})^1(\pi^*_{2p_y})^1$ orbitals with parallel spins (triplet $^3\Sigma_g^-$ ground state)!

Pela teoria das perturbações de segunda ordem (equação de Klopman–Salem), quando duas moléculas se aproximam, a estabilização covalente provém sobretudo dos pares de orbitais ocupadas e virtuais cuja diferença de energia $|\varepsilon_{\text{HOMO}} - \varepsilon_{\text{LUMO}}|$ é mínima e cuja sobreposição $S_{\text{HOMO,LUMO}}$ é máxima: By second-order perturbation theory, when two reactants approach, covalent stabilization is dominated by frontier orbitals with the smallest energy gap $|\varepsilon_{\text{HOMO}} - \varepsilon_{\text{LUMO}}|$ and largest overlap $S$:

💡 Como ler esta equação com intuição física:How to read this equation with physical intuition: Para capturar a correlação eletrónica dinâmica com precisão termoquímica ('precisão química' $\lt 1\text{ kcal mol}^{-1}$), incluímos excitações duplas e triplas para orbitais virtuais no método $\text{CCSD(T)}$, cujo custo computacional cresce com $\mathcal{O}(N^7)$. To achieve 'chemical accuracy' ($\lt 1\text{ kcal mol}^{-1}$), Coupled Cluster $\text{CCSD(T)}$ includes single, double, and perturbative triple excitations into virtual orbitals at $\mathcal{O}(N^7)$ scaling.

A síntese da vitamina $\text{D}_3$ na epiderme humana começa com a abertura eletrocíclica do anel ciclo-hexadieno do 7-desidrocolesterol ($6\text{ eletrões }\pi$). À temperatura corporal no escuro (via térmica, controlada pelo HOMO $\psi_3$), a geometria do esteroide impede a rotação disrotatória necessária, pelo que a reação não ocorre. Mas quando um fotão UVB solar promove um eletrão para o LUMO ($\psi_4^*$), a simetria nos terminais inverte-se e permite a abertura conrotatória fotoquímica instantânea em pré-vitamina $\text{D}_3$, seguida de um rearranjo sigmatrópico $[1,7]$ térmico antarafacial! Vitamin $\text{D}_3$ synthesis in human skin relies directly on Woodward–Hoffmann orbital symmetry: UVB sunlight promotes an electron from $\psi_3$ (HOMO) to $\psi_4^*$ (LUMO) in 7-dehydrocholesterol, unlocking a photochemical conrotatory $6\pi$ electrocyclic ring opening that is sterically impossible in the dark!

Na formulação de Kohn–Sham (1965), calculamos a energia cinética exata $T_s[\rho]$ de um sistema fictício de eletrões não-interagentes que possui exatamente a mesma densidade $\rho(\vec{r})$ do sistema real, concentrando todos os efeitos quânticos de muitos corpos num único termo: o funcional de troca-correlação $E_{\text{XC}}[\rho]$: In Kohn–Sham DFT (1965), a reference system of non-interacting orbitals reproduces the exact density $\rho(\vec{r})$, concentrating all many-body quantum effects into the exchange–correlation functional $E_{\text{XC}}[\rho]$:

💡 Como ler esta equação com intuição física:How to read this equation with physical intuition: Como os funcionais semilocais (LDA/GGA) não anulam perfeitamente a repulsão de Coulomb $J[\rho]$ de um eletrão consigo próprio (Erro de Auto-Interação, SIE), tendem a sobre-deslocalizar eletrões e a subestimar barreiras de reação ($\text{H}_3^\ddagger$), e não capturam forças de London a longa distância ($-C_6/R^6$) sem a correção de dispersão `D3(BJ)` ou funcionais híbridos range-separated! Because semilocal functionals (LDA/GGA) suffer from Self-Interaction Error (SIE), they underestimate reaction barriers and miss long-range London dispersion ($-C_6/R^6$) unless augmented with exact exchange (hybrids) and `D3(BJ)` dispersion corrections!

Dentro de cada folha de grafeno na ponta de um lápis, os átomos de carbono estão unidos por ligações covalentes $\sigma/\pi$ fortíssimas ($\sim 500\text{ kJ mol}^{-1}$). Mas entre folhas adjacentes (separadas por $3.35\text{ Å}$), a coesão provém exclusivamente de forças de dispersão de London! Um cálculo `B3LYP` puro sem correção de dispersão prevê erradamente que a grafite se desfaz sozinha (repulsiva!), enquanto `B3LYP-D3(BJ)` reproduz exatamente a energia de esfoliação experimental ($\sim 52\text{ meV/átomo}$) que permite à ponta do lápis deixar camadas de grafeno no papel! Within each graphene sheet of a pencil, carbon atoms are held by strong covalent bonds, while adjacent sheets ($3.35\text{ Å}$ apart) are held purely by London dispersion! Uncorrected `B3LYP` falsely predicts graphite sheets to repel each other, whereas `B3LYP-D3(BJ)` accurately recovers the interlayer binding energy that lets pencils write on paper!