Rigorous pKa Estimation of Amine Species Using Density-Functional

Dec 5, 2017 - Stepping up gene synthesis. Microarray synthesis of DNA oligonucleotides costs a fraction of a penny per base. But such methods can't...
2 downloads 10 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Rigorous pKa Estimation of Amine Species Using DensityFunctional Tight-Binding-Based Metadynamics Simulations Aditya Wibawa Sakti, Yoshifumi Nishimura, and Hiromi Nakai J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.7b00855 • Publication Date (Web): 05 Dec 2017 Downloaded from http://pubs.acs.org on December 11, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Chemical Theory and Computation is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Rigorous pKa Estimation of Amine Species Using Density-Functional Tight-Binding-Based Metadynamics Simulations

Aditya Wibawa Sakti,a Yoshifumi Nishimura,b and Hiromi Nakaia,b,c,d,*

a

Department of Chemistry and Biochemistry, School of Advanced Science and Engineering,

Waseda University, Tokyo 169-8555, Japan b

c

Research Institute for Science and Engineering, Waseda University, Tokyo 169-8555, Japan

Core Research for Evolutional Science and Technology (CREST), Japan Science and

Technology Agency (JST), Saitama 332-0012, Japan d

Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyoto 615-

8520, Japan

1 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT Predicting pKa values for different types of amine species with high accuracy and efficiency is of critical importance for the design of high performance and economical solvents in carbon capture and storage with aqueous amine solutions. In this study, we demonstrate that density-functional tight-binding-based metadynamics simulations are a promising approach to calculate the free energy difference between the protonated and neutral states of amines in aqueous solution with inexpensive computational cost. The calculated pKa values were in satisfactory agreement with the experimental values, the mean absolute deviation being only 0.08 pKa units for 34 amines commonly used in CO2 scrubbing. Such superior reproducibility and correlation compared to estimations by static quantum mechanical calculations highlight the significant effect of dynamical proton transfer processes in explicit solvent molecules for the improvement of the estimation accuracy.

2 Environment ACS Paragon Plus

Page 2 of 30

Page 3 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

1. INTRODUCTION In the present era, the increasing abundance of CO2 in the atmosphere is arguably leading to global warming.1 To reduce the CO2 levels emitted by fossil fuel combustion sources, consideration efforts have been devoted to the development of carbon capture and storage (CCS) systems.2–4 Aqueous amine solutions play a central role as advanced materials in CCS by chemical absorption.5,6 The pKa values of the conjugate acids of amines are considered to be an indicator of absorption efficiency.7,8 Experimental measurements have thus been performed for a wide range of amine species.8–11 On the theoretical side, quantitative evaluation of pKa values is still challenging. The routinely adopted procedure is a combination of thermodynamic cycles and free energy calculations of associated molecules in individual steps.12,13 Since the pioneering study by da Silva and Svendsen,14 a great deal of attention has been paid to the pKa calculation of amines, in which solvent effects are included implicitly and often explicitly.15–21 In these early theoretical studies, a mean absolute deviation (MAD) of 0.2–0.4 pKa units relative to the selected experimental values was considered an elaborated achievement. Unfortunately, static theoretical calculations are unable to process the various possible solvation patterns in deprotonation reactions. The delicate structural changes and proton transfer processes involved in acid dissociation reactions can be favorably described by full quantum mechanical molecular dynamics (MD) simulations. Density functional theory (DFT)-based MD simulations have provided not only mechanistic insight22–24 but also the ability to predict pKa values25–33 for typical acid–base reactions. Most of those investigations have adopted the metadynamics sampling technique34–36 to study such infrequent deprotonation events within accessible simulation times. The free energy difference between the protonated and deprotonated

3 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

