• We perform ab initio molecular dynamics (AIMD) simulation of liquid water in the canonical ensemble at ambient conditions using the SCAN meta-GGA functional approximation, and carry out systematic comparisons with the results obtained from the GGA-level PBE functional, and Tkatchenko-Scheffler van der Waals (vdW) dispersion correction inclusive PBE functional. We analyze various properties of liquid water including radial distribution functions, oxygen-oxygen-oxygen triplet angular distribution, tetrahedrality, hydrogen bonds, diffusion coefficients, ring statistics, density of states, band gaps, and dipole moments. We find that the SCAN functional is generally more accurate than the other two functionals for liquid water by not only capturing the intermediate-range vdW interactions but also mitigating the overly strong hydrogen bonds prescribed in PBE simulations. We also compare the results of SCAN-based AIMD simulations in the canonical and isothermal-isobaric ensembles. Our results suggest that SCAN provides a reliable description for most structural, electronic, and dynamical properties in liquid water.
  • Proton transfer via hydronium and hydroxide ions in water is ubiquitous. It underlies acid-base chemistry, certain enzyme reactions, and even infection by the flu. Despite two-centuries of investigation, the mechanism underlying why hydronium diffuses faster than hydroxide in water is still not well understood. Herein, we employ state of the art Density Functional Theory based molecular dynamics, with corrections for nonlocal van der Waals interactions, and self-interaction in the electronic ground state, to model water and the hydrated water ions. At this level of theory, structural diffusion of hydronium preserves the previously recognized concerted behavior. However, by contrast, proton transfer via hydroxide is dominated by stepwise events, arising from a stabilized hyper-coordination solvation structure that discourages proton transfer. Specifically, the latter exhibits non-planar geometry, which agrees with neutron scattering results. Asymmetry in the temporal correlation of proton transfer enables hydronium to diffuse faster than hydroxide.
  • Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der Waals interactions that dictates the structure and dynamics of liquid water. Notably, SCAN captures the density difference between water and ice I{\it h} at ambient conditions, as well as many important structural, electronic, and dynamic properties of liquid water. These successful predictions of the versatile SCAN functional open the gates to study complex processes in aqueous phase chemistry and the interactions of water with other materials in an efficient, accurate, and predictive, ab initio manner.
  • Oxygen K-edge X-ray absorption spectra of liquid water are computed based on the configurations from advanced ab initio molecular dynamics simulations, as well as an electron excitation theory from the GW method. One one hand, the molecular structures of liquid water are accurately predicted by including both van der Waals interactions and hybrid functional (PBE0). On the other hand, the dynamic screening effects on electron excitation are approximately described by the recently developed enhanced static Coulomb hole and screened exchange approximation by Kang and Hybertsen [Phys. Rev. B 82, 195108 (2010)]. The resulting spectra of liquid water are in better quantitative agreement with the experimental spectra due to the softened hydrogen bonds and the slightly broadened spectra originating from the better screening model.
  • Originating from a broken spatial inversion symmetry, ferroelectricity is a functionality of materials with an electric dipole that can be switched by external electric fields. Spontaneous polarization is a crucial ferroelectric property, and its amplitude is determined by the strength of polar structural distortions. Density functional theory (DFT) is one of the most widely used theoretical methods to study ferroelectric properties, yet it is limited by the levels of approximations in electron exchange-correlation. On the one hand, the local density approximation (LDA) is considered to be more accurate for the conventional perovskite ferroelectrics such as BaTiO3 and PbTiO3 than the generalized gradient approximation (GGA), which suffers from the so-called super-tetragonality error. On the other hand, GGA is more suitable for hydrogen-bonded ferroelectrics than LDA, which largely overestimates the strength of hydrogen bonding in general. We show here that the recently developed general-purpose strongly constrained and appropriately normed (SCAN) meta-GGA functional significantly improves over the traditional LDA/GGA for structural, electric, and energetic properties of diversely-bonded ferroelectric materials with a comparable computational effort, and thus enhances largely the predictive power of DFT in studies of ferroelectric materials. We also address the observed system-dependent performances of LDA and GGA for ferroelectrics from a chemical bonding point of view.
  • By applying Wannier-based extended Kugel-Khomskii model, we carry out first-principles calculations and electronic structure analysis to understand the spin-phonon coupling effect in transition-metal perovskites. We demonstrate the successful application of our approach to SrMnO$_3$ and BiFeO$_3$. We show that both the electron orbitals under crystal field splitting and the electronic configuration should be taken into account in order to understand the large variances of spin-phonon coupling effects among various phonon modes as well as in different materials.
  • Elastic strain is potentially an important approach in tuning the properties of the improperly multiferroic hexagonal ferrites, the details of which have however been elusive due to the experimental difficulties. Employing the method of restrained thermal expansion, we have studied the effect of isothermal biaxial strain in the basal plane of h-LuFeO3 (001) films. The results indicate that a compressive biaxial strain significantly enhances the ferrodistortion, and the effect is larger at higher temperatures. The compressive biaxial strain and the enhanced ferrodistortion together, cause an increase in the electric polarization and a reduction in the canting of the weak ferromagnetic moments in h-LuFeO3, according to our first principle calculations. These findings are important for understanding the strain effect as well as the coupling between the lattice and the improper multiferroicity in h-LuFeO3. The experimental elucidation of the strain effect in h-LuFeO3 films also suggests that the restrained thermal expansion can be a viable method to unravel the strain effect in many other epitaxial thin film materials.
  • We report a strong thickness dependence of the complex frequency-dependent optical dielectric function in epitaxial CaMnO3(001) thin films on SrTiO3(001), LaAlO3(001), and SrLaAlO4(001) substrates. A doubling of the peak value of the imaginary part of the dielectric function and spectral shifts of 0.5 eV for a given magnitude of absorption are observed. On the basis of the experimental data and first-principles density functional theory calculations of the dielectric function, its evolution with thickness from 4 to 63 nm has several regimes. In the thinnest, strain-coherent films, the response is characterized by a significant contribution from the free surface that dominates strain effects. However, at intermediate and larger thicknesses approaching the bulk-like film, strain coherence and partial strain relaxation persist and in influence the dielectric function.
  • In ABO3 perovskites, oxygen octahedron rotations are common structural distortions that can promote large ferroelectricity in BiFeO3 with an R3c structure [1], but suppress ferroelectricity in CaTiO3 with a Pbnm symmetry [2]. For many CaTiO3-like perovskites, the BiFeO3 structure is a metastable phase. Here, we report the stabilization of the highly-polar BiFeO3-like phase of CaTiO3 in a BaTiO3/CaTiO3 superlattice grown on a SrTiO3 substrate. The stabilization is realized by a reconstruction of oxygen octahedron rotations at the interface from the pattern of nonpolar bulk CaTiO3 to a different pattern that is characteristic of a BiFeO3 phase. The reconstruction is interpreted through a combination of amplitude-contrast sub 0.1nm high-resolution transmission electron microscopy and first-principles theories of the structure, energetics, and polarization of the superlattice and its constituents. We further predict a number of new artificial ferroelectric materials demonstrating that nonpolar perovskites can be turned into ferroelectrics via this interface mechanism. Therefore, a large number of perovskites with the CaTiO3 structure type, which include many magnetic representatives, are now good candidates as novel highly-polar multiferroic materials [3].
  • Kohn-Sham density functional theory (DFT) is a widely-used electronic structure theory for materials as well as molecules. DFT is needed especially for large systems, ab initio molecular dynamics, and high-throughput searches for functional materials. DFT's accuracy and computational efficiency are limited by the approximation to its exchange-correlation energy. Currently, the local density approximation (LDA) and generalized gradient approximations (GGAs) dominate materials computation mainly due to their efficiency. We show here that the recently developed non-empirical strongly constrained and appropriately normed (SCAN) meta-GGA improves significantly over LDA and the standard Perdew-Burke-Ernzerhof GGA for geometries and energies of diversely-bonded materials (including covalent, metallic, ionic, hydrogen, and van der Waals bonds) at comparable efficiency. Thus SCAN may be useful even for soft matter. Often SCAN matches or improves upon the accuracy of a computationally expensive hybrid functional, at almost-GGA cost. SCAN is therefore expected to have a broad impact on materials science.
  • In this work, the solvation and electronic structure of the aqueous chloride ion solution was investigated using Density Functional Theory (DFT) based \textit{ab initio} molecular dynamics (AIMD). From an analysis of radial distribution functions, coordination numbers, and solvation structures, we found that exact exchange ($E_{\rm xx}$) and non-local van der Waals (vdW) interactions effectively \textit{weaken} the interactions between the Cl$^-$ ion and the first solvation shell. With a Cl-O coordination number in excellent agreement with experiment, we found that most configurations generated with vdW-inclusive hybrid DFT exhibit 6-fold coordinated distorted trigonal prism structures, which is indicative of a significantly disordered first solvation shell. By performing a series of band structure calculations on configurations generated from AIMD simulations with varying DFT potentials, we found that the solvated ion orbital energy levels (unlike the band structure of liquid water) strongly depend on the underlying molecular structures. In addition, these orbital energy levels were also significantly affected by the DFT functional employed for the electronic structure; as the fraction of $E_{\rm xx}$ was increased, the gap between the highest occupied molecular orbital of Cl$^-$ and the valence band maximum of liquid water steadily increased towards the experimental value.
  • Magnetoelastic coupling, i.e., the change of crystal lattice induced by a spin order, is not only scientifically interesting, but also technically important. In this work, we propose a general microscopic model from first-principles calculations to describe the magnetoelastic coupling and provide a way to construct the microscopic model from density functional theory calculations. Based on this model, we reveal that there exists a previously unexpected contribution to the electric polarization induced by the spin-order in multiferroics due to the combined effects of magnetoelastic coupling and piezoelectric effect. Interestingly and surprisingly, we find that this lattice deformation contribution to the polarization is even larger than that from the pure electronic and ion-displacement contributions in BiFeO3. This model of magnetoelastic coupling can be generally applied to investigate the other magnetoelastic phenomena.
  • To date, high-resolution electron microscopy has largely relied on using the phase of the exit wave function at the exit surface to form a high-resolution electron microscopic image. We have for the first time used chromatic aberration correction to implement a new imaging mode to achieve amplitude contrast imaging in high-resolution electron microscopy, allowing us to obtain directly interpretable high-resolution electron microscopic images with discrimination between light and heavy atomic columns. Using this imaging approach, we have successfully visualized the atomic structure in a BaTiO3/CaTiO3 superlattice with high spatial accuracy and discrimination between Ba and Ca columns, providing direct visualization of the Ca and Ba associated oxygen octahedral tilt that controls ferroelectric behavior in these superlattice structures. Furthermore, this approach offers new opportunities to unravel the structure in a wide range of materials, especially complex oxides with exotic behaviors based on specific structural arrangements.
  • Using combined theoretical and experimental approaches, we studied the structural and electronic origin of the magnetic structure in hexagonal LuFeO$_3$. Besides showing the strong exchange coupling that is consistent with the high magnetic ordering temperature, the previously observed spin reorientation transition is explained by the theoretically calculated magnetic phase diagram. The structural origin of this spin reorientation that is responsible for the appearance of spontaneous magnetization, is identified by theory and verified by x-ray diffraction and absorption experiments.
  • In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local vdW/dispersion interactions, via a fully self-consistent density-dependent dispersion correction, and approximate nuclear quantum effects (aNQE), via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx, vdW, and aNQE as resulting from a large-scale AIMD simulation of (H$_2$O)$_{128}$ at the PBE0+vdW level of theory, significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, $S_{\rm OO}(Q)$, and corresponding oxygen-oxygen radial distribution function, $g_{\rm OO}(r)$, that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment as demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of the PBE0 hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions.
  • The structural instabilities of ATcO$_3$(A= Ca, Sr, Ba) are investigated by first-principles calculations. In addiction to the large octahedral rotation instability, a weak ferroelectric tendency is identified in CaTcO$_3$. We show that the ferroelectricity in CaTcO$_3$ can be recovered by interface engineering based on CaTcO$_3$/BaTcO$_3$ superlattices, where the octahedral rotation is largely suppressed. The N{\'e}el temperature is found to be $\sim$ 816 K, indicating that CaTcO$_3$ can be engineered into a new room temperature multiferroic material.
  • The ionization potential distributions of hydrated hydroxide and hydronium are computed with many-body approach for electron excitations with configurations generated by {\it ab initio} molecular dynamics. The experimental features are well reproduced and found to be closely related to the molecular excitations. In the stable configurations, the ionization potential is mainly perturbed by water molecules within the first solvation shell. On the other hand, electron excitation is delocalized on both proton receiving and donating complex during proton transfer, which shifts the excitation energies and broadens the spectra for both hydrated ions.
  • Hexagonal LuFeO$_3$ films have been studied using x-ray absorption and optical spectroscopy. The crystal splittings of Fe$^{3+}$ are extracted as $E_{e'}-E_{e"}$=0.7 eV and $E_{a_1'}-E_{e'}$=0.9 eV and a 2.0 eV optical band gap is determined assuming a direct gap. First-principles calculations confirm the experiments that the relative energies of crystal field splitting states do follow $E_{a_1'}>E_{e'}>E_{e"}$ with slightly underestimated values and a band gap of 1.35 eV.
  • We calculate the x-ray absorption spectra of liquid water at ambient conditions and of hexagonal ice close to melting, using a static GW approach that includes approximately local field effects. Quantum dynamics of the nuclei is taken into account by averaging the absorption cross section over molecular configurations generated by path integral simulations. We find that inclusion of quantum disorder is essential to bring the calculated spectra in close agreement with experiment. In particular, the intensity of the pre-edge feature, a spectral signature of broken and distorted hydrogen bonds, is accurately reproduced, in water and ice, only when quantum nuclei are considered. The effect of the local fields is less important but non negligible, particularly in ice.
  • Based on the locality principle of insulating superlattices, we apply the method of Wu {\it et al} [Phys. Rev. Letter {\bf 101}, 087610 (2008)] to the piezoelectric strains of individual layers under fixed displacement field. For a superlattice of arbitrary stacking sequence an accurate model is acquired for predicting piezoelectricity. By applying the model in the superlattices where ferroelectric and antiferrodistortive modes are in competition, functional piezoelectricity can be achieved. A strong nonlinear effect is observed and can be further engineered in the PbTiO$_3$/SrTiO$_3$ superlattice and an interface enhancement of piezoelectricity is found in the BaTiO$_3$/CaTiO$_3$ superlattice.
  • The electric and magnetic properties of (BaTiO$_3$)$_n$/(CaMnO$_3$)$_n$ short-period superlattices are studied by the first-principles calculations. The local electric polarizations in the CaMnO$_3$ layers are significant, comparable to that in the BaTiO$_3$ layers. Remarkably, the electric polarization is almost doubled when the spin configuration changes from antiferromagnetic to ferromagnetic in the superlattices, indicating a giant magnetoelectric coupling. This enhancement of the magnetoelectric coupling is due to the suppression of the antiferrodistortive mode in the CaMnO$_3$ layers at the interfaces.
  • The spinel cobalt oxide Co3O4 is a magnetic semiconductor containing cobalt ions in Co2+ and Co3+ oxidation states. We have studied the electronic, magnetic and bonding properties of Co3O4 using density functional theory (DFT) at the Generalized Gradient Approximation (GGA), GGA+U, and PBE0 hybrid functional levels. The GGA correctly predicts Co3O4 to be a semiconductor, but severely underestimates the band gap. The GGA+U band gap (1.96 eV) agrees well with the available experimental value (~ 1.6 eV), whereas the band gap obtained using the PBE0 hybrid functional (3.42 eV) is strongly overestimated. All the employed exchange-correlation functionals predict 3 unpaired d electrons on the Co2+ ions, in agreement with crystal field theory, but the values of the magnetic moments given by GGA+U and PBE0 are in closer agreement with the experiment than the GGA value, indicating a better description of the cobalt localized d states. Bonding properties are studied by means of Maximally Localized Wannier Functions (MLWFs). We find d-type MLWFs on the cobalt ions, as well as Wannier functions with the character of sp3d bonds between cobalt and oxygen ions. Such hybridized bonding states indicate the presence of a small covalent component in the primarily ionic bonding mechanism of this compound.
  • We carry out first-principles calculations for CaTiO$_3$/BaTiO$_3$ superlattices with epitaxial strain corresponding to growth on a SrTiO$_3$ substrate, and consider octahedral rotations as well as ferroelectric distortions. The calculations are done as a function of electric displacement field, and both a macroscopic and a local electrostatic analysis are carried out. We find that strong octahedral rotations occur for TiO$_6$ octahedra sandwiched between CaO layers on both sides, but are strongly suppressed if either neighboring layer is a BaO layer. Due to the resulting enhancement of the ferroelectric instability in the BaO-neighboring octahedra, we find that overall the ferroelectric instability of the superlattice is enhanced by the interface. Thus, short-period superlattices in this system have a larger ferroelectric polarization than longer-period ones of the same average composition, contrary to the expected trend.
  • The x-ray absorption spectra of water and ice are calculated with a many-body approach for electron-hole excitations. The experimental features, including the small effects of temperature change in the liquid, are quantitatively reproduced from molecular configurations generated by ab-initio molecular dynamics. The spectral difference between the solid and the liquid is due to two major short range order effects. One, due to breaking of hydrogen bonds, enhances the pre-edge intensity in the liquid. The other, due to a non-bonded molecular fraction in the first coordination shell, affects the main spectral edge in the conversion of ice to water. This effect may not involve hydrogen bond breaking as shown by experiment in high-density amorphous ice.
  • Recent theoretical work has provided evidence that hybrid functionals, which include a fraction of exact (Hartree Fock) exchange in the density functional theory (DFT) exchange and correlation terms, significantly improve the description of band gaps of semiconductors compared with local and semilocal approximations. Based on a recently developed order-$N$ method for calculating the exact exchange in extended insulating systems, we have implemented an efficient scheme to determine the hybrid functional band gap. We use this scheme to study the band gap and other electronic properties of the ternary compound In$_{1-x}$Ga$_{x}$N using a 64-atom supercell model.