Search

Data structures and algorithms are essential building blocks for programs, and \emph{distributed data structures}, which automatically partition data across multiple memory locales, are essential to writing high-level parallel programs. While many projects have designed and implemented C++ distributed data structures and algorithms, there has not been widespread adoption of an interoperable model allowing algorithms and data structures from different libraries to work together. This paper introduces distributed ranges, which is a model for building generic data structures, views, and algorithms. A distributed range extends a C++ range, which is an iterable sequence of values, with a concept of segmentation, thus exposing how the distributed range is partitioned over multiple memory locales. Distributed data structures provide this distributed range interface, which allows them to be used with a collection of generic algorithms implemented using the distributed range interface. The modular nature of the model allows for the straightforward implementation of \textit{distributed views}, which are lightweight objects that provide a lazily evaluated view of another range. Views can be composed together recursively and combined with algorithms to implement computational kernels using efficient, flexible, and high-level standard C++ primitives. We evaluate the distributed ranges model by implementing a set of standard concepts and views as well as two execution runtimes, a multi-node, MPI-based runtime and a single-process, multi-GPU runtime. We demonstrate that high-level algorithms implemented using generic, high-level distributed ranges can achieve performance competitive with highly-tuned, expert-written code., Comment: To appear in ACM International Conference on Supercomputing (ICS) 2024

QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org .

QMCPACK has enabled cutting-edge materials research on supercomputers for over a decade. It scales nearly ideally but has low single-node efficiency due to the physics-based abstractions using array-of-structures objects, causing inefficient vectorization. We present a systematic approach to transform QMCPACK to better exploit the new hardware features of modern CPUs in portable and maintainable ways. We develop miniapps for fast prototyping and optimizations. We implement new containers in structure-of-arrays data layout to facilitate vectorizations by the compilers. Further speedup and smaller memory-footprints are obtained by computing data on the fly with the vectorized routines and expanding single-precision use. All these are seamlessly incorporated in production QMCPACK. We demonstrate upto 4.5x speedups on recent Intel processors and IBM Blue Gene/Q for representative workloads. Energy consumption is reduced significantly commensurate to the speedup factor. Memory-footprints are reduced by up-to 3.8x, opening the possibility to solve much larger problems of future., Comment: 12 pages, 10 figures, 2 tables, to be published at SC17

$\alpha$-graphyne is a two-dimensional sheet of $sp$-$sp^2$ hybridized carbon atoms in a honeycomb lattice. While the geometrical structure is similar to that of graphene, the hybridized triple bonds give rise to electronic structure that is different from that of graphene. Similar to graphene, $\alpha$-graphyne can be stacked in bilayers with two stable configurations, but the different stackings have very different electronic structures: one is predicted to have gapless parabolic bands and the other a tunable band gap which is attractive for applications. In order to realize applications, it is crucial to understand which stacking is more stable. This is difficult to model, as the stability is a result of weak interlayer van der Waals interactions which are not well captured by density functional theory (DFT). We have used quantum Monte Carlo simulations that accurately include van der Waals interactions to calculate the interlayer binding energy of bilayer graphyne and to determine its most stable stacking mode. Our results show that interlayer bindings of $sp$- and $sp^{2}$-bonded carbon networks are significantly underestimated in a Kohn-Sham DFT approach, even with an exchange-correlation potential corrected to include, in some approximation, van der Waals interactions. Finally, our quantum Monte Carlo calculations reveal that the interlayer binding energy difference between the two stacking modes is only 0.9(4) meV/atom. From this we conclude that the two stable stacking modes of bilayer $\alpha$-graphyne are almost degenerate with each other, and both will occur with about the same probability at room temperature unless there is a synthesis path that prefers one stacking over the other., Comment: 25 pages, 6 figures

B-spline based orbital representations are widely used in Quantum Monte Carlo (QMC) simulations of solids, historically taking as much as 50% of the total run time. Random accesses to a large four-dimensional array make it challenging to efficiently utilize caches and wide vector units of modern CPUs. We present node-level optimizations of B-spline evaluations on multi/many-core shared memory processors. To increase SIMD efficiency and bandwidth utilization, we first apply data layout transformation from array-of-structures to structure-of-arrays (SoA). Then by blocking SoA objects, we optimize cache reuse and get sustained throughput for a range of problem sizes. We implement efficient nested threading in B-spline orbital evaluation kernels, paving the way towards enabling strong scaling of QMC simulations. These optimizations are portable on four distinct cache-coherent architectures and result in up to 5.6x performance enhancements on Intel Xeon Phi processor 7250P (KNL), 5.7x on Intel Xeon Phi coprocessor 7120P, 10x on an Intel Xeon processor E5v4 CPU and 9.5x on BlueGene/Q processor. Our nested threading implementation shows nearly ideal parallel efficiency on KNL up to 16 threads. We employ roofline performance analysis to model the impacts of our optimizations. This work combined with our current efforts of optimizing other QMC kernels, result in greater than 4.5x speedup of miniQMC on KNL., Comment: 11 pages, 10 figures, 4 tables

