Transferable Classical Force Field for Pure and Mixed Metal Halide Perovskites Parameterized from First-Principles.

Many key features in photovoltaic perovskites occur in relatively long time scales and involve mixed compositions. This requires realistic but also numerically simple models. In this work we present a transferable classical force field to describe the mixed hybrid perovskite MAxFA1-xPb(BryI1-y)3 for variable composition (∀x, y ∈ [0, 1]). The model includes Lennard-Jones and Buckingham potentials to describe the interactions between the atoms of the inorganic lattice and the organic molecule, and the AMBER model to describe intramolecular atomic interactions. Most of the parameters of the force field have been obtained by means of a genetic algorithm previously developed to parametrize the CsPb(BrxI1-x)3 perovskite (Balestra et al. J. Mater. Chem. A. 2020, DOI: 10.1039/d0ta03200j). The algorithm finds the best parameter set that simultaneously fits the DFT energies obtained for several crystalline structures with moderate degrees of distortion with respect to the equilibrium configuration. The resulting model reproduces correctly the XRD patterns, the expansion of the lattice upon I/Br substitution, and the thermal expansion coefficients. We use the model to run classical molecular dynamics simulations with up to 8600 atoms and simulation times of up to 40 ns. From the simulations we have extracted the ion diffusion coefficient of the pure and mixed perovskites, presenting for the first time these values obtained by a fully dynamical method using a transferable model fitted to first-principles calculations. The values here reported can be considered as the theoretical upper limit, that is, without grain boundaries or other defects, for ion migration dynamics induced by halide vacancies in photovoltaic perovskite devices under operational conditions.

Efficient Computation of Structural and Electronic Properties of Halide Perovskites Using Density Functional Tight Binding: GFN1-xTB Method.

In recent years, metal halide perovskites (MHPs) for optoelectronic applications have attracted the attention of the scientific community due to their outstanding performance. The fundamental understanding of their physicochemical properties is essential for improving their efficiency and stability. Atomistic and molecular simulations have played an essential role in the description of the optoelectronic properties and dynamical behavior of MHPs, respectively. However, the complex interplay of the dynamical and optoelectronic properties in MHPs requires the simultaneous modeling of electrons and ions in relatively large systems, which entails a high computational cost, sometimes not affordable by the standard quantum mechanics methods, such as density functional theory (DFT). Here, we explore the suitability of the recently developed density functional tight binding method, GFN1-xTB, for simulating MHPs with the aim of exploring an efficient alternative to DFT. The performance of GFN1-xTB for computing structural, vibrational, and optoelectronic properties of several MHPs is benchmarked against experiments and DFT calculations. In general, this method produces accurate predictions for many of the properties of the studied MHPs, which are comparable to DFT and experiments. We also identify further challenges in the computation of specific geometries and chemical compositions. Nevertheless, we believe that the tunability of GFN1-xTB offers opportunities to resolve these issues and we propose specific strategies for the further refinement of the parameters, which will turn this method into a powerful computational tool for the study of MHPs and beyond.

Phase Transition Induced by Gas Adsorption in Metal-Organic Frameworks.

We present a molecular simulation study with the aim of investigating the structural phase transition of ZJU-198 metal-organic framework. This material has been recently synthetized with the appropriate control of window size, which performs well for the separation of mixtures of gases containing nitrogen and methane. We find that the adsorption of small gases in this structure is unusual, and provide an explanation of the molecular mechanisms involved. Using molecular simulation, we analyze the structural distortions exerted by the adsorption of carbon dioxide, nitrogen, methane, acetylene, and ethene. We found that the separation of mixtures composed of these gases in ZJU-198 is due to the organic linker of the structure. The rotation of this linker causes the expansion of the cavities and enhances gas separation by allowing the adsorption of molecules that a priori are too big to be adsorbed.

Adsorption of Cyclohexane in Pure Silica Zeolites: High-Throughput Computational Screening Validated by Experimental Data.

Adsorption of cyclohexane in pure silica zeolites was studied experimentally and by molecular simulations. Based on the adsorption isobars obtained from the quasi-equilibrated temperature adsorption and desorption (QE-TPDA) measurements and reported adsorption isotherms for high-silica zeolites Y, ZSM-5, and ZSM-11 we refined Lennard-Jones parameters for guest-host interactions available in the literature. Adsorption of cyclohexane from equimolar mixture of twisted-boat and chair conformations has been screened in 171 pure silica zeolitic structures using grand canonical Monte Carlo simulations. Almost 20 frameworks showing extraordinary preference for adsorption of the chair conformation over the twisted boat one or vice versa were found. This selectivity was attributed to the geometry of channels and cavities present in the pore structures, as all t-boat selective structures possess channels or cavities of 8.3-9.1 Å. We also differentiated ways of chair-selectivity depending on the size and shape of the channels or cavities and also on the arrangement of the guest molecules in the pores.

