
We perform ab initio molecular dynamics (AIMD) simulation of liquid water in
the canonical ensemble at ambient conditions using the SCAN metaGGA functional
approximation, and carry out systematic comparisons with the results obtained
from the GGAlevel PBE functional, and TkatchenkoScheffler van der Waals (vdW)
dispersion correction inclusive PBE functional. We analyze various properties
of liquid water including radial distribution functions, oxygenoxygenoxygen
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
intermediaterange vdW interactions but also mitigating the overly strong
hydrogen bonds prescribed in PBE simulations. We also compare the results of
SCANbased AIMD simulations in the canonical and isothermalisobaric 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 acidbase chemistry, certain enzyme reactions, and even infection by
the flu. Despite twocenturies 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
selfinteraction 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 hypercoordination solvation structure that discourages proton
transfer. Specifically, the latter exhibits nonplanar 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 Kedge Xray 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 exchangecorrelation. 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 socalled
supertetragonality error. On the other hand, GGA is more suitable for
hydrogenbonded ferroelectrics than LDA, which largely overestimates the
strength of hydrogen bonding in general. We show here that the recently
developed generalpurpose strongly constrained and appropriately normed (SCAN)
metaGGA functional significantly improves over the traditional LDA/GGA for
structural, electric, and energetic properties of diverselybonded
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 systemdependent performances of LDA
and GGA for ferroelectrics from a chemical bonding point of view.

By applying Wannierbased extended KugelKhomskii model, we carry out
firstprinciples calculations and electronic structure analysis to understand
the spinphonon coupling effect in transitionmetal 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 spinphonon 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 hLuFeO3 (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 hLuFeO3, 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 hLuFeO3. The experimental
elucidation of the strain effect in hLuFeO3 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 frequencydependent
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 firstprinciples density functional theory calculations of the
dielectric function, its evolution with thickness from 4 to 63 nm has several
regimes. In the thinnest, straincoherent 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
bulklike 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 CaTiO3like perovskites, the BiFeO3 structure is a metastable
phase. Here, we report the stabilization of the highlypolar BiFeO3like 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 amplitudecontrast sub 0.1nm highresolution
transmission electron microscopy and firstprinciples 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 highlypolar multiferroic
materials [3].

KohnSham density functional theory (DFT) is a widelyused electronic
structure theory for materials as well as molecules. DFT is needed especially
for large systems, ab initio molecular dynamics, and highthroughput searches
for functional materials. DFT's accuracy and computational efficiency are
limited by the approximation to its exchangecorrelation 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 nonempirical strongly constrained and
appropriately normed (SCAN) metaGGA improves significantly over LDA and the
standard PerdewBurkeErnzerhof GGA for geometries and energies of
diverselybonded 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 almostGGA 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 nonlocal van der Waals (vdW)
interactions effectively \textit{weaken} the interactions between the Cl$^$
ion and the first solvation shell. With a ClO coordination number in excellent
agreement with experiment, we found that most configurations generated with
vdWinclusive hybrid DFT exhibit 6fold 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
firstprinciples 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 spinorder
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 iondisplacement contributions in BiFeO3. This model of
magnetoelastic coupling can be generally applied to investigate the other
magnetoelastic phenomena.

To date, highresolution electron microscopy has largely relied on using the
phase of the exit wave function at the exit surface to form a highresolution
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 highresolution electron microscopy, allowing us to obtain
directly interpretable highresolution 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 xray 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 exchangecorrelation (XC) functionals to investigate
the individual and collective effects of exact exchange (Exx), via the PBE0
hybrid functional, nonlocal vdW/dispersion interactions, via a fully
selfconsistent densitydependent 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 largescale 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 oxygenoxygen structure factor, $S_{\rm OO}(Q)$, and
corresponding oxygenoxygen 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 nonlocal vdW/dispersion
interactions.

The structural instabilities of ATcO$_3$(A= Ca, Sr, Ba) are investigated by
firstprinciples 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 manybody 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 xray 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. Firstprinciples 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 xray 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 preedge
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$
shortperiod superlattices are studied by the firstprinciples 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 exchangecorrelation 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 dtype 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 firstprinciples 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 BaOneighboring octahedra, we find that
overall the ferroelectric instability of the superlattice is enhanced by the
interface. Thus, shortperiod superlattices in this system have a larger
ferroelectric polarization than longerperiod ones of the same average
composition, contrary to the expected trend.

The xray absorption spectra of water and ice are calculated with a manybody
approach for electronhole excitations. The experimental features, including
the small effects of temperature change in the liquid, are quantitatively
reproduced from molecular configurations generated by abinitio 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 preedge intensity in the liquid. The other, due to a nonbonded
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 highdensity 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$_{1x}$Ga$_{x}$N using a 64atom supercell model.