On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment

Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.

Self-consistent field method for open-shell systems within the density-matrix functional theory

The unrestricted Hartree-Fock method is extended to correlation calculation within the density-matrix functional theory. The method is derived from an entropic cumulant functional for the correlation energy. The eigenvalue equations for the spin-orbitals are modified by the orbital occupation numbers. The Euler equation for the occupation numbers results in the Fermi-Dirac distribution, which is very efficient to update as soon as the orbital eigenvalue equations are solved. The method is demonstrated on the ground state of O2 .

Dynamical Screening of Local Spin Moments at Metal-Molecule Interfaces

Transition-metal phthalocyanine molecules have attracted considerable interest in the context of spintronics device development due to their amenability to diverse bonding regimes and their intrinsic magnetism. The latter is highly influenced by the quantum fluctuations that arise at the inevitable metal-molecule interface in a device architecture. In this study, we have systematically investigated the dynamical screening effects in phthalocyanine molecules hosting a series of transition-metal ions (Ti, V, Cr, Mn, Fe, Co, and Ni) in contact with the Cu(111) surface. Using comprehensive density functional theory plus Anderson's Impurity Model calculations, we show that the orbital-dependent hybridization and electron correlation together result in strong charge and spin fluctuations. While the instantaneous spin moments of the transition-metal ions are near atomic-like, we find that screening gives rise to considerable lowering or even quenching of these. Our results highlight the importance of quantum fluctuations in metal-contacted molecular devices, which may influence the results obtained from theoretical or experimental probes, depending on their possibly material-dependent characteristic sampling time-scales.

DFT with corrections for an efficient and accurate description of strong electron correlations in NiO

An efficient and accurate description of the electronic structure of a strongly correlated metal-oxide semiconductor like NiO has been notoriously difficult. Here, we study the capabilities and limitations of two frequently employed correction schemes, a DFT+Uon-site correction and a DFT+1/2 self-energy correction. While both methods individually are unable to provide satisfactory results, in combination they provide a very good description of all relevant physical quantities. Since both methods cope with different shortcomings of common density-functional theory (DFT) methods (using local-density or generalized-gradient approximations), their combination is not mutually dependent and remains broadly applicable. The combined approach retains the computational efficiency of DFT calculations while providing significantly improved predictive power.

Monte Carlo study of cuprate superconductors in a four-bandd-pmodel: role of orbital degrees of freedom

Understanding the various competing phases in cuprate superconductors is a long-standing challenging problem. Recent studies have shown that orbital degrees of freedom, both Cue_gorbitals and Oporbitals, are a key ingredient for a unified understanding of cuprate superconductors, including the material dependence. Here we investigate a four-bandd-pmodel derived from the first-principles calculations with the variational Monte Carlo method, which allows us to elucidate competing phases on an equal footing. The obtained results can consistently explain the doping dependence of superconductivity, antiferromagnetic and stripe phases, phase separation in the underdoped region, and also novel magnetism in the heavily-overdoped region. The presence ofporbitals is critical to the charge-stripe features, which induce two types of stripe phases withs)-wave andd-wave bond stripe. On the other hand, the presence ofdz2orbital is indispensable to material dependence of the superconducting transition temperature (Tc), and enhances local magnetic moment as a source of novel magnetism in the heavily-overdoped region as well. These findings beyond one-band description could provide a major step toward a full explanation of unconventional normal state and highTcin cuprate supercondutors.

On the highest oxidation states of the actinoids in AnO4 molecules (An = Ac - Cm): A DMRG-CASSCF study

Actinoid tetroxide molecules AnO₄ (An = Ac - Cm) are investigated with the ab initio density matrix renormalization group (DMRG) approach. Natural orbital shapes are used to read out the oxidation state (OS) of the f-elements, and the atomic orbital energies and radii are used to explain the trends. The highest OSs reveal a "volcano"-type variation: For An = Ac - Np, the OSs are equal to the number of available valence electrons, that is, Ac^III , Th^IV , Pa^V , U^VI , and Np^VII . Starting with plutonium as the turning point, the highest OSs in the most stable AnO₄ isomers then decrease as Pu^V , Am^V , and Cm^III , indicating that the 5f-electrons are hard to be fully oxidized off from Pu onward. The variations are related to the actinoid contraction and to the 5f-covalency characteristics. Combined with previous work on OSs, we review their general trends throughout the periodic table, providing fundamental understanding of OS-relevant phenomena.

Enhanced Electron Correlation and Significantly Suppressed Thermal Conductivity in Dirac Nodal-Line Metal Nanowires by Chemical Doping

Enhancing electron correlation in a weakly interacting topological system has great potential to promote correlated topological states of matter with extraordinary quantum properties. Here, the enhancement of electron correlation in a prototypical topological metal, namely iridium dioxide (IrO₂ ), via doping with 3d transition metal vanadium is demonstrated. Single-crystalline vanadium-doped IrO₂ nanowires are synthesized through chemical vapor deposition where the nanowire yield and morphology are improved by creating rough surfaces on substrates. Vanadium doping leads to a dramatic decrease in Raman intensity without notable peak broadening, signifying the enhancement of electron correlation. The enhanced electron correlation is further evidenced by transport studies where the electrical resistivity is greatly increased and follows an unusual T $\sqrt T $ dependence on the temperature (T). The lattice thermal conductivity is suppressed by an order of magnitude via doping even at room temperature where phonon-impurity scattering becomes less important. Density functional theory calculations suggest that the remarkable reduction of thermal conductivity arises from the complex phonon dispersion and reduced energy gap between phonon branches, which greatly enhances phase space for phonon-phonon Umklapp scattering. This work demonstrates a unique system combining 3d and 5d transition metals in isostructural materials to enrich the system with various types of interactions.