Molecular Dynamics Analysis of Charge Transport in Ionic-Liquid Electrolytes Containing Added Salt with Mono, Di, and Trivalent Metal Cations.

Among many other applications, room-temperature ionic liquids (ILs) are used as electrolytes for storage and energy-conversion devices. In this work, we investigate, at the microscopic level, the structural and dynamical properties of 1-methyl-1-butyl-pyrrolidinium bis(trifluoromethanesulfonyl) imide [C4 PYR]+ [Tf2 N]- IL-based electrolytes for metal-ion batteries. We carried out molecular dynamics simulations of electrolytes mainly composed of [C4 PYR]+ [Tf2 N]- IL with the addition of Mn+ -[Tf2 N]- metal salts (M=Li+ , Na+ , Ni2+ , Zn2+ , Co2+ , Cd2+ , and Al3+ , n=1, 2, and 3) dissolved in the IL. The addition of low salt concentrations lowers the charge transport and conductivity of the electrolytes. This effect is due to the strong interaction of the metal cations with the [Tf2 N]- anions, which allows for molecular aggregation between them. We analyze how the conformation of the [Tf2 N]- anions surrounding the metal cations determine the charge-transport properties of the electrolyte. We found two main conformations based on the size and charge of the metal cation: monodentate and bidentate (number of oxygen atoms of the anion pointing to the metal atoms). The microscopic local structure of the Mn+ -[Tf2 N]- aggregates influences the microscopic charge transport as well as the macroscopic conductivity of the total electrolyte.

Stepped Propane Adsorption in Pure-Silica ITW Zeolite.

Gas adsorption over zeolites is at the basis of important applications of this class of microporous crystalline solids, notably as separation media and catalysts, but it may also be complex and not straightforward to understand. Here we report that for temperature below 323 K propane adsorption on the small-pore pure-silica zeolite ITW exhibits a clear step (pseudosaturation). This is absent in the case of propene and the other small linear alkanes. An intermediate plateau, clearly observed in the 293 K isotherm, always occurs when one molecule of propane is loaded in every other cage, i.e., at half-saturation. The simulation results show a swelling of the ITW structure upon propane adsorption. The strong dependence of available pore volume on the adsorbate loading level implies that adsorption cannot occur on the void structure while saturation can only be reached on highly loaded structures. To account for this unprecedented adsorption phenomenon, we propose the term "guest-modulated effect".

Potential of polarizable force fields for predicting the separation performance of small hydrocarbons in M-MOF-74.

The separation of light olefins from paraffins via cryogenic distillation is a very energy intensive process. Solid adsorbents and especially metal-organic frameworks with open metal sites have the potential to significantly lower the required energy. Specifically, M-MOF-74 has drawn considerable attention for application in olefin/paraffin separation. To investigate how the separation proceeds on a molecular level and to design better materials, molecular simulation can be a useful tool. Unfortunately, it is still a challenge to model the adsorption behavior of many adsorbates in metal-organic frameworks with open metal sites. Previously, the inclusion of explicit polarization has been suggested to improve the quality of classical force fields for such systems. Here, the potential of polarizable force fields for the description of olefins and paraffins in metal-organic frameworks with open metal sites is investigated. In particular, heats of adsorption, binding geometries, and adsorption isotherms are calculated for C2H4, C2H6, C3H6, and C3H8 in M-MOF-74 (with M = Co, Mn, Fe, and Ni). In this study, no force field parameters are adjusted to improve the model. The results show that including explicit polarization significantly improves the description of the adsorption in comparison to non-polarizable generic force fields which do not consider explicit polarization. The study also reveals that simulation predictions are sensitive to the assigned repulsive potential and framework charges. A fully re-parametrized polarizable force field may have the capability to improve the predictions even further.

Aqueous Solutions of Ionic Liquids: Microscopic Assembly.

