
We propose that the noncentrosymmetric LiGaGetype hexagonal $ABC$ crystal
SrHgPb realizes a new type of topological semimetal that hosts both Dirac and
Weyl points in momentum space. The symmetryprotected Dirac points arise due to
a band inversion and are located on the sixfold rotation $z$axis, whereas the
six pairs of Weyl points related by sixfold symmetry are located on the
perpendicular $k_z=0$ plane. By studying the electronic structure as a function
of the buckling of the HgPb layer, which is the origin of inversion symmetry
breaking, we establish that the coexistence of Dirac and Weyl fermions defines
a phase separating two topologically distinct Dirac semimetals. These two Dirac
semimetals are distinguished by the $\mathbb{Z}_2$ index of the $k_z=0$ plane
and the corresponding presence or absence of 2D Dirac fermions on side
surfaces. We formalize our firstprinciples calculations by deriving and
studying a lowenergy model Hamiltonian describing the DiracWeyl semimetal
phase. We conclude by proposing several other materials in the
noncentrosymmetric $ABC $ material class, in particular SrHgSn and CaHgSn, as
candidates for realizing the DiracWeyl semimetal.

Weyl semimetals are gapless topological states of matter with broken
inversion and/or time reversal symmetry, which can support unconventional
responses to externally applied electrical, optical and magnetic fields. Here
we report a new photogalvanic effect in typeII WSMs, MoTe2 and Mo0.9W0.1Te2,
which are observed to support a circulating photocurrent when illuminated by
circularly polarized light at normal incidence. This effect occurs exclusively
in the inversion broken phase, where crucially we find that it is associated
with a spatially varying beam profile via a new dispersive contribution to the
circular photogalvanic effect (sCPGE). The response functions derived for
sCPGE reveal the microscopic mechanism of this photocurrent, which are
controlled by terms that are allowed in the absence of inversion symmetry,
along with asymmetric carrier excitation and relaxation. By evaluating this
response for a minimal model of a Weyl semimetal, we obtain the frequency
dependent scaling behavior of this form of photocurrent. These results
demonstrate opportunities for controlling photoresponse by patterning optical
fields to store, manipulate and transmit information over a wide spectral
range.

We have investigated the effect of transition metal dopants on the local
structure of the prototypical 0.75 Pb(Mg$_{1/3}$Nb$_{2/3}$)O$_3$0.25 PbTiO$_3$
relaxor ferroelectric. We find that these dopants give rise to very different
local structure and other physical properties. For example, when Mg is
partially substituted by Cu or Zn, the displacement of Cu or Zn is much larger
than that of Mg, and is even comparable to that of Nb. The polarization of
these systems is also increased, especially for the Cudoped solution, due to
the large polarizability of Cu and Zn. As a result, the predicted maximum
dielectric constant temperatures ($T_m$) are increased. On the other hand, the
replacement of a Ti atom with a Mo or Tc dramatically decreases the
displacements of the cations and the polarization, and thus, the $T_m$ values
are also substantially decreased. The higher $T_m$ cannot be explained by the
conventional argument based on the ionic radii of the cations. Furthermore, we
find that Cu, Mo, or Tc doping increase the cations displacement disorder. The
effect of the dopants on the temperature dispersion ${\Delta}T_m$, which is the
change of $T_m$ for different frequencies, is also discussed. Our findings lay
the foundation for further investigations of unexplored dopants.

Controlled chemical transformation of nanostructures is a promising technique
to obtain precisely designed novel materials which are difficult to synthesize
otherwise. We report hightemperature vapor phase anion exchange reactions to
chemically transform IIVI semiconductor nanostructures (100300 nm length
scale) while retaining the single crystallinity, crystal structure, morphology
and even defect distribution of the parent material via atomic templating. The
concept of atomic templating is employed to obtain kinetically controlled,
thermodynamically metastable structural phases such as zincblende CdSe and CdS
from zincblende CdTe upon complete chemical replacement of Te with Se or S. The
underlying transformation mechanisms are explained through firstprinciples
density functional theory calculations. Atomic templating is a unique path to
independently tune materials phase and composition at the nanoscale allowing
synthesis of novel materials.