Enhancing Electron Correlation at a 3d Ferromagnetic Surface

Spin-resolved momentum microscopy and theoretical calculations are combined beyond the one-electron approximation to unveil the spin-dependent electronic structure of the interface formed between iron (Fe) and an ordered oxygen (O) atomic layer, and an adsorbate-induced enhancement of electronic correlations is found. It is demonstrated that this enhancement is responsible for a drastic narrowing of the Fe d-bands close to the Fermi energy (E_F ) and a reduction of the exchange splitting, which is not accounted for in the Stoner picture of ferromagnetism. In addition, correlation leads to a significant spin-dependent broadening of the electronic bands at higher binding energies and their merging with satellite features, which are manifestations of a pure many-electron behavior. Overall, adatom adsorption can be used to vary the material parameters of transition metal surfaces to access different intermediate electronic correlated regimes, which will otherwise not be accessible. The results show that the concepts developed to understand the physics and chemistry of adsorbate-metal interfaces, relevant for a variety of research areas, from spintronics to catalysis, need to be reconsidered with many-particle effects being of utmost importance. These may affect chemisorption energy, spin transport, magnetic order, and even play a key role in the emergence of ferromagnetism at interfaces between non-magnetic systems.

Polycyclic Hydrocarbons from [4n]Annulenes: Correlation versus Hybridization Forces in the Formation of Diradicaloids

The conceptual connections between [4n] Hückel antiaromaticity, disjoint orbitals, correlation energy, pro-aromaticity and diradical character for a variety of extended π-conjugated systems, including some salient recent examples of nanographenes and polycyclic aromatic radicals, are provided based on their [4n]annulene peripheries. The realization of such structure-property relationships has led to a beneficial pedagogic exercise establishing design guidelines for diradicaloids. The antiaromatic fingerprint of the [4n]annulene peripheries upon orbital interactions due to internal covalent connectors gives insights into the diradicaloid property of a diversity of π-conjugated molecules that have fascinated chemists recently.

Photoionization Observables from Multi-Reference Dyson Orbitals Coupled to B-Spline DFT and TD-DFT Continuum

We present a theoretical model to compute the accurate photoionization dynamical parameters (cross-sections, asymmetry parameters and orbital, or cross-section, ratios) from Dyson orbitals obtained with the multi-state complete active space perturbation theory to the second order (MS-CASPT2) method. Our new implementation of Dyson orbitals in OpenMolcas takes advantage of the full Abelian symmetry point group and has the corrected normalization. The Dyson orbitals are coupled to an accurate description of the electronic continuum obtained with a multicentric B-spline basis at the DFT and TD-DFT levels. Two prototype diatomic molecules, i.e., CS and SiS, have been chosen due to their smallness, which hides important correlation effects. These effects manifest themselves in the appearance of well-characterized isolated satellite bands in the middle of the valence region. The rich satellite structures make CS and SiS the perfect candidates for a computational study based on our highly accurate MS-CASPT2/B-spline TD-DFT protocol.

Reparametrization of the Colle-Salvetti formula

We investigate the Colle-Salvetti (CS) formula, the basis of the Lee, Yang and Parr (LYP) correlation functional used in approximate density functional theory. The CS formula is reparametrized using high-accuracy Hartree-Fock (HF) wavefunctions to determine the accuracy of the formula to calculate anions. Fitting to the hydride ion or the two-electron system just prior to electron detachment at the HF level of theory does not, in general, improve the calculated correlation energies using the parameters derived from the CS/LYP method. An analysis of the CS parameters used in the popular LYP functional demonstrates the ingenuity and perhaps fortuitousness of the original formulation by CS.

Energy-based automatic determination of buffer region in the divide-and-conquer second-order Møller-Plesset perturbation theory

In the linear‐scaling divide‐and‐conquer (DC) electronic structure method, each subsystem is calculated together with the neighboring buffer region, the size of which affects the energy error introduced by the fragmentation in the DC method. The DC self‐consistent field calculation utilizes a scheme to automatically determine the appropriate buffer region that is as compact as possible for reducing the computational time while maintaining acceptable accuracy (J. Comput. Chem. 2018, 39, 909). To extend the automatic determination scheme of the buffer region to the DC second‐order Møller–Plesset perturbation (MP2) calculation, a scheme for estimating the subsystem MP2 correlation energy contribution from each atom in the buffer region is proposed. The estimation is based on the atomic orbital Laplace MP2 formalism. Based on this, an automatic buffer determination scheme for the DC‐MP2 calculation is constructed and its performance for several types of systems is assessed.

Fragmentation of natural orbital occupation numbers-based diagnostic of differential multireference character in complexes with hydrogen bonds

