BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrew M. Rappe

Andrew M. Rappe appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.mes-hallphysics.comp-phcond-mat.otherphysics.chem-phphysics.optics

Catalog footprint

What is connected

45works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2022arXiv

Recent Progress in the Theory of Bulk Photovoltaic Effect

The bulk photovoltaic effect (BPVE) occurs in solids with broken inversion symmetry and refers to DC current generation due to uniform illumination, without the need of heterostructures or interfaces, a feature that is distinct from the traditional photovoltaic effect. Its existence has been demonstrated almost 50 years ago, but predictive theories only appeared in the last ten years, allowing for the identification of different mechanisms and the determination of their relative importance in real materials. It is now generally accepted that there is an intrinsic mechanism that is insensitive to scattering, called shift current, where first-principles calculations can now give highly accurate predictions. Another important but more extrinsic mechanism, called ballistic current, is also attracting a lot of attention, but due to the complicated scattering processes, its numerical calculation for real materials is only made possible quite recently. In addition, an intrinsic ballistic current, usually referred to as injection current, will appear under circularly-polarized light and has wide application in experiments. In this article, experiments that are pertinent to the theory development are reviewed, and a significant portion is devoted to discussing the recent progress in the theories of BPVE and their numerical implementations. As a demonstration of the capability of the newly developed theories, a brief review of the materials design strategies enabled by the theory development is given. Finally, remaining questions in the BPVE field and possible future directions are discussed to inspire further investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

First Principles Calculation of Ballistic Current from Electron-Hole Interaction

The bulk photovoltaic effect (BPVE) has attracted an increasing interest due to its potential to overcome the efficiency limit of traditional photovoltaics, and much effort has been devoted to understanding its underlying physics. However, previous work has shown that theoretical models of the shift current and the phonon-assisted ballistic current in real materials do not fully account for the experimental BPVE photocurrent, and so other mechanisms should be investigated in order to obtain a complete picture of BPVE. In this Letter, we demonstrate two approaches that enable the ab initio calculation of the ballistic current originating from the electron-hole interaction in semiconductors. Using BaTiO$_3$ and MoS$_2$ as two examples, we show clearly that for them the asymmetric scattering from electron-hole interaction is less appreciable than that from electron-phonon interaction, indicating more scattering processes need to be included to further improve the BPVE theory. Moreover, our approaches build up a venue for predicting and designing materials with larger ballistic current due to electron-hole interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large Bulk Piezophotovoltaic Effect of Monolayer $2H$-MoS$_2$

The bulk photovoltaic effect in noncentrosymmetric materials is an intriguing physical phenomenon that holds potential for high-efficiency energy harvesting. Here, we study the shift current bulk photovoltaic effect in the transition metal dichalcogenide MoS$_2$. We present a simple automated method to guide materials design and use it to uncover a distortion to monolayer $2H$-MoS$_2$ that dramatically enhances the integrated shift current. Using this distortion, we show that overlap in the Brillouin zone of the distributions of the shift vector (a quantity measuring the net displacement in real space of coherent wave packets during excitation) and the transition intensity is crucial for increasing the shift current. The distortion pattern is related to the material polarization and can be realized through an applied electric field via the converse piezoelectric effect. This finding suggests an additional method to engineer the shift current response of materials to augment previously reported methods using mechanical strain.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The shift photovoltaic current and magnetically-induced bulk photocurrent in piezoelectric sillenite crystals

Recently, it has been shown how shift and ballistic currents in piezoelectric sillenite crystals Bi$_{12}$GeO$_{20}$ and Bi$_{12}$SiO$_{20}$ can be separated experimentally under the assumption that the shift component of the circular current is small. However, it has been claimed that the shift and ballistic currents cannot be quantified by this method, due to the magneto-photovoltaic effect caused either by the change of the crystal's spatial symmetry in the magnetic field or by the breaking of time-reversal symmetry. Presently, we report observations of photovoltaic currents in Bi$_{12}$SiO$_{20}$, excited by linearly- and circularly-polarized light under weak external magnetic field, as well as measurements of the corresponding photo-Hall signals. We demonstrate that the magneto-photovoltaic current constitutes a significant fraction of the measured current in the Hall direction for Bi$_{12}$SiO$_{20}$ under specific experimental conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The tetragonal phase of CH$_{3}$NH$_{3}$PbI$_{3}$ is strongly anharmonic

