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.

Fitting electron density as a physically sound basis for the development of interatomic potentials of complex alloys.

The development of new interatomic potentials to model metallic systems is a difficult task, due in part to the dependence between the parameters that describe the electron density and the short-range interactions. Parameter search methods are prone to false convergence. To solve this problem, we have developed a methodology for obtaining the electron density parameters independently of the short-range interactions, so that physically sound parameters can be obtained to describe the electron density, after which the short-range parameters can be fitted, thus reducing the complexity of the process and yielding better interatomic potentials. With the new method we can develop self-consistent, accurate force fields, using solely calculations, without the need to fit to experimental data. Density functional theory calculations are used to compute the observables with which the potential is fit. We applied the method to a Ni-based Inconel 625 superalloy (IN625), modelled here as Ni, Cr, Mo and Fe solid solution alloys. The capability of the force fields developed using this new method is validated, by comparing the structural and thermo-elastic properties predicted with the force fields, with the corresponding experimental data, both for single crystals and polycrystalline alloys.

On the ionophoric selectivity of nonactin and related macrotetrolide derivatives.

Nonactin and its analogs constitute a central class of macrocycles with an antibiotic activity closely related to their selective ionophoric behavior. In this study, we apply experimental and computational methods to revisit the specificity of cation binding and transport by three nactin variants differing in structural properties, such as the position of the ester linkages, the nature of the side groups, or the flexibility of the backbone. On the one hand, electrospray ionization mass spectrometry and infrared spectroscopy are employed to expose the selectivity of the liquid-liquid (water-chloroform) extraction of alkali cations by nonactin and to demonstrate that the cation complexes are partially hydrated in the organic phase. Furthermore, laser desorption mass spectrometry is employed to determine the intrinsic cation affinities of nonactin under solvent-free conditions. On the other hand, density functional theory calculations are performed to characterize the conformations of the alkali cation complexes of the three nactins, and to assess the role of intermolecular and solvent interactions in determining their relative stability. Depending on the structure of the macrocycle, the cation complexes adopt either a cage-like conformation or a tweezer-like conformation. The computations show that the partial hydration of those different conformations in the organic phase, determine the distinct cation extraction selectivities that are observed experimentally.

Critical Role of Dynamic Flexibility in Ge-Containing Zeolites: Impact on Diffusion.

Incorporation of germanium in zeolites is well known to confer static flexibility to their framework, by stabilizing the formation of small rings. In this work, we show that the flexibility associated to Ge atoms in zeolites goes beyond this static effect, manifesting also a clear dynamic nature, in the sense that it leads to enhanced molecular diffusion. Our study combines experimental and theoretical methods providing evidence for this effect, which has not been described previously, as well as a rationalization for it, based on atomistic grounds. We have used both pure-silica and silico-germanate ITQ-29 (LTA topology) zeolites as a case study. Based on our simulations, we identify the flexibility associated to the pore breathing-like behavior induced by the Ge atoms, as the key factor leading to the enhanced diffusion observed experimentally in Ge-containing zeolites.

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.

Modelling a Linker Mix-and-Match Approach for Controlling the Optical Excitation Gaps and Band Alignment of Zeolitic Imidazolate Frameworks.

Tuning the electronic structure of metal-organic frameworks is the key to extending their functionality to the photocatalytic conversion of absorbed gases. Herein we discuss how the band edge positions in zeolitic imidazolate frameworks (ZIFs) can be tuned by mixing different imidazole-based linkers within the same structure. We present the band alignment for a number of known and hypothetical Zn-based ZIFs with respect to the vacuum level. Structures with a single type of linker exhibit relatively wide band gaps; however, by mixing linkers of a low-lying conduction edge with linkers of a high-lying valence edge, we can predict materials with ideal band positions for visible-light water splitting and CO2 reduction photocatalysis. By introducing copper in the tetrahedral position of the mixed-linker ZIFs, it would be possible to increase both photo-absorption and the electron-hole recombination times.

Binding Selectivity of Macrocycle Ionophores in Ionic Liquids versus Aqueous Solution and Solvent-free Conditions.

The understanding of supramolecular recognition in room-temperature ionic liquids (RTILs) is key to develop the full potential of these materials. In this work, we provide insights into the selectivity of the binding of alkali metal cations by standard cyclodextrin and calixarene macrocycles in RTILs. A direct laser desorption/ionization mass spectrometry approach is employed to determine the relative abundances of the inclusion complexes formed through competitive binding in RTIL solutions. The results are compared with the binding selectivities measured under solvent-free conditions and in water/methanol solutions. Cyclodextrins and calixarenes in which the peripheral OH groups are substituted by bulkier side groups preferentially bind to Cs(+) . Such specific ionophoric behavior is substantially enhanced by solvation effects in the RTIL. This finding is rationalized with the aid of quantum mechanical calculations, in terms of the conformational features and steric interactions that drive the solvation of the inclusion complexes by the bulky RTIL counterions.

