Theoretical investigation of (La4O6)n, (La2Ce2O7)n, and (Ce4O8)n nanoclusters (n = 10, 18): Temperature effects and O-vacancy formation.

We report a theoretical investigation of temperature, size, and composition effects on the structural, energetic, and electronic properties of the (La4O6)n, (La2Ce2O7)n, and (Ce4O8)n nanoclusters (NCs) for n = 10, 18. Furthermore, we investigated the single O vacancy formation energy as a function of the geometric location within the NC. Our calculations are based on the combination of force-field molecular dynamics (MD) simulations and density functional theory calculations. We identified a phase transition from disordered to ordered structures for all NCs via MD simulations and structural analysis, e.g., radius changes, radial distribution function, common neighbor analysis, etc. The transition is sharp for La36Ce36O126, La20Ce20O70, and Ce72O144 due to the crystalline domains in the core and less abrupt for Ce40O80, La40O60, and La72O108. As expected, radius changes are abrupt at the transition temperature, as are morphological differences between NCs located below and above the transition temperature. We found a strong dependence on the O vacancy formation energy (Evac) and its location within the NCs. For example, for La40O60, Evac decreases almost linearly as the distance from the geometric center increases; however, the same trend was not observed for Ce40O80, while there are large deviations from the linear trend for La20Ce20O70. Evac has smaller values for Ce40O80 and higher values for La40O60, that is, almost three times, while Evac has intermediate values for mixed oxides, as expected from weighted averages. Therefore, the mixture of one formula unit of La2O3 with two formula units of CeO2 has the effect of increasing the stability of CeO2 (binding energy), which increases the magnitude of the formation energy of the O vacancy.

Direct Evidence of Reversible Changes in Electrolyte and its Interplay with LiO2 Intermediate in Li-O2 Batteries.

Lithium-oxygen batteries show promising energy storage potential with high theoretical energy density; however, further investigation of chemical reactions is required. In this study, experimental Raman and theoretical analyzes are performed for a Li-O2 battery with LiClO4/dimethyl sulfoxide (DMSO) electrolyte and carbon cathode to understand the role of intermediate species in the reactional mechanism of the cell using a high donor number solvent. Operando Raman results reveal reversible changes in the DMSO bands, in addition to the formation and decomposition of Li2O2. On discharge, a decrease in DMSO polarizability is observed and bands of DMSO-Li+-anion interactions are evidenced and supported by ab initio density functional theory (DFT) calculations. Molecular dynamics (MD) force field simulations and operando Raman show that DMSO interacts with LiO2(sol), highlighting the stability of the electrolyte compared to the interaction with reactive O 2 - ${\rm O}_2^{-}$ . On charging, the presence of Li+ indicates the formation of a lithium-deficient phase, followed by the release of Li+ and oxygen. Therefore, this study contributes to understanding the discharge/charge chemistry of a Li-O2 cell, employing a common carbon cathode and DMSO electrolyte. The combination of a simple characterization technique in operando mode and theoretical studies provides essential information on the mechanism of Li-O2 system.

Theoretical Framework Based on Molecular Dynamics and Data Mining Analyses for the Study of Potential Energy Surfaces of Finite-Size Particles.

Nanoclusters are remarkably promising for the capture and activation of small molecules for fuel production or as precursors for other chemicals of high commercial value. Since this process occurs under a wide variety of experimental conditions, an improved atomistic understanding of the stability and phase transitions of these systems will be key to the development of successful technological applications. In this work, we proposed a theoretical framework to explore the potential energy surface and configuration space of nanoclusters to map the most important morphologies presented by those systems and the phase transitions between them. A fully automated process was developed, which combines global optimization techniques, classical molecular dynamics, and unsupervised machine learning algorithms. To showcase these capabilities of the approach, we explored the example of copper nanoclusters (Cun) where n = 13, 38, 55, 75, 98, 102, and 147. We not only reported a graphical potential energy surface for each size, but also explored the topology of the configuration space via structural and thermodynamic analyses. The effect of size on the potential energy surface and the critical temperature for solid-liquid phase transitions were also reported, highlighting the impact of magic numbers on those quantities.

Screening of the Role of the Chemical Structure in the Electrochemical Stability Window of Ionic Liquids: DFT Calculations Combined with Data Mining.

Ionic liquids have attracted the attention of researchers as possible electrolytes for electrochemical energy storage devices. However, their properties, such as the electrochemical stability window (ESW), ionic conductivity, and diffusivity, are influenced both by the chemical structures of cations and anions and by their combinations. Most studies in the literature focus on the understanding of common ionic liquids, and little effort has been made to find ways to improve our atomistic understanding of those systems. The goal of this paper is to explore the structural characteristics of cations and anions that form ionic liquids that can expand the HOMO/LUMO gap, a property directly linked to the ESW of the electrolyte. For that, we design a framework for randomly generating new ions by combining their fragments. Within this framework, we generate about 104 cations and 104 anions and fully optimize their structures using density functional theory. Our calculations show that aromatic cations are less stable ionic liquids than aliphatic ones, an expected result if chemical rationale is used. More importantly, we can improve the gap by adding electron-donating and electron-withdrawing functional groups to the cations and anions, respectively. The increase can be about 2 V, depending on the case. This improvement is reflected in a wider ESW.

