Showing 1–8 of 8 results for author: Lahnsteiner, J

Tuning Einstein Oscillator Frequencies of Cation Rattlers: A Molecular Dynamics Study of the Lattice Thermal Conductivity of CsPbBr$_3$

Authors: Jonathan Lahnsteiner, Max Rang, Menno Bokdam

Abstract: The pure CsPbBr$_3$ perovskite is an archetypal example of a strongly anharmonic crystal that poses a major challenge for computational methods to describe its thermodynamic properties. Its lattice dynamics exhibits characteristics of a phonon liquid: mode coupling, low lifetimes and `rattlers'. To study the thermal conduction in this crystal, including the effect of dynamic disorder introduced by… ▽ More The pure CsPbBr$_3$ perovskite is an archetypal example of a strongly anharmonic crystal that poses a major challenge for computational methods to describe its thermodynamic properties. Its lattice dynamics exhibits characteristics of a phonon liquid: mode coupling, low lifetimes and `rattlers'. To study the thermal conduction in this crystal, including the effect of dynamic disorder introduced by the Cs rattlers, we apply large-scale molecular dynamics (MD) simulations combined with machine-learning interatomic potentials. We simulate its ultra-low lattice thermal conductivity in the cubic phase and obtain phonon spectra by measuring velocity autocorrelation functions. The thermal conductivity at 500 K is computed to be $0.53\pm0.04\frac{\rm W}{\rm m K}$, which is similar to that of demineralized water under normal indoor conditions. MD based insight into the heat transport mechanism of halide perovskites is presented. In the analysis the Cs cations are interpreted as damped Einstein oscillators. The phonon bandstructure of a system with artificially raised Cs masses demonstrates an increased interference of the Cs rattling with the acoustic phonon modes. We show that the thermal conductivity of the CsPbBr$_3$ perovskite can still be slightly decreased by tuning the cation rattling frequency into the range of the low lying acoustic mode △ Less

Submitted 3 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the J. Phys. Chem. C on 18 Dec. 2023, and in the upcoming Machine Learning in Physical Chemistry Volume 2 VSI Virtual Special Issue

Exploring librational pathways with on-the-fly machine-learning force fields: Methylammonium molecules in MAPbX$_3$ (X=I, Br, Cl) perovskites

Authors: Menno Bokdam, Jonathan Lahnsteiner, D. D. Sarma

Abstract: Two seemingly similar crystal structures of the low-temperature (~100 K) MAPbX$_3$ (X=I,Br,Cl) perovskites, but with different relative Methylammonium (MA) ordering, have appeared as representatives of this orthorhombic phase. Distinguishing them by X-ray diffraction experiments is difficult and conventional first-principles based molecular-dynamics approaches are often too computationally intensi… ▽ More Two seemingly similar crystal structures of the low-temperature (~100 K) MAPbX$_3$ (X=I,Br,Cl) perovskites, but with different relative Methylammonium (MA) ordering, have appeared as representatives of this orthorhombic phase. Distinguishing them by X-ray diffraction experiments is difficult and conventional first-principles based molecular-dynamics approaches are often too computationally intensive to be feasible. Therefore, to determine the thermodynamically stable structure, we use a recently introduced on-the-fly Machine-Learning Force Field method, which reduces the computation time from years to days. The molecules exhibit a large degree of anharmonic motion depending on temperature: i.e. rattling, twisting and tumbling. We observe the crystal's 'librational pathways' while slowly heating it in isothermal-isobaric simulations. Marked differences in the thermal evolution of structural parameters allow us to determine the real structure of the system via a comparison with experimentally determined crystal structures. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Journal ref: J. Phys. Chem. C 2021, 125, 38, 21077-21086

Anharmonic lattice dynamics in large thermodynamic ensembles with machine-learning force fields: CsPbBr$_{3}$ a phonon-liquid with Cs rattlers

Authors: Jonathan Lahnsteiner, Menno Bokdam

