First-Principles Molecular Dynamics in Metal-Halide Perovskites: Contrasting Generalized Gradient Approximation and Hybrid Functionals

First-principles molecular dynamics (FPMD) represents a valuable tool to probe dynamical properties of metal-halide perovskites (MHPs) which are key to their success in optoelectronic devices. Most FPMD studies rely on generalized gradient approximation (GGA) functionals for computational efficiency matters, while hybrid functionals, although computationally demanding, are usually needed to accurately describe structural and electronic properties of MHPs. This Letter reports FPMD simulations on CsPbI₃ based on the hybrid PBE0 functional. Our results demonstrate that PBE0 leads to lattice parameters and phonon modes in excellent agreement with experimental data, while GGA results overestimate the lattice parameter and the electronic band gap and underestimate the phonon energies. Our FPMD results also shed light on anharmonic effects and double-well instabilities in the octahedral tilting, highlighting a lowered free energy barrier for PBE0 and farther separated potential wells. Our results suggest that hybrid functionals are required to accurately describe crystal structure, lattice dynamics, and anharmonicity in MHPs.

Structural and Optical Properties of Solvated PbI2 in γ-Butyrolactone: Insight into the Solution Chemistry of Lead Halide Perovskite Precursors

We employ a fine-tuned theoretical framework, combining ab initio molecular dynamics (AIMD), density functional theory (DFT), and time-dependent (TD) DFT methods, to investigate the interactions and optical properties of the iodoplumbates within the low coordinative γ-butyrolactone (GBL) solvent environment, widely employed in the perovskite synthesis. We uncover the extent of GBL coordination to PbI2 investigating its relation to the solvated PbI2 optical properties. The employed approach has been further validated by comparison with the experimental UV-vis absorption spectrum of PbI2 in GBL solvent. A comparison with other solvents, commonly employed in the perovskite synthesis, such as N,N-dimethylformamide (DMF) and dimethyl sulfoxide (DMSO) is also reported. The methodology developed in this work can be reasonably extended to the investigation of similar systems.

Role of Water on the Rotational Dynamics of the Organic Methylammonium Cation: A First Principles Analysis

Understanding the degradation mechanisms of lead-halide perovskites (CH₃NH₃PbI₃) under exposure to liquid/aerosol water is an essential problem within the photovoltaic community. Herein we investigate both the static and the dynamic properties of the methylammonuim cation (MA) as it coordinates with invading water molecules (MA.(H₂O)_n, n = 1, 2, 3, 4) using both stationary state quantum mechanics and first principle molecular dynamics simulations. Various solvation structures of MA were characterized by their stabilization energies, dipoles, and Maximally-Localized Wannier Function (MLWF) centers. Calculation – and analysis – of vibrational shifts in the IR spectral region were performed for hydrated complexes; the locations of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bf{N}}{{\bf{H}}}_{{\bf{3}}}^{{\boldsymbol{+}}}$$\end{document}NH3+ stretching vibrations allude to significant hydrogen bonding between MA and the water molecules. Through Fourier analysis of the rotational dynamics on several MA · (H₂O)n complexes, we conclude that the water molecules dampen the rotational motion of the MA as the intermolecular bonds formed between the water molecules and the MA act to hinder the rotation of the cation; these findings give explanatory support to earlier computational observations of humidity effects on perovskites (i.e., CH₃NH₃PbI₃) materials. This work is a step toward understanding the water-MA cation interaction in bulk perovskites, thus providing greater understanding of in situ instability/degradation of perovskite bulk materials.

Effects of Electron-Phonon Coupling on Electronic Properties of Methylammonium Lead Iodide Perovskites

Temperature can have a dramatic effect on the solar efficiency of methylammonium lead iodide (CH₃NH₃PbI₃) absorbers due to changes in the electronic structure of the system even within the range of stability of a single phase. Herein, using first-principles density functional theory, we investigate the electron band structure of the tetragonal and orthorhombic phases of CH₃NH₃PbI₃ as a function of temperature. The electron-phonon interactions are computed to all orders using a Monte Carlo approach, which is needed considering that the second-order Allen-Heine-Cardona theory in electron-phonon coupling is not adequate. Our results show that the band gap increases with temperature, in excellent agreement with experimental results. We verified that anharmonic effects are only important near the tetragonal-cubic phase transition temperature. We also found that temperature has a significant effect on the effective masses and Rashba coupling. At room temperature, electron-phonon coupling is found to enhance the band effective mass by a factor of 2 and to diminish the Rashba coupling by the same factor compared to T = 0 K values. Our results underscore the significant impact of electron-phonon coupling on electronic properties of the hybrid perovskites.

New insights into the atomic structure of amorphous TiO2 using tight-binding molecular dynamics

