BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Sandro Sorella

Sandro Sorella contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elphysics.chem-phphysics.comp-phcond-mat.mtrl-scicond-mat.otherhep-latquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 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