BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sandro Sorella

Sandro Sorella appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elphysics.chem-phphysics.comp-phcond-mat.mtrl-scicond-mat.supr-conquant-phcond-mat.otherhep-latcond-mat.mes-hallcond-mat.quant-gasphysics.atom-ph

Catalog footprint

What is connected

42works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

The phase diagram of the Hubbard model by Variational Auxiliary Field quantum Monte Carlo

A systematically improvable wave function is proposed for the numerical solution of strongly correlated systems. With a stochastic optimization method, based on the auxiliary field quantum Monte Carlo technique, an effective temperature Teff is defined, probing the distance of the ground state properties of the model in the thermodynamic limit from the ones of the proposed correlated mean-field ansatz. In this way their uncertainties from the unbiased zero temperature limit may be estimated by simple and stable extrapolations well before the so called sign problem gets prohibitive. At finite Teff the convergence of the energy to the thermodynamic limit is indeed shown to be possible in the Hubbard model already for relatively small square lattices with linear dimension L ~10, thanks to appropriate averages over several twisted boundary conditions. Within the estimated energy accuracy of the proposed variational ansatz, two clear phases are identified, as the energy is lowered by spontaneously breaking some symmetries satisfied by the Hubbard Hamiltonian: a) a stripe phase where both spin and translation symmetries are broken, and b) a strong coupling d-wave superconducting phase when the particle number is not conserved and global U(1) symmetry is broken. On the other hand the symmetric phase is stable in a wide region at large doping and small coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

QMC study of the chiral Heisenberg Gross-Neveu universality class

We investigate a quantum criticality of an antiferromagnetic phase transition in the Hubbard model on a square lattice with a $d$-wave pairing field by large-scale auxiliary-field quantum Monte Carlo simulations. Since the $d$-wave pairing filed induces Dirac cones in the non-interacting single-particle spectrum, the quantum criticality should correspond to the chiral Heisenberg universality class in terms of the Gross-Neveu theory, which is the same as those expected in the Hubbard model on the honeycomb lattice, despite the unit cells being different (e.g., they contain one and two sites, respectively). We show that both the two phase transitions, expected to occur on the square and on the honeycomb lattices, indeed have the same quantum criticality. We also argue that details of the models, i.e., the way of counting the total number $N$ of fermion components and the anisotropy of the Dirac cones, do not change the critical exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Space-warp coordinate transformation for efficient ionic force calculations in quantum Monte Carlo

Ab-initio quantum Monte Carlo (QMC) methods are a state-of-the-art computational approach to obtaining highly accurate many-body wave functions. Although QMC methods are widely used in physics and chemistry to compute ground-state energies, calculation of atomic forces is still under technical/algorithmic development. Very recently, force evaluation has started to become of paramount importance for the generation of machine-learning force-field potentials. Nevertheless, there is no consensus regarding whether an efficient algorithm is available for the QMC force evaluation, namely one that scales well with the number of electrons and the atomic numbers. In this study, we benchmark the accuracy of all-electron variational Monte Carlo (VMC) and lattice-regularized diffusion Monte Carlo (LRDMC) forces for various mono- and heteronuclear dimers. The VMC and LRDMC forces were calculated with and without the so-called space-warp coordinate transformation (SWCT) and appropriate regularization techniques to remove the infinite variance problem. The LRDMC forces were computed with the Reynolds (RE) and variational-drift (VD) approximations. The potential energy surfaces obtained from the LRDMC energies give equilibrium bond lengths ($r_{\rm eq}$) and harmonic frequencies ($ω$) very close to the experimental values for all dimers, improving the corresponding VMC results. The LRDMC forces improve the VMC forces, implying that it is worth computing the DMC forces beyond VMC in spite of the higher computational cost. We find that the ratio of computational costs between QMC energy and forces scales as $Z^{\sim 2.5}$ without the SWCT. In contrast, the application of the SWCT makes the ratio {\it independent} of $Z$. As such, the accessible QMC system size is not affected by the evaluation of ionic forces but governed by the same scaling as the total energy one.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The two-dimensional $t$-$t^{\prime}$ Holstein model