Halide perovskite (HP) semiconductors exhibit unique strong coupling between the electronic and structural dynamics. The high-temperature cubic phase of HPs is known to be entropically stabilized, with imaginary frequencies in the calculated phonon dispersion relation. Similar calculations, based on the static average crystal structure, predict a stable tetragonal phase with no imaginary modes. This work shows that in contrast to standard theory predictions, the room-temperature tetragonal phase of CH$_{3} $NH$_{3} $PbI$_{3}$ is strongly anharmonic. We use Raman polarization-orientation (PO) measurements and \textit{ab initio} molecular dynamics (AIMD) to investigate the origin and temperature evolution of the strong structural anharmonicity throughout the tetragonal phase. Raman PO measurements reveal a new spectral feature that resembles a soft mode. This mode shows an unusual continuous increase in damping with temperature which is indicative of an anharmonic potential surface. The analysis of AIMD trajectories identifies two major sources of anharmonicity: the orientational unlocking of the [CH$_{3} $NH$_{3}$]$^+$ ions and large amplitude octahedral tilting that continuously increases with temperature. Our work suggests that the standard phonon picture cannot describe the structural dynamics of tetragonal CH$_{3} $NH$_{3} $PbI$_{3}$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Halide perovskites under polarized light: Vibrational symmetry analysis using polarized Raman

In the last decade, hybrid organic-inorganic halide perovskites have emerged as a new type of semiconductor for photovoltaics and other optoelectronic applications. Unlike standard, tetrahedrally bonded semiconductors (e.g. Si and GaAs), the ionic thermal fluctuations in the halide perovskites (i.e. structural dynamics) are strongly coupled to the electronic dynamics. Therefore, it is crucial to obtain accurate and detailed knowledge about the nature of atomic motions within the crystal. This has proved to be challenging due to low thermal stability and the complex, temperature dependent structural phase sequence of the halide perovskites. Here, these challenges are overcome and a detailed analysis of the mode symmetries is provided in the low-temperature orthorhombic phase of methylammonium-lead iodide. Raman measurements using linearly- and circularly- polarized light at 1.16 eV excitation are combined with density functional perturbation theory (DFPT). By performing an iterative analysis of Raman polarization-orientation dependence and DFPT mode analysis, the crystal orientation is determined. Subsequently, accounting for birefringence effects detected using circularly polarized light excitation, the symmetries of all the observed Raman-active modes at 10 K are assigned.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Ideal near-Dirac triple-point semimetal in III-V semiconductor alloys

Despite the growing interest in topological materials, the difficulty of experimentally synthesizing and integrating them with other materials has been one of the main barriers restricting access to their unique properties. Recent advances in synthesizing metastable phases of crystalline materials can help to overcome this barrier and offer new platforms to experimentally study and manipulate band topology. Because III-V semiconductors have a wide range of functional material applications (including optoelectronic devices, light-emitting diodes, and highly efficient solar cells), and because Bi-doped III-V materials can be synthesized by ion plantation and ion-cutoff methods, we revisit the effect of bismuth substitution in metastable III-V semiconductors. Through first-principles calculation methods, we show that in wurtzite structure III-V materials, Bi substitution can lead to band inversion phenomena and induce nontrivial topological properties. Specifically, we identify that GaBi and InBi are Dirac-Weyl semimetals, characterized by the coexistence of Dirac points and Weyl points, and $\text{GaAs}_{0.5} \text{Bi}_{0.5}$, $\text{GaSb}_{0.5} \text{Bi}_{0.5}$, $\text{InSb}_{0.5} \text{Bi}_{0.5}$ are triple-point semimetals, characterized by two sets of "near Dirac" triple points on the Fermi level. These experimentally-accessible bismuth-based topological semimetals can be integrated into the large family of functional III-V materials for experimental studies of heterostructures and future optoelectronic applications.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Shift current bulk photovoltaic effect influenced by quasiparticles and excitons

