Analytic calculations of vibrational hyperpolarizabilities in the atomic orbital basis

We present an analytic scheme for the calculation of pure vibrational contributions to linear and nonlinear optical properties such as the polarizability and the first and second hyperpolarizabilities. The formalism is fully expressed in terms of a perturbation- and time-dependent atomic orbital basis, using the elements of the density matrix in the atomic orbital basis as the basic variables. We calculate perturbed densities up to third order with respect to the electric field in accordance with the n + 1 rule, and the approach is therefore applicable for the calculation of pure vibrational contributions involving all vibrational coordinates in large molecular complexes. In the case of static electric fields, we therefore only need to calculate 19 response equations, independent of the size of the molecule. If we can determine the molecular energy and force field, the calculation of pure vibrational contributions to the nonlinear optical properties of the molecule is therefore a rather straightforward task. We illustrate the implementation by calculating pure vibrational contributions to the first and second hyperpolarizabilities of molecules containing up to 66 atoms using basis sets of good quality.     

A hybrid approach combining energy density analysis with the interaction energy decomposition method

We propose a new analysis technique for characterizing molecular interactions that combines an energy decomposition scheme, such as the Kitaura-Morokuma decomposition method, with energy density analysis, which partitions the total energy of the system into atomic contributions. The combined scheme, termed Interaction-EDA, enables us to estimate the local contribution of interaction energy components, such as electrostatic, exchange, polarization, and charge transfer. The evaluation of the local interaction energy is rather important in large systems. For a numerical assessment, the Interaction-EDA method is applied to the process of CO adsorption on Si(100) - (2 x 1) surface.     

Polarizability and second hyperpolarizability evaluation of long molecules by the density functional theory with long-range correction

Polarizabilities and second hyperpolarizabilities of polyacetylene and a hydrogen chain are evaluated by density functional theory (DFT) using a hybrid generalized gradient approximation functional with correct long-range electron-electron interactions. The well known catastrophic overestimate of the hyperpolarizabilities for molecular systems of enhanced length is corrected by the two-electron repulsion operator decomposition technique, integrating the distance-dependent nonlocal exchange effects for long-range interaction, while neither the asymptotically corrected exchange functional for long-range interaction nor ordinary hybrid methods seem to be capable of overcoming the serious drawback of the DFT in polarizability/hyperpolarizability evaluation.     

A perspective on the localizability of Hartree–Fock orbitals

A common perception about molecular systems with a nonlocal electronic structure (as manifested by a nonlocal Hartree–Fock (HF) density matrix), such as conjugated π-systems, is that they can only be described in terms of nonlocal molecular orbitals. This view is mostly founded on chemical intuition, and further, this view is strengthened by traditional approaches for obtaining local occupied and virtual orbital spaces, such as the occupied Pipek–Mezey orbitals, and projected atomic orbitals. In this article, we discuss the limitations for localizability of HF orbitals in terms of restrictions posed by the delocalized character of the underlying density matrix for the molecular system and by the orthogonality constraint on the molecular orbitals. We show that the locality of the orbitals, in terms of nonvanishing charge distributions of orbitals centered far apart, is much more strongly affected by the orthogonality constraint than by the physical requirement that the occupied orbitals must represent the electron density. Thus, the freedom of carrying out unitary transformations among the orbitals provides the flexibility to obtain highly local occupied and virtual molecular orbitals, even for molecular systems with a nonlocal density matrix, provided that a proper localization function is used. As an additional consideration, we clear up the common misconception that projected atomic orbitals in general are more local than localized orthogonal virtual orbitals.     

Solvent and intermolecular effects on first hyperpolarizabilities of organometallic tungsten-carbonyl complexes, a TDDFT study

