Ab Initio Electron-Phonon Interactions in Correlated Electron Systems.

Electron-phonon (e-ph) interactions are pervasive in condensed matter, governing phenomena such as transport, superconductivity, charge-density waves, polarons, and metal-insulator transitions. First-principles approaches enable accurate calculations of e-ph interactions in a wide range of solids. However, they remain an open challenge in correlated electron systems (CES), where density functional theory often fails to describe the ground state. Therefore reliable e-ph calculations remain out of reach for many transition metal oxides, high-temperature superconductors, Mott insulators, planetary materials, and multiferroics. Here we show first-principles calculations of e-ph interactions in CES, using the framework of Hubbard-corrected density functional theory (DFT+U) and its linear response extension (DFPT+U), which can describe the electronic structure and lattice dynamics of many CES. We showcase the accuracy of this approach for a prototypical Mott system, CoO, carrying out a detailed investigation of its e-ph interactions and electron spectral functions. While standard DFPT gives unphysically divergent and short-ranged e-ph interactions, DFPT+U is shown to remove the divergences and properly account for the long-range Fröhlich interaction, allowing us to model polaron effects in a Mott insulator. Our work establishes a broadly applicable and affordable approach for quantitative studies of e-ph interactions in CES, a novel theoretical tool to interpret experiments in this broad class of materials.

The Different Story of π Bonds.

We revisit "classical" issues in multiply bonded systems between main groups elements, namely the structural distortions that may occur at the multiple bonds and that lead, e.g., to trans-bent and bond-length alternated structures. The focus is on the role that orbital hybridization and electron correlation play in this context, here analyzed with the help of simple models for σ- and π-bonds, numerically exact solutions of Hubbard Hamiltonians and first principles (density functional theory) investigations of an extended set of systems.

First-principles Hubbard U approach for small molecule binding in metal-organic frameworks.

We apply first-principles approaches with Hubbard U corrections for calculation of small molecule binding energetics to open-shell transition metal atoms in metal-organic frameworks (MOFs). Using density functional theory with van der Waals dispersion-corrected functionals, we determine Hubbard U values ab initio through an established linear response procedure for M-MOF-74, for a number of different metal centers (M = Ti, V, Cr, Mn, Fe, Co, Ni, and Cu). While our ab initio U values differ from those used in previous work, we show that they result in lattice parameters and electronic contributions to CO2-MOF binding energies that lead to excellent agreement with experiments and previous results, yielding lattice parameters within 3%. In addition, U-dependent calculations for an example system, Co-MOF-74, suggest that the CO2 binding energy grows monotonically with the value of Hubbard U, with the binding energy shifting 4 kJ/mol (or 0.041 eV) over the range of U = 0-5.4 eV. These results provide insight into an approximate but computationally efficient means for calculation of small molecule binding energies to open-shell transition metal atoms in MOFs and suggest that the approach can be predictive with good accuracy, independent of the cations used and the availability of experimental data.

Donor and acceptor levels of organic photovoltaic compounds from first principles.

Accurate and efficient approaches to predict the optical properties of organic semiconducting compounds could accelerate the search for efficient organic photovoltaic materials. Nevertheless, predicting the optical properties of organic semiconductors has been plagued by the inaccuracy or computational cost of conventional first-principles calculations. In this work, we demonstrate that orbital-dependent density-functional theory based upon Koopmans' condition [Phys. Rev. B, 2010, 82, 115121] is apt for describing donor and acceptor levels for a wide variety of organic molecules, clusters, and oligomers within a few tenths of an electron-volt relative to experiment, which is comparable to the predictive performance of many-body perturbation theory methods at a fraction of the computational cost.

Role of electronic localization in the phosphorescence of iridium sensitizing dyes.

In this work we present a systematic study of three representative iridium dyes, namely, Ir(ppy)(3), FIrpic, and PQIr, which are commonly used as sensitizers in organic optoelectronic devices. We show that electronic correlations play a crucial role in determining the excited-state energies in these systems, due to localization of electrons on Ir d orbitals. Electronic localization is captured by employing hybrid functionals within time-dependent density-functional theory and with Hubbard-model corrections within the Δ-SCF approach. The performance of both methods are studied comparatively and shown to be in good agreement with experiment. The Hubbard-corrected functionals provide further insight into the localization of electrons and on the charge-transfer character of excited-states. The gained insight allows us to comment on envisioned functionalization strategies to improve the performance of these systems. Complementary discussions on the Δ-SCF method are also presented in order to fill some of the gaps in the literature.

Energy Level Alignment at the Cobalt Phosphate/Electrolyte Interface: Intrinsic Stability vs Interfacial Chemical Reactions in 5 V Lithium Ion Batteries.