species calculated from reconstructed free energy profiles in metadynamics has shown an error of approximately 0.2 pKa units relative to the experimentally determined values.26–29 The density-functional tight-binding (DFTB)37–39 method is a highly efficient approximation to the Kohn–Shan DFT, which uses atomic and diatomic parameters obtained from DFT calculations. Computational speeds 2–3 orders magnitude faster than that of conventional DFT implementation have been achieved.37,39 Furthermore, the DFTB method retains the robustness of the chemical reaction dynamics by virtue of the self-consistent determination of charge distributions.40 The authors’ group has implemented an efficient program of the DFTB-MD simulation,41 especially, by combining with the linear scaling divideand-conquer technique.42–45 The successes of describing several chemical phenomena in aqueous solutions46–48 have stimulated us to investigate whether the DFTB-based simulations can become an alternative for predicting pKa values of amines that enables fast screening for new compounds. In the present study, we have estimated the pKa values of amine species in solution using metadynamics simulations at the DFTB level. Although the combination of the DFTB method and the metadynamics technique has attracted recent interest for application in complex environments and yielded encouraging results,49–51 evaluation of DFTB accuracy for various types of amine species needs to be addressed for prediction purpose. To this end, the calculated pKa values are compared with available experimental data for 34 amine species. We also examine the effect of system size to the quality of theoretical predictions, which is hard to obtain from computationally more demanding DFT-based MD simulations. The paper is organized as follows. The computational details of the DFTB metadynamics simulations are provided in Section 2. Section 3 presents the accuracy and efficiency of the investigated approach. Finally, concluding remarks are described in Section 4.

4 Environment ACS Paragon Plus

Page 4 of 30

Page 5 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

2. COMPUTATIONAL DETAILS Canonical ensemble simulations were performed for a cubic cell with three-dimensional periodic boundary conditions unless noted otherwise. In the initial condition, the system consisted of one protonated amine, one hydroxide ion, and 50 water molecules. For monoethanolamine (MEA), three larger simulation boxes containing 100, 200, and 500 water molecules were additionally prepared to investigate the system size dependence on the calculated pKa values. The density of each system was determined from classical MD simulations using the general CHARMM force field52 and NAMD package.53 This pre-equilibration procedure was performed for 3 ns with a time step of 1.0 fs under NVE, NVT, and NPT ensembles. The classically equilibrated structure was further equilibrated at the third-order DFTB (DFTB3)54 level with 3ob parameter set55 for 50 ps with a time step of 1.0 fs. The choice of DFTB model and parameter was based on the best estimation of pKa values for MEA at two different temperatures among several DFTB variants considered in preliminary simulations. During the DFTB equilibration, the temperature was controlled using an Andersen thermostat at 298.15 or 313.15 K. The structure and velocity of last equilibration step were adopted as initial conditions of the DFTB metadynamics simulations. In the production run, the time step was 0.5 fs. All DFTB calculations were performed using the DC-DFTB-K program.41 In order to observe the deprotonation reaction given by [Amine(H)]+ + H2O

Amine + H3O+,

(1)

the number of hydrogen atoms coordinated to the nitrogen atom of the amine species (nNH) was chosen as the collective variable (CV) for the DFTB metadynamics simulations. In accordance with previous successful free energy samplings,26–29,56 the rational form of nNH was adopted:

5 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

( )6 ∑ 1 − (rNH r0 )12 i =1 , NH

nNH =

Page 6 of 30

1 − rNHi r0

(2)

i

where rNH is the distance between the nitrogen and hydrogen atoms, r0 is the cut-off distance set i

to 1.6 Å, and the summation runs over the NH hydrogen atoms initially coordinated to the nitrogen atom. In the DFTB metadynamics simulations, the height and width of the Gaussian function to describe the bias potential were set to 1.31 kJ/mol and 0.1, respectively. Note that the CV is dimensionless. The Gaussian bias potential was deposited at time intervals of 0.04 ps. After reaching the convergence of free energy profile, the pKa value at temperature T was calculated using the following relationship: pKa =

∆G 2.303RT ,

(3)

where R is the universal gas constant and ∆G is the free energy difference between the two states in Eq. (1), directly calculated from the sampled free energy profile.

6 Environment ACS Paragon Plus

Page 7 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