We have performed quantum Monte Carlo (QMC) simulations and density functional theory (DFT) calculations to study the equations of state of MgSiO$_3$ perovskite (Pv) and post-perovskite (PPv), up to the pressure and temperature conditions of the base of Earth's lower mantle. The ground state energies were derived using QMC and the temperature dependent Helmholtz free energies were calculated within the quasi-harmonic approximation and density functional perturbation theory. The equations of state for both phases of MgSiO$_3$ agree well with experiments, and better than those from generalized gradient approximation (GGA) calculations. The Pv-PPv phase boundary calculated from our QMC equations of states is also consistent with experiments, and better than previous LDA calculations. We discuss the implications for double crossing of the Pv-PPv boundary in the Earth.

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 implemented 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 shows a quantitative 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., Comment: 9 pages, 5 figures

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 O$_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., Comment: 43 pages, 10 figures

We have evaluated the successes and failures of the Hubbard-corrected density functional theory (DFT+U) approach to study Mg doping of LiCoO$_2$. We computed the effect of the U parameter on the energetic, geometric and electronic properties of two possible doping mechanisms: (1) substitution of Mg onto a Co (or Li) site with an associated impurity state and, (2) formation of impurity-state-free complexes of substitutional Mg and point defects in LiCoO$_2$. We find that formation of impurity states results in changes on the valency of Co in LiCoO$_2$. Variation of the Co U shifts the energy of the impurity state, resulting in energetic, geometric and electronic properties that depend significantly on the specific value of U. In contrast, the properties of the impurity-state-free complexes are insensitive to U. These results identify reasons for the strong dependence on the doping properties on the chosen value of U and for the overall difficulty of achieving agreement with the experimentally known energetic and electronic properties of doped transition metal oxides such as LiCoO$_2$., Comment: 12 pages, 6 figures

In view of the continuous theoretical efforts aimed at an accurate microscopic description of the strongly correlated transition metal oxides and related materials, we show that with continuum quantum Monte Carlo (QMC) calculations it is possible to obtain the value of the spin superexchange coupling constant of a copper oxide in a quantitatively excellent agreement with experiment. The variational nature of the QMC total energy allows us to identify the best trial wave function out of the available pool of wave functions, which makes the approach essentially free from adjustable parameters and thus truly ab initio. The present results on magnetic interactions suggest that QMC is capable of accurately describing ground state properties of strongly correlated materials., Comment: Published in Physical Review X

We have performed quantum Monte Carlo calculations to study the cohesion energetics of carbon allotropes, including $sp^3$-bonded diamond, $sp^2$-bonded graphene, $sp$-$sp^2$ hybridized graphynes, and $sp$-bonded carbyne. The computed cohesive energies of diamond and graphene are found to be in excellent agreement with the corresponding values determined 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 systematically as the ratio of $sp$-bonded carbon atoms increases. The cohesive energy of $\gamma$-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. Experimental difficulty in synthesizing graphynes could be explained by their significantly smaller cohesive energies. Finally we conclude that the cohesive energy of a newly-proposed graphyne can be accurately estimated with the carbon-carbon bond energies determined from the cohesive energies of graphene and three different graphynes considered here., Comment: 18 pages, 4 figures

Solid atomic hydrogen is one of the simplest systems to undergo a metal-insulator transition. Near the transition, the electronic degrees of freedom become strongly correlated and their description provides a difficult challenge for theoretical methods. As a result, the order and density of the phase transition are still subject to debate. In this work we use diffusion quantum Monte Carlo to benchmark the transition between paramagnetic and anti-ferromagnetic body centered cubic atomic hydrogen in its ground state. We locate the density of the transition by computing the equation of state for these two phases and identify the phase transition order by computing the band gap near the phase transition. These benchmark results show that the phase transition is continuous and occurs at a Wigner-Seitz radius of $r_s=2.27(3) a_0$. We compare our results to previously reported density functional theory, Hedin's GW approximation, and dynamical mean field theory results., Comment: 7 pages, 3 figures

