Lithium solvation in dimethyl sulfoxide-acetonitrile mixtures

We present molecular dynamics simulation results pertaining to the solvation of Li+ in dimethyl sulfoxide-acetonitrile binary mixtures. The results are potentially relevant in the design of Li-air batteries that rely on aprotic mixtures as solvent media. To analyze effects derived from differences in ionic size and charge sign, the solvation of Li+ is compared to the ones observed for infinitely diluted K+ and Cl− species, in similar solutions. At all compositions, the cations are preferentially solvated by dimethyl sulfoxide. Contrasting, the first solvation shell of Cl− shows a gradual modification in its composition, which varies linearly with the global concentrations of the two solvents in the mixtures. Moreover, the energetics of the solvation, described in terms of the corresponding solute-solvent coupling, presents a clear non-ideal concentration dependence. Similar nonlinear trends were found for the stabilization of different ionic species in solution, compared to the ones exhibited by their electrically neutral counterparts. These tendencies account for the characteristics of the free energy associated to the stabilization of Li+Cl−, contact-ion-pairs in these solutions. Ionic transport is also analyzed. Dynamical results show concentration trends similar to those recently obtained from direct experimental measurements.Fil: Semino, Rocio. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; ArgentinaFil: Zaldivar, Gervasio. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; ArgentinaFil: Calvo, Ernesto Julio. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; ArgentinaFil: Laria, Daniel Hector. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina. Comisión Nacional de Energía Atómica; Argentin

Combined Effects from Solvation and Nuclear Quantum Fluctuations on Autoionization Mechanisms in Aqueous Clusters

Using path-integral molecular dynamics simulations, we examine isomerization paths involving collective proton transfers in [H2O]5 and [H2O]8 clusters under cryogenic conditions. We focused attention on combined effects derived from solvation and nuclear quantum fluctuations on the characteristics of free energy barriers and relative stabilities of reactants and products. In particular, we analyzed two different processes: the first one involves the exchange of donor-acceptor hydrogen bond roles along cyclic moieties, whereas the second one corresponds to charge separation leading to stable [H3O]+[OH]- ion pairs. In the first case, the explicit incorporation of quantum tunneling introduces important modifications in the classical free energy profile. The resulting quantum profile presents two main contributions, one corresponding to compressions of O-O distances and a second one ascribed to nuclear tunneling of the light protons. Solvation effects promote a moderate polarization of the cyclic structures and a partial loss of concertedness in the collective modes, most notably, at the onset of tunneling. Still, the latter effects are also sufficiently strong to promote the stabilization of ion pairs along the classical trajectories. In contrast, the explicit incorporation of nuclear quantum fluctuations leads to charge separated configurations that are marginally stable. As such, the latter states could also be regarded as short-lived intermediate states along the reactive exchange path.Fil: Turi, Laszlo. Eötvös Loránd University; HungríaFil: Rodriguez, Javier. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Unidad Ejecutora Instituto de Nanociencia y Nanotecnología. Unidad Ejecutora Instituto de Nanociencia y Nanotecnología - Nodo Constituyentes | Comisión Nacional de Energía Atómica. Unidad Ejecutora Instituto de Nanociencia y Nanotecnología. Unidad Ejecutora Instituto de Nanociencia y Nanotecnología - Nodo Constituyentes; ArgentinaFil: Laria, Daniel Hector. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentin

Nuclear quantum effects on the structure and the dynamics of [H2O]8 at low temperatures