3. RESULTS AND DISCUSSION 3.1. pKa Estimation of Aqueous MEA Solution First, we discuss the accuracy and efficiency of DFTB metadynamics simulations taking MEA, which is one of the most commonly used amines in CO2 capture process,5 as representative example. The time-course change of the CV during the deprotonation process of protonated MEA (MEA(P)) at 313.15 K is shown in Figure 1a. The nNH values of ~1.9 and ~2.8 correspond to MEA and MEA(P), respectively. A proton is shuttling between MEA(P) and an adjacent water molecule at 3.5–4.8 ps. The fluctuation at 5.2–5.9 ps is a manifestation of the instantaneous formation of a MEA(P)–water pair. The CV starts to be diffusive from ∼6.5 ps, at which the metadynamics sampling was stopped. In the diffusive region, the proton dissociated from MEA(P) migrates through the hydrogen bond network of water molecules, the so-called Grotthuss shuttling.57,58 The relatively short simulation time to complete the reaction in Eq. (1) is a consequence of the sufficient equilibration achieved before the production run.26–29 Figure 1b shows the time-course change of the free energy difference between the protonated and neutral states of MEA at 313.15 K. Sampling of the diffusive region hardly changes the free energy difference, demonstrating that the reconstructed free energy profile is well converged. Figure 1c shows the free energy profile for the deprotonation of MEA(P) at 313.15 K. Two minima are found at nNH = ~1.9 and ~2.8, representing the MEA and MEA(P) states, respectively. The calculated free energy difference is 54.19 kJ/mol, equivalent to a pKa value of 9.15. This estimation is in good agreement with the experimental value of 9.09.10 We also performed the DFTB metadynamics simulations of MEA at different conditions. One is to treat different unit cell sizes containing a larger number of water molecules. The time-

7 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

course changes of CV and calculated free energy profiles for 100, 200, and 500 water molecules are presented in Figures S1–S3 in the Supporting Information. The estimated pKa values of 100, 200, and 500 water molecules system are 9.19, 9.03, and 9.07, respectively. Considering the insensitivity of calculated pKa value with respect to the system size, 50 water molecules to solvate amine species seems to be sufficient for prediction purpose and was used throughout this study. The other is to calculate the pKa value at 298.15 K (for corresponding simulation results, see Figure S4 in the Supporting Information). The agreement of estimation (9.40) with experimental data (9.47)10 is again excellent, which indicates that the metadynamics sampling is able to capture the effect of slight change in temperature reasonably. The efficiency of the DFTB metadynamics simulation is remarkable. Calculation of a 60 ps long trajectory took less than an hour using an 8-core Intel Xeon E5-2637 v3 (3.50 GHz) workstation. The results demonstrate that the region shown in Figure 1 can be generated within minutes. In contrast, previous DFT-based metadynamics simulations of systems containing a similar number of solvent water molecules have been reported to require more expensive computational resources (i.e., an IBM Blue Gene cluster with 256 processors).26,29

3.2. pKa Estimation of Various Types of Amines The obtained results in Section 3.1 give us confidence to evaluate the performance of DFTB metadynamics simulations for amines being in common use in the CCS technology. Specifically, the number of target amines was extended to 34 in total, including seven primary (10), nine secondary (20), nine tertiary (30), four cyclic secondary (c-20), and five cyclic tertiary (c-30) species. A comparison of the calculated pKa values with the experimental results8–11 is shown in Table 1. The time-course changes of the CV and the simulated free energy profiles are given

8 Environment ACS Paragon Plus

Page 8 of 30

Page 9 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Figures S5–S37 in the Supporting Information, except for MEA at 313.15 K, already shown in Figure 1. In all cases, sampling of reactive region to fulfill the free energy profile with bias potential was continued no longer than 12 ps and protonated and neutral states were clearly distinguished from each other in the free energy profile. Overall consistency was observed between the estimated and experimental pKa values. Moreover, the obtained MAD was only 0.09 pKa units. The largest error was 0.21 pKa units (1.20 kJ/mol) for n-cyclohexylethanolamine. It should be noted that this specific group of amines does not present an obvious trend toward a certain range of pKa values. The averaged statistical error bar based on five MD trajectories for each amine is 0.26 pKa units, with the maximum deviation of 0.53 pKa units for 3-piperidino-1,2propanediol. More details of the estimated pKa values for individual metadynamics runs are given in Table S3 in the Supporting Information. The correlation of pKa values in Table 1 is graphically displayed in Figure 2. An excellent correlation coefficient of R2 = 0.990 was obtained for the 34 amines, with pKa values ranging from 6.6 to 11.1. For comparison with previous studies using static theoretical calculations, see Figures S38–S40 in the Supporting Information. The outperformance of the present approach strongly supports the need to account for the multiple arrangements of the solvent water molecules and the forward/backward proton shuttling processes in order to accurately determine pKa values.