We explore the possible route to approximate natural orbital occupation numbers-based diagnostic of differential multireference character of noncovalent energy differences by techniques based on many-body expansion. It turns out that two-body fragmentation of monomers may lead to a reasonable approximation of such a diagnostic in hydrogen-bonded complexes. The results are useful, for example, for assessment of the expected bias cancellation in energy differences of larger systems obtained by single-reference methods.

A Valence-Bond-Based Multiconfigurational Density Functional Theory: The λ-DFVB Method Revisited

A recently developed valence-bond-based multireference density functional theory, named λ-DFVB, is revisited in this paper. λ-DFVB remedies the double-counting error of electron correlation by decomposing the electron-electron interactions into the wave function term and density functional term with a variable parameter λ. The λ value is defined as a function of the free valence index in our previous scheme, denoted as λ-DFVB(K) in this paper. Here we revisit the λ-DFVB method and present a new scheme based on natural orbital occupation numbers (NOONs) for parameter λ, named λ-DFVB(IS), to simplify the process of λ-DFVB calculation. In λ-DFVB(IS), the parameter λ is defined as a function of NOONs, which are straightforwardly determined from the many-electron wave function of the molecule. Furthermore, λ-DFVB(IS) does not involve further self-consistent field calculation after performing the valence bond self-consistent field (VBSCF) calculation, and thus, the computational effort in λ-DFVB(IS) is approximately the same as the VBSCF method, greatly reduced from λ-DFVB(K). The performance of λ-DFVB(IS) was investigated on a broader range of molecular properties, including equilibrium bond lengths and dissociation energies, atomization energies, atomic excitation energies, and chemical reaction barriers. The computational results show that λ-DFVB(IS) is more robust without losing accuracy and comparable in accuracy to high-level multireference wave function methods, such as CASPT2.

Combining EXAFS and Computer Simulations to Refine the Structural Description of Actinyls in Water

EXAFS spectroscopy is one of the most used techniques to solve the structure of actinoid solutions. In this work a systematic analysis of the EXAFS spectra of four actinyl cations, [UO2]2+, [NpO2]2+, [NpO2]+ and [PuO2]2+ has been carried out by comparing experimental results with theoretical spectra. These were obtained by averaging individual contributions from snapshots taken from classical Molecular Dynamics simulations which employed a recently developed [AnO2]2+/+ -H2O force field based on the hydrated ion model using a quantum-mechanical (B3LYP) potential energy surface. Analysis of the complex EXAFS signal shows that both An-Oyl and An-OW single scattering paths as well as multiple scattering ones involving [AnO2]+/2+ molecular cation and first-shell water molecules are mixed up all together to produce a very complex signal. Simulated EXAFS from the B3LYP force field are in reasonable agreement for some of the cases studied, although the k= 6-8 Å-1 region is hard to be reproduced theoretically. Except uranyl, all studied actinyls are open-shell electron configurations, therefore it has been investigated how simulated EXAFS spectra are affected by minute changes of An-O bond distances produced by the inclusion of static and dynamic electron correlation in the quantum mechanical calculations. A [NpO2]+-H2O force field based on a NEVPT2 potential energy surface has been developed. The small structural changes incorporated by the electron correlation on the actinyl aqua ion geometry, typically smaller than 0.07 Å, leads to improve the simulated spectrum with respect to that obtained from the B3LYP force field. For the other open-shell actinyls, [NpO2]2+ and [PuO2]2+, a simplified strategy has been adopted to improve the simulated EXAFS spectrum. It is computed taking as reference structure the NEVPT2 optimized geometry and including the DW factors of their corresponding MD simulations employing the B3LYP force field. A better agreement between the experimental and the simulated EXAFS spectra is found, confirming the a priori guess that the inclusion of dynamic and static correlation refine the structural description of the open-shell actinyl aqua ions.

First principles investigation of electron correlation and Lifshitz transition within iron polynitrides

Metal poly-nitrogen compounds are gaining great interests as potential high energy density materials. Several iron polynitrides have been recently synthesized and investigated under high pressure (2018Nature Communications92756). In this work the electron correlations within these iron poly-nitrogen compounds were self-consistently determined, benchmarked with those obtained from linear response approach. Along with the increase of the concentration of nitrogen, the Coulomb interaction strengths show a monotonic decrease, where FeN and FeN₂are antiferromagnetic and the others are ferromagnetic. Then the electron correlation is studied along with the pressure, where the electrons are more delocalized as pressure becomes higher. One electronic topological transition was found for FeN₂, owing to a breaking of symmetry of spin and a transition of magnetism induced by a structural change. The band structure, densities of states, Fermi surface and absorption spectra were calculated and discussed.

Electron correlation effect versus spin-orbit coupling for tungsten and impurities

The electron correlation and spin-orbit coupling (SOC) effects are investigated for body-centered-cubic tungsten and intrinsic & irradiative impurities using first principles calculations based upon the density functional theory. It is found that the electron correlation between the localized 5delectrons and the SOC effect are significant in modifying the band structures and the formation energies of defects. For the latter one, the involving of electron correlation always makes the defects stabler than the Perdew-Burke-Ernzerhof results, while the SOC contributes diversely for different defects. Moreover, the migration barrier of single tungsten vacancy moving in ⟨111⟩ direction is explored, where the inclusion of electron correlations remarkably decreases the migration barrier, while the influence of SOC is almost negligible. This study can help to validate the previous studies on irradiative defects in tungsten and improve the further investigations.