Abstract: The phonon dispersion relations of crystal lattices can often be well-described with the harmonic approximation. However, when the potential energy landscape exhibits more anharmonicity, for instance, in case of a weakly bonded crystal or when the temperature is raised, the approximation fails to capture all crystal lattice dynamics properly. Phonon-phonon scattering mechanisms become important an… ▽ More The phonon dispersion relations of crystal lattices can often be well-described with the harmonic approximation. However, when the potential energy landscape exhibits more anharmonicity, for instance, in case of a weakly bonded crystal or when the temperature is raised, the approximation fails to capture all crystal lattice dynamics properly. Phonon-phonon scattering mechanisms become important and limit the phonon lifetimes. We take a novel approach and simulate the phonon dispersion of a complex dynamic solid at elevated temperatures with Machine-Learning Force Fields of near-first-principles accuracy. Through large-scale molecular dynamics simulations the projected velocity autocorrelation function (PVACF) is obtained. We apply this approach to the inorganic perovskite CsPbBr$_{3}$. Imaginary modes in the harmonic picture of this perovskite are absent in the PVACF, indicating a dynamic stabilization of the crystal. The anharmonic nature of the potential makes a decoupling of the system into a weakly interacting phonon gas impossible. The phonon spectra of CsPbBr$_{3}$ show the characteristics of a phonon liquid. Rattling motions of the Cs$^{+}$ cations are studied by self-correlation functions and are shown to be nearly dispersionless motions of the cations with a frequency of $\sim$0.8THz within the lead-bromide framework. △ Less

Submitted 14 December, 2021; v1 submitted 15 January, 2021; originally announced January 2021.

Comments: 16 pages , 13 figures

Long-range order imposed by short-range interactions in methylammonium lead iodide: Comparing point-dipole models to machine-learning force fields

Authors: Jonathan Lahnsteiner, Ryosuke Jinnouchi, Menno Bokdam

Abstract: The crystal structure of the MAPbI$_3$ hybrid perovskite forms an intricate electrostatic puzzle with different ordering patterns of the MA molecules at elevated temperatures. For this perovskite three published model Hamiltonians based on the point-dipole (pd) approximation combined with short-range effective interactions are compared to a recently developed machine-learning force field. A molecu… ▽ More The crystal structure of the MAPbI$_3$ hybrid perovskite forms an intricate electrostatic puzzle with different ordering patterns of the MA molecules at elevated temperatures. For this perovskite three published model Hamiltonians based on the point-dipole (pd) approximation combined with short-range effective interactions are compared to a recently developed machine-learning force field. A molecular order parameter is used to consistently compare the transformation of the anti-ferroelectric ordering in the orthorhombic phase upon raising the temperature. We show that the ground states and the order-disorder transition of the three models are completely different. Our analysis indicates that the long-range order in the low-temperature orthorhombic phase can be captured by pd-based models with a short cutoff radius, including the nearest and next-nearest neighbor molecules. By constructing effective atomic interactions the ordering can already be described within 6 A radius. By extracting the coupling energetics of the molecules from density functional theory calculations on MA$_x$Cs$_{1-x}$PbI$_3$ test systems, we show that the pd-approximation holds at least for static structures. To improve the accuracy of the pd-interaction an Ewald summation is applied combined with a distance dependent electronic screening function. △ Less

Submitted 13 September, 2019; v1 submitted 29 May, 2019; originally announced May 2019.

Journal ref: Phys. Rev. B 100, 094106 (2019)

Phase transitions of hybrid perovskites simulated by machine-learning force fields trained on-the-fly with Bayesian inference

Authors: Ryosuke Jinnouchi, Jonathan Lahnsteiner, Ferenc Karsai, Georg Kresse, Menno Bokdam

Abstract: Realistic finite temperature simulations of matter are a formidable challenge for first principles methods. Long simulation times and large length scales are required, demanding years of compute time. Here we present an on-the-fly machine learning scheme that generates force fields automatically during molecular dynamics simulations. This opens up the required time and length scales, while retaini… ▽ More Realistic finite temperature simulations of matter are a formidable challenge for first principles methods. Long simulation times and large length scales are required, demanding years of compute time. Here we present an on-the-fly machine learning scheme that generates force fields automatically during molecular dynamics simulations. This opens up the required time and length scales, while retaining the distinctive chemical precision of first principles methods and minimizing the need for human intervention. The method is widely applicable to multi-element complex systems. We demonstrate its predictive power on the entropy driven phase transitions of hybrid perovskites, which have never been accurately described in simulations. Using machine learned potentials, isothermal-isobaric simulations give direct insight into the underlying microscopic mechanisms. Finally, we relate the phase transition temperatures of different perovskites to the radii of the involved species, and we determine the order of the transitions in Landau theory. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Journal ref: Phys. Rev. Lett. 122, 225701 (2019)

The Finite Temperature Structure of the MAPbI3 Perovskite: Comparing Density Functional Approximations and Force Fields to Experiment

Authors: Jonathan Lahnsteiner, Georg Kresse, Jurn Heinen, Menno Bokdam