The intrinsic stability of the 5 V LiCoPO4-LiCo2P3O10 thin-film (carbon-free) cathode material coated with MoO3 thin layer is studied using a comprehensive synchrotron electron spectroscopy in situ approach combined with first-principle calculations. The atomic-molecular level study demonstrates fully reversible electronic properties of the cathode after the first electrochemical cycle. The polyanionic oxide is not involved in chemical reactions with the fluoroethylene-containing liquid electrolyte even when charged to 5.1 V vs Li+/Li. The high stability of the cathode is explained on the basis of the developed energy level model. In contrast, the chemical composition of the cathode-electrolyte interface evolves continuously by involving MoO3 in the decomposition reaction with consequent leaching of oxide from the surface. The proposed mechanisms of chemical reactions are attributed to external electrolyte oxidation via charge transfer from the relevant electron level to the MoO3 valence band state and internal electrolyte oxidation via proton transfer to the solvents. This study provides a deeper insight into the development of both a doping strategy to enhance the electronic conductivity of high-voltage cathode materials and an efficient surface coating against unfavorable interfacial chemical reactions.

Spin-state crossover and hyperfine interactions of ferric iron in MgSiO(3) perovskite.

Using density functional theory plus Hubbard U calculations, we show that the ground state of (Mg,Fe)(Si,Fe)O(3) perovskite, the major mineral phase in Earth's lower mantle, has high-spin ferric iron (S=5/2) at both dodecahedral (A) and octahedral (B) sites. With increasing pressure, the B-site iron undergoes a spin-state crossover to the low-spin state (S=1/2) between 40 and 70 GPa, while the A-site iron remains in the high-spin state. This B-site spin-state crossover is accompanied by a noticeable volume reduction and an increase in quadrupole splitting, consistent with recent x-ray diffraction and Mössbauer spectroscopy measurements. The anomalous volume reduction leads to a significant softening in bulk modulus during the crossover, suggesting a possible source of seismic-velocity anomalies in the lower mantle.

Simulation of heme using DFT + U: a step toward accurate spin-state energetics.

We investigate the DFT + U approach as a viable solution to describe the low-lying states of ligated and unligated iron heme complexes. Besides their central role in organometallic chemistry, these compounds represent a paradigmatic case where LDA, GGA, and common hybrid functionals fail to reproduce the experimental magnetic splittings. In particular, the imidazole pentacoordinated heme is incorrectly described as a triplet by all usual DFT flavors. In this study, we show that a U parameter close to 4 eV leads to spin transitions and molecular geometries in quantitative agreement with experiments and that DFT + U represents an appealing tool in the description of iron porphyrin complexes, at a much reduced cost compared to correlated quantum-chemistry methods. The possibility of obtaining the U parameter from first principles is explored through a self-consistent linear-response formulation. We find that this approach, which proved to be successful in other iron systems, produces in this case some overestimation with respect to the optimal values of U.

Searching for high magnetization density in bulk Fe: the new metastable Fe6 phase.

We report the discovery of a new allotrope of iron by first principles calculations. This phase has Pmn2(1) symmetry, a six-atom unit cell (hence the name Fe6), and the highest magnetization density (Ms) among all the known crystalline phases of iron. Obtained from the structural optimizations of the Fe3C-cementite crystal upon carbon removal, Pmn2(1) Fe6 is shown to result from the stabilization of a ferromagnetic FCC phase, further strained along the Bain path. Although metastable from 0 to 50 GPa, the new phase is more stable at low pressures than the other well-known HCP and FCC allotropes and smoothly transforms into the FCC phase under compression. If stabilized to room temperature, for example, by interstitial impurities, Fe6 could become the basis material for high Ms rare-earth-free permament magnets and high-impact applications such as light-weight electric engine rotors or high-density recording media. The new phase could also be key to explaining the enigmatic high Ms of Fe16N2, which is currently attracting intense research activity.

Origin of magnetic interactions and their influence on the structural properties of Ni2MnGa and related compounds.

In this work, we perform first-principles DFT calculations to investigate the interplay between magnetic and structural properties in Ni(2)MnGa. We demonstrate that the relative stability of austenite (cubic) and non-modulated martensite (tetragonal) phases depends critically on the magnetic interactions between Mn atoms. While standard approximate DFT functionals stabilize the latter phase, a more accurate treatment of electronic localization and magnetism, obtained with DFT+U, suppresses the non-modulated tetragonal structure for the stoichiometric compound, in better agreement with experiments. We show that the Anderson impurity model, with Mn atoms treated as magnetic impurities, can explain this observation and that the fine balance between super-exchange RKKY type interactions mediated by Ni d and Ga p orbitals determines the equilibrium structure of the crystal. The Anderson model is also demonstrated to capture the effect of the number of valence electrons per unit cell on the structural properties, often used as an empirical parameter to tune the behavior of Ni(2)MnGa based alloys. Finally, we show that off-stoichiometric compositions with excess Mn promote transitions to a non-modulated tetragonal structure, in agreement with experiments.