Amorphous TiO₂ (a-TiO₂) could offer an attractive alternative to conventional crystalline TiO₂ phases for photocatalytic applications. However, the atomic structure of a-TiO₂ remains poorly understood with respect to that of its crystalline counterparts. Here, we conduct some classical molecular dynamics simulations of a-TiO₂ based on a selection of empirical potentials. We show that, on account of its ability to dynamically assign the charge of each atom based on its local environment, the second-moment tight-binding charge equilibration potential yields an unprecedented agreement with available experimental data. Based on these simulations, we investigate the degree of order and disorder in a-TiO₂. Overall, the results suggest that a-TiO₂ features a large flexibility in its local topology, which may explain the high sensitivity of its structure to the synthesis method being used.

Mapping the Free Energy of Lithium Solvation in the Protic Ionic Liquid Ethylammonuim Nitrate: A Metadynamics Study

Understanding lithium solvation and transport in ionic liquids is important due to their possible application in electrochemical devices. Using first-principles simulations aided by a metadynamics approach we study the free-energy landscape for lithium ions at infinite dilution in ethylammonium nitrate, a protic ionic liquid. We analyze the local structure of the liquid around the lithium cation and obtain a quantitative picture in agreement with experimental findings. Our simulations show that the lowest two free energy minima correspond to conformations with the lithium ion being solvated either by three or four nitrate ions with a transition barrier between them of 0.2 eV. Other less probable conformations having different solvation pattern are also investigated.

Dynamic properties of a hexadecamolybdenum wheel: studies in solution and density functional theory calculations

Variable temperature (1)H NMR studies of the host-guest complex [Mo(16)O(16)S(16)(OH)(16)(H(2)O)(4)(PDA)(2)](4-) (1 ; PDA(2-) = phenylenediacetate) previously carried out in D(2)O have revealed a complex behavior in solution, involving a gliding motion of both parallel phenyl rings of the PDA(2-) ligands. In the present work, we present new NMR spectra carried out in the aprotic solvent CD(3)CN, which allow the observation of the proton signals associated with the bridging hydroxo groups of the inorganic host. The new spectra provide detailed information about the concerted reorganization of the guest components, that is, PDA(2-) and water molecules. The existence of an equilibrium between two distinct isomers differing in the linking mode between the inorganic host and the two equivalent PDA(2-) ligands is evidenced. This equilibrium appears strongly dependent upon the temperature, leading to a complete inversion of the distribution between 300 and 226 K. The thermodynamic data related to the isomerization reaction have been determined (Delta(r)H = -50.5 kJ mol(-1) and Delta(r)S = -215 J mol(-1) K(-1)). Furthermore, at low temperature, one of the isomers exists in two conformations, only differing in the H-bond network involving the inner water molecules. Density functional theory calculations were carried out to push ahead the interpretations obtained from experiment, identify the isomers of 1, and specify the role and the positions of the guest water molecules. Among the various structures that have been calculated for 1, three fall in a narrow energy range and should correspond to the species characterized by variable-temperature (1)H NMR experiments in CD(3)CN. Finally, this study shows how the internal disposition of the ligands affects the ellipticity of the Mo(16) ring which varies from one isomer to the other in the 0.73-1 range and highlights solvation of the ring as one of the key parameters for the conformational design of these flexible host-guest systems.

Host–guest adaptability within oxothiomolybdenum wheels: structures, studies in solution and DFT calculations

The formation of host-guest cyclic architectures, built up through the self-condensation process of [Mo(2)O(2)S(2)](2+) oxothiocations around linear dicarboxylate ions such as adipate (Adip(2-)), suberate (Sub(2-)) and azelaate (Azel(2-)) anions is reported. The complexes [Mo(12)Adip](2-), [Mo(12)Sub](2-) and [Mo(14)Azel](2-) have been characterized in the solid state by X-ray diffraction and in solution by (1)H NMR in different solvents (D(2)O, DMF, DMSO and CD(3)CN). The host-guest dynamics appear to be dependent on the nature of the system and are mainly governed by mutual adaptability between the host and the guest. (1)H NMR DOSY experiments show systematic differences, either positive or negative between the experimental and calculated molecular weights which appear to be correlated with the charge of the anion. The relative stabilities of the twelve-membered rings containing the Adip(2-), Pim(2-) (pimelate) or Sub(2-) anions were determined experimentally and decrease according to the order [Mo(12)Adip](2-) > [Mo(12)Pim](2-) > [Mo(12)Sub](2-). The host-guest adaptability depends on the length of the carbon chain and gives rise to selective encapsulation processes. Finally, theoretical DFT investigations in the gas phase yielded conformations whose symmetry and geometrical parameters proved consistent with X-ray structures and (1)H NMR spectra recorded in DMSO or DMF. Energy calculation highlights the high flexibility of the ring showing that only 3.1 kJ mol(-1) accompanies the conformational change from circular to elliptical. The host-guest bond energy (Delta E) calculated for the Mo(12)-based clusters is consistent with the experimental stability scale, major variations being due to some constraints undergone by the central alkyl chain.