Quantum Correlations in Frustrated Three-Body Systems

Technical Summary: Quantum Correlations in Frustrated Three-Body Systems

Problem Statement
Many physical systems involving multi-body quantum correlations remain difficult to understand from first principles, particularly because their Schrödinger equations lack closed-form solutions. Classical simulations of these systems often fail to capture unique quantum features, leading to discrepancies labeled as "anomalies." The authors focus on three specific frustrated three-body systems where competing interactions (attractive vs. repulsive forces) prevent a simple separation of variables. These systems are:

1. Helium-like atoms: Two electrons bound to a nucleus, featuring two attractive nucleus-electron forces and one repulsive electron-electron force.
2. The hydrogen molecular ion ($H_2^+$): The simplest molecule, consisting of two protons and one electron.
3. The hydrogen bond (O-H-O): Modeled as a proton between two oxygen atoms, structurally analogous to $H_2^+$ but with flipped charge signs and different mass scales.

The central challenge is to characterize the ground state wavefunctions and the nature of quantum entanglement arising from these frustrating interactions using minimal variational parameters and analytical methods.

Methodology
The paper employs a combination of variational methods, analytical series expansions, and numerical fitting:

• Helium-like Atoms: The authors utilize the Rayleigh-Ritz variational principle to estimate ground state energies. They test four specific ansätze (trial wavefunctions) ranging from simple products of hydrogenic orbitals to more complex forms incorporating "permanent" symmetrization (required for identical fermions in a spin-singlet state) and explicit electron-electron repulsion terms ($1 + cr_{12}$). They analyze the system by introducing a frustration parameter $\lambda$ to scale the electron-electron repulsion, allowing them to track the transition from symmetric (non-frustrated) to asymmetric (frustrated) solutions. Entanglement is quantified using the von Neumann entropy derived from the Schmidt decomposition of the spatial wavefunction.
• Hydrogen Molecular Ion ($H_2^+$): The authors solve the Schrödinger equation using prolate spheroidal coordinates, where the equation separates into angular and radial components. They derive series solutions using associated Legendre polynomials and power series, linking the separation constants via recurrence relations. The exact solution is expressed in terms of Confluent Heun functions. To facilitate visualization and approximation, they also fit the exact solution to finite-order polynomial and Chebyshev polynomial expansions.
• Hydrogen Bond: The authors propose a model treating the oxygen atoms as static point charges with an effective charge $Ze$. By rescaling the Hamiltonian to match the units of the $H_2^+$ system, they adapt the analytical solutions derived for the molecular ion. They perform a simultaneous fit of the analytical solution to experimental data regarding the bond length and dissociation energy of a water dimer, varying the effective charge parameter $Ze$ to minimize the difference between theory and experiment.

Key Results

• Helium-like Atoms:

  • The study demonstrates that frustration (the competition between attraction and repulsion) drives the system into an entangled state. When the electron-electron repulsion is active ($\lambda > 0$), the optimal variational parameters break symmetry ($\alpha \neq \beta$), resulting in a wavefunction described by a permanent rather than a simple product state.
  • The entanglement entropy ($S_{vN}$) is zero for non-frustrated systems ($\lambda \le 0$) and increases with frustration. However, it cannot reach the maximum theoretical value because the orbitals cannot become perfectly orthogonal.
  • Ansatz 4 (incorporating both permanent structure and repulsion) provides the most accurate energy estimates, with fractional accuracy improving as the atomic number $Z$ increases. The study finds that for small $Z$, permanent symmetrization and Coulomb repulsion contribute comparably to the ground state energy.

• Hydrogen Molecular Ion ($H_2^+$):

  • The ground state wavefunction is analytically expressed as a product of Confluent Heun functions.
  • The probability distribution reveals a distinct double-peaked structure, indicating a high probability of finding the electron near one of the nuclei rather than in the center. This structure is interpreted as evidence of quantum tunneling between the two nuclei, a fundamental mechanism of chemical bonding.
  • Approximations using Chebyshev polynomials converge faster than standard power series expansions, achieving 6-decimal point accuracy with low-order terms ($N=4$).

• Hydrogen Bond:

  • By mapping the O-H-O system to the $H_2^+$ model, the authors obtain a wavefunction that fits experimental dissociation energy and bond length data reasonably well, despite the simplification of treating oxygen atoms as point charges.
  • The resulting wavefunction exhibits a double-peaked structure similar to $H_2^+$, but with sharper peaks closer to the oxygen atoms. This suggests the proton is significantly more likely to be localized near one oxygen atom than in the intervening region, differing from the more rounded distribution seen in $H_2^+$.

Significance and Claims
The paper claims that studying these simple yet non-trivial three-body systems provides fundamental insights into the nature of quantum correlations and frustration.

• Entanglement Origin: It explicitly links the emergence of quantum entanglement in helium-like atoms to the frustration caused by competing Coulomb forces, showing that symmetry breaking in the wavefunction is a direct consequence of this frustration.
• Bonding Mechanism: It demonstrates that the chemical bond in $H_2^+$ and the hydrogen bond can be understood through the lens of quantum tunneling, visualized via the double-peaked probability distribution of the wavefunction.
• Modeling Utility: The work suggests that efficient modeling of "permanent" structures (as opposed to simple determinants) is crucial for understanding larger quantum systems. Furthermore, it posits that the hydrogen bond can be effectively modeled using the analytical framework of the hydrogen molecular ion, offering a first-principles approach that captures the essential 3D structure of the bond, improving upon models that view the bond merely as tunneling along a single reaction coordinate.

The authors maintain a modest tone, acknowledging limitations such as the point-charge approximation for oxygen atoms and the presence of truncation errors in the numerical fits, while emphasizing that the models yield sensible wavefunctions with a small number of phenomenological parameters.