Congested molecules. Where is the steric repulsion? An analysis of the electron density by the method of interacting quantum atoms

The computed electron density of several congested saturated hydrocarbons and halogenated derivatives has been analyzed by the method of interacting quantum atoms (IQA). For all the molecules studied, the calculations show the existence of a bond path between the congested atoms and which, according to the Quantum Theory of Atoms in Molecules, indicates that there is a stabilizing interaction between these atoms. The bond path is found to exist up to interatomic distances well‐beyond the sum of the van der Waals radii. The IQA results indicate that steric hindrance is not a repulsive force between the congested atoms but that is the result of an increase in the intra‐atomic or self‐energy of the congested atoms. This increase in self‐energy is caused by the deformation of the atomic basin of the congested atoms. © 2013 Wiley Periodicals, Inc.     

Understanding the ring current effects on magnetic shielding of hydrogen and carbon nuclei in naphthalene and anthracene

The local response to an external magnetic field normal to the molecular plane of naphthalene and anthracene was investigated via current density and magnetic shielding density maps. The Biot-Savart law shows that the deshielding caused by pi-ring currents in naphthalene is stronger for alpha- than for beta-protons due to geometrical factors. The shielding tensor of the carbon nuclei in both molecules is strongly anisotropic and its out-of-plane component determines the up-field chemical shift of (13)C in nuclear magnetic resonance spectra. The pi-ring currents flowing beyond the C-skeleton in front of a probe carbon nucleus, and on remote parts of the molecular perimeter, yield positive contributions to the out-of-plane component of carbon shielding as big as approximately 10-15% of the total values. Near Hartree-Fock estimates of magnetizability and magnetic shielding at the nuclei fully consistent with the current model are reported.     

Hydrogen–hydrogen interaction in planar biphenyl: A theoretical study based on the interacting quantum atoms and Hirshfeld atomic energy partitioning methods

The nature of H‐H interaction between ortho‐hydrogen atoms in planar biphenyl is investigated by two different atomic energy partitioning methods, namely fractional occupation iterative Hirshfeld (FOHI) and interacting quantum atoms (IQA), and compared with the traditional virial‐based approach of quantum theory of atoms in molecules (QTAIM). In agreement with Bader's hypothesis of H? H bonding, partitioning the atomic energy into intra‐atomic and interatomic terms reveals that there is a net attractive interaction between the ortho‐hydrogens in the planar biphenyl. This falsifies the classical view of steric repulsion between the hydrogens. In addition, in contrast to the traditional QTAIM energy analysis, both FOHI and IQA show that the total atomic energy of the ortho‐hydrogens remains almost constant when they participate in the H‐H interaction. Although, the interatomic part of atomic energy of the hydrogens plays a stabilizing role during the formation of the H? H bond, it is almost compensated by the destabilizing effects of the intra‐atomic parts and consequently, the total energy of the hydrogens remains constant. The trends in the changes of intra‐atomic and interatomic energy terms of ortho‐hydrogens during H? H bond formation are very similar to those observed for the H₂ molecule. © 2014 Wiley Periodicals, Inc.     

Linear scaling approaches to quantum macromolecular similarity: evaluating the similarity function

The evaluation of the electron density based similarity function scales quadratically with respect to the size of the molecules for simplified, atomic shell densities. Due to the exponential decay of the function's atom-atom terms most interatomic contributions are numerically negligible on large systems. An improved algorithm for the evaluation of the Quantum Molecular Similarity function is presented. This procedure identifies all non-negligible terms without computing unnecessary interatomic squared distances, thus effectively turning to linear scaling the similarity evaluation. Presented also is a minimalist dynamic electron density model. Approximate, single shell densities together with the proposed algorithm facilitate fast electron density based alignments on macromolecules.     

Magnetic ordering and metal‐atom site preference in tetragonal CrMnAs: Electronic correlation effects