The competition and interplay between charge-density wave and superconductivity have become a central subject for quasi-2D compounds. Some of these materials, such as the transition-metal dichalcogenides, exhibit strong electron-phonon coupling, an interaction that may favor both phases, depending on the external parameters, such as hydrostatic pressure. In view of this, here we analyze the single-band $t$-$t^{\prime}$ Holstein model in the square lattice, adding a next-nearest neighbor hopping $t^{\prime}$ in order to play the role of the external pressure. To this end, we perform unbiased quantum Monte Carlo simulations with an efficient inversion sampling technique appropriately devised for this model. Such a methodology drastically reduces the autocorrelation time, and increases the efficiency of the Monte Carlo approach. By investigating the charge-charge correlation functions, we obtain the behavior of the critical temperature as a function of $t^{\prime}$, and from compressibility analysis, we show that a first-order metal-to-insulator phase transition occurs. We also provide a low-temperature phase diagram for the model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

General correlated geminal ansatz for electronic structure calculations: exploiting Pfaffians in place of determinants

We propose here a single Pfaffian correlated variational ansatz, that dramatically improves the accuracy with respect to the single determinant one, while remaining at a similar computational cost. A much larger correlation energy is indeed determined by the most general two electron pairing function, including both singlet and triplet channels, combined with a many-body Jastrow factor, including all possible spin-spin spin-density and density-density terms. The main technical ingredient to exploit this accuracy is the use of the Pfaffian for antisymmetrizing an highly correlated pairing function, thus recovering the Fermi statistics for electrons with an affordable computational cost. Moreover the application of the Diffusion Monte Carlo, within the fixed node approximation, allows us to obtain very accurate binding energies for the first preliminary calculations reported in this study: C$_2$, N$_2$ and O$_2$ and the benzene molecule. This is promising and remarkable, considering that they represent extremely difficult molecules even for computationally demanding multi-determinant approaches, and opens therefore the way for realistic and accurate electronic simulations with an algorithm scaling at most as the fourth power of the number of electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ground-state properties of the hydrogen chain: insulator-to-metal transition, dimerization, and magnetic phases

Accurate and predictive computations of the quantum-mechanical behavior of many interacting electrons in realistic atomic environments are critical for the theoretical design of materials with desired properties, and require solving the grand-challenge problem of the many-electron Schrodinger equation. An infinite chain of equispaced hydrogen atoms is perhaps the simplest realistic model for a bulk material, embodying several central themes of modern condensed matter physics and chemistry, while retaining a connection to the paradigmatic Hubbard model. Here we report a combined application of cutting-edge computational methods to determine the properties of the hydrogen chain in its quantum-mechanical ground state. Varying the separation between the nuclei leads to a rich phase diagram, including a Mott phase with quasi long-range antiferromagnetic order, electron density dimerization with power-law correlations, an insulator-to-metal transition and an intricate set of intertwined magnetic orders.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

TurboRVB: a many-body toolkit for {\it ab initio} electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for {\it ab initio} Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC), and Diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions. The electronic wave function (WF) is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal product with spin-singlet pairing, or as a Pfaffian, including both singlet and triplet correlations. This wave function can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) ansatz, first proposed by L. Pauling and P. W. Anderson in quantum chemistry and condensed matter physics, respectively. The RVB ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization is made possible thanks to state-of-the-art stochastic algorithms for energy minimization. The code has been efficiently parallelized by using a hybrid MPI-OpenMP protocol, that is also an ideal environment for exploiting the computational power of modern GPU accelerators.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

All-electron quantum Monte Carlo with Jastrow single determinant Ansatz: application to the sodium dimer

In this work, we report potential energy surfaces (PESs) of the sodium dimer calculated by variational (VMC) and lattice regularized diffusion Monte Carlo (LRDMC). The VMC calculation is accurate for determining the equilibrium distance and the qualitative shape of the experimental PES. Remarkably, after the application of the LRDMC projection to this single determinant ansatz, namely the Jastrow Antisymmetrized Geminal Product (JAGP), chemical accuracy (~ 1kcal/mol) is reached, and the obtained dissociation energy, equilibrium internuclear distance, and harmonic vibrational frequency are in very good agreement with the experimental ones. This outcome crucially depends on the quality of the optimization used to determine the best possible trial function within the chosen ansatz. The strategy adopted in this work is to minimize the variational energy by initializing the trial function with the DFT single determinant ansatz expanded exactly in the same atomic basis used for the corresponding VMC and LRDMC calculations. This atomic basis is ad-hoc reshaped for QMC calculations. Indeed, we multiply the standard Gaussian type atomic orbitals by a one-body Jastrow factor, satisfying in this way the electron-ion cusp conditions. This allows us to use a very small basis almost converged in the complete basis set limit, by reducing the computational effort as well as the statistical fluctuations on the total energy. In order to achieve these important advantages, we have defined a very efficient DFT algorithm in the mentioned basis, by estimating the corresponding matrix elements on a mesh, and by using a much finer mesh grid in the vicinity of nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Speeding up the ab initio diffusion Monte Carlo by a smart lattice regularization