9 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

4. CONCLUSION In conclusion, DFTB-based metadynamics simulations are able to estimate pKa values of amine species in aqueous solution with remarkable accuracy. The calculated MAD for 34 amines commonly used in CO2 capture was only 0.09 pKa units with respect to the experimental results, superior by a factor of 2–3 to common static theoretical calculations. The observed accuracy may be attributed to the inclusion of Grotthuss shuttling effects and the sampling of solvation shells around the reaction centers during the deprotonation process. In contrast, the single-point calculations at the DFTB3/3OB level underestimate the gas-phase acidities and proton affinities of amines, as shown in Table S1 in the Supporting Information. Incorporation of polarizable continuum model cannot improve the pKa values for single-point calculations, as shown in Table S2 in the Supporting Information. We also found that the pKa values obtained with DFTB metadynamics simulations are less sensitive with respect to the system size while temperature effect can be incorporated adequately. Although possible error sources exist, such as metadynamics sampling errors, DFTB model/parameter dependence, and the neglect quantum effects, those negative factors seem to be cancelled out in the resultant free energy difference results.26–29,59 The required computational demands are small with no deterioration of the prediction accuracy; hence, the present approach presents a great advantage over DFT- and/or wavefunction-based protocols in that prospective candidates for efficient chemical absorption of CO2 can be explored in a feasible manner. The accurate determination of acid dissociation constants for different groups of compounds may be feasible after careful inspection of the DFTB performance.

10 Environment ACS Paragon Plus

Page 10 of 30

Page 11 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

ACKNOWLEDGMENT This study was supported in part by a Grant-in-Aid for Scientific Research (A) “KAKENHI Grant Number JP26248009” from the Japan Society for the Promotion of Science (JSPS), a Grand-in-Aid for Challenging and Exploratory Research “KAKENHI 15K13629” from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT), Japan, and by MEXT as “Priority Issue on Post-K computer” (Development of new fundamental technologies for highefficiency energy creation, conversion/storage, and use). One of the authors (A.W.S.) acknowledges financial support from the Yoshida Scholarship Foundation (YSF). Some of the calculations were performed at the Research Center for Computational Science (RCCS), Okazaki, Japan.

ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publication website. Time-course changes of collective variable and free energy profile for deprotonation reaction of 34 amine species obtained from production run of DFTB metadynamics simulations and correlations between estimated pKa values obtained from previous static theoretical calculations and the experimental results.

11 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR INFORMATION Corresponding Author *Corresponding author: H. Nakai, Phone: +81-3-5286-3452, E-mail: [email protected]

Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

Notes The authors declare no competing financial interest.

12 Environment ACS Paragon Plus

Page 12 of 30

Page 13 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

REFERENCES (1) Stern, N. Stern Review on the Economics of Climate Change; Cambridge University Press: Cambridge, 2006. (2) D’Allessandro, D. M.; Smit, B.; Long, J. R. Carbon Dioxide Capture: Prospects for New Materials. Angew. Chem., Int. Ed. 2010, 49, 6058–6082. (3) MacDowell, N.; Florin, N.; Buchard, A.; Hallett, J.; Galindo, A.; Jackson, G.; Adjiman, C. S.; Williams, C. K.; Shah, N.; Fennell, P. An Overview of CO2 Capture Technologies. Energy Environ. Sci. 2010, 3, 1645–1669. (4) Kenarsari, S. D.; Yang, D.; Jiang, G.; Zhang, S.; Wang, J.; Russell, A. G.; Wei, Q.; Fan, M. Review of Recent Advances in Carbon Dioxide Separation and Capture. RSC Adv. 2013, 3, 22739–22773. (5) Rochelle, G. T. Amine Scrubbing for CO2 Capture. Science 2009, 325, 1652–1654. (6) Mumford, K. A.; Wu, Y.; Smith, K. H.; Stevens, G. W. Review of Solvent Based Carbon-Dioxide Capture Technologies. Front. Chem. Sci. Eng. 2015, 9, 125–141. (7) Versteeg, G. F.; van Dijck, L. A. J.; van Swaaij, W. P. M. On the Kinetics between CO2 and Alkanolamines both in Aqueous and Non-Aqueous Solutions. An Overview. Chem. Eng. Commun. 1996, 144, 113–158. (8) Puxty, G.; Rowland, R.; Allport, A.; Yang, Q.; Bown, M.; Burns, R.; Maeder, M.; Attalla, M. Carbon Dioxide Postcombustion Capture: A Novel Screening Study of the Carbon Dioxide Absorption Performance of 76 Amines. Environ. Sci. Technol. 2009, 43, 6427–6433.