Abstract: Determining the finite temperature structure of the hybrid perovskite MAPbI3 is a challenge for both experimental and theoretical methods. A very powerful computational method that can resolve the atomic structure is molecular dynamics (MD). The resulting structure depends on the density functional approximation (DFA) in the case of ab initio MD and the force field in classical MD. We compare the… ▽ More Determining the finite temperature structure of the hybrid perovskite MAPbI3 is a challenge for both experimental and theoretical methods. A very powerful computational method that can resolve the atomic structure is molecular dynamics (MD). The resulting structure depends on the density functional approximation (DFA) in the case of ab initio MD and the force field in classical MD. We compare the structure between 250K and 400K obtained with different DFAs and force fields in one consistent manner. The symmetry of the PbI3 framework is analyzed as well as the relative ordering of the neighboring organic molecules inside the framework. The distribution function of the molecules is used to map out an effective energy surface for the rotation of a single molecule. This surface is accurately modeled by a pair of cubic harmonics. Available experimental data in literature are discussed and compared to the structure obtained with the different methods. The spread in these data is still too large to uniquely determine the method that 'best' describes the perovskite, however promising candidates and outliers have been identified. △ Less

Submitted 23 April, 2018; originally announced April 2018.

Comments: 15 pages, 8 figures

Journal ref: Phys. Rev. Materials 2, 073604 (2018)

Assessing density functionals using many body theory for hybrid perovskites

Authors: Menno Bokdam, Jonathan Lahnsteiner, Benjamin Ramberger, Tobias Schaefer, Georg Kresse

Abstract: Which density functional is the "best" for structure simulations of a particular material? A concise, first principles, approach to answer this question is presented. The random phase approximation (RPA)--- an accurate many body theory--- is used to evaluate various density functionals. To demonstrate and verify the method, we apply it to the hybrid perovskite MAPbI$_3$, a promising new solar cell… ▽ More Which density functional is the "best" for structure simulations of a particular material? A concise, first principles, approach to answer this question is presented. The random phase approximation (RPA)--- an accurate many body theory--- is used to evaluate various density functionals. To demonstrate and verify the method, we apply it to the hybrid perovskite MAPbI$_3$, a promising new solar cell material. The evaluation is done by first creating finite temperature ensembles for small supercells using RPA molecular dynamics, and then evaluating the variance between the RPA and various approximate density functionals for these ensembles. We find that, contrary to recent suggestions, van der Waals functionals do not improve the description of the material, whereas hybrid functionals and the SCAN (strongly constrained appropriately normed) density functional yield very good agreement with the RPA. Finally, our study shows that in the room temperature tetragonal phase of MAPbI$_3$, the molecules are preferentially parallel to the shorter lattice vectors but reorientation on ps timescales is still possible. △ Less

Submitted 13 October, 2017; v1 submitted 22 August, 2017; originally announced August 2017.

Journal ref: Phys. Rev. Lett. 119, 145501 (2017)

Room temperature dynamic correlation between methylammonium molecules in lead-iodine based perovskites: An ab-initio molecular dynamics perspective

Authors: Jonathan Lahnsteiner, Georg Kresse, Abhinav Kumar, D. D. Sarma, Cesare Franchini, Menno Bokdam

Abstract: The high efficiency of lead organo-metal-halide perovskite solar cells has raised many questions about the role of the methylammonium (MA) molecules in the Pb-I framework. Experiments indicate that the MA molecules are able to 'freely' spin around at room temperature even though they carry an intrinsic dipole moment. We have performed large supercell (2592 atoms) finite temperature ab-initio molec… ▽ More The high efficiency of lead organo-metal-halide perovskite solar cells has raised many questions about the role of the methylammonium (MA) molecules in the Pb-I framework. Experiments indicate that the MA molecules are able to 'freely' spin around at room temperature even though they carry an intrinsic dipole moment. We have performed large supercell (2592 atoms) finite temperature ab-initio molecular dynamics calculations to study the correlation between the molecules in the framework. An underlying long range anti-ferroelectric ordering of the molecular dipoles is observed. The dynamical correlation between neighboring molecules shows a maximum around room temperature in the mid-temperature phase. In this phase, the rotations are slow enough to (partially) couple to neighbors via the Pb-I cage. This results in a collective motion of neighboring molecules in which the cage acts as the mediator. At lower and higher temperatures the motions are less correlated. △ Less

Submitted 2 January, 2017; v1 submitted 17 August, 2016; originally announced August 2016.

Journal ref: Phys. Rev. B 94, 214114 (2016)