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

Crystallization properties of arsenic doped GST alloys.

We present the enhanced properties observed in the phase change memory alloy Ge2Sb2Te5 (GST) when doped with arsenic. Although arsenic is known as a toxic element, our observations show that significant improvement can be obtained in GST systems on thermal stability, transition temperature between amorphous and crystalline phases and switching behaviors when doping with arsenic. Though both the GST and arsenic doped GST are amorphous in the as-deposited state, only GST alloy turns to crystalline NaCl-type structure after annealing at 150 °C for 1 h. Results from the resistance versus temperature study show a systematic increase in the transition temperature and resistivity in the amorphous and crystalline states when the arsenic percentage in the GST alloy increases. The crystallization temperature (Tc) of (GST)0.85As0.15 is higher than the Tc observed in GST. Optical band gap (Eopt) values of the as-deposited films show a clear increasing trend; 0.6 eV for GST to 0.76 eV for (GST)0.85As0.15. The decreases in Eopt for the samples annealed at higher temperatures shows significant optical contrast between the as-deposited and annealed samples. Though all (GST)1-xAsx alloys show memory switching behaviors, threshold switching voltages (VT) of the studied alloys show an increasing trend with arsenic doping. For (GST)0.85As0.15, VT is about 5.2 V, which is higher than GST (4.0 V). Higher transition temperature and higher threshold switching values show arsenic doping in GST can enhance the memory device properties by improving the thermal stability and data readability. Understanding the doping effect on the GST is important to understand its crystallization properties. Structure properties of amorphous GST, Ge2Sb2-0.3As0.3Te5 and (GST)0.85As0.15 models were studied using first principles molecular dynamics simulations, compared their partial radial distribution functions, and q parameter order. Arsenic doping into GST features interesting structural and electronic effects revealed by the radial distribution functions, q order parameter and band gap value, in line with the experimental findings.

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

Understanding the degradation mechanisms of lead-halide perovskites (CH3NH3PbI3) 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.(H2O)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 [Formula: see text] stretching vibrations allude to significant hydrogen bonding between MA and the water molecules. Through Fourier analysis of the rotational dynamics on several MA · (H2O)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., CH3NH3PbI3) 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 (CH3NH3PbI3) 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 CH3NH3PbI3 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 TiO2 (a-TiO2) could offer an attractive alternative to conventional crystalline TiO2 phases for photocatalytic applications. However, the atomic structure of a-TiO2 remains poorly understood with respect to that of its crystalline counterparts. Here, we conduct some classical molecular dynamics simulations of a-TiO2 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-TiO2. Overall, the results suggest that a-TiO2 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.

New insights in the electrocatalytic proton reduction and hydrogen oxidation by bioinspired catalysts: a DFT investigation.

In this paper, we present a DFT study of the proton reduction mechanism catalyzed by the complex [Ni(P2(H)N2(H))2](2+), bioinspired from the hydrogenases. A detailed analysis of the reactive isomers is discussed together with the localizations of the transitions states and energy minima. The reactive catalytic species is a biprotonated Ni(0) complex that can show different conformations and that can be protonated on different sites. The energies of the different conformations and biprotonated species have been calculated and discussed. Energy barriers for two different reaction mechanisms have been identified in solvent and in gas phase. Frequency calculations have been performed to check the nature of the energy minima and for the calculations of entropic energetic terms and zero point energies. We show that only one conformation is mostly reactive. All the others species are nonreactive in their original form, and they have to pass through conformational barriers in order to transform in the reactive species.

Conformational changes in a flexible, encapsulated dicarboxylate: evidence from density functional theory simulations.

The dynamical behavior of the [Mo12O12S12(OH)12{O2C-(CH2)N-CO2}]2- complexes is analyzed via first-principles molecular dynamics simulations. Experimental X-ray data play the role of initial configurations for our search in the configuration space. We show that different stable and metastable conformers are possible, and these are thermally accessible at relatively low temperatures. This is the main outcome of our first-principles molecular dynamics approach in which the temperature allowing for thermal activation has been set to T = 500 K, which is consistent with the variable temperature 1H NMR spectra of these complexes in solution taken at comparable, although moderately lower, temperature. This implies that a relatively large manifold of folding configurations is available to the encapsulated guest species. Moreover, the high flexibility of both the host cage and the inserted guests allows for the accommodation of a rather wide variety of functional groups with potential applications in several fields.

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.

Tuning the thermodynamic stability of oxothiomolybdenum wheels: crystal structures, studies in solution and DFT calculations.

The complexes [Mo12O12S12(OH)12(Muco)]2- (Muco2- = muconate, C6H4O4(2-)) and [Mo12,O12S12(OH)12(TMT)]2 (TMT2- = tetramethylterephthalate, C12H12O4(2-)) have been obtained from the condensation of the [Mo2O2S2]2+ building block in the presence of Muco2- and TMT2-, respectively. Both compounds were structurally characterized, revealing host-guest architectures with one or two encapsulated water molecules. 1H NMR spectra in DMSO and D2O showed that both complexes had an average symmetry higher than that in the solid state, due to changes in the distribution of encapsulated water molecules. The relative stabilities in water of the seven complexes encapsulating various di- or tricarboxylate guests, either rigid or non-rigid, have been determined. The stability scale obtained for the dianionic complexes is interpreted in relation with the rigidity or flexibility of the guests. A DFT study demonstrates that additional stabilization arises from the presence of inner hydrogen bonds involving 1, 2 or 3 water molecules, which even permit the extension of the H-bonds network to the first solvation sphere of the anion. DFT calculations were carried out on all investigated complexes as isolated or solvated anions and provide the sequence of the bond energies between the host and the guests, which is compared to the experimental data.

Changing the oxothiomolybdate ring from an anionic to a cationic receptor.

The ionic recognition properties of neutral oxothiomolybdenum wheels can be changed from anionic to cationic through ionization of the internal aquo ligands. In the solid state, [Mo8O8S8(OH)10(H2O)]2- (1) interacts with two Cs+ cations to give a close supramolecular host-guest arrangement. Such interactions appear to be maintained in solution as a labile association.

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.