Molecular Simulations of Carbon Dioxide and Water: Cation Solvation Louise J

Total Page:16

File Type:pdf, Size:1020Kb

Molecular Simulations of Carbon Dioxide and Water: Cation Solvation Louise J Article pubs.acs.org/est Molecular Simulations of Carbon Dioxide and Water: Cation Solvation Louise J. Criscenti* and Randall T. Cygan Geochemistry Department, P.O. Box 5800 MS 0754, Sandia National Laboratories, Albuquerque, New Mexico 87185-0754, United States ABSTRACT: Proposed carbon dioxide sequestration scenarios in sedimentary reservoirs require investigation into the interactions between supercritical carbon dioxide, brines, and the mineral phases found in the basin and overlying caprock. Molecular simulations can help to understand the partitioning of metal cations between aqueous solutions and supercritical carbon dioxide where limited experimental data exist. In this effort, we used classical molecular dynamics simulations to compare the solvation of alkali and alkaline-earth metal cations in fl water and liquid CO2 at 300 K by combining a exible simple point charge fl fi model for water and an accurate exible force eld for CO2. Solvation energies for these cations are larger in water than in carbon dioxide, suggesting that they will partition preferentially into water. In both aqueous and CO2 solutions, the solvation energies decrease with cation size and increase with cation charge. However, changes in solvation energy with ionic radii are smaller in CO2 than in water suggesting that the partitioning of cations into CO2 will increase with ion size. Simulations of the interface between + + aqueous solution and supercritical CO2 support this suggestion in that some large cations (e.g., Cs and K ) partition into the CO2 phase, often with a partial solvation sphere of water molecules. ■ INTRODUCTION Because of limited experimental data on metal partitioning To develop accurate regional scale models to evaluate the fate between supercritical CO2 and brine solutions, it is valuable to of supercritical CO in various CO sequestration scenarios, it is consider the use of computational chemistry methods to 2 2 explore this subject. Metal partitioning into the supercritical important to assess the exchange of solutes across the 1−7 CO phase may influence chemical reactions such as the supercritical CO −water interface. Experimental studies of 2 2 dissolution of primary phases or precipitation of secondary the solubility of CO in water and simple aqueous solutions 2 minerals, and physical properties such as both the density and have been used to determine the range of CO concentrations 2 viscosity of the CO phase, and the contact angles formed trapped in idealized reservoir fluids in regional scale multiphase 2 between supercritical CO , water, and mineral phases. Super- models. However, the interaction of supercritical CO with the 2 2 critical CO has been considered as a potential fluid for complex brine compositions found in aquifers have typically 2 extracting metal ions from liquid phases. The ions cannot been ignored in such hydrological simulations, most likely readily be extracted in supercritical CO , but the use of metal− because of the lack of data (solubility, density, partition 2 organic complexes greatly increases metal partitioning into the coefficients, etc.) for these interactions at the temperatures and 15 supercritical phase. This suggests that metal speciation in the pressures expected in geological reservoirs. Recent experimental 8−12 brine may also impact species partitioning between the two and simulation studies of the geochemistry of reservoir fluids. fluids and mineral interactions have provided new insights on fi 13 The research presented here provides a rst look at the supercritical CO2 behavior in geological media. Rempel et al. solvation of alkali and alkaline-earth metal cations in CO fluids present experimental data for the fractionation of Na, Fe, Cu, 2 − and the likelihood that these cations will partition between and Zn between brine and CO2 for pressures of 6.6 16 MPa at brines and supercritical CO2 in sedimentary basin environ- 340 K. The data exhibit an increase in Na concentrations with ments. An overall goal of our research is to use molecular increasing CO2 density and suggest that Fe, Cu, and Zn could simulations to model the structure and behavior of injected be transferred to adjacent aquifer systems. However, this supercritical CO in subsurface reservoirs and how it eventually experimental approach is recently developed and difficult to 2 implement. Therefore, geochemical modeling of reactions Special Issue: Carbon Sequestration among reservoir rocks, brine, and supercritical CO2 focus on investigating the impact of high PCO2 on reactions between Received: April 23, 2012 saline brines and minerals under reservoir conditions, but, as a Revised: July 6, 2012 first approximation, neglect the solubility of water, salts, or Accepted: July 10, 2012 metals in the supercritical phase.14 Published: July 10, 2012 © 2012 American Chemical Society 87 dx.doi.org/10.1021/es301608c | Environ. Sci. Technol. 2013, 47, 87−94 Environmental Science & Technology Article interacts with resident aqueous fluids and minerals over time. In Molecular Dynamics Simulations. Classical molecular the present study we use molecular dynamics to investigate the dynamics (MD) was used to evaluate the structure and − structure of the CO2 brine interface and how electrolytes, thermodynamics of the solvation of alkali and alkaline earth fi − speci cally alkali and alkaline earth cations, are partitioned cations by CO2, and to examine the structure of the water 22 between aqueous and CO2 phases. We report the application of CO2 interfaces. The Forcite software was used to evaluate the fl fi an accurate and exible CO2 force eld, compare cation electrostatic and short-range interactions for each atomic fi solvation in water and CO2 systems, and examine the solubility con guration and time step. MD simulations were performed of CO2 in water. using a canonical NVT thermodynamic ensemble, maintaining a fixed cell volume for a fixed number of atoms. Temperature was ■ COMPUTATIONAL METHODS controlled using the Nose−Hoover23 method. One million Carbon Dioxide Force Field. For the molecular time steps of 1 fs were used to obtain 1 ns of simulation time; fl fi simulations of this study, we use the fully exible CO2 force atomic con gurations were saved every 1000 time steps (1 ps) field developed by Cygan et al.16 The three-point force field for efficient data storage and trajectory analysis. MD 17 − − uses the refined Van der Waals parameters of Zhu et al. that simulations of systems involving water CO2 and brine CO2 were derived from the original model of Harris and Young.18 interfaces were completed for 0.5 ns of simulation time using The force field also allows for full intramolecular bond stretch the same frequency of molecular configuration storage. and angle bend and thus correctly predicts the vibrational Cation Solvation in Liquid Carbon Dioxide. MD spectra of CO2. The CO2 potential is compatible with most simulations of alkali and alkaline earth cation solvation in 19 point-charge force fields including Clayff and other related liquid CO2 were performed using the NVT thermodynamic force fields for the simulation of geochemical systems. The ensemble with 512 CO2 molecules in a simulation box (V = 3 3 interatomic potentials should be appropriate for simulating 48.23 nm ; ρ = 0.77 g/cm at 8 MPa) at 300 K for comparison complex CO2 systems involving multiple phases, interfaces, and with cation solvation in liquid water which is typically the addition of dissolved cations in both water and CO phases. determined at ambient temperature and pressure (300 K, 0.1 2 24,25 In the Cygan et al.16 force field, the total potential energy of MPa). A simulation was performed for 1 ns to equilibrate the molecular system is described by Coulombic and van der the liquid CO2 in a cubic periodic cell with a cell length of 36.4 Waals contributions representing nonbonded energies and Å. Subsequently, each cation was placed in the equilibrated bond-stretch and angle-bend terms representing the intra- model CO2 liquid. A similar set of simulations was completed molecular energies. The long-range Coulombic or electrostatic for the same number of CO2 molecules at supercritical energy is given by eq 1, where q and q are the partial charges of conditions (350 K, 20 MPa) using a cubic periodic cell with i j 3 ρ 3 the atoms and rij is the distance between the atoms. a cell length of 39.2 Å (V = 60.24 nm ; = 0.62 g/cm ). For qq both conditions, each solvation model was initially geometry- ij optimized for 500 iterations prior to 1 ns of MD simulation. ECoul = rij (1) Charge balance was established by using a charge-neutralizing background.20,26 The solvation energy for a particular cation The short-range van der Waals energy is given by was estimated by subtracting the potential energy of the ⎡ 12 6⎤ equilibrated liquid CO2 simulation cell from the potential ⎛ σσ⎞ ⎛ ⎞ ⎢⎜ ij ⎟ ⎜ ij ⎟ ⎥ energy of the equilibrated simulation cell containing liquid CO2 EVDW=−4ε ij⎢⎜ ⎟ ⎜ ⎟ ⎥ − (or supercritical CO2) and the cation [i.e., PE(CO2 + cation) ⎣⎝ rrij ⎠ ⎝ ij ⎠ ⎦ (2) PE(CO2 cell)]. The van der Waals parameters used for the + + + + 2+ ε σ alkali (Li ,Na,K, and Cs ) and alkaline earth cations (Mg , where ij and ij are optimized for intermolecular interactions 2+ 2+ 2+ using values determined by Zhu et al.17 The interaction Ca ,Sr ,Ba ) are provided in Table 1. Four cations of each parameters between unlike atoms are calculated according to σ Table 1. Force Field Parameters for Cations the arithmetic mean rule for the distance parameter, ij, and the 20,21 ε a ε σ geometric mean rule for the energy parameter, ij. cation i (kJ/mol) i (Å) Harmonic potentials are used for the bond stretch and angle Li+ 40 0.0766 1.3723 bend terms. Na+ 25 0.5447 2.3500 + 41 1 2 K 0.4187 3.3340 Ekrr=−() + 42 Stretch2 S ij o (3) Cs 0.4187 3.8310 Mg2+ 40 3.6634 1.6444 1 2+ 41 Ek=−()θθ2 Ca 0.4187 2.8720 Bend2 B ijk o (4) Sr2+ 40 0.4187 3.4620 Ba2+ 40 0.1968 3.8166 Both expressions describe the increase in the potential energy aSmith and Dang,25 Koneshan et al.,41 Smith and Dang,42 Åqvist40.
Recommended publications
  • Microscopic Dynamics of Charge Separation at the Aqueous Electrochemical Interface
    Microscopic dynamics of charge separation at the aqueous electrochemical interface John A. Kattirtzia,b, David T. Limmerc,d,e,1, and Adam P. Willarda,1 aDepartment of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02138; bCollege of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China; cDepartment of Chemistry, University of California, Berkeley, CA 94609; dKavili Energy NanoScience Institute, University of California, Berkeley, CA 94609; and eMaterial Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94609 Edited by Michael L. Klein, Temple University, Philadelphia, PA, and approved June 6, 2017 (received for review March 7, 2017) We have used molecular simulation and methods of importance the process of aqueous ion association (forward reaction above) sampling to study the thermodynamics and kinetics of ionic or dissociation (reverse reaction above) is deceptively simple. charge separation at a liquid water–metal interface. We have con- The expression omits the cooperative role of solvent, which sidered this process using canonical examples of two different must restructure to accommodate transitions between the asso- classes of ions: a simple alkali–halide pair, Na+I−, or classical ions, ciated and dissociated states (2, 11). This restructuring, which + − and the products of water autoionization, H3O OH , or water is driven by thermal fluctuations, is both collective (extending ions. We find that for both ion classes, the microscopic mecha- beyond the first solvation shell) and fleeting (12–14). These nism of charge separation, including water’s collective role in the microscopic processes are difficult to probe experimentally (15), process, is conserved between the bulk liquid and the electrode so atomistic simulation has played an important role in reveal- interface.
    [Show full text]
  • Solvent Effects on the Thermodynamic Functions of Dissociation of Anilines and Phenols
    University of Wollongong Research Online University of Wollongong Thesis Collection 1954-2016 University of Wollongong Thesis Collections 1982 Solvent effects on the thermodynamic functions of dissociation of anilines and phenols Barkat A. Khawaja University of Wollongong Follow this and additional works at: https://ro.uow.edu.au/theses University of Wollongong Copyright Warning You may print or download ONE copy of this document for the purpose of your own research or study. The University does not authorise you to copy, communicate or otherwise make available electronically to any other person any copyright material contained on this site. You are reminded of the following: This work is copyright. Apart from any use permitted under the Copyright Act 1968, no part of this work may be reproduced by any process, nor may any other exclusive right be exercised, without the permission of the author. Copyright owners are entitled to take legal action against persons who infringe their copyright. A reproduction of material that is protected by copyright may be a copyright infringement. A court may impose penalties and award damages in relation to offences and infringements relating to copyright material. Higher penalties may apply, and higher damages may be awarded, for offences and infringements involving the conversion of material into digital or electronic form. Unless otherwise indicated, the views expressed in this thesis are those of the author and do not necessarily represent the views of the University of Wollongong. Recommended Citation Khawaja, Barkat A., Solvent effects on the thermodynamic functions of dissociation of anilines and phenols, Master of Science thesis, Department of Chemistry, University of Wollongong, 1982.
    [Show full text]
  • Hydrated Sulfate Clusters SO4^2
    Article Cite This: J. Phys. Chem. B 2019, 123, 4065−4069 pubs.acs.org/JPCB 2− n − Hydrated Sulfate Clusters SO4 (H2O)n ( =1 40): Charge Distribution Through Solvation Shells and Stabilization Maksim Kulichenko,† Nikita Fedik,† Konstantin V. Bozhenko,‡,§ and Alexander I. Boldyrev*,† † Department of Chemistry and Biochemistry, Utah State University, 0300 Old Main Hill, Logan, Utah 84322-0300, United States ‡ Department of Physical and Colloid Chemistry, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St, Moscow 117198, Russian Federation § Institute of Problems of Chemical Physics, Russian Academy of Sciences, Chernogolovka 142432, Moscow Region, Russian Federation *S Supporting Information 2− ABSTRACT: Investigations of inorganic anion SO4 interactions with water are crucial 2− for understanding the chemistry of its aqueous solutions. It is known that the isolated SO4 dianion is unstable, and three H2O molecules are required for its stabilization. In the 2− current work, we report our computational study of hydrated sulfate clusters SO4 (H2O)n (n =1−40) in order to understand the nature of stabilization of this important anion by fi 2− water molecules. We showed that the most signi cant charge transfer from dianion SO4 ≤ 2− to H2O takes place at a number of H2O molecules n 7. The SO4 directly donates its charge only to the first solvation shell and surprisingly, a small amount of electron density of 0.15|e| is enough to be transferred in order to stabilize the dianion. Upon further addition of ff ≤ H2O molecules, we found that the cage e ect played an essential role at n 12, where the fi 2− | | rst solvation shell closes.
    [Show full text]
  • Solvation Structure of Ions in Water
    International Journal of Thermophysics, Vol. 28, No. 2, April 2007 (© 2007) DOI: 10.1007/s10765-007-0154-6 Solvation Structure of Ions in Water Raymond D. Mountain1 Molecular dynamics simulations of ions in water are reported for solutions of varying solute concentration at ambient conditions for six cations and four anions in 10 solutes. The solutes were selected to show trends in properties as the size and charge density of the ions change. The emphasis is on how the structure of water is modified by the presence of the ions and how many water molecules are present in the first solvation shell of the ions. KEY WORDS: anion; cation; molecular dynamics; potential functions; solva- tion; water. 1. INTRODUCTION The behavior of aqueous solutions of salts is a topic of continuing interest [1–3]. In this note we examine how the structure of water,as reflected in the oxy- gen–oxygen pair function, is modified as the concentration of ions increases. The method of generating the pair functions is molecular dynamics using the SPC/E model for water [4] and various models for ion–water interactions taken from the literature. The solutes are LiCl, NaCl, KCl, RbCl, NaF, CaCl2, CaSO4,Na2SO4, NaNO3, and guanidinium chloride (GdmCl) [C(NH2)3Cl]. This work is an extension of earlier work [5] to a larger set of ions. The water molecules and ions interact through site–site pair interac- tions consisting of Lennard–Jones potentials and Coulomb interactions so that the interaction between a pair of sites labeled i, j and separated by an interatomic distance r is 12 6 φij (r) = 4ij (σij /r) − (σij /r) + qiqj /r (1) where qi is the charge on site i.
    [Show full text]
  • Correlation Between Solubility and the Number of Water Molecules in the First Solvation Shell
    Pure &App/. Chem., Vol. 70, No. 10, pp. 1895-1904, 1998. Printed in Great Britain. 0 1998 IUPAC Solubility of gases in water: Correlation between solubility and the number of water molecules in the first solvation shell Pirketta Scharlin,a Rubin Battino,b Estanislao Silla,c Iiiaki Tuii6nc and Juan Luis Pascual-Ahuirc a Department of Chemistry, University of Turku, FIN-20014 Turku, Finland Department of Chemistry, Wright State University, Dayton, Ohio 45435, USA Departamento de Quimica Fisica, Universidad de Valbncia, 46100 Burjassot, Valbncia, Spain Abstract: Using a new version of a program called GEPOL, a consistent set of values for areas of three different kinds of surfaces for 53 gaseous solutes was computed by Silla et al. These surface areas, together with the volumes of space enclosed by the surfaces, are reported in the present paper. The three surfaces are the van der Waals Surface (WS), the Solvent Accessible Surface (SAS) and the Solvent-Excluding Surface (SES). Values for the number of water molecules (N) in the first solvation shell are estimated by a simple surface area approach from the SAS data. Values of N, as well as literature data on solubilities of gases in water, are used to study various semi-empirical correlations between thermodynamic changes on solution and the number of water molecules in the first solvation shell. Dilute aqueous solutions of gases have been of continuous experimental and theoretical interest. Reliable solubility data exist for a large variety of gases in water. '-18 Experimental data have been important for validating the results of theoretical calculations, computer simulations, and various practical applications.
    [Show full text]
  • The Ionization Constant of Water Over Wide Ranges of Temperature and Density
    The Ionization Constant of Water over Wide Ranges of Temperature and Density Andrei V. Bandura Department of Quantum Chemistry, St. Petersburg State University, 26 University Prospect, Petrodvoretz, St. Petersburg 198504, Russia Serguei N. Lvova… The Energy Institute and Department of Energy & Geo-Environmental Engineering, The Pennsylvania State University, 207 Hosler Building, University Park, Pennsylvania 16802 ͑Received 26 April 2004; revised manuscript received 1 March 2005; accepted 5 April 2005; published online 8 December 2005͒ A semitheoretical approach for the ionization constant of water, KW , is used to fit the available experimental data over wide ranges of density and temperature. Statistical ther- modynamics is employed to formulate a number of contributions to the standard state chemical potential of the ionic hydration process. A sorption model is developed for calculating the inner-shell term, which accounts for the ion–water interactions in the immediate ion vicinity. A new analytical expression is derived using the Bragg–Williams approximation that reproduces the dependence of a mean ion solvation number on the solvent chemical potential. The proposed model was found to be correct at the zero- density limit. The final formulation has a simple analytical form, includes seven adjust- able parameters, and provides good fitting of the collected KW data, within experimental uncertainties, for a temperature range of 0 – 800 °C and densities of 0 – 1.2 g cmϪ3.©2006 American Institute of Physics. ͓DOI: 10.1063/1.1928231͔ Key words: Bragg–Williams approximation; high temperatures and low densities; ionization constant of water; statistical thermodynamics solvation model. Contents 3. References on the experimentally obtained KW data used for estimating the empirical parameters Nomenclature.............................
    [Show full text]
  • Structure of Solvated Metal Ions
    Structure of solvated metal ions Solution and crystal structure of Ga3+, In3+, Sc3+, Y3+, La3+ and Ca2+ ions with water and non-aqueous oxygen donor solvents Patric Lindqvist-Reis Doctoral Thesis Department of Chemistry Stockholm 2000 Patric Lindqvist-Reis Department of Chemistry Inorganic Chemistry Royal Institute of Technology S-100 44 Stockholm Sweden Printed by Kista Snabbtryck AB ON THE COVER: 3+ Sc(H2O)8 ions in the crystal structure of [Sc(H2O)8](CF3SO3)3 ROYAL ISBN 91-7170-569-4 INSTITUTE OF ISSN 0348-825X TECHNOLOGY TRITA-OOK-1060 STRUCTURE OF SOLVATED METAL IONS Solution and crystal structure of Ga3+, In3+, Sc3+, Y3+, La3+ and Ca2+ ions with water and non-aqueous oxygen donor solvents Patric Lindqvist-Reis Akademisk avhandling som med tillstånd av Kungliga Tekniska Högskolan i Stockholm framlägges till offentlig granskning för avläggande av filosofie doktorsexamen i oorganisk kemi, onsdagen den 14 juni, kl. 10.00 i kollegiesalen, administrationsbyggnaden, KTH, Vallhallavägen 79, Stockholm Avhandlingen försvaras på engelska. Abstract The structure of the solvated group 3 ions scandium(III), yttrium(III) and lanthanum(III), has been determined in aqueous solution and in some oxygen donor solvents, and compared with the hydrated group 13 ions gallium(III) and indium(III). A combination of X-ray absorption fine structure (XAFS), large angle X-ray scattering (LAXS), crystallography and vibration spectroscopy has been used for the structure studies of the hydrated ions and their dimethylsulfoxide and N,N´-dimethylpropylene urea solvates in solution and in the solid state. For the hexahydrated gallium(III) and indium(III) ions in solution, the metal-oxygen distances were found to be 1.959(6) Å and 2.131(7) Å to the first hydration shell, and 4.05(1) and 4.13(1) Å, respectively, to fairly well-defined second hydration shells containing about twelve water molecules.
    [Show full text]
  • A Computational Study of Mg2+ Dehydration in Aqueous Solution in the Presence of HS and Other Monovalent Anions – Insights to Dolomite Formation
    Available online at www.sciencedirect.com Geochimica et Cosmochimica Acta 88 (2012) 77–87 www.elsevier.com/locate/gca A computational study of Mg2+ dehydration in aqueous solution in the presence of HSÀ and other monovalent anions – Insights to dolomite formation Yang Yang a,b,c, Nita Sahai a,b,d,e,⇑, Christopher S. Romanek f,g, Suvankar Chakraborty f,g a Department of Geoscience, University of Wisconsin, Madison, 1215 West Dayton Street, Madison, WI 53706, United States b NASA Astrobiology Institute, University of Wisconsin, Madison, 1215 West Dayton Street, Madison, WI 53706, United States c Department of Chemistry and Biochemistry, Rowan University, 201 Mullica Hill Rd, Glassboro, NJ 08028-1701, United States d Department of Polymer Science, 170 University Avenue, University of Akron, Akron, OH 44325-3909, United States e NASA Astrobiology Institute, 170 University Avenue, University of Akron, Akron, OH 44325-3909, United States f Department of Earth and Environmental Sciences, University of Kentucky, Lexington, KY 40506, United States g NASA Astrobiology Institute, University of Kentucky, Lexington, KY 40506, United States Received 4 November 2011; accepted in revised form 13 March 2012; available online 29 March 2012 Abstract Massive sedimentary dolomite formed at near-Earth’s surface temperatures is abundant in the ancient geological rock record compared to modern deposition. Extensive experimental work to synthesize dolomite at low temperature and to reveal the for- mation mechanism has been attempted previously. Sulfide, the product of bacterial sulfate reduction, has been proposed in the literature to play an active role in promoting dolomite formation by facilitating desolvation of Mg2+ in the bulk solution and, thus, incorporation into the dolomite crystal structure.
    [Show full text]
  • The Effects of Dissolved Halide Anions on Hydrogen Bonding in Liquid Water Jared D
    Published on Web 10/25/2007 The Effects of Dissolved Halide Anions on Hydrogen Bonding in Liquid Water Jared D. Smith,†,‡ Richard J. Saykally,†,‡ and Phillip L. Geissler*,†,‡,§ Contribution from the Department of Chemistry, UniVersity of California, Berkeley, California 94720, and Chemical and Material Sciences DiVisions, Lawrence Berkeley National Laboratory, Berkeley, California 94720 Received March 19, 2007; E-mail: [email protected] Abstract: It is widely believed that the addition of salts to water engenders structural changes in the hydrogen-bond network well beyond the adjacent shell of solvating molecules. Classification of many ions as “structure makers” and “structure breakers” has been based in part on corresponding changes in the vibrational spectra (Raman and IR). Here we show that changes in O-H vibrational spectra induced by the alkali halides in liquid water result instead from the actions of ions’ electric fields on adjacent water molecules. Computer simulations that accurately reproduce our experimental measurements suggest that the statistics of hydrogen-bond strengths are only weakly modified beyond this first solvation shell. Introduction liquid after the mass and volume associated with the ions have been removed from the solution.4) In contrast, recent ultrafast The physical and chemical properties of ions in aqueous infrared pump-probe measurements appear to indicate that solution are centrally important for many biological, geological, dissolved ions only significantly affect vibrational and orien- environmental, and industrial processes.1 Yet despite enormous tational relaxation times of water molecules in the first hydration efforts, several important issues remain unresolved regarding shell,5,6 implying that structural rearrangements are strongly the microscopic structure of the surrounding solvent and localized.
    [Show full text]
  • Step-Wise Hydration of Magnesium by Four Water Molecules Precedes Phosphate Release in a Myosin Motor
    bioRxiv preprint doi: https://doi.org/10.1101/817254; this version posted October 24, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Step-wise Hydration of Magnesium by Four Water Molecules Precedes Phosphate Release in a Myosin Motor M.L. Mugnai1 and D. Thirumalai1 1Department of Chemistry, The University of Texas at Austin, Austin, TX 78712 October 23, 2019 1 bioRxiv preprint doi: https://doi.org/10.1101/817254; this version posted October 24, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Abstract Molecular motors, such as myosin, kinesin, and dynein, convert the energy re- leased by the hydrolysis of ATP into mechanical work, which allows them to undergo directional motion on cytoskeletal tracks. This process is achieved through synchro- nization between the catalytic activity of the motor and the associated changes in its conformation. A pivotal step in the chemomechanical transduction in myosin motors occurs after they bind to the actin filament, which triggers the release of phosphate (Pi, product of ATP hydrolysis) and the rotation of the lever arm. Here, we investigate the mechanism of phosphate release in myosin VI, which has been debated for over two decades, using extensive molecular dynamics simulations in- volving multiple trajectories each several µs long.
    [Show full text]
  • Enhancing Electrochemical Intermediate Solvation Through Electrolyte Anion Selection to Increase Nonaqueous Li‐O2 Battery Capacity
    Enhancing electrochemical intermediate solvation through electrolyte anion selection to increase nonaqueous Li‐O2 battery capacity 1,2 3 4 Colin M. Burke3 , Vikram Pande , Abhishek Khetan1,2 , Venkatasubramanian Viswanathan1Department of Chemical* and Bryan D. McCloskey and Biomolecular Engineering,* University of California, Berkeley, CA, 94720 2Environmental Energy Technologies Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720 3Department of Mechanical Engineering, Carnegie Mellon University, Pittsburgh, PA, 15213 4Institute for Combustion Technology, RWTH, Aachen 52056, Germany Keywords:*[email protected], [email protected] Significance Statement donor number, solubility, lithium nitrate, NMR, Li‐air battery The Li‐air battery has attracted significant interest as a potential high‐energy alternative to Li‐ion batteries. However, the battery discharge product, lithium peroxide, is both an electronic insulator and insoluble in nonaqueous electrolytes. It therefore passivates the battery cathode as it is uniformly deposited and disallows the battery to achieve even a modest fraction of its potential electrochemical capacity. Our research objective is to circumvent this challenge by enhancing the solubility of electrochemically formed intermediate species. We present a rational basis for electrolyte (i.e. solvent and salt) selection for nonaquoeus Li‐air batteries and + demonstrate a selection criteria for an electrolyte salt that increases the stability of Li in solution, thereby triggering a solution I. Abstractbased process that allows significantly improved battery capacities. 2 Among the ‘beyond Li‐ion’ battery chemistries, nonaqueous Li‐O batteries have the highest theoretical specific energy and as a result have attracted significant research attention over 2 the past decade. A critical scientific challenge facing nonaqueous Li‐O batteries is the electronically insulating nature of the primary discharge product, lithium peroxide, which passivates the battery cathode as it is formed, leading to low ultimate cell capacities.
    [Show full text]
  • Observation of Water Dangling OH Bonds Around Dissolved Nonpolar Groups
    Observation of water dangling OH bonds around dissolved nonpolar groups P. N. Pereraa, K. R. Fegaa, C. Lawrenceb, E. J. Sundstroma, J. Tomlinson-Phillipsa, and Dor Ben-Amotza,1 aDepartment of Chemistry, Purdue University, West Lafayette, IN 47907; and bDepartment of Chemistry, Grand Valley State University, Allendale, MI 49401 Edited by Michael D. Fayer, Stanford University, Stanford, CA, and approved June 12, 2009 (received for review April 3, 2009) We report the experimental observation of water dangling OH spectra (20), and here, we employ it to uncover dangling OH bonds in the hydration shells around dissolved nonpolar (hydro- vibrational bands arising from the hydration shells around dissolved carbon) groups. The results are obtained by combining vibrational nonpolar groups. (Raman) spectroscopy and multivariate curve resolution (MCR), to The remainder of this article begins by describing the results reveal a high-frequency OH stretch peak arising from the hydration of the experiments we performed to confirm our assignment of shell around nonpolar (hydrocarbon) solute groups. The frequency the observed peak to water dangling OH bonds. Changes in the and width of the observed peak is similar to that of dangling OH intensity of the dangling OH peak are then used to estimate the bonds previously detected at macroscopic air–water and oil–water number of water dangling OH bonds around hydrocarbon chains interfaces. The area of the observed peak is used to quantify the of different size. The width of the observed dangling OH peak number of water dangling bonds around hydrocarbon chains of is used to establish a lower bound on the lifetime of the different length.
    [Show full text]