The response functions of a material characterize its behavior under external
stimuli, such as electromagnetic radiation. Such responses may grow linearly
with the amplitude of the incident radiation, as is the case of absorption, or
may be nonlinear. The latter category includes a diverse set of phenomena such
as second harmonic generation (SHG), shift current, sum frequency generation,
and excited state absorption, among others. Despite decades of research into
nonlinear response theory, and the occasional discovery of materials with large
nonlinear responses, there has been no systematic investigation into the
maximum amount of nonlinear optical response attainable in solidstate
materials. In this work, we present an upper bound on the secondorder response
functions of materials, which controls the SHG and shift current responses. We
show that this bound depends on the band gap, band width, and geometrical
properties of the material in question. We find that Kuzyk's bound for the
maximum SHG of isolated molecules can be exceeded by conjugation or
condensation of molecules to form molecular solids, and that strongly coupled
systems generally have larger responses than weakly coupled or isolated ones.
As a proof of principle, we perform firstprinciples calculations of the
response tensors of a wide variety of materials, finding that the materials in
our database do not yet saturate the upper bound. This suggests that new large
SHG and shift current materials will likely be discovered by future materials
research guided by the factors mentioned in this work.

We show that the carrier recombination rate of noncentrosymmetric materials
can be strongly modified by spinorbit coupling. Our proposed mechanism
involves the separation of conduction and valence bands into their respective
spin components, which changes the transition dipole moments between them. The
change in the carrier recombination can be either positive or negative in sign,
or vary depending on the location of carriers in the Brillouin zone. We have
performed a large scale DFT screening study to identify candidate materials
that display this effect. We have selected three materials, Pb$_4$SeBr$_6$,
ReTe$_3$Br$_5$, and CsCu(BiS$_2$)$_2$, which span the range of behaviors, and
discuss their electronic band structure in greater detail. We find transition
dipole moment enhancement factors of up to three orders of magnitude,
reflecting the physical impact of spinorbit coupling on the carrier lifetime.
Therefore, further explorations of the spinorbit coupling and lattice symmetry
could prove to be useful for manipulating the photophysics of materials.

The consistency between the exchangecorrelation functional used in
pseudopotential construction and in the actual density functional theory
calculation is essential for the accurate prediction of fundamental properties
of materials. However, routine hybrid density functional calculations at
present still rely on GGA pseudopotentials due to the lack of hybrid functional
pseudopotentials. Here, we present a scheme for generating hybrid functional
pseudopotentials, and we analyze the importance of pseudopotential density
functional consistency for hybrid functionals. We benchmark our PBE0
pseudopotentials for structural parameters and fundamental electronic gaps of
the G2 molecular dataset and some simple solids. Our results show that using
our new PBE0 pseudopotentials in PBE0 calculations improves agreement with
respect to allelectron calculations.

Tribopolymer formed on the contacts of microelectromechanical and
nanoelectromechanical system (MEMS/NEMS) devices is a major concern hampering
their practical use in information technology. Conductive metal oxides, such as
RuO$_2$ and ReO$_3$, have been regarded as promising candidate materials for
MEMS/NEMS contacts due to their conductivity, hardness, and relatively
chemically inert surfaces. However, recent experimental works demonstrate that
trace amounts of polymer could still form on RuO$_2$ surfaces. In this study,
we demonstrate the mechanism of this class of unexpected tribopolymer formation
by conducting density functional theory based computational compression
experiments with benzene as the contamination gas. First, mechanical force
during compression changes the benzene molecules from slightly physisorbed to
strongly chemisorbed. Further compression causes deformation and chemical
linkage of the benzene molecules. Finally, the two contacts detach, with one
having a complex organic molecule attached and the other with a more reactive
surface. The complex organic molecule, which has an oxabicyclic segment, can be
viewed as the rudiment of tribopolymer, and the more reactive surface can
trigger the next adsorptionreactiontribopolymer formation cycle. Based on
these results, we also predict tribopolymer formation rates by using
transitionstate theory and the secondorder rate law. This study deepens our
understanding of tribopolymer formation (especially on metal oxides) and
provides strategies for suppressing tribopolymerization.