We use ring-polymer-molecular-dynamics (RPMD) techniques and the semi-empirical q-TIP4P/F water model to investigate the relationship between hydrogen bond connectivity and the charac- teristics of nuclear position fluctuations, including explicit incorporation of quantum effects, for the energetically low lying isomers of the prototype cluster [H2O]8 at T = 50 K and at 150 K. Our results reveal that tunneling and zero-point energy effects lead to sensible increments in the magnitudes of the fluctuations of intra and intermolecular distances. The degree of proton spatial delocalization is found to map logically with the hydrogen-bond connectivity pattern of the cluster. Dangling hydro- gen bonds exhibit the largest extent of spatial delocalization and participate in shorter intramolecular O-H bonds. Combined effects from quantum and polarization fluctuations on the resulting individ- ual dipole moments are also examined. From the dynamical side, we analyze the characteristics of the infrared absorption spectrum. The incorporation of nuclear quantum fluctuations promotes red shifts and sensible broadening relative to the classical profile, bringing the simulation results in much more satisfactory agreement with direct experimental information in the mid and high fre- quency range of the stretching band. While RPMD predictions overestimate the peak position of the low frequency shoulder, the overall agreement with that reported using an accurate, parame- terized, many-body potential is reasonable, and far superior to that one obtains by implementing a partially adiabatic centroid molecular dynamics approach. Quantum effects on the collective dynam- ics, as reported by instantaneous normal modes, are also discussedFil: Videla, Pablo Ernesto. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina. Universidad de Buenos Aires; ArgentinaFil: Rossky, Peter J.. University of Texas at Austin; Estados UnidosFil: Laria, Daniel Hector. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina. Comision Nacional de Energia Atomica. Gerencia Quimica. CAC; Argentin

Building up a nanotube in solution: docking of Janus-Cyclodextrin

Using molecular dynamics experiments, we analyze the association of "Janus" 6-amino-6-deoxy-2O-carboxymethyl-β-cyclodextrins (JCD) in aqueous solutions. In JCD dimers, the free energy associated with the primary-rim-secondary-rim docking shows a stable minimum of ∼-45 kcal mol-1. Trimers in solution are also remarkably stable, exhibiting minimal distortions in their spatial and orientational distributions. The resulting geometrical docking shows the incipient characteristics of flexible nanotubes in solution, with eventual water interchange between the central channel and the bulk at the junctions between monomers. Structural and dynamical properties of the trapped water filling the nanotube are dictated to a large exent by the charge density at the rims.Fil: Rodriguez, Javier. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; ArgentinaFil: Semino, Rocio. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; ArgentinaFil: Laria, Daniel Hector. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentin

Equilibrium and dynamical characteristics of Imidazole Langmuir monolayers on graphite sheets

Using molecular dynamics techniques, we examine structural and dynamical characteristics of liquid-like, Imidazole(Im) monolayers physisorbed onto a planar graphite sheet, at T = 384K. Our  simulations reveal that molecular orientations in the saturated monolayer exhibit a bistable distribution, characterized by an inner parallel arrangement of the molecules in close contact with the substrate and a slanted alignment, in those lying in adjacent, outer locations. Compared to the results  found in three dimensional, bulk phases, the analysis of the spatial correlations between sites participating in hydrogen bonding shows a clear enhancement of the intermolecular interactions, which also leads to stronger dipolar correlations. As a result, the gross structural features of the   monolayer can be cast in terms of mesoscopic domains, comprising units articulated via winding hydrogen bonds, that persist along  typical time intervals of a few tens of picoseconds.  On the dynamical side, a similar comparison of the characteristic decorrelation time for orientational motions shows a 4-fold increment.  Contrasting, the reduction of the system dimensionality leads to a larger diffusion constant. Possible substrate-induced anisotropies in the diffusive motions are also investigated.Fil: Rodriguez, Javier. Comisión Nacional de Energía Atómica; Argentina. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Núcleo de Investigación en Educacion Ciencia y Tecnologia; ArgentinaFil: Elola, Maria Dolores. Comisión Nacional de Energía Atómica; ArgentinaFil: Laria, Daniel Hector. Comisión Nacional de Energía Atómica; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentin

Coaxial cross-diffusion through carbon nantoubes