The electronic and magnetic structures of tetragonal, Cu₂Sb‐type CrMnAs were examined using density functional theory. To obtain reasonable agreement with reported atomic and low‐temperature magnetic ordering in this compound, the intra‐atomic electron–electron correlation in term of Hubbard U on Mn atoms are necessary. Using GGA + U, calculations identify four low‐energy antiferromagnetically ordered structures, all of which adopt a magnetic unit cell that contains the same direct Cr Cr and Cr Mn magnetic interaction, as well as the same indirect Mn⋅⋅⋅Mn magnetic interaction across the Cr planes. One of these low‐energy configurations corresponds to the reported case. Effective exchange parameters for metal–metal contacts obtained from SPRKKR calculations indicate both direct and indirect exchange couplings play important roles in tetragonal CrMnAs. © 2018 Wiley Periodicals, Inc.     

Ring current effects on nuclear magnetic shielding of carbon in the benzene molecule

The differential Biot-Savart law of classical electrodynamics was applied to develop a ring current model for the magnetic shielding of the carbon nucleus in benzene. It is shown that the local effect of the pi currents, induced by a magnetic field normal to the molecular plane, on the sigmaC out-of-plane shielding tensor component vanishes. However, approximately 10% of sigmaC is due to the shielding contributions from pi current density in the region of the other carbon atoms. Magnetic shielding density maps obtained via quantum mechanical procedures confirm the predictions of the classical model. Copyright (c) 2004 John Wiley & Sons, Ltd.     

Interatomic magnetizability: a QTAIM-based approach toward deciphering magnetic aromaticity

Interatomic magnetizability provides insight into the extent of electronic current density between two adjacent atomic basins. By studying a number of well-known aromatic, nonaromatic, and antiaromatic molecules, it is demonstrated that interatomic magnetizability (bond magnetizability) not only is able to verify the exact nature of aromaticity/antiaromaticity among different molecules, but also can distinguish the correct aromaticity order among sets of aromatic/antiaromatic molecules. The interatomic magnetizability is a direct measure of the current flux between two adjacent atomic basins and is the first QTAIM-derived index that evaluates aromaticity based on a response property, that is, magnetizability. Bond magnetizability is easy to compute, straightforward to interpret, and can be employed to evaluate the pure π- or σ-orbital contributions to magnetic aromaticity.     

Nuclear magnetic shielding constants of liquid water: insights from hybrid quantum mechanics/molecular mechanics models

We present a gauge-origin independent method for the calculation of nuclear magnetic shielding tensors of molecules in a structured and polarizable environment. The method is based on a combination of density functional theory (DFT) or Hartree-Fock wave functions with molecular mechanics. The method is unique in the sense that it includes three important properties that need to be fulfilled in accurate calculations of nuclear magnetic shielding constants: (i) the model includes electron correlation effects, (ii) the model uses gauge-including atomic orbitals to give gauge-origin independent results, and (iii) the effect of the environment is treated self-consistently using a discrete reaction-field methodology. The authors present sample calculations of the isotropic nuclear magnetic shielding constants of liquid water based on a large number of solute-solvent configurations derived from molecular dynamics simulations employing potentials which treat solvent polarization either explicitly or implicitly. For both the (17)O and (1)H isotropic shielding constants the best predicted results compare fairly well with the experimental data, i.e., they reproduce the experimental solvent shifts to within 4 ppm for the (17)O shielding and 1 ppm for the (1)H shielding.     

NMR shielding constants of CuX,AgX, and AuX (X = F,Cl, Br,and I) investigated by density functional theory based on the Douglas–Kroll–Hess Hamiltonian

Two‐component relativistic density functional theory (DFT) with the second‐order Douglas–Kroll–Hess (DKH2) one‐electron Hamiltonian was applied to the calculation of nuclear magnetic resonance (NMR) shielding constant. Large basis set dependence was observed in the shielding constant of Xe atom. The DKH2‐DFT‐calculated shielding constants of I and Xe in HI, I₂, CuI, AgI, and XeF₂ agree well with those obtained by the four‐component relativistic theory and experiments. The Au NMR shielding constant in AuF is extremely more positive than in AuCl, AuBr, and AuI, as reported recently. This extremely positive shielding constant arises from the much larger Fermi contact (FC) term of AuF than in others. Interestingly, the absolute values of the paramagnetic and the FC terms are considerably larger in CuF and AuF than in others. The large paramagnetic term of AuF arises from the large d‐components in the Au d_π –F p_π and Au sd_σ–F p_σ molecular orbitals (MOs). The large FC term in AuF arises from the small energy difference between the Au sd_σ + F p_σ and Au sd_σ–F p_σ MOs. The second‐order magnetically relativistic effect, which is the effect of DKH2 magnetic operator, is important even in CuF. This effect considerably improves the overestimation of the spin‐orbit effect calculated by the Breit–Pauli magnetic operator. © 2013 Wiley Periodicals, Inc.     