One of the most significant drawbacks of the all-electron ab initio diffusion Monte Carlo (DMC) is that its computational cost drastically increases with the atomic number ($Z$), which typically scales with $Z^{\sim 6}$. In this study, we introduce an algorithm based on a very efficient implementation of the Lattice Regularized Diffusion Monte Carlo (LRDMC), where the conventional time discretization is replaced by its lattice space counterpart. This scheme enables us to conveniently adopt a small lattice space in the vicinity of nuclei, and a large one in the valence region, by which a considerable speedup is achieved, especially for large atomic number $Z$. Indeed, the computational performances of our algorithm can be theoretically established by using the Thomas-Fermi model for heavy atoms, yielding an almost affordable scaling with the atomic number, i.e., $Z^{\sim 5}$. This opens the way for efficient and accurate all-electron ab initio DMC in electronic structure calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Assessing the orbital selective Mott transition with variational wave functions

We study the Mott metal-insulator transition in the two-band Hubbard model with different hopping amplitudes $t_1$ and $t_2$ for the two orbitals on the two-dimensional square lattice by using {\it non-magnetic} variational wave functions, similarly to what has been considered in the limit of infinite dimensions by dynamical mean-field theory. We work out the phase diagram at half filling (i.e., two electrons per site) as a function of $R=t_2/t_1$ and the on-site Coulomb repulsion $U$, for two values of the Hund's coupling $J=0$ and $J/U=0.1$. Our results are in good agreement with previous dynamical mean-field theory calculations, demonstrating that the non-magnetic phase diagram is only slightly modified from infinite to two spatial dimensions. Three phases are present: a metallic one, for small values of $U$, where both orbitals are itinerant; a Mott insulator, for large values of $U$, where both orbitals are localized because of the Coulomb repulsion; and the so-called orbital-selective Mott insulator (OSMI), for small values of $R$ and intermediate $U$'s, where one orbital is localized while the other one is still itinerant. The effect of the Hund's coupling is two-fold: on one side, it favors the full Mott phase over the OSMI; on the other side, it stabilizes the OSMI at larger values of $R$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Boosting the accuracy and speed of quantum Monte Carlo: size-consistency and time-step

Diffusion Monte Carlo (DMC) simulations for fermions are becoming the standard to provide high quality reference data in systems that are too large to be investigated via quantum chemical approaches. DMC with the fixed-node approximation relies on modifications of the Green function to avoid singularities near the nodal surface of the trial wavefunction. We show that these modifications affect the DMC energies in a way that is not size-consistent, resulting in large time-step errors. Building on the modifications of Umrigar {\em et al.} and of DePasquale {\em et al.} we propose a simple Green function modification that restores size-consistency to large values of time-step; substantially reducing the time-step errors. The new algorithm also yields remarkable speedups of up to two orders of magnitude in the calculation of molecule-molecule binding energies and crystal cohesive energies, thus extending the horizons of what is possible with DMC.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Competing collinear magnetic structures in superconducting FeSe by first principles quantum Monte Carlo calculations

Resolving the interplay between magnetic interactions and structural properties in strongly correlated materials through a quantitatively accurate approach has been a major challenge in condensed matter physics. Here we apply highly accurate first principles quantum Monte Carlo (QMC) techniques to obtain structural and magnetic properties of the iron selenide (FeSe) superconductor under pressure. Where comparable, the computed properties are very close to the experimental values. Of potential ordered magnetic configurations, collinear spin configurations are the most energetically favorable over the explored pressure range. They become nearly degenerate in energy with bicollinear spin orderings at around 7 GPa, when the experimental critical temperature $T_c$ is the highest. On the other hand, ferromagnetic, checkerboard, and staggered dimer configurations become relatively higher in energy as the pressure increases. The behavior under pressure is explained by an accurate analysis of the charge compressibility and the orbital occupation as described by the QMC many-body wave function, which reveals how spin, charge and orbital degrees of freedom are strongly coupled in this compound. This remarkable pressure evolution suggests that stripe-like magnetic fluctuations may be responsible for the enhanced $T_c$ in FeSe and that higher T$_c$ is associated with nearness to a crossover between collinear and bicollinear ordering.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hidden Mott transition and large-$U$ superconductivity in the two-dimensional Hubbard model