We present results from nonequilibrium molecular dynamics experiments describing the relaxation of local concentrations at two reservoirs, initially filled with water (W) and acetonitrile (ACN), as they become connected through a membrane composed of (16,16) carbon nanotubes. Within the hydrophobic nanotube cavities, the equilibrium concentrations contrast sharply to those observed at the reservoirs, with a clear enhancement of ACN, in detriment of W. From the dynamical side, the relaxation involves three well-differentiated stages; the first one corresponds to the equilibration of individual concentrations within the nanotubes. An intermediate interval with Fickian characteristics follows, during which the overall transport can be cast in terms of coaxial opposite fluxes, with a central water domain segregated from an external ACN shell, in close contact with the tube walls. We also found evidence of a third, much slower, mechanism to reach equilibration, which involves structural modifications of tightly bound solvation shells, in close contact with the nanotube rims.Fil: Rodriguez, Javier. Comisión Nacional de Energía Atómica. Gerencia del Área de Investigaciones y Aplicaciones no Nucleares. Gerencia de Física (Centro Atómico Constituyentes); Argentina. Universidad Nacional de San Martín; ArgentinaFil: Elola, Maria Dolores. Comisión Nacional de Energía Atómica. Gerencia del Área de Investigaciones y Aplicaciones no Nucleares. Gerencia de Física (Centro Atómico Constituyentes); Argentina. Universidad Nacional de San Martín; ArgentinaFil: Laria, Daniel Hector. Comisión Nacional de Energía Atómica. Gerencia del Área de Investigaciones y Aplicaciones no Nucleares. Gerencia de Física (Centro Atómico Constituyentes); Argentin

Hydrogen bond dynamics at water/Pt interfaces

We present results from computer simulations that shed light on structural and dynamic characteristics of hydrogen bonding of aqueous phases at ambient conditions, at the close vicinity of electrified metal interfaces. Our simulation strategy relied on the consideration of a Hamiltonian that explicitly incorporates effects from polarization fluctuations at the metal surface, induced by the instantaneous local electric field promoted by the partial charges at the solvent molecules. Compared to bulk environments, our results reveal important modifications in the hydrogen bond architectures that critically depend on the atomic arrangements of the interfaces exposed to the liquid phases and the net charges allocated at the metal plates. These modifications have equally important consequences on the characteristic time scales describing the ruptures of hydrogen bonds which are operated by mechanisms which are absent in descriptions that omit atomic detail and polarization fluctuations at the metal plates. We also analyze how the latter modifications are translated into spectral shifts in the stretching bands of infrared spectra of water adlayers.Fil: Videla, Pablo Ernesto. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Química Inorgánica, Analítica y Química Física; ArgentinaFil: Ansourian, Lisandro. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Química Inorgánica, Analítica y Química Física; ArgentinaFil: Laria, Daniel Hector. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Química Inorgánica, Analítica y Química Física; Argentina. Comisión Nacional de Energía Atómica; Argentin

Polar mixtures under nanoconfinement

We present results from molecular dynamics simulations describing structural and dynamical characteristics of equimolar mixtures of water and acetonitrile, confined between two silica walls separated at interplate distances of d = 0.6, 1, and 1.5 nm. Two different environments were investigated: a first one where wall-solvent dispersion forces prevail (hydrophobic confinement) and a second one in which the terminal O atoms at the silica surface are transformed into silanol groups (hydrophilic confinement). For the former case, we found that, at the shortest interplate distance examined, the confined region is devoid of water molecules. At an interplate distance of the order of 1 nm, water moves into the confined region, although, in all cases, there is a clear enhancement of the local concentration of acetonitrile in detriment of that of water. Within hydrophilic environments, we found clear distinctions between a layer of bound water lying in close contact with the silica substrates and a minority of confined water that occupies the inner liquid slab. The bound aqueous layer is fully coordinated to the silanol groups and exhibits minimal hydrogen bonding with the second solvation layer, which exclusively includes acetonitrile molecules. Dynamical characteristics of the solvent mixture are analyzed in terms of diffusive and rotational motions in both environments. Compared to bulk mixtures, we found significant retardations in all dynamical modes, with those ascribed to water molecules bound to the hydrophilic plates being the most dramatic.Fil: Rodriguez, Javier. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Comisión Nacional de Energía Atómica; Argentina. Universidad Nacional de San Martín. Escuela de Ciencia y Tecnología; ArgentinaFil: Elola, Maria Dolores. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Comisión Nacional de Energía Atómica; ArgentinaFil: Laria, Daniel Hector. Comisión Nacional de Energía Atómica; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Química Inorgánica, Analítica y Química Física; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentin

Ionic liquid aqueous solutions under nanoconfinement

We extend our previous molecular dynamics analysis of confined aqueous electrolytes within cylindrical hydrophobic pores of nanometric dimensions [Videla et al. J. Chem. Phys.2011, 135, 104503] to the case of room temperature ionic liquid (RTIL) solutions, with concentrations close to c ∼ 1 M. Equilibrium and dynamical characteristics of two imidazolium-based RTILs, differing in the hydrophobicity of the corresponding anionic species, were considered. The solutions within the pore were modeled in contact with "bulk-like" reservoirs, which served as reference systems to gauge the magnitude of the modifications observed in the global densities and in the transport coefficients. The density fields associated to the ionic species present a marked enhancement near the pore walls; this leads to increments of the global RTIL concentration within the pores, which are intermediate between 2 and 3 times the ones observed in the bulk reservoirs. These modifications are more marked in solutions containing more hydrophobic anionic species. In both cases, selective adsorption of imidazolium groups at the pore walls prevails; these wall-solvation states are characterized by a parallel orientation of the imidazolium ring, with respect to the pore surface. Mass and charge transport were also investigated. The segregation of the ionic species towards the pore wall promotes a sharp drop in the individual ionic diffusion coefficients. Nonuniform trends in the modifications of the ionic conductivity were found. Our results show that charge transport is the result of a complex interplay between competing effects involving modifications in the local concentrations, retardations in the ionic mobility, and dynamical cross-correlations, as well. A physical interpretation of the latter effects is provided in terms of the differences in the spatial correlations of the ionic species within the interior of the pore. © 2012 American Chemical Society.Fil: Rodriguez, Javier. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Comisión Nacional de Energía Atómica; Argentina. Universidad Nacional de San Martín. Escuela de Ciencia y Tecnología; ArgentinaFil: Elola, Maria Dolores. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Comisión Nacional de Energía Atómica; ArgentinaFil: Laria, Daniel Hector. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina. Comisión Nacional de Energía Atómica; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Química Inorgánica, Analítica y Química Física; Argentin

Confinement of polar solvents within β-cyclodextrins

Using molecular dynamics techniques, we examined equilibrium and dynamical characteristics pertaining to the solvation of a single β-cyclodextrin (CD) in water and in dimethylsulfoxide (DMSO). Compared to its global minimum structure, the overall shape of the solute in solution is reasonably well preserved. While in aqueous solutions, the average number of solvent molecules retained within the central cavity of the oligosaccharide is close to 5, for DMSO, that number reduces to ∼1. No evidence of significant orientational correlations of the trapped molecules were found in either solvent. The main contributions to the hydrogen-bond (HB) connectivity between the solute and the bulk phases are due to the more distal HO6-O6 hydroxyl groups, acting as HB donors and acceptors. The average residence time for retained DMSO was found to be in the nanosecond range, and it is, at least, 1 order of magnitude longer that the one observed for water. We also analyzed the characteristics of the solvation of the β-CD in an equimolar water-DMSO mixture. In this environment, we found a preferential localization of a single DMSO molecule in the interior of the CD and a very minor retention of water. In the mixture, the characteristic time of residence of the trapped DMSO molecule increases by a factor of ∼2. The observed difference was rationalized in terms of the fluctuations of the local concentrations of the two species in the vicinity of the CD top and bottom rims.Fil: Rodriguez, Javier. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; ArgentinaFil: Rico, Daniel Hernán. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; ArgentinaFil: Domenianni, Luis Ignacio. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; ArgentinaFil: Laria, Daniel Hector. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química, Física de los Materiales, Medioambiente y Energía. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química, Física de los Materiales, Medioambiente y Energía; Argentina. Comisión Nacional de Energía Atómica; Argentin