Many-Body Convergence of the Electrostatic Properties of Water.

J Chem Theory Comput

Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92093, United States.

Published: November 2013

The many-body convergence of the dipole moment and the dipole-dipole polarizability of water is investigated. It is found that, for systems of low symmetry like the water clusters examined here, simple measures such as dipole magnitudes and average polarizabilities may lead to an incomplete interpretation of the underlying physics. Alternative metrics are introduced that allow for an unambiguous characterization of both properties. The convergence of the many-body decomposition of the total dipole and the polarizability is studied for (H2O)N, with N = 4-6 being minimum energy water clusters and N = 14 being clusters that were extracted from condensed phase simulations. For these clusters, it is demonstrated that both the total dipole and polarizability are almost entirely pairwise additive, with three-body terms contributing less than 4% and all higher-order terms being essentially negligible.

Download full-text PDF

Source
http://dx.doi.org/10.1021/ct400696dDOI Listing

Publication Analysis

Top Keywords

many-body convergence
8
water clusters
8
total dipole
8
dipole polarizability
8
convergence electrostatic
4
electrostatic properties
4
water
4
properties water
4
water many-body
4
dipole
4

Similar Publications

Green's function theory has emerged as a powerful many-body approach not only in condensed matter physics but also in quantum chemistry in recent years. We have developed a new all-electron implementation of the BSE@GW formalism using numeric atom-centered orbital basis sets (Liu, C. 2020, 152, 044105).

View Article and Find Full Text PDF

Electronic structure calculations in enzymes converge very slowly with respect to the size of the model region that is described using quantum mechanics (QM), requiring hundreds of atoms to obtain converged results and exhibiting substantial sensitivity (at least in smaller models) to which amino acids are included in the QM region. As such, there is considerable interest in developing automated procedures to construct a QM model region based on well-defined criteria. However, testing such procedures is burdensome due to the cost of large-scale electronic structure calculations.

View Article and Find Full Text PDF

Fragment-based quantum chemistry methods offer a means to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein-ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost.

View Article and Find Full Text PDF

The emerging field of molecular cavity polaritons has stimulated a surge of experimental and theoretical activities and presents a unique opportunity to develop the many-body simulation methodology. This paper presents a numerical scheme for the extraction of key kinetic information of lossy cavity polaritons based on the transfer tensor method (TTM). Steady state, relaxation timescales, and oscillatory phenomena can all be deduced directly from a set of transfer tensors without the need for long-time simulation.

View Article and Find Full Text PDF

Convergence of the many-body expansion with respect to distance cutoffs in crystals of polar molecules: Acetic acid, formamide, and imidazole.

J Chem Phys

December 2024

Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA.

The many-body expansion, where one computes the total energy of a supersystem as the sum of the dimer, trimer, tetramer, etc., subsystems, provides a convenient approach to compute the lattice energies of molecular crystals. We investigate approximate methods for computing the non-additive three-body contributions to the crystal lattice energy of the polar molecules acetic acid, imidazole, and formamide, comparing to coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] level benchmarks.

View Article and Find Full Text PDF

Want AI Summaries of new PubMed Abstracts delivered to your In-box?

Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!