We consider the one-band Hubbard model on the square lattice by using variational and Green's function Monte Carlo methods, where the variational states contain Jastrow and backflow correlations on top of an uncorrelated wave function that includes BCS pairing and magnetic order. At half filling, where the ground state is antiferromagnetically ordered for any value of the on-site interaction $U$, we can identify a hidden critical point $U_{\rm Mott}$, above which a finite BCS pairing is stabilized in the wave function. The existence of this point is reminiscent of the Mott transition in the paramagnetic sector and determines a separation between a Slater insulator (at small values of $U$), where magnetism induces a potential energy gain, and a Mott insulator (at large values of $U$), where magnetic correlations drive a kinetic energy gain. Most importantly, the existence of $U_{\rm Mott}$ has crucial consequences when doping the system: We observe a tendency to phase separation into a hole-rich and a hole-poor region only when doping the Slater insulator, while the system is uniform by doping the Mott insulator. Superconducting correlations are clearly observed above $U_{\rm Mott}$, leading to the characteristic dome structure in doping. Furthermore, we show that the energy gain due to the presence of a finite BCS pairing above $U_{\rm Mott}$ shifts from the potential to the kinetic sector by increasing the value of the Coulomb repulsion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Toward Accurate Adsorption Energetics on Clay Surfaces

Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universal quantum criticality in the metal-insulator transition of two-dimensional interacting Dirac electrons

The metal-insulator transition has been a subject of intense research since Nevil Mott has first proposed that the metallic behavior of interacting electrons could turn to the insulating one as electron correlations increase. Here, we consider electrons with massless Dirac-like dispersion in two spatial dimensions, described by the Hubbard models on two geometrically different lattices, and perform numerically exact calculations on unprecedentedly large systems that, combined with a careful finite size scaling analysis, allow us to explore the quantum critical behavior in the vicinity of the interaction-driven metal-insulator transition. We find thereby that the transition is continuous and determine the quantum criticality for the corresponding universality class, which is described in the continuous limit by the Gross-Neveu model, a model extensively studied in quantum field theory. We furthermore discuss a fluctuation-driven scenario for the metal-insulator transition in the interacting Dirac electrons: the metal-insulator transition is triggered only by the vanishing of the quasiparticle weight but not the Dirac Fermi velocity, which instead remains finite near the transition. This important feature cannot be captured by a simple mean-field or Gutzwiller-type approximate picture, but is rather consistent with the low energy behavior of the Gross-Neveu model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite size scaling with modified boundary conditions

An efficient scheme is introduced for a fast and smooth convergence to the thermodynamic limit with finite size cluster calculations. This is obtained by modifying the energy levels of the non interacting Hamiltonian in a way consistent with the corresponding one particle density of states in the thermodynamic limit. After this modification exact free electron energies are obtained with finite size calculations and for particular fillings that satisfy the so called "closed shell condition". In this case the "sign problem" is particularly mild in the auxiliary field quantum Monte Carlo technique and therefore, with this technique, it is possible to obtain converged energies for the Hubbard model even for $U>0$. We provide a strong numerical evidence that phase separation occurs in the low doping region and moderate $U\lesssim 4t$ regime of this model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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. Recently it has been reported a complete atomization transition between a mixed-atomic liquid and a completely dissociated fluid in an almost unaccessible pressure range {[\emph{Nat. Commun.} {\bf 5}, 3487 (2014)]}. Here instead, in a much more interesting pressure range, we identify a different transition between the fully molecular liquid and the mixed-atomic fluid at $\sim$ 400 GPa, with numerical evidence supporting its metallic behavior. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lanczos steps to improve variational wave functions

Gutzwiller-projected fermionic states can be efficiently implemented within quantum Monte Carlo calculations to define extremely accurate variational wave functions for Heisenberg models on frustrated two-dimensional lattices, not only for the ground state but also for low-energy excitations. The application of few Lanczos steps on top of these states further improves their accuracy, allowing calculations on large clusters. In addition, by computing both the energy and its variance, it is possible to obtain reliable estimations of exact results. Here, we report the cases of the frustrated Heisenberg models on square and Kagome lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Light-cone effect and supersonic correlations in one- and two-dimensional bosonic superfluids

