Efficient multipole model and linear scaling of NDDO-based methods

Fast growth of computational costs with that of the system's size is a bottleneck for the applications of traditional methods of quantum chemistry to polyatomic molecular systems. This problem is addressed by the development of linear (or almost linear) scaling methods. In the semiempirical domain, it is typically achieved by a series of approximations to the self-consistent field (SCF) solution. By contrast, we propose a route to linear scalability by modifying the trial wave function itself. Our approach is based on variationally determined strictly local one-electron states and a geminal representation of chemical bonds and lone pairs. A serious obstacle previously faced on this route were the numerous transformations of the two-center repulsion integrals characteristic for the neglect of diatomic differential overlap (NDDO) methods. We pass it by replacing the fictitious charge configurations usual for the NDDO scheme by atomic multipoles interacting through semiempirical potentials. It ensures invariance of these integrals and improves the computational efficiency of the whole method. We discuss possible schemes for evaluating the integrals as well as their numerical values. The method proposed is implemented for the most popular modified neglect of diatomic overlap (MNDO), Austin model 1 (AM1), and PM3 parametrization schemes of the NDDO family. Our calculations involving well-justified cutoff procedures for molecular interactions unequivocally show that the proposed scheme provides almost linear scaling of computational costs with the system's size. The numerical results on molecular properties certify that our method is superior with respect to its SCF-based ancestors.     

Fast noniterative orbital localization for large molecules

We use Cholesky decomposition of the density matrix in atomic orbital basis to define a new set of occupied molecular orbital coefficients. Analysis of the resulting orbitals ("Cholesky molecular orbitals") demonstrates their localized character inherited from the sparsity of the density matrix. Comparison with the results of traditional iterative localization schemes shows minor differences with respect to a number of suitable measures of locality, particularly the scaling with system size of orbital pair domains used in local correlation methods. The Cholesky procedure for generating orthonormal localized orbitals is noniterative and may be made linear scaling. Although our present implementation scales cubically, the algorithm is significantly faster than any of the conventional localization schemes. In addition, since this approach does not require starting orbitals, it will be useful in local correlation treatments on top of diagonalization-free Hartree-Fock optimization algorithms.     

Bond-order potentials with split-charge equilibration: application to C-, H-, and O-containing systems

A method for extending charge transfer to bond-order potentials, known as the bond-order potential/split-charge equilibration (BOP/SQE) method [P. T. Mikulski, M. T. Knippenberg, and J. A. Harrison, J. Chem. Phys. 131, 241105 (2009)], is integrated into a new bond-order potential for interactions between oxygen, carbon, and hydrogen. This reactive potential utilizes the formalism of the adaptive intermolecular reactive empirical bond-order potential [S. J. Stuart, A. B. Tutein, and J. A. Harrison, J. Chem. Phys. 112, 6472 (2000)] with additional terms for oxygen and charge interactions. This implementation of the reactive potential is able to model chemical reactions where partial charges change in gas- and condensed-phase systems containing oxygen, carbon, and hydrogen. The BOP/SQE method prevents the unrestricted growth of charges, often observed in charge equilibration methods, without adding significant computational time, because it makes use of a quantity which is calculated as part of the underlying covalent portion of the potential, namely, the bond order. The implementation of this method with the qAIREBO potential is designed to provide a tool that can be used to model dynamics in a wide range of systems without significant computational cost. To demonstrate the usefulness and flexibility of this potential, heats of formation for isolated molecules, radial distribution functions of liquids, and energies of oxygenated diamond surfaces are calculated.     

Polarization and polarizability in extended one-dimensional organic materials

The modern theory of polarization in extended insulators is applied to one-dimensional models for conjugated polymers and charge transfer salts. Closed expressions for the dependence of the polarization on the site and bond energy alternations are presented for uncorrelated models, and results from exact real-space diagonalization are obtained for correlated models. Changes in polarization induced by lattice phonons or molecular vibrations are directly related to the intensity of infrared bands in the far and mid-IR, respectively. We model intensities by introducing linear electron-vibration coupling and show that coupling to delocalized electrons generates a combination band consisting of a lattice phonon and a molecular vibration. The displaced dipole operator is defined on a real-space basis allowing for the finite field calculation of linear polarizability in finite size systems with periodic boundary conditions. Size-consistency arguments are used to demonstrate that the resulting polarizability becomes exact in the thermodynamic limit, and numerical calculations demonstrate that this approach leads to reliable results that converge rapidly to the thermodynamic limit.     