Aqueous solutions of ionic liquids are of special interest, due to the distinctive properties of ionic liquids, in particular, their amphiphilic character. A better understanding of the structure-property relationships of such systems is hence desirable. One of the crucial molecular-level interactions that influences the macroscopic behavior is hydrogen bonding. In this work, we conduct molecular dynamics simulations to investigate the effects of ionic liquids on the hydrogen-bond network of water in dilute aqueous solutions of ionic liquids with various combinations of cations and anions. Calculations are performed for imidazolium-based cations with alkyl chains of different lengths and for a variety of anions, namely, [Br](-), [NO3](-), [SCN](-) [BF4](-), [PF6](-), and [Tf2N](-). The structure of water and the water-ionic liquid interactions involved in the formation of a heterogeneous network are analyzed by using radial distribution functions and hydrogen-bond statistics. To this end, we employ the geometric criterion of the hydrogen-bond definition and it is shown that the structure of water is sensitive to the amount of ionic liquid and to the anion type. In particular, [SCN](-) and [Tf2N](-) were found to be the most hydrophilic and hydrophobic anions, respectively. Conversely, the cation chain length did not influence the results.

Quantum and Classical Molecular Dynamics of Ionic Liquid Electrolytes for Na/Li-based Batteries: Molecular Origins of the Conductivity Behavior.

Compositional effects on the charge-transport properties of electrolytes for batteries based on room-temperature ionic liquids (RTILs) are well-known. However, further understanding is required about the molecular origins of these effects, in particular regarding the replacement of Li by Na. In this work, we investigate the use of RTILs in batteries, by means of both classical molecular dynamics (MD), which provides information about structure and molecular transport, and ab initio molecular dynamics (AIMD), which provides information about structure. The focus has been placed on the effect of adding either Na(+) or Li(+) to 1-methyl-1-butyl-pyrrolidinium [C4 PYR](+) bis(trifluoromethanesulfonyl)imide [Tf2 N](-) . Radial distribution functions show excellent agreement between MD and AIMD, which ensures the validity of the force fields used in the MD. This is corroborated by the MD results for the density, the diffusion coefficients, and the total conductivity of the electrolytes, which reproduce remarkably well the experimental observations for all studied Na/Li concentrations. By extracting partial conductivities, it is demonstrated that the main contribution to the conductivity is that of [C4 PYR](+) and [Tf2 N](-) . However, addition of Na(+) /Li(+) , although not significant on its own, produces a dramatic decrease in the partial conductivities of the RTIL ions. The origin of this indirect effect can be traced to the modification of the microscopic structure of the liquid as observed from the radial distribution functions, owing to the formation of [Na(Tf2 N)n ]((n-1)-) and [Li(Tf2 N)n ]((n-1)-) clusters at high concentrations. This formation hinders the motion of the large ions, hence reducing the total conductivity. We demonstrate that this clustering effect is common to both Li and Na, showing that both ions behave in a similar manner at a microscopic level in spite of their distinct ionic radii. This is an interesting finding for extending Li-ion and Li-air technologies to their potentially cheaper Na-based counterparts.

Compound Defects in Halide Perovskites: A First-Principles Study of CsPbI3.

Lattice defects affect the long-term stability of halide perovskite solar cells. Whereas simple point defects, i.e., atomic interstitials and vacancies, have been studied in great detail, here we focus on compound defects that are more likely to form under crystal growth conditions, such as compound vacancies or interstitials, and antisites. We identify the most prominent defects in the archetype inorganic perovskite CsPbI3, through first-principles density functional theory (DFT) calculations. We find that under equilibrium conditions at room temperature, the antisite of Pb substituting Cs forms in a concentration comparable to those of the most prominent point defects, whereas the other compound defects are negligible. However, under nonequilibrium thermal and operating conditions, other complexes also become as important as the point defects. Those are the Cs substituting Pb antisite, and, to a lesser extent, the compound vacancies of PbI2 or CsPbI3 units, and the I substituting Cs antisite. These compound defects only lead to shallow or inactive charge carrier traps, which testifies to the electronic stability of the halide perovskites. Under operating conditions with a quasi-Fermi level very close to the valence band, deeper traps can develop.

On the Use of Water and Methanol with Zeolites for Heat Transfer.

Reducing carbon dioxide emissions has become a must in society, making it crucial to find alternatives to supply the energy demand. Adsorption-based cooling and heating technologies are receiving attention for thermal energy storage applications. In this paper, we study the adsorption of polar working fluids in hydrophobic and hydrophilic zeolites by means of experimental quasi-equilibrated temperature-programmed desorption and adsorption combined with Monte Carlo simulations. We measured and computed water and methanol adsorption isobars in high-silica HS-FAU, NaY, and NaX zeolites. We use the experimental adsorption isobars to develop a set of parameters to model the interaction between methanol and the zeolite and cations. Once we have the adsorption of these polar molecules, we use a mathematical model based on the adsorption potential theory of Dubinin-Polanyi to assess the performance of the adsorbate-working fluids for heat storage applications. We found that molecular simulations are an excellent tool for investigating energy storage applications since we can reproduce, complement, and extend experimental observations. Our results highlight the importance of controlling the hydrophilic/hydrophobic nature of the zeolites by changing the Al content to maximize the working conditions of the heat storage device.

