Search

Delving into the intricate interplay between spin-orbit coupling and Coulomb correlations in strongly correlated oxides, particularly perovskite compounds, has unveiled a rich landscape of exotic phenomena ranging from unconventional superconductivity to the emergence of topological phases. In this study, we have employed pulsed laser deposition technique to grow SrIrO 3 (SIO) thin films on SrTiO 3 substrates, systematically varying the oxygen content during the post-deposition annealing. X-ray photoelectron spectroscopy (XPS) provided insights into the stoichiometry and spin-orbit splitting energy of Iridium within the SIO film, while high-resolution x-ray studies meticulously examined the structural integrity of the thin films. Remarkably, our findings indicate a decrease in the metallicity of SIO thin films with reduced annealing O 2 partial pressure. Furthermore, we carried out magneto-transport studies on the SIO thin films, the results revealed intriguing insights into spin transport as a function of oxygen content. The tunability of the electronic band structure of SIO films with varying oxygen vacancy is correlated with the density functional theory calculations. Our findings elucidate the intricate mechanisms dictating spin transport properties in SIO thin films, offering invaluable guidance for the design and optimization of spintronic devices based on complex oxide materials. Notably, the ability to tune bandwidth by varying post-annealing oxygen partial pressure in iridate-based spintronic materials holds significant promise for advancing technological applications in the spintronics domain., (© 2024 IOP Publishing Ltd.)

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., (© 2024 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

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 O 2 ., (© 2023 Wiley Periodicals LLC.)

Here, we investigate the structural and electrical properties on SrIrO 3 films grown on LaAlO 3 (100) substrate with varying thickness (18, 25 and 40 nm). The x-ray diffraction shows good quality epitaxial films without any chemical impurity. The out-of-plane lattice parameter increases with the film thickness. All the films show a semiconducting behavior where the resistivity increases with increasing thickness. Our analysis shows at high temperature the charge conduction mechanism follows Mott's two-dimensional variable-range-hopping model. Detailed current-voltage ( I - V ) measurements show a linear Ohmic behavior at room temperature, however, a prominent deviation from linearity has been observed at low temperatures where the exponent n (I∝Vn) increases with decreasing temperature, reachingn ∼1.5 at low temperature. Analysis of I - V data indicates that the charge conduction has dominant contribution of Poole-Frenkel mechanism rather than Schottky behavior. This evolution of charge transport with temperature is quite intriguing which may be related to the development of low temperature magnetism in films of SrIrO 3 ., (© 2023 IOP Publishing Ltd.)

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.

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+ U on-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., (Creative Commons Attribution license.)

Understanding the various competing phases in cuprate superconductors is a long-standing challenging problem. Recent studies have shown that orbital degrees of freedom, both Cu e g orbitals and O p orbitals, 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 of p orbitals is critical to the charge-stripe features, which induce two types of stripe phases with s )-wave and d -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., (© 2023 IOP Publishing Ltd.)