Hybrid leadhalide perovskites have emerged as an excellent class of
photovoltaic materials. Recent reports suggest that the organic molecular
cation is responsible for local polar fluctuations that inhibit carrier
recombination. We combine low frequency Raman scattering with firstprinciples
molecular dynamics (MD) to study the fundamental nature of these local polar
fluctuations. Our observations of a strong central peak in both hybrid
(CH$_3$NH$_3$PbBr$_3$) and allinorganic (CsPbBr$_3$) leadhalide perovskites
show that anharmonic, local polar fluctuations are intrinsic to the general
leadhalide perovskite structure, and not unique to the dipolar organic cation.
MD simulations show that headtohead Cs motion coupled to Br face expansion,
on a few hundred femtosecond time scale, drives the local polar fluctuations in
CsPbBr$_3$.

Thermal expansion in materials can be accurately modeled with careful phonon
calculations within density functional theory. However, because of interest in
controlling thermal expansion and the time consumed evaluating thermal
expansion properties of candidate materials, either theoretically or
experimentally, an approach to rapidly identifying materials with desirable
thermal expansion properties would be of great utility. We show that the
fraction of crystal volume occupied by ions, based upon ionic radii, the
deviation of bond coordination divided by mean bond coordination, and the ratio
of least to greatest atomic mass are descriptors that correlate with the
coefficient of thermal expansion for a variety of materials found in widely
accessible databases. Correlation is greatly improved by combining these
descriptors in a multidimensional quadratic fit. Open space combined with a
range of atomic masses and a variety of local bond coordinations indicates
materials with lower coefficients of thermal expansion. Materials with
singlevalued local coordination and less open space have the highest
coefficients of thermal expansion.

There is a surge of interest in developing environmentally friendly cooling
technology based on the solidstate electrocaloric effect (ECE). Here, we
point out that negative ECE with a fast cooling rate ($\approx$10$^{11}$ K/s)
can be achieved by driving solid crystals to a hightemperature phase with a
properly designed electric field pulse. Specifically, we predict that an
ultrafast electric field pulse can cause a negative ECE up to 35 K in PbTiO$_3$
occurring on few picosecond time scales. We acquire and analyze these results
by clarifying the mechanism of ECE during an adiabatic irreversible process; In
addition to the conventional explanation of the ECE with entropy reallocation,
we simply portray the ECE with the concept of internal energy redistribution.
Electric field does work on a ferroelectric crystal and redistributes its
internal energy. How the kinetic energy is redistributed determines the
temperature change and strongly depends on the electric field temporal profile.
This concept is supported by our allatom molecular dynamics simulations of
PbTiO$_3$ and BaTiO$_3$. Based on this improved understanding of ECE, we
propose strategies for inducing both giant and negative ECE. This work offers a
more general framework to understand the ECE and highlights the opportunities
of electricfield engineering for controlled design of fast and efficient
cooling technology.

Hydrocarbon tribopolymer, a type of polymer formed due to friction between
surfaces, is a major impediment to the development of micro and
nanoelectromechanical systems (MEMS/NEMS) devices for industrial application.
Tribopolymer buildup can prevent MEMS and NEMS from making or breaking
electrical contact. We describe the adsorption of benzene (C6H6) on the
RuO2(110) surface using density functional theory. This adsorption is an
important initial step in the mechanism of hydrocarbon tribopolymer layer
formation on MEMS and NEMS devices. The adsorption interaction is studied by
considering three oxygen coverages of RuO2(110) and all the possible adsorption
sites for benzene. We find that adsorption of benzene on Opoor RuO2(110) via
CRu bonds is stronger than adsorption on the Orich RuO2(110) via HO bonds.
For an indepth study of the adsorption behavior, we include the van der Waals
interaction for a holistic investigation. By incorporating the thermodynamic
chemical potentials into the adsorption simulations, we describe a model that
can provide guidance for realistic situations.

We investigate the shift current bulk photovoltaic response of materials
close to a band inversion topological phase transition. We find that the bulk
photocurrent reverses direction across the band inversion transition, and that
its magnitude is enhanced in the vicinity of the phase transition. These
results are demonstrated with first principles DFT calculations of BiTeI and
CsPbI$_3$ under hydrostatic pressure, and explained with an analytical model,
suggesting that this phenomenon remains robust across disparate material
systems.

