Guglielmo Mazzola

Unexpectedly high pressure for molecular dissociation in liquid hydrogen by electronic simulation.

The study of the high pressure phase diagram of hydrogen has continued with renewed effort for about one century as it remains a fundamental challenge for experimental and theoretical techniques. Here we employ an efficient molecular dynamics based on the quantum Monte Carlo method, which can describe accurately the electronic correlation and treat a large number of hydrogen atoms, allowing a realistic and reliable prediction of thermodynamic properties. We find that the molecular liquid phase is unexpectedly stable, and the transition towards a fully atomic liquid phase occurs at much higher pressure than previously believed. The old standing problem of low-temperature atomization is, therefore, still far from experimental reach.

Finite-temperature electronic simulations without the Born-Oppenheimer constraint

The adiabatic approximation, typically assumed when performing standard Born-Oppenheimer (BO) molecular dynamics, can become unreliable at finite temperature, and specifically when the temperature is larger than the electronic energy gap between the ground state and the low-lying excited states. In this regime, relevant for many important chemical processes, the non-adiabatic couplings between the electronic energy states can produce finite temperature effects in several molecular properties, such as the geometry, the vibrational frequencies, the binding energy, and several chemical reactions. In this work, we introduce a novel finite-temperature non-adiabatic molecular dynamics based on a novel covariant formulation of the electronic partition function. In this framework, the nuclei are not constrained to move in a specific electronic potential energy surface. Then, by using a rigorous variational upper bound to the free energy, we are led to an approximate partition function that can be evaluated numerically. The method can be applied to any technique capable to provide an energy value over a given wave function ansatz depending on several variational parameters and atomic positions. In this work, we have applied the proposed method within a quantum Monte Carlo (QMC) scheme. In particular, we consider in this first application only classical ions, but we explicitly include an electronic correlation (Jastrow) term in the wave function, by extending in this way the standard variational QMC method, from ground state to finite temperature properties. We show that our approximation reduces correctly to the standard ground-state Born-Oppenheimer (gsBO) at zero temperature and to the correct high temperature limit. Moreover, at temperatures large enough, this method improves the upper bound of the free energy obtained with a single BO energy surface, since within our approach it is possible to estimate the electron entropy of a correlated ansatz in an efficient way. We test this new method on the simple hydrogen molecule, where at low temperature we recover the correct gsBO low temperature limit. Moreover, we show that the dissociation of the molecule is possible at a temperature much smaller than the one corresponding to the gsBO energy surface, in good agreement with experimental evidence. Several extensions of the proposed technique are also discussed, as for instance the inclusion of quantum effects for ions and the calculation of critical (magnetic, superconducting) temperatures.

Geminal embedding scheme for optimal atomic basis set construction in correlated calculations

We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.

Quantum Monte Carlo tunneling from quantum chemistry to quantum annealing

Quantum Tunneling is ubiquitous across different fields, from quantum chemical reactions, and magnetic materials to quantum simulators and quantum computers. While simulating the real-time quantum dynamics of tunneling is infeasible for high-dimensional systems, quantum tunneling also shows up in quantum Monte Carlo (QMC) simulations that scale polynomially with system size. Here we extend a recent results obtained for quantum spin models {[{Phys. Rev. Lett.} {\bf 117}, 180402 (2016)]}, and study high-dimensional continuos variable models for proton transfer reactions. We demonstrate that QMC simulations efficiently recover ground state tunneling rates due to the existence of an instanton path, which always connects the reactant state with the product. We discuss the implications of our results in the context of quantum chemical reactions and quantum annealing, where quantum tunneling is expected to be a valuable resource for solving combinatorial optimization problems.

Quantum Monte Carlo annealing with multi-spin dynamics

We introduce a novel Simulated Quantum Annealing (SQA) algorithm which employs a multispin quantum fluctuation operator. At variance with the usual transverse field, short-range two-spin flip interactions are included in the driver Hamiltonian. A Quantum Monte Carlo algorithm, capable of efficiently simulating large disordered systems, is described and tested. A first application to SQA, on a random square lattice Ising spin glass reveals that the multi-spin driver Hamiltonian improves upon the usual transverse field. This work paves the way for more systematic investigations using multi-spin quantum fluctuations on a broader range of problems.