A statistical method is derived for the calculation of thermodynamic properties of many-body systems at low temperatures. This method is based on the self-healing diffusion Monte Carlo method for complex functions [F. A. Reboredo J. Chem. Phys. 136, 204101 (2012)] and some ideas of the correlation function Monte Carlo approach [D. M. Ceperley and B. Bernu, J. Chem. Phys. 89, 6316 (1988)]. In order to allow the evolution in imaginary time to describe the density matrix, we remove the fixed-node restriction using complex antisymmetric guiding wave functions. In the process we obtain a parallel algorithm that optimizes a small subspace of the many-body Hilbert space to have maximum overlap with the subspace spanned by the lowest-energy eigenstates of a many-body Hamiltonian. We show in a model system that the partition function is progressively maximized within this subspace. We show that the subspace spanned by the small basis systematically converges towards the subspace spanned by the lowest energy eigenstates. Possible applications of this method to calculate the thermodynamic properties of many-body systems near the ground state are discussed. The resulting basis can be also used to accelerate the calculation of the ground or excited states with Quantum Monte Carlo.

We establish a physically meaningful representation of a quantum energy density for use in Quantum Monte Carlo calculations. The energy density operator, defined in terms of Hamiltonian components and density operators, returns the correct Hamiltonian when integrated over a volume containing a cluster of particles. This property is demonstrated for a helium-neon "gas," showing that atomic energies obtained from the energy density correspond to eigenvalues of isolated systems. The formation energies of defects or interfaces are typically calculated as total energy differences. Using a model of delta-doped silicon (where dopant atoms form a thin plane) we show how interfacial energies can be calculated more efficiently with the energy density, since the region of interest is small. We also demonstrate how the energy density correctly transitions to the bulk limit away from the interface where the correct energy is obtainable from a separate total energy calculation., Comment: 11 pages, 4 figures

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

We study the electronic structure of a spherical jellium in the presence of a central Gaussian impurity. We test how well the resulting inhomogeneity effects beyond spherical jellium are reproduced by several approximations of density functional theory (DFT). Four rungs of Perdew's ladder of DFT functionals, namely local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA and orbital-dependent hybrid functionals are compared against our quantum Monte Carlo (QMC) benchmarks. We identify several distinct transitions in the ground state of the system as the electronic occupation changes between delocalized and localized states. We examine the parameter space of realistic densities ($1 \le r_s\le 5$) and moderate depths of the Gaussian impurity ($Z<7$). The selected 18 electron system (with closed-shell ground state) presents $1d \to 2s$ transitions while the 30 electron system (with open-shell ground state) exhibits $1f \to 2p$ transitions. For the former system, the accuracy for the transitions is clearly improving with increasing sophistication of functionals with meta-GGA and hybrid functionals having only small deviations from QMC. However, for the latter system, we find much larger differences for the underlying transitions between our pool of DFT functionals and QMC. We attribute this failure to treatment of the exact exchange within these functionals. Additionally, we amplify the inhomogeneity effects by creating the system with spherical shell which leads to even larger errors in DFT approximations., Comment: 8 pages, 4 figures, submitted to PRB as a regular article revisited version after review

Quantum Monte Carlo (QMC) methods are often used to calculate properties of many body quantum systems. The main cost of many QMC methods, for example the variational Monte Carlo (VMC) method, is in constructing a sequence of Slater matrices and computing the ratios of determinants for successive Slater matrices. Recent work has improved the scaling of constructing Slater matrices for insulators so that the cost of constructing Slater matrices in these systems is now linear in the number of particles, whereas computing determinant ratios remains cubic in the number of particles. With the long term aim of simulating much larger systems, we improve the scaling of computing the determinant ratios in the VMC method for simulating insulators by using preconditioned iterative solvers. The main contribution of this paper is the development of a method to efficiently compute for the Slater matrices a sequence of preconditioners that make the iterative solver converge rapidly. This involves cheap preconditioner updates, an effective reordering strategy, and a cheap method to monitor instability of ILUTP preconditioners. Using the resulting preconditioned iterative solvers to compute determinant ratios of consecutive Slater matrices reduces the scaling of QMC algorithms from O(n^3) per sweep to roughly O(n^2), where n is the number of particles, and a sweep is a sequence of n steps, each attempting to move a distinct particle. We demonstrate experimentally that we can achieve the improved scaling without increasing statistical errors. Our results show that preconditioned iterative solvers can dramatically reduce the cost of VMC for large(r) systems., Comment: 24 pages, 10 figures

We present 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 derive its discontinuity at the Fermi wavevector. This yields an accurate measure of the renormalization factor that we compare with quantum Monte Carlo and GW calculations performed both on crystalline sodium and on the homogeneous electron gas. Our calculated results are in good agreement with the experiment., Comment: 5 pages, 3 figures, 1 table