The first hyperpolarizability of two tungsten-carbonyl complexes, tungsten pentacarbonyl pyridine and tungsten pentacarbonyl trans-1,2-bis(4-pyridyl)-ethylene, has been studied by the high-level TDDFT method. The consideration of the solvent effect and intermolecular pi-pi weak interaction in the calculations quantitatively improve the final result of both the electronic excitations and the first hyperpolarizabilities. By using the orbital decomposition scheme (J. Phys. Chem. A 2006, 110, 1014-1021), the NLO mechanisms of these two complexes have been ascribed to the dominant contribution from the metal-to-ligand charge transfer, with HOMO --> LUMO character, and the indispensable contribution from the intraligand charge transfer as well. A supplementary formula has been proposed to implement the orbital-pair transition analysis. This study reports the significant influences of solvation and intermolecular interactions on the first hyperpolarizabilities of organometallic NLO chromophores.     

Theory of intermolecular interactions: The long range terms in the dipole–dipole,monopoles–dipole,and monopoles–bond polarizabilities approximations

The problem of evaluating the long range terms (electrostatic, polarization, dispersion) of the interaction energy between molecules at intermediate distances (i.e. distances of the order of magnitude of the molecular dimensions) is considered. Instead of being approximated by its dipole part, the exact interaction Hamiltonian is treated as proposed by Longuet-Higgins [11], i.e. the matrix elements are interpreted as electrostatic interactions between state and transition charge distributions. These charge distributions are approximated in a systematic way by sets of point charges (localized on the atoms) or sets of dipoles (localized on the bonds). The various contributions to the energy may then be expressed in terms of atomic net charges and bond polarizabilities. More refined approximations of the charge distributions could be used and correspondingly improved formulae could be derived: as an example, a formula for the σ-π dispersion energy is derived, where the σ charge distributions are approximated by bond transition dipoles (leading to σ bond polarizabilities in the final formula) while the π charge distributions are approximated by atomic charges.     

The kinetic energy of molecular charge distributions and molecular stability

This paper examines the relationship between the topographical features of a molecular charge distribution and the kinetic energy of the system. Specifically, the spatial contributions to the kinetic energy are related to the Laplacian of the total charge density and to the gradients of the natural-orbital densities. It is concluded that a necessary requirement for molecular stability is the existence of a net negative curvature for the molecular charge distribution in the internuclear region. It is shown that the charge density accumulated in the internuclear region of a stable molecule is distributed in such a way as to keep the accompanying increase in the kinetic energy to a minimum. A comparison of the contributions to the kinetic energy from the atomic and molecular charge distributions indicates that in the formation of a stable molecule the contribution from the molecular charge density in the binding region is decreased relative to that of the atoms.     

A study of the partitioning of the first-order reduced density matrix according to the theory of atoms in molecules

This work describes a simple spatial decomposition of the first-order reduced density matrix corresponding to an N-electron system into first-order density matrices, each of them associated to an atomic domain defined in the theory of atoms in molecules. A study of the representability of the density matrices arisen from this decomposition is reported and analyzed. An appropriate treatment of the eigenvectors of the matrices defined over atomic domains or over unions of these domains allows one to describe satisfactorily molecular properties and chemical bondings within a determined molecule and among its fragments. Numerical determinations, performed in selected molecules, confirm the reliability of our proposal.     

Cumulative atomic multipole representation of the molecular charge distribution and its basis set dependence

A simple procedure to decompose the theoretical molecular charge distribution into cumulative atomic multipoles supplementing any population analysis scheme has been described and tested for a number of molecules in extended basis sets. This approach may be applied to describe local charge distributions in neutral as well as charged systems and also leads to a simplified point-charge model conserving the local anisotropy of the atomic charge distribution in molecules. Such an approach may be useful in estimating intermolecular interactions, representing the molecular environment in solvent effect or enzyme catalytic activity studies, evaluation of molecular electrostatic potentials or tracing the quality of basis set functions.     

Local and nonlocal density functional calculations of the molecular structure of isomeric thiadiazoles