We study the spreading of density-density correlations in Bose-Hubbard models after a quench of the interaction strength, using time-dependent variational Monte Carlo simulations. It gives access to unprecedented long propagation times and to dimensions higher than one. In both one and two dimensions, we find ballistic light-cone spreading of correlations and extract accurate values of the light-cone velocity in the superfluid regime. We show that the spreading of correlations is generally supersonic, with a light-cone propagating faster than sound modes but slower than the maximum group velocity of density excitations, except at the Mott transition, where all the characteristic velocities are equal. Further, we show that in two dimensions the correlation spreading is highly anisotropic and presents nontrivial interference effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Static and dynamical correlation in diradical molecules by Quantum Monte Carlo using the Jastrow Antisymmetrized Geminal Power ansatz

Diradical molecules are essential species involved in many organic and inorganic chemical reactions. The computational study of their electronic structure is often challenging, because a reliable description of the correlation, and in particular of the static one, requires multi-reference techniques. The Jastrow correlated Antisymmetrized Geminal Power (JAGP) is a compact and efficient wave function ansatz, based on the valence-bond representation, which can be used within Quantum Monte Carlo (QMC) approaches. The AGP part can be rewritten in terms of molecular orbitals, obtaining a multi-determinant expansion with zero-seniority number. In the present work we demonstrate the capability of the JAGP ansatz to correctly describe the electronic structure of two diradical prototypes: the orthogonally twisted ethylene, C2H4, and the methylene, CH2, representing respectively a homosymmetric and heterosymmetric system. On the other hand, we show that the simple ansatz of a Jastrow correlated Single Determinant (JSD) wave function is unable to provide an accurate description of the electronic structure in these diradical molecules, both at variational level and, more remarkably, in the fixed-nodes projection schemes showing that a poor description of the static correlation yields an inaccurate nodal surface. The suitability of JAGP to correctly describe diradicals with a computational cost comparable with that of a JSD calculation, in combination with a favorable scalability of QMC algorithms with the system size, opens new perspectives in the ab initio study of large diradical systems, like the transition states in cycloaddition reactions and the thermal isomerization of biological chromophores.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unexpectedly high pressure for molecular dissociation in liquid hydrogen by a reliable 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 roperties. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Direct evidence for a gapless $Z_2$ spin liquid by frustrating Néel antiferromagnetism

By direct calculations of the spin gap in the frustrated Heisenberg model on the square lattice, with nearest- ($J_1$) and next-nearest-neighbor ($J_2$) super-exchange couplings, we provide a solid evidence that the spin-liquid phase in the frustrated regime $0.45 \lesssim J_2/J_1 \lesssim 0.6$ is gapless. Our numerical method is based on a variational wave function that is {\it systematically} improved by the application of few Lanczos steps and allows us to obtain reliable extrapolations in the thermodynamic limit. The peculiar nature of the non-magnetic state is unveiled by the existence of S=1 gapless excitations at $k=(π,0)$ and $(0,π)$. The magnetic transition can be described and interpreted by a variational state that is built from Abrikosov fermions having a $Z_2$ gauge structure and four Dirac points in the spinon spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gapless spin-liquid phase in the kagome spin-1/2 Heisenberg antiferromagnet

We study the energy and the static spin structure factor of the ground state of the spin-1/2 quantum Heisenberg antiferromagnetic model on the kagome lattice. By the iterative application of a few Lanczos steps on accurate projected fermionic wave functions and the Green's function Monte Carlo technique, we find that a gapless (algebraic) U(1) Dirac spin liquid is competitive with previously proposed gapped (topological) Z2 spin liquids. By performing a finite-size extrapolation of the ground-state energy, we obtain an energy per site E/J=-0.4365(2), which is equal, within three error bars, to the estimates given by the density-matrix renormalization group (DMRG). Our estimate is obtained for a translationally invariant system, and, therefore, does not suffer from boundary effects, like in DMRG. Moreover, on finite toric clusters at the pure variational level, our energies are lower compared to those from DMRG calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Improper s-wave symmetry for the electronic pairing in iron-based superconductors by first-principles calculation

