Electric field-gradient contributions to the chemical shifts of liquid water

MOLSIM: A modular molecular simulation software

The modular software MOLSIM for all-atom molecular and coarse-grained simulations is presented with focus on the underlying concepts used. The software possesses four unique features: (1) it is an integrated software for molecular dynamic, Monte Carlo, and Brownian dynamics simulations; (2) simulated objects are constructed in a hierarchical fashion representing atoms, rigid molecules and colloids, flexible chains, hierarchical polymers, and cross-linked networks; (3) long-range interactions involving charges, dipoles and/or anisotropic dipole polarizabilities are handled either with the standard Ewald sum, the smooth particle mesh Ewald sum, or the reaction-field technique; (4) statistical uncertainties are provided for all calculated observables. In addition, MOLSIM supports various statistical ensembles, and several types of simulation cells and boundary conditions are available. Intermolecular interactions comprise tabulated pairwise potentials for speed and uniformity and many-body interactions involve anisotropic polarizabilities. Intramolecular interactions include bond, angle, and crosslink potentials. A very large set of analyses of static and dynamic properties is provided. The capability of MOLSIM can be extended by user-providing routines controlling, for example, start conditions, intermolecular potentials, and analyses. An extensive set of case studies in the field of soft matter is presented covering colloids, polymers, and crosslinked networks.

The coupling constant polarizability and hyperpolarizabilty of 1J(NH) in N-methylacetamide, and its application for the multipole spin-spin coupling constant polarizability/reaction field approach to solvation

We present a benchmark study of a combined multipole spin-spin coupling constant (SSCC) polarizability/reaction field (MJP/RF) approach to the calculation of both specific and bulk solvation effects on SSCCs of solvated molecules. The MJP/RF scheme is defined by an expansion of the SSCCs of the solvated molecule in terms of coupling constant dipole and quadrupole polarizabilities and hyperpolarizabilities derived from single molecule ab initio calculations. The solvent electric field and electric field gradient are calculated based on data derived from molecular dynamics (MD) simulations thereby accounting for solute-solvent dynamical effects. The MJP/RF method is benchmarked against polarizable QM/MM calculations for the one-bond N-H coupling constant in N-methylacetamide. The best agreement between the MJP/RF and QM/MM approaches is found by truncating the electric field expansion in the MJP/RF approach at the linear electric field level. In addition, we investigate the sensitivity of the results due to the choice of one-electron basis set in the ab initio calculations of the coupling constant (hyper-)polarizabilities and find that they are affected by the basis set in a way similar to the coupling constants themselves.

Electric quadrupole polarizabilities of nuclear magnetic shielding in some small molecules

Computational procedures, based on (i) the Ramsey common origin approach and (ii) the continuous transformation of the origin of the quantum mechanical current density-diamagnetic zero (CTOCD-DZ), were applied at the Hartree-Fock level to determine electric quadrupole polarizabilities of nuclear magnetic shielding for molecules in the presence of a nonuniform electric field with a uniform gradient. The quadrupole polarizabilities depend on the origin of the coordinate system, but values of the magnetic field induced at a reference nucleus, determined via the CTOCD-DZ approach, are origin independent for any calculations relying on the algebraic approximation, irrespective of size and quality of the (gaugeless) basis set employed. On the other hand, theoretical estimates of the induced magnetic field obtained by single-origin methods are translationally invariant only in the limit of complete basis sets. Calculations of electric quadrupole polarizabilities of nuclear magnetic shielding are reported for H(2), HF, H(2)O, NH(3), and CH(4) molecules.

Refractive index of liquid water in different solvent models

We present a combined molecular dynamics/quantum chemical perturbation method for calculating the refractive index of liquid water at different temperatures. We compare results of this method with the refractive index obtained from other solvent models. The best agreement with the experimental refractive index of liquid water and its temperature dependence is obtained using correlated gas-phase polarizabilities in the classical Lorentz-Lorenz expression. Also, the iterative self-consistent reaction field approach in the semicontinuum implementation matches the experimental refractive index reasonably well.

Zero-point vibrational effects on proton shieldings: functional-group contributions from ab initio calculations

We investigate the effects of zero-point vibrational motion on the nuclear magnetic shielding constants of a large number of organic molecules. The vibrational corrections include anharmonic contributions from the potential energy surface and harmonic contributions from the curvature of the property surface. Particular attention is paid to vibrational corrections to hydrogen shielding constants where we show that vibrational corrections may be substantial, ranging from about +0.50 to -0.70 ppm, and thus demonstrating that ignoring these effects may give errors in the chemical shifts by more than 1 ppm in certain extreme cases. These effects can therefore not be neglected when comparing calculated results with experiment, not even for the chemical shifts. However, we also demonstrate that the vibrational corrections to the hydrogen shieldings are to a large extent transferable from one molecule to another. We have tabulated functional vibrational corrections to the hydrogen shieldings, based on results for more than 35 molecules. Unfortunately, no similar transferability has been observed for the vibrational corrections to shielding constants of other nuclei such as carbon, nitrogen, or oxygen.