Dispersible exfoliated zeolite nanosheets and their application as a selective membrane.

Thin zeolite films are attractive for a wide range of applications, including molecular sieve membranes, catalytic membrane reactors, permeation barriers, and low-dielectric-constant materials. Synthesis of thin zeolite films using high-aspect-ratio zeolite nanosheets is desirable because of the packing and processing advantages of the nanosheets over isotropic zeolite nanoparticles. Attempts to obtain a dispersed suspension of zeolite nanosheets via exfoliation of their lamellar precursors have been hampered because of their structure deterioration and morphological damage (fragmentation, curling, and aggregation). We demonstrated the synthesis and structure determination of highly crystalline nanosheets of zeolite frameworks MWW and MFI. The purity and morphological integrity of these nanosheets allow them to pack well on porous supports, facilitating the fabrication of molecular sieve membranes.

Extended DFT + U + V method with on-site and inter-site electronic interactions.

In this paper we introduce a generalization of the popular DFT + U method based on the extended Hubbard model that includes on-site and inter-site electronic interactions. The novel corrective Hamiltonian is designed to study systems for which electrons are not completely localized on atomic states (according to the general scheme of Mott localization) and hybridization between orbitals from different sites plays an important role. The application of the extended functional to archetypal Mott-charge-transfer (NiO) and covalently bonded insulators (Si and GaAs) demonstrates its accuracy and versatility and the possibility to obtain a unifying and equally accurate description for a broad range of very diverse systems.

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials.

QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

Realistic quantitative descriptions of electron transfer reactions: diabatic free-energy surfaces from first-principles molecular dynamics.

A general approach to calculate the diabatic surfaces for electron-transfer reactions is presented, based on first-principles molecular dynamics of the active centers and their surrounding medium. The excitation energy corresponding to the transfer of an electron at any given ionic configuration (the Marcus energy gap) is accurately assessed within ground-state density-functional theory via a novel penalty functional for oxidation-reduction reactions that appropriately acts on the electronic degrees of freedom alone. The self-interaction error intrinsic to common exchange-correlation functionals is also corrected by the same penalty functional. The diabatic free-energy surfaces are then constructed from umbrella sampling on large ensembles of configurations. As a paradigmatic case study, the self-exchange reaction between ferrous and ferric ions in water is studied in detail.

Density functional theory in transition-metal chemistry: a self-consistent Hubbard U approach.

Transition-metal centers are the active sites for a broad variety of biological and inorganic chemical reactions. Notwithstanding this central importance, density-functional theory calculations based on generalized-gradient approximations often fail to describe energetics, multiplet structures, reaction barriers, and geometries around the active sites. We suggest here an alternative approach, derived from the Hubbard U correction to solid-state problems, that provides an excellent agreement with correlated-electron quantum chemistry calculations in test cases that range from the ground state of Fe2 and Fe2- to the addition elimination of molecular hydrogen on FeO+. The Hubbard U is determined with a novel self-consistent procedure based on a linear-response approach.

A unified electrostatic and cavitation model for first-principles molecular dynamics in solution.

The electrostatic continuum solvent model developed by [Fattebert and Gygi J. Comput. Chem. 23, 662 (2002); Int. J. Quantum Chem. 93, 139 (2003)] is combined with a first-principles formulation of the cavitation energy based on a natural quantum-mechanical definition for the surface of a solute. Despite its simplicity, the cavitation contribution calculated by this approach is found to be in remarkable agreement with that obtained by more complex algorithms relying on a large set of parameters. Our model allows for very efficient Car-Parrinello simulations of finite or extended systems in solution and demonstrates a level of accuracy as good as that of established quantum-chemistry continuum solvent methods. We apply this approach to the study of tetracyanoethylene dimers in dichloromethane, providing valuable structural and dynamical insights on the dimerization phenomenon.

Electronic-enthalpy functional for finite systems under pressure.

We introduce the notion of electronic enthalpy for first-principles structural and dynamical calculations of finite systems under pressure. An external pressure field is allowed to act directly on the electronic structure of the system studied via the ground-state minimization of the functional E+PV(q), where V(q) is the quantum volume enclosed by a charge isosurface. The Hellmann-Feynman theorem applies, and assures that the ionic equations of motion follow an isoenthalpic dynamics. No pressurizing medium is explicitly required, while coatings of environmental ions or ligands can be introduced if chemically relevant. We apply this novel approach to the study of group-IV nanoparticles during a shock wave, highlighting the significant differences in the plastic or elastic response of the diamond cage under load, and their potential use as novel nanostructured impact-absorbing materials.