Development of an Advanced Force Field for Water using Variational Energy Decomposition Analysis

Given the piecewise approach to modeling intermolecular interactions for force fields, they can be difficult to parameterize since they are fit to data like total energies that only indirectly connect to their separable functional forms. Furthermore, by neglecting certain types of molecular interactions such as charge penetration and charge transfer, most classical force fields must rely on, but do not always demonstrate, how cancellation of errors occurs among the remaining molecular interactions accounted for such as exchange repulsion, electrostatics, and polarization. In this work we present the first generation of the (many-body) MB-UCB force field that explicitly accounts for the decomposed molecular interactions commensurate with a variational energy decomposition analysis, including charge transfer, with force field design choices that reduce the computational expense of the MB-UCB potential while remaining accurate. We optimize parameters using only single water molecule and water cluster data up through pentamers, with no fitting to condensed phase data, and we demonstrate that high accuracy is maintained when the force field is subsequently validated against conformational energies of larger water cluster data sets, radial distribution functions of the liquid phase, and the temperature dependence of thermodynamic and transport water properties. We conclude that MB-UCB is comparable in performance to MB-Pol, but is less expensive and more transferable by eliminating the need to represent short-ranged interactions through large parameter fits to high order polynomials

Parallel algorithms for solving linear equations using Givens transformations

AbstractThe use of the Givens method to solve linear equations on a parallel computer is reviewed, and a new algorithm which requires fewer time steps in the infinite processor case is presented

Cluster decomposition of full configuration interaction wave functions: a tool for chemical interpretation of systems with strong correlation

Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation amplitudes contain information on the complexity of the electronic wave function, but this information is contaminated by contributions from disconnected excitations, i.e. those excitations that are just products of independent lower-level excitations. The unwanted contributions can be removed via a cluster decomposition procedure, making it possible to examine the importance of connected excitations in complicated multireference molecules which are outside the reach of conventional algorithms. We present an implementation of the cluster decomposition analysis and apply it to both true FCI wave functions, as well as wave functions generated from the adaptive sampling CI (ASCI) algorithm. The cluster decomposition is useful for interpreting calculations in chemical studies, as a diagnostic for the convergence of various excitation manifolds, as well as as a guidepost for polynomially scaling electronic structure models. Applications are presented for (i) the double dissociation of water, (ii) the carbon dimer, (iii) the {\pi} space of polyacenes, as well as (iv) the chromium dimer. While the cluster amplitudes exhibit rapid decay with increasing rank for the first three systems, even connected octuple excitations still appear important in Cr$_2$, suggesting that spin-restricted single-reference coupled-cluster approaches may not be tractable for some problems in transition metal chemistry.Comment: 15 pages, 5 figure

A deterministic alternative to the full configuration interaction quantum Monte Carlo method

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 to achieve chemical accuracy for a wide range of systems, including the difficult Cr$_{2}$ dimer. In addition our method also allows efficient calculation of excited state energies, for which we illustrate with benchmark results for the excited states of C$_{2}$.Comment: 4 pages, 2 figure

Electron-vibration interaction in transport through atomic gold wires

We calculate the effect of electron-vibration coupling on conduction through atomic gold wires, which was measured in the experiments of Agra\"it et al. [Phys. Rev. Lett. 88, 216803 (2002)]. The vibrational modes, the coupling constants, and the inelastic transport are all calculated using a tight-binding parametrization and the non-equilibrium Green function formalism. The electron-vibration coupling gives rise to small drops in the conductance at voltages corresponding to energies of some of the vibrational modes. We study systematically how the position and height of these steps vary as a linear wire is stretched and more atoms are added to it, and find a good agreement with the experiments. We also consider two different types of geometries, which are found to yield qualitatively similar results. In contrast to previous calculations, we find that typically there are several close-lying drops due to different longitudinal modes. In the experiments, only a single drop is usually visible, but its width is too large to be accounted for by temperature. Therefore, to explain the experimental results, we find it necessary to introduce a finite broadening to the vibrational modes, which makes the separate drops merge into a single, wide one. In addition, we predict how the signatures of vibrational modes in the conductance curves differ between linear and zigzag-type wires.Comment: 19 pages, 12 figure

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method.

Recent advances in selected configuration interaction methods have made them competitive with the most accurate techniques available and, hence, creating an increasingly powerful tool for solving quantum Hamiltonians. In this work, we build on recent advances from the adaptive sampling configuration interaction (ASCI) algorithm. We show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The ASCI search algorithm can be used in several different modes, which includes an integral driven search and a coefficient driven search. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques, we present ASCI results applied to a large range of systems and basis sets to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-ζ benchmark data for the G1 data set. The largest of these calculations is Si2H6 which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets

The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients.

We present the theory and implementation of a Poisson-Boltzmann implicit solvation model for electrolyte solutions. This model can be combined with arbitrary electronic structure methods that provide an accurate charge density of the solute. A hierarchy of approximations for this model includes a linear approximation for weak electrostatic potentials, finite size of the mobile electrolyte ions, and a Stern-layer correction. Recasting the Poisson-Boltzmann equations into Euler-Lagrange equations then significantly simplifies the derivation of the free energy of solvation for these approximate models. The parameters of the model are either fit directly to experimental observables-e.g., the finite ion size-or optimized for agreement with experimental results. Experimental data for this optimization are available in the form of Sechenov coefficients that describe the linear dependence of the salting-out effect of solutes with respect to the electrolyte concentration. In the final part, we rationalize the qualitative disagreement of the finite ion size modification to the Poisson-Boltzmann model with experimental observations by taking into account the electrolyte concentration dependence of the Stern layer. A route toward a revised model that captures the experimental observations while including the finite ion size effects is then outlined. This implementation paves the way for the study of electrochemical and electrocatalytic processes of molecules and cluster models with accurate electronic structure methods

Zeolite-Catalyzed Isobutene Amination: Mechanism and Kinetics

Amination of isobutene with NH was investigated over Brønsted acidic zeolites at 1 atm and 453-483 K. To compare catalytic activities over different zeolites, the measured reaction rates are normalized by the number of active sites determined by tert-butylamine temperature-programmed desorption (TPD). Small- A nd medium-pore zeolites with one-dimensional channels exhibit low activity because of pore blockage by adsorbed tert-butylammonium ions. However, turnover frequencies and activation energies are not sensitive to framework identity, as long as the active site is accessible to isobutene and tert-butylamine. Kinetic measurements and FTIR spectroscopy reveal that the Brønsted acid sites in MFI are covered predominantly with tert-butylammonium ions under reaction conditions. The desorption of tert-butylamine is assisted by the concurrent adsorption of isobutene. DFT simulations show that at very low tert-butylamine partial pressures, for example, at the inlet to the reactor, tert-butylamine desorption is rate-limiting. However, at sufficiently high tert-butylamine partial pressures (>0.03 kPa), protonation of isobutene to the corresponding carbenium ion limits the rate of amination.