We report density functional theory (DFT) investigation of $B$site doped
CaFeO$_3$, a prototypical chargeordered perovskite. At 290 K, CaFeO$_3$
undergoes a metalinsulator transition and a charge disproportionation reaction
2Fe$^{4+}$$\rightarrow$Fe$^{5+}$+Fe$^{3+}$. We observe that when Zr dopants
occupy a (001) layer, the band gap of the resulting solid solution increases to
0.93 eV due to a 2D JahnTeller type distortion, where FeO$_6$ cages on the
$xy$ plane elongate along $x$ and $y$ alternatively between neighboring Fe
sites. Furthermore, we show that the rocksalt ordering of the Fe$^{5+}$ and
Fe$^{3+}$ cations can be enhanced when the $B$site dopants are arranged in a
(111) plane due to a collective steric effect that facilitates the size
discrepancy between the Fe$^{5+}$O$_6$ and Fe$^{3+}$O$_6$ octahedra and
therefore gives rise to a larger band gap. The enhanced charge
disproportionation in these solid solutions is verified by rigorously
calculating the oxidation states of the Fe cations with different octahedral
cage sizes. We therefore predict that the corresponding transition temperature
will increase due to the enhanced charge ordering and larger band gap. The
compositional, structural and electrical relationships exploited in this paper
can be extended to a variety of perovskites and nonperovskite oxides providing
guidance in structurally manipulating electrical properties of functional
materials.

The dielectric function is one of the most important quantities that
describes the electrical and optical properties of solids. Accurate modeling of
the frequencydependent dielectric function has great significance in the study
of the longrange van der Waals (vdW) interaction for solids and adsorption. In
this work, we calculate the frequencydependent dielectric functions of
semiconductors and insulators using the $GW$ method with and without exciton
effects, as well as efficient semilocal density functional theory (DFT), and
compare these calculations with a model frequencydependent dielectric
function. We find that for semiconductors with moderate band gaps, the model
dielectric functions, $GW$ values, and DFT calculations all agree well with
each other. However, for insulators with strong exciton effects, the model
dielectric functions have a better agreement with accurate $GW$ values than the
DFT calculations, particularly in highfrequency region. To understand this, we
repeat the DFT calculations with scissors correction, by shifting DFT KohnSham
energy gap to match the experimental band gap. We find that scissors correction
only moderately improves the DFT dielectric function in lowfrequency region.
Based on the dielectric functions calculated with different methods, we make a
comparative study by applying these dielectric functions to calculate the vdW
coefficients ($C_3$ and $C_5$) for adsorption of raregas atoms on a variety of
surfaces. We find that the vdW coefficients obtained with the nearlyfree
electron gasbased model dielectric function agree quite well with those
obtained from the $GW$ dielectric function, in particular for adsorption on
semiconductors, leading to an overall error of less than 7% for $C_3$ and 5%
for $C_5$. This demonstrates the reliability of the model dielectric function
for the study of physisorption.

We study negative thermal expansion (NTE) in model lattices with multiple
atoms per cell and first and secondnearest neighbor interactions using the
(anharmonic) Morse potential. By exploring the phase space of neighbor
distances and thermal expansion rates of the bonds, we determine the conditions
under which NTE emerges. By permitting all bond lengths to expand at different
rates, we find that NTE is possible without appealing to fully rigid units.
Nearly constant, largeamplitude, isotropic NTE is observed up to the melting
temperature in a classical molecular dynamics model of a $\mathrm{ReO}_3$like
structure when the rigidity of octahedral units is almost completely
eliminated. Only weak NTE, changing over to positive expansion is observed when
the cornerlinked octahedra are rigid, with flexible secondneighbor bonds
between neighboring octahedra permitting easy rotation. We observe similar
changes to thermal expansion behavior for the diamond lattice: NTE when
secondneighbor interactions are weak to positive thermal expansion when
secondneighbor interactions are strong. From these observations, we suggest
that the only essential local conditions for NTE are atoms with low
coordination numbers along with very low energies for changing bond angles
relative to bondstretching energies.

The noncollinear cycloidal magnetic order breaks the inversion symmetry in
CaMn$_{7}$O$_{12}$, generating one of the largest spinorbit driven
ferroelectric polarizations measured to date. In this Letter, the microscopic
origin of the polarization, including its direction, charge density
redistribution, magnetic exchange interactions, and its coupling to the spin
helicity, is explored via first principles calculations. The Berry phase
computed polarization exhibits almost pure electronic behavior, as the Mn
displacements are negligible, $\approx$~0.7~m\textrm{\AA}. The polarization
magnitude and direction are both determined by the Mn spin current, where the
\emph{p}\emph{d} orbital mixing is driven by the inequivalent exchange
interactions within the \emph{B}site Mn cycloidal spiral chains along each
Cartesian direction. We employ the generalized spincurrent model with
Heisenbergexchange DzyaloshinskiiMoriya interaction energetics to provide
insight into the underlying physics of this spindriven polarization.
Persistent electronic polarization induced by helical spin order in nearly
inversionsymmetric ionic crystal lattices suggests opportunities for ultrafast
magnetoelectric response.