We compute the shift current bulk photovoltaic effect (BPVE) in bulk BaTiO$_3$ and two-dimensional monochalcogenide SnSe considering quasi-particle corrections and exciton effects. We explore changes in shift current peak position and magnitude reduction due to band renormalization. For BaTiO$_3$, we demonstrate that shift current is reduced near the band edge due to exciton effects. Comparison of these results with experiments on BaTiO$_3$ indicate that mechanisms other than shift current may be contributing to BPVE. Additionally, we reveal that the shift current near the band gap shows only a small change due to excitons in two-dimensional SnSe, suggesting that the thin film geometry provides a feasible way to reduce the exciton effect on the shift current. These results suggest that many-body corrections are important for accurate assessments of bulk photovoltaic materials and to understand the mechanisms behind the BPVE

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Atomistic Description for Temperature-Driven Phase Transitions in BaTiO$_3$

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Classical model of negative thermal expansion in solids with expanding bonds

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Enhanced charge ordering transition in doped CaFeO3 through steric templating

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Enhancement of bulk photovoltaic effect in topological insulators

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bulk photovoltaic effect enhancement via electrostatic control in layered ferroelectrics

The correlation between the shift current mechanism for the bulk photovoltaic effect (BPVE) and the structural and electronic properties of ferroelectric perovskite oxides is not well understood. Here, we study and engineer the shift current photovoltaic effect using a visible-light-absorbing ferroelectric Pb(Ni$_{x}$Ti$_{1-x}$)O$_{3-x}$ solid solution from first principles. We show that the covalent orbital character dicates the direction, magnitude, and onset energy of shift current in a predictable fashion. In particular, we find that the shift current response can be enhanced via electrostatic control in layered ferroelectrics, as bound charges face a stronger impetus to screen the electric field in a thicker material, delocalizing electron densities. This heterogeneous layered structure with alternative photocurrent generating and insulating layers is ideal for BPVE applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dirac Line Nodes in Inversion Symmetric Crystals

We propose and characterize a new $\mathbb{Z}_2$ class of topological semimetals with a vanishing spin--orbit interaction. The proposed topological semimetals are characterized by the presence of bulk one-dimensional (1D) Dirac Line Nodes (DLNs) and two-dimensional (2D) nearly-flat surface states, protected by inversion and time--reversal symmetries. We develop the $\mathbb{Z}_2$ invariants dictating the presence of DLNs based on parity eigenvalues at the parity--invariant points in reciprocal space. Moreover, using first-principles calculations, we predict DLNs to occur in Cu$_3$N near the Fermi energy by doping non-magnetic transition metal atoms, such as Zn and Pd, with the 2D surface states emerging in the projected interior of the DLNs. This paper includes a brief discussion of the effects of spin--orbit interactions and symmetry-breaking as well as comments on experimental implications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

First-Principles Calculation of the Bulk Photovoltaic Effect in KNbO$_{3}$ and (K,Ba)(Ni,Nb)O$_{3-δ}$

The connection between noncentrosymmetric materials' structure, electronic structure, and bulk photovoltaic performance remains not well understood. In particular, it is still unclear which photovoltaic (PV) mechanism(s) are relevant for the recently demonstrated visible-light ferroelectric photovoltaic (K,Ba)(Ni,Nb)O$_{3-δ}$. In this paper, we study the bulk photovoltaic effect (BPVE) of (K,Ba)(Ni,Nb)O$_{3-δ}$ and KNbO$_{3}$ by calculating the shift current from first principles. The effects of structural phase, lattice distortion, oxygen vacancies, cation arrangement, composition, and strain on BPVE are systematically studied. We find that (K,Ba)(Ni,Nb)O$_{3-δ}$ has a comparable BPVE with that of the broadly explored BiFeO$_{3}$, but for a much lower photon energy. In particular, the Glass coefficient of (K,Ba)(Ni,Nb)O$_{5}$ in a simply layered structure can be as large as 12 times that of BiFeO$_{3}$. Furthermore, the nature of the wavefunctions dictates the eventual shift current yield, which can be significantly affected and engineered by changing the O vacancy location, cation arrangement, and strain. This is not only helpful for understanding other PV mechanisms that relate to the motion of the photocurrent carriers, but also provides guidelines for the design and optimization of PV converters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

First-principles materials design of high-performing bulk photovoltaics with the LiNbO$_3$ structure

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Layered Topological Crystalline Insulators