A Comparison of the Interacting Quantum Atoms (IQA) Analysis of the Two-Particle Density-Matrices of MP4SDQ and CCSD

Within the quantum topological energy partitioning method called Interacting Quantum Atoms (IQA) we transition from Møller-Plesset (MP4SDQ) to CCSD in calculating intra- and interatomic electron correlation energies for a set of hydrides, diatomics, a few simple molecules and non-covalently bonded complexes, using the uncontracted basis set 6-31++G(2d,2p). CCSD-IQA allows a more rigorous analysis of atomic electron correlation than that offered by Møller-Plesset, which returns IQA contributions that are identical to Hartree–Fock counterparts except for two-electron terms. The CCSD-IQA analysis returns bond and other interatomic correlation energies that are typically much larger in magnitude than the MP4SDQ values. Crisp patterns of energy transferability are detected in water clusters, both for intra-atomic and interatomic correlation energies. CCSD determines that the intra-atomic correlation energy of an oxygen drops by 15 kJ · mol^–1 for donating a hydrogen and by 25 kJ · mol^–1 for accepting a hydrogen.     

Is the Conventional Interpretation of the Anisotropic Effects of CC Double Bonds and Aromatic Rings in NMR Spectra in Terms of the π‐Electron Shielding/Deshielding Contributions Correct?

Based on the nucleus‐independent chemical shift (NICS) concept, isotropic magnetic shielding values have been computed along the three Cartesian axes for ethene, cyclobutadiene, benzene, naphthalene, and benzocyclobutadiene, starting from the molecular/ring center up to 10 Å away. These through‐space NMR spectroscopic shielding (TSNMRS) values, which reflect the anisotropic effects, have been broken down into contributions from localized‐ and canonical molecular orbitals (LMOs and CMOs); these contributions revealed that the proton NMR spectroscopic chemical shifts of nuclei that are spatially close to the C?C double bond or the aromatic ring should not be explained in terms of the conventionally accepted π‐electron shielding/deshielding effects. In fact, these effects followed the predictions only for the antiaromatic cyclobutadiene ring.     

Counterintuitive deshielding on the 13C NMR chemical shift for the trifluoromethyl anion

The trifluoromethyl anion (CF₃⁻) displays ¹³C NMR chemical shift (175.0 ppm) surprisingly larger than neutral (CHF₃, 122.2 ppm) and cation (CF₃⁺, 150.7 ppm) compounds. This unexpected deshielding effect for a carbanion is investigated by density functional theory calculations and decomposition analyses of the ¹³C shielding tensor into localized molecular orbital contributions. The present work determines the shielding mechanisms involved in the observed behaviour of the fluorinated anion species, shedding light on the experimental NMR data and demystify the classical correlation between electron density and NMR chemical shift. The presence of fluorine atoms induces the carbon lone pair to create a paramagnetic shielding on the carbon nucleus.     

Plasmon resonances and plasmon‐induced charge transport in linear atomic chains

In linear hydrogen atomic chains, plasmon resonances and plasmon‐induced charge transport are studied by time‐dependent density functional theory. For the large linear chain, it is a general phenomenon that, in the longitudinal excitation, there are high‐energy resonances and a large low‐energy resonance. The energy of the large low‐energy resonance conforms to the results calculated by the classical Drude model. In order to explain the formation mechanism of the high‐energy resonances, we present a simple harmonic oscillator model. This model may reasonably account for the relationship between low‐energy and high‐energy resonances, and has a certain degree of universality. As the interatomic distance decreases, the current shows a gradual transition from insulator to metal. The current enhancement mainly depends on the local field enhancement associated with plasmon excitation, and the enhanced electron delocalization effect as a result of the decrease of the interatomic distance. © 2013 Wiley Periodicals, Inc.     