Electronic, magnetic and optical properties of transition-metal and hydroxides doped monolayer g-C3N4: a first principles investigation

The graphitic carbon nitride (g-C3N4) is a promising layered two-dimension material with an opened bandgap. It is of interest to explore the tunability of the bandgap together with the magnetism by doping transition metal atoms. In this work, we investigated the transition metals (Mn, Fe, Co, Ni) and their hydroxides doped g-C3N4monolayers. The electron correlations between the 3delectrons of the doped transition metal atoms are self-consistently calculated and analyzed based on the density functional theory. The magnetism, electronic band structures and optical properties are systematically investigated. It reveals that the transition metal doped g-C3N4is ferromagnetic (FM) state at small doping concentration, where the two spins show different bandgaps. When the doping is high enough, it turns to metallic antiferromagnetic (AFM) state except that Mn doped g-C3N4is metallic FM state. On another hand, the system shows variable absorption spectra at different doping level. When the vacancy sites are fully occupied, a large absorption peak appears around 1.5 eV suitable for visible light. Moreover, within the transition metal hydroxides doped g-C3N4, the global ground state shows as AFM, and the absorption spectra within low energy range is distinct due to the presence of hydroxyl group. Therefore, doping with transition metal atoms and hydroxides can effectively tune the bandgap, magnetism and optical properties of g-C3N4so as to promote its applications.

Spin-dependent transient current in transistor-like nanostructures

Transient current in transistor-like nanostructures has been studied by a model of a few electrons confined in a one-dimensional effective potential consisting of three quantum wells, 'source', 'gate', and 'drain'. The time-dependent Schrödinger equation for the electrons has been integrated relying on the symplectic integrator method and the transient current has been calculated as the flux of the probability density of electrons absorbed by the complex absorbing potential placed at the far edge of the drain region. The electrons are initially placed in the source domain as their lowest-energy state for a given spin multiplicity and the source-drain current has been calculated for different gate potential heights. The current for different spin configurations has shown strong emission at different values of the gate potential, suggesting use of the studied nanostructures for extracting current with a specific spin configuration from spin-unpolarized normal current. Dependence of the current emission on electron correlation has also been studied by changing the size of the source domain. The current has shown appreciable differences for different spin configurations for the medium and strong confinement regimes, while these differences become smaller for smaller confinement and tend to diminish in the weak limit of confinement. This observed trend has been rationalized on the basis of the formation of the Wigner lattice states.

Cyclotron Resonance Study of Monolayer Graphene under Double Moiré Potentials

We report the first cyclotron resonance study of monolayer graphene under double-moiré potentials in which the crystal axis of graphene is nearly aligned to those of both the top and bottom hexagonal boron nitride (h-BN) layers. Under mid-infrared light irradiation, we observe cyclotron resonance absorption with the following unique features: (1) cyclotron resonance magnetic field B_CR is entirely different from that of nonaligned monolayer graphene, (2) B_CR exhibits strong electron-hole asymmetry, and (3) splitting of B_CR is observed for |ν| < 1, with the split maximum at |ν| = 1, resulting in eyeglass-shaped trajectories. These features are well explained by considering the large bandgap induced by the double moiré potentials, the electron-hole asymmetry in the Fermi velocity, and the Fermi-level-dependent enhancement of spin gaps, which suggests a large electron-electron correlation contribution in this system.

An Interacting Quantum Atoms (IQA) and Relative Energy Gradient (REG) Study of the Halogen Bond with Explicit Analysis of Electron Correlation

Energy profiles of seven halogen-bonded complexes were analysed with the topological energy partitioning called Interacting Quantum Atoms (IQA) at MP4(SDQ)/6-31+G(2d,2p) level of theory. Explicit interatomic electron correlation energies are included in the analysis. Four complexes combine X₂ (X = Cl or F) with HCN or NH₃, while the remaining three combine ClF with HCN, NH₃ or N₂. Each complex was systematically deformed by translating the constituent molecules along its central axis linking X and N, and reoptimising its remaining geometry. The Relative Energy Gradient (REG) method (Theor. Chem. Acc. 2017, 136, 86) then computes which IQA energies most correlate with the total energy during the process of complex formation and further compression beyond the respective equilibrium geometries. It turns out that the covalent energy (i.e., exchange) of the halogen bond, X…N, itself drives the complex formation. When the complexes are compressed from their equilibrium to shorter X…N distance then the intra-atomic energy of N is in charge. When the REG analysis is restricted to electron correlation then the interatomic correlation energy between X and N again drives the complex formation, and the complex compression is best described by the destabilisation of the through-space correlation energy between N and the "outer" halogen.

Electron-Doping Mottronics in Strongly Correlated Perovskite