By means of space-group symmetry arguments, we argue that the electronic pairing in iron-based high temperature superconductors shows a structure which is a linear combination of planar s-wave and d-wave symmetry channels, both preserving the 3-dimensional A_1g irreducible representation of the corresponding crystal point-group. We demonstrate that the s- and d-wave channels are determined by the parity under reflection of the electronic orbitals through the iron planes, and by improper rotations around the iron sites. We provide evidence of these general properties by performing accurate quantum Monte Carlo ab-initio calculations of the pairing function, for a FeSe lattice with tetragonal experimental geometry at ambient pressure. We find that this picture survives even in the FeSe under pressure and at low temperatures, when the tetragonal point-group symmetry is slightly broken. In order to achieve a higher resolution in momentum space we introduce a BCS model that faithfully describes our QMC variational pairing function on the simulated 4x4 FeSe unit cell. This allows us to provide a k-resolved image of the pairing function, and show that non-isotropic contributions in the BCS gap function are related to the improper s-wave symmetry. Our theory can rationalize and explain a series of contradictory experimental findings, such as the observation of twofold symmetry in the FeSe superconducting phase, the anomalous drop of T_c with Co-impurity in LaFeAsO_(1-x)F_x, the s-to-d-wave gap transition in BaFe_2As_2 under K doping, and the nodes appearing in the LiFeAs superconducting gap upon P isovalent substitution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum Monte Carlo study of the protonated water dimer

We report an extensive theoretical study of the protonated water dimer (Zundel ion) by means of the highly correlated variational Monte Carlo and lattice regularized Monte Carlo approaches. This system represents the simplest model for proton transfer (PT) and a correct description of its properties is essential in order to understand the PT mechanism in more complex acqueous systems. Our Jastrow correlated AGP wave function ensures an accurate treatment of electron correlations. Exploiting the advantages of contracting the primitive basis set over atomic hybrid orbitals, we are able to limit dramatically the number of variational parameters with a systematic control on the numerical precision, crucial in order to simulate larger systems. We investigate energetics and geometrical properties of the Zundel ion as a function of the oxygen-oxygen distance, taken as reaction coordinate. In both cases, our QMC results are found in excellent agreement with coupled cluster CCSD(T) technique, the quantum chemistry "gold standard". Calculations on proton transfer static barriers and dissociation energies display the same agreement. A comparison with density functional theory (DFT) results in the PBE approximation points out the crucial role of electron correlations for a correct description of the PT in the dimer. Our method is able to resolve the tiny energy differences (~ 0.1 Kcal/mol) and the corresponding structural variations at which the proton transfer takes place; it combines these features with a N^3-N^4 scaling with number of particles (favorable with respect to other post-DFT methods), as proven by a benchmark calculation on a larger protonated water cluster. The QMC approach used here represents a promising candidate to provide the first high-level ab initio description of PT in water.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Absence of a Spin Liquid Phase in the Hubbard Model on the Honeycomb Lattice

A spin liquid is a novel quantum state of matter with no conventional order parameter where a finite charge gap exists even though the band theory would predict metallic behavior. Finding a stable spin liquid in two or higher spatial dimensions is one of the most challenging and debated issues in condensed matter physics. Very recently, it has been reported that a model of graphene, i.e., the Hubbard model on the honeycomb lattice, can show a spin liquid ground state in a wide region of the phase diagram, between a semi-metal (SM) and an antiferromagnetic insulator (AFMI). Here, by performing numerically exact quantum Monte Carlo simulations, we extend the previous study to much larger clusters (containing up to 2592 sites), and find, if any, a very weak evidence of this spin liquid region. Instead, our calculations strongly indicate a direct and continuous quantum phase transition between SM and AFMI.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Absence of static stripes in the two-dimensional $t{-}J$ model by an accurate and systematic quantum Monte Carlo approach

We examine the two-dimensional $t{-}J$ model by using variational approach combined with well established quantum Monte Carlo techniques [S. Sorella {\it et al.}, \prl {\bf 88}, 117002 (2002)] that are used to improve systematically the accuracy of the variational ansatz. Contrary to recent density-matrix renormalization group and projected entangled-pair state calculations [P. Corboz {\it et al.}, \prb {\bf 84}, 041108(R) (2011)], a uniform phase is found for $J/t=0.4$, even when the calculation is biased with an ansatz that explicitly contains stripe order. Moreover, in the small hole doping regime, i.e., $δ\lesssim 0.1$, our results support the coexistence of antiferromagnetism and superconductivity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite temperature electronic simulations beyond the Born-Oppenheimer approximation

