Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2511.01132

Screening in the Heitler-London Model: Revisiting the Bonding and Antibonding States of the Hydrogen Molecule

Washington P. da Silva, Daniel Vieira, Jonas Maziero, Edgard P. M. Amorim

physics.chem-phquant-phSubmitted Nov 3, 2025• Announced Nov 4, 2025
arXivHTMLPDF

Abstract

The present manuscript revisits one of the earliest approaches to treating molecular systems within the Schrödinger formalism of quantum mechanics: the Heitler-London (HL) model. Originally proposed in 1927 and based on a linear combination of atomic orbitals, the HL model provided a foundational description of covalent bonds and has served as the basis for numerous variational methods. Focusing on the hydrogen molecule, we begin by revisiting the analytical calculations of the original HL model, from which the qualitative physics of bonding and antibonding states can be obtained. Subsequently, we propose including electronic screening effects directly in the original HL wave function. We then compare our proposal with variational quantum Monte Carlo (VQMC) calculations, whose trial wave function allows us to optimize the electronic screening potential as a function of the inter-proton distance. We obtain the bond length, binding energy, and vibrational frequency of the H$_2$ molecule. Beyond revisiting this foundational approach in quantum mechanics, our proposal can serve as improved input for constructing new, but still analytically simple, variational wave functions to describe dissociation or bond formation.

Journal reference
Braz. J. Phys. 56, 9 (2026)
DOI
10.1007/s13538-025-01935-2
Comments
18 pages, 5 figures, one column

Related papers

Molecular Properties by Quantum Monte Carlo: An Investigation on the Role of the Wave Function Ansatz and the Basis Set in the Water Molecule

Andrea Zen, Luo, Y, Sorella, S +1 more

Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely: the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudo potential, and the basis set for QMC calculations. We also introduce a new strategy for the definition of the atomic orbitals involved in the Jastrow - Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.

•Oct 1, 2013•37 citations

Molecular Properties by Quantum Monte Carlo: An Investigation on the Role of the Wave Function Ansatz and the Basis Set in the Water Molecule

Andrea Zen, Ye Luo, Sandro Sorella +1 more

Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.

arXiv:1309.0427•Journal of Chemical Theory and Computation•Aug 28, 2013•44 citations

Developing density functional theory for Bose-Einstein condensates. The case of chemical bonding

Mihai V. Putz

Since the nowadays growing interest in Bose-Einstein condensates due to the expanded experimental evidence on various atomic systems within optical lattices in weak and strong coupling regimes, the connection with Density Functional Theory is firstly advanced within the mean field framework at three levels of comprehension: the many-body normalization condition, Thomas-Fermi limit, and the chemical hardness closure with the inter-bosonic strength and universal Hohenberg-Kohn functional. As an application the traditional Heitler-London quantum mechanical description of the chemical bonding for homopolar atomic systems is reloaded within the non-linear Schrödinger (Gross-Pitaevsky) Hamiltonian; the results show that a two-fold energetic solution is registered either for bonding and antibonding states, with the bosonic contribution being driven by the square of the order parameter for the Bose-Einstein condensate density in free (gas) motion, while the associate wave functions remain as in classical molecular orbital model.

AIP conference proceedings•Jan 1, 2015•0 citations

Chemical bond and entanglement of electrons in the hydrogen molecule

Nikos Iliopoulos, Andreas Terzis

We theoretically investigate the quantum correlations (in terms of concurrence of indistinguishable electrons) in a prototype molecular system (hydrogen molecule). With the assistance of the standard approximations of the linear combination of atomic orbitals and the configuration interaction methods we describe the electronic wavefunction of the ground state of the H 2 molecule. Moreover, we managed to find a rather simple analytic expression for the concurrence (the most used measure of quantum entanglement) of the two electrons when the molecule is in its lowest energy. We have found that concurrence does not really show any relation to the construction of the chemical bond.

arXiv:1408.0099•International Journal of Quantum Information•Aug 1, 2014•1 citation

Exchange

Jean‐Philippe Ansermet

In 1927, while working in Zurich with F. London, W. Heitler used the newly established quantum mechanics to investigate the electronic structure of the hydrogen molecule. In 1941, W. Heitler became a professor at the Dublin Institute for Advanced Studies, where E. Schroedinger served as director. By 1949, W. Heitler had transitioned to the University of Zurich, where he continued his career in theoretical physics.

Spintronics•Jun 28, 2024•0 citations

Wrong Potential Energy Term in Schrödinger’s Equation for Hydrogen Atom

Gurcharn Singh Sandhu