The discovery of hydrogen-induced electron localization and highly insulating states in d-band electron correlated perovskites has opened a new paradigm for exploring novel electronic phases of condensed matters and applications in emerging field-controlled electronic devices (e.g., Mottronics). Although a significant understanding of doping-tuned transport properties of single crystalline correlated materials exists, it has remained unclear how doping-controlled transport properties behave in the presence of planar defects. The discovery of an unexpected high-concentration doping effect in defective regions is reported for correlated nickelates. It enables electronic conductance by tuning the Fermi-level in Mott-Hubbard band and shaping the lower Hubbard band state into a partially filled configuration. Interface engineering and grain boundary designs are performed for H_x SmNiO₃ /SrRuO₃ heterostructures, and a Mottronic device is achieved. The interfacial aggregation of hydrogen is controlled and quantified to establish its correlation with the electrical transport properties. The chemical bonding between the incorporated hydrogen with defective SmNiO₃ is further analyzed by the positron annihilation spectroscopy. The present work unveils new materials physics in correlated materials and suggests novel doping strategies for developing Mottronic and iontronic devices via hydrogen-doping-controlled orbital occupancy in perovskite heterostructures.

Atomic Partitioning of the MPn (n = 2, 3, 4) Dynamic Electron Correlation Energy by the Interacting Quantum Atoms Method: A Fast and Accurate Electrostatic Potential Integral Approach

Recently, the quantum topological energy partitioning method called interacting quantum atoms (IQA) has been extended to MPn (n = 2, 3, 4) wave functions. This enables the extraction of chemical insight related to dynamic electron correlation. The large computational expense of the IQA-MPn approach is compensated by the advantages that IQA offers compared to older nontopological energy decomposition schemes. This expense is problematic in the construction of a machine learning training set to create kriging models for topological atoms. However, the algorithm presented here markedly accelerates the calculation of atomically partitioned electron correlation energies. Then again, the algorithm cannot calculate pairwise interatomic energies because it applies analytical integrals over whole space (rather than over atomic volumes). However, these pairwise energies are not needed in the quantum topological force field FFLUX, which only uses the energy of an atom interacting with all remaining atoms of the system that it is part of. Thus, it is now feasible to generate accurate and sizeable training sets at MPn level of theory. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Relativistic quantum crystallography of diphenyl- and dicyanomercury. Theoretical structure factors and Hirshfeld atom refinement

Quantum crystallographic refinement of heavy-element-containing compounds is a challenge, because many physical effects have to be accounted for adequately. Here, the impact and magnitude of relativistic effects are compared with those of electron correlation, polarization through the environment, choice of basis set and treatment of thermal motion effects on the structure factors of diphenylmercury(II) [Hg(Ph)₂] and dicyanomercury(II) [Hg(CN)₂]. Furthermore, the individual atomic contributions to the structure factors are explored in detail (using Mulliken population analysis and the exponential decay of atomic displacement parameters) to compare the contributions of lighter atoms, especially hydrogen atoms, against mercury. Subsequently, relativistic Hirshfeld atom refinement (HAR) is validated against theoretical structure factors of Hg(Ph)₂ and Hg(CN)₂, starting from perturbed geometries, to test if the relativistic variant of HAR leads to multiple solutions. Generally, relativistic HAR is successful, leading to a perfect match with the reference geometries, but some limitations are pointed out.

The Coulomb Hole of the Ne Atom

We analyze the Coulomb hole of Ne from highly-accurate CISD wave functions obtained from optimized even-tempered basis sets. Using a two-fold extrapolation procedure we obtain highly accurate results that recover 97 % of the correlation energy. We confirm the existence of a shoulder in the short-range region of the Coulomb hole of the Ne atom, which is due to an internal reorganization of the K-shell caused by electron correlation of the core electrons. The feature is very sensitive to the quality of the basis set in the core region and it is not exclusive to Ne, being also present in most of second-row atoms, thus confirming that it is due to K-shell correlation effects.

Electron Transport in Carbon Nanotubes with Adsorbed Chromium Impurities

We employ Green’s function method for describing multiband models with magnetic impurities and apply the formalism to the problem of chromium impurities adsorbed onto a carbon nanotube. Density functional theory is used to determine the bandstructure, which is then fit to a tight-binding model to allow for the subsequent Green’s function description. Electron–electron interactions, electron–phonon coupling, and disorder scattering are all taken into account (perturbatively) with a theory that involves a cluster extension of the coherent potential approximation. We show how increasing the cluster size produces more accurate results and how the final calculations converge as a function of the cluster size. We examine the spin-polarized electrical current on the nanotube generated by the magnetic impurities adsorbed onto the nanotube surface. The spin polarization increases with both increasing concentration of chromium impurities and with increasing magnetic field. Its origin arises from the strong electron correlations generated by the Cr impurities.

Electron correlation in Li+, He, H- and the critical nuclear charge system Z C : energies, densities and Coulomb holes

This paper presents high-accuracy correlation energies, intracule densities and Coulomb hole(s) for the lithium cation, helium, hydride ion and the system with the critical nuclear charge, Z C , for binding two electrons. The fully correlated (FC) wave function and the Hartree-Fock (HF) wave function are both determined using a Laguerre-based wave function. It is found that for the lithium cation and the helium atom a secondary Coulomb hole is present, in agreement with a previous literature finding, confirming a counterintuitive conclusion that electron correlation can act to bring distant electrons closer together. However, no evidence for a tertiary Coulomb hole is found. For the hydride anion and the system just prior to electron detachment only a single Coulomb hole is present and electron correlation decreases the probability of finding the electrons closer together at all radial distances. The emergence of a secondary Coulomb hole is investigated and found to occur between Z = 1.15 and Z = 1.20. The FC and HF energies and intracule densities (in atomic units) used to calculate the correlation energy and Coulomb hole, respectively, are accurate to at least the nano-scale for helium and the cation and at least the micro-scale for the anions.