Topological crystalline insulators (TCIs) are insulating materials whose topological property relies on generic crystalline symmetries. Based on first-principles calculations, we study a three-dimensional (3D) crystal constructed by stacking two-dimensional TCI layers. Depending on the inter-layer interaction, the layered crystal can realize diverse 3D topological phases characterized by two mirror Chern numbers (MCNs) ($μ_1,μ_2$) defined on inequivalent mirror-invariant planes in the Brillouin zone. As an example, we demonstrate that new TCI phases can be realized in layered materials such as a PbSe (001) monolayer/h-BN heterostructure and can be tuned by mechanical strain. Our results shed light on the role of the MCNs on inequivalent mirror-symmetric planes in reciprocal space and open new possibilities for finding new topological materials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low Operation Voltage Ferroelectric Field-Effect Transistor Based on Polarization Rotation Effect

The effect of polarization rotation on the performance of metal oxide semiconductor field-effect transistors was investigated with a Landau-Ginzburg-Devonshire theory based model. In this analytical model, depolarization field, polarization rotations and the electrostatic properties of the doped silicon substrate are considered to illustrate the size effect of ferroelectric oxides and the stability of polarization in each direction. Based on this model, we demonstrate that MOSFET operation could be achieved by polarization reorientation with a low operating voltage, if the thickness of ferroelectric oxide is properly selected. Polarization reorientation can boost the surface potential of the silicon substrate, leading to a subthreshold swing S lower than 60 mV/decade. We believe that this model could provide guidance in designing electronic logic devices with low operating voltages and low active energy consumption.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Photoferroelectric and Photopiezoelectric Properties of Organometal Halide Perovskites

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Polarization Dependence of Water Adsorption to CH$_3$NH$_3$PbI$_3$ (001) Surfaces

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rashba Spin-Orbit Coupling Enhanced Carrier Lifetime in Organometal Halide Perovskites

Organometal halide perovskites are promising solar-cell materials for next-generation photovoltaic applications. The long carrier lifetime and diffusion length of these materials make them very attractive for use in light absorbers and carrier transporters. While these aspects of organometal halide perovskites have attracted the most attention, the consequences of the Rashba effect, driven by strong spin-orbit coupling, on the photovoltaic properties of these materials are largely unexplored. In this work, taking the electronic structure of methylammonium lead iodide as an example, we propose an intrinsic mechanism for enhanced carrier lifetime in 3D Rashba materials. Based on first-principles calculations and a Rashba spin-orbit model, we demonstrate that the recombination rate is reduced due to the spin-forbidden transition. These results are important for understanding the fundamental physics of organometal halide perovskites and for optimizing and designing the materials with better performance. The proposed mechanism including spin degrees of freedom offers a new paradigm of using 3D Rashba materials for photovoltaic applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strain Induced Ferroelectric Topological Insulator

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structural and Magnetic Reconstruction of Thermodynamically Stable CaMnO$_{3}$ (001) Surfaces

The relative thermodynamic stability of surface reconstructions including vacancies, adatoms and additional layers on both CaO and MnO$_{2}$ terminations is calculated to predict the surface phase diagram of ($\sqrt{2}\times\sqrt{2}$)$\it{R}$45$^{\circ}$ CaMnO$_{3}$ (001) using $\it{ab}$ $\it{initio}$ thermodynamics. Stoichiometric and nonstoichiometric reconstructions are considered. A set of boundary conditions driven by the binary and ternary sub-phases from CaMnO$_{3}$ defines its bulk stability region, enclosing the stable surface reconstructions that can be in equilibrium with the bulk. Most of the surfaces take the magnetic ordering of the bulk ground state, G-type antiferromagnetic (AFM); however, some of the MnO$_{2}$-terminated surfaces are more stable when the surface layer spins flip. At 573 K, the stoichiometric CaO- and MnO$_{2}$-terminated surfaces are predicted to be thermodynamically stable, as well as a CaO-terminated surface reconstruction where half the Ca are replaced by Mn. The MnO$_{2}$-terminated surface reconstructions dominate the phase diagrams at high temperatures, including phases with MnO and MnO$_{2}$ adatoms per surface unit cell.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Substantial optical dielectric enhancement by volume compression in LiAsSe$_2$

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Theoretical Modeling of Tribochemical Reaction on Pt and Au Contacts: Mechanical Load and Catalysis

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coupling between octahedral rotations and local polar displacements in WO3/ReO3 superlattices