Dielectric and Viscoelastic Study of Entanglement Dynamics: A Review of Recent Findings

This article gives a review of the results of recent dielectric and viscoelastic studies for entangled binary blends of linear cis-polyisoprenes to explain the current understanding of the equilibrium entanglement dynamics on the basis of the molecular picture of dynamic tube dilation (DTD). Comparison of dielectric and viscoelastic properties reveals that the full-DTD picture regarding the relaxed portions of the chains as a solvent fails for the high molecular weight component chain in the blends at intermediate times. This failure is related to insufficient constraint release (CR) equilibration of the entanglement segments of this chain. A partial-DTD picture properly considering this CR equilibration successfully describes the linear relaxation behavior of the blends. The dielectric and viscoelastic properties of PI under fast flow, being affected by the flow-activated CR/DTD mechanism, are also presented in order to demonstrate the usefulness of the comparison of these properties in both equilibrium and non-equilibrium states.     

Improvement of semiempirical response properties with charge-dependent response density

The present work outlines a new method for treatment of charge-dependent polarizability in semiempirical quantum models for use in combined quantum-mechanical/molecular mechanical simulations of biological reactions. The method addresses a major shortcoming in the performance of conventional semiempirical models for these simulations that is tied to the use of a localized minimal atomic-orbital basis set. The present approach has the advantages that it uses a density basis that retains a set of linear-response equations, does not increase the atomic-orbital basis, and avoids the problem of artificial charge transfer and scaling of the polarizability seen in related models that allow atomic charges to fluctuate. The model introduces four new atom-based parameters and has been tested with the modified neglect of differential overlap d-orbital Hamiltonian against 1132 molecules and ions and shown to decrease the dipole moment and polarizability errors by factors of 2 and 10, respectively, with respect to density-functional results. The method performs impressively for a variety of charge states (from 2+ to 2-), and offers a potentially powerful extension in the design of next generation semiempirical quantum models for accurate simulations of highly charged biological reactions.     

Molecule‐specific determination of atomic polarizabilities with the polarizable atomic multipole model

Recently, many polarizable force fields have been devised to describe induction effects between molecules. In popular polarizable models based on induced dipole moments, atomic polarizabilities are the essential parameters and should be derived carefully. Here, we present a parameterization scheme for atomic polarizabilities using a minimization target function containing both molecular and atomic information. The main idea is to adopt reference data only from quantum chemical calculations, to perform atomic polarizability parameterizations even when relevant experimental data are scarce as in the case of electronically excited molecules. Specifically, our scheme assigns the atomic polarizabilities of any given molecule in such a way that its molecular polarizability tensor is well reproduced. We show that our scheme successfully works for various molecules in mimicking dipole responses not only in ground states but also in valence excited states. The electrostatic potential around a molecule with an externally perturbing nearby charge also exhibits a near‐quantitative agreement with the reference data from quantum chemical calculations. The limitation of the model with isotropic atoms is also discussed to examine the scope of its applicability. © 2012 Wiley Periodicals, Inc.     

A Hirshfeld partitioning of polarizabilities of water clusters

A new Hirshfeld partitioning of cluster polarizability into intrinsic polarizabilities and charge delocalization contributions is presented. For water clusters, density-functional theory calculations demonstrate that the total polarizability of a water molecule in a cluster depends upon the number and type of hydrogen bonds the molecule makes with its neighbors. The intrinsic contribution to the molecular polarizability is transferable between water molecules displaying the same H-bond scheme in clusters of different sizes, and geometries, while the charge delocalization contribution also depends on the cluster size. These results could be used to improve the existing force fields.     

A generalization of the charge equilibration method for nonmetallic materials