Water-Gas Shift Reaction to Capture Carbon Dioxide and Separate Hydrogen on Single-Walled Carbon Nanotubes.

In view of the increasingly severe global warming and ocean acidification caused by CO2 emissions, we report a new procedure, named "reactive separation", to capture CO2. We used advanced Monte Carlo and molecular dynamics methods to simulate the water-gas shift reaction in single-walled carbon nanotubes. We found that (11,11) carbon nanotubes with a diameter of 0.75 nm have the best ability to capture CO2 generated in the water-gas shift reaction. When the feed water-gas ratio is 1:1, the pressure is 3 MPa, and the temperature is 473 K, the storage capacity of CO2 reaches 2.18 mmol/g, the molar fraction of CO2 and H2 inside the carbon nanotube is 0.87 and 0.09, respectively, the conversion of CO in the pore is as high as 97.6%, and the CO2/H2 separation factor is 10.3. Therefore, utilizing the reaction and separation coupling effect of carbon nanotubes to adsorb and store the product CO2 formed in the water-gas shift reaction, while separating the generated clean energy gas H2, is a promising strategy for developing novel CO2 capture technologies.

Atomistic Insights Into the Degradation of Inorganic Halide Perovskite CsPbI3: A Reactive Force Field Molecular Dynamics Study.

Halide perovskites make efficient solar cells but suffer from several stability issues. The characterization of these degradation processes is challenging because of the limited spatiotemporal resolution in experiments and the absence of efficient computational methods to study these reactive processes. Here, we present the first reactive force field for molecular dynamics simulations of the phase instability and the defect-induced degradation in CsPbI3. We find that the phase transitions are driven by the anharmonic fluctuations of the atoms in the perovskite lattice. At low temperatures, the Cs cations tend to move away from their preferential positions, resulting in worse contacts with the surrounding metal halide framework which initiates the conversion to a nonperovskite phase. Moreover, our simulations of defective structures reveal that, although both iodine vacancies and interstitials are mobile in the perovskite lattice, the vacancies have a detrimental effect on the stability, leading to the decomposition of perovskites to PbI2.

Separation of CF4/N2, C2F6/N2, and SF6/N2 Mixtures in Amorphous Activated Carbons Using Molecular Simulations.

The capture and separation of CF4, C2F6, and SF6 and their mixtures containing nitrogen is a challenging process. To solve this, we propose the use of saccharose coke-based carbons as membranes for the adsorption and separation of these gases. By means of advanced techniques of Monte Carlo and molecular dynamics simulations, we have studied the adsorption and diffusion of CF4, C2F6, and SF6 as well as their mixtures with nitrogen in three HRMC carbon models, namely, CS400, CS1000, and CS1000a. We have computed the adsorption isotherms of the single components and the heat of adsorption as a function of the adsorbed concentration. We have also calculated the competitive adsorption of fluoride molecules and nitrogen at two different molar fractions, 0.1 and 0.9. We have computed the transport properties of the adsorbed gases in terms of the self-diffusivities and corrected diffusivities. The performance of the membranes for the targeted separations has been characterized by the calculation of the permselectivity. Our results indicate that the activated amorphous carbon CS1000a is an efficient adsorbent for the capture of the fluoride adsorbates as well as their purification from nitrogen-based mixtures.

Aqueous Co-Solvent in Zwitterionic-based Protic Ionic Liquids as Electrolytes in 2.0†V Supercapacitors.

High-performance energy-storage devices are receiving great interest in sustainable terms as a required complement to renewable energy sources to level out the imbalances between supply and demand. Besides electrode optimization, a primary objective is also the judicious design of high-performance electrolytes combining novel ionic liquids (ILs) and mixtures of aqueous solvents capable of offering "à la carte" properties. Herein, it is described the stoichiometric addition of a zwitterion such as betaine (BET) to protic ILs (PILs) such as those formed between methane sulfonic acid (MSAH) or p-toluenesulfonic acid (PTSAH) with ethanolamine (EOA). This addition resulted in the formation of zwitterionic-based PILs (ZPILs) containing the original anion and cation as well as the zwitterion. The ZPILs prepared in this work ([EOAH]+ [BET][MSA]- and [EOAH]+ [BET][PTSA]- ) were liquid at room temperature even though the original PILs ([EOAH]+ [MSA]- and [EOAH]+ [PTSA]- ) were not. Moreover, ZPILs exhibited a wide electrochemical stability window, up to 3.7 V vs. Ag wire for [EOAH]+ [BET][MSA]- and 4.0 V vs. Ag wire for [EOAH]+ [BET][PTSA]- at room temperature, and a high miscibility with both water and aqueous co-solvent (WcS) mixtures. In particular, "WcS-in-ZPIL" mixtures of [EOAH]+ [BET][MSA]- in 2 H2 O/ACN/DMSO provided specific capacitances of approximately 83 F g-1 at current densities of 1 A g-1 , and capacity retentions of approximately 90 % after 6000 cycles when operating at a voltage of 2.0 V and a current density of 4 A g-1 .