13 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 30

(9) Rayer, A. V.; Sumon, K. Z.; Jaffari, L.; Henni, A. Dissociation Constants (pKa) of Tertiary and Cyclic Amines: Structural and Temperature Dependences. J. Chem. Eng. Data 2014, 59, 3805–3813. (10) Tagiuri,

A.;

Mohamedali, M.;

Henni,

A.

Dissociation Constant

(pKa) and

Thermodynamic Properties of Some Tertiary and Cyclic Amines from (298 to 333) K. J. Chem. Eng. Data 2016, 61, 247–254. (11) Hamborg, E. S.; Versteeg, G. F. Dissociation Constants and Thermodynamic Properties of Amines and Alkanolamines from (293 to 353) K. J. Chem. Eng. Data 2009, 54, 1318– 1328. (12) Alongi, K. S.; Shields, G. C. Theoretical Calculations of Acid Dissociation Constants: A Review Article. Annu. Rep. Comput. Chem. 2010, 6, 113–138. (13) Ho, J.; Coote, M. L. A Universal Approach for Continuum Solvent pKa Calculations: Are We There Yet? Theor. Chem. Acc. 2010, 125, 3–21. (14) Da Silva, E. F.; Svendsen, H. F. Prediction of the pKa Values of Amines Using Ab Initio Methods and Free-Energy Perturbations. Ind. Eng. Chem. Res. 2003, 42, 4414–4421. (15) Bryantsev, V. S.; Diallo, M. S.; Goddard, W. A., III pKa Calculations of Aliphatic Amines, Diamines, and Aminoamides via Density Functional Theory with a PoissonBoltzmann Continuum Solvent Model. J. Phys. Chem. A 2007, 111, 4422–4430. (16) Khalili, F.; Henni, A.; East, A. L. L. Entropy Contributions in pKa Computation: Application to Alkanolamines and Piperazines. J. Mol. Struct.: THEOCHEM 2009, 916, 1–9. (17) Yamada, H.; Shimizu, S.; Okabe, H.; Matsuzaki, Y.; Chowdhury, F. A.; Fujioka, Y. Prediction of the Basicity of Aqueous Amine Solutions and the Species Distribution in

14 Environment ACS Paragon Plus

Page 15 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

the Amine-H2O-CO2 System Using the COSMO-RS Method. Ind. Eng. Chem. Res. 2010, 49, 2449–2455. (18) Gupta, M.; da Silva, E. F.; Svendsen, H. F. Modeling Temperature Dependency of Amine Basicity Using PCM and SM8T Implicit Solvation Models. J. Phys. Chem. B 2012, 116, 1865–1875. (19) Zhang, S. A Reliable and Efficient First Principles-Based Method for Predicting pKa Values. 4. Organic Bases. J. Comput. Chem. 2012, 33, 2469–2482. (20) Gangarapu, S.; Marcelis, A. T. M.; Zuilhof, H. Accurate pKa Calculation of the Conjugate Acids of Alkanolamines, Alkaloids and Nucleotide Bases by Quantum Chemical Methods. ChemPhysChem 2013, 14, 990–995. (21) Gupta, M.; da Silva, E. F.; Svendsen, H. F. Postcombustion CO2 Capture Solvent Characterization Employing the Explicit Solvation Shell Model and Continuum Solvation Models. J. Phys. Chem. B 2016, 120, 9034–9050. (22) Park, J. M.; Laio, A.; Iannuzzi, M.; Parrinello, M. Dissociation Mechanism of Acetic Acid in Water. J. Am. Chem. Soc. 2006, 128, 11318–11319. (23) Lee, J.-G.; Asciutto, E.; Babin, V.; Sagui, C.; Darden, T.; Roland, C. Deprotonation of Solvated Formic Acid: Car-Parrinello and Metadynamics Simulations. J. Phys. Chem. B 2006, 110, 2325–2331. (24) Galib, M.; Hanna, G. Mechanistic Insights into the Dissociation and Decomposition of Carbonic Acid in Water via the Hydroxide Route: An Ab Initio Metadynamics Study. J. Phys. Chem. B 2011, 115, 15024–15035.