Actinoid tetroxide molecules AnO 4 (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 4 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., (© 2022 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

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., (© 2022 The Authors. Advanced Materials published by Wiley-VCH GmbH.)

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 2 ), via doping with 3d transition metal vanadium is demonstrated. Single-crystalline vanadium-doped IrO 2 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., (© 2022 The Authors. Advanced Science published by Wiley-VCH GmbH.)

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., (© 2022 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Correlated states have emerged in low-dimensional systems owing to enhanced Coulomb interactions. Elucidating these states requires atomic-scale characterization and delicate control capabilities. Herein, spectroscopic imaging-scanning tunneling microscopy was employed to investigate the correlated states residing in 1D electrons of the monolayer and bilayer MoSe 2 mirror twin boundary (MTB). The Coulomb energies, determined by the wire length, drive the MTB into two types of ground states with distinct respective out-of-phase and in-phase charge orders. The two ground states can be reversibly converted through a metastable zero-energy state with in situ voltage pulses, which tune the electron filling of the MTB via a polaronic process, substantiated by first-principles calculations. Our Hubbard model calculation with an exact diagonalization method reveals the ground states as correlated insulators from an on-site U-originated Coulomb interaction, dubbed the Hubbard-type Coulomb blockade effect. Our study lays a foundation for understanding and tailoring correlated physics in complex systems., (© The Author(s) 2022. Published by Oxford University Press on behalf of China Science Publishing & Media Ltd.)

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.

The discrepancy among one-electron and two-electron densities for diverse N -electron atomss, enclosing neutral systems (with nuclear charge Z=N) and charge-one ions (|N-Z|=1), is quantified by means of mutual information, I , and Quantum Similarity Index, QSI, in the conjugate spaces position/momentum. These differences can be interpreted as a measure of the electron correlation of the system. The analysis is carried out by considering systems with a nuclear charge up to Z=103 and singly charged ions (cations and anions) as far as N=54. The interelectronic correlation, for any given system, is quantified through the comparison of its double-variable electron pair density and the product of the respective one-particle densities. An in-depth study along the Periodic Table reveals the importance, far beyond the weight of the systems considered, of their shell structure.

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., (© 2022 The Authors.)

The interplay of electron correlations and topological phases gives rise to various exotic phenomena including fractionalization, excitonic instability and axionic excitation. Recently discovered transition-metal pentatellurides can reach the ultra-quantum limit in low magnetic fields and serve as good candidates for achieving such a combination. Here, we report evidence of density wave and metal-insulator transition in HfTe 5 induced by intense magnetic fields. Using the non-linear transport technique, we detect a distinct non-linear conduction behavior in the longitudinal resistivity within the a-c plane, corresponding to the formation of a density wave induced by magnetic fields. In high fields, the onset of non-linear conduction in the Hall resistivity indicates an impurity-pinned magnetic freeze-out as the possible origin of the insulating behavior. These frozen electrons can be gradually reactivated into mobile states above a threshold of electric field. This experimental evidence calls for further investigation into the underlying mechanism of the bulk quantum Hall effect and field-induced phase transitions in pentatellurides., (© The Author(s) 2021. Published by Oxford University Press on behalf of China Science Publishing & Media Ltd.)

The ground-state structural transition in small lithium clusters Li n (n = 4 - 6) is analyzed based on the many-body expansion of the interaction energy using the total energy calculated by the fixed-node diffusion Monte Carlo (FN-DMC) simulations. The results show that the transition from 2D to 3D structure occurs through an intricate competition of attractive and repulsive interaction energies. As the structure dimensionality increases from 2D to 3D, the electron-correlation contribution to the interaction energy in the isomer of the ground-state structure is always the largest.

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., (© 2021 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

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., (© 2020 Wiley Periodicals LLC.)

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.

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, [UO 2 ] 2+ , [NpO 2 ] 2+ , [NpO 2 ] + and [PuO 2 ] 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 [AnO 2 ] 2+/+ -H 2 O 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 [AnO 2 ] +/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 [NpO 2 ] + -H 2 O 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, [NpO 2 ] 2+ and [PuO 2 ] 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.

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 (2018 Nature Communications 9 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 2 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 2 , 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., (© 2020 IOP Publishing Ltd.)

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 5 d electrons 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., (© 2020 IOP Publishing Ltd.)

The graphitic carbon nitride (g-C 3 N 4 ) 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-C 3 N 4 monolayers. The electron correlations between the 3 d electrons 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-C 3 N 4 is 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-C 3 N 4 is 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-C 3 N 4 , 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-C 3 N 4 so as to promote its applications., (© 2020 IOP Publishing Ltd.)

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., (© 2020 IOP Publishing Ltd.)

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 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.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.

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 2 (X = Cl or F) with HCN or NH 3 , while the remaining three combine ClF with HCN, NH 3 or N 2 . 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 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.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.

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 3 /SrRuO 3 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 3 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., (© 2019 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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., (© 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.)

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) 2 ] and dicyanomercury(II) [Hg(CN) 2 ]. 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) 2 and Hg(CN) 2 , 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 driven similarity renormalization group (DSRG) provides an alternative way to address the intruder state problem in quantum chemistry. In this review, we discuss recent developments of multireference methods based on the DSRG. We provide a pedagogical introduction to the DSRG and its various extensions and discuss its formal properties in great detail. In addition, we report several illustrative applications of the DSRG to molecular systems.

We elucidate the approaches used to incorporate electron correlation in existing semiempirical molecular orbital theory (SEMO) methods and compare them with the techniques used in other quantum chemical methods. After analyzing expressions for electron correlation in ab initio wavefunction theory, density functional theory, and density functional-based tight-binding (TB) methods, we suggest a framework for developing hybrid TB-SEMO methods. We provide a numerical proof-of-concept for such a method based on the OM2 method.

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., Competing Interests: The authors declare no conflict of interest.

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.

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 , 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 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., Competing Interests: The authors declare that they have no competing interests.

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., (© 2018 Wiley Periodicals, Inc.)

Undecacene, the largest member of the acene family, was photogenerated in a polymer matrix under cryogenic conditions from a precursor with two α-diketone bridges. The electronic absorption of undecacene extends into the NIR region, but the spectrum is dominated by two strong absorptions in the UV/Vis range. The HOMO-LUMO transition is continuously shifted to lower energies as the number of rings of the acene increases (E vs. 1/N, N=number of rings), which is in line with the strong electron correlation calculated by DFT/MRCI., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

When electronic correlation energy is partitioned topologically, a detailed picture of its distribution emerges, both within atoms and between any two atoms. This methodology allows one to study dispersion beyond its more narrow definition in long-range Rayleigh-Schrödinger perturbation theory. The interacting quantum atoms (IQA) method was applied to MP2/6-31G(d,p) (uncontracted) wave functions of a wide variety of systems: glycine…water (hydration), the ethene dimer (π-π interactions), benzene (aromaticity), cyclobutadiene (antiaromaticity), and NH 3 BH 3 (dative bond). Through the study of molecular complexes it turns out that dispersion energy is either important to a system's stabilization (for the C 2 H 4 dimer) or not important (for Gly…H 2 O). We have also discovered that the delocalization in benzene lowers the strength of Coulomb repulsion in the bonds, which has been quantified for the first time through IQA. Finally, we showed that the nature of the dative bond is much different from that of a regular covalent bond as it is not destabilized by electronic correlation. Finally, the conclusions obtained for these archetypical systems have implications for the future of the quantum topological force field FFLUX in the simulation of larger systems. Graphical abstract Atomic and bond dynamic correlation energies are now available thanks to IQA. Larger molecules can now be accessed to include resonance and solvation of FFLUX force field.

The description of the electronic structure of molecules in terms of molecular orbitals is a highly successful concept in chemistry. However, it commonly fails if the electrons in a molecule are strongly correlated and cannot be treated as independent particles. Electron correlation is essential to understand inner-valence X-ray spectroscopies, it can drive ultrafast charge migration in molecules, and it is responsible for many exotic properties of strongly correlated materials. Time-resolved spectroscopy with attosecond resolution is generally capable of following electronic motion in real time and can thus provide experimental access to electron-correlation-driven phenomena. High-harmonic spectroscopy in particular uses the precisely timed laser-driven recollision of electrons to interrogate the electronic structure and dynamics of the investigated system on a sub-femtosecond timescale. In this Review, the capabilities of high-harmonic spectroscopy to follow electronic motion in molecules are discussed. Both qualitative and quantitative approaches to unraveling the detailed dynamical responses of molecular systems following ionization are presented. A new theoretical formalism for the reconstruction of correlation-driven charge migration is introduced. The importance of electron-ion entanglement and electronic coherence in the reconstruction of attosecond hole dynamics are discussed. These advances make high-harmonic spectroscopy a promising technique to decode fundamental electron correlations and to provide experimental data on the complex manifestations of multi-electron dynamics., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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 3 , in the B3LYP functional is analyzed for archetype ionic and covalent systems as well as for the Fe 2+ 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 3 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 3 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 ) 0 (dz2 ) 2 and (M-Lx2-y2 ) 0 (dx2-y2 ) 2 configurations into the dominant (M-Lx2-y2 ) 2 (dx2-y2 ) 0 and (M-Lx2-y2 ) 2 (dx2-y2 ) 0 ones in the low(er) spin states is put forward. That is, in the effort to mimic such electron-electron interactions, E x 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.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., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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'., (© 2018 The Author(s).)

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., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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 (F 16 CoPc/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.

The geometric and electronic ground-state structures of 30 isomers of six MS 4 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 4 species were compared to analogous MO 4 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 4 and (Ru,Os,Hs)O 4 , whereas a low MOS of two appeared in the high-spin septet D 2d species Fe(S 2 ) 2 and (slightly excited) metastable Fe(O 2 ) 2 . The ground states of all other molecules had intermediate MOS values, with S 2- , S 2 2- , S 2 1- (and O 2- , O 1- , O 2 2- , O 2 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 4 and MS 4 species provides insight into the periodicity of oxidation states and bonding., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Random-phase approximation (RPA) methods are rapidly emerging as cost-effective validation tools for semilocal density functional computations. We present the theoretical background of RPA in an intuitive rather than formal fashion, focusing on the physical picture of screening and simple diagrammatic analysis. A new decomposition of the RPA correlation energy into plasmonic modes leads to an appealing visualization of electron correlation in terms of charge density fluctuations. Recent developments in the areas of beyond-RPA methods, RPA correlation potentials, and efficient algorithms for RPA energy and property calculations are reviewed. The ability of RPA to approximately capture static correlation in molecules is quantified by an analysis of RPA natural occupation numbers. We illustrate the use of RPA methods in applications to small-gap systems such as open-shell d- and f-element compounds, radicals, and weakly bound complexes, where semilocal density functional results exhibit strong functional dependence.

The X-ray constrained wavefunction (XC-WF) method proposed by Jayatilaka [Jayatilaka & Grimwood (2001) ▸, Acta Cryst. A 57 , 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.

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'., (© 2016 The Author(s).)

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., (© 2015 Wiley Periodicals, Inc.)

Although the organic superconductors κ-(BEDT-TTF)2X (X=Cu[N(CN)2]Br and Cu(NCS)2) have been studied by NMR spectroscopy, hyperfine coupling tensors are required to quantify NMR spectra. Angle dependences of NMR spectra were measured to determine hyperfine coupling tensors applicable to further NMR assessments of attractive physical phenomena on κ-salts. The tensors of κ-(BEDT-TTF)2X and β(')-(BEDT-TTF)2ICl2 salts were compared to determine the hyperfine coupling mechanism in organic metals, with the results indicating that off-site dimer contribution should be considered. We also report the electron correlation of these salts and further application of the tensors., (Copyright © 2015 Elsevier Inc. All rights reserved.)

We propose a set of design rules with a model Hamiltonian that allows electrons to form attracting pairs through the exploitation of a new combination of resonant band alignment and Coulombic repulsion. The pair bands and single particle bands in various lattices are calculated and compared in energy, and regions of net attraction are identified. This work provides guidelines for the construction of molecular systems, nanocrystals, and nanoparticle arrays with the potential for superconductivity.