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

Miguel A. Morales

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

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phcond-mat.str-elphysics.comp-phcond-mat.othercond-mat.stat-mechcond-mat.mes-hallphysics.geo-ph

Catalog footprint

What is connected

18works
8topics
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

18 published item(s)

preprint2025arXiv

Critical gate distance for Wigner crystallization in the two-dimensional electron gas

We report on the properties of the two-dimensional electron gas in a dual-gate geometry, using quantum Monte Carlo methods to obtain aspects of the phase diagram as a function of electron density and gate distance. We identify the critical gate distance below which the Wigner crystal phase disappears. For larger gate distances, the system undergoes a re-entrant transition from crystal to liquid at sufficiently low density. We also present preliminary evidence for a fully polarized ferromagnetic liquid state at low electron density and intermediate gate distances. The quantum Monte Carlo results are compared with simpler approximate methods, which are shown to be semi-quantitatively reliable for determining key features of the phase diagram. These methods are then used to obtain the phase boundary between the Wigner crystal and liquid in the single-gate geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Phaseless Auxiliary-Field Quantum Monte Carlo Perspective on the Uniform Electron Gas at Finite Temperatures: Issues, Observations, and Benchmark Study

We investigate the viability of the phaseless finite temperature auxiliary field quantum Monte Carlo (ph-FT-AFQMC) method for ab initio systems using the uniform electron gas as a model. Through comparisons with exact results and finite temperature coupled cluster theory, we find that ph-FT-AFQMC is sufficiently accurate at high to intermediate electronic densities. We show both analytically and numerically that the phaseless constraint at finite temperature is fundamentally different from its zero temperature counterpart (i.e., ph-ZT-AFQMC) and generally one should not expect ph-FT-AFQMC to agree with ph-ZT-AFQMC in the low temperature limit. With an efficient implementation, we are able to compare exchange-correlation energies to existing results in the thermodynamic limit and find that existing parameterizations are highly accurate. In particular, we found that ph-FT-AFQMC exchange-correlation energies are in a better agreement with a known parametrization than is restricted path-integral Monte Carlo in the regime of $Θ\le0.5$ and $r_s \le 2$, which highlights the strength of ph-FT-AFQMC.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nature of the bonded-to-atomic transition in liquid silica to TPa pressures

First-principles calculations and analysis of the thermodynamic, structural, and electronic properties of liquid SiO$_2$ characterize the bonded-to-atomic transition at 0.1--1.6 TPa and 10$^4$--10$^5$ K (1--7 eV), the high-energy-density regime relevant to understanding planetary interiors. We find strong ionic bonds that become short-lived due to high kinetics during the transition, with sensitivity of the transition temperature to pressure, and our calculated Hugoniots agree with past experimental data. These results reconcile previous experimental and theoretical findings by clarifying the nature of the bond dissociation process in early Earth and "rocky" (oxide) constituents of large planets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerating Auxiliary-Field Quantum Monte Carlo Simulations of Solids with Graphical Processing Unit

We outline how auxiliary-field quantum Monte Carlo (AFQMC) can leverage graphical processing units (GPUs) to accelerate the simulation of solid state sytems. By exploiting conservation of crystal momentum in the one- and two-electron integrals we show how to efficiently formulate the algorithm to best utilize current GPU architectures. We provide a detailed description of different optimization strategies and profile our implementation relative to standard approaches, demonstrating a factor of 40 speed up over a CPU implementation. With this increase in computational power we demonstrate the ability of AFQMC to systematically converge solid state calculations with respect to basis set and system size by computing the cohesive energy of Carbon in the diamond structure to within 0.02 eV of the experimental result.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerating the Convergence of Auxiliary-Field Quantum Monte Carlo inSolids with Optimized Gaussian Basis Sets

We investigate the use of optimized correlation consistent gaussian basis sets for the study of insulating solids with auxiliary-field quantum Monte Carlo (AFQMC). The exponents of the basis set are optimized through the minimization of the second order Møller--Plesset perturbation theory (MP2) energy in a small unit cell of the solid. We compare against other alternative basis sets proposed in the literature, namely calculations in the Kohn--Sham basis and in the natural orbitals of an MP2 calculation. We find that our optimized basis sets accelerate the convergence of the AFQMC correlation energy compared to a Kohn--Sham basis, and offer similar convergence to MP2 natural orbitals at a fraction of the cost needed to generate them. We also suggest the use of an improved, method independent, MP2-based basis set correction that significantly reduces the required basis set sizes needed to converge the correlation energy. With these developments, we study the relative performance of these basis sets in LiH, Si and MgO, and determine that our optimized basis sets yield the most consistent results as a function of volume. Using these optimized basis sets, we systematically converge the AFQMC calculations to the complete basis set and thermodynamic limit and find excellent agreement with experiment for systems studied. Although we focus on AFQMC, our basis set generation procedure is independent of the subsequent correlated wavefunction method used.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C$_{36}$ Fullerene and an Iron Porphyrin Model Complex