Ultraviolet Pretreatment of Titanium Dioxide and Tin-Doped Indium Oxide Surfaces as a Promoter of the Adsorption of Organic Molecules in Dry Deposition Processes: Light Patterning of Organic Nanowires.

In this article we present the preactivation of TiO2 and ITO by UV irradiation under ambient conditions as a tool to enhance the incorporation of organic molecules on these oxides by evaporation at low pressures. The deposition of π-stacked molecules on TiO2 and ITO at controlled substrate temperature and in the presence of Ar is thoroughly followed by SEM, UV-vis, XRD, RBS, and photoluminescence spectroscopy, and the effect is exploited for the patterning formation of small-molecule organic nanowires (ONWs). X-ray photoelectron spectroscopy (XPS) in situ experiments and molecular dynamics simulations add critical information to fully elucidate the mechanism behind the increase in the number of adsorption centers for the organic molecules. Finally, the formation of hybrid organic/inorganic semiconductors is also explored as a result of the controlled vacuum sublimation of organic molecules on the open thin film microstructure of mesoporous TiO2.

Thermostructural behaviour of Ni-Cr materials: modelling of bulk and nanoparticle systems.

The thermostructural properties of Ni-Cr materials, as bulk and nanoparticle (NP) systems, have been predicted with a newly developed interatomic potential, for Ni/Cr ratios from 100/0 to 60/40. The potential, which has been fitted using experimental data and further validated using Density Functional Theory (DFT), describes correctly the variation with temperature of lattice parameters and the coefficient of thermal expansion, from 100 K to 1000 K. Using this potential, we have performed Molecular Dynamics (MD) simulations on bulk Ni-Cr alloys of various compositions, for which no experimental data are available. Similarly, NPs with diameters of 3, 5, 7, and 10 nm were studied. We found a very rapid convergence of NP properties with the size of the systems, showing already the 5 nm NPs with a thermostructural behaviour similar to the bulk. MD simulations of two 5 nm NPs show very little sintering and thermally induced damage, for temperatures between 300 K and 1000 K, suggesting that materials formed by agglomeration of Ni-Cr NPs meet the thermostructural stability requirements for catalysis applications.

Nucleation of zeolitic imidazolate frameworks: from molecules to nanoparticles.

We have studied the clusters involved in the initial stages of nucleation of Zeolitic Imidazolate Frameworks, employing a wide range of computational techniques. In the pre-nucleating solution, the prevalent cluster is the ZnIm4 cluster (formed by a zinc cation, Zn2+, and four imidazolate anions, Im-), although clusters such as ZnIm3, Zn2Im7, Zn2Im7, Zn3Im9, Zn3Im10, or Zn4Im12 have energies that are not much higher, so they would also be present in solution at appreciable quantities. All these species, except ZnIm3, have a tetrahedrally coordinated Zn2+ cation. Small ZnxImy clusters are less stable than the ZnIm4 cluster. The first cluster that is found to be more stable than ZnIm4 is the Zn41Im88 cluster, which is a disordered cluster with glassy structure. Bulk-like clusters do not begin to be more stable than glassy clusters until much larger sizes, since the larger cluster we have studied (Zn144Im288) is still less stable than the glassy Zn41Im88 cluster, suggesting that Ostwald's rule (the less stable polymorph crystallizes first) could be fulfilled, not for kinetic, but for thermodynamic reasons. Our results suggest that the first clusters formed in the nucleation process would be glassy clusters, which then undergo transformation to any of the various crystal structures possible, depending on the kinetic routes provided by the synthesis conditions. Our study helps elucidate the way in which the various species present in solution interact, leading to nucleation and crystal growth.

Multipodal coordination of a tetracarboxylic crown ether with NH4(+): a vibrational spectroscopy and computational study.