Barium titanate (BaTiO$_3$) is a prototypical ferroelectric perovskite that
undergoes the rhombohedralorthorhombictetragonalcubic phase transitions as
the temperature increases. In this work, we develop a classical interatomic
potential for BaTiO$_3$ within the framework of the bondvalence theory. The
force field is parameterized from firstprinciples results, enabling accurate
largescale molecular dynamics (MD) simulations at finite temperatures. Our
model potential for BaTiO$_3$ reproduces the temperaturedriven phase
transitions in isobaricisothermal ensemble (NPT) MD simulations. This
potential allows the analysis of BaTiO$_3$ structures with atomic resolution.
By analyzing the local displacements of Ti atoms, we demonstrate that the phase
transitions of BaTiO$_3$ exhibit a mix of orderdisorder and displacive
characters. Besides, from detailed observation of structural dynamics during
phase transition, we discover that the global phase transition is associated
with changes in the equilibrium value and fluctuations of each polarization
component, including the ones already averaging to zero, Contrary to the
conventional understanding that temperature increase generally causes
bondsoftening transition, the x polarization component exhibits a
bondhardening character during the orthorhombic to tetragonal transition.
These results provide further insights about the temperaturedriven phase
transitions in BaTiO$_3$.

Piezoelectrics play a critical role in various applications. The permanent
dipole associated with the molecular cations in organometalhalide perovskites
(OMHPs) may lead to spontaneous polarization and thus piezoelectricity. Here,
we explore the piezoelectric properties of OMHPs with density functional
theory. We find that the piezoelectric coefficient depends sensitively on the
molecular ordering, and that the experimentally observed lightenhanced
piezoelectricity is due to to a nonpolar to polar structural transition. By
comparing OMHPs with different atomic substitutions in the $ABX_3$
architecture, we find that the displacement of the $B$site cation contributes
to nearly all the piezoelectric response, and that the competition between
$A$$X$ hydrogen bond and $B$$X$ metalhalide bond in OMHPs controls the
piezoelectric properties. These results highlight the potential of the OMHP
architecture for designing new functional photoferroelectrics and
photopiezoelectrics.

Microelectromechanical system and nanoelectromechanical system (MEMS and
NEMS) transistors are considered promising for sizereducing and
powermaximizing electronic devices. However, the tribopolymer which forms due
to the mechanical load to the surface contacts affects the conductivity between
the contacts dramatically. This is one of the challenging problems that prevent
widespread practical use of these otherwise promising devices. Here, we use
density functional theory (DFT) to investigate the mechanisms of tribopolymer
formation, including normal mechanical loading, the catalytic effect, as well
as the electrochemical effect of the metal contacts. We select benzene select
as the background gas, because it is one of the most common and severe
hydrocarbon contaminants. Two adsorption cases are considered: one is benzene
on the reactive metal surface, Pt(111), and the other is benzene on the noble
metal, Au(111). We demonstrate that the formation of tribopolymer is induced
both by the mechanical load and by the catalytic effect of the contact. First,
benzene molecules are adsorbed on the Pt surfaces. Then, due to the closure of
the Pt contacts, stress is applied to the adsorbates, making the CH bonds more
fragile. As the stress increases further, H atoms are pressed close to the Pt
substrate and begin to bond with Pt atoms. Thus Pt acts as a catalyst,
accelerating the dehydrogenation process. When there is voltage applied across
the contacts, the catalytic effect is enhanced by electrochemistry. Finally,
due to the loss of H atoms, C atoms become more reactive and link together or
pile up to form tribopolymer. By understanding these mechanisms, we provide
guidance on design strategies for suppressing tribopolymer formation.

