Transferability of topologically partitioned polarizabilities: the case of n-alkanes

A fast approximate extension of the interacting quantum atoms energy decomposition to excited states

An approach is developed for the fast calculation of the interacting quantum atoms energy decomposition (IQA) from the information contained in the first order reduced density matrix only. The proposed methodology utilizes an approximate exchange-correlation density from Density Matrix Functional Theory without the need to evaluate the correlation-exchange contribution directly. Instead, weight factors are estimated to decompose the exact V_xc into atomic and pairwise contributions. In this way, the sum of the IQA contributions recovers the energy obtained from the electronic structure calculation. This method can, hence, be applied to obtain atomic contributions in excited states on the same footing as in their ground states using any method that delivers the reduced first-order density matrix. In this way, one can locate chromophores from first principles quantum chemical calculations. Test calculations on the ground and excited states of a set of small molecules indicate that the scaled atomic contributions reproduce vertical electronic transition energies calculated exactly. This approach may be useful to extend the applicability of the IQA approach in the study of large photochemical systems especially when the calculations of the second order reduced density matrices is prohibitive or not possible.

Origin and control of superlinear polarizability scaling in chemical potential equalization methods

Many common chemical potential equalization (muEq) methods are known to suffer from a superlinear scaling of the polarizability with increasing molecular size that interferes with model transferability and prevents the straightforward application of these methods to large, biochemically relevant molecules. In the present work, we systematically investigate the origins of this scaling and the mechanisms whereby some existing methods successfully temper the scaling. We demonstrate several types of topological charge constraints distinct from the usual single molecular charge constraint that can successfully achieve linear polarizability scaling in atomic charge based equilibration models. We find the use of recently employed charge conservation constraints tied to small molecular units to be an effective and practical approach for modulating the polarizability scaling in atomic muEq schemes. We also analyze the scaling behavior of several muEq schemes in the bond representation and derive closed-form expressions for the polarizability scaling in a linear atomic chain model; for a single molecular charge constraint these expressions demonstrate a cubic dependence of the polarizability on molecular size compared with linear scaling obtainable in the case of the atom-atom charge transfer (AACT) and split-charge equilibration (SQE) schemes. Application of our results to the trans N-alkane series reveals that in certain situations, the AACT and SQE schemes can become unstable due to an indefinite Hessian matrix. Consequently, we discuss sufficient criteria for ensuring stability within these schemes.

Distributed polarizability analysis for para-nitroaniline and meta-nitroaniline: Functional group and charge-transfer contributions

Topological partitioning of electronic properties is used to investigate the polarizability of para-nitroaniline and meta-nitroaniline. The distributed polarizabilities for atoms are combined into total local or generalized distributed contributions for the amino, ring, and nitro functional groups; generalized distributed group contributions have not been calculated before. The local group contributions are transferable between the two molecules only when charge transfer is suppressed, but the generalized distributed contributions prove surprisingly similar in the two molecules, apparently because they treat charge-transfer contributions explicitly.

Microscopic calculation of the energetics of charged states in amorphous polyethylene

Polarization energies are calculated for a single excess charge on a polyethylene chain in amorphous polyethylene using (i) local segment and nonlocal distributed molecular polarizabilities, (ii) material structures simulated by both general-purpose and specialist Monte Carlo software, and (iii) uniform and Gaussian distributions of charge with different extents of charge delocalization. Local and distributed response lead to results that are essentially the same except that they correspond to different mean polarizabilities. With increasing delocalization of the charge along the chain, the polarization energies shift to higher values and the width of their distribution decreases, the differences being more pronounced for the uniform distribution. The polarization energies for charges delocalized over 10-20 methylene units form a distribution some 14 eV wide centered around 1 eV, narrowing significantly for more homogeneous polymer melts. The calculations are relevant to trapping of charge in polyethylene. They also yield the microscopic variation in the potential along the polymer chain caused by the polarization energy difference, and so may provide useful inputs to theories of electronic conduction in polymer materials.

OPEP: a tool for the optimal partitioning of electric properties

OPEP is a suite of FORTRAN programs targeted at the optimal partitioning of molecular electric properties. It includes an interactive module for the construction of Cartesian grids of points, on which either the molecular electrostatic potential or the induction energy is mapped. The generation of distributed multipoles and polarizabilities is achieved using either the formalism of the normal equations of the least-squares problem, which restates the fitting procedure in terms of simple matrix operations, or a statistical approach, which provides a pictorial description of the distributed models of multipoles and polarizabilities, thereby allowing the pinpointing of pathological cases. Molecular symmetry is accounted for by means of local atomic frames, which are generated in an automated fashion. A Tcl/Tk graphical user interface wraps the suite of programs, thereby making OPEP a user-friendly package for building models of distributed multipoles and polarizabilities. OPEP is an open-source suite of programs distributed free of charge under the GNU general public license (GPL) at http://www.lctn.uhp-nancy.fr/Opep.

Application of the quantum theory of atoms in molecules to selected physico-chemical and biophysical problems: focus on correlation with experiment

This article reviews how the quantum theory of atoms in molecules (QTAIM) can be used to predict experimental physico-chemical properties of molecules of biologic interest: the amino acids, the polycyclic aromatic hydrocarbons (PAH), and the opiates, for example, morphine and PEO. The predicted experimental properties are as diverse as the partial molar volumes, the free energies of hydration, the second code-letter in the genetic code, the resonance energies, and the proton spin-spin coupling constants. Recent examples of the utilization of QTAIM to construct excellent statistical models (with squared correlation coefficients (r(2)) > 0.9) correlating properties of the electron density and of the pair density to experiment are reviewed. Some new results on the solvent effects on electron delocalization are also presented.

An atoms-in-molecules study of the genetically-encoded amino acids: I. Effects of conformation and of tautomerization on geometric, atomic, and bond properties

The theory of Atoms-In-Molecules (AIM) is a partitioning of the real space of a molecule into disjoint atomic constituents as determined by the topology of the electron density, rho(r). This theory identifies an atom in a molecule with a quantum mechanical open system and, consequently, all of the atom's properties are unambiguously defined. AIM recovers the basic empirical cornerstone of chemistry: that atoms and functional groups possess characteristic and additive properties that in many cases exhibit a remarkable transferability between different molecules. As a result, the theory enables the theoretical synthesis of a large molecule and the prediction of its properties by joining fragments that are predetermined as open systems. The present article is the first of a series (in preparation) that explore this possibility for polypeptides by determining the transferability of the building blocks: the amino acid residues. Transferability of group properties requires transferability of the electron density rho(r), which in turn requires the transferability of the geometric parameters. This article demonstrates that these parameters are conformation-insensitive for a representative amino acid, leucine, and that the atomic and bond properties exhibit a corresponding transferability. The effects of hydrogen bonding are determined and a set of geometrical conditions for the occurrence of such bonding is identified. The effects of transforming neutral leucine into its zwitter-ionic form on its atomic and bond properties are shown to be localized primarily to the sites of ionization.