Finite-temperature electronic simulations without the Born-Oppenheimer constraint

Quantum computing for chemistry and physics applications from a Monte Carlo perspective

This Perspective focuses on the several overlaps between quantum algorithms and Monte Carlo methods in the domains of physics and chemistry. We will analyze the challenges and possibilities of integrating established quantum Monte Carlo solutions into quantum algorithms. These include refined energy estimators, parameter optimization, real and imaginary-time dynamics, and variational circuits. Conversely, we will review new ideas for utilizing quantum hardware to accelerate the sampling in statistical classical models, with applications in physics, chemistry, optimization, and machine learning. This review aims to be accessible to both communities and intends to foster further algorithmic developments at the intersection of quantum computing and Monte Carlo methods. Most of the works discussed in this Perspective have emerged within the last two years, indicating a rapidly growing interest in this promising area of research.

Electronic Free Energy Surface of the Nitrogen Dimer Using First-Principles Finite Temperature Electronic Structure Methods

We use full configuration interaction and density matrix quantum Monte Carlo methods to calculate the electronic free energy surface of the nitrogen dimer within the free-energy Born-Oppenheimer approximation. As the temperature is raised from T = 0, we find a temperature regime in which the internal energy causes bond strengthening. At these temperatures, adding in the entropy contributions is required to cause the bond to gradually weaken with increasing temperature. We predict a thermally driven dissociation for the nitrogen dimer between 22,000 to 63,200 K depending on symmetries and basis set. Inclusion of more spatial and spin symmetries reduces the temperature required. The origin of these observations is explored using the structure of the density matrix at various temperatures and bond lengths.

Classical Quantum Friction at Water-Carbon Interfaces

Friction at water-carbon interfaces remains a major puzzle with theories and simulations unable to explain experimental trends in nanoscale waterflow. A recent theoretical framework─quantum friction (QF)─proposes to resolve these experimental observations by considering nonadiabatic coupling between dielectric fluctuations in water and graphitic surfaces. Here, using a classical model that enables fine-tuning of the solid's dielectric spectrum, we provide evidence from simulations in general support of QF. In particular, as features in the solid's dielectric spectrum begin to overlap with water's librational and Debye modes, we find an increase in friction in line with that proposed by QF. At the microscopic level, we find that this contribution to friction manifests more distinctly in the dynamics of the solid's charge density than that of water. Our findings suggest that experimental signatures of QF may be more pronounced in the solid's response rather than liquid water's.

Accelerating ab initio Molecular Dynamics and Probing the Weak Dispersive Forces in Dense Liquid Hydrogen

We propose an ab initio molecular dynamics method, capable of dramatically reducing the autocorrelation time required for the simulation of classical and quantum particles at finite temperatures. The method is based on an efficient implementation of a first order Langevin dynamics modified by means of a suitable, position dependent acceleration matrix S. Here, we apply this technique to both Lennard-Jones models, to demonstrate the accuracy and speeding-up of the sampling, and within a quantum Monte Carlo based wave function approach, for determining the phase diagram of high-pressure hydrogen with simulations much longer than the autocorrelation time. With the proposed method, we are able to equilibrate in a few hundred steps even close to the liquid-liquid phase transition (LLT). Within our approach, we find that the LLT transition is consistent with recent density functionals predicting a much larger transition pressure when the long range dispersive forces are taken into account.

Correlation Effects in Scanning Tunneling Microscopy Images of Molecules Revealed by Quantum Monte Carlo

Scanning tunneling microscopy (STM) and spectroscopy probe the local density of states of single molecules electrically insulated from the substrate. The experimental images, although usually interpreted in terms of single-particle molecular orbitals, are associated with quasiparticle wave functions dressed by the whole electron-electron interaction. Here we propose an ab initio approach based on quantum Monte Carlo to calculate the quasiparticle wave functions of molecules. Through the comparison between Monte Carlo wave functions and their uncorrelated Hartree-Fock counterparts we visualize the electronic correlation embedded in the simulated STM images, highlighting the many-body features that might be observed.

Quantum Monte Carlo Treatment of the Charge Transfer and Diradical Electronic Character in a Retinal Chromophore Minimal Model

The penta-2,4-dieniminium cation (PSB3) displays similar ground state and first excited state potential energy features as those of the retinal protonated Schiff base (RPSB) chromophore in rhodopsin. Recently, PSB3 has been used to benchmark several electronic structure methods, including highly correlated multireference wave function approaches, highlighting the necessity to accurately describe the electronic correlation in order to obtain reliable properties even along the ground state (thermal) isomerization paths. In this work, we apply two quantum Monte Carlo approaches, the variational Monte Carlo and the lattice regularized diffusion Monte Carlo, to study the energetics and electronic properties of PSB3 along representative minimum energy paths and scans related to its thermal cis–trans isomerization. Quantum Monte Carlo is used in combination with the Jastrow antisymmetrized geminal power ansatz, which guarantees an accurate and balanced description of the static electronic correlation thanks to the multiconfigurational nature of the antisymmetrized geminal power term, and of the dynamical correlation, due to the presence of the Jastrow factor explicitly depending on electron–electron distances. Along the two ground state isomerization minimum energy paths of PSB3, CASSCF calculations yield wave functions having either charge transfer or diradical character in proximity of the two transition state configurations. Here, we observe that at the quantum Monte Carlo level of theory, only the transition state with charge transfer character can be located. The conical intersection, which becomes highly sloped, is observed only if the path connecting the two original CASSCF transition states is extended beyond the diradical one, namely by increasing the bond-length-alternation (BLA). These findings are in good agreement with the results obtained by MRCISD+Q calculations, and they demonstrate the importance of having an accurate description of the static and dynamical correlation when studying isomerization and transition states of conjugated systems.