We present an experimental and theoretical viewpoint on the electronic
carrier mobilities of typical hybrid organicinorganic perovskites (HOIPs).
While these mobilities are often quoted as high, a review of them shows that
although otherwise the semiconducting properties of HOIPs are impressively
good, mobilities of HOIPs used in most solar cells are actually not that high.
This is especially apparent if they are compared to those of inorganic
semiconductors used in other high efficiency solar cells. We critically examine
possible causes and focus on electronlattice coupling mechanisms that are
active at room temperature, and can lead to carrier scattering. From this, we
propose scattering due to acoustic phonons or polarons as possible causes, but
also point out the difficulties with each of these in view of additional
experimental and theoretical findings in the literature. Further research in
this direction will contribute to making HOIP solar cells even more efficient
than they already are.

Based on firstprinciples calculations, we predict a substantial increase in
the optical dielectric function of LiAsSe$_2$ under pressure. We find that the
optical dielectric constant is enhanced threefold under volume compression.
This enhancement is mainly due to the dimerization strength reduction of the
onedimensional (1D) AsSe chains in LiAsSe$_2$, which significantly alters
the wavefunction phase mismatch between two neighboring chains and changes the
transition intensity. By developing a tightbinding model of the interacting 1D
chains, the essential features of the lowenergy electronic structure of
LiAsSe$_2$ are captured. Our findings are important for understanding the
fundamental physics of LiAsSe$_2$ and provide a feasible way to enhance the
material optical response that can be applied to light harvesting for energy
applications.

The bulk photovoltaic effect is a longknown but poorly understood
phenomenon. Recently, however, the multiferroic bismuth ferrite has been
observed to produce strong photovoltaic response to visible light, suggesting
that the effect has been underexploited as well. Here we present three polar
oxides in the LiNbO$_3$ structure that we predict to have band gaps in the 12
eV range and very high bulk photovoltaic response: PbNiO$_3$,
Mg$_{1/2}$Zn$_{1/2}$PbO$_3$, and LiBiO$_3$. All three have band gaps determined
by cations with $d^{10}s^0$ electronic configurations, leading to conduction
bands composed of cation $s$orbitals and O $p$orbitals. This both
dramatically lowers the band gap and increases the bulk photovoltaic response
by as much as an order of magnitude over previous materials, demonstrating the
potential for highperforming bulk photovoltaics.

The simultaneous presence of seemingly incompatible properties of solids
often provides a unique opportunity to address questions of fundamental and
practical importance. The coexistence of ferroelectric and topological orders
is one such example. Ferroelectrics, which have a spontaneous macroscopic
polarization switchable by an applied electric field, usually are
semiconductors with a welldeveloped wide band gap with a few exceptions. On
the other hand, timereversal symmetric $Z_2$ topological insulators (TI),
characterized by robust metallic surface states protected by the topology of
the bulk, usually are narrowgap semiconductors ($< 0.7$ eV) which allow band
inversion induced by the spinorbit interaction. To date, a ferroelectric
topological insulator (FETI) has remained elusive, owing to the seemingly
contradictory characters of the ferroelectric and topological orders. Here, we
report that the FETI can be realized in halide perovskite CsPbI$_3$ under
strain. Our firstprinciples study reveals that a noncentrosymmetric
ferroelectric structure of CsPbI$_3$ is energetically favored under a wide
range of pressures, while maintaining its topological order. The proposed FETI
is characterized by switchable polar surfaces with spinmomentum locked Dirac
cones, which allows for electricfield control of topological surface states
(TSSs) and the surface spin current. Our demonstration of a FETI in a feasible
material opens doors for future studies combining ferroelectric and topological
orders, and offers a new paradigm for diverse applications in electronics,
spintronics, and quantum information.

The instability of organometal halide perovskites when in contact with water
is a serious challenge to their feasibility as solar cell materials. Although
studies of moisture exposure have been conducted, an atomistic understanding of
the degradation mechanism is required. Toward this goal, we study the
interaction of water with the (001) surfaces of CH$_3$NH$_3$PbI$_3$ under both
low and high water concentrations using density functional theory. We find that
water adsorption is heavily influenced by the orientation of the methylammonium
cations close to the surface. It is demonstrated that, depending on
methylammonium orientation, the water molecule can infiltrate into the hollow
site of the surface and get trapped. Controlling dipole orientation via poling
or interfacial engineering could thus enhance its moisture stability. We do not
see a direct reaction between the water and methylammonium molecules.
Furthermore, calculations with an implicit solvation model indicate that a
higher water concentration may facilitate degradation.