RI-MP3 calculations of biomolecules based on the fragment molecular orbital method

In this study, the third-order Møller-Plesset perturbation (MP3) theory using the resolution of the identity (RI) approximation was combined with the fragment molecular orbital (FMO) method to efficiently calculate a high-order electron correlation energy of biomolecular systems. We developed a new algorithm for the RI-MP3 calculation, which can be used with the FMO scheme. After test calculations using a small molecule, the FMO-RI-MP3 calculations were performed for two biomolecular systems comprising a protein and a ligand. The computational cost of these calculations was only around 5 and 4 times higher than those of the FMO-RHF calculations. The error associated with the RI approximation was around 2.0% of the third-order correlation contribution to the total energy. However, the RI approximation error in the interaction energy between the protein and ligand molecule was insignificantly small, which reflected the negligible error in the inter fragment interaction energy. © 2018 Wiley Periodicals, Inc.

Hartree-Fock implementation using a Laguerre-based wave function for the ground state and correlation energies of two-electron atoms

An implementation of the Hartree-Fock (HF) method using a Laguerre-based wave function is described and used to accurately study the ground state of two-electron atoms in the fixed nucleus approximation, and by comparison with fully correlated (FC) energies, used to determine accurate electron correlation energies. A variational parameter A is included in the wave function and is shown to rapidly increase the convergence of the energy. The one-electron integrals are solved by series solution and an analytical form is found for the two-electron integrals. This methodology is used to produce accurate wave functions, energies and expectation values for the helium isoelectronic sequence, including at low nuclear charge just prior to electron detachment. Additionally, the critical nuclear charge for binding two electrons within the HF approach is calculated and determined to be Z^HF_C=1.031 177 528.This article is part of the theme issue 'Modern theoretical chemistry'.

Conceptual Insights into DFT Spin-State Energetics of Octahedral Transition-Metal Complexes through a Density Difference Analysis

In this study, an intuitive concept is derived, which explains the characteristic dependence of spin-state energetics on the exact exchange admixture of DFT functionals in the case of octahedral transition metal complexes. The change in electron density distributions upon varying the admixture, c₃ , in the B3LYP functional is analyzed for archetype ionic and covalent systems as well as for the Fe²⁺ ion in an ideal octahedral field. An understanding of how the DFT description of the electronic structure of octahedral complexes changes as a function of c₃ is sought. A systematic spin-state energy analysis of 50 octahedral complexes of various metals and ligands with consistent experimental data is presented, allowing the derivation, in theory, of an optimal c₃ value for each system. The notion that the admixture dependence of DFT spin-state energetics stems from the treatment of nondynamic electrons arising from the mixing of (M-Lz2 )⁰ (dz2 )² and (M-Lx2-y2 )⁰ (dx2-y2 )² configurations into the dominant (M-Lx2-y2 )² (dx2-y2 )⁰ and (M-Lx2-y2 )² (dx2-y2 )⁰ ones in the low(er) spin states is put forward. That is, in the effort to mimic such electron-electron interactions, E_x^LDA overestimates, whereas exact exchange downplays the contribution of this type of electron correlation to the stability of low(er) spin states, leading to the widespread practical observation that the higher the exact exchange admixture, the more stable the high-spin-state configuration.

The Transferability of Topologically Partitioned Electron Correlation Energies in Water Clusters

The electronic effects that govern the cohesion of water clusters are complex, demanding the inclusion of N-body, Coulomb, exchange and correlation effects. Here we present a much needed quantitative study of the effect of correlation (and hence dispersion) energy on the stabilization of water clusters. For this purpose we used a topological energy partitioning method called Interacting Quantum Atoms (IQA) to partition water clusters into topological atoms, based on a MP2/6-31G(d,p) wave function, and modified versions of GAUSSIAN09 and the Quantum Chemical Topology (QCT) program MORFI. Most of the cohesion in the water clusters provided by electron correlation comes from intramolecular energy stabilization. Hydrogen bond-related interactions tend to largely cancel each other. Electron correlation energies are transferable in almost all instances within 1 kcal mol^-1 . This observed transferability is very important to the further development of the QCT force field FFLUX, especially to the future modelling of liquid water.

Electronic structure, transport, and collective effects in molecular layered systems

The great potential of organic heterostructures for organic device applications is exemplified by the targeted engineering of the electronic properties of phthalocyanine-based systems. The transport properties of two different phthalocyanine systems, a pure copper phthalocyanine (CoPc) and a flourinated copper phthalocyanine-manganese phthalocyanine (F16CoPc/MnPc) heterostructure, are investigated by means of density functional theory (DFT) and the non-equilibrium Green's function (NEGF) approach. Furthermore, a master-equation-based approach is used to include electronic correlations beyond the mean-field-type approximation of DFT. We describe the essential theoretical tools to obtain the parameters needed for the master equation from DFT results. Finally, an interacting molecular monolayer is considered within a master-equation approach.

Diversity of Chemical Bonding and Oxidation States in MS4 Molecules of Group 8 Elements