We model short-period superlattices of WO$_3$ and ReO$_3$ with first-principles calculations. In fully-relaxed superlattices, we observe that octahedral tilts about an axis in the planes of the superlattices do not propagate from one material, despite the presence of the corner-shared oxygen atoms. However, we find that octahedral rotation is enhanced within WO$_3$ layers in cases in which strain couples with native antiferroelectric displacements of tungsten within their octahedral cages. Resulting structures remain antiferroelectric with low net global polarization. Thermodynamic analysis reveals that superlattices with sufficiently thick ReO$_3$ layers, the absolute number being three or more layers and the Re fraction $\geq 50\%$, tend to be more stable than the separated material phases and also show enhanced octahedral rotations in the WO$_3$ layers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Density Functional Theory Study Of Hypothetical PbTiO3-Based Oxysulfides

Using density functional theory (DFT) within the local density approximation (LDA), we calculate the physical and electronic properties of PbTiO3 (PTO) and a series of hypothetical compounds PbTiO3-xSx x = 0.2, 0.25, 0.33, 0.5, 1, 2, and 3 arranged in the corner-sharing cubic perovskite structure. We determine that replacing the apical oxygen atom in the PTO tetragonal unit cell with a sulfur atom reduces the x = 0 LDA calculated band gap of 1.47 eV to 0.43 - 0.67 eV for x = 0.2 - 1 and increases the polarization. PBE0 and GW methods predict that the hypothetical compositions x = 0.2 to x = 2 will have band gaps in the visible range. For all values of x < 2, the oxysulfide perovskite retains the tetragonal phase of PbTiO3, and the a-lattice parameter remains within 2.5% of the oxide. Thermodynamic analysis indicates that chemical routes using high temperature gas, such as H2S and CS2, can be used to substitute O for S in PTO for the compositions x = 0.2 - 0.5.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First-principles calculation of bulk photovoltaic effect in CH$_3$NH$_3$PbI$_3$ and CH$_3$NH$_3$PbI$_{3-x}$Cl$_{x}$

Hybrid halide perovskites exhibit nearly 20% power conversion efficiency, but the origin of their high efficiency is still unknown. Here, we compute the shift current, a dominant mechanism of bulk photovoltaic (PV) effect for ferroelectric photovoltaics, in CH$_3$NH$_3$PbI$_3$ and CH$_3$NH$_3$PbI$_{3-x}$Cl$_{x}$ from first principles. We find that these materials give approximately three times larger shift current PV response to near-IR and visible light than the prototypical ferroelectric photovoltaic BiFeO$_3$. The molecular orientations of CH$_3$NH$_3^{+}$ can strongly affect the corresponding PbI$_3$ inorganic frame so as to alter the magnitude of the shift current response. Specifically, configurations with dipole moments aligned in parallel distort the inorganic PbI$_3$ frame more significantly than configurations with near net zero dipole, yielding a larger shift current response. Furthermore, we explore the effect of Cl substitution on shift current, and find that Cl substitution at the equatorial site induces a larger response than does substitution at the apical site.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First-Principles Calculation of the Bulk Photovoltaic Effect in the Polar Compounds LiAsS$_\text{2}$, LiAsSe$_\text{2}$, and NaAsSe$_\text{2}$

We calculate the shift current response, which has been identified as the dominant mechanism for the bulk photovoltaic effect, for the polar compounds LiAsS$_\text{2}$, LiAsSe$_\text{2}$, and NaAsSe$_\text{2}$. We find that the magnitudes of the photovoltaic responses in the visible range for these compounds exceed the maximum response obtained for BiFeO$_\text{3}$ by 10 - 20 times. We correlate the high shift current response with the existence of $p$ states at both the valence and conduction band edges, as well as the dispersion of these bands, while also showing that high polarization is not a requirement. With low experimental band gaps of less than 2 eV and high shift current response, these materials have potential for use as bulk photovoltaics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Structural Diversity of ABS3 Compounds with d0 Electronic Configuration for the B-cation