The elucidation of the structural requirements for molecular recognition by the crown ether (18-crown-6)-2,3,11,12-tetracarboxylic acid (18c6H(4)) and its cationic complexes constitutes a topic of current fundamental and practical interest in catalysis and analytical sciences. The flexibility of the central ether ring and its four carboxyl side arms poses important challenges to experimental and theoretical approaches. In this study, infrared action vibrational spectroscopy and quantum mechanical computations are employed to characterize the conformational structure of the isolated gas phase complex formed by the 18c6H(4) host with NH(4)(+) as guest. The results show that the most stable gas-phase structure is a barrel-like conformation sustained by tetrapodal H-bonding of the ammonia cation with two C=O side groups and with four oxygen atoms of the ether ring in a bifurcated arrangement. Interestingly, a similar structure had been proposed in previous crystallographic studies. The experiment also provides evidence for a significant contribution of a higher energy bowl-like conformer with features resembling those adopted by 18c6H(4) in the analogous complexes with secondary amines. Such a conformation displays H-bonding between confronted side carboxyl groups and tetrapodal binding of the NH(4)(+) with the ether ring and with one C=O group. Structures involving even more extensive intramolecular H-bonding in the 18c6H(4) substrate are found to lie higher in energy and are ruled out by the experiment.

Incorporating zeolitic-imidazolate framework-8 nanoparticles into kidney scaffolds: a first step towards innovative renal therapies.

We demonstrate for the first time the potential of zeolitic-imidazolate framework-8 nanoparticles to be incorporated within a renal scaffold while retaining their ability to remove uremic toxins (mainly hydrophobic toxins like p-cresol) under flow conditions. This work may pave the way for the future development of novel adsorbents for dialysis and/or artificial kidneys.

Unravelling the key factors in the chlorine-promoted epoxidation of ethylene over a silver-copper oxide nanocatalyst.

Ethylene oxide is one of the most important raw materials in the chemical industry, with an annual production close to 35 million metric tons. Despite its importance, to date, no metal has been found that can compete with the original silver bulk material catalyst discovered in 1931. Recently, a few copper and copper-silver based nanostructures have demonstrated remarkable selectivity and activity, especially when coupled with an industrial chlorine promoter. The present work evaluates the mechanistic role of chlorine as an active promoter of the selective oxidation of ethylene to ethylene oxide in the presence of a silver-copper oxide hybrid nanocatalyst (AgCuO). Experimental kinetic studies combined with density functional theory (DFT) calculations provide insight into the influence that Ag/CuO-supported chlorine atoms have over the ethylene epoxidation reaction. Remarkably, the typically described indirect route via the formation of an oxametallacycle (OMC) is also accompanied by a direct route. Furthermore, the presence of chlorine seems to facilitate a more favorable adsorption energy for ethylene oxide (EO) than for acetaldehyde (AA), the main reaction by-product. As a result, complete oxidation of EO can be further prevented in the presence of this AgCuO hybrid heteronanostructure.

Effect of air humidity on the removal of carbon tetrachloride from air using Cu-BTC metal-organic framework.

We have used interatomic potential-based simulations to study the removal of carbon tetrachloride from air at 298 K, using Cu-BTC metal organic framework. We have developed new sets of Lennard-Jones parameters that accurately describe the vapour-liquid equilibrium curves of carbon tetrachloride and the main components from air (oxygen, nitrogen, and argon). Using these parameters we performed Monte Carlo simulations for the following systems: (a) single component adsorption of carbon tetrachloride, oxygen, nitrogen, and argon molecules, (b) binary Ar/CCl(4), O(2)/CCl(4), and N(2)/CCl(4) mixtures with bulk gas compositions 99 : 1 and 99.9 : 0.1, (c) ternary O(2)/N(2)/Ar mixtures with both, equimolar and 21 : 78 : 1 bulk gas composition, (d) quaternary mixture formed by 0.1% of CCl(4) pollutant, 20.979% O(2), 77.922% N(2), and 0.999% Ar, and (e) five-component mixtures corresponding to 0.1% of CCl(4) pollutant in air with relative humidity ranging from 0 to 100%. The carbon tetrachloride adsorption selectivity and the self-diffusivity and preferential sitting of the different molecules in the structure are studied for all the systems.

Crown ether complexes with H3O⁺ and NH4⁺: proton localization and proton bridge formation.

The complexes formed by crown ethers with hydronium and ammonium cations are of key relevance for the understanding of their supramolecular behavior in protic solvents. In this work, the complexes of the 15-crown-5 (15c5) and 18-crown-6 (18c6) ethers with H3O⁺ and NH4⁺ and their deuterated variants are investigated under isolated conditions. The study employs infrared multiple photon dissociation (IRMPD) vibrational spectroscopy and DFT B3LYP/6-31++G(d,p) calculations for conformational assignment. The 18c6 ether provides two energetically nearby C(3v) conformations with commensurate linear O-H···O and N-H···O bonds. The 15c5 ether ring adopts partially folded asymmetric pyramidal geometries, yielding one shorter linear H bond and two longer non-linear H bonds. Remarkably, an appreciable broadening of the IRMPD vibrational bands is observed for the 15c5-H3O⁺/D3O⁺ complexes. This can be interpreted as a signature for partial sharing of the proton (or deuteron) between the water and the crown ether along the linear O-H···O intermolecular H bond, which is indeed particularly short for this complex.