The geometric and electronic ground-state structures of 30 isomers of six MS₄ molecules (M=Group 8 metals Fe, Ru, Os, Hs, Sm, and Pu) have been studied by using quantum-chemical density functional theory and correlated wavefunction approaches. The MS₄ species were compared to analogous MO₄ species recently investigated (W. Huang, W.-H. Xu, W. H. E. Schwarz, J. Li, Inorg. Chem. 2016, 55, 4616). A metal oxidation state (MOS) with a high value of eight appeared in the low-spin singlet T_d geometric species (Os,Hs)S₄ and (Ru,Os,Hs)O₄ , whereas a low MOS of two appeared in the high-spin septet D_2d species Fe(S₂ )₂ and (slightly excited) metastable Fe(O₂ )₂ . The ground states of all other molecules had intermediate MOS values, with S^2- , S₂^2- , S₂^1- (and O^2- , O^1- , O₂^2- , O₂^1- ) ligands bonded by ionic, covalent, and correlative contributions. The known tendencies toward lower MOS on going from oxides to sulfides, from Hs to Os to Ru, and from Pu to Sm, and the specific behavior of Fe, were found to arise from the different atomic orbital energies and radii of the (n-1)p core and (n-1)d and (n-2)f valence shells of the metal atoms in row n of the periodic table. The comparative results of the electronic and geometric structures of the MO₄ and MS₄ species provides insight into the periodicity of oxidation states and bonding.

Can X-ray constrained Hartree-Fock wavefunctions retrieve electron correlation?

The X-ray constrained wavefunction (XC-WF) method proposed by Jayatilaka [Jayatilaka & Grimwood (2001) ▸, Acta Cryst. A57, 76-86] has attracted much attention because it represents a possible third way of theoretically studying the electronic structure of atoms and molecules, combining features of the more popular wavefunction- and DFT-based approaches. In its original formulation, the XC-WF technique extracts statistically plausible wavefunctions from experimental X-ray diffraction data of molecular crystals. A weight is used to constrain the pure Hartree-Fock solution to the observed X-ray structure factors. Despite the wavefunction being a single Slater determinant, it is generally assumed that its flexibility could guarantee the capture, better than any other experimental model, of electron correlation effects, absent in the Hartree-Fock Hamiltonian but present in the structure factors measured experimentally. However, although the approach has been known for long time, careful testing of this fundamental hypothesis is still missing. Since a formal demonstration is impossible, the validation can only be done heuristically and, to accomplish this task, X-ray constrained Hartree-Fock calculations have been performed using structure factor amplitudes computed at a very high correlation level (coupled cluster) for selected molecules in isolation, in order to avoid the perturbations due to intermolecular interactions. The results show that a single-determinant XC-WF is able to capture the electron correlation effects only partially. The largest amount of electron correlation is extracted when: (i) a large external weight is used (much larger than what has normally been used in XC-WF calculations using experimental data); and (ii) the high-order reflections, which carry less information on the electron correlation, are down-weighted (or even excluded), otherwise they would bias the fitting towards the unconstrained Hartree-Fock wavefunction.

Unconventional high-Tc superconductivity in fullerides

A3C60 molecular superconductors share a common electronic phase diagram with unconventional high-temperature superconductors such as the cuprates: superconductivity emerges from an antiferromagnetic strongly correlated Mott-insulating state upon tuning a parameter such as pressure (bandwidth control) accompanied by a dome-shaped dependence of the critical temperature, Tc However, unlike atom-based superconductors, the parent state from which superconductivity emerges solely by changing an electronic parameter-the overlap between the outer wave functions of the constituent molecules-is controlled by the C60 (3-) molecular electronic structure via the on-molecule Jahn-Teller effect influence of molecular geometry and spin state. Destruction of the parent Mott-Jahn-Teller state through chemical or physical pressurization yields an unconventional Jahn-Teller metal, where quasi-localized and itinerant electron behaviours coexist. Localized features gradually disappear with lattice contraction and conventional Fermi liquid behaviour is recovered. The nature of the underlying (correlated versus weak-coupling Bardeen-Cooper-Schrieffer theory) s-wave superconducting states mirrors the unconventional/conventional metal dichotomy: the highest superconducting critical temperature occurs at the crossover between Jahn-Teller and Fermi liquid metal when the Jahn-Teller distortion melts.This article is part of the themed issue 'Fullerenes: past, present and future, celebrating the 30th anniversary of Buckminster Fullerene'.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

The diverse manifold of electronic states generated by a single carbon defect in a graphene sheet: multireference calculations using a pyrene defect model

Detailed calculations have been performed on the electronic states occurring in a single vacancy defect model based on pyrene from which one of the central carbon atoms has been removed. Complete active space self-consistent field and multireference configuration interaction with singles and doubles calculations have been performed using the 6-31G and 6-31G* basis sets. Two types of defect geometries have been defined: 1) The unrelaxed defect structure based on pyrene and 2) a relaxed structure. In total 12 electronic states have been computed for the unrelaxed structure at C2v symmetry, comprising four singlets, triplets and quintets each. The lowest six states are formed from singlet and triplet states and appear in a rather narrow gap of ∼0.6 eV. The lowest quintet state is found 1.43 eV above the (3) B1 ground state. As predicted from Jahn-Teller distortions, a CC bond is formed between dangling carbon bonds in the (1, 3) B1 states, leading to the formation of a five-membered ring. The (1, 3) A2 states show initial repulsive behavior along the bond formation coordinate until an avoided crossing is reached by which these states are furnished with CC bonding character so that finally also in these cases a CC bond is established. Linear interpolation curves between the initial unrelaxed defect structure and the final optimized structure are used to give an overview of the evolution of electronic states and the occurrence of avoided crossings. Out-of-plane structures are investigated with special emphasis on the carbon atom containing a dangling bond in the relaxed structure. Unpaired electron densities are used to characterize the electronic structure of the different states.