We develop an all-electron quantum Monte Carlo (QMC) method for solids that does not rely on pseudopotentials, and use it to construct a primary ultra-high pressure calibration based the equation of state of cubic boron nitride(c-BN). We compute the static contribution to the free energy with QMC, and obtain the phonon contribution from density functional theory, yielding a high-accuracy calibration up to 900 GPa usable directly in experiment. Furthermore, we compute the anharmonic Raman frequency shift with QMC as a function of pressure and temperature, allowing optical pressure calibration in table-top experiments. 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 methodology is generally applicable to first-row solids, and can be used to provide a new reference for ab initio calculations of solids and to benchmark pseudopotential accuracy., Comment: 4 pages, 2 figures, 1 table

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 mHa, suggesting 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 empirical models and previous calculations, and we discuss the physical mechanisms of the interaction, the role of diffuse basis functions, and the charge redistribution in the bond., Comment: 12 pages, 4 figures

The morphology of nanoscopic Ag grains significantly affects the phonons. Atomistic simulations show that realistic nanograin models display complex vibrational properties. (1) Single-crystalline grains. Nearly-pure torsional and radial phonons appear at low frequencies. For low-energy, faceted models, the breathing mode and acoustic gap (lowest frequency) are about 10% lower than predicted by elasticity theory (ET) for a continuum sphere of the same volume. The sharp edges and the atomic lattice split the ET-acoustic-gap quintet into a doublet and triplet. The surface protrusions associated with nearly spherical, high-energy models produce a smaller acoustic gap and a higher vibrational density of states (DOS) at frequencies \nu<2 THz. (2) Twined icosahedra. In contrast to the single-crystal case, the inherent strain produce a larger acoustic gap, while the core atoms yield a DOS tail extending beyond the highest frequency of single-crystalline grains. (3) Mark's decahedra, in contrast to (1) and (2), do not have a breathing mode; although twined and strained, do not exhibit a high-frequency tail in the DOS. (4) Irregular nanograins. Grain boundaries and surface disorder yield non-degenerate phonon frequencies, and significantly smaller acoustic gap. Only these nanograins exhibit a low-frequency \nu^2 DOS in the interval 1-2 THz., Comment: Version published in Phys. Rev. B

We present {\it ab initio} density-functional calculations of molecular solids formed from C$_{28}$-derived closed-shell fullerenes. Solid C$_{28}$H$_4$ is found to bind weakly and exhibits many of the electronic structure features of solid C$_{60}$ with an enhanced electron-phonon interaction potential. We show that chemical doping of this structure is feasible, albeit more restrictive than its C$_{60}$ counterpart, with an estimated superconducting transition temperature exceeding those of the alkali-doped C$_{60}$ solids., Comment: Lower quality postscript file for Figure 1 is used in the manuscript in order to meet submission quota for pre-print server. Higher quality postscript file available from author: naromero@alum.mit.edu This article has been updated to reflect changes incorporated during the peer review process. It is published in PRB 70, 140504(R) 2004

We have contrasted the quantum confinement of (i) multiple quantum wells of flat GaAs and AlAs layers, i.e. $(\GaAs)_{m}/(\AlAs)_n/(\GaAs)_p/(\AlAs)_q$, with (ii) ``cylindrical Russian Dolls'' -- an equivalent sequence of wells and barriers arranged as concentric wires. Using a pseudopotential plane-wave calculation, we identified theoretically a set of numbers ($m,n,p$ and $q$) such that charge separation can exist in ``cylindrical Russian Dolls'': the CBM is localized in the inner GaAs layer, while the VBM is localized in the outer GaAs layer., Comment: latex, 8 pages

We perform total energy calculations based on the tight-binding Hamiltonian scheme (i) to study the structural properties and energetics of the extended {311} defects depending upon their dimensions and interstitial concentrations and (ii) to find possible mechanisms of interstitial capture by and release from the {311} defects. The generalized orbital-based linear-scaling method implemented on Cray-T3D is used for supercell calculations of large scale systems containing more than 1000 Si atoms., Comment: 31 pages, LaTeX file, and 8 figures in postscript files, submitted to PRB

An energy functional for orbital based $O(N)$ calculations is proposed, which depends on a number of non orthogonal, localized orbitals larger than the number of occupied states in the system, and on a parameter, the electronic chemical potential, determining the number of electrons. We show that the minimization of the functional with respect to overlapping localized orbitals can be performed so as to attain directly the ground state energy, without being trapped at local minima. The present approach overcomes the multiple minima problem present within the original formulation of orbital based $O(N)$ methods; it therefore makes it possible to perform $O(N)$ calculations for an arbitrary system, without including any information about the system bonding properties in the construction of the input wavefunctions. Furthermore, while retaining the same computational cost as the original approach, our formulation allows one to improve the variational estimate of the ground state energy, and the energy conservation during a molecular dynamics run. Several numerical examples for surfaces, bulk systems and clusters are presented and discussed., Comment: 24 pages, RevTex file, 5 figures available upon request