15 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(25) Ivanov, I.; Chen, B.; Raugei, S.; Klein, M. L. Relative pKa Values from First-Principles Molecular Dynamics: The Case of Histidine Deprotonation. J. Phys. Chem. B 2006, 110, 6365–6371. (26) Tummanapelli, A. K.; Vasudevan, S. Dissociation Constants of Weak Acids from Ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pKa Values. J. Phys. Chem. B 2014, 118, 13651–13657. (27) Tummanapelli, A. K.; Vasudevan, S. Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pKa Values from Metadynamics Sampling. J. Phys. Chem. B 2015, 119, 12249–12255. (28) Tummanapelli, A. K.; Vasudevan, S. Ab Initio MD Simulations of the Brønsted Acidity of Glutathione in Aqueous Solutions: Predicting pKa Shifts of the Cysteine Residue. J. Phys. Chem. B 2015, 119, 15353–15358. (29) Tummanapelli, A. K.; Vasudevan, S. Estimating Successive pKa Values of Polyprotic Acids from Ab Initio Molecular Dynamics Using Metadynamics: The Dissociation of Phthalic Acid and Its Isomers. Phys. Chem. Chem. Phys. 2015, 17, 6383–6388. (30) Leung, K.; Nielsen, I. M. B.; Criscenti, L. J. Elucidating the Bimodal Acid-Base Behavior of the Water-Silica Interface from First Principles. J. Am. Chem. Soc. 2009, 131, 18358–18365. (31) Sulpizi, M.; Sprik, M. Acidity Constants from DFT-Based Molecular Dynamics Simulations. J. Phys.: Condens. Matter 2010, 22, 284116. (32) 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–1961.

16 Environment ACS Paragon Plus

Page 16 of 30

Page 17 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

(33) Gaigeot, M.-P.; Sprik, M.; Sulpizi, M. J. Phys.: Condens. Matter Oxide/Water Interfaces: How the Surface Chemistry Modifies Interfacial Water Properties. 2012, 24, 124106. (34) Laio, A.; Parrinello, M. Escaping Free-Energy Minima. Proc. Natl. Acad. Sci. U. S. A. 2002, 99, 12562–12566. (35) Ensing, B.; De Vivo, M.; Liu, Z.; Moore, P.; Klein, M. L. Metadynamics as a Tool for Exploring Free Energy Landscapes of Chemical Reactions. Acc. Chem. Res. 2006, 39, 73–81. (36) Barducci, A.; Bonomi, M.; Parrinello, M. Metadynamics. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2011, 1, 826–843. (37) Seifert, G.; Joswig, J.-O. Density-Functional Tight Binding–An Approximate DensityFunctional Theory Method. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 456–465. (38) Elstner, M.; Seifert, G. Density Functional Tight Binding. Philos. Trans. R. Soc., A 2014, 372, 20120483. (39) Gaus, M.; Cui, Q.; Elstner, M. Density Functional Tight Binding: Application to Organic and Biological Molecules. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2014, 4, 49–61. (40) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-Consistent-Charge Density-Functional Tight-Binding Method for Simulations of Complex Materials Properties. Phys. Rev. B 1998, 58, 7260–7268. (41) Nishizawa, H.; Nishimura, Y.; Kobayashi, M.; Irle, S.; Nakai, H. Three Pillars for Achieving Quantum Mechanical Molecular Dynamics Simulations of Huge Systems: Divide-and-Conquer,

Density-Functional

Tight-Binding,

Computation. J. Comput. Chem. 2016, 37, 1983–1992.

17 Environment ACS Paragon Plus

and

Massively