Cyclobutadiene automerization and rotation of ethylene: energetics of the barriers by using spin-polarized wave functions

Spin-projected spin polarized Møller-Plesset and spin polarized coupled clusters calculations have been made to estimate the cyclobutadiene automerization, the ethylene torsion barriers in their ground state, and the gap between the singlet and triplet states of ethylene. The results have been obtained optimizing the geometries at MP4 and/or CCSD levels, by an extensive Gaussian basis set. A comparative analysis with more complex calculations, up to MP5 and CCSDTQP, together with others from the literature, have also been made, showing the efficacy of using spin-polarized wave functions as a reference wave function for Møller-Plesset and coupled clusters calculations, in such problems.

Electrostatic Domination of the Effect of Electron Correlation in Intermolecular Interactions

The electron-electron correlation energy is negative, and attractive dispersion interactions are entirely a correlation effect; therefore, the contribution of correlation to intermolecular binding is commonly assumed to be negative, or binding in nature. However, there are many cases where the long-range correlation binding energy is positive, with certain geometries of the water dimer as a prominent example. Geometries with dipoles misaligned can also have an electrostatically dominated, though negative, long-range correlation binding. In either case, the interaction decays as R(-3). This has its origin in the systematic overestimation of dipole moments by Hartree-Fock theory, leading to a reduction in the calculated electrostatic attraction upon inclusion of correlation. Thus, energy decomposition analyses that include correlation but do not correct mean field electrostatic terms are suboptimal. Attenuated second-order Møller-Plesset theory, which smoothly truncates long-range electron correlation effects to zero, can, paradoxically, have the correct long-range behavior for many intermolecular interactions.

On the basis set convergence of electron-electron entanglement measures: helium-like systems

A systematic investigation of three different electron–electron entanglement measures, namely the von Neumann, the linear and the occupation number entropy at full configuration interaction level has been performed for the four helium-like systems hydride, helium, Li⁺ and Be²⁺ using a large number of different basis sets. The convergence behavior of the resulting energies and entropies revealed that the latter do in general not show the expected strictly monotonic increase upon increase of the one–electron basis. Overall, the three different entanglement measures show good agreement among each other, the largest deviations being observed for small basis sets. The data clearly demonstrates that it is important to consider the nature of the chemical system when investigating entanglement phenomena in the framework of Gaussian type basis sets: while in case of hydride the use of augmentation functions is crucial, the application of core functions greatly improves the accuracy in case of cationic systems such as Li⁺ and Be²⁺. In addition, numerical derivatives of the entanglement measures with respect to the nucleic charge have been determined, which proved to be a very sensitive probe of the convergence leading to qualitatively wrong results (i.e., the wrong sign) if too small basis sets are used.

Multidimensional attosecond resonant X-ray spectroscopy of molecules: lessons from the optical regime

New free-electron laser and high-harmonic generation X-ray light sources are capable of supplying pulses short and intense enough to perform resonant nonlinear time-resolved experiments in molecules. Valence-electron motions can be triggered impulsively by core excitations and monitored with high temporal and spatial resolution. We discuss possible experiments that employ attosecond X-ray pulses to probe the quantum coherence and correlations of valence electrons and holes, rather than the charge density alone, building on the analogy with existing studies of vibrational motions using femtosecond techniques in the visible regime.

Recent progress in high-pressure studies on organic conductors

Recent high-pressure studies of organic conductors and superconductors are reviewed. The discovery of the highest T_c superconductivity among organics under high pressure has triggered the further progress of the high-pressure research. Owing to this finding, various organic conductors with the strong electron correlation were investigated under high pressures. This review includes the pressure techniques using the cubic anvil apparatus, as well as high-pressure studies of the organic conductors up to 10 GPa showing extraordinary temperature and pressure dependent transport phenomena.

Iron pnictides as a new setting for quantum criticality

Two major themes in the physics of condensed matter are quantum critical phenomena and unconventional superconductivity. These usually occur in the context of competing interactions in systems of strongly correlated electrons. All this interesting physics comes together in the behavior of the recently discovered iron pnictide compounds that have generated enormous interest because of their moderately high-temperature superconductivity. The ubiquity of antiferromagnetic ordering in their phase diagrams naturally raises the question of the relevance of magnetic quantum criticality, but the answer remains uncertain both theoretically and experimentally. Here, we show that the undoped iron pnictides feature a unique type of magnetic quantum critical point, which results from a competition between electronic localization and itinerancy. Our theory provides a mechanism to understand the experimentally observed variation of the ordered moment among the undoped iron pnictides. We suggest P substitution for As in the undoped iron pnictides as a means to access this example of magnetic quantum criticality in an unmasked fashion. Our findings point to the iron pnictides as a much-needed setting for quantum criticality, one that offers a unique set of control parameters.