Steric and Electrostatic Effects on the Diffusion of CH4/CH3OH in Copper-Exchanged Zeolites: Insights from Enhanced Sampling Molecular Dynamics and Free Energy Calculations.

Copper-exchanged zeolites have demonstrated high selectivity in methane-to-methanol conversion carried out on copper-oxo centers. Nevertheless, the reaction can only occur if the methane molecules reach the active site while the methanol molecules must leave the material without high energetic cost for the migration. In this context, we have used force field-based molecular dynamics simulations with the potential of mean force method to estimate the energy barrier in cage to cage diffusion of methane and methanol molecules in the chabazite framework type zeolite. The results show considerably higher energy barrier for methanol diffusion. The steric effect of the active site and the electrostatic environment favors the CH3OH diffusion toward nonactive cages where it tends to accumulate due to the strong interactions with the zeolite. The same behavior is observed in the water molecules distribution, which emphasizes the control of the electrostatic potential over the polar molecules migration. For high concentration of polar molecules, the electrostatic effect is shielded and the driving force is reduced for CH3OH diffusion. The results show that if the electrostatic environment can be controlled, the product migration may be facilitated, which can improve the catalytic process.

Interfacial Structures in Ionic Liquid-Based Ternary Electrolytes for Lithium-Metal Batteries: A Molecular Dynamics Study.

Lithium-metal batteries are promising candidates to fulfill the future performance requirements for energy storage applications. However, the tendency to form metallic dendrites and the undesirable side reactions between the electrolyte and the Li electrode lead to poor performance and safety issues in these batteries. Therefore, understanding the interfacial properties and the Li-metal surface/electrolyte interactions is crucial to resolve the remaining obstacles and make these devices feasible. Here, we report a computational study on the interface effects in ternary polymer electrolytes composed by poly(ethylene oxide) (PEO), lithium salts, and different ionic liquids (ILs) confined between two Li-metal slabs. Atomistic simulations are used to characterize the local environment of the Li+ ions and the transport properties in the bulk and at the interface regions. Aggregation of ions at the metal surface is seen in all investigated systems; the structure and composition are directly correlated to the IL components. The strong interactions between the electrolyte species and the Li-metal atoms result in the structuring of the electrolyte at the interface region, in which comparatively small and flat ions result in a well-defined region with extensive Li+ populations and high self-diffusion coefficients. In contrast, large ions such as [P222mom]+ increase the PEO density in the bulk due to large steric effects at the interface. Therefore, the choice of specific ILs in ternary polymer electrolytes can tune the structure-dynamic properties at the Li-metal surface/electrolyte interface, controlling the SEI formation at the electrode surface, and thereby improve battery performance.

Solvation Structure and Dynamics of Li+ in Ternary Ionic Liquid-Lithium Salt Electrolytes.

The structural and dynamical changes in the solvation shell surrounding Li+ in a multianion environment are studied by Raman spectroscopy and molecular dynamics (MD) simulations. The ternary electrolyte is composed of a mixture of two ionic liquids (ILs), n-methyl- n-propylpyrrolidinium bis(trifluoromethanesulfonyl)imide ([PYR13][TFSI]) and 1-ethyl-3-methylimidazolium dicyanamide ([EMIM][DCA]), and a lithium bis(trifluoromethanesulfonyl)imide ([Li][TFSI]) salt (0-1 M). A 1:9 volumetric mixture of [PYR13][TFSI]/[EMIM][DCA] formed an eutectic that exhibited a lower melting point than that of either parent IL. The local structure of Li+ in this eutectic is found to be heterogenous and preferentially solvated by [DCA], which is the smaller and more abundant anion. Whereas [TFSI] is able to bridge multiple Li+ at high salt concentrations and form both monodentate and bidentate conformations through its oxygen atoms, [DCA] is capable of forming only monodentate coordination with Li+ through either of its end nitrogen atoms. The Raman and MD analyses suggest a wide distribution of solvation structures in the form of [Li(TFSI) m(DCA) n]( m+ n-1)- where m = 0-1 and n = 3-4. The computations showed increased ion pair lifetime for Li+-[DCA] and decreased lifetimes for Li+-[TFSI] in the ternary mixture with the increase in the [Li][TFSI] concentration. These results show that the solvation and transport properties of charge carriers in ILs can be modified via the presence of multiple ions with varying degree of coordination, which provides an approach to impact the performance in electrochemical processes.

A molecular dynamics study of lithium-containing aprotic heterocyclic ionic liquid electrolytes.