Even after nine decades of successful run of the Quantum Mechanics (QM), different viewpoints on foundational problems of Quantum Physics are still being actively debated. That is because mathematical logic of QM often defies the physical intuition which constitutes the main spirit of Physics. De Broglie’s hypothesis of matter waves implied that the dynamic characteristics of a micro particle in motion, can be ascribed to the wave characteristics of the wavelet accompanying the particle. The Schrödinger equation models the matter-wave interactions through wavefunction ψ and effectively serves as the foundation of QM. Even though mathematical structure of the Schrödinger equation is sound and elegant, here we show a conceptual mistake in the development of this equation wherein the physical situation has not been correctly modeled in the equation. The Coulomb potential energy of the proton electron pair in Hydrogen atom is essentially the negative interaction energy between their superposed electrostatic fields which is inversely proportional to their instantaneous separation distance. Assuming the proton to be relatively fixed at the origin of an appropriate coordinate system, the potential energy of the orbiting electron will be a function of instantaneous position coordinates of the electron. This has not been properly modeled in the Schrödinger equation. The resulting errors in the solution have been quantitatively demonstrated in this paper. We have stressed the necessity of incorporating a specific correction in the potential energy term of the Schrödinger equation, after which it may facilitate the adoption of Bohmian QM.

Journal of Modern Physics•Jan 1, 2018•3 citations
The present manuscript revisits one of the earliest approaches to treating molecular systems within the Schrödinger formalism of quantum mechanics: the Heitler-London (HL) model. Originally proposed in 1927 and based on a linear combination of atomic orbitals, the HL model provided a foundational description of covalent bonds and has served as the basis for numerous variational methods. Focusing on the hydrogen molecule, we begin by revisiting the analytical calculations of the original HL model, from which the qualitative physics of bonding and antibonding states can be obtained. Subsequently, we propose including electronic screening effects directly in the original HL wave function. We then compare our proposal with variational quantum Monte Carlo (VQMC) calculations, whose trial wave function allows us to optimize the electronic screening potential as a function of the inter-proton distance. We obtain the bond length, binding energy, and vibrational frequency of the H2_22​ molecule. Beyond revisiting this foundational approach in quantum mechanics, our proposal can serve as improved input for constructing new, but still analytically simple, variational wave functions to describe dissociation or bond formation.
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely: the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudo potential, and the basis set for QMC calculations. We also introduce a new strategy for the definition of the atomic orbitals involved in the Jastrow - Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.
Since the nowadays growing interest in Bose-Einstein condensates due to the expanded experimental evidence on various atomic systems within optical lattices in weak and strong coupling regimes, the connection with Density Functional Theory is firstly advanced within the mean field framework at three levels of comprehension: the many-body normalization condition, Thomas-Fermi limit, and the chemical hardness closure with the inter-bosonic strength and universal Hohenberg-Kohn functional. As an application the traditional Heitler-London quantum mechanical description of the chemical bonding for homopolar atomic systems is reloaded within the non-linear Schrödinger (Gross-Pitaevsky) Hamiltonian; the results show that a two-fold energetic solution is registered either for bonding and antibonding states, with the bosonic contribution being driven by the square of the order parameter for the Bose-Einstein condensate density in free (gas) motion, while the associate wave functions remain as in classical molecular orbital model.
We theoretically investigate the quantum correlations (in terms of concurrence of indistinguishable electrons) in a prototype molecular system (hydrogen molecule). With the assistance of the standard approximations of the linear combination of atomic orbitals and the configuration interaction methods we describe the electronic wavefunction of the ground state of the H 2 molecule. Moreover, we managed to find a rather simple analytic expression for the concurrence (the most used measure of quantum entanglement) of the two electrons when the molecule is in its lowest energy. We have found that concurrence does not really show any relation to the construction of the chemical bond.
Even after nine decades of successful run of the Quantum Mechanics (QM), different viewpoints on foundational problems of Quantum Physics are still being actively debated. That is because mathematical logic of QM often defies the physical intuition which constitutes the main spirit of Physics. De Broglie’s hypothesis of matter waves implied that the dynamic characteristics of a micro particle in motion, can be ascribed to the wave characteristics of the wavelet accompanying the particle. The Schrödinger equation models the matter-wave interactions through wavefunction ψ and effectively serves as the foundation of QM. Even though mathematical structure of the Schrödinger equation is sound and elegant, here we show a conceptual mistake in the development of this equation wherein the physical situation has not been correctly modeled in the equation. The Coulomb potential energy of the proton electron pair in Hydrogen atom is essentially the negative interaction energy between their superposed electrostatic fields which is inversely proportional to their instantaneous separation distance. Assuming the proton to be relatively fixed at the origin of an appropriate coordinate system, the potential energy of the orbiting electron will be a function of instantaneous position coordinates of the electron. This has not been properly modeled in the Schrödinger equation. The resulting errors in the solution have been quantitatively demonstrated in this paper. We have stressed the necessity of incorporating a specific correction in the potential energy term of the Schrödinger equation, after which it may facilitate the adoption of Bohmian QM.