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