We use first-principles density functional theory (DFT) within the local density approx- imation (LDA) to ascertain the ground state structure of real and theoretical compounds with the formula ABS3 (A = K, Rb, Cs, Ca, Sr, Ba, Tl, Sn, Pb, and Bi; and B = Sc, Y, Ti, Zr, V, and Nb) under the constraint that B must have a d0 electronic configuration. Our findings indicate that none of these AB combinations prefer a perovskite ground state with corner-sharing BS6 octahedra, but that they prefer phases with either edge- or face-sharing motifs. Further, a simple two-dimensional structure field map created from A and B ionic radii provides a neat demarcation between combinations preferring face-sharing versus edge- sharing phases for most of these combinations. We then show that by modifying the common Goldschmidt tolerance factor with a multiplicative term based on the electronegativity dif- ference between A and S, the demarcation between predicted edge-sharing and face-sharing ground state phases is enhanced. We also demonstrate that, by calculating the free energy contribution of phonons, some of these compounds may assume multiple phases as synthesis temperatures are altered, or as ambient temperatures rise or fall.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bulk Dirac points in distorted spinels

We report on a Dirac-like Fermi surface in three-dimensional bulk materials in a distorted spinel structure on the basis of density functional theory (DFT) as well as tight-binding theory. The four examples we provide in this paper are BiZnSiO4, BiCaSiO4, BiMgSiO4, and BiAlInO4. A necessary characteristic of these structures is that they contain a Bi lattice which forms a hierarchy of chain-like substructures, with consequences for both fundamental understanding and materials design.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exploration of the momentum of ferroelectric domain walls via molecular dynamics simulations

The motion of ferroelectric domain walls (DW) is critical for various technological applications of ferroelectric materials. One important question that is of interest both scientifically and technologically is whether the ferroelectric DW has momentum. To address this problem, we performed canonical ensemble molecular dynamics simulations of 180$^\circ$ and 90$^\circ$ DWs under applied electric field. Examination of the evolution of the polarization and local structure of DWs reveals that they stop moving after the removal of electric field. Thus, our computational study shows that ferroelectric domain walls do not have momentum.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reinterpretation of bond-valence model with bond-order formalism: an improved bond-valence based interatomic potential for PbTiO$_3$

We present a modified bond-valence model of PbTiO$_3$ based on the principles of bond-valence and bond-valence vector conservation. The relationship between the bond-valence model and the bond-order potential is derived analytically in the framework of a tight-binding model. A new energy term, bond-valence vector energy, is introduced into the atomistic model and the potential parameters are re-optimized. The new model potential can be applied both to canonical ensemble ($NVT$) and isobaric-isothermal ensemble ($NPT$) molecular dynamics (MD) simulations. This model reproduces the experimental phase transition in $NVT$ MD simulations and also exhibits the experimental sequence of temperature-driven and pressure-driven phase transitions in $NPT$ simulations. We expect that this improved bond-valence model can be applied to a broad range of inorganic materials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anisotropic Local Correlations and Dynamics in a Relaxor Ferroelectric

Relaxor ferroelectrics have been a focus of intense attention due to their anomalous dielectric characteristics, diffuse phase transitions, and strong piezoelectricity. Understanding the structure and dynamics of relaxors has been one of the long-standing challenges in solid-state physics, with the current model of polar nanoregions in a non-polar matrix providing only a qualitative description of the relaxor phase transitions. In this paper, we investigate the local structure and dynamics in 75%PbMg$_{1/3}$Nb$_{2/3}$O$_3$-25%PbTiO$_3$ (PMN-PT) using molecular dynamics simulations and the dynamic pair distribution function technique. We show for the first time that relaxor transitions can be described by local order parameters. We find that structurally, the relaxor phase is characterized by the presence of highly anisotropic correlations between the local cation displacements. These correlations resemble the hydrogen bond network in water. Our findings contradict the current polar nanoregion model; instead, we suggest a new model of a homogeneous random network of anisotropically coupled dipoles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

First principles calculations of the Shift Current Bulk Photovoltaic Effect in Ferroelectrics