Parallel

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(42) Yang, W.; Lee, T.-S. Density-Matrix Divide-and-Conquer Approach for Electronic Structure Calculations of Large Molecules. J. Chem. Phys. 1995, 103, 5674–5678. (43) Akama, T.; Kobayashi, M.; Nakai, H. Implementation of Divide-and-Conquer Method Including Hartree-Fock Exchange Interaction. J. Comput. Chem. 2007, 28, 2003–2012. (44) Kobayashi, M.; Nakai, H. In Linear-Scaling Techniques in Computational Chemistry and Physics; Springer: Dordrecht, 2011; pp 97–127. (45) Kobayashi, M.; Nakai, H. How Does It Become Possible to Treat Delocalized and/or Open-Shell Systems in Fragmentation-Based Linear-Scaling Electronic Structure Calculations? The Case of the Divide-and-Conquer Method. Phys. Chem. Chem. Phys. 2012, 14, 7629–7639. (46) Nakai, H.; Nishimura, Y.; Kaiho, T.; Kubota, T.; Sato, H. Contrasting Mechanisms for CO2 Absorption and Regeneration Processes in Aqueous Amine Solutions: Insights from Density-Functional Tight-Binding Molecular Dynamics Simulations. Chem. Phys. Lett. 2016, 647, 127–131. (47) Nakai, H.; Sakti, A. W.; Nishimura, Y. Divide-and-Conquer-Type Density-Functional Tight-Binding Molecular Dynamics Simulations of Proton Diffusion in a Bulk Water System. J. Phys. Chem. B 2016, 120, 217–221. (48) Sakti, A. W.; Nishimura, Y.; Nakai, H. Divide-and-Conquer-Type Density-Functional Tight-Binding Simulations of Hydroxide Ion Diffusion in Bulk Water. J. Phys. Chem. B 2017, 121, 1362–1371. (49) Selli, D.; Baburin, I. A.; Martoňák, R.; Leoni, S. Novel Metastable Metallic and Semiconducting Germaniums. Sci. Rep. 2013, 3, 1466.

18 Environment ACS Paragon Plus

Page 18 of 30

Page 19 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

(50) Knaup, J. M.; Wehlau, M.; Frauenheim, T. Permutation-Invariant Collective Variable to Track and Drive Vacancy Dynamics in Simulations of Solids. Phys. Rev. B 2013, 88, 220101. (51) Cuny, J.; Korchagina, K.; Menakbi, C.; Mineva, T. Metadynamics Combined with Auxiliary Density Functional and Density Functional Tight-Binding Methods: Alanine Dipeptide as a Case Study. J. Mol. Model. 2017, 23, 72. (52) Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; Mackerell, A. D., Jr. CHARMM General Force Field: A Force Field for Drug-Like Molecules Compatible with the CHARMM All-Atom Additive Biological Force Fields. J. Comput. Chem. 2009, 31, 671–690. (53) Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kalé, L; Schulten, K. Scalable Molecular Dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. (54) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge DensityFunctional Tight-Binding Method (SCC-DFTB). J. Chem. Theory Comput. 2011, 7, 931– 948. (55) Gaus, M.; Goez, A.; Elstner, M. Parametrization and Benchmark of DFTB3 for Organic Molecules. J. Chem. Theory Comput. 2013, 9, 338–354. (56) Iannuzzi, M.; Laio, A.; Parrinello, M. Efficient Exploration of Reactive Potential Energy Surfaces Using Car-Parrinello Molecular Dynamics. Phys. Rev. Lett. 2003, 90, 238302. (57) Agmon, N. The Grotthuss Mechanism. Chem. Phys. Lett. 1995, 244, 456–462. (58) Marx, D. Proton Transfer 200 Years after von Grotthuss: Insights from Ab Initio Simulations. ChemPhysChem 2006, 7, 1848–1870.

19 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(59) Striling, A.; Pápai, I. H2CO3 Forms via HCO3– in Water. J. Phys. Chem. B 2010, 114, 16854–16859.

20 Environment ACS Paragon Plus

Page 20 of 30

Page 21 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

TABLES Table 1. Comparison of the pKa values calculated from DFTB metadynamics simulations with the experimental results for 34 amine species in aqueous solution Entry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 MAD