The chemistry of thiadiazoles and their derivatives is of considerable interest in chemistry owing to their pharmacological and potential industrial applications. In this context, a detailed study of isomeric thiadiazole molecules has been done using local (SVWN; Slater, and Vosko, Wilk and Nusair) and nonlocal (BLYP; Becke, and Lee, Yang and Parr) density functionals and optimizing the molecular geometries by means of the gradient technique. A charge sensitivity analysis of the studied molecule has been performed by resorting to density functional theory, obtaining several sensitivity coefficients such as the molecular energy, net atomic charges, global and local hardness, global and local softness and Fukui functions. With these results and the analysis of the dipole moments, the molecular electrostatic potentials and the total electron density maps, several conclusions have been inferred about the preferred sites of chemical reaction of the studied compounds. The condensed Fukui functions are shown to be one of the best criteria for predicting chemical reactivity.     

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.     

Net charge changes in the calculation of relative ligand‐binding free energies via classical atomistic molecular dynamics simulation

The calculation of binding free energies of charged species to a target molecule is a frequently encountered problem in molecular dynamics studies of (bio‐)chemical thermodynamics. Many important endogenous receptor‐binding molecules, enzyme substrates, or drug molecules have a nonzero net charge. Absolute binding free energies, as well as binding free energies relative to another molecule with a different net charge will be affected by artifacts due to the used effective electrostatic interaction function and associated parameters (e.g., size of the computational box). In the present study, charging contributions to binding free energies of small oligoatomic ions to a series of model host cavities functionalized with different chemical groups are calculated with classical atomistic molecular dynamics simulation. Electrostatic interactions are treated using a lattice‐summation scheme or a cutoff‐truncation scheme with Barker–Watts reaction‐field correction, and the simulations are conducted in boxes of different edge lengths. It is illustrated that the charging free energies of the guest molecules in water and in the host strongly depend on the applied methodology and that neglect of correction terms for the artifacts introduced by the finite size of the simulated system and the use of an effective electrostatic interaction function considerably impairs the thermodynamic interpretation of guest‐host interactions. Application of correction terms for the various artifacts yields consistent results for the charging contribution to binding free energies and is thus a prerequisite for the valid interpretation or prediction of experimental data via molecular dynamics simulation. Analysis and correction of electrostatic artifacts according to the scheme proposed in the present study should therefore be considered an integral part of careful free‐energy calculation studies if changes in the net charge are involved. © 2013 The Authors Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Atoms-in-molecules in momentum space: a Hirshfeld partitioning of electron momentum densities

A direct application of the Hirshfeld atomic partitioning (HAP) scheme is implemented for molecular electron momentum densities (EMDs). The momentum density contributions of individual atoms in diverse molecular systems are analyzed along with their topographical features and the kinetic energies of the atomic partitions. The proposed p-space HAP-based charge scheme does seem to possess the desirable attributes expected of any atoms in molecules partitioning. In addition to this, the main strength of the p-space HAP is the exact knowledge of the kinetic energy functional and the inherent ease in computing the kinetic energy. The charges derived from HAP in momentum space are found to match chemical intuition and the generally known chemical characteristics such as electronegativity, etc.     

A density functional theory-based chemical potential equalisation approach to molecular polarizability

The electron density changes in molecular systems in the presence of external electric fields are modeled for simplicity in terms of the induced charges and dipole moments at the individual atomic sites. A chemical potential equalisation scheme is proposed for the calculation of these quantities and hence the dipole polarizability within the framework of density functional theory based linear response theory. The resulting polarizability is expressed in terms of the contributions from individual atoms in the molecule. A few illustrative numerical calculations are shown to predict the molecular polarizabilities in good agreement with available results. The usefulness of the approach to the calculation of intermolecular interaction needed for computer simulation is highlighted.     

Enhancement of the third-order nonlinear optical properties in open-shell singlet transition-metal dinuclear systems: effects of the group, of the period, and of the charge of the metal atom