Vibronic Theory of Ultrafast Intersystem Crossing Dynamics in a Single Spin-Crossover Molecule at Finite Temperature beyond the Born--Oppenheimer Approximation

Quantum density matrix theory is carried out to study the ultrafast dynamics of the photoinduced state in a spin-crossover (SC) molecule interacting with a heat bath. The investigations are realized at finite temperature and beyond the usual Born-Oppenheimer (BO) approach. We found that the SC molecule experiences in the photoexcited state (PES) a huge internal pressure, estimated at several gigapascals, partly released in an "explosive" way within ∼100 fs, causing large bond length oscillations, which dampen in the picosecond time scale because of internal conversion processes. During this regime, the BO approximation is not valid. Depending on the tunneling strength, the ultrafast relaxation may proceed through the thermodynamic metastable high-spin state or prevent it. Interestingly, we demonstrate that final relaxation toward the low-spin state always follows a local equilibrium pathway, where the BO approach is valid. Our formulation reconciles the nonequilibrium and the equilibrium properties of this fascinating phenomenon and opens the way to quantum studies on cluster molecules.

Ground State Geometries of Polyacetylene Chains from Many-Particle Quantum Mechanics

Due to the crucial role played by electron correlation, the accurate determination of ground state geometries of π-conjugated molecules is still a challenge for many quantum chemistry methods. Because of the high parallelism of the algorithms and their explicit treatment of electron correlation effects, Quantum Monte Carlo calculations can offer an accurate and reliable description of the electronic states and of the geometries of such systems, competing with traditional quantum chemistry approaches. Here, we report the structural properties of polyacetylene chains H-(C₂H₂)(N)-H up to N = 12 acetylene units, by means of Variational Monte Carlo (VMC) calculations based on the multi-determinant Jastrow Antisymmetrized Geminal Power (JAGP) wave function. This compact ansatz can provide for such systems an accurate description of the dynamical electronic correlation as recently detailed for the 1,3-butadiene molecule [J. Chem. Theory Comput. 2015 11 (2), 508-517]. The calculated Bond Length Alternation (BLA), namely the difference between the single and double carbon bonds, extrapolates, for N → ∞, to a value of 0.0910(7) Å, compatible with the experimental data. An accurate analysis was able to distinguish between the influence of the multi-determinantal AGP expansion and of the Jastrow factor on the geometrical properties of the fragments. Our size-extensive and self-interaction-free results provide new and accurate ab initio references for the structures of the ground state of polyenes.

Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo

Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.

Distinct metallization and atomization transitions in dense liquid hydrogen

We perform molecular dynamics simulations driven by accurate quantum Monte Carlo forces on dense liquid hydrogen. There is a recent report of a complete atomization transition between a mixed molecular-atomic liquid and a completely dissociated fluid in an almost unaccessible pressure range [Nat. Commun. 5, 3487 (2014)]. Here, instead, we identify a different transition between the fully molecular liquid and the mixed-atomic fluid at ∼400  GPa, i.e., in a much more interesting pressure range. We provide numerical evidence supporting the metallic behavior of this intermediate phase. Therefore, we predict that the metallization at finite temperature occurs in this partially dissociated molecular fluid, well before the complete atomization of the liquid. At high temperature this first-order transition becomes a crossover, in very good agreement with the experimental observation. Several systematic tests supporting the quality of our large scale calculations are also reported.

Ab Initio Geometry and Bright Excitation of Carotenoids: Quantum Monte Carlo and Many Body Green's Function Theory Calculations on Peridinin

In this letter, we report the singlet ground state structure of the full carotenoid peridinin by means of variational Monte Carlo (VMC) calculations. The VMC relaxed geometry has an average bond length alternation of 0.1165(10) Å, larger than the values obtained by DFT (PBE, B3LYP, and CAM-B3LYP) and shorter than that calculated at the Hartree-Fock (HF) level. TDDFT and EOM-CCSD calculations on a reduced peridinin model confirm the HOMO-LUMO major contribution of the Bu(+)-like (S2) bright excited state. Many Body Green's Function Theory (MBGFT) calculations of the vertical excitation energy of the Bu(+)-like state for the VMC structure (VMC/MBGFT) provide an excitation energy of 2.62 eV, in agreement with experimental results in n-hexane (2.72 eV). The dependence of the excitation energy on the bond length alternation in the MBGFT and TDDFT calculations with different functionals is discussed.

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

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.

Reaction pathways by quantum Monte Carlo: insight on the torsion barrier of 1,3-butadiene, and the conrotatory ring opening of cyclobutene

Quantum Monte Carlo (QMC) methods are used to investigate the intramolecular reaction pathways of 1,3-butadiene. The ground state geometries of the three conformers s-trans, s-cis, and gauche, as well as the cyclobutene structure are fully optimised at the variational Monte Carlo (VMC) level, obtaining an excellent agreement with the experimental results and other quantum chemistry high level calculations. Transition state geometries are also estimated at the VMC level for the s-trans to gauche torsion barrier of 1,3-butadiene and for the conrotatory ring opening of cyclobutene to the gauche-1,3-butadiene conformer. The energies of the conformers and the reaction barriers are calculated at both variational and diffusional Monte Carlo levels providing a precise picture of the potential energy surface of 1,3-butadiene and supporting one of the two model profiles recently obtained by Raman spectroscopy [Boopalachandran et al., J. Phys. Chem. A 115, 8920 (2011)]. Considering the good scaling of QMC techniques with the system's size, our results also demonstrate how variational Monte Carlo calculations can be applied in the future to properly investigate the reaction pathways of large and correlated molecular systems.