• We propose that the noncentrosymmetric LiGaGe-type hexagonal $ABC$ crystal SrHgPb realizes a new type of topological semimetal that hosts both Dirac and Weyl points in momentum space. The symmetry-protected 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 first-principles calculations by deriving and studying a low-energy model Hamiltonian describing the Dirac-Weyl semimetal phase. We conclude by proposing several other materials in the non-centrosymmetric $ABC $ material class, in particular SrHgSn and CaHgSn, as candidates for realizing the Dirac-Weyl 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 type-II 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 (s-CPGE). The response functions derived for s-CPGE 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 Cu-doped 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 high-temperature vapor phase anion exchange reactions to chemically transform II-VI semiconductor nanostructures (100-300 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 first-principles 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 solid-state materials. In this work, we present an upper bound on the second-order 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 first-principles 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 spin-orbit 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 spin-orbit coupling on the carrier lifetime. Therefore, further explorations of the spin-orbit coupling and lattice symmetry could prove to be useful for manipulating the photophysics of materials.
  • The consistency between the exchange-correlation 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 all-electron 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 adsorption--reaction--tribopolymer formation cycle. Based on these results, we also predict tribopolymer formation rates by using transition--state theory and the second--order rate law. This study deepens our understanding of tribopolymer formation (especially on metal oxides) and provides strategies for suppressing tribopolymerization.
  • Hybrid lead-halide 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 first-principles 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 all-inorganic (CsPbBr$_3$) lead-halide perovskites show that anharmonic, local polar fluctuations are intrinsic to the general lead-halide perovskite structure, and not unique to the dipolar organic cation. MD simulations show that head-to-head 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 single-valued 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 solid--state 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 high--temperature 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 all--atom 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 electric--field 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 O-poor RuO2(110) via C-Ru bonds is stronger than adsorption on the O-rich RuO2(110) via H-O bonds. For an in-depth 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 charge-ordered perovskite. At 290 K, CaFeO$_3$ undergoes a metal-insulator 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 Jahn-Teller 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 rock-salt 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 non-perovskite 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 frequency-dependent dielectric function has great significance in the study of the long-range van der Waals (vdW) interaction for solids and adsorption. In this work, we calculate the frequency-dependent 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 frequency-dependent 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 high-frequency region. To understand this, we repeat the DFT calculations with scissors correction, by shifting DFT Kohn-Sham energy gap to match the experimental band gap. We find that scissors correction only moderately improves the DFT dielectric function in low-frequency 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 rare-gas atoms on a variety of surfaces. We find that the vdW coefficients obtained with the nearly-free electron gas-based 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 second-nearest 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, large-amplitude, 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 corner-linked octahedra are rigid, with flexible second-neighbor bonds between neighboring octahedra permitting easy rotation. We observe similar changes to thermal expansion behavior for the diamond lattice: NTE when second-neighbor interactions are weak to positive thermal expansion when second-neighbor 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 bond-stretching energies.
  • The noncollinear cycloidal magnetic order breaks the inversion symmetry in CaMn$_{7}$O$_{12}$, generating one of the largest spin-orbit 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 spin-current model with Heisenberg-exchange Dzyaloshinskii-Moriya interaction energetics to provide insight into the underlying physics of this spin-driven polarization. Persistent electronic polarization induced by helical spin order in nearly inversion-symmetric ionic crystal lattices suggests opportunities for ultrafast magnetoelectric response.
  • Barium titanate (BaTiO$_3$) is a prototypical ferroelectric perovskite that undergoes the rhombohedral-orthorhombic-tetragonal-cubic phase transitions as the temperature increases. In this work, we develop a classical interatomic potential for BaTiO$_3$ within the framework of the bond-valence theory. The force field is parameterized from first-principles results, enabling accurate large-scale molecular dynamics (MD) simulations at finite temperatures. Our model potential for BaTiO$_3$ reproduces the temperature-driven phase transitions in isobaric-isothermal 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 order-disorder 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 bond-softening transition, the x polarization component exhibits a bond-hardening character during the orthorhombic to tetragonal transition. These results provide further insights about the temperature-driven phase transitions in BaTiO$_3$.
  • Piezoelectrics play a critical role in various applications. The permanent dipole associated with the molecular cations in organometal-halide 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 light-enhanced piezoelectricity is due to to a non-polar 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$ metal-halide bond in OMHPs controls the piezoelectric properties. These results highlight the potential of the OMHP architecture for designing new functional photoferroelectrics and photopiezoelectrics.
  • Micro-electro-mechanical system and nano-electro-mechanical system (MEMS and NEMS) transistors are considered promising for size-reducing and power-maximizing 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 C-H 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 organic-inorganic 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 electron-lattice 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 first-principles 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 one-dimensional (1D) As--Se chains in LiAsSe$_2$, which significantly alters the wavefunction phase mismatch between two neighboring chains and changes the transition intensity. By developing a tight-binding model of the interacting 1D chains, the essential features of the low-energy 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 long-known 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 1-2 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 high-performing 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 well-developed wide band gap with a few exceptions. On the other hand, time-reversal symmetric $Z_2$ topological insulators (TI), characterized by robust metallic surface states protected by the topology of the bulk, usually are narrow-gap semiconductors ($< 0.7$ eV) which allow band inversion induced by the spin-orbit 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 first-principles study reveals that a non-centrosymmetric 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 spin-momentum locked Dirac cones, which allows for electric-field 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.