Amine Monoethanolamine 2-Amino-2-methyl-1-propanol 2-(2-Aminoethoxy)ethanol 2-Amino-2-ethyl-1,3-propanediol 3-Amino-1-propanol Serinol(2-aminopropane-1,3-diol) Tris(hydroxymethyl)aminomethane Methylmonoethanolamine 2-(Ethylamino)ethanol 2-(tert-Butylamino)ethanol Diisopropanolamine n-Cyclohexylethanolamine Diethanolamine n-Cyclopentylethanolamine Tricine n-Cyclopropylethanolamine Methyldiethanolamine 2-(Diisopropylamino)ethanol Triethanolamine Tris[2-(2-methoxyethoxy)ethyl]amine Triethylamine 3-Dimethylamino-1-propanol Ethyldiethanolamine N,N-Dimethylisopropanolamine tert-Butyldiethanolamine Piperazine 2-Piperidinemethanol 2-Piperidineethanol 3-Piperidinemethanol 1-Piperidinepropanolamine 3-Piperidino-1,2-propanediol 3-Quinoclidinol 4,2-Hydroxyethylmorpholine 3-Morpholinopropylamine

Abbreviation MEA AMP 2-AEE AEPD MPA SAPD THMAM MMEA EAE TBAE DIPA n-CHEA DEA n-CPEA TRC n-CPREA MDEA 2-DIPA TEA TMEEA TREA 3-DMAP EDEA DMIPA t-BDEA PZ 2-PPM 2-PPE 3-PPM 1-PPP 3-PPPD 3-QCD 4,2-HEMO 3-MOPA

Typea 10 10 10 10 10 10 10 20 20 20 20 20 20 20 20 20 30 30 30 30 30 30 30 30 30 c-20 c-20 c-20 c-20 c-30 c-30 c-30 c-30 c-30

21 Environment ACS Paragon Plus

T /K 313.15 313.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 313.15 313.15 313.15 313.15 313.15 298.15 298.15 298.15 298.15 313.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15

pKa (Calc.) 9.14 9.16 9.49 8.78 9.86 8.58 8.16 9.94 10.00 9.76 8.80 9.89 8.80 10.25 8.00 8.50 8.43 9.27 7.41 6.67 10.30 9.30 8.66 9.45 8.92 9.30 10.47 11.08 10.98 9.32 8.79 10.22 6.82 9.89 0.09

pKa (Exptl.) 9.09b 9.23c 9.42c 8.82c 9.96c 8.55d 8.08d 9.85c 10.00d 9.70d 8.88d 10.10d 8.88d 10.10d 8.10d 8.40d 8.31b 9.11b 7.45c 6.69b 10.32c 9.27d 8.80e 9.47e 9.03e 9.38b 10.60d 10.90d 10.80d 9.43d 8.91d 10.10d 6.90e 9.95e

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

a

Primary, secondary, tertiary, cyclic secondary, and cyclic tertiary amines are denoted as 10, 20,

30, c-20, and c-30, respectively. bRef. 10. cRef. 11. dRef. 8. eRef. 9. FIGURE CAPTIONS Figure 1. Time-course changes of the (a) CV, nNH, and (b) free energy difference, and (c) the free energy profile as a function of the CV for the deprotonation of MEA(P) at 313.15 K.

Figure 2. Correlation between the calculated pKa values from DFTB metadynamics simulations and the experimental results8–11 for 34 amine species in aqueous solution.

22 Environment ACS Paragon Plus

Page 22 of 30

Page 23 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

FIGURES

(a)

(b)

(c)

Figure 1. Time-course changes of the (a) CV, nNH, and (b) free energy difference, and (c) the free energy profile as a function of the CV for the deprotonation of MEA(P) at 313.15 K.

23 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. Correlation between the calculated pKa values from DFTB metadynamics simulations and the experimental results8–11 for 34 amine species in aqueous solution.

24 Environment ACS Paragon Plus

Page 24 of 30

Page 25 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

TABLE OF CONTENTS GRAPHICS

25 Environment ACS Paragon Plus

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Fig. 1a 85x59mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Fig. 1b 85x59mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Fig. 1c 88x62mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Fig. 2 635x635mm (72 x 72 DPI)

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table of contents 338x190mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 30 of 30