We introduce a general technique to compute finite temperature electronic properties by a novel covariant formulation of the electronic partition function. By using a rigorous variational upper bound to the free energy we are led to the evaluation of a partition function that can be computed stochastically by sampling electronic wave functions and atomic positions (assumed classical). In order to achieve this target we show that it is extremely important to consider the non trivial geometry of the space defined by the wave function ansatz. The method can be extended to any technique capable to provide an energy value over a given wave function ansatz depending on several variational parameters and atomic positions. In particular we can take into account electronic correlation, by using the standard variational quantum Monte Carlo method, that has been so far limited to zero temperature ground state properties. We show that our approximation reduces correctly to the standard Born-Oppenheimer (BO) one at zero temperature and to the correct high temperature limit. At large enough temperatures this method allows to improve the BO, providing lower values of the electronic free energy, because within this method it is possible to take into account the electron entropy. We test this new method on the simple hydrogen molecule, where at low temperature we recover the correct BO low temperature limit. Moreover, we show that the dissociation of the molecule is possible at a temperature much smaller than the BO prediction. Several extension of the proposed technique are also discussed, as for instance the calculation of critical (magnetic, superconducting) temperatures, or transition rates in chemical reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Linearized Auxiliary fields Monte Carlo: efficient sampling of the fermion sign

We introduce a method that combines the power of both the lattice Green function Monte Carlo (LGFMC) with the auxiliary field techniques (AFQMC), and allows us to compute exact ground state properties of the Hubbard model for U<~ 4t on finite clusters. Thanks to LGFMC one obtains unbiased zero temperature results, not affected by the so called Trotter approximation of the imaginary time propagator exp(- H t). On the other hand the AFQMC formalism yields a remarkably fast convergence in t before the fermion sign problem becomes prohibitive. As a first application we report ground state energies in the Hubbard model at U/t=4 with up to one hundred sites.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strong electronic correlation in the Hydrogen chain: a variational Monte Carlo study

In this article, we report a fully ab initio variational Monte Carlo study of the linear, and periodic chain of Hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density matrix renormalization group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in details the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The fate of the resonating valence bond in graphene

We apply a variational wave function capable of describing qualitatively and quantitatively the so called "resonating valence bond" in realistic materials, by improving standard ab initio calculations by means of quantum Monte Carlo methods. In this framework we clearly identify the Kekulé and Dewar contributions to the chemical bond of the benzene molecule, and we establish the corresponding resonating valence bond energy of these well known structures ($\simeq 0.01$eV/atom). We apply this method to unveil the nature of the chemical bond in undoped graphene and show that this picture remains only within a small "resonance length" of few atomic units.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab-initio calculations for the beta-tin diamond transition in Silicon: comparing theories with experiments

We investigate the pressure-induced metal-insulator transition from diamond to beta-tin in bulk Silicon, using quantum Monte Carlo (QMC) and density functional theory (DFT) approaches. We show that it is possible to efficiently describe many-body effects, using a variational wave function with an optimized Jastrow factor and a Slater determinant. Variational results are obtained with a small computational cost and are further improved by performing diffusion Monte Carlo calculations and an explicit optimization of molecular orbitals in the determinant. Finite temperature corrections and zero point motion effects are included by calculating phonon dispersions in both phases at the DFT level. Our results indicate that the theoretical QMC (DFT) transition pressure is significantly larger (smaller) than the accepted experimental value. We discuss the limitation of DFT approaches due to the choice of the exchange and correlation functionals and the difficulty to determine consistent pseudopotentials within the QMC framework, a limitation that may significantly affect the accuracy of the technique.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Algorithmic differentiation and the calculation of forces by quantum Monte Carlo

We describe an efficient algorithm to compute forces in quantum Monte Carlo using adjoint algorithmic differentiation. This allows us to apply the space warp coordinate transformation in differential form, and compute all the 3M force components of a system with M atoms with a computational effort comparable with the one to obtain the total energy. Few examples illustrating the method for an electronic system containing several water molecules are presented. With the present technique, the calculation of finite-temperature thermodynamic properties of materials with quantum Monte Carlo will be feasible in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Size-consistent variational approaches to non-local pseudopotentials: standard and lattice regularized diffusion Monte Carlo methods revisited

