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 wavefunctions 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 and graphical processing unit 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://qmcpack.org.

Reaction rates of common organic reactions have been reported to increase by one to six orders of magnitude in aqueous microdroplets compared to bulk solution, but the reasons for the rate acceleration are poorly understood. Using a coarse-grained electron model that describes structural organization and electron densities for water droplets without the expense of ab initio methods, we investigate the electric field distributions at the air-water interface to understand the origin of surface reactivity. We find that electric field alignments along free O–H bonds at the surface are ~16 MV/cm larger on average than that found for O–H bonds in the interior of the water droplet. Furthermore, electric field distributions can be an order of magnitude larger than the average due to non-linear coupling of intramolecular solvent polarization with intermolecular solvent modes which may contribute to even greater surface reactivity for weakening or breaking chemical bonds at the droplet surface.

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

Nanoconfined aqueous environments and the recent advent of accelerated chemistry in microdroplets are increasingly being investigated for catalysis. The mechanisms underlying the enhanced reactivity in alternate solvent environments, and whether the enhanced reactivity due to nanoconfinement is a universal phenomenon, are not fully understood. Here, we use ab initio molecular dynamics simulations to characterize the free energy of a retro-Diels-Alder reaction in bulk water at very different densities and in water nanoconfined by parallel graphene sheets. We find that the broadly different global solvation environments accelerate the reactions to a similar degree with respect to the gas phase reaction, with activation free energies that do not differ by more than kbT from each other.The reason for the same acceleration factor in the extremely different solvation environments is that it is the microsolvation of the dienophile's carbonyl group that governs the transition state stabilization and mechanism, which is not significantly disrupted by either the lower density in bulk water or the strong nanoconfinement conditions used here. Our results also suggest that significant acceleration of Diels Alder reactions in microdroplets or on-water conditions can't arise from local microsolvation when water is present, but instead must come from highly altered reaction environments that drastically change the reaction mechanisms.

Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to non-reactive force fields. In this Perspective we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid non-reactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.

We report a new deep learning message passing network that takes inspiration from Newton's equations of motion to learn interatomic potentials and forces. With the advantage of directional information from...

There is accumulating evidence that many chemical reactions are accelerated by several orders of magnitude in micrometer-sized aqueous or organic liquid droplets compared to their corresponding bulk liquid phase. However, the molecular origin of the enhanced rates remains unclear as in the case of spontaneous appearance of 1 μM hydrogen peroxide in water microdroplets. In this Letter, we consider the range of ionization energies and whether interfacial electric fields of a microdroplet can feasibly overcome the high energy step from hydroxide ions (OH − ) to hydroxyl radicals (OH • ) in a primary H 2 O 2 mechanism. We find that the vertical ionization energies (VIEs) of partially solvated OH − ions are greatly lowered relative to the average VIE in the bulk liquid, unlike the case of the Cl − anion which shows no reduction in the VIEs regardless of solvation environment. Overall reduced hydrogen-bonding and undercoordination of OH − are structural features that are more readily present at the air−water interface, where the energy scale for ionization can be matched by statistically probable electric field values.

Non-reactive force fields are defined by perturbations of electron density that are relatively small, whereas chemical reactivity involves wholesale electronic rearrangements that make and break bonds.Thus reactive force fields are incredibly difficult to develop compared to non-reactive force fields, and yet at the same time they fill a critical need when ab initio molecular dynamics methods are not affordable. We introduce a new reactive force field model for water that combines modified nonbonded terms of the ReaxFF model and its embedding in the electrostatic interactions described by our recently introduced coarse-grained electron model (C-GeM). The ReaxFF/C-GeM force field is characterized for many energetic and dissociative water properties for water clusters, structure and dynamical properties at ambient conditions in the condensed phase, as well the temperature dependence of density and water diffusion, with very good agreement with experiment. The ReaxFF/C-GeM should be more transferable and more broadly applicable to a range of reactive systems involving both proton and electron transfer in the condensed phase.