Classical molecular dynamics simulations were performed on twelve different ionic liquids containing aprotic heterocyclic anions doped with Li+. These ionic liquids have been shown to be promising electrolytes for lithium ion batteries. Self-diffusivities, lithium transference numbers, densities, and free volumes were computed as a function of lithium concentration. The dynamics and free volume decreased with increasing lithium concentration, and the trends were rationalized by examining the changes to the liquid structure. Of those examined in the present work, it was found that (methyloxymethyl)triethylphosphonium triazolide ionic liquids have the overall best performance.

Local environment structure and dynamics of CO2 in the 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide and related ionic liquids.

Despite the innumerous papers regarding the study of the ionic liquids as a potential candidate for CO2 capture, many details concerning the structure and dynamics of CO2 in the system are still to be revealed, i.e., the correlation between the local environment structure and the dynamic properties of the substance. This present work relied on the performance of molecular dynamics both for the neat [C2mim][Tf2N] and [C2mim][Tf2N]/CO2 mixtures in an attempt to elucidate the local environment of CO2 and their effects on the dynamic properties of [C2mim][Tf2N]. A slight change in the orientation of the cation and anion could be observed, which was correlated to the cation and anion moving away from each other in order to receive the carbon dioxide. The gas molecules pushed both the cation and the anion away to create sufficient void to its accommodation. The diffusion coefficient of [C2mim]+ is higher than [Tf2N]- regardless the increase of the CO2 concentration. The addition of CO2 in the ionic liquid has shown an increase of 4-5 times for the diffusivity of ions, which was related to the decrease of cation-anion interaction strength. The transport properties' results showed that the addition of CO2 in the ionic liquid generates the fluidization of the system, decreasing the viscosity as a consequence of the local environment structure changing. Likewise, the effect of the type of anion and cation on the system properties was studied considering [Ac]- and [BMpyr]+ ions, showing large effects by the change of anion to [Ac]- which rise from the strong [C2mim]+-[Ac]- interaction, which conditions the solvation of ions by CO2 molecules.

Understanding the Solubility of Acetaminophen in 1-n-Alkyl-3-methylimidazolium-Based Ionic Liquids Using Molecular Simulation.

During the manufacturing of pharmaceutical compounds, solvent mixtures are commonly used, where the addition of a cosolvent allows for the tuning of the intermolecular interactions present in the system. Here we demonstrate how a similar effect can be accomplished using a room temperature ionic liquid. The pharmaceutical compound acetaminophen is studied in 21 common ionic liquids composed of a 1-n-alkyl-3-methylimidazolium cation with 1 of 7 anions. Using the acetate anion, we predict a large enhancement in solubility of acetaminophen relative to water. We show how this is caused by a synergistic effect of favorable interactions between the ionic liquid and the phenyl, hydroxyl and amide groups of acetaminophen, demonstrating how the ionic liquid cation and anion may be chosen to preferentially solvate different functional groups of complex pharmaceutical compounds. Additionally, while the use of charge scaling in ionic liquid force fields has previously been found to have a minute effect on ionic liquid structural properties, we find it appreciably affects the computed solvation free energy of acetaminophen, which in turn affects the predicted solubility.

Insights on the solubility of CO2 in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide from the microscopic point of view.

Emissions of greenhouse gases due to human activities have been well documented as well as the effects on global warming resulting from it. Efforts to reduce greenhouse gases at the source are crucial to curb climate change, but due to insignificant economic incentives to reduce usage of fossil fuels, not a lot of progress has been made by this route. This necessitates additional measures to reduce the occurrence of greenhouse gases in the atmosphere. Here we used theoretical methods to study the solubility of carbon dioxide in ionic liquids (ILs) since sequestration of CO2 in ILs has been proposed as a possible technology for reducing the emissions of CO2 to the atmosphere. Ionic liquids form a class of solvents with melting temperatures below 100 °C and, due to very low vapor pressures, which are not volatile. We have performed molecular dynamics (MD) simulations of 1-ethyl-3-methylimidazolium (C2mim) bis(trifluoromethylsulfonyl)imide (Tf2N) and its mixtures with carbon dioxide in order to investigate the CO2 concentration effect on the CO2-cation and CO2-anion interactions. A systematic investigation of CO2 concentration effects on resulting equilibrium liquid structure, and the local environment of the ions is provided. The Quantum Theory of Atoms in Molecules (QTAIM) was used to determine the interaction energy for CO2-cation and CO2-anion complexes from uncorrelated structures derived from MD simulations. A spatial distribution function analysis demonstrates the specific interactions between CO2 and the ionic liquid. Our findings indicate that the total volume of the system increases with the CO2 concentration, with a molar volume of CO2 of about 0.038 L/mol, corresponding to liquid CO2 under a pressure of 100 bar. In other words, the IL effectively pressurizes the CO2 inside its matrix. The thermodynamics of CO2 solvation in C2 min-Tf2N were computed using free energy techniques, and the solubility of CO2 is found to be higher in this IL (-3.7 ± 1 kcal/mol) than in water (+0.2 kJ/mol), predominantly due to anion-CO2 interactions.