1
|
Curk T, Yuan J, Luijten E. Accelerated simulation method for charge regulation effects. J Chem Phys 2022; 156:044122. [PMID: 35105090 DOI: 10.1063/5.0066432] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The net charge of solvated entities, ranging from polyelectrolytes and biomolecules to charged nanoparticles and membranes, depends on the local dissociation equilibrium of individual ionizable groups. Incorporation of this phenomenon, charge regulation (CR), in theoretical and computational models requires dynamic, configuration-dependent recalculation of surface charges and is therefore typically approximated by assuming constant net charge on particles. Various computational methods exist that address this. We present an alternative, particularly efficient CR Monte Carlo method (CR-MC), which explicitly models the redistribution of individual charges and accurately samples the correct grand-canonical charge distribution. In addition, we provide an open-source implementation in the large-scale Atomic/Molecular Massively Parallel Simulator molecular dynamics (MD) simulation package, resulting in a hybrid MD/CR-MC simulation method. This implementation is designed to handle a wide range of implicit-solvent systems that model discreet ionizable groups or surface sites. The computational cost of the method scales linearly with the number of ionizable groups, thereby allowing accurate simulations of systems containing thousands of individual ionizable sites. By matter of illustration, we use the CR-MC method to quantify the effects of CR on the nature of the polyelectrolyte coil-globule transition and on the effective interaction between oppositely charged nanoparticles.
Collapse
Affiliation(s)
- Tine Curk
- Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, USA
| | - Jiaxing Yuan
- School of Physics and Astronomy and Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Erik Luijten
- Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|
2
|
Rivard U, Thomas V, Bruhacs A, Siwick B, Iftimie R. Donor-Bridge-Acceptor Proton Transfer in Aqueous Solution. J Phys Chem Lett 2014; 5:3200-3205. [PMID: 26276332 DOI: 10.1021/jz501378d] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We use ab initio molecular dynamics to study proton transfer in a donor-bridge-acceptor system in which the bridge is a single water molecule and the entire system is embedded in aqueous solution. The results, based on a large number of proton transfer trajectories, demonstrate that the dominant charge-transfer pathway is a subpicosecond "through bridge" event in which the bridge adopts an Eigen-like (hydronium) structure. We also identify another state in which the bridge forms a Zundel-like configuration with the acceptor that appears to be a dead end for the charge transfer. The reaction coordinate is inherently multidimensional and, as we demonstrate, cannot be given in terms of either local structural parameters of the donor-bridge-acceptor system or local solvent coordination numbers.
Collapse
Affiliation(s)
- Ugo Rivard
- †Département de Chimie, Université de Montréal, CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| | - Vibin Thomas
- †Département de Chimie, Université de Montréal, CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| | - Andrew Bruhacs
- ‡Departments of Chemistry and Physics, Center for the Physics of Materials, McGill University, 801 Sherbrooke Street West, Montréal H3A2K6, Canada
| | - Bradley Siwick
- ‡Departments of Chemistry and Physics, Center for the Physics of Materials, McGill University, 801 Sherbrooke Street West, Montréal H3A2K6, Canada
| | - Radu Iftimie
- †Département de Chimie, Université de Montréal, CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| |
Collapse
|
3
|
Schmidt TC, Paasche A, Grebner C, Ansorg K, Becker J, Lee W, Engels B. QM/MM investigations of organic chemistry oriented questions. Top Curr Chem (Cham) 2014; 351:25-101. [PMID: 22392477 DOI: 10.1007/128_2011_309] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.
Collapse
Affiliation(s)
- Thomas C Schmidt
- Institut für Phys. und Theor. Chemie, Emil-Fischer-Strasse 42, Campus Hubland Nord, 97074, Würzburg, Germany
| | | | | | | | | | | | | |
Collapse
|
4
|
Iftimie R, Tremblay MH, Thomas V, Hétu S, de Lasalle F, Rivard U. Moderately strong phenols dissociate by forming an ion-pair kinetic intermediate. J Phys Chem A 2013; 117:13976-87. [PMID: 24299203 DOI: 10.1021/jp410858d] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We show computational evidence that ground-state moderately strong hydroxyarenes (Ar-OH, pKa ∼ 0) dissociate by forming an ion-pair intermediate that lives for 3-5 ps. The concentration of this intermediate is approximately 2 times smaller than that of the un-ionized acid at pH ∼ 0.6 and is characterized by average C-O bond lengths (1.30 Å) that are intermediate between those of un-ionized (1.29 Å) and fully dissociated (1.34 Å) species. During the lifetime of the ion-pair intermediate the excess proton fluctuates between the oxygen atom of the phenolic moiety and those of water molecules in the first and second solvation shells on a subpicosecond time scale (∼100-300 fs).
Collapse
Affiliation(s)
- Radu Iftimie
- Département de Chimie, Université de Montréal , CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| | | | | | | | | | | |
Collapse
|
5
|
Thomas V, Rivard U, Maurer P, Bruhács A, Siwick BJ, Iftimie R. Concerted and Sequential Proton Transfer Mechanisms in Water-Separated Acid-Base Encounter Pairs. J Phys Chem Lett 2012; 3:2633-2637. [PMID: 26295883 DOI: 10.1021/jz3012639] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The proton transfer mechanisms involved inside aqueous, solvent-separated encounter complexes between phenol and carboxyl moieties are studied using ab initio molecular dynamics and computational time-resolved vibrational spectroscopy. This model framework can be viewed as a ground-state analog of the excited-state proton transfer reactions that have been actively investigated using ultrafast spectroscopy. Three qualitatively distinct proton transfer pathways are observed in the simulations. These can be described as direct concerted, direct sequential, and through bulk transfers. The primary difference between the sequential and concerted mechanism is the involvement of a reaction intermediate in which the proton fluctuates for several picoseconds through the hydrogen bonds connecting donor and acceptor but resides primarily on an intervening water molecule in the encounter complex. These results contribute to our molecular level understanding of the diverse processes involved in proton transfer within water-separated encounter complexes.
Collapse
Affiliation(s)
- Vibin Thomas
- †Département de Chimie, Université de Montréal,CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| | - Ugo Rivard
- †Département de Chimie, Université de Montréal,CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| | - Patrick Maurer
- †Département de Chimie, Université de Montréal,CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| | - Andrew Bruhács
- ‡Departments of Chemistry and Physics, Center for the Physics of Materials, McGill University, 801 Sherbrooke Street West, Montréal, Canada
| | - Bradley J Siwick
- ‡Departments of Chemistry and Physics, Center for the Physics of Materials, McGill University, 801 Sherbrooke Street West, Montréal, Canada
| | - Radu Iftimie
- †Département de Chimie, Université de Montréal,CP 6128, succursale Centre-Ville, Montréal H3C3J7, Canada
| |
Collapse
|
6
|
Glaser R, Jost M. Disproportionation of bromous acid HOBrO by direct O-transfer and via anhydrides O(BrO)2 and BrO-BrO2. An ab initio study of the mechanism of a key step of the Belousov-Zhabotinsky oscillating reaction. J Phys Chem A 2012; 116:8352-65. [PMID: 22871057 DOI: 10.1021/jp301329g] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The results are reported of an ab initio study of the thermochemistry and of the kinetics of the HOBrO disproportionation reaction 2HOBrO (2) ⇄ HOBr (1) + HBrO(3) (3), reaction ( R4' ), in gas phase (MP2(full)/6-311G*) and aqueous solution (SMD(MP2(full)/6-311G*)). The reaction energy of bromous acid disproportionation is discussed in the context of the coupled reaction system R2-R4 of the FKN mechanism of the Belousov-Zhabotinsky reaction and considering the acidities of HBr and HOBrO(2). The structures were determined of ten dimeric aggregates 4 of bromous acid, (HOBrO)(2), of eight mixed aggregates 5 formed between the products of disproportionation, (HOBr)(HOBrO(2)), and of four transition states structures 6 for disproportionation by direct O-transfer. It was found that the condensation of two HOBrO molecules provides facile access to bromous acid anhydride 7, O(BrO)(2). A discussion of the potential energy surface of Br(2)O(3) shows that O(BrO)(2) is prone to isomerization to the mixed anhydride 8, BrO-BrO(2), and to dissociation to 9, BrO, and 10, BrO(2), and their radical pair 11. Hence, three possible paths from O(BrO)(2) to the products of disproportionation, HOBr and HOBrO(2), are discussed: (1) hydrolysis of O(BrO)(2) along a path that differs from its formation, (2) isomerization of O(BrO)(2) to BrO-BrO(2) followed by hydrolysis, and (3) O(BrO)(2) dissociation to BrO and BrO(2) and their reactions with water. The results of the potential energy surface analysis show that the rate-limiting step in the disproportionation of HOBrO consists of the formation of the hydrate 12a of bromous acid anhydride 7 via transition state structure 14a. The computed activation free enthalpy ΔG(act)(SMD) = 13.6 kcal/mol for the process 2·2a → [14a](‡) → 12a corresponds to the reaction rate constant k(4) = 667.5 M(-1) s(-1) and is in very good agreement with experimental measurements. The potential energy surface analysis further shows that anhydride 7 is kinetically and thermodynamically unstable with regard to hydrolysis to HOBr and HOBrO(2) via transition state structure 14b. The transition state structure 14b is much more stable than 14a, and, hence, the formation of the "symmetrical anhydride" from bromous acid becomes an irreversible reaction for all practical purposes because 7 will instead be hydrolyzed as a "mixed anhydride" to afford HOBr and HOBrO(2). The mixed anhydride 8, BrO-BrO(2), does not play a significant role in bromous acid disproportionation.
Collapse
Affiliation(s)
- Rainer Glaser
- Department of Chemistry, University of Missouri, Columbia, Missouri 65211, United States.
| | | |
Collapse
|
7
|
Mangold M, Rolland L, Costanzo F, Sprik M, Sulpizi M, Blumberger J. Absolute pKa Values and Solvation Structure of Amino Acids from Density Functional Based Molecular Dynamics Simulation. J Chem Theory Comput 2011; 7:1951-61. [PMID: 26596456 DOI: 10.1021/ct100715x] [Citation(s) in RCA: 72] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Absolute pKa values of the amino acid side chains of arginine, aspartate, cysteine, histidine, and tyrosine; the C- and N-terminal group of tyrosine; and the tryptophan radical cation are calculated using a revised density functional based molecular dynamics simulation technique introduced previously [ Cheng , J. ; Sulpizi , M. ; Sprik , M. J. Chem. Phys. 2009 , 131 , 154504 ]. In the revised scheme, acid deprotonation is considered as a dissociation rather than a proton transfer reaction, and a correction term for treating the proton as a hydronium ion is suggested. The acidity constants of the amino acids are obtained from the vertical energy gaps for removal or insertion of the acidic proton and the computed solvation free energy of the proton. The unsigned mean error relative to experimental results is 2.1 pKa units with a maximum error of 4.0 pKa units. The estimated mean statistical uncertainty due to the finite length of the trajectories is ±1.1 pKa units. The solvation structures of the protonated and deprotonated amino acids are analyzed in terms of radial distribution functions, which can serve as reference data for future force field developments.
Collapse
Affiliation(s)
- Martina Mangold
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Leslie Rolland
- Departement de Chimie, Ecole Normale Superieure , 24 rue Lhomond, 75231 Paris Cedex 05, France
| | - Francesca Costanzo
- Dipartimento di Chimica Fisica e Inorganica, Universita di Bologna , Viale Risorgimento 4, I-40136 Bologna, Italy
| | - Michiel Sprik
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Marialore Sulpizi
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London , London WC1E 6BT, United Kingdom
| |
Collapse
|
8
|
Maurer P, Thomas V, Iftimie R. A computational study of ultrafast acid dissociation and acid–base neutralization reactions. II. The relationship between the coordination state of solvent molecules and concerted versus sequential acid dissociation. J Chem Phys 2011; 134:094505. [DOI: 10.1063/1.3554654] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
|
9
|
Maurer P, Thomas V, Rivard U, Iftimie R. A computational study of ultrafast acid dissociation and acid-base neutralization reactions. I. The model. J Chem Phys 2010; 133:044108. [PMID: 20687634 DOI: 10.1063/1.3461162] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Ultrafast, time-resolved investigations of acid-base neutralization reactions have recently been performed using systems containing the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS) and various Bronsted bases. Two conflicting neutralization mechanisms have been formulated by Mohammed et al. [Science 310, 83 (2005)] and Siwick et al. [J. Am. Chem. Soc. 129, 13412 (2007)] for the same acid-base system. Herein an ab initio molecular dynamics based computational model is formulated, which is able to investigate the validity of the proposed mechanisms in the general context of ground-state acid-base neutralization reactions. Our approach consists of using 2,4,6-tricyanophenol (exp. pKa congruent with 1) as a model for excited-state HPTS( *) (pKa congruent with 1.4) and carboxylate ions for the accepting base. We employ our recently proposed dipole-field/quantum mechanics (QM) treatment [P. Maurer and R. Iftimie, J. Chem. Phys. 132, 074112 (2010)] of the proton donor and acceptor molecules. This approach allows one to tune the free energy of neutralization to any desired value as well as model initial nonequilibrium hydration effects caused by a sudden increase in acidity, making it possible to achieve a more realistic comparison with experimental data than could be obtained via a full-QM treatment of the entire system. It is demonstrated that the dipole-field/QM model reproduces correctly key properties of the 2,4,6-tricyanophenol acid molecule including gas-phase proton dissociation energies and dipole moments, and condensed-phase hydration structure and pKa values.
Collapse
Affiliation(s)
- Patrick Maurer
- Département de Chimie, Université de Montréal, CP 6128, succursale Centre-Ville, Canada
| | | | | | | |
Collapse
|