We calculate the bulk photovoltaic response of the ferroelectrics BaTiO$_3$ and PbTiO$_3$ from first principles by applying "shift current" theory to the electronic structure from density functional theory. The first principles results for BaTiO$_3$ reproduce eperimental photocurrent direction and magnitude as a function of light frequency, as well as the dependence of current on light polarization, demonstrating that shift current is the dominant mechanism of the bulk photovoltaic effect in BaTiO$_3$. Additionally, we analyze the relationship between response and material properties in detail. The photocurrent does not depend simply or strongly on the magnitude of material polarization, as has been previously assumed; instead, electronic states with delocalized, covalent bonding that is highly asymmetric along the current direction are required for strong shift current enhancements. The complexity of the response dependence on both external and material parameters suggests applications not only in solar energy conversion, but to photocatalysis and sensor and switch type devices as well.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rigorous definition of oxidation states of ions in solids

We present justification and rigorous procedure for electron partitioning among atoms in extended systems. The method is based on wavefunction topology and the modern theory of polarization, rather than charge density partitioning or wavefunction projection, and, as such, re-formulates the concept of oxidation state without assuming real-space charge transfer between atoms. This formulation provides rigorous electrostatics of finite extent solids, including films and nanowires.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spin texture on the Fermi surface of tensile strained HgTe

We present ab initio and k.p calculations of the spin texture on the Fermi surface of tensile strained HgTe, which is obtained by stretching the zincblende lattice along the (111) axis. Tensile strained HgTe is a semimetal with pointlike accidental degeneracies between a mirror symmetry protected twofold degenerate band and two nondegenerate bands near the Fermi level. The Fermi surface consists of two ellipsoids which contact at the point where the Fermi level crosses the twofold degenerate band along the (111) axis. However, the spin texture of occupied states indicates that neither ellipsoid carries a compensating Chern number. Consequently, the spin texture is locked in the plane perpendicular to the (111) axis, exhibits a nonzero winding number in that plane, and changes winding number from one end of the Fermi ellipsoids to the other. The change in the winding of the spin texture suggests the existence of singular points. An ordered alloy of HgTe with ZnTe has the same effect as stretching the zincblende lattice in the (111) direction. We present ab initio calculations of ordered Hg_xZn_1-xTe that confirm the existence of a spin texture locked in a 2D plane on the Fermi surface with different winding numbers on either end.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lattice normal modes and electronic properties of the correlated metal LaNiO$_3$

We use density functional theory (DFT) calculations to study the lattice vibrations and electronic properties of the correlated metal LaNiO$_3$. To characterize the rhombohedral to cubic structural phase transition of perovskite LaNiO$_3$, we examine the evolution of the Raman-active phonon modes with temperature. We find that the $A_{1g}$ Raman mode, whose frequency is sensitive to the electronic band structure, is a useful signature to characterize the octahedral rotations in rhombohedral LaNiO$_3$. We also study the importance of electron--electron correlation effects on the atomic structure with two approaches which go beyond the conventional band theory (local spin density approximation): the local spin density+Hubbard $U$ method (LSDA$+U$) and hybrid exchange-correlation density functionals which include portions of exact Fock-exchange. We find the conventional LSDA accurately reproduces the delocalized nature of the valence states in LaNiO$_3$ and gives the best structural and vibrational agreement to the available experimental data. Based on our calculations, we show that the electronic screening effect from the delocalized Ni 3$d$ and O-2$p$ states mitigate the electronic correlations of the $d^7$ Ni cations, making LaNiO$_3$ a weakly correlated metal.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Polarization-dependence of palladium deposition on ferroelectric lithium niobate (0001) surfaces

We investigate the effect of ferroelectric polarization direction on the geometric properties of Pd deposited on the positive and negative surfaces of LiNbO$_3$ (0001). We predict preferred geometries and diffusion properties of small Pd clusters using density functional theory, and use these calculations as the basis for kinetic Monte Carlo simulations of Pd deposition on a larger scale. Our results show that on the positive surface, Pd atoms favor a clustered configuration, while on the negative surface, Pd atoms are adsorbed in a more dispersed pattern due to suppression of diffusion and agglomeration. This suggests that the effect of LiNbO$_3$ polarization direction on the catalytic activity of Pd [J. Phys. Chem. \textbf{88}, 1148 (1984)] is due, at least in part, to differences in adsorption geometry. Further investigations using these methods can aid the search for catalysts whose activities switch reversibly with the polarization of their ferroelectric substrates.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Post density functional theoretical studies of highly polar semiconductive Pb(Ti$_{1-x}$Ni$_{x}$)O$_{3-x}$ solid solutions: The effects of cation arrangement on band gap

