1
|
Lazaridis T, Sepehri A. Amino acid deprotonation rates from classical force fields. J Chem Phys 2022; 157:085101. [DOI: 10.1063/5.0101960] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Acid ionization constants (pKa's) of titratable amino acid side chains have received a large amount of experimental and theoretical attention. In many situations, however, the rates of protonation and deprotonation, kon & koff, may also be important, for example in understanding the mechanism of action of proton channels or membrane proteins that couple proton transport to other processes. Protonation and deprotonation involve making and breaking of covalent bonds, which cannot be studied by classical force fields. However, environment effects on the rates should be captured by such methods. Here we present an approach for estimating deprotonation rates based on Warshel's extension of Marcus's theory of electron transfer, with input from molecular simulations. The missing bond dissociation energy is represented by a constant term determined by fitting the pKa value in solution. The statistics of the energy gap between protonated and deprotonated states is used to compute free energy curves of the two states and thus free energy barriers, from which the rate can be deduced. The method is applied to Glu, Asp, and His in bulk solution and select membrane proteins: the M2 proton channel, bacteriorhodopsin, and cytochrome c oxidase.
Collapse
Affiliation(s)
- Themis Lazaridis
- Department of Chemistry, City College of New York, United States of America
| | | |
Collapse
|
2
|
Joutsuka T, Ando K. Efficient Free-energy Calculation of Proton Transfer by Constrained Density Functional Theory and Geometrically Restrained Molecular Dynamics Simulation. CHEM LETT 2021. [DOI: 10.1246/cl.210132] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Tatsuya Joutsuka
- Institute of Quantum Beam Science, Graduate School of Science and Engineering, Ibaraki University, Hitachi, Ibaraki 316-8511, Japan
| | - Koji Ando
- Department of Information and Sciences, Tokyo Woman’s Christian University, 2-6-1 Zenpukuji, Suginami-ku, Tokyo 167-8585, Japan
| |
Collapse
|
3
|
Joutsuka T, Ando K. Constrained Density Functional Theory Molecular Dynamics Simulation of Deprotonation in Aqueous Silicic Acid. J Phys Chem B 2020; 124:8323-8330. [PMID: 32856917 DOI: 10.1021/acs.jpcb.0c05096] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The solubility of silicic acid in water is so low that the molecular mechanism behind the physical properties such as pKa remains poorly understood, despite the importance in fields such as chemistry and geology. Theoretical calculations provide the molecular-level information on such a rare chemical species yet face difficulties in the selection of reaction coordinates and the rare-event sampling of proton transfers (PTs) in condensed phases. We propose a straightforward calculation scheme of pKa with ab initio molecular dynamics (MD) simulation and the constrained density functional theory (CDFT), which provides structural and dynamical properties such as radial distribution functions, vibrational spectra, and reaction paths. The calculated reaction free energies of deprotonations agreed with experiments within a few kcal/mol. Analysis of the solvation structure shows that, after deprotonation, the hydronium ion D3O+ repels away from the deprotonated silicic acid SiO(OD)3- without forming a stable contact-ion pair. The calculated vibrational spectra are consistent with the spectroscopic measurements, and the dynamical analysis of the reaction path quantifies the couplings of the OD stretch vibrations of silicic acid and water to the PT reaction in terms of the vertical energy-gap coordinate that includes both the solute and the solvent components.
Collapse
Affiliation(s)
- Tatsuya Joutsuka
- Institute of Quantum Beam Science, Graduate School of Science and Engineering, Ibaraki University, Hitachi 316-8511 Ibaraki, Japan
| | - Koji Ando
- Department of Information and Sciences, Tokyo Womans Christian University, 2-6-1 Zenpukuji, Suginami-ku, Tokyo 167-8585, Japan
| |
Collapse
|
4
|
Prampolini G, Ingrosso F, Segalina A, Caramori S, Foggi P, Pastore M. Dynamical and Environmental Effects on the Optical Properties of an Heteroleptic Ru(II)–Polypyridine Complex: A Multilevel Approach Combining Accurate Ground and Excited State QM-Derived Force Fields, MD and TD-DFT. J Chem Theory Comput 2018; 15:529-545. [DOI: 10.1021/acs.jctc.8b01031] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Francesca Ingrosso
- Université de Lorraine, CNRS, Laboratoire de Physique et Chimie Théoriques, F-54000 Nancy, France
| | - Alekos Segalina
- Université de Lorraine, CNRS, Laboratoire de Physique et Chimie Théoriques, F-54000 Nancy, France
| | - Stefano Caramori
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Ferrara, Via Luigi Borsari 46, I-44100, Ferrara, Italy
| | - Paolo Foggi
- European Laboratory for Non Linear Spectroscopy (LENS), Università di Firenze, Via Nello Carrara 1, I-50019 Sesto Fiorentino Florence, Italy
- INO−CNR, Istituto Nazionale di Ottica, Consiglio Nazionale delle Ricerche, Largo Fermi 6, I-50125 Florence, Italy
- Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, Via Elce di Sotto 8, I-06123 Perugia, Italy
| | - Mariachiara Pastore
- Université de Lorraine, CNRS, Laboratoire de Physique et Chimie Théoriques, F-54000 Nancy, France
| |
Collapse
|
5
|
Mushrif SH, Varghese JJ, Krishnamurthy CB. Solvation dynamics and energetics of intramolecular hydride transfer reactions in biomass conversion. Phys Chem Chem Phys 2015; 17:4961-9. [DOI: 10.1039/c4cp05063k] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Solvent dynamics and non-equilibrium solvation alter the energetics of the hydride transfer step and thus, can significantly affect reaction kinetics.
Collapse
Affiliation(s)
- Samir H. Mushrif
- School of Chemical and Biomedical Engineering
- Nanyang Technological University
- Singapore 637459
| | - Jithin J. Varghese
- School of Chemical and Biomedical Engineering
- Nanyang Technological University
- Singapore 637459
| | | |
Collapse
|
6
|
Nibbering ET, Chudoba C, Elsaesser T. Hydrogen-Bond Dynamics and Solvation of Electronically Excited States as Determined by Femtosecond Vibrational Spectroscopy. Isr J Chem 2013. [DOI: 10.1002/ijch.199900040] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
7
|
Iida K, Sato H. Theoretical study on ionization process in aqueous solution. J Chem Phys 2012; 136:144510. [DOI: 10.1063/1.3700225] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|
8
|
YOSHIMORI AKIRA. TIME DEPENDENT DENSITY FUNCTIONAL METHODS AND THEIR APPLICATION TO CHEMICAL PHYSICS. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2011. [DOI: 10.1142/s0219633604000878] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
This article reviews microscopic development of time dependent functional method and its application to chemical physics. It begins with the formulation of density functional theory. The time dependent extension is discussed after the equilibrium formulation. Its application is explained by solvation dynamics. In addition, it reviews studies of nonlinear effects on polar liquids and simple mixtures.
Collapse
Affiliation(s)
- AKIRA YOSHIMORI
- Department of Physics, Kyushu University, Fukuoka, 812-8581, Japan
| |
Collapse
|
9
|
Joutsuka T, Ando K. Dynamics of Proton Transfer and Vibrational Relaxation in Dilute Hydrofluoric Acid. J Phys Chem A 2011; 115:678-84. [DOI: 10.1021/jp108413p] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tatsuya Joutsuka
- Department of Chemistry, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan
| | - Koji Ando
- Department of Chemistry, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
10
|
Nishihara Y, Kato S, Hayashi S. Protein collective motions coupled to ligand migration in myoglobin. Biophys J 2010; 98:1649-57. [PMID: 20409486 DOI: 10.1016/j.bpj.2009.12.4318] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2009] [Revised: 12/06/2009] [Accepted: 12/18/2009] [Indexed: 11/19/2022] Open
Abstract
Ligand migration processes inside myoglobin and protein dynamics coupled to the migration were theoretically investigated with molecular dynamics simulations. Based on a linear response theory, we identified protein motions coupled to the transient migration of ligand, carbon monoxide (CO), through channels. The result indicates that the coupled protein motions involve collective motions extended over the entire protein correlated with local gating motions at the channels. Protein motions, coupled to opening of a channel from the distal pocket to a neighboring xenon site, were found to share the collective motion with experimentally observed protein motions coupled to a doming motion of the heme Fe atom upon photodissociation of the ligand. Analysis based on generalized Langevin dynamics elucidated slow and diffusive features of the protein response motions. Remarkably small transmission coefficients for rates of the CO migrations through myoglobin were found, suggesting that the CO migration dynamics are characterized as motions governed by the protein dynamics involving the collective motions, rather than as thermally activated transitions across energy barriers of well-structured channels.
Collapse
Affiliation(s)
- Yasutaka Nishihara
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto Japan
| | | | | |
Collapse
|
11
|
Nagata Y, Lennartz C. Atomistic simulation on charge mobility of amorphous tris(8-hydroxyquinoline) aluminum (Alq3): Origin of Poole–Frenkel–type behavior. J Chem Phys 2008; 129:034709. [DOI: 10.1063/1.2949506] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
|
12
|
Suzuki Y, Tanimura Y. Exploring a free energy landscape by means of multidimensional infrared and terahertz spectroscopies. J Chem Phys 2008; 128:164501. [PMID: 18447453 DOI: 10.1063/1.2897982] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
A model for the dipolar crystal system is employed to explore a role of free energy landscape (FEL), in which dipolar molecules are posted on two-dimensional lattice sites with two-state libratinal dynamics. All dipole-dipole interactions are included to have frustrated interactions among the dipoles. For the regular and distorted lattice cases, the FEL is calculated from the interaction energies and the total polarizations for all possible dipolar states at various temperatures. At high temperatures, the shape of the calculated FEL is smooth and parabolic, while it becomes bumpy at low temperatures exhibiting multiple local minima. To study dynamical aspects of the system, the single flip dynamics and the single-double mixed flips dynamics of dipoles are examined from a master equation approach. As the observables of linear absorption and two-dimensional (2D) infrared, far infrared, and terahertz spectroscopies, the first- and third-order response functions of polarization are calculated for different physical conditions characterized by the FEL. While the linear absorption signals decay in time in a similar manner regardless of the FEL profiles, the 2D signals exhibit prominent differences for those profiles. This indicates that we may differentiate the FEL profiles by changing two-time valuables in 2D spectroscopy. As illustrated in the single-double flips case, the FEL study by means of 2D spectroscopy, however, relies on the dynamics which is set independently from the FEL. The Smoluchowski equation is applied to examine the description of the collective dynamics on the microscopically calculated FEL. We found that the one-dimensional and 2D signals calculated from the Smoluchowski equation agree with those from master equation only at temperatures where the FEL becomes parabolic shape.
Collapse
Affiliation(s)
- Yohichi Suzuki
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa, Sakyoku, Kyoto 606-8502, Japan.
| | | |
Collapse
|
13
|
Ando K, Hynes JT. Acid-Base Proton Transfer and Ion Pair Formation in Solution. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470141694.ch6] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
14
|
Suzuki Y, Tanimura Y. Free energy landscapes of electron transfer system in dipolar environment below and above the rotational freezing temperature. J Chem Phys 2007; 126:054504. [PMID: 17302482 DOI: 10.1063/1.2431172] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Electron transfer reaction in a polar solvent is modeled by a solute dipole surrounded by dipolar molecules with simple rotational dynamics posted on the three-dimensional distorted lattice sites. The interaction energy between the solute and solvent dipoles as a reaction coordinate is adopted and free energy landscapes are calculated by generating all possible states for a 26 dipolar system and by employing Wang-Landau sampling algorithm for a 92 dipolar system. For temperatures higher than the energy scale of dipole-dipole interactions, the free energy landscapes for the small reaction coordinate region have quadratic shape as predicted by Marcus [Rev. Mod. Phys. 65, 599 (1993)] whereas for the large reaction coordinate region, the landscapes exhibit a nonquadratic shape. When the temperature drops, small notched structures appear on the free energy profiles because of the frustrated interactions among dipoles. The formation of notched structure is analyzed with statistical approach and it is shown that the amplitude of notched structure depend upon the segment size of the reaction coordinate and is characterized by the interaction energy among the dipoles. Using simulated free energy landscapes, the authors calculate the reaction rates as a function of the energy gap for various temperatures. At high temperature, the reactions rates follow a bell shaped (inverted parabolic) energy gap law in the small energy gap regions, while it becomes steeper than the parabolic shape in a large energy gap regions due to the nonquadratic shape of the free energy landscape. The peak position of parabola also changes as the function of temperature. At low temperature, the profile of the reaction rates is no longer smooth because of the many local minima of the free energy landscape.
Collapse
Affiliation(s)
- Yohichi Suzuki
- Department of Chemistry, Kyoto University, Oiwakecho, Kitashirakawa, Sakyoku, Kyoto 606-8502, Japan
| | | |
Collapse
|
15
|
Fernandez-Ramos A, Miller JA, Klippenstein SJ, Truhlar DG. Modeling the kinetics of bimolecular reactions. Chem Rev 2007; 106:4518-84. [PMID: 17091928 DOI: 10.1021/cr050205w] [Citation(s) in RCA: 393] [Impact Index Per Article: 23.1] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Antonio Fernandez-Ramos
- Departamento de Quimica Fisica, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | | | | | | |
Collapse
|
16
|
Suzuki Y, Tanimura Y. Free energy landscape analysis of two-dimensional dipolar solvent model at temperatures below and above the rotational freezing point. J Chem Phys 2006; 124:124508. [PMID: 16599698 DOI: 10.1063/1.2178785] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Ionic solvation in a polar solvent is modeled by a central charge surrounded by dipolar molecules posted on two-dimensional distorted lattice sites with simple rotational dynamics. Density of states is calculated by applying the Wang-Landau algorithm to both the energy and polarization states. The free energy landscapes of solvent molecules as a function of polarization are depicted to explore the competition between the thermal fluctuation and solvation energy. Without a central charge, for temperatures higher than the energy scale of the dipole-dipole interactions, the energy landscape for the small polarization region exhibits a parabolic shape as predicted by Marcus [Rev. Mod. Phys. 65, 599 (1993)] for electron transfer reaction, while there is an additional quartic contribution to the landscape for the large polarization region. When the temperature drops, the simulated free energy landscapes are no longer smooth due to the presence of multiple local minima arising from the frustrated interaction among the dipoles. The parabolic contribution becomes negligible and the energy landscape becomes quartic in shape. For a strong central charge, the energy landscape exhibits an asymmetric profile due to the contributions of linear and cubic terms that arise from the charge-dipole interactions.
Collapse
Affiliation(s)
- Yohichi Suzuki
- Department of Chemistry, Kyoto University, Oiwakecho, Kitashirakawa, Sakyoku, Kyoto 606-8502, Japan.
| | | |
Collapse
|
17
|
Quantum versus classical electron transfer energy as reaction coordinate for the aqueous Ru2+/Ru3+ redox reaction. Theor Chem Acc 2005. [DOI: 10.1007/s00214-005-0058-0] [Citation(s) in RCA: 71] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
18
|
Higashi M, Kato S. Theoretical Study on Electronic and Spin Structures of [Fe2S2]2+,+ Cluster: Reference Interaction Site Model Self-Consistent Field (RISM-SCF) and Multireference Second-Order Møller−Plesset Perturbation Theory (MRMP) Approach. J Phys Chem A 2005; 109:9867-74. [PMID: 16833302 DOI: 10.1021/jp0581429] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Electronic structures of [Fe(2)S(2)(SCH(3))(4)](2-,3-) in DMSO solution are calculated using reference interaction site model complete active space self-consistent field (RISM-CASSCF)/multireference second-order Møller-Plesset perturbation theory (MRMP) method. For the reduced state, we obtain both the low-spin Fe(3+)Fe(2+) localized and high-spin Fe(2.5+)Fe(2.5+) delocalized forms, which are very close in energy. The spin interaction constants obtained from the energies of states with various spin multiplicities are in good agreement with the available experimental estimates both for the oxidized and for the reduced states. The dynamic electron correlation effect is found to be important in estimating the spin interaction between the Fe ions. The redox potentials are calculated to be 2.87 and 2.78 eV for the localized and delocalized reduced states, respectively, which are close to the experimental values. We devise a simple model for calculating the free energy curves of the reduction process based on the RISM-SCF theory. The activation barrier height is calculated to be 7.4 kcal/mol at the equilibrium geometry of oxidized state, indicating that the reduction reaction will occur efficiently in DMSO solvent. The effect of solvent fluctuation on the free energy profiles is discussed on the basis of the present calculations.
Collapse
Affiliation(s)
- Masahiro Higashi
- Department of Chemistry, Graduate School of Science, Kyoto Unicersity, Kitashirakawa, Sakyo-ku, Kyoto 606-8502, Japan
| | | |
Collapse
|
19
|
Iuchi S, Morita A, Kato S. Electronic relaxation dynamics of Ni2+-ion aqueous solution: Molecular-dynamics simulation. J Chem Phys 2005; 123:24505. [PMID: 16050757 DOI: 10.1063/1.1949212] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Electronic relaxation dynamics of Ni2+-ion aqueous solution is investigated using molecular-dynamics (MD) simulations with the model-effective Hamiltonian developed previously. The nonadiabatic transition rates from the first three excited states to the ground state are evaluated by the golden rule formula with the adiabatic MD simulations. The MD simulations with the fewest-switch surface-hopping method are also carried out to obtain a more detailed description of the electronic relaxation dynamics among the excited states. We found out that the transitions among the three excited states are very fast, in the order of 10 fs, while the transition between the excited and ground states is slow, about 800 ps. These findings are consistent with the time scales of energy dissipation detected by the transient lens experiment. In both simulations, we explore the effects of the quantum decoherence, where the decoherence functions are derived by the energy-gap dynamics with the displaced harmonic-oscillator model.
Collapse
Affiliation(s)
- Satoru Iuchi
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa, Sakyo-ku, Kyoto 606-8502, Japan
| | | | | |
Collapse
|
20
|
Bernardi E, Marques Martins M, Stassen H. The breakdown of linear response theory in non-polar solvation dynamics. Chem Phys Lett 2005. [DOI: 10.1016/j.cplett.2005.03.067] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
21
|
Li S, Thompson WH. Proton Transfer in Nanoconfined Polar Solvents. 1. Free Energies and Solute Position. J Phys Chem B 2005; 109:4941-6. [PMID: 16863151 DOI: 10.1021/jp045036i] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The reaction free energy curves for a model phenol-amine proton-transfer system in a confined CH3Cl solvent have been calculated by Monte Carlo simulations. The free energy curves, as a function of a collective solvent coordinate, have been obtained for several fixed reaction complex radial positions (based on the center-of-mass). A smooth, hydrophobic spherical cavity was used to confine the solvent, and radii of 10 and 15 A have been considered. Quantum effects associated with the transferring proton have been included by adding the proton zero-point energy to the classical free energy. The results indicate the reaction complex position can be an important component of the reaction coordinate for proton-transfer reactions in nanoconfined solvents.
Collapse
Affiliation(s)
- Shenmin Li
- Department of Chemistry and Chemical Engineering, Dalian University, Dalian 116622, PR China
| | | |
Collapse
|
22
|
Faeder J, Ladanyi BM. Solvation Dynamics in Reverse Micelles: The Role of Headgroup−Solute Interactions. J Phys Chem B 2005; 109:6732-40. [PMID: 16851757 DOI: 10.1021/jp045202m] [Citation(s) in RCA: 107] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
We present molecular dynamics simulation results for solvation dynamics in the water pool of anionic-surfactant reverse micelles (RMs) of varying water content, w(0). The model RMs are designed to represent water/aerosol-OT/oil systems, where aerosol-OT is the common name for sodium bis(2-ethylhexyl)sulfosuccinate. To determine the effects of chromophore-headgroup interactions on solvation dynamics, we compare the results for charge localization in model ionic diatomic chromophores that differ only in charge sign. Electronic excitation in both cases is modeled as charge localization on one of the solute sites. We find dramatic differences in the solvation responses for anionic and cationic chromophores. Solvation dynamics for the cationic chromophore are considerably slower and more strongly w(0)-dependent than those for the anionic chromophore. Further analysis indicates that the difference in the responses can be ascribed in part to the different initial locations of the two chromophores relative to the surfactant interface. In addition, slow motion of the cationic chromophore relative to the interface is the main contributor to the longer-time decay of the solvation response to charge localization in this case.
Collapse
Affiliation(s)
- James Faeder
- Theoretical Biology and Biophysics Group, MS K710, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | |
Collapse
|
23
|
Ingrosso F, Ladanyi BM, Mennucci B, Elola MD, Tomasi J. Solvation Dynamics in Acetonitrile: A Study Incorporating Solute Electronic Response and Nuclear Relaxation. J Phys Chem B 2005; 109:3553-64. [PMID: 16851393 DOI: 10.1021/jp0456032] [Citation(s) in RCA: 56] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.
Collapse
Affiliation(s)
- Francesca Ingrosso
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa, Italy
| | | | | | | | | |
Collapse
|
24
|
Blumberger J, Sprik M. Ab Initio Molecular Dynamics Simulation of the Aqueous Ru2+/Ru3+ Redox Reaction: The Marcus Perspective. J Phys Chem B 2005; 109:6793-804. [PMID: 16851765 DOI: 10.1021/jp0455879] [Citation(s) in RCA: 71] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A well-behaved (low spin) transition metal aqua ion, Ru(aq)(2+), is used as a model system in an ab initio molecular dynamics study of a redox half reaction to which the Marcus theory of electron transfer is assumed to apply. Using constraint methods, we show that aqueous Ru(2+) can be reversibly transformed to Ru(3+) under the control of the classical solvent electrostatic potential as order parameter. The mean force is found to vary linearly with the order parameter in accordance with the Marcus theory. As can be expected for a half reaction, the slope in the oxidized and reduced states are asymmetric differing by approximately a factor of two. As a further test, we verify that the corresponding quadratic potential of mean force is in excellent agreement with the free energy profile obtained from the Gaussian distribution of potential fluctuations sampled from free (unconstrained) runs of the reduced and oxidized system. Similar to experimental electrochemical methods, our simulation scheme enables us to manipulate the effective thermodynamic driving force and align the free energy minima of product and reactant state. The activation energy and reaction entropy computed under these conditions are discussed and analyzed from the Marcus perspective.
Collapse
Affiliation(s)
- Jochen Blumberger
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | | |
Collapse
|
25
|
Martins LR, Skaf MS, Ladanyi BM. Solvation Dynamics at the Water/Zirconia Interface: Molecular Dynamics Simulations. J Phys Chem B 2004. [DOI: 10.1021/jp0470896] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lucimara R. Martins
- Instituto de Química, Universidade Estadual de Campinas, Cx. P. 6154, Campinas-SP 13083-970, Brazil
| | - Munir S. Skaf
- Instituto de Química, Universidade Estadual de Campinas, Cx. P. 6154, Campinas-SP 13083-970, Brazil
| | - Branka M. Ladanyi
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523
| |
Collapse
|
26
|
Grabowski ZR, Rotkiewicz K, Rettig W. Structural Changes Accompanying Intramolecular Electron Transfer: Focus on Twisted Intramolecular Charge-Transfer States and Structures. Chem Rev 2003; 103:3899-4032. [PMID: 14531716 DOI: 10.1021/cr940745l] [Citation(s) in RCA: 2251] [Impact Index Per Article: 107.2] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
27
|
Ohmiya K, Kato S. Solution reaction path Hamiltonian based on reference interaction site model self-consistent field method: Application to Menshutkin-type reactions. J Chem Phys 2003. [DOI: 10.1063/1.1580803] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Vieceli J, Benjamin I. Solvation Dynamics at the Interface between Water and Self-assembled Monolayers. J Phys Chem B 2003. [DOI: 10.1021/jp0222043] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- John Vieceli
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, California 95064
| | - Ilan Benjamin
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, California 95064
| |
Collapse
|
29
|
|
30
|
Ladanyi BM, Perng BC. Solvation Dynamics in Dipolar−Quadrupolar Mixtures: A Computer Simulation Study of Dipole Creation in Mixtures of Acetonitrile and Benzene. J Phys Chem A 2002. [DOI: 10.1021/jp013580a] [Citation(s) in RCA: 54] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Branka M. Ladanyi
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523
| | - Baw-Ching Perng
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523
| |
Collapse
|
31
|
Ohmiya K, Kato S. RISM-SCF study for the rate constant of SN2 reaction CH3Cl+Cl− in aqueous solution. Chem Phys Lett 2001. [DOI: 10.1016/s0009-2614(01)01091-0] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
32
|
Faeder J, Ladanyi BM. Solvation Dynamics in Aqueous Reverse Micelles: A Computer Simulation Study. J Phys Chem B 2001. [DOI: 10.1021/jp010632n] [Citation(s) in RCA: 150] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- James Faeder
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523
| | - Branka M. Ladanyi
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523
| |
Collapse
|
33
|
Ando K. A stable fluctuating-charge polarizable model for molecular dynamics simulations: Application to aqueous electron transfers. J Chem Phys 2001. [DOI: 10.1063/1.1394923] [Citation(s) in RCA: 40] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
Ando K. Solvent nuclear quantum effects in electron transfer reactions. III. Metal ions in water. Solute size and ligand effects. J Chem Phys 2001. [DOI: 10.1063/1.1369128] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Ando K. Solvent nuclear quantum effects in electron transfer reactions. II. Molecular dynamics study on methanol solution. J Chem Phys 2001. [DOI: 10.1063/1.1367384] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
Kobayashi C, Saito S, Ohmine I. Mechanism of fast proton transfer in ice: Potential energy surface and reaction coordinate analyses. J Chem Phys 2000. [DOI: 10.1063/1.1319636] [Citation(s) in RCA: 59] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Changenet-Barret P, Choma CT, Gooding EF, DeGrado WF, Hochstrasser RM. Ultrafast Dielectric Response of Proteins from Dynamics Stokes Shifting of Coumarin in Calmodulin. J Phys Chem B 2000. [DOI: 10.1021/jp001634v] [Citation(s) in RCA: 105] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
38
|
Aherne D, Tran V, Schwartz BJ. Nonlinear, Nonpolar Solvation Dynamics in Water: The Roles of Electrostriction and Solvent Translation in the Breakdown of Linear Response. J Phys Chem B 2000. [DOI: 10.1021/jp000326u] [Citation(s) in RCA: 83] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
39
|
Hilczer M, Tachiya M. Computer simulation of electron transfer reactions between two molecular species in solvents of various polarity. J Mol Liq 2000. [DOI: 10.1016/s0167-7322(99)00130-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
40
|
Rettig W, Zietz B. Do twisting and pyramidalization contribute to the reaction coordinate of charge-transfer formation in DMABN and derivatives? Chem Phys Lett 2000. [DOI: 10.1016/s0009-2614(99)01323-8] [Citation(s) in RCA: 72] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
41
|
Ando K, Hynes JT. Molecular Mechanism of HF Acid Ionization in Water: An Electronic Structure−Monte Carlo Study. J Phys Chem A 1999. [DOI: 10.1021/jp992481i] [Citation(s) in RCA: 99] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Koji Ando
- Institute of Materials Science, University of Tsukuba, Tsukuba, Ibaraki 305-8573, Japan
| | - James T. Hynes
- Department of Chemistry and Biochemistry, University of Colorado, Boulder, Colorado 80309-0215
| |
Collapse
|
42
|
Kovalenko A, Hirata F. Potential of Mean Force between Two Molecular Ions in a Polar Molecular Solvent: A Study by the Three-Dimensional Reference Interaction Site Model. J Phys Chem B 1999. [DOI: 10.1021/jp991300+] [Citation(s) in RCA: 162] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
43
|
Mercer IP, Gould IR, Klug DR. A Quantum Mechanical/Molecular Mechanical Approach to Relaxation Dynamics: Calculation of the Optical Properties of Solvated Bacteriochlorophyll-a. J Phys Chem B 1999. [DOI: 10.1021/jp990284d] [Citation(s) in RCA: 48] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ian P. Mercer
- Department of Chemistry, Imperial College, South Kensington, London, SW7 2AY, U.K
| | - Ian R. Gould
- Department of Chemistry, Imperial College, South Kensington, London, SW7 2AY, U.K
| | - David R. Klug
- Department of Chemistry, Imperial College, South Kensington, London, SW7 2AY, U.K
| |
Collapse
|
44
|
Ishida T, Hirata F, Kato S. Solvation dynamics of benzonitrile excited state in polar solvents: A time-dependent reference interaction site model self-consistent field approach. J Chem Phys 1999. [DOI: 10.1063/1.479083] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Free energy profiles of electron transfer at water–electrode interface studied by the reference interaction site model theory. Chem Phys Lett 1999. [DOI: 10.1016/s0009-2614(99)00372-3] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
46
|
Kiefer AM, Kast SM, Wasielewski MR, Laukenmann K, Kothe G. Exploring the Structure of a Photosynthetic Model by Quantum-Chemical Calculations and Time-Resolved Q-Band Electron Paramagnetic Resonance. J Am Chem Soc 1999. [DOI: 10.1021/ja981930+] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
47
|
Hilczer M, Tachiya M. Effect of solvent polarity on the potential of mean force between two molecular ions: MD simulation. Chem Phys Lett 1998. [DOI: 10.1016/s0009-2614(98)00970-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
48
|
Ladanyi BM, Maroncelli M. Mechanisms of solvation dynamics of polyatomic solutes in polar and nondipolar solvents: A simulation study. J Chem Phys 1998. [DOI: 10.1063/1.476911] [Citation(s) in RCA: 121] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
49
|
Georgievskii Y, Hsu CP, Marcus RA. Dynamic Stokes shift in solution: Effect of finite pump pulse duration. J Chem Phys 1998. [DOI: 10.1063/1.476155] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
|