Chemically accurate description of aromatic rings interaction using quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Azadi, Sam
We present an accurate study of interactions between benzene molecules using wave function based quantum Monte Carlo (QMC) methods. We compare our QMC results with density functional theory (DFT) using various van der Waals (vdW) functionals. This comparison enables us to tune vdW functionals. We show that highly optimizing the wave function and introducing more dynamical correlation into the wave function are crucial to calculate the weak chemical binding energy between benzene molecules. The good agreement among our results, experiments and quantum chemistry methods, is an important sign of the capability of the wave function based QMC methods to provide accurate description of very weak intermolecular interactions based on vdW dispersive forces.
Accurate band gaps of semiconductors and insulators from Quantum Monte Carlo calculations
NASA Astrophysics Data System (ADS)
Nazarov, Roman; Hood, Randolph; Morales, Miguel
2015-03-01
Ab initio calculations are useful tools in developing materials with targeted band gaps for semiconductor industry. Unfortunately, the main workhorse of ab initio calculations - density functional theory (DFT) in local density approximation (LDA) or generalized gradient approximation (GGA) underestimates band gaps. Several approaches have been proposed starting from empirical corrections to more elaborate exchange-correlation functionals to deal with this problem. But none of these work well for the entire range of semiconductors and insulators. Deficiencies of DFT as a mean field method can be overcome using many-body techniques. Quantum Monte Carlo (QMC) methods can obtain a nearly exact numerical solutions of both total energies and spectral properties. Diffusion Monte Carlo (DMC), the most widely used QMC method, has been shown to provide gold standard results for different material properties, including spectroscopic constants of dimers and clusters, equation of state for solids, accurate descriptions of defects in metals and insulators. To test DMC's accuracy in a wider range of semiconductors and insulators we have computed band gaps of several semiconductors and insulators. We show that DMC can provide superior agreement with experiment compared with more traditional DFT approaches including high level exchange-correlation functionals (e.g. HSE).
Properties of Solar Thermal Fuels by Accurate Quantum Monte Carlo Calculations
NASA Astrophysics Data System (ADS)
Saritas, Kayahan; Ataca, Can; Grossman, Jeffrey C.
2014-03-01
Efficient utilization of the sun as a renewable and clean energy source is one of the major goals of this century due to increasing energy demand and environmental impact. Solar thermal fuels are materials that capture and store the sun's energy in the form of chemical bonds, which can then be released as heat on demand and charged again. Previous work on solar thermal fuels faced challenges related to the cyclability of the fuel over time, as well as the need for higher energy densities. Recently, it was shown that by templating photoswitches onto carbon nanostructures, both high energy density as well as high stability can be achieved. In this work, we explore alternative molecules to azobenzene in such a nano-templated system. We employ the highly accurate quantum Monte Carlo (QMC) method to predict the energy storage potential for each molecule. Our calculations show that in many cases the level of accuracy provided by density functional theory (DFT) is sufficient. However, in some cases, such as dihydroazulene, the drastic change in conjugation upon light absorption causes the DFT predictions to be inconsistent and incorrect. For this case, we compare our QMC results for the geometric structure, band gap and reaction enthalpy with different DFT functionals.
Kanai, Y; Takeuchi, N
2009-10-14
We revisit the molecular line growth mechanism of styrene on the hydrogenated Si(001) 2x1 surface. In particular, we investigate the energetics of the radical chain reaction mechanism by means of diffusion quantum Monte Carlo (QMC) and density functional theory (DFT) calculations. For the exchange correlation (XC) functional we use the non-empirical generalized-gradient approximation (GGA) and meta-GGA. We find that the QMC result also predicts the intra dimer-row growth of the molecular line over the inter dimer-row growth, supporting the conclusion based on DFT results. However, the absolute magnitudes of the adsorption and reaction energies, and the heights of the energy barriers differ considerably between the QMC and DFT with the GGA/meta-GGA XC functionals.
NASA Astrophysics Data System (ADS)
Luo, Ye; Sorella, Sandro
2014-03-01
We introduce a general and efficient method for the calculation of vibrational frequencies of electronic systems, ranging from molecules to solids. By performing damped molecular dynamics with ab initio forces, we show that quantum vibrational frequencies can be evaluated by diagonalizing the time averaged position-position or force-force correlation matrices, although the ionic motion is treated on the classical level within the Born-Oppenheimer approximation. The novelty of our approach is to evaluate atomic forces with QMC by means of a highly accurate and correlated variational wave function which is optimized simultaneously during the dynamics. QMC is an accurate and promising many-body technique for electronic structure calculation thanks to massively parallel computers. However, since infinite statistics is not feasible, property evaluation may be affected by large noise that is difficult to harness. Our approach controls the QMC stochastic bias systematically and gives very accurate results with moderate computational effort, namely even with noisy forces. We prove the accuracy and efficiency of our method on the water monomer[A. Zen et al., JCTC 9 (2013) 4332] and dimer. We are currently working on the challenging problem of simulating liquid water at ambient conditions.
Accurate quantum chemical calculations
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Taylor, Peter R.
1989-01-01
An important goal of quantum chemical calculations is to provide an understanding of chemical bonding and molecular electronic structure. A second goal, the prediction of energy differences to chemical accuracy, has been much harder to attain. First, the computational resources required to achieve such accuracy are very large, and second, it is not straightforward to demonstrate that an apparently accurate result, in terms of agreement with experiment, does not result from a cancellation of errors. Recent advances in electronic structure methodology, coupled with the power of vector supercomputers, have made it possible to solve a number of electronic structure problems exactly using the full configuration interaction (FCI) method within a subspace of the complete Hilbert space. These exact results can be used to benchmark approximate techniques that are applicable to a wider range of chemical and physical problems. The methodology of many-electron quantum chemistry is reviewed. Methods are considered in detail for performing FCI calculations. The application of FCI methods to several three-electron problems in molecular physics are discussed. A number of benchmark applications of FCI wave functions are described. Atomic basis sets and the development of improved methods for handling very large basis sets are discussed: these are then applied to a number of chemical and spectroscopic problems; to transition metals; and to problems involving potential energy surfaces. Although the experiences described give considerable grounds for optimism about the general ability to perform accurate calculations, there are several problems that have proved less tractable, at least with current computer resources, and these and possible solutions are discussed.
Linear-scaling quantum Monte Carlo calculations.
Williamson, A J; Hood, R Q; Grossman, J C
2001-12-10
A method is presented for using truncated, maximally localized Wannier functions to introduce sparsity into the Slater determinant part of the trial wave function in quantum Monte Carlo calculations. When combined with an efficient numerical evaluation of these localized orbitals, the dominant cost in the calculation, namely, the evaluation of the Slater determinant, scales linearly with system size. This technique is applied to accurate total energy calculation of hydrogenated silicon clusters and carbon fullerenes containing 20-1000 valence electrons. PMID:11736525
Quantum Monte Carlo for vibrating molecules
Brown, W.R. |
1996-08-01
Quantum Monte Carlo (QMC) has successfully computed the total electronic energies of atoms and molecules. The main goal of this work is to use correlation function quantum Monte Carlo (CFQMC) to compute the vibrational state energies of molecules given a potential energy surface (PES). In CFQMC, an ensemble of random walkers simulate the diffusion and branching processes of the imaginary-time time dependent Schroedinger equation in order to evaluate the matrix elements. The program QMCVIB was written to perform multi-state VMC and CFQMC calculations and employed for several calculations of the H{sub 2}O and C{sub 3} vibrational states, using 7 PES`s, 3 trial wavefunction forms, two methods of non-linear basis function parameter optimization, and on both serial and parallel computers. In order to construct accurate trial wavefunctions different wavefunctions forms were required for H{sub 2}O and C{sub 3}. In order to construct accurate trial wavefunctions for C{sub 3}, the non-linear parameters were optimized with respect to the sum of the energies of several low-lying vibrational states. In order to stabilize the statistical error estimates for C{sub 3} the Monte Carlo data was collected into blocks. Accurate vibrational state energies were computed using both serial and parallel QMCVIB programs. Comparison of vibrational state energies computed from the three C{sub 3} PES`s suggested that a non-linear equilibrium geometry PES is the most accurate and that discrete potential representations may be used to conveniently determine vibrational state energies.
Quantum Monte Carlo for atoms and molecules
Barnett, R.N.
1989-11-01
The diffusion quantum Monte Carlo with fixed nodes (QMC) approach has been employed in studying energy-eigenstates for 1--4 electron systems. Previous work employing the diffusion QMC technique yielded energies of high quality for H{sub 2}, LiH, Li{sub 2}, and H{sub 2}O. Here, the range of calculations with this new approach has been extended to include additional first-row atoms and molecules. In addition, improvements in the previously computed fixed-node energies of LiH, Li{sub 2}, and H{sub 2}O have been obtained using more accurate trial functions. All computations were performed within, but are not limited to, the Born-Oppenheimer approximation. In our computations, the effects of variation of Monte Carlo parameters on the QMC solution of the Schroedinger equation were studied extensively. These parameters include the time step, renormalization time and nodal structure. These studies have been very useful in determining which choices of such parameters will yield accurate QMC energies most efficiently. Generally, very accurate energies (90--100% of the correlation energy is obtained) have been computed with single-determinant trail functions multiplied by simple correlation functions. Improvements in accuracy should be readily obtained using more complex trial functions.
Accurate Evaluation of Quantum Integrals
NASA Technical Reports Server (NTRS)
Galant, David C.; Goorvitch, D.
1994-01-01
Combining an appropriate finite difference method with Richardson's extrapolation results in a simple, highly accurate numerical method for solving a Schr\\"{o}dinger's equation. Important results are that error estimates are provided, and that one can extrapolate expectation values rather than the wavefunctions to obtain highly accurate expectation values. We discuss the eigenvalues, the error growth in repeated Richardson's extrapolation, and show that the expectation values calculated on a crude mesh can be extrapolated to obtain expectation values of high accuracy.
Optimized trial functions for quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Huang, Sheng-Yu; Sun, Zhiwei; Lester, William A., Jr.
1990-01-01
An algorithm to optimize trial functions for fixed-node quantum Monte Carlo calculations has been developed based on variational random walks. The approach is applied to wave functions that are products of a simple Slater determinant and correlation factor explicitly dependent on interelectronic distance, and is found to provide improved ground-state total energies. A modification of the method for ground-states that makes use of a projection operator technique is shown to make possible the calculation of more accurate excited-state energies. In this optimization method the Young tableaux of the permutation group is used to facilitate the treatment of fermion properties and multiplets. Application to ground states of H2, Li2, H3, H+3, and to the first-excited singlets of H2, H3, and H4 are presented and discussed.
Optimized trial functions for quantum Monte Carlo
Huang, S.; Sun, Z.; Lester, W.A. Jr. )
1990-01-01
An algorithm to optimize trial functions for fixed-node quantum Monte Carlo calculations has been developed based on variational random walks. The approach is applied to wave functions that are products of a simple Slater determinant and correlation factor explicitly dependent on interelectronic distance, and is found to provide improved ground-state total energies. A modification of the method for ground-states that makes use of a projection operator technique is shown to make possible the calculation of more accurate excited-state energies. In this optimization method the Young tableaux of the permutation group is used to facilitate the treatment of fermion properties and multiplets. Application to ground states of H{sub 2}, Li{sub 2}, H{sub 3}, H{sup +}{sub 3}, and to the first-excited singlets of H{sub 2}, H{sub 3}, and H{sub 4} are presented and discussed.
Quantum Monte Carlo simulations of tunneling in quantum adiabatic optimization
NASA Astrophysics Data System (ADS)
Brady, Lucas T.; van Dam, Wim
2016-03-01
We explore to what extent path-integral quantum Monte Carlo methods can efficiently simulate quantum adiabatic optimization algorithms during a quantum tunneling process. Specifically we look at symmetric cost functions defined over n bits with a single potential barrier that a successful quantum adiabatic optimization algorithm will have to tunnel through. The height and width of this barrier depend on n , and by tuning these dependencies, we can make the optimization algorithm succeed or fail in polynomial time. In this article we compare the strength of quantum adiabatic tunneling with that of path-integral quantum Monte Carlo methods. We find numerical evidence that quantum Monte Carlo algorithms will succeed in the same regimes where quantum adiabatic optimization succeeds.
Quantum speedup of Monte Carlo methods
Montanaro, Ashley
2015-01-01
Monte Carlo methods use random sampling to estimate numerical quantities which are hard to compute deterministically. One important example is the use in statistical physics of rapidly mixing Markov chains to approximately compute partition functions. In this work, we describe a quantum algorithm which can accelerate Monte Carlo methods in a very general setting. The algorithm estimates the expected output value of an arbitrary randomized or quantum subroutine with bounded variance, achieving a near-quadratic speedup over the best possible classical algorithm. Combining the algorithm with the use of quantum walks gives a quantum speedup of the fastest known classical algorithms with rigorous performance bounds for computing partition functions, which use multiple-stage Markov chain Monte Carlo techniques. The quantum algorithm can also be used to estimate the total variation distance between probability distributions efficiently. PMID:26528079
Quantum Monte Carlo Endstation for Petascale Computing
Lubos Mitas
2011-01-26
NCSU research group has been focused on accomplising the key goals of this initiative: establishing new generation of quantum Monte Carlo (QMC) computational tools as a part of Endstation petaflop initiative for use at the DOE ORNL computational facilities and for use by computational electronic structure community at large; carrying out high accuracy quantum Monte Carlo demonstration projects in application of these tools to the forefront electronic structure problems in molecular and solid systems; expanding the impact of QMC methods and approaches; explaining and enhancing the impact of these advanced computational approaches. In particular, we have developed quantum Monte Carlo code (QWalk, www.qwalk.org) which was significantly expanded and optimized using funds from this support and at present became an actively used tool in the petascale regime by ORNL researchers and beyond. These developments have been built upon efforts undertaken by the PI's group and collaborators over the period of the last decade. The code was optimized and tested extensively on a number of parallel architectures including petaflop ORNL Jaguar machine. We have developed and redesigned a number of code modules such as evaluation of wave functions and orbitals, calculations of pfaffians and introduction of backflow coordinates together with overall organization of the code and random walker distribution over multicore architectures. We have addressed several bottlenecks such as load balancing and verified efficiency and accuracy of the calculations with the other groups of the Endstation team. The QWalk package contains about 50,000 lines of high quality object-oriented C++ and includes also interfaces to data files from other conventional electronic structure codes such as Gamess, Gaussian, Crystal and others. This grant supported PI for one month during summers, a full-time postdoc and partially three graduate students over the period of the grant duration, it has resulted in 13
Chemical accuracy from quantum Monte Carlo for the benzene dimer
Azadi, Sam; Cohen, R. E.
2015-09-14
We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of −2.3(4) and −2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is −2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.
Chemical accuracy from quantum Monte Carlo for the benzene dimer.
Azadi, Sam; Cohen, R E
2015-09-14
We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods. PMID:26374029
Interaction picture density matrix quantum Monte Carlo
Malone, Fionn D. Lee, D. K. K.; Foulkes, W. M. C.; Blunt, N. S.; Shepherd, James J.; Spencer, J. S.
2015-07-28
The recently developed density matrix quantum Monte Carlo (DMQMC) algorithm stochastically samples the N-body thermal density matrix and hence provides access to exact properties of many-particle quantum systems at arbitrary temperatures. We demonstrate that moving to the interaction picture provides substantial benefits when applying DMQMC to interacting fermions. In this first study, we focus on a system of much recent interest: the uniform electron gas in the warm dense regime. The basis set incompleteness error at finite temperature is investigated and extrapolated via a simple Monte Carlo sampling procedure. Finally, we provide benchmark calculations for a four-electron system, comparing our results to previous work where possible.
Fast quantum Monte Carlo on a GPU
NASA Astrophysics Data System (ADS)
Lutsyshyn, Y.
2015-02-01
We present a scheme for the parallelization of quantum Monte Carlo method on graphical processing units, focusing on variational Monte Carlo simulation of bosonic systems. We use asynchronous execution schemes with shared memory persistence, and obtain an excellent utilization of the accelerator. The CUDA code is provided along with a package that simulates liquid helium-4. The program was benchmarked on several models of Nvidia GPU, including Fermi GTX560 and M2090, and the Kepler architecture K20 GPU. Special optimization was developed for the Kepler cards, including placement of data structures in the register space of the Kepler GPUs. Kepler-specific optimization is discussed.
Novel Quantum Monte Carlo Approaches for Quantum Liquids
NASA Astrophysics Data System (ADS)
Rubenstein, Brenda M.
Quantum Monte Carlo methods are a powerful suite of techniques for solving the quantum many-body problem. By using random numbers to stochastically sample quantum properties, QMC methods are capable of studying low-temperature quantum systems well beyond the reach of conventional deterministic techniques. QMC techniques have likewise been indispensible tools for augmenting our current knowledge of superfluidity and superconductivity. In this thesis, I present two new quantum Monte Carlo techniques, the Monte Carlo Power Method and Bose-Fermi Auxiliary-Field Quantum Monte Carlo, and apply previously developed Path Integral Monte Carlo methods to explore two new phases of quantum hard spheres and hydrogen. I lay the foundation for a subsequent description of my research by first reviewing the physics of quantum liquids in Chapter One and the mathematics behind Quantum Monte Carlo algorithms in Chapter Two. I then discuss the Monte Carlo Power Method, a stochastic way of computing the first several extremal eigenvalues of a matrix too memory-intensive to be stored and therefore diagonalized. As an illustration of the technique, I demonstrate how it can be used to determine the second eigenvalues of the transition matrices of several popular Monte Carlo algorithms. This information may be used to quantify how rapidly a Monte Carlo algorithm is converging to the equilibrium probability distribution it is sampling. I next present the Bose-Fermi Auxiliary-Field Quantum Monte Carlo algorithm. This algorithm generalizes the well-known Auxiliary-Field Quantum Monte Carlo algorithm for fermions to bosons and Bose-Fermi mixtures. Despite some shortcomings, the Bose-Fermi Auxiliary-Field Quantum Monte Carlo algorithm represents the first exact technique capable of studying Bose-Fermi mixtures of any size in any dimension. In Chapter Six, I describe a new Constant Stress Path Integral Monte Carlo algorithm for the study of quantum mechanical systems under high pressures. While
Quantum Monte Carlo calculations of light nuclei
Pieper, S.C.
1998-12-01
Quantum Monte Carlo calculations using realistic two- and three-nucleon interactions are presented for nuclei with up to eight nucleons. We have computed the ground and a few excited states of all such nuclei with Greens function Monte Carlo (GFMC) and all of the experimentally known excited states using variational Monte Carlo (VMC). The GFMC calculations show that for a given Hamiltonian, the VMC calculations of excitation spectra are reliable, but the VMC ground-state energies are significantly above the exact values. We find that the Hamiltonian we are using (which was developed based on {sup 3}H,{sup 4}He, and nuclear matter calculations) underpredicts the binding energy of p-shell nuclei. However our results for excitation spectra are very good and one can see both shell-model and collective spectra resulting from fundamental many-nucleon calculations. Possible improvements in the three-nucleon potential are also be discussed. {copyright} {ital 1998 American Institute of Physics.}
Quantum Monte Carlo calculations of light nuclei
Pieper, Steven C.
1998-12-21
Quantum Monte Carlo calculations using realistic two- and three-nucleon interactions are presented for nuclei with up to eight nucleons. We have computed the ground and a few excited states of all such nuclei with Greens function Monte Carlo (GFMC) and all of the experimentally known excited states using variational Monte Carlo (VMC). The GFMC calculations show that for a given Hamiltonian, the VMC calculations of excitation spectra are reliable, but the VMC ground-state energies are significantly above the exact values. We find that the Hamiltonian we are using (which was developed based on {sup 3}H,{sup 4}He, and nuclear matter calculations) underpredicts the binding energy of p-shell nuclei. However our results for excitation spectra are very good and one can see both shell-model and collective spectra resulting from fundamental many-nucleon calculations. Possible improvements in the three-nucleon potential are also be discussed.
Quantum Monte Carlo calculations of light nuclei.
Pieper, S. C.
1998-08-25
Quantum Monte Carlo calculations using realistic two- and three-nucleon interactions are presented for nuclei with up to eight nucleons. We have computed the ground and a few excited states of all such nuclei with Greens function Monte Carlo (GFMC) and all of the experimentally known excited states using variational Monte Carlo (VMC). The GFMC calculations show that for a given Hamiltonian, the VMC calculations of excitation spectra are reliable, but the VMC ground-state energies are significantly above the exact values. We find that the Hamiltonian we are using (which was developed based on {sup 3}H, {sup 4}He, and nuclear matter calculations) underpredicts the binding energy of p-shell nuclei. However our results for excitation spectra are very good and one can see both shell-model and collective spectra resulting from fundamental many-nucleon calculations. Possible improvements in the three-nucleon potential are also be discussed.
Interaction picture density matrix quantum Monte Carlo.
Malone, Fionn D; Blunt, N S; Shepherd, James J; Lee, D K K; Spencer, J S; Foulkes, W M C
2015-07-28
The recently developed density matrix quantum Monte Carlo (DMQMC) algorithm stochastically samples the N-body thermal density matrix and hence provides access to exact properties of many-particle quantum systems at arbitrary temperatures. We demonstrate that moving to the interaction picture provides substantial benefits when applying DMQMC to interacting fermions. In this first study, we focus on a system of much recent interest: the uniform electron gas in the warm dense regime. The basis set incompleteness error at finite temperature is investigated and extrapolated via a simple Monte Carlo sampling procedure. Finally, we provide benchmark calculations for a four-electron system, comparing our results to previous work where possible. PMID:26233116
Discovering correlated fermions using quantum Monte Carlo.
Wagner, Lucas K; Ceperley, David M
2016-09-01
It has become increasingly feasible to use quantum Monte Carlo (QMC) methods to study correlated fermion systems for realistic Hamiltonians. We give a summary of these techniques targeted at researchers in the field of correlated electrons, focusing on the fundamentals, capabilities, and current status of this technique. The QMC methods often offer the highest accuracy solutions available for systems in the continuum, and, since they address the many-body problem directly, the simulations can be analyzed to obtain insight into the nature of correlated quantum behavior. PMID:27518859
Discovering correlated fermions using quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Wagner, Lucas K.; Ceperley, David M.
2016-09-01
It has become increasingly feasible to use quantum Monte Carlo (QMC) methods to study correlated fermion systems for realistic Hamiltonians. We give a summary of these techniques targeted at researchers in the field of correlated electrons, focusing on the fundamentals, capabilities, and current status of this technique. The QMC methods often offer the highest accuracy solutions available for systems in the continuum, and, since they address the many-body problem directly, the simulations can be analyzed to obtain insight into the nature of correlated quantum behavior.
Applications of Maxent to quantum Monte Carlo
Silver, R.N.; Sivia, D.S.; Gubernatis, J.E. ); Jarrell, M. . Dept. of Physics)
1990-01-01
We consider the application of maximum entropy methods to the analysis of data produced by computer simulations. The focus is the calculation of the dynamical properties of quantum many-body systems by Monte Carlo methods, which is termed the Analytical Continuation Problem.'' For the Anderson model of dilute magnetic impurities in metals, we obtain spectral functions and transport coefficients which obey Kondo Universality.'' 24 refs., 7 figs.
Approaching chemical accuracy with quantum Monte Carlo.
Petruzielo, F R; Toulouse, Julien; Umrigar, C J
2012-03-28
A quantum Monte Carlo study of the atomization energies for the G2 set of molecules is presented. Basis size dependence of diffusion Monte Carlo atomization energies is studied with a single determinant Slater-Jastrow trial wavefunction formed from Hartree-Fock orbitals. With the largest basis set, the mean absolute deviation from experimental atomization energies for the G2 set is 3.0 kcal/mol. Optimizing the orbitals within variational Monte Carlo improves the agreement between diffusion Monte Carlo and experiment, reducing the mean absolute deviation to 2.1 kcal/mol. Moving beyond a single determinant Slater-Jastrow trial wavefunction, diffusion Monte Carlo with a small complete active space Slater-Jastrow trial wavefunction results in near chemical accuracy. In this case, the mean absolute deviation from experimental atomization energies is 1.2 kcal/mol. It is shown from calculations on systems containing phosphorus that the accuracy can be further improved by employing a larger active space. PMID:22462844
Quantum Monte Carlo Methods for First Principles Simulation of Liquid Water
ERIC Educational Resources Information Center
Gergely, John Robert
2009-01-01
Obtaining an accurate microscopic description of water structure and dynamics is of great interest to molecular biology researchers and in the physics and quantum chemistry simulation communities. This dissertation describes efforts to apply quantum Monte Carlo methods to this problem with the goal of making progress toward a fully "ab initio"…
Quantum Monte Carlo calculations of light nuclei.
Pieper, S. C.; Physics
2008-01-01
Variational Monte Carlo and Green's function Monte Carlo are powerful tools for cal- culations of properties of light nuclei using realistic two-nucleon (NN) and three-nucleon (NNN) potentials. Recently the GFMC method has been extended to multiple states with the same quantum numbers. The combination of the Argonne v18 two-nucleon and Illinois-2 three-nucleon potentials gives a good prediction of many energies of nuclei up to 12 C. A number of other recent results are presented: comparison of binding energies with those obtained by the no-core shell model; the incompatibility of modern nuclear Hamiltonians with a bound tetra-neutron; difficulties in computing RMS radii of very weakly bound nuclei, such as 6He; center-of-mass effects on spectroscopic factors; and the possible use of an artificial external well in calculations of neutron-rich isotopes.
Quantum Monte Carlo calculations for carbon nanotubes
NASA Astrophysics Data System (ADS)
Luu, Thomas; Lähde, Timo A.
2016-04-01
We show how lattice quantum Monte Carlo can be applied to the electronic properties of carbon nanotubes in the presence of strong electron-electron correlations. We employ the path-integral formalism and use methods developed within the lattice QCD community for our numerical work. Our lattice Hamiltonian is closely related to the hexagonal Hubbard model augmented by a long-range electron-electron interaction. We apply our method to the single-quasiparticle spectrum of the (3,3) armchair nanotube configuration, and consider the effects of strong electron-electron correlations. Our approach is equally applicable to other nanotubes, as well as to other carbon nanostructures. We benchmark our Monte Carlo calculations against the two- and four-site Hubbard models, where a direct numerical solution is feasible.
Quantum Monte Carlo for the Spectroscopy of Core Excited States
NASA Astrophysics Data System (ADS)
Zubarev, Dmitry
2012-02-01
X-ray absorption spectroscopy is a powerful experimental tool that is capable of delivering valuable information about very delicate aspects of electronic structure and reveals details of the local chemical environment in many systems of fundamental and applied importance. The rigorous interpretation of core-level spectra requires very accurate quantum chemical simulations. The trade-off between feasibility of treatment of large systems and consistency in description of electron correlation tremendously hinders the generation of accurate theoretical results for many experimental studies. We show that the fixed-node diffusion Monte Carlo (FN-DMC) approach can be used straightforwardly for the accurate simulation of core-level spectra. Basic methodological aspects are addressed, including the strategy for the construction of adequate trial wave functions. Examples of FN-DMC calculations of core-level spectra of water and pyrrole are presented. The possibility of the simulation of X-ray absorption spectra of solvent-solute systems is discussed.
Quantum Monte Carlo calculations for light nuclei.
Wiringa, R. B.
1998-10-23
Quantum Monte Carlo calculations of ground and low-lying excited states for nuclei with A {le} 8 are made using a realistic Hamiltonian that fits NN scattering data. Results for more than 40 different (J{pi}, T) states, plus isobaric analogs, are obtained and the known excitation spectra are reproduced reasonably well. Various density and momentum distributions and electromagnetic form factors and moments have also been computed. These are the first microscopic calculations that directly produce nuclear shell structure from realistic NN interactions.
Noncovalent Interactions by Quantum Monte Carlo.
Dubecký, Matúš; Mitas, Lubos; Jurečka, Petr
2016-05-11
Quantum Monte Carlo (QMC) is a family of stochastic methods for solving quantum many-body problems such as the stationary Schrödinger equation. The review introduces basic notions of electronic structure QMC based on random walks in real space as well as its advances and adaptations to systems with noncovalent interactions. Specific issues such as fixed-node error cancellation, construction of trial wave functions, and efficiency considerations that allow for benchmark quality QMC energy differences are described in detail. Comprehensive overview of articles covers QMC applications to systems with noncovalent interactions over the last three decades. The current status of QMC with regard to efficiency, applicability, and usability by nonexperts together with further considerations about QMC developments, limitations, and unsolved challenges are discussed as well. PMID:27081724
Chemical application of diffusion quantum Monte Carlo
NASA Technical Reports Server (NTRS)
Reynolds, P. J.; Lester, W. A., Jr.
1984-01-01
The diffusion quantum Monte Carlo (QMC) method gives a stochastic solution to the Schroedinger equation. This approach is receiving increasing attention in chemical applications as a result of its high accuracy. However, reducing statistical uncertainty remains a priority because chemical effects are often obtained as small differences of large numbers. As an example, the single-triplet splitting of the energy of the methylene molecule CH sub 2 is given. The QMC algorithm was implemented on the CYBER 205, first as a direct transcription of the algorithm running on the VAX 11/780, and second by explicitly writing vector code for all loops longer than a crossover length C. The speed of the codes relative to one another as a function of C, and relative to the VAX, are discussed. The computational time dependence obtained versus the number of basis functions is discussed and this is compared with that obtained from traditional quantum chemistry codes and that obtained from traditional computer architectures.
Quantum Monte Carlo Endstation for Petascale Computing
David Ceperley
2011-03-02
CUDA GPU platform. We restructured the CPU algorithms to express additional parallelism, minimize GPU-CPU communication, and efficiently utilize the GPU memory hierarchy. Using mixed precision on GT200 GPUs and MPI for intercommunication and load balancing, we observe typical full-application speedups of approximately 10x to 15x relative to quad-core Xeon CPUs alone, while reproducing the double-precision CPU results within statistical error. We developed an all-electron quantum Monte Carlo (QMC) method for solids that does not rely on pseudopotentials, and used it to construct a primary ultra-high-pressure calibration based on the equation of state of cubic boron nitride. We computed the static contribution to the free energy with the QMC method and obtained the phonon contribution from density functional theory, yielding a high-accuracy calibration up to 900 GPa usable directly in experiment. We computed the anharmonic Raman frequency shift with QMC simulations as a function of pressure and temperature, allowing optical pressure calibration. In contrast to present experimental approaches, small systematic errors in the theoretical EOS do not increase with pressure, and no extrapolation is needed. This all-electron method is applicable to first-row solids, providing a new reference for ab initio calculations of solids and benchmarks for pseudopotential accuracy. We compared experimental and theoretical results on the momentum distribution and the quasiparticle renormalization factor in sodium. From an x-ray Compton-profile measurement of the valence-electron momentum density, we derived its discontinuity at the Fermi wavevector finding an accurate measure of the renormalization factor that we compared with quantum-Monte-Carlo and G0W0 calculations performed both on crystalline sodium and on the homogeneous electron gas. Our calculated results are in good agreement with the experiment. We have been studying the heat of formation for various Kubas complexes of molecular
Computing Entanglement Entropy in Quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Melko, Roger
2012-02-01
The scaling of entanglement entropy in quantum many-body wavefunctions is expected to be a fruitful resource for studying quantum phases and phase transitions in condensed matter. However, until the recent development of estimators for Renyi entropy in quantum Monte Carlo (QMC), we have been in the dark about the behaviour of entanglement in all but the simplest two-dimensional models. In this talk, I will outline the measurement techniques that allow access to the Renyi entropies in several different QMC methodologies. I will then discuss recent simulation results demonstrating the richness of entanglement scaling in 2D, including: the prevalence of the ``area law''; topological entanglement entropy in a gapped spin liquid; anomalous subleading logarithmic terms due to Goldstone modes; universal scaling at critical points; and examples of emergent conformal-like scaling in several gapless wavefunctions. Finally, I will explore the idea that ``long range entanglement'' may complement the notion of ``long range order'' for quantum phases and phase transitions which lack a conventional order parameter description.
Quantum Monte Carlo methods for nuclear physics
Carlson, J.; Gandolfi, S.; Pederiva, F.; Pieper, Steven C.; Schiavilla, R.; Schmidt, K. E.; Wiringa, R. B.
2015-09-01
Quantum Monte Carlo methods have proved valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments, and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. The nuclear interactions and currents are reviewed along with a description of the continuum quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. A variety of results are presented, including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. Low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars are also described. Furthermore, a coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.
Quantum Monte Carlo methods for nuclear physics
Carlson, J.; Gandolfi, S.; Pederiva, F.; Pieper, Steven C.; Schiavilla, R.; Schmidt, K. E.; Wiringa, R. B.
2015-09-01
Quantum Monte Carlo methods have proved valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments, and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. The nuclear interactions and currents are reviewed along with a description of the continuum quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit,more » and three-body interactions. A variety of results are presented, including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. Low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars are also described. Furthermore, a coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.« less
Quantum Monte Carlo methods for nuclear physics
Carlson, Joseph A.; Gandolfi, Stefano; Pederiva, Francesco; Pieper, Steven C.; Schiavilla, Rocco; Schmidt, K. E,; Wiringa, Robert B.
2014-10-19
Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-bodymore » interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.« less
Metallic lithium by quantum Monte Carlo
Sugiyama, G.; Zerah, G.; Alder, B.J.
1986-12-01
Lithium was chosen as the simplest known metal for the first application of quantum Monte Carlo methods in order to evaluate the accuracy of conventional one-electron band theories. Lithium has been extensively studied using such techniques. Band theory calculations have certain limitations in general and specifically in their application to lithium. Results depend on such factors as charge shape approximations (muffin tins), pseudopotentials (a special problem for lithium where the lack of rho core states requires a strong pseudopotential), and the form and parameters chosen for the exchange potential. The calculations are all one-electron methods in which the correlation effects are included in an ad hoc manner. This approximation may be particularly poor in the high compression regime, where the core states become delocalized. Furthermore, band theory provides only self-consistent results rather than strict limits on the energies. The quantum Monte Carlo method is a totally different technique using a many-body rather than a mean field approach which yields an upper bound on the energies. 18 refs., 4 figs., 1 tab.
Chemical application of diffusion quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Reynolds, P. J.; Lester, W. A., Jr.
1983-10-01
The diffusion quantum Monte Carlo (QMC) method gives a stochastic solution to the Schroedinger equation. As an example the singlet-triplet splitting of the energy of the methylene molecule CH2 is given. The QMC algorithm was implemented on the CYBER 205, first as a direct transcription of the algorithm running on our VAX 11/780, and second by explicitly writing vector code for all loops longer than a crossover length C. The speed of the codes relative to one another as a function of C, and relative to the VAX is discussed. Since CH2 has only eight electrons, most of the loops in this application are fairly short. The longest inner loops run over the set of atomic basis functions. The CPU time dependence obtained versus the number of basis functions is discussed and compared with that obtained from traditional quantum chemistry codes and that obtained from traditional computer architectures. Finally, preliminary work on restructuring the algorithm to compute the separate Monte Carlo realizations in parallel is discussed.
Quantum Monte Carlo methods for nuclear physics
Carlson, J.; Gandolfi, S.; Pederiva, F.; Pieper, Steven C.; Schiavilla, R.; Schmidt, K. E.; Wiringa, R. B.
2015-09-09
Quantum Monte Carlo methods have proved valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments, and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. The nuclear interactions and currents are reviewed along with a description of the continuum quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. A variety of results are presented, including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. Low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars are also described. Furthermore, a coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.
Quantum Monte Carlo methods for nuclear physics
Carlson, Joseph A.; Gandolfi, Stefano; Pederiva, Francesco; Pieper, Steven C.; Schiavilla, Rocco; Schmidt, K. E,; Wiringa, Robert B.
2014-10-19
Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.
Performance of quantum Monte Carlo for calculating molecular bond lengths
NASA Astrophysics Data System (ADS)
Cleland, Deidre M.; Per, Manolo C.
2016-03-01
This work investigates the accuracy of real-space quantum Monte Carlo (QMC) methods for calculating molecular geometries. We present the equilibrium bond lengths of a test set of 30 diatomic molecules calculated using variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC) methods. The effect of different trial wavefunctions is investigated using single determinants constructed from Hartree-Fock (HF) and Density Functional Theory (DFT) orbitals with LDA, PBE, and B3LYP functionals, as well as small multi-configurational self-consistent field (MCSCF) multi-determinant expansions. When compared to experimental geometries, all DMC methods exhibit smaller mean-absolute deviations (MADs) than those given by HF, DFT, and MCSCF. The most accurate MAD of 3 ± 2 × 10-3 Å is achieved using DMC with a small multi-determinant expansion. However, the more computationally efficient multi-determinant VMC method has a similar MAD of only 4.0 ± 0.9 × 10-3 Å, suggesting that QMC forces calculated from the relatively simple VMC algorithm may often be sufficient for accurate molecular geometries.
A pure-sampling quantum Monte Carlo algorithm
Ospadov, Egor; Rothstein, Stuart M.
2015-01-14
The objective of pure-sampling quantum Monte Carlo is to calculate physical properties that are independent of the importance sampling function being employed in the calculation, save for the mismatch of its nodal hypersurface with that of the exact wave function. To achieve this objective, we report a pure-sampling algorithm that combines features of forward walking methods of pure-sampling and reptation quantum Monte Carlo (RQMC). The new algorithm accurately samples properties from the mixed and pure distributions simultaneously in runs performed at a single set of time-steps, over which extrapolation to zero time-step is performed. In a detailed comparison, we found RQMC to be less efficient. It requires different sets of time-steps to accurately determine the energy and other properties, such as the dipole moment. We implement our algorithm by systematically increasing an algorithmic parameter until the properties converge to statistically equivalent values. As a proof in principle, we calculated the fixed-node energy, static α polarizability, and other one-electron expectation values for the ground-states of LiH and water molecules. These quantities are free from importance sampling bias, population control bias, time-step bias, extrapolation-model bias, and the finite-field approximation. We found excellent agreement with the accepted values for the energy and a variety of other properties for those systems.
Understanding Quantum Tunneling through Quantum Monte Carlo Simulations
NASA Astrophysics Data System (ADS)
Boixo, Sergio; Isakov, Sergei; Mazzola, Guglielmo; Smelyanskiy, Vadim; Jiang, Zhang; Neven, Hartmut; Troyer, Matthias
The tunneling between the two ground states of an Ising ferromagnet is a typical example of many-body tunneling processes between two local minima, as they occur during quantum annealing. Performing quantum Monte Carlo (QMC) simulations we find that the QMC tunneling rate displays the same scaling (in the exponent) with system size, as the rate of incoherent tunneling. The scaling in both cases is O (Δ2) , where Δ is the tunneling splitting. An important consequence is that QMC simulations can be used to predict the performance of a quantum annealer for tunneling through a barrier. Furthermore, by using open instead of periodic boundary conditions in imaginary time, equivalent to a projector QMC algorithm, we obtain a quadratic speedup for QMC, and achieve linear scaling in Δ. We provide a physical understanding of these results and their range of applicability based on an instanton picture.
Quantum Monte Carlo simulations in novel geometries
NASA Astrophysics Data System (ADS)
Iglovikov, Vladimir
Quantum Monte Carlo simulations are giving increasing insight into the physics of strongly interacting bosons, spins, and fermions. Initial work focused on the simplest geometries, like a 2D square lattice. Increasingly, modern research is turning to more rich structures such as honeycomb lattice of graphene, the Lieb lattice of the CuO2 planes of cuprate superconductors, the triangular lattice, and coupled layers. These new geometries possess unique features which affect the physics in profound ways, eg a vanishing density of states and relativistic dispersion ("Dirac point'') of a honeycomb lattice, frustration on a triangular lattice, and a flat bands on a Lieb lattice. This thesis concerns both exploring the performance of QMC algorithms on different geometries(primarily via the "sign problem'') and also applying those algorithms to several interesting open problems.
Pseudopotentials for quantum Monte Carlo studies of transition metal oxides
NASA Astrophysics Data System (ADS)
Krogel, Jaron T.; Santana, Juan A.; Reboredo, Fernando A.
2016-02-01
Quantum Monte Carlo (QMC) calculations of transition metal oxides are partially limited by the availability of high-quality pseudopotentials that are both accurate in QMC and compatible with major plane-wave electronic structure codes. We have generated a set of neon-core pseudopotentials with small cutoff radii for the early transition metal elements Sc to Zn within the local density approximation of density functional theory. The pseudopotentials have been directly tested for accuracy within QMC by calculating the first through fourth ionization potentials of the isolated transition metal (M) atoms and the binding curve of each M-O dimer. We find the ionization potentials to be accurate to 0.16(1) eV, on average, relative to experiment. The equilibrium bond lengths of the dimers are within 0.5(1)% of experimental values, on average, and the binding energies are also typically accurate to 0.18(3) eV. The level of accuracy we find for atoms and dimers is comparable to what has recently been observed for bulk metals and oxides using the same pseudopotentials. Our QMC pseudopotential results also compare well with the findings of previous QMC studies and benchmark quantum chemical calculations.
Pseudopotentials for quantum Monte Carlo studies of transition metal oxides
Krogel, Jaron T.; Santana Palacio, Juan A.; Reboredo, Fernando A.
2016-02-22
Quantum Monte Carlo (QMC) calculations of transition metal oxides are partially limited by the availability of high-quality pseudopotentials that are both accurate in QMC and compatible with major plane-wave electronic structure codes. We have generated a set of neon-core pseudopotentials with small cutoff radii for the early transition metal elements Sc to Zn within the local density approximation of density functional theory. The pseudopotentials have been directly tested for accuracy within QMC by calculating the first through fourth ionization potentials of the isolated transition metal (M) atoms and the binding curve of each M-O dimer. We find the ionization potentialsmore » to be accurate to 0.16(1) eV, on average, relative to experiment. The equilibrium bond lengths of the dimers are within 0.5(1)% of experimental values, on average, and the binding energies are also typically accurate to 0.18(3) eV. The level of accuracy we find for atoms and dimers is comparable to what has recently been observed for bulk metals and oxides using the same pseudopotentials. Our QMC pseudopotential results compare well with the findings of previous QMC studies and benchmark quantum chemical calculations.« less
Quantum Monte Carlo studies on small molecules
NASA Astrophysics Data System (ADS)
Galek, Peter T. A.; Handy, Nicholas C.; Lester, William A., Jr.
The Variational Monte Carlo (VMC) and Fixed-Node Diffusion Monte Carlo (FNDMC) methods have been examined, through studies on small molecules. New programs have been written which implement the (by now) standard algorithms for VMC and FNDMC. We have employed and investigated throughout our studies the accuracy of the common Slater-Jastrow trial wave function. Firstly, we have studied a range of sizes of the Jastrow correlation function of the Boys-Handy form, obtained using our optimization program with analytical derivatives of the central moments in the local energy. Secondly, we have studied the effects of Slater-type orbitals (STOs) that display the exact cusp behaviour at nuclei. The orbitals make up the all important trial determinant, which determines the fixed nodal surface. We report all-electron calculations for the ground state energies of Li2, Be2, H2O, NH3, CH4 and H2CO, in all cases but one with accuracy in excess of 95%. Finally, we report an investigation of the ground state energies, dissociation energies and ionization potentials of NH and NH+. Recent focus paid in the literature to these species allow for an extensive comparison with other ab initio methods. We obtain accurate properties for the species and reveal a favourable tendency for fixed-node and other systematic errors to cancel. As a result of our accurate predictions, we are able to obtain a value for the heat of formation of NH, which agrees to within less than 1 kcal mol-1 to other ab initio techniques and 0.2 kcal mol-1 of the experimental value.
Properties of reactive oxygen species by quantum Monte Carlo
Zen, Andrea; Trout, Bernhardt L.; Guidoni, Leonardo
2014-07-07
The electronic properties of the oxygen molecule, in its singlet and triplet states, and of many small oxygen-containing radicals and anions have important roles in different fields of chemistry, biology, and atmospheric science. Nevertheless, the electronic structure of such species is a challenge for ab initio computational approaches because of the difficulties to correctly describe the statical and dynamical correlation effects in presence of one or more unpaired electrons. Only the highest-level quantum chemical approaches can yield reliable characterizations of their molecular properties, such as binding energies, equilibrium structures, molecular vibrations, charge distribution, and polarizabilities. In this work we use the variational Monte Carlo (VMC) and the lattice regularized Monte Carlo (LRDMC) methods to investigate the equilibrium geometries and molecular properties of oxygen and oxygen reactive species. Quantum Monte Carlo methods are used in combination with the Jastrow Antisymmetrized Geminal Power (JAGP) wave function ansatz, which has been recently shown to effectively describe the statical and dynamical correlation of different molecular systems. In particular, we have studied the oxygen molecule, the superoxide anion, the nitric oxide radical and anion, the hydroxyl and hydroperoxyl radicals and their corresponding anions, and the hydrotrioxyl radical. Overall, the methodology was able to correctly describe the geometrical and electronic properties of these systems, through compact but fully-optimised basis sets and with a computational cost which scales as N{sup 3} − N{sup 4}, where N is the number of electrons. This work is therefore opening the way to the accurate study of the energetics and of the reactivity of large and complex oxygen species by first principles.
Properties of reactive oxygen species by quantum Monte Carlo.
Zen, Andrea; Trout, Bernhardt L; Guidoni, Leonardo
2014-07-01
The electronic properties of the oxygen molecule, in its singlet and triplet states, and of many small oxygen-containing radicals and anions have important roles in different fields of chemistry, biology, and atmospheric science. Nevertheless, the electronic structure of such species is a challenge for ab initio computational approaches because of the difficulties to correctly describe the statical and dynamical correlation effects in presence of one or more unpaired electrons. Only the highest-level quantum chemical approaches can yield reliable characterizations of their molecular properties, such as binding energies, equilibrium structures, molecular vibrations, charge distribution, and polarizabilities. In this work we use the variational Monte Carlo (VMC) and the lattice regularized Monte Carlo (LRDMC) methods to investigate the equilibrium geometries and molecular properties of oxygen and oxygen reactive species. Quantum Monte Carlo methods are used in combination with the Jastrow Antisymmetrized Geminal Power (JAGP) wave function ansatz, which has been recently shown to effectively describe the statical and dynamical correlation of different molecular systems. In particular, we have studied the oxygen molecule, the superoxide anion, the nitric oxide radical and anion, the hydroxyl and hydroperoxyl radicals and their corresponding anions, and the hydrotrioxyl radical. Overall, the methodology was able to correctly describe the geometrical and electronic properties of these systems, through compact but fully-optimised basis sets and with a computational cost which scales as N(3) - N(4), where N is the number of electrons. This work is therefore opening the way to the accurate study of the energetics and of the reactivity of large and complex oxygen species by first principles. PMID:25005287
Properties of reactive oxygen species by quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Zen, Andrea; Trout, Bernhardt L.; Guidoni, Leonardo
2014-07-01
The electronic properties of the oxygen molecule, in its singlet and triplet states, and of many small oxygen-containing radicals and anions have important roles in different fields of chemistry, biology, and atmospheric science. Nevertheless, the electronic structure of such species is a challenge for ab initio computational approaches because of the difficulties to correctly describe the statical and dynamical correlation effects in presence of one or more unpaired electrons. Only the highest-level quantum chemical approaches can yield reliable characterizations of their molecular properties, such as binding energies, equilibrium structures, molecular vibrations, charge distribution, and polarizabilities. In this work we use the variational Monte Carlo (VMC) and the lattice regularized Monte Carlo (LRDMC) methods to investigate the equilibrium geometries and molecular properties of oxygen and oxygen reactive species. Quantum Monte Carlo methods are used in combination with the Jastrow Antisymmetrized Geminal Power (JAGP) wave function ansatz, which has been recently shown to effectively describe the statical and dynamical correlation of different molecular systems. In particular, we have studied the oxygen molecule, the superoxide anion, the nitric oxide radical and anion, the hydroxyl and hydroperoxyl radicals and their corresponding anions, and the hydrotrioxyl radical. Overall, the methodology was able to correctly describe the geometrical and electronic properties of these systems, through compact but fully-optimised basis sets and with a computational cost which scales as N3 - N4, where N is the number of electrons. This work is therefore opening the way to the accurate study of the energetics and of the reactivity of large and complex oxygen species by first principles.
Theory and Applications of Quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Deible, Michael John
With the development of peta-scale computers and exa-scale only a few years away, the quantum Monte Carlo (QMC) method, with favorable scaling and inherent parrallelizability, is poised to increase its impact on the electronic structure community. The most widely used variation of QMC is the diffusion Monte Carlo (DMC) method. The accuracy of the DMC method is only limited by the trial wave function that it employs. The effect of the trial wave function is studied here by initially developing correlation-consistent Gaussian basis sets for use in DMC calculations. These basis sets give a low variance in variance Monte Carlo calculations and improved convergence in DMC. The orbital type used in the trial wave function is then investigated, and it is shown that Brueckner orbitals result in a DMC energy comparable to a DMC energy with orbitals from density functional theory and significantly lower than orbitals from Hartree-Fock theory. Three large weakly interacting systems are then studied; a water-16 isomer, a methane clathrate, and a carbon dioxide clathrate. The DMC method is seen to be in good agreement with MP2 calculations and provides reliable benchmarks. Several strongly correlated systems are then studied. An H4 model system that allows for a fine tuning of the multi-configurational character of the wave function shows when the accuracy of the DMC method with a single Slater-determinant trial function begins to deviate from multi-reference benchmarks. The weakly interacting face-to-face ethylene dimer is studied with and without a rotation around the pi bond, which is used to increase the multi-configurational nature of the wave function. This test shows that the effect of a multi-configurational wave function in weakly interacting systems causes DMC with a single Slater-determinant to be unable to achieve sub-chemical accuracy. The beryllium dimer is studied, and it is shown that a very large determinant expansion is required for DMC to predict a binding
Instantons in Quantum Annealing: Thermally Assisted Tunneling Vs Quantum Monte Carlo Simulations
NASA Technical Reports Server (NTRS)
Jiang, Zhang; Smelyanskiy, Vadim N.; Boixo, Sergio; Isakov, Sergei V.; Neven, Hartmut; Mazzola, Guglielmo; Troyer, Matthias
2015-01-01
Recent numerical result (arXiv:1512.02206) from Google suggested that the D-Wave quantum annealer may have an asymptotic speed-up than simulated annealing, however, the asymptotic advantage disappears when it is compared to quantum Monte Carlo (a classical algorithm despite its name). We show analytically that the asymptotic scaling of quantum tunneling is exactly the same as the escape rate in quantum Monte Carlo for a class of problems. Thus, the Google result might be explained in our framework. We also found that the transition state in quantum Monte Carlo corresponds to the instanton solution in quantum tunneling problems, which is observed in numerical simulations.
Quantum Monte Carlo Calculations of Nucleon-Nucleus Scattering
NASA Astrophysics Data System (ADS)
Wiringa, R. B.; Nollett, Kenneth M.; Pieper, Steven C.; Brida, I.
2009-10-01
We report recent quantum Monte Carlo (variational and Green's function) calculations of elastic nucleon-nucleus scattering. We are adding the cases of proton-^4He, neutron-^3H and proton-^3He scattering to a previous GFMC study of neutron-^4He scattering [1]. To do this requires generalizing our methods to include long-range Coulomb forces and to treat coupled channels. The two four-body cases can be compared to other accurate four-body calculational methods such as the AGS equations and hyperspherical harmonic expansions. We will present results for the Argonne v18 interaction alone and with Urbana and Illinois three-nucleon potentials. [4pt] [1] K.M. Nollett, S. C. Pieper, R.B. Wiringa, J. Carlson, and G.M. Hale, Phys. Rev. Lett. 99, 022502 (2007)
Quantum Monte Carlo calculations on positronium compounds
NASA Astrophysics Data System (ADS)
Jiang, Nan
The stability of compounds containing one or more positrons in addition to electrons and nuclei has been the focus of extensive scientific investigations. Interest in these compounds stems from the important role they play in the process of positron annihilation, which has become a useful technique in material science studies. Knowledge of these compounds comes mostly from calculations which are presently less difficult than laboratory experiments. Owing to the small binding energies of these compounds, quantum chemistry methods beyond the molecular orbital approximation must be used. Among them, the quantum Monte Carlo (QMC) method is most appealing because it is easy to implement, gives exact results within the fixed nodes approximation, and makes good use of existing approximate wavefunctions. Applying QMC to small systems like PsH for binding energy calculation is straightforward. To apply it to systems with heavier atoms, to systems for which the center-of-mass motion needs to be separated, and to calculate annihilation rates, special techniques must be developed. In this project a detailed study and several advancements to the QMC method are carried out. Positronium compounds PsH, Ps2, PsO, and Ps2O are studied with algorithms we developed. Results for PsH and Ps2 agree with the best accepted to date. Results for PsO confirm the stability of this compound, and are in fair agreement with an earlier calculation. Results for Ps2O establish the stability of this compound and give an approximate annihilation rate for the first time. Discussions will include an introduction to QMC methods, an in-depth discussion on the QMC formalism, presentation of new algorithms developed in this study, and procedures and results of QMC calculations on the above mentioned positronium compounds.
Quantum Monte Carlo : not just for energy levels.
Nollett, K. M.; Physics
2007-01-01
Quantum Monte Carlo and realistic interactions can provide well-motivated vertices and overlaps for DWBA analyses of reactions. Given an interaction in vaccum, there are several computational approaches to nuclear systems, as you have been hearing: No-core shell model with Lee-Suzuki or Bloch-Horowitz for Hamiltonian Coupled clusters with G-matrix interaction Density functional theory, granted an energy functional derived from the interaction Quantum Monte Carlo - Variational Monte Carlo Green's function Monte Carlo. The last two work directly with a bare interaction and bare operators and describe the wave function without expanding in basis functions, so they have rather different sets of advantages and disadvantages from the others. Variational Monte Carlo (VMC) is built on a sophisticated Ansatz for the wave function, built on shell model like structure modified by operator correlations. Green's function Monte Carlo (GFMC) uses an operator method to project the true ground state out of a reasonable guess wave function.
Systematic study of finite-size effects in quantum Monte Carlo calculations of real metallic systems
Azadi, Sam Foulkes, W. M. C.
2015-09-14
We present a systematic and comprehensive study of finite-size effects in diffusion quantum Monte Carlo calculations of metals. Several previously introduced schemes for correcting finite-size errors are compared for accuracy and efficiency, and practical improvements are introduced. In particular, we test a simple but efficient method of finite-size correction based on an accurate combination of twist averaging and density functional theory. Our diffusion quantum Monte Carlo results for lithium and aluminum, as examples of metallic systems, demonstrate excellent agreement between all of the approaches considered.
Spin-orbit interactions in electronic structure quantum Monte Carlo methods
NASA Astrophysics Data System (ADS)
Melton, Cody A.; Zhu, Minyi; Guo, Shi; Ambrosetti, Alberto; Pederiva, Francesco; Mitas, Lubos
2016-04-01
We develop generalization of the fixed-phase diffusion Monte Carlo method for Hamiltonians which explicitly depends on particle spins such as for spin-orbit interactions. The method is formulated in a zero-variance manner and is similar to the treatment of nonlocal operators in commonly used static-spin calculations. Tests on atomic and molecular systems show that it is very accurate, on par with the fixed-node method. This opens electronic structure quantum Monte Carlo methods to a vast research area of quantum phenomena in which spin-related interactions play an important role.
Linear Scaling Quantum Monte Carlo Calculations
NASA Astrophysics Data System (ADS)
Williamson, Andrew
2002-03-01
New developments to the quantum Monte Carlo approach are presented that improve the scaling of the time required to calculate the total energy of a configuration of electronic coordinates from N^3 to nearly linear[1]. The first factor of N is achieved by applying a unitary transform to the set of single particle orbitals used to construct the Slater determinant, creating a set of maximally localized Wannier orbitals. These localized functions are then truncated beyond a given cutoff radius to introduce sparsity into the Slater determinant. The second factor of N is achieved by evaluating the maximally localized Wannier orbitals on a cubic spline grid, which removes the size dependence of the basis set (e.g. plane waves, Gaussians) typically used to expand the orbitals. Application of this method to the calculation of the binding energy of carbon fullerenes and silicon nanostructures will be presented. An extension of the approach to deal with excited states of systems will also be presented in the context of the calculation of the excitonic gap of a variety of systems. This work was performed under the auspices of the U.S. Dept. of Energy at the University of California/LLNL under contract no. W-7405-Eng-48. [1] A.J. Williamson, R.Q. Hood and J.C. Grossman, Phys. Rev. Lett. 87 246406 (2001)
Quantum Monte Carlo Simulations of Correlated-Electron Models
NASA Astrophysics Data System (ADS)
Zhang, Shiwei
1996-05-01
We briefly review quantum Monte Carlo simulation methods for strongly correlated fermion systems and the well-known ``sign'' problem that plagues these methods. We then discuss recent efforts to overcome the problem in the context of simulations of lattice models of electron correlations. In particular, we describe a new algorithm^1, called the constrained path Monte Carlo (CPMC), for studying ground-state (T=0K) properties. It has the form of a random walk in a space of mean-field solutions (Slater determinants); the exponential decay of ``sign'' or signal-to-noise ratio is eliminated by constraining the paths of the random walk according to a known trial wave function. Applications of this algorithm to the Hubbard model have enabled accurate and systematic studies of correlation functions, including s- and d-wave pairings, and hence the long-standing problem of the model's relevance to superconductivity. The method is directly applicable to a variety of other models important to understand high-Tc superconductors and heavy-fermion compounds. In addition, it is expected to be useful to simulations of nuclei, atoms, molecules, and solids. We also comment on possible extensions of the algorithm to finite-temperature calculations. Work supported in part by the Department of Energy's High Performance Computing and Communication Program at Los Alamos National Laboratory, and at OSU by DOE-Basic Energy Sciences, Division of Materials Sciences. ^1 Shiwei Zhang, J. Carlson, and J. E. Gubernatis, Phys. Rev. Lett. 74, 3652 (1995).
Improving light propagation Monte Carlo simulations with accurate 3D modeling of skin tissue
Paquit, Vincent C; Price, Jeffery R; Meriaudeau, Fabrice; Tobin Jr, Kenneth William
2008-01-01
In this paper, we present a 3D light propagation model to simulate multispectral reflectance images of large skin surface areas. In particular, we aim to simulate more accurately the effects of various physiological properties of the skin in the case of subcutaneous vein imaging compared to existing models. Our method combines a Monte Carlo light propagation model, a realistic three-dimensional model of the skin using parametric surfaces and a vision system for data acquisition. We describe our model in detail, present results from the Monte Carlo modeling and compare our results with those obtained with a well established Monte Carlo model and with real skin reflectance images.
NASA Astrophysics Data System (ADS)
Park, Sungjin; Shin, Hyeondeok; Kwon, Yongkyung
2012-08-01
The recently-proposed fourth-order propagator based on the multi-product expansion has been applied to path-integral Monte Carlo calculations for asymmetric quantum quadruploar rotors fixed at face-centered cubic lattice sites. The rotors are observed to undergo an orientational orderdisorder phase transition at a low temperature when the electric quadrupole-quadrupole interaction is strong enough. At intermediate interaction strength, a further decrease of temperature after the first transition to the ordered phase results in a reentrant transition back to the disordered phase. The theoretical phase diagram of these asymmetric rotors determined by using fourth-order path-integral Monte Carlo calculations is found to be in good quantitative agreement with the experimental one for solid hydrogen deuteride. This leads us to conclude that the fourth-order propagator can be effectively implemented for an accurate path-integral Monte Carlo calculation of a quantum many-body system with rotational degrees of freedom.
Quantum Monte Carlo Calculations Applied to Magnetic Molecules
Larry Engelhardt
2006-08-09
We have calculated the equilibrium thermodynamic properties of Heisenberg spin systems using a quantum Monte Carlo (QMC) method. We have used some of these systems as models to describe recently synthesized magnetic molecules, and-upon comparing the results of these calculations with experimental data-have obtained accurate estimates for the basic parameters of these models. We have also performed calculations for other systems that are of more general interest, being relevant both for existing experimental data and for future experiments. Utilizing the concept of importance sampling, these calculations can be carried out in an arbitrarily large quantum Hilbert space, while still avoiding any approximations that would introduce systematic errors. The only errors are statistical in nature, and as such, their magnitudes are accurately estimated during the course of a simulation. Frustrated spin systems present a major challenge to the QMC method, nevertheless, in many instances progress can be made. In this chapter, the field of magnetic molecules is introduced, paying particular attention to the characteristics that distinguish magnetic molecules from other systems that are studied in condensed matter physics. We briefly outline the typical path by which we learn about magnetic molecules, which requires a close relationship between experiments and theoretical calculations. The typical experiments are introduced here, while the theoretical methods are discussed in the next chapter. Each of these theoretical methods has a considerable limitation, also described in Chapter 2, which together serve to motivate the present work. As is shown throughout the later chapters, the present QMC method is often able to provide useful information where other methods fail. In Chapter 3, the use of Monte Carlo methods in statistical physics is reviewed, building up the fundamental ideas that are necessary in order to understand the method that has been used in this work. With these
Overy, Catherine; Blunt, N. S.; Shepherd, James J.; Booth, George H.; Cleland, Deidre; Alavi, Ali
2014-12-28
Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the full configuration interaction quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself) and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems.
Accurately modeling Gaussian beam propagation in the context of Monte Carlo techniques
NASA Astrophysics Data System (ADS)
Hokr, Brett H.; Winblad, Aidan; Bixler, Joel N.; Elpers, Gabriel; Zollars, Byron; Scully, Marlan O.; Yakovlev, Vladislav V.; Thomas, Robert J.
2016-03-01
Monte Carlo simulations are widely considered to be the gold standard for studying the propagation of light in turbid media. However, traditional Monte Carlo methods fail to account for diffraction because they treat light as a particle. This results in converging beams focusing to a point instead of a diffraction limited spot, greatly effecting the accuracy of Monte Carlo simulations near the focal plane. Here, we present a technique capable of simulating a focusing beam in accordance to the rules of Gaussian optics, resulting in a diffraction limited focal spot. This technique can be easily implemented into any traditional Monte Carlo simulation allowing existing models to be converted to include accurate focusing geometries with minimal effort. We will present results for a focusing beam in a layered tissue model, demonstrating that for different scenarios the region of highest intensity, thus the greatest heating, can change from the surface to the focus. The ability to simulate accurate focusing geometries will greatly enhance the usefulness of Monte Carlo for countless applications, including studying laser tissue interactions in medical applications and light propagation through turbid media.
Cohesion Energetics of Carbon Allotropes: Quantum Monte Carlo Study
Shin, Hyeondeok; Kang, Sinabro; Koo, Jahyun; Lee, Hoonkyung; Kim, Jeongnim; Kwon, Yongkyung
2014-01-01
We have performed quantum Monte Carlo calculations to study the cohesion energetics of carbon allotropes, including sp3-bonded diamond, sp2-bonded graphene, sp-sp2 hybridized graphynes, and sp-bonded carbyne. The comput- ed cohesive energies of diamond and graphene are found to be in excellent agreement with the corresponding values de- termined experimentally for diamond and graphite, respectively, when the zero-point energies, along with the interlayer binding in the case of graphite, are included. We have also found that the cohesive energy of graphyne decreases system- atically as the ratio of sp-bonded carbon atoms increases. The cohesive energy of -graphyne, the most energetically- stable graphyne, turns out to be 6.766(6) eV/atom, which is smaller than that of graphene by 0.698(12) eV/atom. Experi- mental difficulty in synthesizing graphynes could be explained by their significantly smaller cohesive energies. Finally we conclude that the cohesive energy of a newly-proposed two-dimensional carbon network can be accurately estimated with the carbon-carbon bond energies determined from the cohesive energies of graphene and three different graphynes.
Pseudopotentials for quantum Monte Carlo calculations of transition metal oxides
NASA Astrophysics Data System (ADS)
Krogel, Jaron; Santana, Juan; Kent, Paul; Reboredo, Fernando
2015-03-01
Quantum Monte Carlo calculations of transition metal oxides are partially limited by the availability of high quality pseudopotentials that are both accurate in QMC and compatible with major electronic structure codes, e.g. by not being overly hard in the standard planewave basis. Following insight gained from recent GW calculations, a set of neon core pseudopotentials with small cutoff radii have been created for the early transition metal elements Sc to Zn within the local density approximation of DFT. The pseudopotentials have been tested for energy consistency within QMC by calculating the first through fourth ionization potentials of the isolated transition metal (TM) atoms and the binding curve of each TM-O dimer. The vast majority of the ionization potentials fall within 0.3 eV of the experimental values, with exceptions occurring mainly for atoms with multiple unpaired d electrons where multireference effects are the strongest. The equilibrium bond lengths of the dimers are within 1% of experimental values and the binding energy errors are typically less than 0.3 eV. Given the uniform treatment of the core, the larger deviations occasionally observed may primarily reflect the limitations of a Slater-Jastrow trial wavefunction. This work is supported by the Materials Sciences & Engineering Division of the Office of Basic Energy Sciences, U.S. DOE. Research by PRCK was conducted at the Center for Nanophase Materials Sciences, which is a DOE Office of Science User Facility.
Quantum Monte Carlo Simulation of condensed van der Waals Systems
NASA Astrophysics Data System (ADS)
Benali, Anouar; Shulenburger, Luke; Romero, Nichols A.; Kim, Jeongnim; Anatole von Lilienfeld, O.
2012-02-01
Van der Waals forces are as ubiquitous as infamous. While post-Hartree-Fock methods enable accurate estimates of these forces in molecules and clusters, they remain elusive for dealing with many-electron condensed phase systems. We present Quantum Monte Carlo [1,2] results for condensed van der Waals systems. Interatomic many-body contributions to cohesive energies and bulk modulus will be discussed. Numerical evidence is presented for crystals of rare gas atoms, and compared to experiments and methods [3]. Sandia National Laboratories is a multiprogram laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. DoE's National Nuclear Security Administration under Contract No. DE-AC04-94AL85000.[4pt] [1] J. Kim, K. Esler, J. McMinis and D. Ceperley, SciDAC 2010, J. of Physics: Conference series, Chattanooga, Tennessee, July 11 2011 [0pt] [2] QMCPACK simulation suite, http://qmcpack.cmscc.org (unpublished)[0pt] [3] O. A. von Lillienfeld and A. Tkatchenko, J. Chem. Phys. 132 234109 (2010)
Quantum Monte Carlo Simulation of Tunneling Devices Using Bohm Trajectories
NASA Astrophysics Data System (ADS)
Oriols, X.; García-García, J. J.; Martín, F.; Suñé, J.; González, T.; Mateos, J.; Pardo, D.
1997-11-01
A generalization of the classical Monte Carlo (MC) device simulation technique is proposed to simultaneously deal with quantum-mechanical phase-coherence effects and scattering interactions in tunneling devices. The proposed method restricts the quantum treatment of transport to the regions of the device where the potential profile significantly changes in distances of the order of the de Broglie wavelength of the carriers (the quantum window). Bohm trajectories associated to time-dependent Gaussian wavepackets are used to simulate the electron transport in the quantum window. Outside this window, the classical ensemble simulation technique is used. Classical and quantum trajectories are smoothly matched at the boundaries of the quantum window according to a criterium of total energy conservation. A simple one-dimensional simulator for resonant tunneling diodes is presented to demonstrate the feasibility of our proposal.
Monte Carlo simulation of quantum Zeno effect in the brain
NASA Astrophysics Data System (ADS)
Georgiev, Danko
2015-12-01
Environmental decoherence appears to be the biggest obstacle for successful construction of quantum mind theories. Nevertheless, the quantum physicist Henry Stapp promoted the view that the mind could utilize quantum Zeno effect to influence brain dynamics and that the efficacy of such mental efforts would not be undermined by environmental decoherence of the brain. To address the physical plausibility of Stapp's claim, we modeled the brain using quantum tunneling of an electron in a multiple-well structure such as the voltage sensor in neuronal ion channels and performed Monte Carlo simulations of quantum Zeno effect exerted by the mind upon the brain in the presence or absence of environmental decoherence. The simulations unambiguously showed that the quantum Zeno effect breaks down for timescales greater than the brain decoherence time. To generalize the Monte Carlo simulation results for any n-level quantum system, we further analyzed the change of brain entropy due to the mind probing actions and proved a theorem according to which local projections cannot decrease the von Neumann entropy of the unconditional brain density matrix. The latter theorem establishes that Stapp's model is physically implausible but leaves a door open for future development of quantum mind theories provided the brain has a decoherence-free subspace.
Monte Carlo studies of nuclei and quantum liquid drops
Pandharipande, V.R.; Pieper, S.C.
1989-01-01
The progress in application of variational and Green's function Monte Carlo methods to nuclei is reviewed. The nature of single-particle orbitals in correlated quantum liquid drops is discussed, and it is suggested that the difference between quasi-particle and mean-field orbitals may be of importance in nuclear structure physics. 27 refs., 7 figs., 2 tabs.
Reagents for Electrophilic Amination: A Quantum Monte CarloStudy
Amador-Bedolla, Carlos; Salomon-Ferrer, Romelia; Lester Jr.,William A.; Vazquez-Martinez, Jose A.; Aspuru-Guzik, Alan
2006-11-01
Electroamination is an appealing synthetic strategy toconstruct carbon-nitrogen bonds. We explore the use of the quantum MonteCarlo method and a proposed variant of the electron-pair localizationfunction--the electron-pair localization function density--as a measureof the nucleophilicity of nitrogen lone-pairs as a possible screeningprocedure for electrophilic reagents.
Quantum Monte Carlo simulation with a black hole
NASA Astrophysics Data System (ADS)
Benić, Sanjin; Yamamoto, Arata
2016-05-01
We perform quantum Monte Carlo simulations in the background of a classical black hole. The lattice discretized path integral is numerically calculated in the Schwarzschild metric and in its approximated metric. We study spontaneous symmetry breaking of a real scalar field theory. We observe inhomogeneous symmetry breaking induced by an inhomogeneous gravitational field.
Monte Carlo techniques for real-time quantum dynamics
Dowling, Mark R. . E-mail: dowling@physics.uq.edu.au; Davis, Matthew J.; Drummond, Peter D.; Corney, Joel F.
2007-01-10
The stochastic-gauge representation is a method of mapping the equation of motion for the quantum mechanical density operator onto a set of equivalent stochastic differential equations. One of the stochastic variables is termed the 'weight', and its magnitude is related to the importance of the stochastic trajectory. We investigate the use of Monte Carlo algorithms to improve the sampling of the weighted trajectories and thus reduce sampling error in a simulation of quantum dynamics. The method can be applied to calculations in real time, as well as imaginary time for which Monte Carlo algorithms are more-commonly used. The Monte-Carlo algorithms are applicable when the weight is guaranteed to be real, and we demonstrate how to ensure this is the case. Examples are given for the anharmonic oscillator, where large improvements over stochastic sampling are observed.
Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes
Sjostrom, Travis; Daligault, Jerome
2014-10-10
Here, we develop and implement a new quantum molecular dynamics approximation that allows fast and accurate simulations of dense plasmas from cold to hot conditions. The method is based on a carefully designed orbital-free implementation of density functional theory. The results for hydrogen and aluminum are in very good agreement with Kohn-Sham (orbital-based) density functional theory and path integral Monte Carlo calculations for microscopic features such as the electron density as well as the equation of state. The present approach does not scale with temperature and hence extends to higher temperatures than is accessible in the Kohn-Sham method and lower temperatures than is accessible by path integral Monte Carlo calculations, while being significantly less computationally expensive than either of those two methods.
Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes
Sjostrom, Travis; Daligault, Jerome
2014-10-10
Here, we develop and implement a new quantum molecular dynamics approximation that allows fast and accurate simulations of dense plasmas from cold to hot conditions. The method is based on a carefully designed orbital-free implementation of density functional theory. The results for hydrogen and aluminum are in very good agreement with Kohn-Sham (orbital-based) density functional theory and path integral Monte Carlo calculations for microscopic features such as the electron density as well as the equation of state. The present approach does not scale with temperature and hence extends to higher temperatures than is accessible in the Kohn-Sham method and lowermore » temperatures than is accessible by path integral Monte Carlo calculations, while being significantly less computationally expensive than either of those two methods.« less
NASA Astrophysics Data System (ADS)
García-García, J.; Martín, F.; Oriols, X.; Suñé, J.
1998-12-01
A tool for the simulation of resonant tunneling diodes (RTDs) has been developed. This is based on the solution of the quantum Liouville equation in the active region of the device and the Boltzman transport equation in the regions adjacent to the contacts by means of a Monte Carlo algorithm. By accurately coupling both approaches to current transport, we have developed a quantum simulation tool that allows the use of simulation domains much larger and realistic than those previously considered, without a significant increase in computational burden. The main characteristics expected for the considered devices are clearly obtained, thus supporting the validity of our tool for the simulation of RTDs.
NASA Astrophysics Data System (ADS)
Lonardoni, D.; Pederiva, F.; Gandolfi, S.
2014-01-01
Background: An accurate assessment of the hyperon-nucleon interaction is of great interest in view of recent observations of very massive neutron stars. The challenge is to build a realistic interaction that can be used over a wide range of masses and in infinite matter starting from the available experimental data on the binding energy of light hypernuclei. To this end, accurate calculations of the hyperon binding energy in a hypernucleus are necessary. Purpose: We present a quantum Monte Carlo study of Λ and ΛΛ hypernuclei up to A =91. We investigate the contribution of two- and three-body Λ-nucleon forces to the Λ binding energy. Method: Ground state energies are computed solving the Schrödinger equation for nonrelativistic baryons by means of the auxiliary field diffusion Monte Carlo algorithm extended to the hypernuclear sector. Results: We show that a simple adjustment of the parameters of the ΛNN three-body force yields a very good agreement with available experimental data over a wide range of hypernuclear masses. In some cases no experiments have been performed yet, and we give new predictions. Conclusions: The newly fitted ΛNN force properly describes the physics of medium-heavy Λ hypernuclei, correctly reproducing the saturation property of the hyperon separation energy.
Accuracy of electronic wave functions in quantum Monte Carlo: The effect of high-order correlations
NASA Astrophysics Data System (ADS)
Huang, Chien-Jung; Umrigar, C. J.; Nightingale, M. P.
1997-08-01
Compact and accurate wave functions can be constructed by quantum Monte Carlo methods. Typically, these wave functions consist of a sum of a small number of Slater determinants multiplied by a Jastrow factor. In this paper we study the importance of including high-order, nucleus-three-electron correlations in the Jastrow factor. An efficient algorithm based on the theory of invariants is used to compute the high-body correlations. We observe significant improvements in the variational Monte Carlo energy and in the fluctuations of the local energies but not in the fixed-node diffusion Monte Carlo energies. Improvements for the ground states of physical, fermionic atoms are found to be smaller than those for the ground states of fictitious, bosonic atoms, indicating that errors in the nodal surfaces of the fermionic wave functions are a limiting factor.
Molecular hydrogen adsorbed on benzene: Insights from a quantum Monte Carlo study.
Beaudet, Todd D; Casula, Michele; Kim, Jeongnim; Sorella, Sandro; Martin, Richard M
2008-10-28
We present a quantum Monte Carlo study of the hydrogen-benzene system where binding is very weak. We demonstrate that the binding is well described at both variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC) levels by a Jastrow correlated single determinant geminal wave function with an optimized compact basis set that includes diffuse orbitals. Agreement between VMC and fixed-node DMC binding energies is found to be within 0.18 mhartree, suggesting that the calculations are well converged with respect to the basis. Essentially the same binding is also found in independent DMC calculations using a different trial wave function of a more conventional Slater-Jastrow form, supporting our conclusion that the binding energy is accurate and includes all effects of correlation. We compare with previous calculations, and we discuss the physical mechanisms of the interaction, the role of diffuse basis functions, and the charge redistribution in the bond. PMID:19045302
Low-pressure phase diagram of crystalline benzene from quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Azadi, Sam; Cohen, R. E.
2016-08-01
We studied the low-pressure (0-10 GPa) phase diagram of crystalline benzene using quantum Monte Carlo and density functional theory (DFT) methods. We performed diffusion quantum Monte Carlo (DMC) calculations to obtain accurate static phase diagrams as benchmarks for modern van der Waals density functionals. Using density functional perturbation theory, we computed the phonon contributions to the free energies. Our DFT enthalpy-pressure phase diagrams indicate that the Pbca and P21/c structures are the most stable phases within the studied pressure range. The DMC Gibbs free-energy calculations predict that the room temperature Pbca to P21/c phase transition occurs at 2.1(1) GPa. This prediction is consistent with available experimental results at room temperature. Our DMC calculations give 50.6 ± 0.5 kJ/mol for crystalline benzene lattice energy.
Ab initio quantum Monte Carlo calculations of ground-state properties of manganese's oxides
NASA Astrophysics Data System (ADS)
Sharma, Vinit; Krogel, Jaron T.; Kent, P. R. C.; Reboredo, Fernando A.
One of the critical scientific challenges of contemporary research is to obtain an accurate theoretical description of the electronic properties of strongly correlated systems such as transition metal oxides and rare-earth compounds, since state-of-art ab-initio methods based on approximate density functionals are not always sufficiently accurate. Quantum Monte Carlo (QMC) methods, which use statistical sampling to evaluate many-body wave functions, have the potential to answer this challenge. Owing to the few fundamental approximations made and the direct treatment of electron correlation, QMC methods are among the most accurate electronic structure methods available to date. We assess the accuracy of the diffusion Monte Carlo method in the case of rocksalt manganese oxide (MnO). We study the electronic properties of this strongly-correlated oxide, which has been identified as a suitable candidate for many applications ranging from catalysts to electronic devices. ``This work was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division.'' Ab initio quantum Monte Carlo calculations of ground-state properties of manganese's oxides.
Quantum Monte Carlo calculations for point defects in semiconductors
NASA Astrophysics Data System (ADS)
Hennig, Richard
2010-03-01
Point defects in silicon have been studied extensively for many years. Nevertheless the mechanism for self diffusion in Si is still debated. Direct experimental measurements of the selfdiffusion in silicon are complicated by the lack of suitable isotopes. Formation energies are either obtained from theory or indirectly through the analysis of dopant and metal diffusion experiments. Density functional calculations predict formation energies ranging from 3 to 5 eV depending on the approximations used for the exchange-correlation functional [1]. Analysis of dopant and metal diffusion experiments result in similar broad range of diffusion activation energies of 4.95 [2], 4.68 [3], 2.4 eV [4]. Assuming a migration energy barrier of 0.1-0.3 eV [5], the resulting experimental interstitial formation energies range from 2.1 - 4.9 eV. To answer the question of the formation energy of Si interstitials we resort to a many-body description of the wave functions using quantum Monte Carlo (QMC) techniques. Previous QMC calculations resulted in formation energies for the interstitials of around 5 eV [1,6]. We present a careful analysis of all the controlled and uncontrolled approximations that affect the defect formation energies in variational and diffusion Monte Carlo calculations. We find that more accurate trial wave functions for QMC using improved Jastrow expansions and most importantly a backflow transformation for the electron coordinates significantly improve the wave functions. Using zero-variance extrapolation, we predict interstitial formation energies in good agreement with hybrid DFT functionals [1] and recent GW calculations [7]. [4pt] [1] E. R. Batista, J. Heyd, R. G. Hennig, B. P. Uberuaga, R. L. Martin, G. E. Scuseria, C. J. Umrigar, and J. W. Wilkins. Phys. Rev. B 74, 121102(R) (2006).[0pt] [2] H. Bracht, E. E. Haller, and R. Clark-Phelps, Phys. Rev. Lett. 81, 393 (1998). [0pt] [3] A. Ural, P. B. Griffin, and J. D. Plummer, Phys. Rev. Lett. 83, 3454 (1999). [0pt
An accurate and simple quantum model for liquid water.
Paesani, Francesco; Zhang, Wei; Case, David A; Cheatham, Thomas E; Voth, Gregory A
2006-11-14
The path-integral molecular dynamics and centroid molecular dynamics methods have been applied to investigate the behavior of liquid water at ambient conditions starting from a recently developed simple point charge/flexible (SPC/Fw) model. Several quantum structural, thermodynamic, and dynamical properties have been computed and compared to the corresponding classical values, as well as to the available experimental data. The path-integral molecular dynamics simulations show that the inclusion of quantum effects results in a less structured liquid with a reduced amount of hydrogen bonding in comparison to its classical analog. The nuclear quantization also leads to a smaller dielectric constant and a larger diffusion coefficient relative to the corresponding classical values. Collective and single molecule time correlation functions show a faster decay than their classical counterparts. Good agreement with the experimental measurements in the low-frequency region is obtained for the quantum infrared spectrum, which also shows a higher intensity and a redshift relative to its classical analog. A modification of the original parametrization of the SPC/Fw model is suggested and tested in order to construct an accurate quantum model, called q-SPC/Fw, for liquid water. The quantum results for several thermodynamic and dynamical properties computed with the new model are shown to be in a significantly better agreement with the experimental data. Finally, a force-matching approach was applied to the q-SPC/Fw model to derive an effective quantum force field for liquid water in which the effects due to the nuclear quantization are explicitly distinguished from those due to the underlying molecular interactions. Thermodynamic and dynamical properties computed using standard classical simulations with this effective quantum potential are found in excellent agreement with those obtained from significantly more computationally demanding full centroid molecular dynamics
Energies of the first row atoms from quantum Monte Carlo.
Brown, M D; Trail, J R; Ríos, P López; Needs, R J
2007-06-14
All-electron variational and diffusion quantum Monte Carlo calculations of the ground state energies of the first row atoms (from Li to Ne) are reported. The authors use trial wave functions of four types: single-determinant Slater-Jastrow wave functions, multideterminant Slater-Jastrow wave functions, single-determinant Slater-Jastrow wave functions with backflow transformations, and multideterminant Slater-Jastrow wave functions with backflow transformations. At the diffusion quantum Monte Carlo level and using their multideterminant Slater-Jastrow wave functions with backflow transformations, they recover 99% or more of the correlation energies for Li, Be, B, C, N, and Ne, 97% for O, and 98% for F. PMID:17581047
Takahashi, F; Endo, A
2007-01-01
A system utilising radiation transport codes has been developed to derive accurate dose distributions in a human body for radiological accidents. A suitable model is quite essential for a numerical analysis. Therefore, two tools were developed to setup a 'problem-dependent' input file, defining a radiation source and an exposed person to simulate the radiation transport in an accident with the Monte Carlo calculation codes-MCNP and MCNPX. Necessary resources are defined by a dialogue method with a generally used personal computer for both the tools. The tools prepare human body and source models described in the input file format of the employed Monte Carlo codes. The tools were validated for dose assessment in comparison with a past criticality accident and a hypothesized exposure. PMID:17510203
Valence-bond quantum Monte Carlo algorithms defined on trees.
Deschner, Andreas; Sørensen, Erik S
2014-09-01
We present a class of algorithms for performing valence-bond quantum Monte Carlo of quantum spin models. Valence-bond quantum Monte Carlo is a projective T=0 Monte Carlo method based on sampling of a set of operator strings that can be viewed as forming a treelike structure. The algorithms presented here utilize the notion of a worm that moves up and down this tree and changes the associated operator string. In quite general terms, we derive a set of equations whose solutions correspond to a whole class of algorithms. As specific examples of this class of algorithms, we focus on two cases. The bouncing worm algorithm, for which updates are always accepted by allowing the worm to bounce up and down the tree, and the driven worm algorithm, where a single parameter controls how far up the tree the worm reaches before turning around. The latter algorithm involves only a single bounce where the worm turns from going up the tree to going down. The presence of the control parameter necessitates the introduction of an acceptance probability for the update. PMID:25314561
GORRAM: Introducing accurate operational-speed radiative transfer Monte Carlo solvers
NASA Astrophysics Data System (ADS)
Buras-Schnell, Robert; Schnell, Franziska; Buras, Allan
2016-06-01
We present a new approach for solving the radiative transfer equation in horizontally homogeneous atmospheres. The motivation was to develop a fast yet accurate radiative transfer solver to be used in operational retrieval algorithms for next generation meteorological satellites. The core component is the program GORRAM (Generator Of Really Rapid Accurate Monte-Carlo) which generates solvers individually optimized for the intended task. These solvers consist of a Monte Carlo model capable of path recycling and a representative set of photon paths. Latter is generated using the simulated annealing technique. GORRAM automatically takes advantage of limitations on the variability of the atmosphere. Due to this optimization the number of photon paths necessary for accurate results can be reduced by several orders of magnitude. For the shown example of a forward model intended for an aerosol satellite retrieval, comparison with an exact yet slow solver shows that a precision of better than 1% can be achieved with only 36 photons. The computational time is at least an order of magnitude faster than any other type of radiative transfer solver. Merely the lookup table approach often used in satellite retrieval is faster, but on the other hand suffers from limited accuracy. This makes GORRAM-generated solvers an eligible candidate as forward model in operational-speed retrieval algorithms and data assimilation applications. GORRAM also has the potential to create fast solvers of other integrable equations.
Monte Carlo modeling provides accurate calibration factors for radionuclide activity meters.
Zagni, F; Cicoria, G; Lucconi, G; Infantino, A; Lodi, F; Marengo, M
2014-12-01
Accurate determination of calibration factors for radionuclide activity meters is crucial for quantitative studies and in the optimization step of radiation protection, as these detectors are widespread in radiopharmacy and nuclear medicine facilities. In this work we developed the Monte Carlo model of a widely used activity meter, using the Geant4 simulation toolkit. More precisely the "PENELOPE" EM physics models were employed. The model was validated by means of several certified sources, traceable to primary activity standards, and other sources locally standardized with spectrometry measurements, plus other experimental tests. Great care was taken in order to accurately reproduce the geometrical details of the gas chamber and the activity sources, each of which is different in shape and enclosed in a unique container. Both relative calibration factors and ionization current obtained with simulations were compared against experimental measurements; further tests were carried out, such as the comparison of the relative response of the chamber for a source placed at different positions. The results showed a satisfactory level of accuracy in the energy range of interest, with the discrepancies lower than 4% for all the tested parameters. This shows that an accurate Monte Carlo modeling of this type of detector is feasible using the low-energy physics models embedded in Geant4. The obtained Monte Carlo model establishes a powerful tool for first instance determination of new calibration factors for non-standard radionuclides, for custom containers, when a reference source is not available. Moreover, the model provides an experimental setup for further research and optimization with regards to materials and geometrical details of the measuring setup, such as the ionization chamber itself or the containers configuration. PMID:25195174
Fast Monte Carlo Electron-Photon Transport Method and Application in Accurate Radiotherapy
NASA Astrophysics Data System (ADS)
Hao, Lijuan; Sun, Guangyao; Zheng, Huaqing; Song, Jing; Chen, Zhenping; Li, Gui
2014-06-01
Monte Carlo (MC) method is the most accurate computational method for dose calculation, but its wide application on clinical accurate radiotherapy is hindered due to its poor speed of converging and long computation time. In the MC dose calculation research, the main task is to speed up computation while high precision is maintained. The purpose of this paper is to enhance the calculation speed of MC method for electron-photon transport with high precision and ultimately to reduce the accurate radiotherapy dose calculation time based on normal computer to the level of several hours, which meets the requirement of clinical dose verification. Based on the existing Super Monte Carlo Simulation Program (SuperMC), developed by FDS Team, a fast MC method for electron-photon coupled transport was presented with focus on two aspects: firstly, through simplifying and optimizing the physical model of the electron-photon transport, the calculation speed was increased with slightly reduction of calculation accuracy; secondly, using a variety of MC calculation acceleration methods, for example, taking use of obtained information in previous calculations to avoid repeat simulation of particles with identical history; applying proper variance reduction techniques to accelerate MC method convergence rate, etc. The fast MC method was tested by a lot of simple physical models and clinical cases included nasopharyngeal carcinoma, peripheral lung tumor, cervical carcinoma, etc. The result shows that the fast MC method for electron-photon transport was fast enough to meet the requirement of clinical accurate radiotherapy dose verification. Later, the method will be applied to the Accurate/Advanced Radiation Therapy System ARTS as a MC dose verification module.
Continuous-time quantum Monte Carlo impurity solvers
NASA Astrophysics Data System (ADS)
Gull, Emanuel; Werner, Philipp; Fuchs, Sebastian; Surer, Brigitte; Pruschke, Thomas; Troyer, Matthias
2011-04-01
Continuous-time quantum Monte Carlo impurity solvers are algorithms that sample the partition function of an impurity model using diagrammatic Monte Carlo techniques. The present paper describes codes that implement the interaction expansion algorithm originally developed by Rubtsov, Savkin, and Lichtenstein, as well as the hybridization expansion method developed by Werner, Millis, Troyer, et al. These impurity solvers are part of the ALPS-DMFT application package and are accompanied by an implementation of dynamical mean-field self-consistency equations for (single orbital single site) dynamical mean-field problems with arbitrary densities of states. Program summaryProgram title: dmft Catalogue identifier: AEIL_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEIL_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: ALPS LIBRARY LICENSE version 1.1 No. of lines in distributed program, including test data, etc.: 899 806 No. of bytes in distributed program, including test data, etc.: 32 153 916 Distribution format: tar.gz Programming language: C++ Operating system: The ALPS libraries have been tested on the following platforms and compilers: Linux with GNU Compiler Collection (g++ version 3.1 and higher), and Intel C++ Compiler (icc version 7.0 and higher) MacOS X with GNU Compiler (g++ Apple-version 3.1, 3.3 and 4.0) IBM AIX with Visual Age C++ (xlC version 6.0) and GNU (g++ version 3.1 and higher) compilers Compaq Tru64 UNIX with Compq C++ Compiler (cxx) SGI IRIX with MIPSpro C++ Compiler (CC) HP-UX with HP C++ Compiler (aCC) Windows with Cygwin or coLinux platforms and GNU Compiler Collection (g++ version 3.1 and higher) RAM: 10 MB-1 GB Classification: 7.3 External routines: ALPS [1], BLAS/LAPACK, HDF5 Nature of problem: (See [2].) Quantum impurity models describe an atom or molecule embedded in a host material with which it can exchange electrons. They are basic to nanoscience as
Barborini, Matteo; Sorella, Sandro; Guidoni, Leonardo
2014-01-01
We present full structural optimizations of the ground state and of the low lying triplet state of the ethylene molecule by means of Quantum Monte Carlo methods. Using the efficient structural optimization method based on renormalization techniques and on adjoint differentiation algorithms recently proposed [Sorella, S.; Capriotti, L. J. Chem. Phys. 2010, 133, 234111], we present the variational convergence of both wave function parameters and atomic positions. All of the calculations were done using an accurate and compact wave function based on Pauling’s resonating valence bond representation: the Jastrow Antisymmetrized Geminal Power (JAGP). All structural and wave function parameters are optimized, including coefficients and exponents of the Gaussian primitives of the AGP and the Jastrow atomic orbitals. Bond lengths and bond angles are calculated with a statistical error of about 0.1% and are in good agreement with the available experimental data. The Variational and Diffusion Monte Carlo calculations estimate vertical and adiabatic excitation energies in the ranges 4.623(10)–4.688(5) eV and 3.001(5)–3.091(5) eV, respectively. The adiabatic gap, which is in line with other correlated quantum chemistry methods, is slightly higher than the value estimated by recent photodissociation experiments. Our results demonstrate how Quantum Monte Carlo calculations have become a promising and computationally affordable tool for the structural optimization of correlated molecular systems. PMID:24634617
Sign Learning Kink-based (SiLK) Quantum Monte Carlo for molecular systems.
Ma, Xiaoyao; Hall, Randall W; Löffler, Frank; Kowalski, Karol; Bhaskaran-Nair, Kiran; Jarrell, Mark; Moreno, Juana
2016-01-01
The Sign Learning Kink (SiLK) based Quantum Monte Carlo (QMC) method is used to calculate the ab initio ground state energies for multiple geometries of the H2O, N2, and F2 molecules. The method is based on Feynman's path integral formulation of quantum mechanics and has two stages. The first stage is called the learning stage and reduces the well-known QMC minus sign problem by optimizing the linear combinations of Slater determinants which are used in the second stage, a conventional QMC simulation. The method is tested using different vector spaces and compared to the results of other quantum chemical methods and to exact diagonalization. Our findings demonstrate that the SiLK method is accurate and reduces or eliminates the minus sign problem. PMID:26747795
Sign Learning Kink-based (SiLK) Quantum Monte Carlo for molecular systems
NASA Astrophysics Data System (ADS)
Ma, Xiaoyao; Hall, Randall W.; Löffler, Frank; Kowalski, Karol; Bhaskaran-Nair, Kiran; Jarrell, Mark; Moreno, Juana
2016-01-01
The Sign Learning Kink (SiLK) based Quantum Monte Carlo (QMC) method is used to calculate the ab initio ground state energies for multiple geometries of the H2O, N2, and F2 molecules. The method is based on Feynman's path integral formulation of quantum mechanics and has two stages. The first stage is called the learning stage and reduces the well-known QMC minus sign problem by optimizing the linear combinations of Slater determinants which are used in the second stage, a conventional QMC simulation. The method is tested using different vector spaces and compared to the results of other quantum chemical methods and to exact diagonalization. Our findings demonstrate that the SiLK method is accurate and reduces or eliminates the minus sign problem.
Semi-stochastic full configuration interaction quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Holmes, Adam; Petruzielo, Frank; Khadilkar, Mihir; Changlani, Hitesh; Nightingale, M. P.; Umrigar, C. J.
2012-02-01
In the recently proposed full configuration interaction quantum Monte Carlo (FCIQMC) [1,2], the ground state is projected out stochastically, using a population of walkers each of which represents a basis state in the Hilbert space spanned by Slater determinants. The infamous fermion sign problem manifests itself in the fact that walkers of either sign can be spawned on a given determinant. We propose an improvement on this method in the form of a hybrid stochastic/deterministic technique, which we expect will improve the efficiency of the algorithm by ameliorating the sign problem. We test the method on atoms and molecules, e.g., carbon, carbon dimer, N2 molecule, and stretched N2. [4pt] [1] Fermion Monte Carlo without fixed nodes: a Game of Life, death and annihilation in Slater Determinant space. George Booth, Alex Thom, Ali Alavi. J Chem Phys 131, 050106, (2009).[0pt] [2] Survival of the fittest: Accelerating convergence in full configuration-interaction quantum Monte Carlo. Deidre Cleland, George Booth, and Ali Alavi. J Chem Phys 132, 041103 (2010).
a Quantum Monte Carlo Study of the High Pressure Phases of Solid Hydrogen
NASA Astrophysics Data System (ADS)
Natoli, Vincent Dominic
1994-01-01
Variational and Diffusion Monte Carlo are powerful computational methods which can afford accurate estimates of the ground state properties of quantum many-body problems. We have applied these Monte Carlo methods to the high pressure phases of solid hydrogen to elucidate those parts of the phase diagram where experimental results are inconclusive or lacking. The method allows us to treat both electrons and protons as quantum particles by incorporating them in the trial wavefunction and avoids the Born-Oppenheimer and harmonic approximations. Our trial wavefunction uses single-body solutions from a mean-field calculation coupled with standard pair potential terms to achieve the most accurate results to date. Equally accurate results were realized for calculations in the disparate insulating molecular and metallic atomic regime. We performed a study of the possible ground state structures of the atomic metallic phase of hydrogen which identifies a new family of low energy atomic structures. Another study was done on the molecular phase over the range of pressures (40-180GPa) where recent experiments have observed spectral discontinuities and other interesting new phenomena. Particular attention was directed to determining the equation of state and orientational ordering. We find that molecular hydrogen adopts a lower symmetry insulating structure over a wide range of pressure. The results of the atomic and molecular studies are combined to draw conclusions about the molecular-atomic transition and other details about the high pressure phase diagram.
Hood, R Q; Williamson, A J; Dubois, J L; Reboredo, F A
2008-02-07
We have developed a highly accurate computational capability to calculate the equation of state (EOS) and defect formation energies of metallic systems. We are using a newly developed algorithm that enables the study of metallic systems with quantum Monte Carlo (QMC) methods. To date, technical limitations have restricted the application of QMC methods to semiconductors, insulators and the homogeneous electron gas. Using this new 'QMC for metals' we can determine, for the first time, the significance of correlation effects in the EOS and in the formation energies of point defects, impurities, surfaces and interfaces in metallic systems. These calculations go beyond the state-of-the-art accuracy which is currently obtained with Density Functional Theory approaches. Such benchmark calculations can provide more accurate predictions for the EOS and the formation energies of vacancies and interstitials in simple metals. These are important parameters in determining the mechanical properties as well as the micro-structural evolution of metals in irradiated materials or under extreme conditions. We describe the development of our 'QMC for metals' code, which has been adapted to run efficiently on a variety of computer architectures including BG/L. We present results of the first accurate quantum Monte Carlo calculation of an EOS of a realistic metallic system that goes beyond the homogeneous electron gas.
Infinite variance in fermion quantum Monte Carlo calculations
NASA Astrophysics Data System (ADS)
Shi, Hao; Zhang, Shiwei
2016-03-01
For important classes of many-fermion problems, quantum Monte Carlo (QMC) methods allow exact calculations of ground-state and finite-temperature properties without the sign problem. The list spans condensed matter, nuclear physics, and high-energy physics, including the half-filled repulsive Hubbard model, the spin-balanced atomic Fermi gas, and lattice quantum chromodynamics calculations at zero density with Wilson Fermions, and is growing rapidly as a number of problems have been discovered recently to be free of the sign problem. In these situations, QMC calculations are relied on to provide definitive answers. Their results are instrumental to our ability to understand and compute properties in fundamental models important to multiple subareas in quantum physics. It is shown, however, that the most commonly employed algorithms in such situations have an infinite variance problem. A diverging variance causes the estimated Monte Carlo statistical error bar to be incorrect, which can render the results of the calculation unreliable or meaningless. We discuss how to identify the infinite variance problem. An approach is then proposed to solve the problem. The solution does not require major modifications to standard algorithms, adding a "bridge link" to the imaginary-time path integral. The general idea is applicable to a variety of situations where the infinite variance problem may be present. Illustrative results are presented for the ground state of the Hubbard model at half-filling.
Minimising biases in full configuration interaction quantum Monte Carlo.
Vigor, W A; Spencer, J S; Bearpark, M J; Thom, A J W
2015-03-14
We show that Full Configuration Interaction Quantum Monte Carlo (FCIQMC) is a Markov chain in its present form. We construct the Markov matrix of FCIQMC for a two determinant system and hence compute the stationary distribution. These solutions are used to quantify the dependence of the population dynamics on the parameters defining the Markov chain. Despite the simplicity of a system with only two determinants, it still reveals a population control bias inherent to the FCIQMC algorithm. We investigate the effect of simulation parameters on the population control bias for the neon atom and suggest simulation setups to, in general, minimise the bias. We show a reweight ing scheme to remove the bias caused by population control commonly used in diffusion Monte Carlo [Umrigar et al., J. Chem. Phys. 99, 2865 (1993)] is effective and recommend its use as a post processing step. PMID:25770522
Machine learning of parameters for accurate semiempirical quantum chemical calculations
Dral, Pavlo O.; von Lilienfeld, O. Anatole; Thiel, Walter
2015-04-14
We investigate possible improvements in the accuracy of semiempirical quantum chemistry (SQC) methods through the use of machine learning (ML) models for the parameters. For a given class of compounds, ML techniques require sufficiently large training sets to develop ML models that can be used for adapting SQC parameters to reflect changes in molecular composition and geometry. The ML-SQC approach allows the automatic tuning of SQC parameters for individual molecules, thereby improving the accuracy without deteriorating transferability to molecules with molecular descriptors very different from those in the training set. The performance of this approach is demonstrated for the semiempiricalmore » OM2 method using a set of 6095 constitutional isomers C7H10O2, for which accurate ab initio atomization enthalpies are available. The ML-OM2 results show improved average accuracy and a much reduced error range compared with those of standard OM2 results, with mean absolute errors in atomization enthalpies dropping from 6.3 to 1.7 kcal/mol. They are also found to be superior to the results from specific OM2 reparameterizations (rOM2) for the same set of isomers. The ML-SQC approach thus holds promise for fast and reasonably accurate high-throughput screening of materials and molecules.« less
Machine learning of parameters for accurate semiempirical quantum chemical calculations
Dral, Pavlo O.; von Lilienfeld, O. Anatole; Thiel, Walter
2015-04-14
We investigate possible improvements in the accuracy of semiempirical quantum chemistry (SQC) methods through the use of machine learning (ML) models for the parameters. For a given class of compounds, ML techniques require sufficiently large training sets to develop ML models that can be used for adapting SQC parameters to reflect changes in molecular composition and geometry. The ML-SQC approach allows the automatic tuning of SQC parameters for individual molecules, thereby improving the accuracy without deteriorating transferability to molecules with molecular descriptors very different from those in the training set. The performance of this approach is demonstrated for the semiempirical OM2 method using a set of 6095 constitutional isomers C_{7}H_{10}O_{2}, for which accurate ab initio atomization enthalpies are available. The ML-OM2 results show improved average accuracy and a much reduced error range compared with those of standard OM2 results, with mean absolute errors in atomization enthalpies dropping from 6.3 to 1.7 kcal/mol. They are also found to be superior to the results from specific OM2 reparameterizations (rOM2) for the same set of isomers. The ML-SQC approach thus holds promise for fast and reasonably accurate high-throughput screening of materials and molecules.
Constrained Path Quantum Monte Carlo Method for Fermion Ground States
NASA Astrophysics Data System (ADS)
Zhang, Shiwei; Carlson, J.; Gubernatis, J. E.
1995-05-01
We propose a new quantum Monte Carlo algorithm to compute fermion ground-state properties. The ground state is projected from an initial wave function by a branching random walk in an over-complete basis space of Slater determinants. By constraining the determinants according to a trial wave function \\|ΨT>, we remove the exponential decay of signal-to-noise ratio characteristic of the sign problem. The method is variational and is exact if \\|ΨT> is exact. We report results on the two-dimensional Hubbard model up to size 16×16, for various electron fillings and interaction strengths.
Monte Carlo simulation of a noisy quantum channel with memory.
Akhalwaya, Ismail; Moodley, Mervlyn; Petruccione, Francesco
2015-10-01
The classical capacity of quantum channels is well understood for channels with uncorrelated noise. For the case of correlated noise, however, there are still open questions. We calculate the classical capacity of a forgetful channel constructed by Markov switching between two depolarizing channels. Techniques have previously been applied to approximate the output entropy of this channel and thus its capacity. In this paper, we use a Metropolis-Hastings Monte Carlo approach to numerically calculate the entropy. The algorithm is implemented in parallel and its performance is studied and optimized. The effects of memory on the capacity are explored and previous results are confirmed to higher precision. PMID:26565361
Quantum Monte Carlo calculations of neutron matter with chiral three-body forces
NASA Astrophysics Data System (ADS)
Tews, I.; Gandolfi, S.; Gezerlis, A.; Schwenk, A.
2016-02-01
Chiral effective field theory (EFT) enables a systematic description of low-energy hadronic interactions with controlled theoretical uncertainties. For strongly interacting systems, quantum Monte Carlo (QMC) methods provide some of the most accurate solutions, but they require as input local potentials. We have recently constructed local chiral nucleon-nucleon (NN) interactions up to next-to-next-to-leading order (N2LO ). Chiral EFT naturally predicts consistent many-body forces. In this paper, we consider the leading chiral three-nucleon (3N) interactions in local form. These are included in auxiliary field diffusion Monte Carlo (AFDMC) simulations. We present results for the equation of state of neutron matter and for the energies and radii of neutron drops. In particular, we study the regulator dependence at the Hartree-Fock level and in AFDMC and find that present local regulators lead to less repulsion from 3N forces compared to the usual nonlocal regulators.
Spinor path integral Quantum Monte Carlo for fermions
NASA Astrophysics Data System (ADS)
Shin, Daejin; Yousif, Hosam; Shumway, John
2007-03-01
We have developed a continuous-space path integral method for spin 1/2 fermions with fixed-phase approximation. The internal spin degrees of freedom of each particle is represented by four extra dimensions. This effectively maps each spinor onto two of the excited states of a four dimensional harmonic oscillator. The phases that appear in the problem can be treated within the fixed-phase approximation. This mapping preserves rotational invariance and allows us to treat spin interactions and fermionic exchange on equal footing, which may lead to new theoretical insights. The technique is illustrated for a few simple models, including a spin in a magnetic field and interacting electrons in a quantum dot in a magnetic field at finite temperature. We will discuss possible extensions of the method to molecules and solids using variational and diffusion Quantum Monte Carlo.
Quantum Monte Carlo simulations with tensor-network states
NASA Astrophysics Data System (ADS)
Song, Jeong Pil; Clay, R. T.
2011-03-01
Matrix-product states, generated by the density-matrix renormalization group method, are among the most powerful methods for simulation of quasi-one dimensional quantum systems. Direct application of a matrix-product state representation fails for two dimensional systems, although a number of tensor-network states have been proposed to generalize the concept for two dimensions. We introduce a useful approximate method replacing a 4-index tensor by two matrices in order to contract tensors in two dimensions. We use this formalism as a basis for variational quantum Monte Carlo, optimizing the matrix elements stochastically. We present results on a two dimensional spinless fermion model including nearest- neighbor Coulomb interactions, and determine the critical Coulomb interaction for the charge density wave state by finite size scaling. This work was supported by the Department of Energy grant DE-FG02-06ER46315.
Correlated wavefunction quantum Monte Carlo approach to solids
Louie, S.G.
1992-10-01
A method for calculating the electronic and structural properties of solids using correlated wavefunctions together with quantum Monte Carlo techniques is described. The approach retains the exact Coulomb interaction between the electrons and employs a many-electron wavefunction of the Jastrow-Slater form. Several examples are given to illustrate the utility of the method. Topics discussed include the cohesive properties of bulk semiconductors, the magnetic-field- induced Wigner crystal in two dimensions, and the magnetic structure of bcc hydrogen. Landau level mixing is shown to be important in determining the transition between the fractional quantum Hall liquid and the Wigner crystal. Information on electron correlations such as the pair correlation functions which are not accessible to one- electron theories is also obtained. 24 refs, 5 figs, 1 tab.
Two-Dimensional Ferromagnet: Quantum Monte Carlo results
NASA Astrophysics Data System (ADS)
Henelius, Patrik; Timm, Carsten; Girvin, Steven M.; Sandvik, Anders
1997-03-01
In the quantum Hall system the Zeeman interaction between electronic spins and the external magnetic field is typically weak compared to both the Landau-level splitting and the exchange interaction. Therefore, quantum Hall systems at integer filling factors can be ferromagnets. The magnetization and, recently, the nuclear magnetic relaxation rate 1/T1 have been measured for these magnets.(S.E. Barrett et al.), Phys. Rev. Lett. 72, 1368 (1994); 74, 5112 (1995) These quantities have been calculated in a Schwinger-boson mean-field approach.(N. Read and S. Sachdev, Phys. Rev. Lett. 75), 3509 (1995) We have calculated these same quantities using a Stochastic Series Expansion Monte Carlo Method. The results are compared with the experimental data, the mean-field results and with 1/N corrections for the mean-field results, calculated by our group.
Excited states of methylene from quantum Monte Carlo.
Zimmerman, Paul M; Toulouse, Julien; Zhang, Zhiyong; Musgrave, Charles B; Umrigar, C J
2009-09-28
The ground and lowest three adiabatic excited states of methylene are computed using the variational Monte Carlo and diffusion Monte Carlo (DMC) methods using progressively larger Jastrow-Slater multideterminant complete active space (CAS) wave functions. The highest of these states has the same symmetry, (1)A(1), as the first excited state. The DMC excitation energies obtained using any of the CAS wave functions are in excellent agreement with experiment, but single-determinant wave functions do not yield accurate DMC energies of the states of (1)A(1) symmetry, indicating that it is important to include in the wave function Slater determinants that describe static (strong) correlation. Excitation energies obtained using recently proposed pseudopotentials [Burkatzki et al., J. Chem. Phys. 126, 234105 (2007)] differ from the all-electron excitation energies by at most 0.04 eV. PMID:19791848
Quantum Monte Carlo theory and applications for molecular systems
NASA Astrophysics Data System (ADS)
Kollias, Alexander C.
New directions for the quantum Monte Carlo (QMC) electronic structure method are discussed. Diffusion Monte Carlo (DMC) results for the atomization energy and heats for formation of CO+2 are presented, while the bonding character is examined using the electron localization function. DMC all-electron and effective-core potential trial functions are used to obtain the atomization energies, heats of formation, and energy differences of the C2H4 singlet and triplet states. In addition, DMC is applied to obtain the heat of reaction and barrier height of the proton extraction reaction, CH3OH + Cl → CH 2OH + HCl. The results of the barrier height and heat of reaction are verified by examining the atomization energies and heats for formation of the reactants and products. DMC calculations were carried out on 22 small hydrocarbons. In this benchmark study the DMC atomization and bond dissociation energies, and heats of formation of these hydrocarbons are presented and compared to other ab initio methods. Methods for geometry optimization and calculating forces for QMC are discussed. The response surface methodology is applied to variational Monte Carlo (VMC) and DMC methods to obtain an optimized geometry, force constants and vibrational frequencies of CH2O. Finally, the zero-variance principle is applied to obtain VMC and DMC effective-core potential force estimators. These estimators are used to obtain a force curve for LiH.
Quantum Monte Carlo calculations of {Alpha} = 8 nuclei.
Wiringa, R. B.; Pieper, S. C.; Carlson, J.; Pandharipande, V. R.; Physics; LANL; Univ. of Illinois
2000-07-01
We report quantum Monte Carlo calculations of ground and low-lying excited states for {Alpha}=8 nuclei using a realistic Hamiltonian containing the Argonne v{sub 18} two-nucleon and Urbana IX three-nucleon potentials. The calculations begin with correlated eight-body wave functions that have a filled {alpha}-like core and four p-shell nucleons LS coupled to the appropriate (J{sup {pi}},T) quantum numbers for the state of interest. After optimization, these variational wave functions are used as input to a Green's function Monte Carlo calculation made with a new constrained path algorithm. We find that the Hamiltonian produces a {sup 8}Be ground state that is within 2 MeV of the experimental resonance, but the other eight-body energies are progressively worse as the neutron-proton asymmetry increases. The {sup 8}Li ground state is stable against breakup into subclusters, but the {sup 8}He ground state is not. The excited state spectra are in fair agreement with experiment, with both the single-particle behavior of {sup 8}He and {sup 8}Li and the collective rotational behavior of {sup 8}Be being reproduced. We also examine energy differences in the T=1,2 isomultiplets and isospin-mixing matrix elements in the excited states of {sup 8}Be. Finally, we present densities, momentum distributions, and studies of the intrinsic shapes of these nuclei, with {sup 8}Be exhibiting a definite 2{alpha} cluster structure.
Variational quantum Monte Carlo calculations for solid surfaces
Bahnsen, R.; Eckstein, H.; Schattke, W.; Fitzer, N.; Redmer, R.
2001-06-15
Quantum Monte Carlo methods have proven to predict atomic and bulk properties of light and nonlight elements with high accuracy. Here we report on variational quantum Monte Carlo (VMC) calculations for solid surfaces. Taking the boundary condition for the simulation from a finite-layer geometry, the Hamiltonian, including a nonlocal pseudopotential, is cast in a layer-resolved form and evaluated with a two-dimensional Ewald summation technique. The exact cancellation of all jellium contributions to the Hamiltonian is ensured. The many-body trial wave function consists of a Slater determinant with parametrized localized orbitals and a Jastrow factor with a common two-body term plus an additional confinement term representing further variational freedom to take into account the existence of the surface. We present results for the ideal (110) surface of gallium arsenide for different system sizes. With the optimized trial wave function, we determine some properties related to a solid surface to illustrate that VMC techniques provide reasonable results under full inclusion of many-body effects at solid surfaces.
H2 Adsorption in a Porous Crystal: Accurate First-Principles Quantum Simulation.
D'Arcy, Jordan H; Jordan, Meredith J T; Frankcombe, Terry J; Collins, Michael A
2015-12-17
A general method is presented for constructing, from ab initio quantum chemistry calculations, the potential energy surface (PES) for H2 absorbed in a porous crystalline material. The method is illustrated for the metal-organic framework material MOF-5. Rigid body quantum diffusion Monte Carlo simulations are used in the construction of the PES and to evaluate the quantum ground state of H2 in MOF-5, the zero-point energy, and the enthalpy of adsorption at 0 K. PMID:26322374
NASA Astrophysics Data System (ADS)
Garrison, Stephen L.
2005-07-01
The combination of molecular simulations and potentials obtained from quantum chemistry is shown to be able to provide reasonably accurate thermodynamic property predictions. Gibbs ensemble Monte Carlo simulations are used to understand the effects of small perturbations to various regions of the model Lennard-Jones 12-6 potential. However, when the phase behavior and second virial coefficient are scaled by the critical properties calculated for each potential, the results obey a corresponding states relation suggesting a non-uniqueness problem for interaction potentials fit to experimental phase behavior. Several variations of a procedure collectively referred to as quantum mechanical Hybrid Methods for Interaction Energies (HM-IE) are developed and used to accurately estimate interaction energies from CCSD(T) calculations with a large basis set in a computationally efficient manner for the neon-neon, acetylene-acetylene, and nitrogen-benzene systems. Using these results and methods, an ab initio, pairwise-additive, site-site potential for acetylene is determined and then improved using results from molecular simulations using this initial potential. The initial simulation results also indicate that a limited range of energies important for accurate phase behavior predictions. Second virial coefficients calculated from the improved potential indicate that one set of experimental data in the literature is likely erroneous. This prescription is then applied to methanethiol. Difficulties in modeling the effects of the lone pair electrons suggest that charges on the lone pair sites negatively impact the ability of the intermolecular potential to describe certain orientations, but that the lone pair sites may be necessary to reasonably duplicate the interaction energies for several orientations. Two possible methods for incorporating the effects of three-body interactions into simulations within the pairwise-additivity formulation are also developed. A low density
NASA Astrophysics Data System (ADS)
Dinpajooh, Mohammadhasan; Bai, Peng; Allan, Douglas A.; Siepmann, J. Ilja
2015-09-01
Since the seminal paper by Panagiotopoulos [Mol. Phys. 61, 813 (1997)], the Gibbs ensemble Monte Carlo (GEMC) method has been the most popular particle-based simulation approach for the computation of vapor-liquid phase equilibria. However, the validity of GEMC simulations in the near-critical region has been questioned because rigorous finite-size scaling approaches cannot be applied to simulations with fluctuating volume. Valleau [Mol. Simul. 29, 627 (2003)] has argued that GEMC simulations would lead to a spurious overestimation of the critical temperature. More recently, Patel et al. [J. Chem. Phys. 134, 024101 (2011)] opined that the use of analytical tail corrections would be problematic in the near-critical region. To address these issues, we perform extensive GEMC simulations for Lennard-Jones particles in the near-critical region varying the system size, the overall system density, and the cutoff distance. For a system with N = 5500 particles, potential truncation at 8σ and analytical tail corrections, an extrapolation of GEMC simulation data at temperatures in the range from 1.27 to 1.305 yields Tc = 1.3128 ± 0.0016, ρc = 0.316 ± 0.004, and pc = 0.1274 ± 0.0013 in excellent agreement with the thermodynamic limit determined by Potoff and Panagiotopoulos [J. Chem. Phys. 109, 10914 (1998)] using grand canonical Monte Carlo simulations and finite-size scaling. Critical properties estimated using GEMC simulations with different overall system densities (0.296 ≤ ρt ≤ 0.336) agree to within the statistical uncertainties. For simulations with tail corrections, data obtained using rcut = 3.5σ yield Tc and pc that are higher by 0.2% and 1.4% than simulations with rcut = 5 and 8σ but still with overlapping 95% confidence intervals. In contrast, GEMC simulations with a truncated and shifted potential show that rcut = 8σ is insufficient to obtain accurate results. Additional GEMC simulations for hard-core square-well particles with various ranges of the
Dinpajooh, Mohammadhasan; Bai, Peng; Allan, Douglas A.; Siepmann, J. Ilja
2015-09-21
Since the seminal paper by Panagiotopoulos [Mol. Phys. 61, 813 (1997)], the Gibbs ensemble Monte Carlo (GEMC) method has been the most popular particle-based simulation approach for the computation of vapor–liquid phase equilibria. However, the validity of GEMC simulations in the near-critical region has been questioned because rigorous finite-size scaling approaches cannot be applied to simulations with fluctuating volume. Valleau [Mol. Simul. 29, 627 (2003)] has argued that GEMC simulations would lead to a spurious overestimation of the critical temperature. More recently, Patel et al. [J. Chem. Phys. 134, 024101 (2011)] opined that the use of analytical tail corrections would be problematic in the near-critical region. To address these issues, we perform extensive GEMC simulations for Lennard-Jones particles in the near-critical region varying the system size, the overall system density, and the cutoff distance. For a system with N = 5500 particles, potential truncation at 8σ and analytical tail corrections, an extrapolation of GEMC simulation data at temperatures in the range from 1.27 to 1.305 yields T{sub c} = 1.3128 ± 0.0016, ρ{sub c} = 0.316 ± 0.004, and p{sub c} = 0.1274 ± 0.0013 in excellent agreement with the thermodynamic limit determined by Potoff and Panagiotopoulos [J. Chem. Phys. 109, 10914 (1998)] using grand canonical Monte Carlo simulations and finite-size scaling. Critical properties estimated using GEMC simulations with different overall system densities (0.296 ≤ ρ{sub t} ≤ 0.336) agree to within the statistical uncertainties. For simulations with tail corrections, data obtained using r{sub cut} = 3.5σ yield T{sub c} and p{sub c} that are higher by 0.2% and 1.4% than simulations with r{sub cut} = 5 and 8σ but still with overlapping 95% confidence intervals. In contrast, GEMC simulations with a truncated and shifted potential show that r{sub cut} = 8σ is insufficient to obtain accurate results. Additional GEMC simulations for hard
Quantum states of confined hydrogen plasma species: Monte Carlo calculations
NASA Astrophysics Data System (ADS)
Micca Longo, G.; Longo, S.; Giordano, D.
2015-12-01
The diffusion Monte Carlo method with symmetry-based state selection is used to calculate the quantum energy states of \\text{H}2+ confined into potential barriers of atomic dimensions (a model for these ions in solids). Special solutions are employed, permitting one to obtain satisfactory results with rather simple native code. As a test case, {{}2}{{\\Pi}u} and {{}2}{{\\Pi}g} states of \\text{H}2+ ions under spherical confinement are considered. The results are interpreted using the correlation of \\text{H}2+ states to atomic orbitals of H atoms lying on the confining surface and perturbation calculations. The method is straightforwardly applied to cavities of any shape and different hydrogen plasma species (at least one-electron ones, including H) for future studies with real crystal symmetries.
Itinerant scenario for Fe pnictides: Comparison with quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Chubukov, Andrey V.; Xing, Rui-Qi
2016-04-01
Recent applications of quantum Monte Carlo (QMC) technique to Fe-based superconductors opened a way to directly verify the applicability of the itinerant scenario for these systems. Fe-based superconductors undergo various instabilities upon lowering temperature (magnetism, superconductivity, nematicity/orbital order), and one can check whether the hierarchy of instabilities obtained within the itinerant approach is the same as in unbiased QMC simulations. In a recent paper [arXiv:1512.08523] the authors considered the simplest two-band model with interaction tailored to favor orbital order. The type of the orbital order found in QMC is different from the one found in earlier itinerant analysis. We report the results of our calculations within the itinerant scenario and argue that they are in perfect agreement with QMC.
Quantum Monte Carlo Calculations in Solids with Downfolded Hamiltonians
NASA Astrophysics Data System (ADS)
Ma, Fengjie; Purwanto, Wirawan; Zhang, Shiwei; Krakauer, Henry
2015-06-01
We present a combination of a downfolding many-body approach with auxiliary-field quantum Monte Carlo (AFQMC) calculations for extended systems. Many-body calculations operate on a simpler Hamiltonian which retains material-specific properties. The Hamiltonian is systematically improvable and allows one to dial, in principle, between the simplest model and the original Hamiltonian. As a by-product, pseudopotential errors are essentially eliminated using frozen orbitals constructed adaptively from the solid environment. The computational cost of the many-body calculation is dramatically reduced without sacrificing accuracy. Excellent accuracy is achieved for a range of solids, including semiconductors, ionic insulators, and metals. We apply the method to calculate the equation of state of cubic BN under ultrahigh pressure, and determine the spin gap in NiO, a challenging prototypical material with strong electron correlation effects.
Quantum Monte Carlo study of bilayer ionic Hubbard model
NASA Astrophysics Data System (ADS)
Jiang, Mi
The interaction-driven insulator-to-metal transition has been reported in the ionic Hubbard model (IHM) for intermediate interaction U, which poses fundamental interest in the correlated electronic systems. Here we use determinant quantum Monte Carlo to study the interplay of interlayer hybridization V and two types of intralayer staggered potentials: one with the same (in-phase) and the other with a π-phase shift (anti-phase) potential in two layers termed as ``bilayer ionic Hubbard model''. We demonstrate that the interaction-driven Insulator-Metal transition extends to bilayer IHM with finite V for both types of staggered potentials. Besides, the system with in-phase potential is prone to metallic phase with turning on interlayer hybridization while that with anti-phase potential tends to insulators with stronger charge density order. The author thanks CSCS, Lugano, Switzerland for computing facilities.
Quantum Monte Carlo study of bilayer ionic Hubbard model
NASA Astrophysics Data System (ADS)
Jiang, M.; Schulthess, T. C.
2016-04-01
The interaction-driven insulator-to-metal transition has been reported in the ionic Hubbard model (IHM) for moderate interaction U , while its metallic phase only occupies a narrow region in the phase diagram. To explore the enlargement of the metallic regime, we extend the ionic Hubbard model to two coupled layers and study the interplay of interlayer hybridization V and two types of intralayer staggered potentials Δ : one with the same (in-phase) and the other with a π -phase shift (antiphase) potential between layers. Our determinant quantum Monte Carlo (DQMC) simulations at lowest accessible temperatures demonstrate that the interaction-driven metallic phase between Mott and band insulators expands in the Δ -V phase diagram of bilayer IHM only for in-phase ionic potentials; while antiphase potential always induces an insulator with charge density order. This implies possible further extension of the ionic Hubbard model from the bilayer case here to a realistic three-dimensional model.
Continuous-time quantum Monte Carlo using worm sampling
NASA Astrophysics Data System (ADS)
Gunacker, P.; Wallerberger, M.; Gull, E.; Hausoel, A.; Sangiovanni, G.; Held, K.
2015-10-01
We present a worm sampling method for calculating one- and two-particle Green's functions using continuous-time quantum Monte Carlo simulations in the hybridization expansion (CT-HYB). Instead of measuring Green's functions by removing hybridization lines from partition function configurations, as in conventional CT-HYB, the worm algorithm directly samples the Green's function. We show that worm sampling is necessary to obtain general two-particle Green's functions which are not of density-density type and that it improves the sampling efficiency when approaching the atomic limit. Such two-particle Green's functions are needed to compute off-diagonal elements of susceptibilities and occur in diagrammatic extensions of the dynamical mean-field theory and in efficient estimators for the single-particle self-energy.
Quantum Monte Carlo study of magnetism in the Lieb Lattice
NASA Astrophysics Data System (ADS)
Costa, Natanael; Santos, Tiago; Paiva, Thereza; Dos Santos, Raimundo; Scalettar, Richard
The Hubbard model on the `Lieb lattice' provides an important example of how flat band systems may lead to ferromagnetism: at half filling Lieb proved that a ferrimagnetic ground state can be achieved. Since a rigorous proof that long range order does indeed emerge is still lacking, here we report Determinant Quantum Monte Carlo (DQMC) simulations for this model. We found that the spin correlation between nearest neighbors are always antiferromagnetic, and that for small U ferromagnetic long range order does set in in the ground state. However, spatial spin correlations weaken as U is increased, and we established that long range order is suppressed above Uc ~ 4 . 5 . We obtain the dependence of the magnetization with the on-site repulsion U, and show that it displays a maximum at U ~ 3 . The behavior of the compressibility and of the double occupancy across this transition is also discussed. Also at Department of Physics, UC Davis.
Quantum Monte Carlo calculations of A=8 nuclei
Wiringa, R. B.; Pieper, Steven C.; Carlson, J.; Pandharipande, V. R.
2000-07-01
We report quantum Monte Carlo calculations of ground and low-lying excited states for A=8 nuclei using a realistic Hamiltonian containing the Argonne v{sub 18} two-nucleon and Urbana IX three-nucleon potentials. The calculations begin with correlated eight-body wave functions that have a filled {alpha}-like core and four p-shell nucleons LS coupled to the appropriate (J{sup {pi}};T) quantum numbers for the state of interest. After optimization, these variational wave functions are used as input to a Green's function Monte Carlo calculation made with a new constrained path algorithm. We find that the Hamiltonian produces a {sup 8}Be ground state that is within 2 MeV of the experimental resonance, but the other eight-body energies are progressively worse as the neutron-proton asymmetry increases. The {sup 8}Li ground state is stable against breakup into subclusters, but the {sup 8}He ground state is not. The excited state spectra are in fair agreement with experiment, with both the single-particle behavior of {sup 8}He and {sup 8}Li and the collective rotational behavior of {sup 8}Be being reproduced. We also examine energy differences in the T=1,2 isomultiplets and isospin-mixing matrix elements in the excited states of {sup 8}Be. Finally, we present densities, momentum distributions, and studies of the intrinsic shapes of these nuclei, with {sup 8}Be exhibiting a definite 2{alpha} cluster structure. (c) 2000 The American Physical Society.
Fast evaluation of multideterminant wavefunctions in quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Morales, Miguel A.; Clark, Bryan K.; McMinis, Jeremy; Kim, Jeongnim; Scuseria, Gustavo
2011-03-01
Quantum Monte Carlo (QMC) methods such as variational and 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, 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 talk we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient, and easily parallelized. The algorithm computes the multi determinant ratios of a series of particle hole excitations in time O(n 2) + O(n s n)+O(n e) where n, n s and n e are the number of particles, single particle excitations, and total number of excitations, respectively. This is accomplished by producing a (relatively) compact table that contains all the information required to read off the excitation ratios. In addition we describe how to compute the gradients and laplacians of these multi determinant terms. This work was performed under the auspices of: the US DOE by LLNL under Contract DE-AC52-07NA27344, the US DOE under Contract DOE-DE-FG05-08OR23336 and by NSF under No.0904572.
2015-01-01
The penta-2,4-dieniminium cation (PSB3) displays similar ground state and first excited state potential energy features as those of the retinal protonated Schiff base (RPSB) chromophore in rhodopsin. Recently, PSB3 has been used to benchmark several electronic structure methods, including highly correlated multireference wave function approaches, highlighting the necessity to accurately describe the electronic correlation in order to obtain reliable properties even along the ground state (thermal) isomerization paths. In this work, we apply two quantum Monte Carlo approaches, the variational Monte Carlo and the lattice regularized diffusion Monte Carlo, to study the energetics and electronic properties of PSB3 along representative minimum energy paths and scans related to its thermal cis–trans isomerization. Quantum Monte Carlo is used in combination with the Jastrow antisymmetrized geminal power ansatz, which guarantees an accurate and balanced description of the static electronic correlation thanks to the multiconfigurational nature of the antisymmetrized geminal power term, and of the dynamical correlation, due to the presence of the Jastrow factor explicitly depending on electron–electron distances. Along the two ground state isomerization minimum energy paths of PSB3, CASSCF calculations yield wave functions having either charge transfer or diradical character in proximity of the two transition state configurations. Here, we observe that at the quantum Monte Carlo level of theory, only the transition state with charge transfer character can be located. The conical intersection, which becomes highly sloped, is observed only if the path connecting the two original CASSCF transition states is extended beyond the diradical one, namely by increasing the bond-length-alternation (BLA). These findings are in good agreement with the results obtained by MRCISD+Q calculations, and they demonstrate the importance of having an accurate description of the static and
Zen, Andrea; Coccia, Emanuele; Gozem, Samer; Olivucci, Massimo; Guidoni, Leonardo
2015-03-10
The penta-2,4-dieniminium cation (PSB3) displays similar ground state and first excited state potential energy features as those of the retinal protonated Schiff base (RPSB) chromophore in rhodopsin. Recently, PSB3 has been used to benchmark several electronic structure methods, including highly correlated multireference wave function approaches, highlighting the necessity to accurately describe the electronic correlation in order to obtain reliable properties even along the ground state (thermal) isomerization paths. In this work, we apply two quantum Monte Carlo approaches, the variational Monte Carlo and the lattice regularized diffusion Monte Carlo, to study the energetics and electronic properties of PSB3 along representative minimum energy paths and scans related to its thermal cis–trans isomerization. Quantum Monte Carlo is used in combination with the Jastrow antisymmetrized geminal power ansatz, which guarantees an accurate and balanced description of the static electronic correlation thanks to the multiconfigurational nature of the antisymmetrized geminal power term, and of the dynamical correlation, due to the presence of the Jastrow factor explicitly depending on electron–electron distances. Along the two ground state isomerization minimum energy paths of PSB3, CASSCF calculations yield wave functions having either charge transfer or diradical character in proximity of the two transition state configurations. Here, we observe that at the quantum Monte Carlo level of theory, only the transition state with charge transfer character can be located. The conical intersection, which becomes highly sloped, is observed only if the path connecting the two original CASSCF transition states is extended beyond the diradical one, namely by increasing the bond-length-alternation (BLA). These findings are in good agreement with the results obtained by MRCISD+Q calculations, and they demonstrate the importance of having an accurate description of the static and
Nielsen, Jens; D’Avezac, Mayeul; Hetherington, James; Stamatakis, Michail
2013-12-14
Ab initio kinetic Monte Carlo (KMC) simulations have been successfully applied for over two decades to elucidate the underlying physico-chemical phenomena on the surfaces of heterogeneous catalysts. These simulations necessitate detailed knowledge of the kinetics of elementary reactions constituting the reaction mechanism, and the energetics of the species participating in the chemistry. The information about the energetics is encoded in the formation energies of gas and surface-bound species, and the lateral interactions between adsorbates on the catalytic surface, which can be modeled at different levels of detail. The majority of previous works accounted for only pairwise-additive first nearest-neighbor interactions. More recently, cluster-expansion Hamiltonians incorporating long-range interactions and many-body terms have been used for detailed estimations of catalytic rate [C. Wu, D. J. Schmidt, C. Wolverton, and W. F. Schneider, J. Catal. 286, 88 (2012)]. In view of the increasing interest in accurate predictions of catalytic performance, there is a need for general-purpose KMC approaches incorporating detailed cluster expansion models for the adlayer energetics. We have addressed this need by building on the previously introduced graph-theoretical KMC framework, and we have developed Zacros, a FORTRAN2003 KMC package for simulating catalytic chemistries. To tackle the high computational cost in the presence of long-range interactions we introduce parallelization with OpenMP. We further benchmark our framework by simulating a KMC analogue of the NO oxidation system established by Schneider and co-workers [J. Catal. 286, 88 (2012)]. We show that taking into account only first nearest-neighbor interactions may lead to large errors in the prediction of the catalytic rate, whereas for accurate estimates thereof, one needs to include long-range terms in the cluster expansion.
Aneesu-Rahman Prize Lecture: The ``sign problem'' in Quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Ceperley, D. M.
1998-03-01
Quantum simulation methods have been quite successful in giving exact results for certain systems, primarily bosons(Ceperley, D.M. , Rev. Mod. Phys. 67), 279 (1995).. Use of the same techniques in general quantum systems leads to the so-called ``sign problem''; the results are correct but the methods are very inefficient. There are two important questions to ask of a proposed method. Given enough computer time can arbitrarily accurate results be obtained? If so, how long does it take to achieve a given error? There are several methods (released-node or transient estimate) that are exact; the difficulty is in finding a method which also scales well with the number of quantum degrees of freedom. Exact methods, in general, scale exponentially with the number of fermions and in the inverse temperature (or accuracy). At root, the fact that wavefunction is complex or changes sign, gives rise to the poor scaling and the ``sign problem.'' It is not the fermion nature of the system, per se, that causes the difficulty. The desired state is not the absolute ground state. Methods which cancel random walks from positive and negative regions have also been limited to quite small systems because they scale poorly. There are a variety of approximate simulation methods which do scale well, such as variational Monte Carlo, and a variety of fixed-node methods (restricted Path Integral Monte Carlo at non-zero temperature and constrained path methods for lattice models) which fix only boundary conditions not the sampling function. For many systems, the variational and fixed-node methods can be very accurate. The lecture notes and references are on my group's homepage.
Generalized directed loop method for quantum Monte Carlo simulations.
Alet, Fabien; Wessel, Stefan; Troyer, Matthias
2005-03-01
Efficient quantum Monte Carlo update schemes called directed loops have recently been proposed, which improve the efficiency of simulations of quantum lattice models. We propose to generalize the detailed balance equations at the local level during the loop construction by accounting for the matrix elements of the operators associated with open world-line segments. Using linear programming techniques to solve the generalized equations, we look for optimal construction schemes for directed loops. This also allows for an extension of the directed loop scheme to general lattice models, such as high-spin or bosonic models. The resulting algorithms are bounce free in larger regions of parameter space than the original directed loop algorithm. The generalized directed loop method is applied to the magnetization process of spin chains in order to compare its efficiency to that of previous directed loop schemes. In contrast to general expectations, we find that minimizing bounces alone does not always lead to more efficient algorithms in terms of autocorrelations of physical observables, because of the nonuniqueness of the bounce-free solutions. We therefore propose different general strategies to further minimize autocorrelations, which can be used as supplementary requirements in any directed loop scheme. We show by calculating autocorrelation times for different observables that such strategies indeed lead to improved efficiency; however, we find that the optimal strategy depends not only on the model parameters but also on the observable of interest. PMID:15903632
Quantum Monte Carlo simulation of spin-polarized H
Markic, L. Vranjes; Boronat, J.; Casulleras, J.
2007-02-01
The ground-state properties of spin polarized hydrogen H{down_arrow} are obtained by means of diffusion Monte Carlo calculations. Using the most accurate to date ab initio H{down_arrow}-H{down_arrow} interatomic potential we have studied its gas phase, from the very dilute regime until densities above its freezing point. At very small densities, the equation of state of the gas is very well described in terms of the gas parameter {rho}a{sup 3}, with a the s-wave scattering length. The solid phase has also been studied up to high pressures. The gas-solid phase transition occurs at a pressure of 173 bar, a much higher value than suggested by previous approximate descriptions.
Polynomial-time-scaling quantum dynamics with time-dependent quantum Monte Carlo.
Christov, Ivan P
2009-05-21
Here we study the dynamics of many-body quantum systems using the time-dependent quantum Monte Carlo method where the evolution is described by ensembles of particles and guide waves. The exponential time scaling inherent to the quantum many-body problem is reduced to polynomial-time computation by solving concurrently a set of coupled Schrodinger equations for the guide waves in physical space and a set of first-order equations for the Monte Carlo walkers. We use effective potentials to account for the local and nonlocal quantum correlations in time-varying fields, where for fermionic states an exchange "hole" is introduced explicitly through screened Coulomb potentials. The walker distributions for the ground states of para- and ortho-helium reproduce well the statistical properties, such as the electron-pair density function, of the real atoms. Our predictions for the dipole response and the ionization of an atom exposed to strong ultrashort optical pulse are in good agreement with the exact results. PMID:19391581
Quantum Monte Carlo with very large multideterminant wavefunctions.
Scemama, Anthony; Applencourt, Thomas; Giner, Emmanuel; Caffarel, Michel
2016-07-01
An algorithm to compute efficiently the first two derivatives of (very) large multideterminant wavefunctions for quantum Monte Carlo calculations is presented. The calculation of determinants and their derivatives is performed using the Sherman-Morrison formula for updating the inverse Slater matrix. An improved implementation based on the reduction of the number of column substitutions and on a very efficient implementation of the calculation of the scalar products involved is presented. It is emphasized that multideterminant expansions contain in general a large number of identical spin-specific determinants: for typical configuration interaction-type wavefunctions the number of unique spin-specific determinants Ndetσ ( σ=↑,↓) with a non-negligible weight in the expansion is of order O(Ndet). We show that a careful implementation of the calculation of the Ndet -dependent contributions can make this step negligible enough so that in practice the algorithm scales as the total number of unique spin-specific determinants, Ndet↑+Ndet↓, over a wide range of total number of determinants (here, Ndet up to about one million), thus greatly reducing the total computational cost. Finally, a new truncation scheme for the multideterminant expansion is proposed so that larger expansions can be considered without increasing the computational time. The algorithm is illustrated with all-electron fixed-node diffusion Monte Carlo calculations of the total energy of the chlorine atom. Calculations using a trial wavefunction including about 750,000 determinants with a computational increase of ∼400 compared to a single-determinant calculation are shown to be feasible. © 2016 Wiley Periodicals, Inc. PMID:27302337
Quantum Monte Carlo Algorithms for Diagrammatic Vibrational Structure Calculations
NASA Astrophysics Data System (ADS)
Hermes, Matthew; Hirata, So
2015-06-01
Convergent hierarchies of theories for calculating many-body vibrational ground and excited-state wave functions, such as Møller-Plesset perturbation theory or coupled cluster theory, tend to rely on matrix-algebraic manipulations of large, high-dimensional arrays of anharmonic force constants, tasks which require large amounts of computer storage space and which are very difficult to implement in a parallel-scalable fashion. On the other hand, existing quantum Monte Carlo (QMC) methods for vibrational wave functions tend to lack robust techniques for obtaining excited-state energies, especially for large systems. By exploiting analytical identities for matrix elements of position operators in a harmonic oscillator basis, we have developed stochastic implementations of the size-extensive vibrational self-consistent field (MC-XVSCF) and size-extensive vibrational Møller-Plesset second-order perturbation (MC-XVMP2) theories which do not require storing the potential energy surface (PES). The programmable equations of MC-XVSCF and MC-XVMP2 take the form of a small number of high-dimensional integrals evaluated using Metropolis Monte Carlo techniques. The associated integrands require independent evaluations of only the value, not the derivatives, of the PES at many points, a task which is trivial to parallelize. However, unlike existing vibrational QMC methods, MC-XVSCF and MC-XVMP2 can calculate anharmonic frequencies directly, rather than as a small difference between two noisy total energies, and do not require user-selected coordinates or nodal surfaces. MC-XVSCF and MC-XVMP2 can also directly sample the PES in a given approximation without analytical or grid-based approximations, enabling us to quantify the errors induced by such approximations.
Accurate energies of the He atom with undergraduate quantum mechanics
NASA Astrophysics Data System (ADS)
Massé, Robert C.; Walker, Thad G.
2015-08-01
Estimating the energies and splitting of the 1s2s singlet and triplet states of helium is a classic exercise in quantum perturbation theory but yields only qualitatively correct results. Using a six-line computer program, the 1s2s energies calculated by matrix diagonalization using a seven-state basis improve the results to 0.4% error or better. This is an effective and practical illustration of the quantitative power of quantum mechanics, at a level accessible to undergraduate students.
Auxiliary-field quantum Monte Carlo calculations of the molybdenum dimer
NASA Astrophysics Data System (ADS)
Purwanto, Wirawan; Zhang, Shiwei; Krakauer, Henry
2016-06-01
Chemical accuracy is difficult to achieve for systems with transition metal atoms. Third row transition metal atoms are particularly challenging due to strong electron-electron correlation in localized d-orbitals. The Cr2 molecule is an outstanding example, which we previously treated with highly accurate auxiliary-field quantum Monte Carlo (AFQMC) calculations [W. Purwanto et al., J. Chem. Phys. 142, 064302 (2015)]. Somewhat surprisingly, computational description of the isoelectronic Mo2 dimer has also, to date, been scattered and less than satisfactory. We present high-level theoretical benchmarks of the Mo2 singlet ground state (X1Σg+) and first triplet excited state (a3Σu+), using the phaseless AFQMC calculations. Extrapolation to the complete basis set limit is performed. Excellent agreement with experimental spectroscopic constants is obtained. We also present a comparison of the correlation effects in Cr2 and Mo2.
An auxiliary-field quantum Monte Carlo study of the chromium dimer
Purwanto, Wirawan Zhang, Shiwei; Krakauer, Henry
2015-02-14
The chromium dimer (Cr{sub 2}) presents an outstanding challenge for many-body electronic structure methods. Its complicated nature of binding, with a formal sextuple bond and an unusual potential energy curve (PEC), is emblematic of the competing tendencies and delicate balance found in many strongly correlated materials. We present an accurate calculation of the PEC and ground state properties of Cr{sub 2}, using the auxiliary-field quantum Monte Carlo (AFQMC) method. Unconstrained, exact AFQMC calculations are first carried out for a medium-sized but realistic basis set. Elimination of the remaining finite-basis errors and extrapolation to the complete basis set limit are then achieved with a combination of phaseless and exact AFQMC calculations. Final results for the PEC and spectroscopic constants are in excellent agreement with experiment.
Recent developments in auxiliary-field quantum Monte Carlo methods for cold atoms
NASA Astrophysics Data System (ADS)
Shi, Hao; Rosenberg, Peter; Vitali, Ettore; Chiesa, Simone; Zhang, Shiwei
Exact calculations are performed on the two-dimensional strongly interacting, unpolarized, uniform Fermi gas with a zero-range attractive interaction. We describe recent advances in auxiliary-field quantum Monte Carlo techniques, which eliminate an infinite variance problem in the standard algorithm, and improve both acceptance ratio and efficiency. The new methods enable calculations on large enough lattices to reliably compute ground-state properties in the thermodynamic limit. An equation of state is obtained, with a parametrization provided, which can serve as a benchmark and allow accurate comparisons with experiments. The pressure, contact parameter, condensate fraction, and pairing gap will be presented. The same methods are also applied to obtain exact results on the two-dimensional strongly interacting Fermi gas in the presence of Rashba spin-orbit (SOC), providing insights on the interplay between pairing and SOC. Supported by NSF, DOE, and the Simons Foundation.
Constrained-path quantum Monte Carlo approach for non-yrast states within the shell model
NASA Astrophysics Data System (ADS)
Bonnard, J.; Juillet, O.
2016-04-01
The present paper intends to present an extension of the constrained-path quantum Monte Carlo approach allowing to reconstruct non-yrast states in order to reach the complete spectroscopy of nuclei within the interacting shell model. As in the yrast case studied in a previous work, the formalism involves a variational symmetry-restored wave function assuming two central roles. First, it guides the underlying Brownian motion to improve the efficiency of the sampling. Second, it constrains the stochastic paths according to the phaseless approximation to control sign or phase problems that usually plague fermionic QMC simulations. Proof-of-principle results in the sd valence space are reported. They prove the ability of the scheme to offer remarkably accurate binding energies for both even- and odd-mass nuclei irrespective of the considered interaction.
Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures.
Drummond, N D; Monserrat, Bartomeu; Lloyd-Williams, Jonathan H; López Ríos, P; Pickard, Chris J; Needs, R J
2015-01-01
Establishing the phase diagram of hydrogen is a major challenge for experimental and theoretical physics. Experiment alone cannot establish the atomic structure of solid hydrogen at high pressure, because hydrogen scatters X-rays only weakly. Instead, our understanding of the atomic structure is largely based on density functional theory (DFT). By comparing Raman spectra for low-energy structures found in DFT searches with experimental spectra, candidate atomic structures have been identified for each experimentally observed phase. Unfortunately, DFT predicts a metallic structure to be energetically favoured at a broad range of pressures up to 400 GPa, where it is known experimentally that hydrogen is non-metallic. Here we show that more advanced theoretical methods (diffusion quantum Monte Carlo calculations) find the metallic structure to be uncompetitive, and predict a phase diagram in reasonable agreement with experiment. This greatly strengthens the claim that the candidate atomic structures accurately model the experimentally observed phases. PMID:26215251
Clay, Raymond C.; Mcminis, Jeremy; McMahon, Jeffrey M.; Pierleoni, Carlo; Ceperley, David M.; Morales, Miguel A.
2014-05-01
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 outperform local density approximation and Perdew-Burke-Ernzerhof. 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.
Variational quantum Monte Carlo calculation of electronic and structural properties of crystals
Louie, S.G.
1989-09-01
Calculation of the electronic and structural properties of solids using a variational quantum Monte Carlo nonlocal pseudopotential approach is described. Ionization potentials and electron affinities for atoms, and binding energies and structural properties for crystals are found to be in very good agreement with experiment. The approach employs a correlated many-electron wavefunction of the Jastrow-Slater form and the exact Coulomb interaction between valence electrons. One- and two-body terms in the Jastrow factor are used and found necessary for an accurate description of the electron-electron energy for the systems considered. The method has further been applied to compute various single-particle properties for solids including the single-particle orbital occupancy, electron pair correlation functions, and quasiparticle excitation energies. 23 refs., 3 figs., 3 tabs.
A fast and efficient algorithm for Slater determinant updates in quantum Monte Carlo simulations
Nukala, Phani K. V. V.; Kent, P. R. C.
2009-05-28
We present an efficient low-rank updating algorithm for updating the trial wave functions used in quantum Monte Carlo (QMC) simulations. The algorithm is based on low-rank updating of the Slater determinants. In particular, the computational complexity of the algorithm is O(kN) during the kth step compared to traditional algorithms that require O(N{sup 2}) computations, where N is the system size. For single determinant trial wave functions the new algorithm is faster than the traditional O(N{sup 2}) Sherman-Morrison algorithm for up to O(N) updates. For multideterminant configuration-interaction-type trial wave functions of M+1 determinants, the new algorithm is significantly more efficient, saving both O(MN{sup 2}) work and O(MN{sup 2}) storage. The algorithm enables more accurate and significantly more efficient QMC calculations using configuration-interaction-type wave functions.
A fast and efficient algorithm for Slater determinant updates in quantum Monte Carlo simulations
NASA Astrophysics Data System (ADS)
Nukala, Phani K. V. V.; Kent, P. R. C.
2009-05-01
We present an efficient low-rank updating algorithm for updating the trial wave functions used in quantum Monte Carlo (QMC) simulations. The algorithm is based on low-rank updating of the Slater determinants. In particular, the computational complexity of the algorithm is O(kN) during the kth step compared to traditional algorithms that require O(N2) computations, where N is the system size. For single determinant trial wave functions the new algorithm is faster than the traditional O(N2) Sherman-Morrison algorithm for up to O(N ) updates. For multideterminant configuration-interaction-type trial wave functions of M +1 determinants, the new algorithm is significantly more efficient, saving both O(MN2) work and O(MN2) storage. The algorithm enables more accurate and significantly more efficient QMC calculations using configuration-interaction-type wave functions.
A fast and efficient algorithm for Slater determinant updates in quantum Monte Carlo simulations.
Nukala, Phani K V V; Kent, P R C
2009-05-28
We present an efficient low-rank updating algorithm for updating the trial wave functions used in quantum Monte Carlo (QMC) simulations. The algorithm is based on low-rank updating of the Slater determinants. In particular, the computational complexity of the algorithm is O(kN) during the kth step compared to traditional algorithms that require O(N(2)) computations, where N is the system size. For single determinant trial wave functions the new algorithm is faster than the traditional O(N(2)) Sherman-Morrison algorithm for up to O(N) updates. For multideterminant configuration-interaction-type trial wave functions of M+1 determinants, the new algorithm is significantly more efficient, saving both O(MN(2)) work and O(MN(2)) storage. The algorithm enables more accurate and significantly more efficient QMC calculations using configuration-interaction-type wave functions. PMID:19485435
A Fast and efficient Algorithm for Slater Determinant Updates in Quantum Monte Carlo Simulations
Nukala, Phani K; Kent, Paul R
2009-01-01
We present an efficient low-rank updating algorithm for updating the trial wavefunctions used in Quantum Monte Carlo (QMC) simulations. The algorithm is based on low-rank updating of the Slater determinants. In particular, the computational complexity of the algorithm is $\\mathcal{O}(k N)$ during the $k$-th step compared with traditional algorithms that require $\\mathcal{O}(N^2)$ computations, where $N$ is the system size. For single determinant trial wavefunctions the new algorithm is faster than the traditional $\\mathcal{O}(N^2)$ Sherman-Morrison algorithm for up to $\\mathcal{O}(N)$ updates. For multideterminant configuration-interaction type trial wavefunctions of $M+1$ determinants, the new algorithm is significantly more efficient, saving both $\\mathcal{O}(MN^2)$ work and $\\mathcal{O}(MN^2)$ storage. The algorithm enables more accurate and significantly more efficient QMC calculations using configuration interaction type wavefunctions.
Quantum Monte Carlo study of charged transition-metal organometallic cluster systems
NASA Astrophysics Data System (ADS)
Tokar, Kamil; Derian, Rene; Stich, Ivan
2015-03-01
Using accurate fixed-node quantum Monte Carlo (QMC) methods we study 1D clusters formed by transition metal atoms separated by benzene molecules (TMBz), both positively and negatively charged. TMBz are among the most important π-bonded organometallics, which, however, often require charged states for experimental studies. We have performed a systematic study of ground-sate spin multiplets, ionization potentials, electron affinities, and dissociation energies of vanadium-benzene cationic and anionic half- and full-sandwiches. By comparison of QMC and DFT results, we find a very strong impact of electronic correlation on properties of these systems, such as dissociation energies, where ~1 eV energy corrections are found. In particular, the anions are unstable at the DFT level and are stabilized only at the QMC level after sophisticated optimization of the trial wavefunction. Supported by APVV-0207-11 and VEGA (2/0007/12) projects.
Energy density matrix formalism for interacting quantum systems: a quantum Monte Carlo study
Krogel, Jaron T; Kim, Jeongnim; Reboredo, Fernando A
2014-01-01
We develop an energy density matrix that parallels the one-body reduced density matrix (1RDM) for many-body quantum systems. Just as the density matrix gives access to the number density and occupation numbers, the energy density matrix yields the energy density and orbital occupation energies. The eigenvectors of the matrix provide a natural orbital partitioning of the energy density while the eigenvalues comprise a single particle energy spectrum obeying a total energy sum rule. For mean-field systems the energy density matrix recovers the exact spectrum. When correlation becomes important, the occupation energies resemble quasiparticle energies in some respects. We explore the occupation energy spectrum for the finite 3D homogeneous electron gas in the metallic regime and an isolated oxygen atom with ground state quantum Monte Carlo techniques imple- mented in the QMCPACK simulation code. The occupation energy spectrum for the homogeneous electron gas can be described by an effective mass below the Fermi level. Above the Fermi level evanescent behavior in the occupation energies is observed in similar fashion to the occupation numbers of the 1RDM. A direct comparison with total energy differences demonstrates a quantita- tive connection between the occupation energies and electron addition and removal energies for the electron gas. For the oxygen atom, the association between the ground state occupation energies and particle addition and removal energies becomes only qualitative. The energy density matrix provides a new avenue for describing energetics with quantum Monte Carlo methods which have traditionally been limited to total energies.
NASA Astrophysics Data System (ADS)
Clay, Raymond C.; Morales, Miguel A.
2015-06-01
Multideterminant wavefunctions, while having a long history in quantum chemistry, are increasingly being used in highly accurate quantum Monte Carlo calculations. Since the accuracy of QMC is ultimately limited by the quality of the trial wavefunction, multi-Slater determinants wavefunctions offer an attractive alternative to Slater-Jastrow and more sophisticated wavefunction ansatz for several reasons. They can be efficiently calculated, straightforwardly optimized, and systematically improved by increasing the number of included determinants. In spite of their potential, however, the convergence properties of multi-Slater determinant wavefunctions with respect to orbital set choice and excited determinant selection are poorly understood, which hinders the application of these wavefunctions to large systems and solids. In this paper, by performing QMC calculations on the equilibrium and stretched carbon dimer, we find that convergence of the recovered correlation energy with respect to number of determinants can depend quite strongly on basis set and determinant selection methods, especially where there is strong correlation. We demonstrate that properly chosen orbital sets and determinant selection techniques from quantum chemistry methods can dramatically reduce the required number of determinants (and thus the computational cost) to reach a given accuracy, which we argue shows clear need for an automatic QMC-only method for selecting determinants and generating optimal orbital sets.
Clay, Raymond C.; Morales, Miguel A.
2015-06-21
Multideterminant wavefunctions, while having a long history in quantum chemistry, are increasingly being used in highly accurate quantum Monte Carlo calculations. Since the accuracy of QMC is ultimately limited by the quality of the trial wavefunction, multi-Slater determinants wavefunctions offer an attractive alternative to Slater-Jastrow and more sophisticated wavefunction ansatz for several reasons. They can be efficiently calculated, straightforwardly optimized, and systematically improved by increasing the number of included determinants. In spite of their potential, however, the convergence properties of multi-Slater determinant wavefunctions with respect to orbital set choice and excited determinant selection are poorly understood, which hinders the application of these wavefunctions to large systems and solids. In this paper, by performing QMC calculations on the equilibrium and stretched carbon dimer, we find that convergence of the recovered correlation energy with respect to number of determinants can depend quite strongly on basis set and determinant selection methods, especially where there is strong correlation. We demonstrate that properly chosen orbital sets and determinant selection techniques from quantum chemistry methods can dramatically reduce the required number of determinants (and thus the computational cost) to reach a given accuracy, which we argue shows clear need for an automatic QMC-only method for selecting determinants and generating optimal orbital sets.
Clay, Raymond C; Morales, Miguel A
2015-06-21
Multideterminant wavefunctions, while having a long history in quantum chemistry, are increasingly being used in highly accurate quantum Monte Carlo calculations. Since the accuracy of QMC is ultimately limited by the quality of the trial wavefunction, multi-Slater determinants wavefunctions offer an attractive alternative to Slater-Jastrow and more sophisticated wavefunction ansatz for several reasons. They can be efficiently calculated, straightforwardly optimized, and systematically improved by increasing the number of included determinants. In spite of their potential, however, the convergence properties of multi-Slater determinant wavefunctions with respect to orbital set choice and excited determinant selection are poorly understood, which hinders the application of these wavefunctions to large systems and solids. In this paper, by performing QMC calculations on the equilibrium and stretched carbon dimer, we find that convergence of the recovered correlation energy with respect to number of determinants can depend quite strongly on basis set and determinant selection methods, especially where there is strong correlation. We demonstrate that properly chosen orbital sets and determinant selection techniques from quantum chemistry methods can dramatically reduce the required number of determinants (and thus the computational cost) to reach a given accuracy, which we argue shows clear need for an automatic QMC-only method for selecting determinants and generating optimal orbital sets. PMID:26093546
Quantum Monte Carlo calculations applied to magnetic molecules
NASA Astrophysics Data System (ADS)
Engelhardt, Larry Paul
In this dissertation, we have implemented a quantum Monte Carlo (QMC) algorithm, and have used it to perform calculations for a variety of finite Heisenberg spin systems. A detailed description of the QMC method has been provided, which is followed by applications of the method to various systems. These applications begin with a detailed analysis of the (calculated) equilibrium magnetization and magnetic susceptibility for a number of Heisenberg Hamiltonians. In particular, we have studied the dependence of these quantities on intrinsic spin s, and have quantified the approach to the classical (s → infinity) limit. These results are not specific to a particular physical system, but are potentially applicable to many systems. We have also analyzed four recently synthesized species of magnetic molecules, each of which is theoretically challenging to the methods that are normally used for such analyses. Using the QMC method, we have distinguished the microscopic (exchange) parameters that describe the interactions in each of these magnetic molecules, and, based upon these parameters, we have made predictions for future experiments. The well-known "negative sign problem" (NSP) can be problematic for QMC calculations. However, for some systems, our analysis was able to proceed despite the NSP. For other systems, this is not the cases, so we have clearly indicated when the NSP is, and is not, insurmountable for these types of calculations.
Quantum Monte Carlo investigations of adsorption energetics on graphene.
Hsing, C R; Wei, C M; Chou, M Y
2012-10-01
We have performed calculations of adsorption energetics on the graphene surface using the state-of-the-art diffusion quantum Monte Carlo method. Two types of configurations are considered in this work: the adsorption of a single O, F, or H atom on the graphene surface and the H-saturated graphene system (graphane). The adsorption energies are compared with those obtained from density functional theory with various exchange-correlation functionals. The results indicate that the approximate exchange-correlation functionals significantly overestimate the binding of O and F atoms on graphene, although the preferred adsorption sites are consistent. The energy errors are much less for atomic hydrogen adsorbed on the surface. We also find that a single O or H atom on graphene has a higher energy than in the molecular state, while the adsorption of a single F atom is preferred over the gas phase. In addition, the energetics of graphane is reported. The calculated equilibrium lattice constant turns out to be larger than that of graphene, at variance with a recent experimental suggestion. PMID:22909778
Quantum Monte Carlo simulations for disordered Bose systems
Trivedi, N.
1992-03-01
Interacting bosons in a random potential can be used to model {sup 3}He adsorbed in porous media, universal aspects of the superconductor-insulator transition in disordered films, and vortices in disordered type II superconductors. We study a model of bosons on a 2D square lattice with a random potential of strength V and on-site repulsion U. We first describe the path integral Monte Carlo algorithm used to simulate this system. The 2D quantum problem (at T=0) gets mapped onto a classical problem of strings or directed polymers moving in 3D with each string representing the world line of a boson. We discuss efficient ways of sampling the polymer configurations as well as the permutations between the bosons. We calculate the superfluid density and the excitation spectrum. Using these results we distinguish between a superfluid, a localized or Bose glass'' insulator with gapless excitations and a Mott insulator with a finite gap to excitations (found only at commensurate densities). We discover novel effects arising from the interpaly between V and U and present preliminary results for the phase diagram at incommensurate and commensurate densities.
Quantum Monte Carlo simulations for disordered Bose systems
Trivedi, N.
1992-03-01
Interacting bosons in a random potential can be used to model {sup 3}He adsorbed in porous media, universal aspects of the superconductor-insulator transition in disordered films, and vortices in disordered type II superconductors. We study a model of bosons on a 2D square lattice with a random potential of strength V and on-site repulsion U. We first describe the path integral Monte Carlo algorithm used to simulate this system. The 2D quantum problem (at T=0) gets mapped onto a classical problem of strings or directed polymers moving in 3D with each string representing the world line of a boson. We discuss efficient ways of sampling the polymer configurations as well as the permutations between the bosons. We calculate the superfluid density and the excitation spectrum. Using these results we distinguish between a superfluid, a localized or ``Bose glass`` insulator with gapless excitations and a Mott insulator with a finite gap to excitations (found only at commensurate densities). We discover novel effects arising from the interpaly between V and U and present preliminary results for the phase diagram at incommensurate and commensurate densities.
Bond breaking with auxiliary-field quantum Monte Carlo.
Al-Saidi, W A; Zhang, Shiwei; Krakauer, Henry
2007-10-14
Bond stretching mimics different levels of electron correlation and provides a challenging test bed for approximate many-body computational methods. Using the recently developed phaseless auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well-studied molecules BH and N(2) and in the H(50) chain. To control the sign/phase problem, the phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. With single Slater determinants from unrestricted Hartree-Fock as trial wave function, the phaseless AF QMC method generally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) method in mapping the potential-energy curve. In both BH and N(2), we also study the use of multiple-determinant trial wave functions from multiconfiguration self-consistent-field calculations. The increase in computational cost versus the gain in statistical and systematic accuracy are examined. With such trial wave functions, excellent results are obtained across the entire region between equilibrium and the dissociation limit. PMID:17935380
Doblhoff-Dier, Katharina; Meyer, Jörg; Hoggan, Philip E; Kroes, Geert-Jan; Wagner, Lucas K
2016-06-14
Transition metals and transition metal compounds are important to catalysis, photochemistry, and many superconducting systems. We study the performance of diffusion Monte Carlo (DMC) applied to transition metal containing dimers (TMCDs) using single-determinant Slater-Jastrow trial wavefunctions and investigate the possible influence of the locality and pseudopotential errors. We find that the locality approximation can introduce nonsystematic errors of up to several tens of kilocalories per mole in the absolute energy of Cu and CuH if Ar or Mg core pseudopotentials (PPs) are used for the 3d transition metal atoms. Even for energy differences such as binding energies, errors due to the locality approximation can be problematic if chemical accuracy is sought. The use of the Ne core PPs developed by Burkatzki et al. (J. Chem. Phys. 2008, 129, 164115), the use of linear energy minimization rather than unreweighted variance minimization for the optimization of the Jastrow function, and the use of large Jastrow parametrizations reduce the locality errors. In the second section of this article, we study the general performance of DMC for 3d TMCDs using a database of binding energies of 20 TMCDs, for which comparatively accurate experimental data is available. Comparing our DMC results to these data for our results that compare best with experiment, we find a mean unsigned error (MUE) of 4.5 kcal/mol. This compares well with the achievable accuracy in CCSDT(2)Q (MUE = 4.6 kcal/mol) and the best all-electron DFT results (MUE = 4.5 kcal/mol) for the same set of systems (Truhlar et al. J. Chem. Theory Comput. 2015, 11, 2036-2052). The mean errors in DMC depend less on the exchange-correlation functionals used to generate the trial wavefunction than the corresponding mean errors in the underlying DFT calculations. Furthermore, the QMC results obtained for each molecule individually vary less with the functionals used. These observations are relevant for systems such as
Quantum dynamics at finite temperature: Time-dependent quantum Monte Carlo study
NASA Astrophysics Data System (ADS)
Christov, Ivan P.
2016-08-01
In this work we investigate the ground state and the dissipative quantum dynamics of interacting charged particles in an external potential at finite temperature. The recently devised time-dependent quantum Monte Carlo (TDQMC) method allows a self-consistent treatment of the system of particles together with bath oscillators first for imaginary-time propagation of Schrödinger type of equations where both the system and the bath converge to their finite temperature ground state, and next for real time calculation where the dissipative dynamics is demonstrated. In that context the application of TDQMC appears as promising alternative to the path-integral related techniques where the real time propagation can be a challenge.
NASA Astrophysics Data System (ADS)
Clay, Raymond C.; Holzmann, Markus; Ceperley, David M.; Morales, Miguel A.
2016-01-01
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 density-functional-theory-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 (QMC) -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 insight 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 we estimate the magnitudes of these effects when possible.
Clay, Raymond C.; Holzmann, Markus; Ceperley, David M.; Morales, Maguel A.
2016-01-19
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 rst principles methods have the potential to provide the accuracy and computational e ciency required for this task, recent benchmarking in hydrogen has shown that achieving this accuracy requires a judicious choice of functional, and a quanti cation 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 atmore » 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 di erent density functionals. We nd 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 di erences exhibited by the major classes of functionals, and estimate the magnitudes of these e ects when possible.« less
Quantum Monte Carlo for electronic structure: Recent developments and applications
Rodriquez, M. M.S.
1995-04-01
Quantum Monte Carlo (QMC) methods have been found to give excellent results when applied to chemical systems. The main goal of the present work is to use QMC to perform electronic structure calculations. In QMC, a Monte Carlo simulation is used to solve the Schroedinger equation, taking advantage of its analogy to a classical diffusion process with branching. In the present work the author focuses on how to extend the usefulness of QMC to more meaningful molecular systems. This study is aimed at questions concerning polyatomic and large atomic number systems. The accuracy of the solution obtained is determined by the accuracy of the trial wave function`s nodal structure. Efforts in the group have given great emphasis to finding optimized wave functions for the QMC calculations. Little work had been done by systematically looking at a family of systems to see how the best wave functions evolve with system size. In this work the author presents a study of trial wave functions for C, CH, C{sub 2}H and C{sub 2}H{sub 2}. The goal is to study how to build wave functions for larger systems by accumulating knowledge from the wave functions of its fragments as well as gaining some knowledge on the usefulness of multi-reference wave functions. In a MC calculation of a heavy atom, for reasonable time steps most moves for core electrons are rejected. For this reason true equilibration is rarely achieved. A method proposed by Batrouni and Reynolds modifies the way the simulation is performed without altering the final steady-state solution. It introduces an acceleration matrix chosen so that all coordinates (i.e., of core and valence electrons) propagate at comparable speeds. A study of the results obtained using their proposed matrix suggests that it may not be the optimum choice. In this work the author has found that the desired mixing of coordinates between core and valence electrons is not achieved when using this matrix. A bibliography of 175 references is included.
NASA Astrophysics Data System (ADS)
Das, Arnab; Chakrabarti, Bikas K.
2008-12-01
Here we discuss the annealing behavior of an infinite-range ±J Ising spin glass in the presence of a transverse field using a zero-temperature quantum Monte Carlo method. Within the simulation scheme, we demonstrate that quantum annealing not only helps finding the ground state of a classical spin glass, but can also help simulating the ground state of a quantum spin glass, in particular, when the transverse field is low, much more efficiently.
Das, Arnab; Chakrabarti, Bikas K
2008-12-01
Here we discuss the annealing behavior of an infinite-range +/-J Ising spin glass in the presence of a transverse field using a zero-temperature quantum Monte Carlo method. Within the simulation scheme, we demonstrate that quantum annealing not only helps finding the ground state of a classical spin glass, but can also help simulating the ground state of a quantum spin glass, in particular, when the transverse field is low, much more efficiently. PMID:19256816
Correlated electron dynamics with time-dependent quantum Monte Carlo: three-dimensional helium.
Christov, Ivan P
2011-07-28
Here the recently proposed time-dependent quantum Monte Carlo method is applied to three dimensional para- and ortho-helium atoms subjected to an external electromagnetic field with amplitude sufficient to cause significant ionization. By solving concurrently sets of up to 20,000 coupled 3D time-dependent Schrödinger equations for the guide waves and corresponding sets of first order equations of motion for the Monte Carlo walkers we obtain ground state energies in close agreement with the exact values. The combined use of spherical coordinates and B-splines along the radial coordinate proves to be especially accurate and efficient for such calculations. Our results for the dipole response and the ionization of an atom with un-correlated electrons are in good agreement with the predictions of the conventional time-dependent Hartree-Fock method while the calculations with correlated electrons show enhanced ionization that is due to the electron-electron repulsion. PMID:21806103
Quantum Monte Carlo Study of the Reactions of CH with Acrolein: Major and Minor Channels.
Pakhira, Srimanta; Singh, Ravi I; Olatunji-Ojo, Olayinka; Frenklach, Michael; Lester, William A
2016-05-26
Acrolein is an important unsaturated hydrocarbon, containing both C═O and C═C bonds, and responsible for atmospheric pollution. A recent study of major reactions of CH with acrolein has been supplemented with computations of other reactions of the system. Similar to the previous approach, the quantum Monte Carlo (QMC) method in the accurate diffusion Monte Carlo (DMC) method was implemented. Single determinant wave functions were used as trial functions for the random walks. Rate coefficients and product branching ratios were computed by solving master equations using the MultiWell software suite. At room temperature, the dominant product channels are 2-methylvinyl + CO (P6), 1,3-butadienal + H (P2), and furan + H (P1). At elevated temperatures, 2,3-butadienal + H (P10) is also a major product. The chain decomposition pathway to form C3H4 + HCO was not competitive with the cyclization pathway at any of the temperatures studied. The DMC branching fractions of the products formed in the subject reaction are in reasonable accord with previous experimental and theoretical values. The computed rate coefficients were found to be independent of pressure at temperatures relevant to combustion (1500-2500 K). PMID:27046018
Dynamic load balancing for petascale quantum Monte Carlo applications: The Alias method
NASA Astrophysics Data System (ADS)
Sudheer, C. D.; Krishnan, S.; Srinivasan, A.; Kent, P. R. C.
2013-02-01
Diffusion Monte Carlo is a highly accurate Quantum Monte Carlo method for electronic structure calculations of materials, but it requires frequent load balancing or population redistribution steps to maintain efficiency on parallel machines. This step can be a significant factor affecting performance, and will become more important as the number of processing elements increases. We propose a new dynamic load balancing algorithm, the Alias Method, and evaluate it theoretically and empirically. An important feature of the new algorithm is that the load can be perfectly balanced with each process receiving at most one message. It is also optimal in the maximum size of messages received by any process. We also optimize its implementation to reduce network contention, a process facilitated by the low messaging requirement of the algorithm: a simple renumbering of the MPI ranks based on proximity and a space filling curve significantly improves the MPI Allgather performance. Empirical results on the petaflop Cray XT Jaguar supercomputer at ORNL show up to 30% improvement in performance on 120,000 cores. The load balancing algorithm may be straightforwardly implemented in existing codes. The algorithm may also be employed by any method with many near identical computational tasks that require load balancing.
Dynamic load balancing for petascale quantum Monte Carlo applications: The Alias method
Sudheer, C. D.; Krishnan, S.; Srinivasan, A.; Kent, P. R. C.
2013-02-01
Diffusion Monte Carlo is the most accurate widely used Quantum Monte Carlo method for the electronic structure of materials, but it requires frequent load balancing or population redistribution steps to maintain efficiency and avoid accumulation of systematic errors on parallel machines. The load balancing step can be a significant factor affecting performance, and will become more important as the number of processing elements increases. We propose a new dynamic load balancing algorithm, the Alias Method, and evaluate it theoretically and empirically. An important feature of the new algorithm is that the load can be perfectly balanced with each process receiving at most one message. It is also optimal in the maximum size of messages received by any process. We also optimize its implementation to reduce network contention, a process facilitated by the low messaging requirement of the algorithm. Empirical results on the petaflop Cray XT Jaguar supercomputer at ORNL showing up to 30% improvement in performance on 120,000 cores. The load balancing algorithm may be straightforwardly implemented in existing codes. The algorithm may also be employed by any method with many near identical computational tasks that requires load balancing.
Quantum Monte Carlo Studies of Interaction-Induced Localization in Quantum Dots and Wires
NASA Astrophysics Data System (ADS)
Devrim Güçlü, A.
2009-03-01
We investigate interaction-induced localization of electrons in both quantum dots and inhomogeneous quantum wires using variational and diffusion quantum Monte Carlo methods. Quantum dots and wires are highly tunable systems that enable the study of the physics of strongly correlated electrons. With decreasing electronic density, interactions become stronger and electrons are expected to localize at their classical positions, as in Wigner crystallization in an infinite 2D system. (1) Dots: We show that the addition energy shows a clear progression from features associated with shell structure to those caused by commensurability of a Wigner crystal. This cross-over is, then, a signature of localization; it occurs near rs˜20. For higher values of rs, the configuration symmetry of the quantum dot becomes fully consistent with the classical ground state. (2) Wires: We study an inhomogeneous quasi-one-dimensional system -- a wire with two regions, one at low density and the other high. We find that strong localization occurs in the low density quantum point contact region as the gate potential is increased. The nature of the transition from high to low density depends on the density gradient -- if it is steep, a barrier develops between the two regions, causing Coulomb blockade effects. We find no evidence for ferromagnetic spin polarization for the range of parameters studied. The picture emerging here is in good agreement with the experimental measurements of tunneling between two wires. Collaborators: C. J. Umrigar (Cornell), Hong Jiang (Fritz Haber Institut), Amit Ghosal (IISER Calcutta), and H. U. Baranger (Duke).
Quantum Monte Carlo methods and lithium cluster properties
Owen, R.K.
1990-12-01
Properties of small lithium clusters with sizes ranging from n = 1 to 5 atoms were investigated using quantum Monte Carlo (QMC) methods. Cluster geometries were found from complete active space self consistent field (CASSCF) calculations. A detailed development of the QMC method leading to the variational QMC (V-QMC) and diffusion QMC (D-QMC) methods is shown. The many-body aspect of electron correlation is introduced into the QMC importance sampling electron-electron correlation functions by using density dependent parameters, and are shown to increase the amount of correlation energy obtained in V-QMC calculations. A detailed analysis of D-QMC time-step bias is made and is found to be at least linear with respect to the time-step. The D-QMC calculations determined the lithium cluster ionization potentials to be 0.1982(14) [0.1981], 0.1895(9) [0.1874(4)], 0.1530(34) [0.1599(73)], 0.1664(37) [0.1724(110)], 0.1613(43) [0.1675(110)] Hartrees for lithium clusters n = 1 through 5, respectively; in good agreement with experimental results shown in the brackets. Also, the binding energies per atom was computed to be 0.0177(8) [0.0203(12)], 0.0188(10) [0.0220(21)], 0.0247(8) [0.0310(12)], 0.0253(8) [0.0351(8)] Hartrees for lithium clusters n = 2 through 5, respectively. The lithium cluster one-electron density is shown to have charge concentrations corresponding to nonnuclear attractors. The overall shape of the electronic charge density also bears a remarkable similarity with the anisotropic harmonic oscillator model shape for the given number of valence electrons.
Quantum Monte Carlo methods and lithium cluster properties. [Atomic clusters
Owen, R.K.
1990-12-01
Properties of small lithium clusters with sizes ranging from n = 1 to 5 atoms were investigated using quantum Monte Carlo (QMC) methods. Cluster geometries were found from complete active space self consistent field (CASSCF) calculations. A detailed development of the QMC method leading to the variational QMC (V-QMC) and diffusion QMC (D-QMC) methods is shown. The many-body aspect of electron correlation is introduced into the QMC importance sampling electron-electron correlation functions by using density dependent parameters, and are shown to increase the amount of correlation energy obtained in V-QMC calculations. A detailed analysis of D-QMC time-step bias is made and is found to be at least linear with respect to the time-step. The D-QMC calculations determined the lithium cluster ionization potentials to be 0.1982(14) (0.1981), 0.1895(9) (0.1874(4)), 0.1530(34) (0.1599(73)), 0.1664(37) (0.1724(110)), 0.1613(43) (0.1675(110)) Hartrees for lithium clusters n = 1 through 5, respectively; in good agreement with experimental results shown in the brackets. Also, the binding energies per atom was computed to be 0.0177(8) (0.0203(12)), 0.0188(10) (0.0220(21)), 0.0247(8) (0.0310(12)), 0.0253(8) (0.0351(8)) Hartrees for lithium clusters n = 2 through 5, respectively. The lithium cluster one-electron density is shown to have charge concentrations corresponding to nonnuclear attractors. The overall shape of the electronic charge density also bears a remarkable similarity with the anisotropic harmonic oscillator model shape for the given number of valence electrons.
i QIST: An open source continuous-time quantum Monte Carlo impurity solver toolkit
NASA Astrophysics Data System (ADS)
Huang, Li; Wang, Yilin; Meng, Zi Yang; Du, Liang; Werner, Philipp; Dai, Xi
2015-10-01
Quantum impurity solvers have a broad range of applications in theoretical studies of strongly correlated electron systems. Especially, they play a key role in dynamical mean-field theory calculations of correlated lattice models and realistic materials. Therefore, the development and implementation of efficient quantum impurity solvers is an important task. In this paper, we present an open source interacting quantum impurity solver toolkit (dubbed i QIST). This package contains several highly optimized quantum impurity solvers which are based on the hybridization expansion continuous-time quantum Monte Carlo algorithm, as well as some essential pre- and post-processing tools. We first introduce the basic principle of continuous-time quantum Monte Carlo algorithm and then discuss the implementation details and optimization strategies. The software framework, major features, and installation procedure for i QIST are also explained. Finally, several simple tutorials are presented in order to demonstrate the usage and power of i QIST.
Assaraf, Roland; Caffarel, Michel; Kollias, A C
2011-04-15
We present a method to efficiently evaluate small energy differences of two close N-body systems by employing stochastic processes having a stability versus chaos property. By using the same random noise, energy differences are computed from close trajectories without reweighting procedures. The approach is presented for quantum systems but can be applied to classical N-body systems as well. It is exemplified with diffusion Monte Carlo simulations for long chains of hydrogen atoms and molecules for which it is shown that the long-standing problem of computing energy derivatives is solved. PMID:21568537
Spin-orbit coupling in the strongly interacting Fermi gas: an exact quantum Monte Carlo study
NASA Astrophysics Data System (ADS)
Rosenberg, Peter; Shi, Hao; Chiesa, Simone; Zhang, Shiwei
Spin-orbit coupling (SOC) plays an essential role in a variety of intriguing condensed matter phenomena, including the quantum Hall effect, and topological insulators and superconductors. The recent experimental realization of spin-orbit coupled Fermi gases provides a unique opportunity to study the effects of SOC in a tunable, disorder-free system. Motivated by this experimental progress, we present here the first exact numerical results on the two-dimensional, unpolarized, uniform Fermi gas with attractive interactions and Rashba SOC. Using auxiliary-field quantum Monte Carlo and incorporating recent algorithmic advances, we carry out exact calculations on sufficiently large system sizes to provide accurate results systematically as a function of experimental parameters. We obtain the equation of state, study the spin behavior and momentum distribution, and examine the interplay of SOC and pairing in real and momentum space. Our results help illuminate the rich pairing structure induced by SOC, and provide important guidance to future experimental efforts. Supported by DOE SciDAC and NSF.
Han, Min Cheol; Yeom, Yeon Soo; Kim, Chan Hyeong; Kim, Seonghoon; Sohn, Jason W
2015-02-21
In the present study, to achieve accurate 4D Monte Carlo dose calculation in radiation therapy, we devised a new approach that combines (1) modeling of the patient body using tetrahedral-mesh geometry based on the patient's 4D CT data, (2) continuous movement/deformation of the tetrahedral patient model by interpolation of deformation vector fields acquired through deformable image registration, and (3) direct transportation of radiation particles during the movement and deformation of the tetrahedral patient model. The results of our feasibility study show that it is certainly possible to construct 4D patient models (= phantoms) with sufficient accuracy using the tetrahedral-mesh geometry and to directly transport radiation particles during continuous movement and deformation of the tetrahedral patient model. This new approach not only produces more accurate dose distribution in the patient but also replaces the current practice of using multiple 3D voxel phantoms and combining multiple dose distributions after Monte Carlo simulations. For routine clinical application of our new approach, the use of fast automatic segmentation algorithms is a must. In order to achieve, simultaneously, both dose accuracy and computation speed, the number of tetrahedrons for the lungs should be optimized. Although the current computation speed of our new 4D Monte Carlo simulation approach is slow (i.e. ~40 times slower than that of the conventional dose accumulation approach), this problem is resolvable by developing, in Geant4, a dedicated navigation class optimized for particle transportation in tetrahedral-mesh geometry. PMID:25615567
Ganesh, P.; Kim, Jeongnim; Park, Changwon; Yoon, Mina; Reboredo, Fernando A.; Kent, Paul R. C.
2014-11-03
In highly accurate diffusion quantum Monte Carlo (QMC) studies of the adsorption and diffusion of atomic lithium in AA-stacked graphite are compared with van der Waals-including density functional theory (DFT) calculations. Predicted QMC lattice constants for pure AA graphite agree with experiment. Pure AA-stacked graphite is shown to challenge many van der Waals methods even when they are accurate for conventional AB graphite. Moreover, the highest overall DFT accuracy, considering pure AA-stacked graphite as well as lithium binding and diffusion, is obtained by the self-consistent van der Waals functional vdW-DF2, although errors in binding energies remain. Empirical approaches based on point charges such as DFT-D are inaccurate unless the local charge transfer is assessed. Our results demonstrate that the lithium carbon system requires a simultaneous highly accurate description of both charge transfer and van der Waals interactions, favoring self-consistent approaches.
Ganesh, P; Kim, Jeongnim; Park, Changwon; Yoon, Mina; Reboredo, Fernando A; Kent, Paul R C
2014-12-01
Highly accurate diffusion quantum Monte Carlo (QMC) studies of the adsorption and diffusion of atomic lithium in AA-stacked graphite are compared with van der Waals-including density functional theory (DFT) calculations. Predicted QMC lattice constants for pure AA graphite agree with experiment. Pure AA-stacked graphite is shown to challenge many van der Waals methods even when they are accurate for conventional AB graphite. Highest overall DFT accuracy, considering pure AA-stacked graphite as well as lithium binding and diffusion, is obtained by the self-consistent van der Waals functional vdW-DF2, although errors in binding energies remain. Empirical approaches based on point charges such as DFT-D are inaccurate unless the local charge transfer is assessed. The results demonstrate that the lithium-carbon system requires a simultaneous highly accurate description of both charge transfer and van der Waals interactions, favoring self-consistent approaches. PMID:26583215
NASA Astrophysics Data System (ADS)
Warren, Gary Lee, Jr.
2005-11-01
Quantum Monte Carlo (QMC) methods are a class of powerful computer simulation techniques for solving the many-body Schrodinger equation. These techniques deliver essentially exact results and boast favorable computational scaling with system size. Calculations provide a full quantum mechanical treatment and may be carried to arbitrary precision. These characteristics make QMC a promising choice for the investigation of doped helium clusters, where quantum effects are substantial. Stochastic in nature, QMC methods are susceptible to statistical bias and error, which must be carefully controlled. Moreover, the relationship between the finite sampling error and the statistical uncertainty in observables has never been systematically investigated. Estimates of arbitrary observables are often substandard and can be plagued by statistical uncertainties an order of magnitude or greater than those for corresponding estimates of the energy. In this work, we present an analysis of how finite populations, importance sampling, and dimensionality affect the statistical uncertainties in QMC estimates of arbitrary observables. We find that the uncertainty depends exponentially on the dimensionality of the system, independent of the observable or nature of the system. This provides insight into the minimal population sizes and importance sampling requirements necessary to obtain useful QMC estimates of properties in high-dimensional systems. With this understanding, we develop new, more robust energy optimization procedures for cluster wavefunctions. We also implement a high quality eight parameter ansatz for the investigation of both pure and doped helium cluster systems. Compared to exact DMC results, the optimized wavefunctions recover over 90% of the total energy for clusters of size n ≤ 20. Finally, we apply this knowledge directly to the study of the solvation behavior of neutral calcium and magnesium impurities in helium nanodroplets. Diffusion Monte Carlo calculations
Quantum spin models with long-range interactions and tunnelings: a quantum Monte Carlo study
NASA Astrophysics Data System (ADS)
Maik, Michał; Hauke, Philipp; Dutta, Omjyoti; Zakrzewski, Jakub; Lewenstein, Maciej
2012-11-01
We use a quantum Monte Carlo method to investigate various classes of two-dimensional spin models with long-range interactions at low temperatures. In particular, we study a dipolar XXZ model with U(1) symmetry that appears as a hard-core boson limit of an extended Hubbard model describing polarized dipolar atoms or molecules in an optical lattice. Tunneling, in such a model, is short-range, whereas density-density couplings decay with distance following a cubic power law. We also investigate an XXZ model with long-range couplings of all three spin components—such a model describes a system of ultracold ions in a lattice of microtraps. We describe an approximate phase diagram for such systems at zero and at finite temperature, and compare their properties. In particular, we compare the extent of crystalline, superfluid and supersolid phases. Our predictions apply directly to current experiments with mesoscopic numbers of polar molecules and trapped ions.
Variational quantum Monte Carlo ground state of lithium on a Slater orbital basis
NASA Astrophysics Data System (ADS)
Eckstein, H.; Schattke, W.
1995-02-01
The ground state of bulk lithium at zero temperature is simulated by the variational quantum Monte Carlo algorithm. The total energy and its constituents are determined for two parametrized sets of trial wave functions. Including correlation by a Jastrow factor the one-determinant ansatz consists of either plane waves or a linear combination of Slater orbitals for the Li 2 s states. The latter yields results near those of the diffusion Monte Carlo algorithm.
NASA Astrophysics Data System (ADS)
Dornheim, T.; Groth, S.; Schoof, T.; Hann, C.; Bonitz, M.
2016-05-01
In a recent publication [S. Groth et al., Phys. Rev. B 93, 085102 (2016), 10.1103/PhysRevB.93.085102], we have shown that the combination of two complementary quantum Monte Carlo approaches, namely configuration path integral Monte Carlo [T. Schoof et al., Phys. Rev. Lett. 115, 130402 (2015), 10.1103/PhysRevLett.115.130402] and permutation blocking path integral Monte Carlo [T. Dornheim et al., New J. Phys. 17, 073017 (2015), 10.1088/1367-2630/17/7/073017], allows for the accurate computation of thermodynamic properties of the spin-polarized uniform electron gas over a wide range of temperatures and densities without the fixed-node approximation. In the present work, we extend this concept to the unpolarized case, which requires nontrivial enhancements that we describe in detail. We compare our simulation results with recent restricted path integral Monte Carlo data [E. W. Brown et al., Phys. Rev. Lett. 110, 146405 (2013), 10.1103/PhysRevLett.110.146405] for different energy contributions and pair distribution functions and find, for the exchange correlation energy, overall better agreement than for the spin-polarized case, while the separate kinetic and potential contributions substantially deviate.
Bohm trajectories for the Monte Carlo simulation of quantum-based devices
NASA Astrophysics Data System (ADS)
Oriols, X.; García-García, J. J.; Martín, F.; Suñé, J.; González, T.; Mateos, J.; Pardo, D.
1998-02-01
A generalization of the classical ensemble Monte Carlo (MC) device simulation technique is proposed to simultaneously deal with quantum-mechanical phase-coherence effects and scattering interactions in quantum-based devices. The proposed method restricts the quantum treatment of transport to the regions of the device where the potential profile significantly changes in distances of the order of the de Broglie wavelength of the carriers (the quantum window). Bohm trajectories associated to time-dependent Gaussian wave packets are used to simulate the electron transport in the quantum window. Outside this window, the classical ensemble MC simulation technique is used. Classical and quantum trajectories are smoothly matched at the boundaries of the quantum window according to a criterium of total-energy conservation. A self-consistent one-dimensional simulator for resonant tunneling diodes has been developed to demonstrate the feasibility of our proposal.
Ng, Yee-Hong; Bettens, Ryan P A
2016-03-01
Using the method of modified Shepard's interpolation to construct potential energy surfaces of the H2O, O3, and HCOOH molecules, we compute vibrationally averaged isotropic nuclear shielding constants ⟨σ⟩ of the three molecules via quantum diffusion Monte Carlo (QDMC). The QDMC results are compared to that of second-order perturbation theory (PT), to see if second-order PT is adequate for obtaining accurate values of nuclear shielding constants of molecules with large amplitude motions. ⟨σ⟩ computed by the two approaches differ for the hydrogens and carbonyl oxygen of HCOOH, suggesting that for certain molecules such as HCOOH where big displacements away from equilibrium happen (internal OH rotation), ⟨σ⟩ of experimental quality may only be obtainable with the use of more sophisticated and accurate methods, such as quantum diffusion Monte Carlo. The approach of modified Shepard's interpolation is also extended to construct shielding constants σ surfaces of the three molecules. By using a σ surface with the equilibrium geometry as a single data point to compute isotropic nuclear shielding constants for each descendant in the QDMC ensemble representing the ground state wave function, we reproduce the results obtained through ab initio computed σ to within statistical noise. Development of such an approach could thereby alleviate the need for any future costly ab initio σ calculations. PMID:26835785
A Variable Coefficient Method for Accurate Monte Carlo Simulation of Dynamic Asset Price
NASA Astrophysics Data System (ADS)
Li, Yiming; Hung, Chih-Young; Yu, Shao-Ming; Chiang, Su-Yun; Chiang, Yi-Hui; Cheng, Hui-Wen
2007-07-01
In this work, we propose an adaptive Monte Carlo (MC) simulation technique to compute the sample paths for the dynamical asset price. In contrast to conventional MC simulation with constant drift and volatility (μ,σ), our MC simulation is performed with variable coefficient methods for (μ,σ) in the solution scheme, where the explored dynamic asset pricing model starts from the formulation of geometric Brownian motion. With the method of simultaneously updated (μ,σ), more than 5,000 runs of MC simulation are performed to fulfills basic accuracy of the large-scale computation and suppresses statistical variance. Daily changes of stock market index in Taiwan and Japan are investigated and analyzed.
Danshita, Ippei; Polkovnikov, Anatoli
2010-09-01
We study the quantum dynamics of supercurrents of one-dimensional Bose gases in a ring optical lattice to verify instanton methods applied to coherent macroscopic quantum tunneling (MQT). We directly simulate the real-time quantum dynamics of supercurrents, where a coherent oscillation between two macroscopically distinct current states occurs due to MQT. The tunneling rate extracted from the coherent oscillation is compared with that given by the instanton method. We find that the instanton method is quantitatively accurate when the effective Planck's constant is sufficiently small. We also find phase slips associated with the oscillations.
Quantum Monte Carlo calculation of the properties of atomic carbon and diamond
Fahy, S.; Wang, X.W.; Louie, S.G.
1988-06-01
A new method of calculating total energies of solids using non-local pseudopotentials in conjunction with the variational quantum Monte Carlo approach is presented. By using pseudopotentials, the large fluctuations of the energies in the core region of the atoms which occur in quantum Monte Carlo all-electron schemes are avoided. The method is applied to calculate the cohesive energy and structural properties of diamond and the first ionization energy and electron affinity of the carbon atom. Results are in excellent agreement with experiment. 8 refs., 1 fig., 2 tabs.
Quantum Monte Carlo Benchmark of Exchange-Correlation Functionals for Bulk Water.
Morales, Miguel A; Gergely, John R; McMinis, Jeremy; McMahon, Jeffrey M; Kim, Jeongnim; Ceperley, David M
2014-06-10
The accurate description of the thermodynamic and dynamical properties of liquid water from first-principles is a very important challenge to the theoretical community. This represents not only a critical test of the predictive capabilities of first-principles methods, but it will also shed light into the microscopic properties of such an important substance. Density Functional Theory, the main workhorse in the field of first-principles methods, has been so far unable to properly describe water and its unusual properties in the liquid state. With the recent introduction of exact exchange and an improved description of dispersion interaction, the possibility of an accurate description of the liquid is finally within reach. Unfortunately, there is still no way to systematically improve exchange-correlation functionals, and the number of available functionals is very large. In this article we use highly accurate quantum Monte Carlo calculations to benchmark a selection of exchange-correlation functionals typically used in Density Functional Theory simulations of bulk water. This allows us to test the predictive capabilities of these functionals in water, giving us a way to choose optimal functionals for first-principles simulations. We compare and contrast the importance of different features of functionals, including the hybrid component, the vdW component, and their importance within different aspects of the PES. In addition, in order to correct the inaccuracies in the description of short-range interactions in the liquid, we test a recently introduced scheme that combines Density Functional Theory with Coupled Cluster calculations through a Many-Body expansion of the energy. PMID:26580755
An introduction to applied quantum mechanics in the Wigner Monte Carlo formalism
NASA Astrophysics Data System (ADS)
Sellier, J. M.; Nedjalkov, M.; Dimov, I.
2015-05-01
The Wigner formulation of quantum mechanics is a very intuitive approach which allows the comprehension and prediction of quantum mechanical phenomena in terms of quasi-distribution functions. In this review, our aim is to provide a detailed introduction to this theory along with a Monte Carlo method for the simulation of time-dependent quantum systems evolving in a phase-space. This work consists of three main parts. First, we introduce the Wigner formalism, then we discuss in detail the Wigner Monte Carlo method and, finally, we present practical applications. In particular, the Wigner model is first derived from the Schrödinger equation. Then a generalization of the formalism due to Moyal is provided, which allows to recover important mathematical properties of the model. Next, the Wigner equation is further generalized to the case of many-body quantum systems. Finally, a physical interpretation of the negative part of a quasi-distribution function is suggested. In the second part, the Wigner Monte Carlo method, based on the concept of signed (virtual) particles, is introduced in detail for the single-body problem. Two extensions of the Wigner Monte Carlo method to quantum many-body problems are introduced, in the frameworks of time-dependent density functional theory and ab-initio methods. Finally, in the third and last part of this paper, applications to single- and many-body problems are performed in the context of quantum physics and quantum chemistry, specifically focusing on the hydrogen, lithium and boron atoms, the H2 molecule and a system of two identical Fermions. We conclude this work with a discussion on the still unexplored directions the Wigner Monte Carlo method could take in the next future.
Grover search algorithm with Rydberg-blockaded atoms: quantum Monte Carlo simulations
NASA Astrophysics Data System (ADS)
Petrosyan, David; Saffman, Mark; Mølmer, Klaus
2016-05-01
We consider the Grover search algorithm implementation for a quantum register of size N={2}k using k (or k+1) microwave- and laser-driven Rydberg-blockaded atoms, following the proposal by Mølmer et al (2011 J. Phys. B 44 184016). We suggest some simplifications for the microwave and laser couplings, and analyze the performance of the algorithm for up to k = 4 multilevel atoms under realistic experimental conditions using quantum stochastic (Monte Carlo) wavefunction simulations.
Kersten, J A F; Booth, George H; Alavi, Ali
2016-08-01
The Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method has proved able to provide near-exact solutions to the electronic Schrödinger equation within a finite orbital basis set, without relying on an expansion about a reference state. However, a drawback to the approach is that being based on an expansion of Slater determinants, the FCIQMC method suffers from a basis set incompleteness error that decays very slowly with the size of the employed single particle basis. The FCIQMC results obtained in a small basis set can be improved significantly with explicitly correlated techniques. Here, we present a study that assesses and compares two contrasting "universal" explicitly correlated approaches that fit into the FCIQMC framework: the [2]R12 method of Kong and Valeev [J. Chem. Phys. 135, 214105 (2011)] and the explicitly correlated canonical transcorrelation approach of Yanai and Shiozaki [J. Chem. Phys. 136, 084107 (2012)]. The former is an a posteriori internally contracted perturbative approach, while the latter transforms the Hamiltonian prior to the FCIQMC simulation. These comparisons are made across the 55 molecules of the G1 standard set. We found that both methods consistently reduce the basis set incompleteness, for accurate atomization energies in small basis sets, reducing the error from 28 mEh to 3-4 mEh. While many of the conclusions hold in general for any combination of multireference approaches with these methodologies, we also consider FCIQMC-specific advantages of each approach. PMID:27497549
Jorgensen, Wiliiam L.
2014-01-01
A recent review (Acc. Chem. Res. 2010, 43:142–151) examined our use and development of a combined quantum and molecular mechanical (QM/MM) technique for modelling organic and enzymatic reactions. Advances included the PDDG/PM3 semiempirical QM (SQM) method, computation of multi-dimensional potentials of mean force (PMF), incorporation of on-the-fly QM in Monte Carlo simulations, and a polynomial quadrature method for rapidly treating proton-transfer reactions. The current article serves as a follow up on our progress. Highlights include new reactions, alternative SQM methods, a polarizable OPLS force field, and novel solvent environments, e.g., “on water” and room temperature ionic liquids. The methodology is strikingly accurate across a wide range of condensed-phase and antibody-catalyzed reactions including substitution, decarboxylation, elimination, isomerization, and pericyclic classes. Comparisons are made to systems treated with continuum-based solvents and ab initio or density functional theory (DFT) methods. Overall, the QM/MM methodology provides detailed characterization of reaction paths, proper configurational sampling, several advantages over implicit solvent models, and a reasonable computational cost. PMID:25431625
A deterministic alternative to the full configuration interaction quantum Monte Carlo method
NASA Astrophysics Data System (ADS)
Tubman, Norm M.; Lee, Joonho; Takeshita, Tyler Y.; Head-Gordon, Martin; Whaley, K. Birgitta
2016-07-01
Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such technique, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, along with a stochastic projected wave function, to find the important parts of Hilbert space. However, the stochastic representation of the wave function is not required to search Hilbert space efficiently, and here we describe a highly efficient deterministic method that can achieve chemical accuracy for a wide range of systems, including the difficult Cr2 molecule. We demonstrate for systems like Cr2 that such calculations can be performed in just a few cpu hours which makes it one of the most efficient and accurate methods that can attain chemical accuracy for strongly correlated systems. In addition our method also allows efficient calculation of excited state energies, which we illustrate with benchmark results for the excited states of C2.
Diffusion Quantum Monte Carlo predictions for bulk MnNiO3
NASA Astrophysics Data System (ADS)
Mitra, Chandrima; Krogel, Jaron; Reboredo, Fernando A.
MnNiO3 is a strongly correlated transition metal oxide that has recently been investigated theoretically for its potential application as an oxygen-evolution photo-catalyst. However, there is no experimental report on critical quantities like its band gap or its bulk modulus. Recent theoretical predictions with standard functionals, such as PBE +U and HSE show large discrepancies in the band-gaps (about 1.23 eV), depending on the nature of the functional used. Hence, there is clearly a need for an accurate quantitative prediction of the band-gap in order to decide its usefulness as a photo-catalyst. In this work, we present Diffusion Quantum Monte Carlo (DMC) study of the bulk properties of MnNiO3. This includes the quasiparticle band gap for the two spin channels, the equilibrium lattice parameter and the bulk modulus. The DMC approach has already been shown to achieve excellent agreement with experimental results for other oxides such as ZnO NiO and Fe2O3. To our knowledge, MnNiO3 is the first case where this theory is applied before experiments are done. This work was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division.
NASA Astrophysics Data System (ADS)
Kersten, J. A. F.; Booth, George H.; Alavi, Ali
2016-08-01
The Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method has proved able to provide near-exact solutions to the electronic Schrödinger equation within a finite orbital basis set, without relying on an expansion about a reference state. However, a drawback to the approach is that being based on an expansion of Slater determinants, the FCIQMC method suffers from a basis set incompleteness error that decays very slowly with the size of the employed single particle basis. The FCIQMC results obtained in a small basis set can be improved significantly with explicitly correlated techniques. Here, we present a study that assesses and compares two contrasting "universal" explicitly correlated approaches that fit into the FCIQMC framework: the [2]R12 method of Kong and Valeev [J. Chem. Phys. 135, 214105 (2011)] and the explicitly correlated canonical transcorrelation approach of Yanai and Shiozaki [J. Chem. Phys. 136, 084107 (2012)]. The former is an a posteriori internally contracted perturbative approach, while the latter transforms the Hamiltonian prior to the FCIQMC simulation. These comparisons are made across the 55 molecules of the G1 standard set. We found that both methods consistently reduce the basis set incompleteness, for accurate atomization energies in small basis sets, reducing the error from 28 mEh to 3-4 mEh. While many of the conclusions hold in general for any combination of multireference approaches with these methodologies, we also consider FCIQMC-specific advantages of each approach.
Bond Breaking of Simple Molecules in Auxiliary-Field Quantum Monte Carlo with GVB Wave Functions
NASA Astrophysics Data System (ADS)
Purwanto, Wirawan; Krakauer, Henry; Zhang, Shiwei; Al-Saidi, Wissam
2007-03-01
Accurate potential energy curves are an essential ingredient in understanding chemical reactions. This problem spans a wide range of correlations, with correlation effects being the most important in the bond-breaking regime. We report potential energy curves of simple molecules, including water and the carbon dimer, within the framework of the auxiliary-field quantum Monte Carlo (AFQMC) method. AFQMC projects the many-body ground-state from a trial wave function, which is also used to control the sign/phase problem. A previous calculation showed that AFQMC could provide a fairly uniform description of the bond stretching of a water molecule, even with a simple unrestricted Hartree-Fock (UHF) trial wave function. We investigate the use of Generalized Valence Bond (GVB). GVB gives a better description of the molecule than UHF; so it is a simple yet efficient alternative to using a single Slater determinant trial wave function. We will compare AFQMC results with other correlated methods and the exact configuration interaction calculations. Al-Saidi, Zhang, Krakauer, J. Chem. Phys. 124, 224101 (2006)
Quantum Monte Carlo Simulations of Adulteration Effect on Bond Alternating Spin=1/2 Chain
NASA Astrophysics Data System (ADS)
Zhang, Peng; Xu, Zhaoxin; Ying, Heping; Dai, Jianhui; Crompton, Peter
The S=1/2 Heisenberg chain with bond alternation and randomness of antiferromagnetic (AFM) and ferromagnetic (FM) interactions is investigated by quantum Monte Carlo simulations of loop/cluster algorithm. Our results have shown interesting finite temperature magnetic properties of this model. The relevance of our study to former investigation results is discussed.
Quantum-trajectory Monte Carlo method for study of electron-crystal interaction in STEM.
Ruan, Z; Zeng, R G; Ming, Y; Zhang, M; Da, B; Mao, S F; Ding, Z J
2015-07-21
In this paper, a novel quantum-trajectory Monte Carlo simulation method is developed to study electron beam interaction with a crystalline solid for application to electron microscopy and spectroscopy. The method combines the Bohmian quantum trajectory method, which treats electron elastic scattering and diffraction in a crystal, with a Monte Carlo sampling of electron inelastic scattering events along quantum trajectory paths. We study in this work the electron scattering and secondary electron generation process in crystals for a focused incident electron beam, leading to understanding of the imaging mechanism behind the atomic resolution secondary electron image that has been recently achieved in experiment with a scanning transmission electron microscope. According to this method, the Bohmian quantum trajectories have been calculated at first through a wave function obtained via a numerical solution of the time-dependent Schrödinger equation with a multislice method. The impact parameter-dependent inner-shell excitation cross section then enables the Monte Carlo sampling of ionization events produced by incident electron trajectories travelling along atom columns for excitation of high energy knock-on secondary electrons. Following cascade production, transportation and emission processes of true secondary electrons of very low energies are traced by a conventional Monte Carlo simulation method to present image signals. Comparison of the simulated image for a Si(110) crystal with the experimental image indicates that the dominant mechanism of atomic resolution of secondary electron image is the inner-shell ionization events generated by a high-energy electron beam. PMID:26082190
Bayesian inference and the analytic continuation of imaginary-time quantum Monte Carlo data
Gubernatis, J.E.; Bonca, J.; Jarrell, M.
1995-12-31
We present brief description of how methods of Bayesian inference are used to obtain real frequency information by the analytic continuation of imaginary-time quantum Monte Carlo data. We present the procedure we used, which is due to R. K. Bryan, and summarize several bottleneck issues.
Monte-Carlo Quantum Chemistry of Biogene Amines. Laser and Neutron Capture Effects
Glushkov, A. V.; Malinovskaya, S. V.; Khetselius, O. Yu.; Loboda, A. V.
2009-03-09
Monte-Carlo quantum calculation of the cluster consisting of the serotonine ST (histamine HM) molecules and 100 molecules of water is carried out. It is found that the zwitterion appears as expected to be strongly favoured with respect to neutral molecule. The perspective possibilities of laser and neutron capture action on different biomolecules are indicated.
Monte-Carlo Quantum Chemistry of Biogene Amines. Laser and Neutron Capture Effects
NASA Astrophysics Data System (ADS)
Glushkov, A. V.; Malinovskaya, S. V.; Khetselius, O. Yu.; Loboda, A. V.
2009-03-01
Monte-Carlo quantum calculation of the cluster consisting of the serotonine ST (histamine HM) molecules and 100 molecules of water is carried out. It is found that the zwitterion appears as expected to be strongly favoured with respect to neutral molecule. The perspective possibilities of laser and neutron capture action on different biomolecules are indicated.
Multi-Jastrow trial wavefunctions for electronic structure calculations with quantum Monte Carlo.
Bouabça, Thomas; Braïda, Benoît; Caffarel, Michel
2010-07-28
A new type of electronic trial wavefunction suitable for quantum Monte Carlo calculations of molecular systems is presented. In contrast with the standard Jastrow-Slater form built with a unique global Jastrow term, it is proposed to introduce individual Jastrow factors attached to molecular orbitals. Such a form is expected to be more physical since it allows to describe differently the local electronic correlations associated with various molecular environments (1s-core orbitals, 3d-magnetic orbitals, localized two-center sigma-orbitals, delocalized pi-orbitals, atomic lone pairs, etc.). In contrast with the standard form, introducing different Jastrow terms allows us to change the nodal structure of the wavefunction, a point which is important in the context of building better nodes for more accurate fixed-node diffusion Monte Carlo (FN-DMC) calculations. Another important aspect resulting from the use of local Jastrow terms is the possibility of defining and preoptimizing local and transferable correlated units for building complex trial wavefunctions from simple parts. The practical aspects associated with the computation of the intricate derivatives of the multi-Jastrow trial function are presented in detail. Some first illustrative applications for atoms of increasing size (O, S, and Cu) and for the potential energy curve and spectroscopic constants of the FH molecule are presented. In the case of the copper atom, the use of the multi-Jastrow form at the variational Monte Carlo level has allowed us to improve significantly the value of the total ground-state energy (about 75% of the correlation energy with only one determinant and three atomic orbital Jastrow factors). In the case of the FH molecule (fluorine hydride), it has been found that the multi-Jastrow nodes lead to an almost exact FN-DMC value of the dissociation energy [D(0)=-140.7(4) kcal/mol instead of the estimated nonrelativistic Born-Oppenheimer exact value of -141.1], which is not the case
A High-Accurate and High-Efficient Monte Carlo Code by Improved Molière Functions with Ionization
NASA Astrophysics Data System (ADS)
Nakatsuka, Takao; Okei, Kazuhide
2003-07-01
Although the Molière theory of multiple Coulomb scattering is less accue rate in tracing solid angles than the Goudsmit and Saunderson theory due to the small angle approximation, it still acts very important roles in developments of high-efficient simulation codes of relativistic charged particles like cosmic-ray particles. Molière expansion is well explained by the physical model, that is the e normal distribution attributing to the high-frequent moderate scatterings and subsequent correction terms attributing to the additive large-angle scatterings. Based on these physical concepts, we have improved a high-accurate and highefficient Monte Carlo code taking account of ionization loss.
Communication: Variation after response in quantum Monte Carlo.
Neuscamman, Eric
2016-08-28
We present a new method for modeling electronically excited states that overcomes a key failing of linear response theory by allowing the underlying ground state ansatz to relax in the presence of an excitation. The method is variational, has a cost similar to ground state variational Monte Carlo, and admits both open and periodic boundary conditions. We present preliminary numerical results showing that, when paired with the Jastrow antisymmetric geminal power ansatz, the variation-after-response formalism delivers accuracies for valence and charge transfer single excitations on par with equation of motion coupled cluster, while surpassing coupled cluster's accuracy for excitations with significant doubly excited character. PMID:27586897
Luo, Ye Sorella, Sandro; Zen, Andrea
2014-11-21
We present a systematic study of a recently developed ab initio simulation scheme based on molecular dynamics and quantum Monte Carlo. In this approach, a damped Langevin molecular dynamics is employed by using a statistical evaluation of the forces acting on each atom by means of quantum Monte Carlo. This allows the use of an highly correlated wave function parametrized by several variational parameters and describing quite accurately the Born-Oppenheimer energy surface, as long as these parameters are determined at the minimum energy condition. However, in a statistical method both the minimization method and the evaluation of the atomic forces are affected by the statistical noise. In this work, we study systematically the accuracy and reliability of this scheme by targeting the vibrational frequencies of simple molecules such as the water monomer, hydrogen sulfide, sulfur dioxide, ammonia, and phosphine. We show that all sources of systematic errors can be controlled and reliable frequencies can be obtained with a reasonable computational effort. This work provides convincing evidence that this molecular dynamics scheme can be safely applied also to realistic systems containing several atoms.
Rota, R; Casulleras, J; Mazzanti, F; Boronat, J
2015-03-21
We present a method based on the path integral Monte Carlo formalism for the calculation of ground-state time correlation functions in quantum systems. The key point of the method is the consideration of time as a complex variable whose phase δ acts as an adjustable parameter. By using high-order approximations for the quantum propagator, it is possible to obtain Monte Carlo data all the way from purely imaginary time to δ values near the limit of real time. As a consequence, it is possible to infer accurately the spectral functions using simple inversion algorithms. We test this approach in the calculation of the dynamic structure function S(q, ω) of two one-dimensional model systems, harmonic and quartic oscillators, for which S(q, ω) can be exactly calculated. We notice a clear improvement in the calculation of the dynamic response with respect to the common approach based on the inverse Laplace transform of the imaginary-time correlation function. PMID:25796238
Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method
NASA Astrophysics Data System (ADS)
Qin, Mingpu; Shi, Hao; Zhang, Shiwei
2016-08-01
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U /t from 2 to 8, for a range of densities including half-filling and n =0.3 ,0.5 ,0.6 ,0.75 , and 0.875 . At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to 16 ×16 are also provided for benchmark purposes.
Quantum Monte Carlo calculations of neutron and nuclear matter
NASA Astrophysics Data System (ADS)
Gandolfi, Stefano
2014-09-01
Recent advances in experiments of the symmetry energy of nuclear matter and in neutron star observations yield important new insights on the equation of state of neutron matter at nuclear densities. In this regime the EOS of neutron matter plays a critical role in determining the mass-radius relationship for neutron stars. We show how microscopic calculations of neutron matter, based on realistic two- and three-nucleon forces, make clear predictions for the relation between the isospin-asymmetry energy of nuclear matter and its density dependence, and the maximum mass and radius for a neutron star. We will also discuss the recent extension of the Auxiliary Field Diffusion Monte Carlo method to study the equation of state of nuclear matter using two-body nucleon interactions. The equation of state of isospin-asymmetric nuclear matter will also be discussed.
NASA Astrophysics Data System (ADS)
Infantino, Angelo; Marengo, Mario; Baschetti, Serafina; Cicoria, Gianfranco; Longo Vaschetto, Vittorio; Lucconi, Giulia; Massucci, Piera; Vichi, Sara; Zagni, Federico; Mostacci, Domiziano
2015-11-01
Biomedical cyclotrons for production of Positron Emission Tomography (PET) radionuclides and radiotherapy with hadrons or ions are widely diffused and established in hospitals as well as in industrial facilities and research sites. Guidelines for site planning and installation, as well as for radiation protection assessment, are given in a number of international documents; however, these well-established guides typically offer analytic methods of calculation of both shielding and materials activation, in approximate or idealized geometry set up. The availability of Monte Carlo codes with accurate and up-to-date libraries for transport and interactions of neutrons and charged particles at energies below 250 MeV, together with the continuously increasing power of nowadays computers, makes systematic use of simulations with realistic geometries possible, yielding equipment and site specific evaluation of the source terms, shielding requirements and all quantities relevant to radiation protection. In this work, the well-known Monte Carlo code FLUKA was used to simulate two representative models of cyclotron for PET radionuclides production, including their targetry; and one type of proton therapy cyclotron including the energy selection system. Simulations yield estimates of various quantities of radiological interest, including the effective dose distribution around the equipment, the effective number of neutron produced per incident proton and the activation of target materials, the structure of the cyclotron, the energy degrader, the vault walls and the soil. The model was validated against experimental measurements and comparison with well-established reference data. Neutron ambient dose equivalent H*(10) was measured around a GE PETtrace cyclotron: an average ratio between experimental measurement and simulations of 0.99±0.07 was found. Saturation yield of 18F, produced by the well-known 18O(p,n)18F reaction, was calculated and compared with the IAEA recommended
Quantum Monte Carlo calculations of neutron-alpha scattering.
Nollett, K. M.; Pieper, S. C.; Wiringa, R. B.; Carlson, J.; Hale, G. M.; Physics
2007-07-13
We describe a new method to treat low-energy scattering problems in few-nucleon systems, and we apply it to the five-body case of neutron-alpha scattering. The method allows precise calculations of low-lying resonances and their widths. We find that a good three-nucleon interaction is crucial to obtain an accurate description of neutron-alpha scattering.
Quantum Monte Carlo Calculations of Neutron-{alpha} Scattering
Nollett, Kenneth M.; Pieper, Steven C.; Wiringa, R. B.; Carlson, J.; Hale, G. M.
2007-07-13
We describe a new method to treat low-energy scattering problems in few-nucleon systems, and we apply it to the five-body case of neutron-alpha scattering. The method allows precise calculations of low-lying resonances and their widths. We find that a good three-nucleon interaction is crucial to obtain an accurate description of neutron-alpha scattering.
NASA Astrophysics Data System (ADS)
Luitz, David J.; Alet, Fabien; Laflorencie, Nicolas
2014-03-01
Shannon-Renyi entropies are measures of the participation of basis states in a wave function. Previous work for one dimensional systems showed that they exhibit a subleading scaling behavior with system size that contains universal information, such as e.g. the Luttinger Liquid parameter. Here, we introduce quantum Monte Carlo schemes to calculate these quantities and the related participation spectra for unfrustrated quantum many body systems in any dimension and apply them to interacting spin systems. Our results demonstrate the universality of subleading scaling terms for different kinds of phase transitions with a spontaneous breaking of discrete or continuous symmetries and at quantum critical points. Aditionally, we also discuss the signature of quantum phase transitions in the participation spectra of subsystems.
Structural stability and defect energetics of ZnO from diffusion quantum Monte Carlo
Santana, Juan A.; Krogel, Jaron T.; Kim, Jeongnim; Reboredo, Fernando A.; Kent, Paul R. C.
2015-04-28
We have applied the many-body ab initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure and the energetics of the oxygen vacancy, zinc interstitial, and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory (DFT) approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O{sub 2}, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.
Towards prediction of correlated material properties using quantum Monte Carlo methods
NASA Astrophysics Data System (ADS)
Wagner, Lucas
Correlated electron systems offer a richness of physics far beyond noninteracting systems. If we would like to pursue the dream of designer correlated materials, or, even to set a more modest goal, to explain in detail the properties and effective physics of known materials, then accurate simulation methods are required. Using modern computational resources, quantum Monte Carlo (QMC) techniques offer a way to directly simulate electron correlations. I will show some recent results on a few extremely challenging materials including the metal-insulator transition of VO2, the ground state of the doped cuprates, and the pressure dependence of magnetic properties in FeSe. By using a relatively simple implementation of QMC, at least some properties of these materials can be described truly from first principles, without any adjustable parameters. Using the QMC platform, we have developed a way of systematically deriving effective lattice models from the simulation. This procedure is particularly attractive for correlated electron systems because the QMC methods treat the one-body and many-body components of the wave function and Hamiltonian on completely equal footing. I will show some examples of using this downfolding technique and the high accuracy of QMC to connect our intuitive ideas about interacting electron systems with high fidelity simulations. The work in this presentation was supported in part by NSF DMR 1206242, the U.S. Department of Energy, Office of Science, Office of Advanced Scientific Computing Research, Scientific Discovery through Advanced Computing (SciDAC) program under Award Number FG02-12ER46875, and the Center for Emergent Superconductivity, Department of Energy Frontier Research Center under Grant No. DEAC0298CH1088. Computing resources were provided by a Blue Waters Illinois grant and INCITE PhotSuper and SuperMatSim allocations.
Structural Stability and Defect Energetics of ZnO from Diffusion Quantum Monte Carlo
Santana Palacio, Juan A.; Krogel, Jaron T.; Kim, Jeongnim; Kent, Paul R.; Reboredo, Fernando A.
2015-04-28
We have applied the many-body ab-initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure, and the energetics of the oxygen vacancy, zinc interstitial and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O2, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy asmore » a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type, and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.« less
Structural Stability and Defect Energetics of ZnO from Diffusion Quantum Monte Carlo
Santana Palacio, Juan A.; Krogel, Jaron T.; Kim, Jeongnim; Kent, Paul R.; Reboredo, Fernando A.
2015-04-28
We have applied the many-body ab-initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure, and the energetics of the oxygen vacancy, zinc interstitial and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O2, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type, and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.
Quantum Monte Carlo studies of relativistic effects in light nuclei
J. L. Forest; V. R. Pandharipande; A. Arriaga
1998-05-01
Relativistic Hamiltonians are defined as the sum of relativistic one-body kinetic energy, two- and three-body potentials and their boost corrections. In this work the authors use the variational Monte Carlo method to study two kinds of relativistic effects in the binding energy of {sup 3}H and {sup 4}He. The first is due to the nonlocalities in the relativistic kinetic energy and relativistic one-pion exchange potential (OPEP), and the second is from boost interaction. The OPEP contribution is reduced by about 15% by the relativistic nonlocality, which may also have significant effects on pion exchange currents. However, almost all of this reduction is canceled by changes in the kinetic energy and other interaction terms, and the total effect of the nonlocalities on the binding energy is very small. The boost interactions, on the other hand, give repulsive contributions of 0.4 (1.9) MeV in {sup 3}H ({sup 4}He) and account for 37% of the phenomenological part of the three-nucleon interaction needed in the nonrelativistic Hamiltonians.
Hyperon Puzzle: Hints from Quantum Monte Carlo Calculations
NASA Astrophysics Data System (ADS)
Lonardoni, Diego; Lovato, Alessandro; Gandolfi, Stefano; Pederiva, Francesco
2015-03-01
The onset of hyperons in the core of neutron stars and the consequent softening of the equation of state have been questioned for a long time. Controversial theoretical predictions and recent astrophysical observations of neutron stars are the grounds for the so-called hyperon puzzle. We calculate the equation of state and the neutron star mass-radius relation of an infinite systems of neutrons and Λ particles by using the auxiliary field diffusion Monte Carlo algorithm. We find that the three-body hyperon-nucleon interaction plays a fundamental role in the softening of the equation of state and for the consequent reduction of the predicted maximum mass. We have considered two different models of three-body force that successfully describe the binding energy of medium mass hypernuclei. Our results indicate that they give dramatically different results on the maximum mass of neutron stars, not necessarily incompatible with the recent observation of very massive neutron stars. We conclude that stronger constraints on the hyperon-neutron force are necessary in order to properly assess the role of hyperons in neutron stars.
Hyperon puzzle: hints from quantum Monte Carlo calculations.
Lonardoni, Diego; Lovato, Alessandro; Gandolfi, Stefano; Pederiva, Francesco
2015-03-01
The onset of hyperons in the core of neutron stars and the consequent softening of the equation of state have been questioned for a long time. Controversial theoretical predictions and recent astrophysical observations of neutron stars are the grounds for the so-called hyperon puzzle. We calculate the equation of state and the neutron star mass-radius relation of an infinite systems of neutrons and Λ particles by using the auxiliary field diffusion Monte Carlo algorithm. We find that the three-body hyperon-nucleon interaction plays a fundamental role in the softening of the equation of state and for the consequent reduction of the predicted maximum mass. We have considered two different models of three-body force that successfully describe the binding energy of medium mass hypernuclei. Our results indicate that they give dramatically different results on the maximum mass of neutron stars, not necessarily incompatible with the recent observation of very massive neutron stars. We conclude that stronger constraints on the hyperon-neutron force are necessary in order to properly assess the role of hyperons in neutron stars. PMID:25793808
Majorana Positivity and the Fermion Sign Problem of Quantum Monte Carlo Simulations
NASA Astrophysics Data System (ADS)
Wei, Z. C.; Wu, Congjun; Li, Yi; Zhang, Shiwei; Xiang, T.
2016-06-01
The sign problem is a major obstacle in quantum Monte Carlo simulations for many-body fermion systems. We examine this problem with a new perspective based on the Majorana reflection positivity and Majorana Kramers positivity. Two sufficient conditions are proven for the absence of the fermion sign problem. Our proof provides a unified description for all the interacting lattice fermion models previously known to be free of the sign problem based on the auxiliary field quantum Monte Carlo method. It also allows us to identify a number of new sign-problem-free interacting fermion models including, but not limited to, lattice fermion models with repulsive interactions but without particle-hole symmetry, and interacting topological insulators with spin-flip terms.
Majorana Positivity and the Fermion Sign Problem of Quantum Monte Carlo Simulations.
Wei, Z C; Wu, Congjun; Li, Yi; Zhang, Shiwei; Xiang, T
2016-06-24
The sign problem is a major obstacle in quantum Monte Carlo simulations for many-body fermion systems. We examine this problem with a new perspective based on the Majorana reflection positivity and Majorana Kramers positivity. Two sufficient conditions are proven for the absence of the fermion sign problem. Our proof provides a unified description for all the interacting lattice fermion models previously known to be free of the sign problem based on the auxiliary field quantum Monte Carlo method. It also allows us to identify a number of new sign-problem-free interacting fermion models including, but not limited to, lattice fermion models with repulsive interactions but without particle-hole symmetry, and interacting topological insulators with spin-flip terms. PMID:27391709
Charged vanadium-benzene multidecker clusters: DFT and quantum Monte Carlo study
NASA Astrophysics Data System (ADS)
Tokár, K.; Derian, R.; Mitas, L.; Štich, I.
2016-02-01
Using explicitly correlated fixed-node quantum Monte Carlo and density functional theory (DFT) methods, we study electronic properties, ground-state multiplets, ionization potentials, electron affinities, and low-energy fragmentation channels of charged half-sandwich and multidecker vanadium-benzene systems with up to 3 vanadium atoms, including both anions and cations. It is shown that, particularly in anions, electronic correlations play a crucial role; these effects are not systematically captured with any commonly used DFT functionals such as gradient corrected, hybrids, and range-separated hybrids. On the other hand, tightly bound cations can be described qualitatively by DFT. A comparison of DFT and quantum Monte Carlo provides an in-depth understanding of the electronic structure and properties of these correlated systems. The calculations also serve as a benchmark study of 3d molecular anions that require a balanced many-body description of correlations at both short- and long-range distances.
Charged vanadium-benzene multidecker clusters: DFT and quantum Monte Carlo study.
Tokár, K; Derian, R; Mitas, L; Štich, I
2016-02-14
Using explicitly correlated fixed-node quantum Monte Carlo and density functional theory (DFT) methods, we study electronic properties, ground-state multiplets, ionization potentials, electron affinities, and low-energy fragmentation channels of charged half-sandwich and multidecker vanadium-benzene systems with up to 3 vanadium atoms, including both anions and cations. It is shown that, particularly in anions, electronic correlations play a crucial role; these effects are not systematically captured with any commonly used DFT functionals such as gradient corrected, hybrids, and range-separated hybrids. On the other hand, tightly bound cations can be described qualitatively by DFT. A comparison of DFT and quantum Monte Carlo provides an in-depth understanding of the electronic structure and properties of these correlated systems. The calculations also serve as a benchmark study of 3d molecular anions that require a balanced many-body description of correlations at both short- and long-range distances. PMID:26874484
Ferromagnetism of a Repulsive Atomic Fermi Gas in an Optical Lattice: A Quantum Monte Carlo Study
NASA Astrophysics Data System (ADS)
Pilati, Sebastiano; Zintchenko, Ilia; Troyer, Matthias
2015-05-01
We investigate the ferromagnetic behavior of a two-component repulsive Fermi gas under the influence of a periodic potential that describes the effect of a 3D optical lattice, using continuous-space quantum Monte Carlo simulations. We find that a shallow optical lattice below half-filling strongly favors the ferromagnetic instability compared to the homogeneous Fermi gas. Instead, in the regime of deep optical lattices and weak interactions, where the conventional description in terms of single-band tight-binding models is reliable, our results indicate that the paramagnetic state is stable, in agreement with previous quantum Monte Carlo simulations of the Hubbard model. Our findings shed light on the important role played by multi-band effects and by interaction-induced hopping in the physics of atomic gases trapped in optical lattices.
Rasch, Kevin M.; Hu, Shuming; Mitas, Lubos
2014-01-28
We elucidate the origin of large differences (two-fold or more) in the fixed-node errors between the first- vs second-row systems for single-configuration trial wave functions in quantum Monte Carlo calculations. This significant difference in the valence fixed-node biases is studied across a set of atoms, molecules, and also Si, C solid crystals. We show that the key features which affect the fixed-node errors are the differences in electron density and the degree of node nonlinearity. The findings reveal how the accuracy of the quantum Monte Carlo varies across a variety of systems, provide new perspectives on the origins of the fixed-node biases in calculations of molecular and condensed systems, and carry implications for pseudopotential constructions for heavy elements.
A sparse algorithm for the evaluation of the local energy in quantum Monte Carlo.
Aspuru-Guzik, Alán; Salomón-Ferrer, Romelia; Austin, Brian; Lester, William A
2005-05-01
A new algorithm is presented for the sparse representation and evaluation of Slater determinants in the quantum Monte Carlo (QMC) method. The approach, combined with the use of localized orbitals in a Slater-type orbital basis set, significantly extends the size molecule that can be treated with the QMC method. Application of the algorithm to systems containing up to 390 electrons confirms that the cost of evaluating the Slater determinant scales linearly with system size. PMID:15761862
NASA Astrophysics Data System (ADS)
Seth, Priyanka; Krivenko, Igor; Ferrero, Michel; Parcollet, Olivier
2016-03-01
We present TRIQS/CTHYB, a state-of-the art open-source implementation of the continuous-time hybridisation expansion quantum impurity solver of the TRIQS package. This code is mainly designed to be used with the TRIQS library in order to solve the self-consistent quantum impurity problem in a multi-orbital dynamical mean field theory approach to strongly-correlated electrons, in particular in the context of realistic electronic structure calculations. It is implemented in C++ for efficiency and is provided with a high-level Python interface. The code ships with a new partitioning algorithm that divides the local Hilbert space without any user knowledge of the symmetries and quantum numbers of the Hamiltonian. Furthermore, we implement higher-order configuration moves and show that such moves are necessary to ensure ergodicity of the Monte Carlo in common Hamiltonians even without symmetry-breaking.
Introducing QMC/MMpol: Quantum Monte Carlo in Polarizable Force Fields for Excited States.
Guareschi, Riccardo; Zulfikri, Habiburrahman; Daday, Csaba; Floris, Franca Maria; Amovilli, Claudio; Mennucci, Benedetta; Filippi, Claudia
2016-04-12
We present for the first time a quantum mechanics/molecular mechanics scheme which combines quantum Monte Carlo with the reaction field of classical polarizable dipoles (QMC/MMpol). In our approach, the optimal dipoles are self-consistently generated at the variational Monte Carlo level and then used to include environmental effects in diffusion Monte Carlo. We investigate the performance of this hybrid model in describing the vertical excitation energies of prototypical small molecules solvated in water, namely, methylenecyclopropene and s-trans acrolein. Two polarization regimes are explored where either the dipoles are optimized with respect to the ground-state solute density (polGS) or different sets of dipoles are separately brought to equilibrium with the states involved in the electronic transition (polSS). By comparing with reference supermolecular calculations where both solute and solvent are treated quantum mechanically, we find that the inclusion of the response of the environment to the excitation of the solute leads to superior results than the use of a frozen environment (point charges or polGS), in particular, when the solute-solvent coupling is dominated by electrostatic effects which are well recovered in the polSS condition. QMC/MMpol represents therefore a robust scheme to treat important environmental effects beyond static point charges, combining the accuracy of QMC with the simplicity of a classical approach. PMID:26959751
Ganesh, P.; Kim, Jeongnim; Park, Changwon; Yoon, Mina; Reboredo, Fernando A.; Kent, Paul R. C.
2014-11-03
In highly accurate diffusion quantum Monte Carlo (QMC) studies of the adsorption and diffusion of atomic lithium in AA-stacked graphite are compared with van der Waals-including density functional theory (DFT) calculations. Predicted QMC lattice constants for pure AA graphite agree with experiment. Pure AA-stacked graphite is shown to challenge many van der Waals methods even when they are accurate for conventional AB graphite. Moreover, the highest overall DFT accuracy, considering pure AA-stacked graphite as well as lithium binding and diffusion, is obtained by the self-consistent van der Waals functional vdW-DF2, although errors in binding energies remain. Empirical approaches based onmore » point charges such as DFT-D are inaccurate unless the local charge transfer is assessed. Our results demonstrate that the lithium carbon system requires a simultaneous highly accurate description of both charge transfer and van der Waals interactions, favoring self-consistent approaches.« less
Quantum Monte Carlo algorithms for electronic structure at the petascale; the endstation project.
Kim, J; Ceperley, D M; Purwanto, W; Walter, E J; Krakauer, H; Zhang, S W; Kent, P.R. C; Hennig, R G; Umrigar, C; Bajdich, M; Kolorenc, J; Mitas, L; Srinivasan, A
2008-10-01
Over the past two decades, continuum quantum Monte Carlo (QMC) has proved to be an invaluable tool for predicting of the properties of matter from fundamental principles. By solving the Schrodinger equation through a stochastic projection, it achieves the greatest accuracy and reliability of methods available for physical systems containing more than a few quantum particles. QMC enjoys scaling favorable to quantum chemical methods, with a computational effort which grows with the second or third power of system size. This accuracy and scalability has enabled scientific discovery across a broad spectrum of disciplines. The current methods perform very efficiently at the terascale. The quantum Monte Carlo Endstation project is a collaborative effort among researchers in the field to develop a new generation of algorithms, and their efficient implementations, which will take advantage of the upcoming petaflop architectures. Some aspects of these developments are discussed here. These tools will expand the accuracy, efficiency and range of QMC applicability and enable us to tackle challenges which are currently out of reach. The methods will be applied to several important problems including electronic and structural properties of water, transition metal oxides, nanosystems and ultracold atoms.
Fermionic Quantum Monte Carlo simulations without fixed nodes
NASA Astrophysics Data System (ADS)
Dornheim, Tobias; Schoof, Tim; Groth, Simon; Bonitz, Michael
Recent restricted PIMC (RPIMC) simulations [PRL 110, 146405 (2013)] of the uniform electron gas (UEG) at finite temperature have turned out to be surprisingly inaccurate [PRL 115, 130402 (2015)]. Therefore, there exists a high need for alternative approaches which circumvent the fermion sign problem (FSP). In this work, we present two independent approaches which exhibit a complementary behavior. The configuration PIMC (CPIMC) method [Contrib. Plasma Phys. 51, 687-697 (2011)], which operates in Fock space, excels at high density and allows for cutting edge results at strong degeneracy. In contrast, the permutation blocking PIMC (PB-PIMC) approach [New J. Phys. 17, 073017 (2015)] is formulated in coordinate space and combines antisymmetric imaginary time propagators (determinants) with a higher order factorization of the density matrix. This leads to a significant reduction of the sign problem and extends the range of applicability of standard PIMC towards higher density and lower temperature [arXiv:1508.03221 (2015)]. Joining these two complementary methods allows us to present accurate thermodynamic results for the uniform electron gas over a broad parameter range and, therefore, to partly avoid the FSP.
Ising nematic quantum critical point in a metal: a Monte Carlo study
NASA Astrophysics Data System (ADS)
Lederer, Samuel
The Ising nematic quantum critical point (QCP) associated with the zero temperature transition from a symmetric to a nematic metal is an exemplar of metallic quantum criticality. We have carried out a minus sign-free quantum Monte Carlo study of this QCP for a two dimensional lattice model with sizes up to 24 × 24 sites. The system remains non-superconducting down to the lowest accessible temperatures. The results exhibit critical scaling behavior over the accessible ranges of temperature, (imaginary) time, and distance. This scaling behavior has remarkable similarities with recently measured properties of the Fe-based superconductors proximate to their putative nematic QCP. With Yoni Schattner, Steven A. Kivelson, and Erez Berg.
Sign-problem-free quantum Monte Carlo of the onset of antiferromagnetism in metals.
Berg, Erez; Metlitski, Max A; Sachdev, Subir
2012-12-21
The quantum theory of antiferromagnetism in metals is necessary for our understanding of numerous intermetallic compounds of widespread interest. In these systems, a quantum critical point emerges as external parameters (such as chemical doping) are varied. Because of the strong coupling nature of this critical point and the "sign problem" plaguing numerical quantum Monte Carlo (QMC) methods, its theoretical understanding is still incomplete. Here, we show that the universal low-energy theory for the onset of antiferromagnetism in a metal can be realized in lattice models, which are free from the sign problem and hence can be simulated efficiently with QMC. Our simulations show Fermi surface reconstruction and unconventional spin-singlet superconductivity across the critical point. PMID:23258893
NASA Astrophysics Data System (ADS)
Lu, Shih-I.
2004-12-01
For a test set of 17 first-row small molecules, the equilibrium structures are calculated with Ornstein-Uhlenbeck diffusion quantum Monte Carlo simulations guiding by trial wave functions constructed from floating spherical Gaussian orbitals and spherical Gaussian geminals. To measure performance of the Monte Carlo calculations, the mean deviation, the mean absolute deviation, the maximum absolute deviation, and the standard deviation of Monte Carlo calculated equilibrium structures with respect to empirical equilibrium structures are given. This approach is found to yield results having a uniformly high quality, being consistent with empirical equilibrium structures and surpassing calculated values from the coupled cluster model with single, double, and noniterative triple excitations [CCSD(T)] with the basis sets of cc-pCVQZ and cc-pVQZ. The nonrelativistic equilibrium atomization energies are also presented to assess performance of the calculated methods. The mean absolute deviations regarding experimental atomization energy are 0.16 and 0.21 kcal/mol for the Monte Carlo and CCSD(T)/cc-pCV(56)Z calculations, respectively.
Efficient continuous-time quantum Monte Carlo method for the ground state of correlated fermions
NASA Astrophysics Data System (ADS)
Wang, Lei; Iazzi, Mauro; Corboz, Philippe; Troyer, Matthias
2015-06-01
We present the ground state extension of the efficient continuous-time quantum Monte Carlo algorithm for lattice fermions of M. Iazzi and M. Troyer, Phys. Rev. B 91, 241118 (2015), 10.1103/PhysRevB.91.241118. Based on continuous-time expansion of an imaginary-time projection operator, the algorithm is free of systematic error and scales linearly with projection time and interaction strength. Compared to the conventional quantum Monte Carlo methods for lattice fermions, this approach has greater flexibility and is easier to combine with powerful machinery such as histogram reweighting and extended ensemble simulation techniques. We discuss the implementation of the continuous-time projection in detail using the spinless t -V model as an example and compare the numerical results with exact diagonalization, density matrix renormalization group, and infinite projected entangled-pair states calculations. Finally we use the method to study the fermionic quantum critical point of spinless fermions on a honeycomb lattice and confirm previous results concerning its critical exponents.
Geometrically Constructed Markov Chain Monte Carlo Study of Quantum Spin-phonon Complex Systems
NASA Astrophysics Data System (ADS)
Suwa, Hidemaro
2013-03-01
We have developed novel Monte Carlo methods for precisely calculating quantum spin-boson models and investigated the critical phenomena of the spin-Peierls systems. Three significant methods are presented. The first is a new optimization algorithm of the Markov chain transition kernel based on the geometric weight allocation. This algorithm, for the first time, satisfies the total balance generally without imposing the detailed balance and always minimizes the average rejection rate, being better than the Metropolis algorithm. The second is the extension of the worm (directed-loop) algorithm to non-conserved particles, which cannot be treated efficiently by the conventional methods. The third is the combination with the level spectroscopy. Proposing a new gap estimator, we are successful in eliminating the systematic error of the conventional moment method. Then we have elucidated the phase diagram and the universality class of the one-dimensional XXZ spin-Peierls system. The criticality is totally consistent with the J1 -J2 model, an effective model in the antiadiabatic limit. Through this research, we have succeeded in investigating the critical phenomena of the effectively frustrated quantum spin system by the quantum Monte Carlo method without the negative sign. JSPS Postdoctoral Fellow for Research Abroad
Quantum Monte-Carlo simulation of spin-one antiferromagnets with single-ion anisotropy
NASA Astrophysics Data System (ADS)
Kato, Yasuyuki; Wierschem, Keola; Nishida, Yusuke; Batista, Cristian; Sengupta, Pinaki
2013-03-01
We study a spin-one Heisenberg model with uniaxial single-ion anisotropy, D, and Zeeman coupling to a magnetic field, B, parallel to the symmetry axis. We compute the (D / J , B / J) quantum phase diagram for square and simple cubic lattices by combining analytical and Quantum Monte Carlo approaches, and find a transition between XY-antiferromagnetic and ferronematic phases that spontaneously break the U(1) symmetry of the model. In the language of bosonic gases, this is a transition between a Bose-Einstein condensate (BEC) of single bosons and a BEC of pairs. For the efficient simulation of ferronematic phase, we developed and implemented a new multi-discontinuity algorithm based on the directed-loop algorithm. The ordinary quantum Monte-Carlo methods fall into freezing problems when we apply them to this system at large D / J and finite B / J ~ 1 . The new method does not suffer from the freezing problems. This research used resources of the NERSCC (DOE Contract No. DE-AC02-05CH11231). Work at LANL was performed under the auspices of a J. Robert Oppenheimer Fellowship and the U.S. DOE contract No. DE-AC52-06NA25396 through the LDRD program.
Quantum Monte Carlo calculations of excited states in A=6-8 nuclei
Pieper, Steven C.; Wiringa, R.B.; Carlson, J.
2004-11-01
A variational Monte Carlo method is used to generate sets of orthogonal trial functions, {psi}{sub T}(J{sup {pi}};T), for given quantum numbers in various light p-shell nuclei. These {psi}{sub T} are then used as input to Green's function Monte Carlo (GFMC) calculations of first, second, and higher excited (J{sup {pi}};T) states. Realistic two- and three-nucleon interactions are used. We find that if the physical excited state is reasonably narrow, the GFMC energy converges to a stable result. With the combined Argonne v{sub 18} two-nucleon and Illinois-2 three-nucleon interactions, the results for many second and higher states in A=6-8 nuclei are close to the experimental values.
Electron density of states of Fe-based superconductors: Quantum trajectory Monte Carlo method
NASA Astrophysics Data System (ADS)
Kashurnikov, V. A.; Krasavin, A. V.; Zhumagulov, Ya. V.
2016-03-01
The spectral and total electron densities of states in two-dimensional FeAs clusters, which simulate iron-based superconductors, have been calculated using the generalized quantum Monte Carlo algorithm within the full two-orbital model. Spectra have been reconstructed by solving the integral equation relating the Matsubara Green's function and spectral density by the method combining the gradient descent and Monte Carlo algorithms. The calculations have been performed for clusters with dimensions up to 10 × 10 FeAs cells. The profiles of the Fermi surface for the entire Brillouin zone have been presented in the quasiparticle approximation. Data for the total density of states near the Fermi level have been obtained. The effect of the interaction parameter, size of the cluster, and temperature on the spectrum of excitations has been studied.
Twist-averaged boundary conditions in continuum quantum Monte Carlo algorithms
NASA Astrophysics Data System (ADS)
Lin, C.; Zong, F. H.; Ceperley, D. M.
2001-07-01
We develop and test Quantum Monte Carlo algorithms that use a``twist'' or a phase in the wave function for fermions in periodic boundary conditions. For metallic systems, averaging over the twist results in faster convergence to the thermodynamic limit than periodic boundary conditions for properties involving the kinetic energy and has the same computational complexity. We determine exponents for the rate of convergence to the thermodynamic limit for the components of the energy of coulomb systems. We show results with twist averaged variational Monte Carlo on free particles, the Stoner model and the electron gas using Hartree-Fock, Slater-Jastrow, and three-body and backflow wave function. We also discuss the use of twist averaging in the grand canonical ensemble, and numerical methods to accomplish the twist averaging.
NASA Astrophysics Data System (ADS)
Ge, Rong-Chun; Hughes, Stephen
2015-11-01
We study the quantum dynamics of two quantum dots (QDs) or artificial atoms coupled through the fundamental localized plasmon of a gold nanorod resonator. We derive an intuitive and efficient time-local master equation, in which the effect of the metal nanorod is taken into consideration self-consistently using a quasinormal mode (QNM) expansion technique of the photon Green function. Our efficient QNM technique offers an alternative and more powerful approach over the standard Jaynes-Cummings model, where the radiative decay, nonradiative decay, and spectral reshaping effect of the electromagnetic environment is rigorously included in a clear and transparent way. We also show how one can use our approach to compliment the approximate Jaynes-Cummings model in certain spatial regimes where it is deemed to be valid. We then present a study of the quantum dynamics and photoluminescence spectra of the two plasmon-coupled QDs. We first explore the non-Markovian regime, which is found to be important only on the ultrashort time scale of the plasmon mode which is about 40 fs. For the field free evolution case of excited QDs near the nanorod, we demonstrate how spatially separated QDs can be effectively coupled through the plasmon resonance and we show how frequencies away from the plasmon resonance can be more effective for coherently coupling the QDs. Despite the strong inherent dissipation of gold nanoresonators, we show that qubit entanglements as large as 0.7 can be achieved from an initially separate state, which has been limited to less than 0.5 in previous work for weakly coupled reservoirs. We also study the superradiance and subradiance decay dynamics of the QD pair. Finally, we investigate the rich quantum dynamics of QDs that are incoherently pumped, and study the polarization dependent behavior of the emitted photoluminescence spectrum where a double-resonance structure is observed due to the strong photon exchange interactions. Our general quantum plasmonics
NASA Astrophysics Data System (ADS)
Busemeyer, Brian; Dagrada, Mario; Sorella, Sandro; Casula, Michele; Wagner, Lucas K.
2016-07-01
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 Tc 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 analysis of the local 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 stripelike magnetic fluctuations may be responsible for the enhanced Tc in FeSe and that higher Tc is associated with nearness to a crossover between collinear and bicollinear ordering.
A Monte Carlo-quantum mechanics study of a solvatochromic π* probe.
Domínguez, Moisés; Rezende, Marcos Caroli
2016-09-01
The solvation and the solvatochromic behavior of 5-(dimethylamino)-5'-nitro-2,2'-bithiophene 1, the basis of a π* scale of solvent polarities, was investigated theoretically in toluene, dichloromethane, methanol and formamide with a Monte Carlo and quantum mechanics (QM/MM) iterative approach. The calculated transition energies of the solvatochromic band of 1, obtained as averages of statistically uncorrelated configurations, including the solute and explicit solvent molecules of the first solvation layer, besides showing good agreement with the experimental transitions, reproduced very well the positive solvatochromism of this probe in various solvents. PMID:27553303
Neutron matter with Quantum Monte Carlo: chiral 3N forces and static response
NASA Astrophysics Data System (ADS)
Buraczynski, M.; Gandolfi, S.; Gezerlis, A.; Schwenk, A.; Tews, I.
2016-03-01
Neutron matter is related to the physics of neutron stars and that of neutron-rich nuclei. Quantum Monte Carlo (QMC) methods offer a unique way of solving the many-body problem non-perturbatively, providing feedback on features of nuclear interactions and addressing scenarios that are inaccessible to other approaches. In this contribution we go over two recent accomplishments in the theory of neutron matter: a) the fusing of QMC with chiral effective field theory interactions, focusing on local chiral 3N forces, and b) the first attempt to find an ab initio solution to the problem of static response.
NASA Astrophysics Data System (ADS)
Calcavecchia, Francesco; Holzmann, Markus
2016-04-01
We use the shadow wave function formalism as a convenient model to study the fermion sign problem affecting all projector quantum Monte Carlo methods in continuum space. We demonstrate that the efficiency of imaginary-time projection algorithms decays exponentially with increasing number of particles and/or imaginary-time propagation. Moreover, we derive an analytical expression that connects the localization of the system with the magnitude of the sign problem, illustrating this behavior through numerical results. Finally, we discuss the computational complexity of the fermion sign problem and methods for alleviating its severity.
Calcavecchia, Francesco; Holzmann, Markus
2016-04-01
We use the shadow wave function formalism as a convenient model to study the fermion sign problem affecting all projector quantum Monte Carlo methods in continuum space. We demonstrate that the efficiency of imaginary-time projection algorithms decays exponentially with increasing number of particles and/or imaginary-time propagation. Moreover, we derive an analytical expression that connects the localization of the system with the magnitude of the sign problem, illustrating this behavior through numerical results. Finally, we discuss the computational complexity of the fermion sign problem and methods for alleviating its severity. PMID:27176442
Theory of melting at high pressures: Amending density functional theory with quantum Monte Carlo
NASA Astrophysics Data System (ADS)
Shulenburger, L.; Desjarlais, M. P.; Mattsson, T. R.
2014-10-01
We present an improved first-principles description of melting under pressure based on thermodynamic integration comparing density functional theory (DFT) and quantum Monte Carlo (QMC) treatments. The method is applied to address the longstanding discrepancy between DFT calculations and diamond anvil cell (DAC) experiments on the melting curve of xenon, a noble gas solid where van der Waals binding is challenging for traditional DFT methods. The calculations show agreement with data below 20 GPa and that the high-pressure melt curve is well described by a Lindemann behavior up to at least 80 GPa, in contrast to DAC data.
Neutron matter with Quantum Monte Carlo: chiral 3N forces and static response
Buraczynski, M.; Gandolfi, S.; Gezerlis, A.; Schwenk, A.; Tews, I.
2016-03-01
Neutron matter is related to the physics of neutron stars and that of neutron-rich nuclei. Moreover, Quantum Monte Carlo (QMC) methods offer a unique way of solving the many-body problem non-perturbatively, providing feedback on features of nuclear interactions and addressing scenarios that are inaccessible to other approaches. Our contribution goes over two recent accomplishments in the theory of neutron matter: a) the fusing of QMC with chiral effective field theory interactions, focusing on local chiral 3N forces, and b) the first attempt to find an ab initio solution to the problem of static response.
Theory of finite size effects for electronic quantum Monte Carlo calculations of liquids and solids
NASA Astrophysics Data System (ADS)
Holzmann, Markus; Clay, Raymond C.; Morales, Miguel A.; Tubman, Norm M.; Ceperley, David M.; Pierleoni, Carlo
2016-07-01
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 effective procedures, such as using the potential and wave function split up into long and short range functions to simplify the method, and we discuss how to treat backflow wave functions. 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.