Assigning effective atomic charges that properly reproduce the electrostatic fields of molecules is a crucial step in the construction of accurate interatomic potentials. We propose a new approach to calculate these charges, which as previous approaches are, is based on the idea of charge equilibration. However, we only allow charge to flow between covalently bonded neighbors by using the concept of so-called split charges. The semiempirical fit parameters in our approach do not only reflect atomic properties (electronegativity and atomic hardness) but also bond-dependent properties. The new method contains two popular but hitherto disjunct approaches as limiting cases. We apply our methodology to a set of molecules containing the elements silicon, carbon, oxygen, and hydrogen. Effective charges derived from electrostatic potential surfaces can be predicted more than twice as accurately as with previous works, at the expense of one additional fit parameter per bond type controlling the polarizability between two bonded atoms. Additional bond-type parameters can be introduced, but barely improve the results. An increase in accuracy of only 30% over existing techniques is achieved when predicting Mulliken charges. However, this could be improved with additional bond-type parameters.     

Linear scaling multireference singles and doubles configuration interaction

A linear scaling multireference singles and doubles configuration interaction (MRSDCI) method has been developed. By using localized bases to span the occupied and virtual subspace, local truncation schemes can be applied in tandem with integral screening to reduce the various bottlenecks in a MRSDCI calculation. Among these, the evaluation of electron repulsion integrals and their subsequent transformation, together with the diagonalization of the large CI Hamiltonian matrix, correspond to the most computationally intensive steps in a MRSDCI calculation. We show that linear scaling is possible within each step. The scaling of the method with system size is explored with a system of linear alkane chains and we proceed to demonstrate this method can produce smooth potential energy surfaces via calculating the dissociation of trans-6-dodecene (C(12)H(24)) along the central C[Double Bond]C bond.     

Dipole preserving and polarization consistent charges

A method for estimating dipole preserving and polarization consistent (DPPC) charges is described, which reproduces exactly the molecular dipole moment as well as the local, atomic hybridization dipoles determined from the corresponding wave function and can yield accurate molecular polarization. The method is based on a model described by Thole and van Duijnen and a new feature is introduced to treat molecular polarization. Thus, the DPPC method offers a convenient procedure to describe molecular polarization in applications using semiempirical models and ab initio molecular orbital theory with relatively small basis functions such as 6‐31+G(d,p) or without inclusion of electron correlation; these methods tend to underestimate molecular polarizability. The trends of the DPPC partial atomic charges are found to be in good accord with those of the CM2 model, a class IV charge analysis method that has been used in a variety of applications. The DPPC method is illustrated to mimic the correct molecular polarizability in a water dimer test case and in water‐halide ion complexes using the explicit polarization (X‐Pol) potential with the Austin model 1 Hamiltonian. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Charge-dependent model for many-body polarization, exchange, and dispersion interactions in hybrid quantum mechanical/molecular mechanical calculations

This work explores a new charge-dependent energy model consisting of van der Waals and polarization interactions between the quantum mechanical (QM) and molecular mechanical (MM) regions in a combined QMMM calculation. van der Waals interactions are commonly treated using empirical Lennard-Jones potentials, whose parameters are often chosen based on the QM atom type (e.g., based on hybridization or specific covalent bonding environment). This strategy for determination of QMMM nonbonding interactions becomes tedious to parametrize and lacks robust transferability. Problems occur in the study of chemical reactions where the "atom type" is a complex function of the reaction coordinate. This is particularly problematic for reactions, where atoms or localized functional groups undergo changes in charge state and hybridization. In the present work we propose a new model for nonelectrostatic nonbonded interactions in QMMM calculations that overcomes many of these problems. The model is based on a scaled overlap model for repulsive exchange and attractive dispersion interactions that is a function of atomic charge. The model is chemically significant since it properly correlates atomic size, softness, polarizability, and dispersion terms with minimal one-body parameters that are functions of the atomic charge. Tests of the model are examined for rare-gas interactions with neutral and charged atoms in order to demonstrate improved transferability. The present work provides a new framework for modeling QMMM interactions with improved accuracy and transferability.     

原子电荷计算方法的对比

原子电荷是对化学体系中电荷分布最简单、最直观的描述形式之一,在理论和实际应用中都有重要意义.本文介绍了12种重要的原子电荷计算方法的原理和特点,通过大量实例从不同角度比较了它们的优缺点.这些方法包括Mulliken、分子环境中的原子轨道(AOIM)、Hirshfeld、原子偶极矩校正的Hirshfeld布居(ADCH)、自然布居分析(NPA)、Merz-Kollmann(MK)、分子中的原子(AIM)、Merck分子力场94(MMFF94)、AM1-BCC、Gasteiger、电荷模型2(CM2)以及电荷均衡(QEq)方法.最后本文对如何在实际应用中选择合适的计算方法给出了建议.     