Assessment of sigma-diatropicity of the cyclopropane molecule

Spatial models of the current density field induced in the cyclopropane molecule by stationary, homogeneous magnetic fields, parallel to either the C3 or the C2 symmetry axis, have been constructed. A compact, abridged representation of the models is given via stagnation graphs that convey essential information. Maps of streamlines and moduli are also reported to complete current models that have proven useful to rationalize magnetic tensor properties, that is, magnetizability, 1H and 13C nuclear shieldings, and magnetic shielding along the C3 symmetry axis. Plots of Biot-Savart magnetic shielding density combined with current density visualization yield an accurate, detailed account of the shielding mechanisms. The magnetropicity of the system described by the current density model is fully consistent with the magnitude of magnetic tensors calculated at near Hartree-Fock level. In a field perpendicular to the molecular plane, cyclopropane sustains a diatropic sigma-ring current with the following peculiar features: (i) it follows the molecular periphery rather than the CC framework; (ii) it bifurcates in the proximity of the methylene moieties flowing along the CH bonds, both above and below the sigma(h) plane; (iii) it has an effect on the values of response properties, although it is not as large as expected from naive arguments (e.g., the center-of-mass value of the magnetic shielding constant is dominated by in-plane components rather than the out-of-plane component, which is in contrast to pi-aromatic systems such as benzene); (iv) it has a negligible effect on the strong anisotropy of carbon magnetic shielding, which is shown to arise from local currents. No evidence for strong diatropism, and therefore sigma-aromaticity of the cyclopropane molecule, was found on the magnetic criterion.     

NMR shielding tensors for density fitted local second-order Mo̸ller-Plesset perturbation theory using gauge including atomic orbitals

An efficient method for the calculation of nuclear magnetic resonance (NMR) shielding tensors is presented, which treats electron correlation at the level of second-order Mo?ller-Plesset perturbation theory. It uses spatially localized functions to span occupied and virtual molecular orbital spaces, respectively, which are expanded in a basis of gauge including atomic orbitals (GIAOs or London atomic orbitals). Doubly excited determinants are restricted to local subsets of the virtual space and pair energies with an interorbital distance beyond a certain threshold are omitted. Furthermore, density fitting is employed to factorize the electron repulsion integrals. Ordinary Gaussians are employed as fitting functions. It is shown that the errors in the resulting NMR shielding constant, introduced (i) by the local approximation and (ii) by density fitting, are very small or even negligible. The capabilities of the new program are demonstrated by calculations on some extended molecular systems, such as the cyclobutane pyrimidine dimer photolesion with adjacent nucleobases in the native intrahelical DNA double strand (ATTA sequence). Systems of that size were not accessible to correlated ab initio calculations of NMR spectra before. The presented method thus opens the door to new and interesting applications in this area.     

Charge distribution and hydrogen bonding of a collagen α2‐chain in vacuum,hydrated, neutral,and charged structural models

A challenging task in computational biophysics is to ascertain the solvent effect on the electronic structure and interatomic bonding at the atomistic level. Simulations must be carried out on reasonably large biomolecules for accurate calculations to yield valid results. We report the results of a calculation on collagen model in the form of a peptide under three different environments: vacuum, solvated and with neutral and charged sites. Quantitative results and analysis of the partial charge (PC) distribution on each amino acid are discussed. A significant charge transfer of more than 1 electron from protein to water molecules is found with similar results when the model contains charged sites. The main contributions to the interatomic bonding are from hydrogen bonds (HBs) between water‐water and water‐protein pairs. A connection between PC and HBs can be established since the nonpolar amino acids form no HBs and have the smallest PC and vice versa. The ab initio PC obtained are used in the NAMD simulation showing significant improvement over the default values as reflected in the root mean square deviation of atomic positions in the MD steps and the total free energy in energy minimization. These results could facilitate the interpretation of data on interaction of various ligands in charged proteins in relation to isoelectric points. © 2016 Wiley Periodicals, Inc.     