We use a combination of conventional density functional theory (DFT) and post-DFT methods, including the local density approximation plus Hubbard $U$ (LDA+$U$), PBE0, and self-consistent $GW$ to study the electronic properties of Ni-substituted PbTiO$_{3}$ (Ni-PTO) solid solutions. We find that LDA calculations yield unreasonable band structures, especially for Ni-PTO solid solutions that contain an uninterrupted NiO$_{2}$ layer. Accurate treatment of localized states in transition-metal oxides like Ni-PTO requires post-DFT methods. $B$-site Ni/Ti cation ordering is also investigated. The $B$-site cation arrangement alters the bonding between Ni and O, and therefore strongly affects the band gap ($E_{\rm g}$) of Ni-PTO. We predict that Ni-PTO solid solutions should have a direct band gap in the visible light energy range, with polarization similar to the parent PbTiO$_{3}$. This combination of properties make Ni-PTO solid solutions promising candidate materials for solar energy conversion devices.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theoretical investigation of the evolution of the topological phase of Bi$_{2}$Se$_{3}$ under mechanical strain

The topological insulating phase results from inversion of the band gap due to spin-orbit coupling at an odd number of time-reversal symmetric points. In Bi$_2$Se$_3$, this inversion occurs at the $Γ$ point. For bulk Bi$_2$Se$_3$, we have analyzed the effect of arbitrary strain on the $Γ$ point band gap using Density Functional Theory. By computing the band structure both with and without spin-orbit interactions, we consider the effects of strain on the gap via Coulombic interaction and spin-orbit interaction separately. While compressive strain acts to decrease the Coulombic gap, it also increases the strength of the spin-orbit interaction, increasing the inverted gap. Comparison with Bi$_2$Te$_3$ supports the conclusion that effects on both Coulombic and spin-orbit interactions are critical to understanding the behavior of topological insulators under strain, and we propose that the topological insulating phase can be effectively manipulated by inducing strain through chemical substitution.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hybrid density functional calculations of the band gap of Ga$_x$In$_{1-x}$N

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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Structure and dielectric response in the high $T_c$ ferroelectric Bi(Zn,Ti)O$_3$-PbTiO$_3$ solid solutions

Theoretical {\em ab initio} and experimental methods were used to investigate the $x$Bi(Zn,Ti)O$_3$-(1-$x$)PbTiO$_3$ (BZT-PT) solid solution. We find that hybridization between Zn 4$p$ and O 2$p$ orbitals allows the formation of short, covalent Zn-O bonds, enabling favorable coupling between A-site and B-site displacements. This leads to large polarization, strong tetragonality and an elevated ferroelectric to paraelectric phase transition temperature. nhomogeneities in local structure near the 90$^\circ$ domain boundaries can be deduced from the asymetric peak broadening in the neutron and x-ray diffraction spectra. These extrinsic effects make the ferroelectric to paraelectric phase transition diffuse in BZT-PT solid solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Support-mediated activation and deactivation of Pt thin films

Using ab initio methods, we examine the the charge distribution at the interface of alpha-alumina-supported Pt films, and we consider the influence of this interface on CO adsorption. We demonstrate that a combination of electrostatic charge transfer and covalent bonding governs the interfacial interactions, and that these interactions play an important role in the metal reactivity. By modifying the interface and varying the Pt film thickness over a nanoscale range, CO adsorption can be significantly enhanced or diminished. These observations could be used to tune the reactivity of Pt particles.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Optimization of Quantum Monte Carlo Wave Functions Using Analytical Energy Derivatives

An algorithm is proposed to optimize quantum Monte Carlo (QMC) wave functions based on New ton's method and analytical computation of the first and second derivatives of the variati onal energy. This direct application of the variational principle yields significantly low er energy than variance minimization methods when applied to the same trial wave function. Quadratic convergence to the local minimum of the variational parameters is achieved. A g eneral theorem is presented, which substantially simplifies the analytic expressions of de rivatives in the case of wave function optimization. To demonstrate the method, the ground state energies of the first-row elements are calculated.

0 citations0 reviewsNo code linkedOpen access