Bandgap Engineering in the Configurational Space of Solid Solutions via Machine Learning: (Mg,Zn)O Case Study.

Computer simulations of alloys' properties often require calculations in a large space of configurations in a supercell of the crystal structure. A common approach is to map density functional theory results into a simplified interaction model using so-called cluster expansions, which are linear on the cluster correlation functions. Alternative descriptors have not been sufficiently explored so far. We show here that a simple descriptor based on the Coulomb matrix eigenspectrum clearly outperforms the cluster expansion for both total energy and bandgap energy predictions in the configurational space of a MgO-ZnO solid solution, a prototypical oxide alloy for bandgap engineering. Bandgap predictions can be further improved by introducing non-linearity via gradient-boosted decision trees or neural networks based on the Coulomb matrix descriptor.

Exploration of multiple energy landscapes for zirconia nanoclusters.

We have predicted the stable and low-energy metastable structures for (ZrO(2))(n) clusters, where n = 1 to 12, employing Density Functional Theory (DFT) at the PBEsol0 level. A process of data mining and the application of an evolutionary algorithm to three different energy landscapes, as defined by interatomic potentials, for each cluster size, was used to generated the plausible structures for refinement using DFT at the PBEsol level. The structures for zirconia were found to be similar to that predicted for titania except that the order, with respect to the binding energies, of the configurations for the two compounds were different for the larger sized clusters. The energies of the local minima configurations and their respective HOMO and LUMO are discussed, as well as the insights gained from employing more than one set of interatomic potentials in the initial global optimisations.

Gas-phase complexes of cyclic and linear polyethers with alkali cations.

The flexibility of polymer backbones constitutes one key aspect in their molecular recognition properties. This investigation characterizes the structure of the gas-phase complexes formed by the cyclic and linear polyethers with the heavier alkali metal cations. In particular, the cyclic 15-crown-5 ether (15c5), (OCH(2)CH(2))(5), and the polyethylene glycol linear chains PEG4 and PEG9, H(OCH(2)CH(2))(n=4,9)OH, are considered. Infrared multiple photon dissociation spectroscopy is applied to probe the polymer vibrational modes within the spectral range 800-1500 cm(-1), in combination with Density Functional Theory calculations. The experimental spectra of the 15c5-M(+) (M = K, Rb, Cs) complexes correlate with distorted asymmetric backbone structures in which the cation lies above the ether ring, and four oxygens are oriented toward the cation. In contrast, the PEG4-K(+) complex features an inclusion-like fivefold coordination structure in which the cation and four oxygens are quasi-coplanar and one terminal oxygen lies out of plane. For the PEG9-K(+) complex, the ether chain builds stable cages involving a coordination of eight oxygens with the cation, sustained by hydrogen bonds between the terminal hydroxyl groups.

Modelling nano-clusters and nucleation.

We review the growing role of computational techniques in modelling the structures and properties of nano-particulate oxides and sulphides. We describe the main methods employed, including those based on both electronic structure and interatomic potential approaches. Particular attention is paid to the techniques used in searching for global minima in the energy landscape defined by the nano-particle cluster. We summarise applications to the widely studied ZnO and ZnS systems, to silica nanochemistry and to group IV oxides including TiO(2). We also consider the special case of silica cluster chemistry in solution and its importance in understanding the hydrothermal synthesis of microporous materials. The work summarised, together with related experimental studies, demonstrates a rich and varied nano-cluster chemistry for these materials.

Emergence of symmetry and chirality in crown ether complexes with alkali metal cations.

Crown ethers provide a valuable benchmark for the comprehension of molecular recognition mediated by inclusion complexes. One of the most relevant crown ethers, 18-crown-6 (18c6), features a flexible six-oxygen cyclic backbone that is well-known for its selective cation binding. This study employs infrared spectroscopy and quantum mechanical calculations to elucidate the structure of the gas-phase complexes formed by the 18c6 ether with the alkali metal cations. It is shown that symmetric and chiral arrangements play a dominant role in the conformational landscape of the 18c6-alkali system. Most stable 18c6-M(+) conformers are found to have symmetries C(3v) and C(2) for Cs(+), D(3d) for K(+), C(1) and D(3d) for Na(+), and D(2) for Li(+). Remarkably, whereas the bare 18c6 ether is achiral, chirality emerges in the C(2) and D(2) 18c6-M(+) conformations, both of which involve pairs of stable atropoisomers capable of acting as enantiomeric selective substrates.