Counter‐Rotating Spin‐Polarised Ring Currents in Odd‐Electron Carbocycles

We propose a molecular‐orbital model to explain how majority and minority spins in odd‐π‐electron carbocycles sustain counter‐rotating magnetic‐field‐induced ring currents. The model is based on the ipsocentric approach to magnetic response, in which ring currents are dominated by frontier‐orbital contributions obeying angular‐momentum selection rules. Coupled unrestricted Hartree–Fock ab initio calculations of the ring‐current responses for singly charged benzene and planarised cyclo‐octatetraene ions confirm the predictions of the qualitative model, and are consistent with correlated MP2 spin‐polarised current calculations.     

Relativistic Spin–Orbit Electronegativity and the Chemical Bond Between a Heavy Atom and a Light Atom

Relativistic effects are known to alter the chemical bonds and spectroscopic properties of heavy-element compounds. In this work, we introduce the concept of spin–orbit (SO) electronegativity of a heavy atom, as reflected by an SO-induced change in the interatomic distance between the heavy atom (HA) and a neighboring light atom (LA). We provide a transparent interpretation of these SO effects by using the concept of spin–orbit electron deformation density (SO-EDD). Spin–orbit coupling at the HA induces rearrangement of the electron density for the scalar-relativistically optimized geometry that, in turn, exerts a new force on the LA. The resulting expansion or contraction of the HA−LA bond depends on the nature and electron configuration of the HA. In addition, we quantify the change in atomic electronegativity induced by SO coupling for a series of hydrides, thereby complementing the SO-EDD picture. The trends in the SO-induced electronegativity and the HA−LA bond length across the periodic table of elements are demonstrated and interpreted, and also linked, intuitively, with the SO-induced NMR shielding at the LA.     

Local and nonlocal contributions to molecular first-order hyperpolarizability: a Hirshfeld partitioning analysis

Based on first-principles calculations, a decomposition scheme is proposed to investigate the molecular site-specific first-order hyperpolarizability (β) responses by means of Hirshfeld population analysis and finite field method. For a molecule, its β is decomposed into local and nonlocal contributions of individual atoms or groups. The former describes the response within the atomic sphere, while the latter describes the contributions from interatomic charge transfer. This scheme is then applied to six prototypical donor-acceptor (D-A) or D-π-A molecules for which the local and nonlocal hyperpolarizabilities are evaluated based on their MP2 density. Both the local and nonlocal parts exhibit site-specific characteristics, but vary differently with molecular structures. The local part depends mainly on the atomic attributes such as electronegativity and charge state, as well as its location in the molecule, while the nonlocal part relates to the ability and distance of charge delocalization within the molecule, increasing rapidly with molecular size. The proposed decomposition scheme provides a way to distinguish atomic or group contributions to molecular hyperpolarizabilities, which is useful in the molecular design for organic nonlinear optical materials.     

A reference‐free stockholder partitioning method based on the force on electrons

We argue that when one divides a molecular property into atom‐in‐a‐molecule contributions, one should perform the division based on the property density of the quantity being partitioned. This is opposition to the normal approach, where the electron density is given a privileged role in defining the properties of atoms‐in‐a‐molecule. Because partitioning each molecular property based on its own property density is inconvenient, we design a reference‐free approach that does not (directly) refer atomic property densities. Specifically, we propose a stockholder partitioning method based on relative influence of a molecule's atomic nuclei on the electrons at a given point in space. The resulting method does not depend on an “arbitrary” choice of reference atoms and it has some favorable properties, including the fact that all of the electron density at an atomic nucleus is assigned to that nucleus and the fact all the atoms in a molecule decay at a uniform asymptotic rate. Unfortunately, the resulting model is not easily applied to spatially degenerate ground states. Furthermore, the practical realizations of this strategy that we tried here gave disappointing numerical results. © 2017 Wiley Periodicals, Inc.