Comparing reduced partial charge models with polarizable simulations of ionic liquids

Molecular ionic liquids are typically characterized by strong electrostatic interactions resulting in a charge ordering and retardation of their translational and rotational behaviour. Unfortunately, this effect is often overestimated in classical molecular dynamics simulations. This can be circumvented in a twofold way: the easiest way is to reduce the partial charges of the ions to sub-integer values of ±0.7-0.9 e. The more realistic model is to include polarizable forces, e.g. Drude-oscillators, but it comes along with an increasing computational effort. On the other hand, charge-scaled models are claimed to take an average polarizability into account. But do both models have the same impact on structure and dynamics of molecular ionic liquids? In the present study several molecular dynamics simulations of 1-ethyl-3-methylimidazolium trifluoromethanesulfonate are performed with different levels of polarization as well as with varying charge scaling factors of 0.74 to 0.90. The analysis of the structural and dynamical results are performed in different levels: from the atomic point of view over the molecular level to collective properties determined by the complete sample.     

Quantitative structure-property relationships for longitudinal, transverse, and molecular static polarizabilities in polyynes

The present work reports for the first time quantitative structure-property relationships, derived at the benchmark CCSD(T)/cc-PVTZ level of theory that estimate the static longitudinal, transverse, and molecular polarizability in polyynes (C2nH2), as a function of their length (L). In the case of independent electron models, regardless of the form of the nuclei potential that the electrons experience, the polarizability increases strongly with system size, scaling as L(4). In contrast, the static longitudinal polarizability in polyynes have a considerably weaker length-dependence (L(1.64)). This is shown to predominantly arise from electron-electron repulsion rather than electron correlation by a systematic study of the polarizability length dependence in several simple quantum mechanical systems (e.g., particle-in-box, simple harmonic oscillator) and other molecular systems (e.g., H2, H2(+), polyynes). Decrease of the electron-electron repulsion term is suggested to be the key factor in enhancing nonlinear polarizability characteristics of linear oligomeric and polymeric materials.     

Fast tools for calculation of atomic charges well suited for drug design1

Two novel approaches to construct empirical schemes for partial atomic charge calculation were proposed. The charge schemes possess important benefits. First, they produce both topologically symmetrical and environment dependent charges. Second, they can be parameterised to reasonably reproduce ab initio molecular electrostatic potential (MEP), which guarantees their successful use in molecular modelling. To validate the approaches, the parameters of the proposed charge schemes were fitted to best reproduce MEP simultaneously on grids around a set of 227 diverse organic compounds. The residual errors in MEP reproduction due to calculated atomic charges were compared to those due to charges from known charge schemes.     

A reformulation of the coupled perturbed self-consistent field equations entirely within a local atomic orbital density matrix-based scheme

To exploit the exponential decay found in numerical studies for the density matrix and its derivative with respect to nuclear displacements, we reformulate the coupled perturbed self-consistent field (CPSCF) equations and a quadratically convergent SCF (QCSCF) method for Hartree-Fock and density functional theory within a local density matrix-based scheme. Our D-CPSCF (density matrix-based CPSCF) and D-QCSCF schemes open the way for exploiting sparsity and to achieve asymptotically linear scaling of computational complexity with molecular size (M), in case of D-CPSCF for all (M) derivative densities. Furthermore, these methods are even for small molecules strongly competitive to conventional algorithms.     

Predication of second-order optical nonlinearity of [(Bu2t Im)AuX] (X=halogen) using time-dependent density-functional theory combined with sum-over-states method

The dipole polarizability and second-order polarizability of recently synthesized (1,3-di-ter-butylimidazol-2-ylidine) gold complexes [(Bu2t Im)AuX] (X=halogen) were investigated by using time-dependent density-functional theory combined with sum-over-states method. We have discovered that these complexes possess remarkably larger molecular second-order polarizability compared with the organometallic and organic complexes. The value of the second-order polarizability increases in the order of F相似文献