We present three distinct examples where phaseless auxiliary-field Quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wavefunction with essential symmetry breaking. We first utilized essential time-reversal symmetry breaking with ph-AFQMC to compute the triplet-singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the single-triplet gap of a known biradicaloid, C$_{36}$. ph-AFQMC with spin-unrestricted Hartree-Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller-Plesset perturbation theory ($κ$-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C$_{36}$. Lastly, we applied ph-AFQMC to the computation of the quintet-triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, $κ$-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and $κ$-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single determinant trial wavefunction with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Auxiliary-field based trial wave functions in quantum Monte Carlo calculations

Quantum Monte Carlo (QMC) algorithms have long relied on Jastrow factors to incorporate dynamic correlation into trial wave functions. While Jastrow-type wave functions have been widely employed in real-space algorithms, they have seen limited use in second-quantized QMC methods, particularly in projection methods that involve a stochastic evolution of the wave function in imaginary time. Here we propose a scheme for generating Jastrow-type correlated trial wave functions for auxiliary-field QMC methods. The method is based on decoupling the two-body Jastrow into one-body projectors coupled to auxiliary fields, which then operate on a single determinant to produce a multi-determinant trial wave function. We demonstrate that intelligent sampling of the most significant determinants in this expansion can produce compact trial wave functions that reduce errors in the calculated energies. Our technique may be readily generalized to accommodate a wide range of two-body Jastrow factors and applied to a variety of model and chemical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice

The MB-pol many-body potential has recently emerged as an accurate molecular model for water simulations from the gas to the condensed phase. In this study, the accuracy of MB-pol is systematically assessed across the three phases of water through extensive comparisons with experimental data and high-level ab initio calculations. Individual many-body contributions to the interaction energies as well as vibrational spectra of water clusters calculated with MB-pol are in excellent agreement with reference data obtained at the coupled cluster level. Several structural, thermodynamic, and dynamical properties of the liquid phase at atmospheric pressure are investigated through classical molecular dynamics simulations as a function of temperature. The structural properties of the liquid phase are in nearly quantitative agreement with X-ray diffraction data available over the temperature range from 268 to 368 K. The analysis of other thermodynamic and dynamical quantities emphasizes the importance of explicitly including nuclear quantum effects in the simulations, especially at low temperature, for a physically correct description of the properties of liquid water. Furthermore, both densities and lattice energies of several ice phases are also correctly reproduced by MB-pol. Following a recent study of DFT models for water, a score is assigned to each computed property, which demonstrates the high and, in many respects, unprecedented accuracy of MB-pol in representing all three phases of water.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Theory of Finite Size Effects for Electronic Quantum Monte Carlo Calculations of Liquids and Solids

Concentrating on zero temperature Quantum Monte Carlo calculations of electronic systems, we give a general description of the theory of finite size extrapolations of energies to the thermodynamic limit based on one and two-body correlation functions. We introduce new effective procedures, such as using the potential and wavefunction split-up into long and short range functions to simplify the method and we discuss how to treat backflow wavefunctions. Then we explicitly test the accuracy of our method to correct finite size errors on example hydrogen and helium many-body systems and show that the finite size bias can be drastically reduced for even small systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Benchmarking Hydrogen-Helium Mixtures with QMC: Energetics, Pressures, and Forces

An accurate understanding of the phase diagram of dense hydrogen and helium mixtures is a crucial component in the construction of accurate models of Jupiter, Saturn, and Jovian extrasolar planets. Though DFT based first principles methods have the potential to provide the accuracy and computational efficiency required for this task, recent benchmarking in hydrogen has shown that achieving this accuracy requires a judicious choice of functional, and a quantification of the errors introduced. In this work, we present a quantum Monte Carlo based benchmarking study of a wide range of density functionals for use in hydrogen-helium mixtures at thermodynamic conditions relevant for Jovian planets. Not only do we continue our program of benchmarking energetics and pressures, but we deploy QMC based force estimators and use them to gain insights into how well the local liquid structure is captured by different density functionals. We find that TPSS, BLYP and vdW-DF are the most accurate functionals by most metrics, and that the enthalpy, energy, and pressure errors are very well behaved as a function of helium concentration. Beyond this, we highlight and analyze the major error trends and relative differences exhibited by the major classes of functionals, and estimate the magnitudes of these effects when possible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the representation of many-body interactions in water

Recent work has shown that the many-body expansion of the interaction energy can effectively be used to develop analytical representations of global potential energy surfaces (PESs) for water. In this study, the role of short- and long-range contri- butions at different orders is investigated by analyzing water potentials that treat the leading terms of the many-body expansion through implicit (i.e., TTM3-F and TTM4-F PESs) and explicit (i.e., WHBB and MB-pol PESs) representations. It is found that explicit short-range representations of 2-body and 3-body interactions along with a physically correct integration of short- and long-range contributions are necessary for an accurate representation of the water interactions from the gas to the condensed phase. Similarly, a complete many-body representation of the dipole moment surface is found to be crucial to reproducing the correct intensities of the infrared spectrum of liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Benchmark of Exchange-Correlation Functionals for High Pressure Hydrogen using Quantum Monte Carlo

The ab-initio phase diagram of dense hydrogen is very sensitive to errors in the treatment of electronic correlation. Recently, it has been shown that the choice of the density functional has a large effect on the predicted location of both the liquid-liquid phase transition and the solid insulator-to-metal transition in dense hydrogen. To identify the most accurate functional for dense hydrogen applications, we systematically benchmark some of the most commonly used functionals using Quantum Monte Carlo. By considering several measures of functional accuracy, we conclude that the van der Waals and hybrid functionals significantly out perform LDA and PBE. We support these conclusions by analyzing the impact of functional choice on structural optimization in the molecular solid, and on the location of the liquid-liquid phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Combined Influence of Nuclear Quantum Effects and van der Waals Interactions on the Structure of Ambient Water

Path-integral molecular dynamics simulations based on density functional theory employing exchange-correlation density functionals capable of treating nonlocal van der Waals (vdW) interactions self-consistently provide a remarkably accurate description of ambient water. Moreover, they suggest that water's structure may be impacted by a combined influence between nuclear quantum effects and vdW interactions. The latter strongly favor the formation of a high-density liquid, whereas the inclusion of the former mitigates this by decreasing the mean hydrogen-bond (H-bond) distance. Examining the structure of water reveals that while the major fraction of molecules do in fact exhibit the traditional picture of near-tetrahedral coordination, the liquid considerably softer than previously simulations have suggested, including a much lower proportion of molecules double-donating H-bonds as well as a much larger distribution of their angles.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computing the energy of a water molecule using MultiDeterminants: A simple, efficient algorithm

Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave-functions are critical to ascertaining new physics. One such wave function is the multiSlater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wavefunctions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally we implement this method and use it to compute the ground state energy of a water molecule.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Free energy methods in Coupled Electron Ion Monte Carlo

Recent progress in simulation methodologies and in computer power allow first principle simulations of condensed systems with Born-Oppenheimer electronic energies obtained by Quantum Monte Carlo methods. Computing free energies and therefore getting a quantitative determination of phase diagrams is one step more demanding in terms of computer resources. In this paper we derive a general relation to compute the free energy of an ab-initio model with Reptation Quantum Monte Carlo (RQMC) energies from the knowledge of the free energy of the same ab-initio model in which the electronic energies are computed by the less demanding but less accurate Variational Monte Carlo (VMC) method. Moreover we devise a procedure to correct transition lines based on the use of the new relation. In order to illustrate the procedure, we consider the liquid-liquid phase transition in hydrogen, a first order transition between a lower pressure, molecular and insulating phase and a higher pressure, partially dissociated and conducting phase. We provide new results along the T = 600K isotherm across the phase transition and find good agreement between the transition pressure and specific volumes at coexistence for the model with RQMC accuracy between the prediction of our procedure and the values that can be directly inferred from the observed plateau in the pressure-volume curve along the isotherm. This work paves the way for future use of VMC in first principle simulations of high pressure hydrogen, an essential simplification when considering larger system sizes or quantum proton effects by Path Integral Monte Carlo methods.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electrical conductivity of high-pressure liquid hydrogen by quantum Monte Carlo methods

We compute the electrical conductivity for liquid hydrogen at high pressure using quantum Monte Carlo. The method uses Coupled Electron-Ion Monte Carlo to generate configurations of liquid hydrogen. For each configuration correlated sampling of electrons is performed in order to calculate a set of lowest many-body eigenstates and current-current correlation functions of the system, which are summed over in the many-body Kubo formula to give AC electrical conductivity directly. The extrapolated DC conductivity at 3000 K for several densities shows a liquid semiconductor to liquid-metal transition at high pressure. Our results are in good agreement with shock-wave data.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Equation of state of metallic hydrogen from Coupled Electron-Ion Monte Carlo simulations

We present a study of hydrogen at pressures higher than molecular dissociation using the Coupled Electron-Ion Monte Carlo method. These calculations use the accurate Reptation Quantum Monte Carlo method to estimate the electronic energy and pressure while doing a Monte Carlo simulation of the protons. In addition to presenting simulation results for the equation of state over a large region of phase space, we report the free energy obtained by thermodynamic integration. We find very good agreement with DFT calculations for pressures beyond 600 GPa and densities above $ρ=1.4 g/cm^3$. Both thermodynamic as well as structural properties are accurately reproduced by DFT calculations. This agreement gives a strong support to the different approximations employed in DFT, specifically the approximate exchange-correlation potential and the use of pseudopotentials for the range of densities considered. We find disagreement with chemical models, which suggests a reinvestigation of planetary models, previously constructed using the Saumon-Chabrier-Van Horn equations of state.

0 citations0 reviewsNo code linkedOpen access