We propose improved versions of the standard diffusion Monte Carlo (DMC) and the lattice regularized diffusion Monte Carlo (LRDMC) algorithms. For the DMC method, we refine a scheme recently devised to treat non-local pseudopotential in a variational way. We show that such scheme --when applied to large enough systems-- maintains its effectiveness only at correspondingly small enough time-steps, and we present two simple upgrades of the method which guarantee the variational property in a size-consistent manner. For the LRDMC method, which is size-consistent and variational by construction, we enhance the computational efficiency by introducing (i) an improved definition of the effective lattice Hamiltonian which remains size-consistent and entails a small lattice-space error with a known leading term, and (ii) a new randomization method for the positions of the lattice knots which requires a single lattice-space.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The nature and strength of inter-layer binding in graphite

We computed the inter-layer bonding properties of graphite using an ab-initio many body theory. We carried out variational and diffusion quantum Monte Carlo calculations and found an equilibrium inter-layer binding energy in good agreement with most recent experiments. We also analyzed the behavior of the total energy as a function of interlayer separation at large distances comparing the results with the predictions of the random phase approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Metal-insulator transition and strong-coupling spin liquid in the $t{-}t^\prime$ Hubbard model

We study the phase diagram of the frustrated $t{-}t^\prime$ Hubbard model on the square lattice by using a novel variational wave function. Taking the clue from the backflow correlations that have been introduced long-time ago by Feynman and Cohen and have been used for describing various interacting systems on the continuum (like liquid $^3$He, the electron jellium, and metallic Hydrogen), we consider many-body correlations to construct a suitable approximation for the ground state of this correlated model on the lattice. In this way, a very accurate {\it ansatz} can be achieved both at weak and strong coupling. We present the evidence that an insulating and non-magnetic phase can be stabilized at strong coupling and sufficiently large frustrating ratio $t^\prime/t$.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Ground state properties of the one dimensional Coulomb gas

We study the ground state properties of a quasi one dimensional electron gas, interacting via an effective potential with a harmonic transversal confinement and long range Coulomb tail. The exact correlation energy has been calculated for a wide range of electron densities by using the lattice regularized diffusion Monte Carlo method, which is a recent development of the standard projection Monte Carlo technique. In this case it is particularly useful as it allows to sample the exact ground state of the system, even in the low density regime when the exchange between electrons is extremely small. For different values of the width parameter b (0.1 a*_0 <= b <= 4 a*_0), we give a simple parametrization of the correlation energy, which provides an accurate local density energy functional for quasi one dimensional systems. Moreover we show that static correlations are in qualitative agreement with those obtained for the Luttinger liquid model with long range interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Two spin liquid phases in the spatially anisotropic triangular Heisenberg model

The quantum spin-1/2 antiferromagnetic Heisenberg model on a two dimensional triangular lattice geometry with spatial anisotropy is relevant to describe materials like ${\rm Cs_2 Cu Cl_4}$ and organic compounds like {$κ$-(ET)$_2$Cu$_2$(CN)$_3$}. The strength of the spatial anisotropy can increase quantum fluctuations and can destabilize the magnetically ordered state leading to non conventional spin liquid phases. In order to understand these intriguing phenomena, quantum Monte Carlo methods are used to study this model system as a function of the anisotropic strength, represented by the ratio $J'/J$ between the intra-chain nearest neighbor coupling $J$ and the inter-chain one $J'$. We have found evidence of two spin liquid regions. The first one is stable for small values of the coupling $J'/J \alt 0.65$, and appears gapless and fractionalized, whereas the second one is a more conventional spin liquid with a small spin gap and is energetically favored in the region $0.65\alt J'/J \alt 0.8$. We have also shown that in both spin liquid phases there is no evidence of broken translation symmetry with dimer or spin-Peirls order or any broken spatial reflection symmetry of the lattice. The various phases are in good agreement with the experimental findings, thus supporting the existence of spin liquid phases in two dimensional quantum spin-1/2 systems.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Stripes and spin-incommensurabilities are favored by lattice anisotropies

Structural distortions in cuprate materials give a natural origin for anisotropies in electron properties. We study a modified one-band t-J model in which we allow for different hoppings and antiferromagnetic couplings in the two spatial directions ($t_x \ne t_y$ and $J_x \ne J_y$). Incommensurate peaks in the spin structure factor show up only in the presence of a lattice anisotropy, whereas charge correlations, indicating enhanced fluctuations at incommensurate wave vectors, are almost unaffected with respect to the isotropic case.

0 citations0 reviewsNo code linkedOpen access