Metal-metal multiply bonded complexes in their singlet state have been predicted to form a novel class of "σ-dominant" third-order nonlinear optical compounds based on the results of dichromium(II) and dimolybdenum(II) systems (H. Fukui et al. J. Phys. Chem. Lett.2011, 2, 2063) whose second hyperpolarizabilities (γ) are enhanced by the contribution of the dσ electrons with an intermediate diradical character. In this study, using the spin-unrestricted coupled-cluster method with singles and doubles as well as with perturbative triples, we investigate the dependences of γ on the group and on the period of the transition metals as well as on their atomic charges in several open-shell singlet dimetallic systems. A significant enhancement of γ is observed in those dimetallic systems composed of (i) transition metals with a small group number, (ii) transition metals with a large periodic number, and (iii) transition metals with a small positive charge. From the decomposition of the γ values into the contributions of dσ, dπ, and dδ electrons, the γ enhancements are shown to originate from the dσ contribution, because it corresponds to the intermediate diradical character region. Furthermore, the amplitude of dσ contribution turns out to be related to the size of the d(z(2)) atomic orbital of the transition metal, which accounts for the dependence of γ on the group, on the period, and on the charge of the metal atoms. These dependences provide a guideline for an effective molecular design of highly efficient third-order nonlinear optical (NLO) systems based on the metal-metal bonded systems.     

A numerically stable restrained electrostatic potential charge fitting method

Inspired by the idea of charge decomposition in calculation of the dipole preserving and polarization consistent charges (Zhang et al., J. Comput. Chem. 2011, 32, 2127), we have proposed a numerically stable restrained electrostatic potential (ESP)‐based charge fitting method for protein. The atomic charge is composed of two parts. The dominant part is fixed to a predefined value (e.g., AMBER charge), and the residual part is to be determined by restrained fitting to residual ESP on grid points around the molecule. Nonuniform weighting factors as a function of the dominant charge are assigned to the atoms. Because the residual part is several folds to several orders smaller than the dominant part, the impact of ill‐conditioning is alleviated. This charge fitting method can be used in quantum mechanical/molecular mechanical (QM/MM) simulations and similar studies, where QM calculated electronic properties are frequently mapped to partial atomic charges. © 2012 Wiley Periodicals, Inc.     

First-principles calculation of local atomic polarizabilities

Common methods of determining atomic polarizabilities suffer from the inclusion of nonlocal effects such as charge polarization. A new method is described for determining fully ab initio atomic polarizabilities based on calculating the response of atomic multipoles to the local electrostatic potential. The localized atomic polarizabilities are then used to calculate induction energies that are compared to ab initio induction energies to test their usefulness in practical applications. These polarizabilities are shown to be an improvement over the corresponding molecular polarizabilities, in terms of both absolute accuracy and the convergence of the multipolar induction series. The transferability of localized polarizabilities for the alkane series is also discussed.     

In-Situ Electronegativity and the Bridging of Chemical Bonding Concepts

One challenge in chemistry is the plethora of often disparate models for rationalizing the electronic structure of molecules. Chemical concepts abound, but their connections are often frail. This work describes a quantum-mechanical framework that enables a combination of ideas from three approaches common for the analysis of chemical bonds: energy decomposition analysis (EDA), quantum chemical topology, and molecular orbital (MO) theory. The glue to our theory is the electron energy density, interpretable as one part electrons and one part electronegativity. We present a three-dimensional analysis of the electron energy density and use it to redefine what constitutes an atom in a molecule. Definitions of atomic partial charge and electronegativity follow in a way that connects these concepts to the total energy of a molecule. The formation of polar bonds is predicted to cause inversion of electronegativity, and a new perspective of bonding in diborane and guanine−cytosine base-pairing is presented. The electronegativity of atoms inside molecules is shown to be predictive of pK_a.     

Modified hybridization displacement charge scheme for π‐electron systems: Study of molecular electrostatic potential maps

A modified form of the hybridization displacement charge scheme (HDC‐PI) suitable for application to molecules where π‐electrons play the dominant role in controlling molecular electrostatic potential (MEP) features is presented. This modified scheme does not only preserve the atomic contributions to molecular dipole moment but also it represents the two opposing components of atomic hybridization dipole moment arising due to π‐electrons, one due to the π‐electrons located above the molecular plane and the other due to those located below the plane, so as to reproduce MEP features appropriately. The method has been applied within the framework of the AMl and MNDO methods to several different types of molecules containing π‐electrons demonstrating its very satisfactory performance. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 191–200, 1999