Adsorption and Diffusion of Benzene in Mg-MOF-74 with Open Metal Sites.

We performed molecular simulations to investigate the adsorption and diffusion of benzene in metal-organic framework Mg-MOF-74. At 300 K and 20 Pa, the saturated loading of benzene reaches 8.2 mmol/g, almost twice of (12,12) single-walled carbon nanotube with a similar pore size, and 93% of the benzene molecules in Mg-MOF-74 can desorb at 390 K. The energy analysis indicates that the van der Waals contribution still dominates 70-80% of the total fluid-wall interaction energy compared with the Coulombic contribution. We further analyzed the structure of benzene confined in Mg-MOF-74 by the molecular snapshots, pair correlation functions, orientational order parameters, and local density profiles. It is found that low temperature and high pressure make the structure of adsorbed benzene more similar to that of the liquid benzene. Moreover, the benzene molecules in the contact adsorption layer lie flat on the surface of adsorbent, whereas those molecules near the pore center have no particular orientations. Due to the existence of open metal sites, the structures of adsorbed benzene are more compact and ordered than those of bulk liquid benzene. Consequently, the self-diffusion coefficient of saturated benzene in Mg-MOF-74 at 300 K is significantly lower than that of bulk liquid benzene and confined liquid benzene in slit pores and disordered carbons by 4-5 orders of magnitude. We investigated the separation and diffusion of benzene/cyclohexane in the mixture in Mg-MOF-74 and found that the pores almost completely adsorbed benzene, although its self-diffusion coefficient was slightly lower than that of cyclohexane.

Role of Ionic Liquid [EMIM]+[SCN]- in the Adsorption and Diffusion of Gases in Metal-Organic Frameworks.

We study the adsorption performance of metal-organic frameworks (MOFs) impregnated of ionic liquids (ILs). To this aim we calculated adsorption and diffusion of light gases (CO2, CH4, N2) and their mixtures in hybrid composites using molecular simulations. The hybrid composites consist of 1-ethyl-3-methylimidazolium thiocyanate impregnated in IRMOF-1, HMOF-1, MIL-47, and MOF-1. We found that the increase of the amount of IL enhances the adsorption selectivity in favor of carbon dioxide for the mixtures CO2/CH4 and CO2/N2 and in favor of methane in the mixture CH4/N2. We also provide detailed analysis of the microscopic organization of ILs and adsorbates via radial distribution functions and average occupation profiles and study the impact of the ILs in the diffusion of the adsorbates inside the pores of the MOFs. Based on our findings, we discuss the advantages of using IL/MOF composites for gas adsorption to increase the adsorption of gases and to control the pore sizes of the structures to foster selective adsorption.

Electrochemical Reduction of Oxygen in Aprotic Ionic Liquids Containing Metal Cations: A Case Study on the Na-O2 system.

Metal-air batteries are intensively studied because of their high theoretical energy-storage capability. However, the fundamental science of electrodes, electrolytes, and reaction products still needs to be better understood. In this work, the ionic liquid N-butyl-N-methylpyrrolidinium bis(trifluoromethanesulfonyl)imide (PYR14TFSI) was chosen to study the influence of a wide range of metal cations (Mn+ ) on the electrochemical behavior of oxygen. The relevance of the theory of Lewis hard and soft acids and bases to predict satisfactorily the reduction potential of oxygen in electrolytes containing metal cations is demonstrated. Systems with soft and intermediate Mn+ acidity are shown to facilitate oxygen reduction and metal oxide formation, whereas oxygen reduction is hampered by hard acid cations such as sodium and lithium. Furthermore, DFT calculations on the energy of formation of the resulting metal oxides rationalize the effect of Mn+ on oxygen reduction. A case study on the Na-O2 system is described in detail. Among other things, the Na+ concentration of the electrolyte is shown to control the electrochemical pathway (solution precipitation vs. surface deposition) by which the discharge product grows. All in all, fundamental insights for the design of advanced electrolytes for metal-air batteries, and Na-air batteries in particular, are provided.