1
|
Yan S, Wang B, Lin H. Reshaping the QM Region On-the-Fly: Adaptive-Shape QM/MM Dynamic Simulations of a Hydrated Proton in Bulk Water. J Chem Theory Comput 2024; 20:3462-3472. [PMID: 38671391 DOI: 10.1021/acs.jctc.4c00164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2024]
Abstract
Adaptive quantum mechanics/molecular mechanics (QM/MM) reclassifies on-the-fly a molecule or molecular fragment as QM or MM during dynamics simulations without abrupt changes in the energy or forces. Notably, the permuted adaptive-partitioning (PAP) algorithms have been applied to simulate a hydrated proton, with a mobile QM zone anchored at a pseudoatom called a proton indicator. The position of the proton indicator approximates the location of the delocalized excess proton, yielding a smooth trajectory of the proton diffusing via the Grotthuss mechanism in aqueous solutions. The mobile QM zone, which has been taken to be a sphere with a preset radius, follows the proton wherever it goes. Although the simulations are successful, the use of a spherical QM zone has one disadvantage: A large preset radius must be utilized to minimize the chance of missing water molecules that are important to proton translocation. A large radius leads to a large QM zone, which is computationally expensive. In this work, we report a new way to set up the QM zone, where one includes only the water molecules important to proton transfer. The importance of a given water molecule is quantified by its "weight" that depends on its relation to the reaction path of proton transfer. The weight varies smoothly, ensuring that a water molecule gradually appears in or disappears from the QM zone without abrupt changes, as required by the PAP method. Consequently, the shape of the QM zone evolves on-the-fly, keeping the QM zone as small as possible and as large as necessary. Test simulations demonstrate that the new algorithm significantly improves the computation efficiency while maintaining the proper descriptions of proton transfer in bulk water.
Collapse
Affiliation(s)
- Shengheng Yan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| | - Hai Lin
- Department of Chemistry, CB 194, University of Colorado Denver, Denver, P.O. Box 173364, Colorado 80217, United States
| |
Collapse
|
2
|
Luan B, McDonagh JL. Developing semi-empirical water model for efficiently simulating temperature-dependent chemisorption of CO 2 in amine solvents. Phys Chem Chem Phys 2024; 26:3540-3547. [PMID: 38214052 DOI: 10.1039/d3cp05874c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
Classical molecular dynamics (MD) simulations without bond forming/breaking cannot be used to model chemical reactions (CRs) among small molecules. Although the first-principle MD simulation can adequately describe CRs with explicit water molecules, such simulation is normally too costly for most researchers to afford. Generally, water molecules in a solvent can exert hydrophobic forces on reacting molecules, which yields a so-called caging effect that cannot be ignored when constructing a free energy landscape for reacting molecules. Many recently developed semi-empirical methods (such as DFTB, PM6 and xTB) are highly efficient for modeling CRs, however none of them can be directly used to model bulk water properly. Here, we developed a modified xTB approach that enables the simulation of CRs in explicit water. Using the chemisorption of CO2 by amines in water as an example application, we demonstrate that our approach yielded results comparable with the first-principle ones, while only using a limited computing resource. Potentially, our proposed semi-empirical water model can be utilized for the computational study of any CR in water.
Collapse
Affiliation(s)
- Binquan Luan
- IBM Thomas J. Watson Research, Yorktown Heights, NY 10598, USA.
| | - James L McDonagh
- IBM Research Europe, Hartree Centre, SciTech Daresbury, Warrington, Chesire WA4 4AD, UK
| |
Collapse
|
3
|
Cinq N, Simon A, Louisnard F, Cuny J. Accurate SCC-DFTB Parametrization of Liquid Water with Improved Atomic Charges and Iterative Boltzmann Inversion. J Phys Chem B 2023; 127:7590-7601. [PMID: 37603798 DOI: 10.1021/acs.jpcb.3c03479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
This work presents improvements of the description of liquid water within the self-consistent-charge density-functional based tight-binding scheme combining the use of Weighted Mulliken (WMull) charges and optimized O-H repulsive potential through the iterative Boltzmann inversion (IBI) process. The quality of the newly developed models is validated considering pair radial distribution functions (RDFs), as well as other structural, energetic, thermodynamic, and dynamic properties. The use of WMull charges certainly improves the agreement with experimental data, however leading to over-structured RDFs at short distance, that can be further improved by considering an optimized O-H repulsive potential obtained by the IBI process. Three different schemes were used to optimize this potential: (i) optimization including short O-H distances. This led to accurate RDFs as well as improved self-diffusion coefficient and heat of vaporization, while the proton transfer energy barrier is severely deteriorated; (ii) optimization starting at long distance. The proton transfer energy barrier is recovered while the heat of vaporization is deteriorated and the O-H RDF is less accurate at short distance; (iii) optimization within the path-integral molecular dynamics scheme which allows us to exclude nuclear quantum effects from the repulsive potential. The latter potential, in conjunction with the WMull improved atomic charges, provides similar results as (i) for structural, dynamic, and thermodynamic properties while recovering a large part of the proton transfer energy barrier. It therefore offers a good compromise to study both dynamic properties and chemistry within liquid water at a quantum chemical level.
Collapse
Affiliation(s)
- Nicolas Cinq
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Aude Simon
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Fernand Louisnard
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| |
Collapse
|
4
|
Yan S, Wang B, Lin H. Tracking the Delocalized Proton in Concerted Proton Transfer in Bulk Water. J Chem Theory Comput 2023; 19:448-459. [PMID: 36630655 DOI: 10.1021/acs.jctc.2c01097] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
A solvated proton in water is often characterized as a charge or structural defect, and it is important to track its evolution on-the-fly in certain dynamics simulations. Previously, we introduced the proton indicator, a pseudo-atom, whose position approximates the location of the excess proton modeled as a structural defect. The proton indicator generally yields a smooth trajectory of a hydrated proton diffusing in aqueous solutions, including in the events of stepwise proton transfer via the Grotthuss mechanism; however, the proton indicator did not perform well in the events of concerted proton transfer, for which it occasionally yielded large position displacements between two successive time steps. To overcome this hurdle, we develop a new algorithm of a proton indicator with greatly enhanced performance for concerted proton transfer in bulk water. A protocol is proposed to exhaustively explore the hydrogen-bonding network of the water wires over which the excess proton is delocalized and to properly account for the contributions of the water molecules in this network as the geometry evolves. The new proton indicator (called Indicator 2.0) is assessed in dynamics simulations of an excess proton in bulk water and in specially constructed model systems of more complex architectures. The results demonstrate that the new indicator yields a smooth trajectory in both stepwise and concerted proton transfers.
Collapse
Affiliation(s)
- Shengheng Yan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen360015P. R. China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen360015P. R. China
| | - Hai Lin
- Department of Chemistry, CB 194, University of Colorado Denver, P.O. Box 173364, Denver, Colorado80217, United States
| |
Collapse
|
5
|
Benayad Z, Bova Saint-André M, Stirnemann G. Molecular Mechanisms of Phosphoester Bond Formation in Water Using Tight-Binding Ab Initio Molecular Dynamics. J Phys Chem B 2022; 126:8251-8265. [PMID: 36201374 DOI: 10.1021/acs.jpcb.2c04259] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Phosphate groups are ubiquitous in biomolecules and are usually incorporated through phosphoester bonds between alcohol groups and orthophosphate. The formation of this bond is exceptionally difficult, with associated barriers of 30-45 kcal/mol in the absence of catalysts. In abiotic conditions, polymerizing nucleic acids without enzymes remains very challenging and is still a partly unsolved problem that severely questions the RNA World hypothesis for the origins of life. Offering a solution to this problem would involve a detailed knowledge of the reaction energetics and mechanisms, yet these remain not fully understood at a molecular level, especially because of the very slow reaction rates that represent a significant challenge for the experiments. The number of involved reaction coordinates and the possible role of the solvent in assisting the reaction are challenging for computational studies. Here, we use extensive ab initio molecular dynamics simulations using semiempirical tight-binding methods and enhanced sampling to address these issues. We first show that the choice of the tight-binding method is greatly limited by the instability of the water liquid phase for most DFTB generations and parameter sets that are widely available. We then focus on a model reaction involving methanol and orthophosphate, for which the two protonation states (mono- and dianionic) that are dominant around neutral pH are considered. We compare different proton coordinates that enable (or not) the participation of solvent water molecules. Our simulations suggest that in all cases, a dissociative associative mechanism, with an intermediate metaphosphate, is favored. The main difference between the two phosphate species is that reaction with the monoanion is assisted by the substrate, while that with the dianion involves solvent water molecules. Our results are in agreement with early experimental measurements, but the reaction barriers are underestimated in our framework. We believe that our approach provides an interesting perspective on how to sample the reaction phase space efficiently, but it calls for future studies using more accurate descriptions of chemical reactivity.
Collapse
Affiliation(s)
- Zakarya Benayad
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, PSL University, Université de Paris, 13 rue Pierre et Marie Curie, 75005Paris, France
| | - Matthias Bova Saint-André
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, PSL University, Université de Paris, 13 rue Pierre et Marie Curie, 75005Paris, France
| | - Guillaume Stirnemann
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, PSL University, Université de Paris, 13 rue Pierre et Marie Curie, 75005Paris, France
| |
Collapse
|
6
|
Ignatov I, Huether F, Neshev N, Kiselova-Kaneva Y, Popova TP, Bankova R, Valcheva N, Ignatov AI, Angelcheva M, Angushev I, Baiti S. Research of Water Molecules Cluster Structuring during Haberlea rhodopensis Friv. Hydration. PLANTS (BASEL, SWITZERLAND) 2022; 11:2655. [PMID: 36235522 PMCID: PMC9572004 DOI: 10.3390/plants11192655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/01/2022] [Accepted: 10/04/2022] [Indexed: 06/16/2023]
Abstract
Gesneriaceae plant family is comprised of resurrection species, namely Boea hygrometrica and Paraboea rufescens, that are native to the Southeast Asia and Haberlea rhodopensis, Ramonda myconi, and Ramonda serbica, which are mainly found in the Balkan Peninsula. Haberlea rhodopensis is known to be able to survive extreme and prolonged dehydration. Study was carried out after the dried plant Haberlea rhodopensis Friv. had been hydrated and had reached its fresh state. Two juice samples were collected from the plant blossom: The first sample was prepared with 1% filtered water through a patented EVOdrop device. Then the sample was saturated with hydrogen with EVOdrop booster to a concentration of 1.2 ppm, pH = 7.3, ORP = -390 mV. This first sample was prepared with filtered tap water from Sofia, Bulgaria. The second sample, which was a control one, was developed with tap water from Sofia, Bulgaria, consisting of 1% solutions of Haberlea rhodopensis. A study revealed that during the drying process in H. rhodopensis the number of free water molecules decreases, and water dimers are formed. The aim of our study was to determine the number of water molecules in clusters in 1% solutions of hydrated H. rhodopensis plants. Results were analyzed according to the two types of water used in the experiment. Th EVOdrop device is equipped with an ultranano membrane and rotating jet nozzle to create a vortex water and saturation thanks to a second device EVObooster to obtain hydrogen-rich water. In the current study Hydrogen-rich water is referred to as Hydrogen EVOdrop Water (HEW). Research was conducted using the following methods-spectral methods non-equilibrium energy spectrum (NES) and differential non-equilibrium energy spectrum (DNES), mathematical models, and study of the distribution of water molecules in water clusters. In a licensed Eurotest Laboratory, the research of tap water before and after flowing through the EVOdrop device was proven. Studies have been carried out on the structuring of water molecule clusters after change of hydrogen bond energies. The restructuring comes with rearrangement of water molecules by the energy levels of hydrogen bonds. Local extrema can be observed in the spectrum with largest amount of water molecules. The structural changes were tested using the NES and DNES spectral methods. The conducted research proved that the application of EVOdrop device and EVObooster changes the parameters of water to benefit hydration and health.
Collapse
Affiliation(s)
- Ignat Ignatov
- Scientific Research Center of Medical Biophysics (SRCMB), 1111 Sofia, Bulgaria
| | | | - Nikolai Neshev
- Faculty of Physics, Sofia University “St. Kliment Ohridski”, 1000 Sofia, Bulgaria
| | - Yoana Kiselova-Kaneva
- Department of Biochemistry, Molecular Medicine and Nutrigenomics, Medical University-Varna, 9002 Varna, Bulgaria
| | - Teodora P. Popova
- Faculty of Veterinary Medicine, University of Forestry, 10 Kl. Ohridski Blvd., 1756 Sofia, Bulgaria
| | - Ralitsa Bankova
- Department of Internal Noncommunicable Diseases, Pathology and Pharmacology, Faculty of Veterinary Medicine, University of Forestry, 10 Kl. Ohridski Blvd., 1756 Sofia, Bulgaria
| | - Nedyalka Valcheva
- Faculty of Agriculture, Department Biochemistry, Microbiology, Physics, Trakia University, 6000 Stara Zagora, Bulgaria
| | | | - Mariana Angelcheva
- Department of Kinesitherapy and Rehabilitation, National Sports Academy “B. Levski”, 1700 Sofia, Bulgaria
| | | | | |
Collapse
|
7
|
Lee IS, Min SK. Generalized Formulation of the Density Functional Tight Binding-Based Restricted Ensemble Kohn-Sham Method with Onsite Correction to Long-Range Correction. J Chem Theory Comput 2022; 18:3391-3409. [PMID: 35549266 DOI: 10.1021/acs.jctc.2c00037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a generalized formulation for the combination of the density functional tight binding (DFTB) approach and the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method by considering onsite correction (OC) as well as the long-range corrected (LC) functional. The OC contribution provides more accurate energies and analytic gradients for individual microstates, while the multireference character of the SSR provides the correct description for conical intersections. We benchmark the LC-OC-DFTB/SSR method against various DFTB calculation methods for excitation energies and conical intersection structures with π/π* or n/π* characters. Furthermore, we perform excited-state molecular dynamics simulations with a molecular rotary motor with variations of LC-OC-DFTB/SSR approaches. We show that the OC contribution to the LC functional is crucial to obtain the correct geometry of conical intersections.
Collapse
Affiliation(s)
- In Seong Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| |
Collapse
|
8
|
Li C, Voth GA. Using Machine Learning to Greatly Accelerate Path Integral Ab Initio Molecular Dynamics. J Chem Theory Comput 2022; 18:599-604. [PMID: 34982562 PMCID: PMC8864787 DOI: 10.1021/acs.jctc.1c01085] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Ab initio molecular dynamics (AIMD) has become one of the most popular and robust approaches for modeling complicated chemical, liquid, and material systems. However, the formidable computational cost often limits its widespread application in simulations of the largest-scale systems. The situation becomes even more severe in cases where the hydrogen nuclei may be better described as quantized particles using a path integral representation. Here, we present a computational approach that combines machine learning with recent advances in path integral contraction schemes, and we achieve a 2 orders of magnitude acceleration over direct path integral AIMD simulation while at the same time maintaining its accuracy.
Collapse
|
9
|
Kitheka M, Redington M, Zhang J, Yao Y, Goyal P. BENCHMARKS OF THE DENSITY FUNCTIONAL TIGHT-BINDING METHOD FOR REDOX, PROTONATION AND ELECTRONIC PROPERTIES OF QUINONES. Phys Chem Chem Phys 2022; 24:6742-6756. [DOI: 10.1039/d1cp05333g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Organic materials with controllable molecular design and sustainable resources are promising electrode materials. Crystalline quinones have been investigated in a variety of rechargeable battery chemistries due to their ubiquitous nature,...
Collapse
|
10
|
Priyadarsini A, Mallik BS. Comparative first principles-based molecular dynamics study of catalytic mechanism and reaction energetics of water oxidation reaction on 2D-surface. J Comput Chem 2021; 42:1138-1149. [PMID: 33851446 DOI: 10.1002/jcc.26528] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Revised: 03/22/2021] [Accepted: 03/26/2021] [Indexed: 01/02/2023]
Abstract
The study of the water-splitting process, which can proceed in 2e- as well as 4e- pathway, reveals that the process is entirely an uphill process, and the third step, that is, the oxooxo bond formation is the rate-determining step. The kinetic barrier of the oxygen evolution reaction (OER) on the 2D material catalysts in the presence of explicit solvents is scarcely studied. Here, we investigate the dynamics of the OER on the undoped graphene and the activation energy barrier of each step using first principles molecular dynamics simulations. Here we provide a detailed analysis of the kinetics of all the 4e- transfer steps of OER on the graphene surface. We also compare the accuracy of one of the density functional theory (DFT) functionals and density functional based tight binding (DFTB) method in explaining the OER steps. The comparative study reveals that DFTB can be used for performing metadynamics simulations quipped with much less computational cost than DFT functionals. By both Perdew-Burke-Ernzerhof and DFTB methods, the third step is revealed to be the rate-determining step with an energy barrier of 21.19 ± 0.51 and 20.23 ± 0.20 kcal mol-1 , respectively. DFTB gives an impression of being successful in predicting the energy barriers of OER in 4e- transfer pathway and comparable to the DFT method, and we would like to extend the use of DFTB for further studies with a sizable and complex system.
Collapse
Affiliation(s)
- Adyasa Priyadarsini
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy, Telangana, India
| | - Bhabani S Mallik
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy, Telangana, India
| |
Collapse
|
11
|
Watanabe HC, Yamada M, Suzuki Y. Proton transfer in bulk water using the full adaptive QM/MM method: integration of solute- and solvent-adaptive approaches. Phys Chem Chem Phys 2021; 23:8344-8360. [PMID: 33875999 DOI: 10.1039/d1cp00116g] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The quantum mechanical/molecular mechanical (QM/MM) method is a hybrid molecular simulation technique that increases the accessibility of local electronic structures of large systems. The technique combines the benefit of accuracy found in the QM method and that of cost efficiency found in the MM method. However, it is difficult to directly apply the QM/MM method to the dynamics of solution systems, particularly for proton transfer. As explained in the Grotthuss mechanism, proton transfer is a structural interconversion between hydronium ions and solvent water molecules. Hence, when the QM/MM method is applied, an adaptive treatment, namely on-the-fly revisions on molecular definitions, is required for both the solute and solvent. Although several solvent-adaptive methods have been proposed, a full adaptive framework, which is an approach that also considers adaptation for solutes, remains untapped. In this paper, we propose a new numerical expression for the coordinates of the excess proton and its control algorithm. Furthermore, we confirm that this method can stably and accurately simulate proton transfer dynamics in bulk water.
Collapse
Affiliation(s)
- Hiroshi C Watanabe
- Quantum Computing Center, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan.
| | | | | |
Collapse
|
12
|
Cha J. Morphological effect of side chain on H 3O + transfer inside polymer electrolyte membranes across polymeric chain via molecular dynamics simulation. Sci Rep 2020; 10:22014. [PMID: 33328487 PMCID: PMC7745029 DOI: 10.1038/s41598-020-77971-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Accepted: 11/12/2020] [Indexed: 11/25/2022] Open
Abstract
Performance and durability of polymer electrolyte membrane are critical to fuel cell quality. As fuel cell vehicles become increasingly popular, membrane fundamentals must be understood in detail. Here, this study used molecular dynamic simulations to explore the morphological effects of perfluorosulfonic acid (PFSA)-based membranes on ionic conductivity. In particular, I developed an intuitive quantitative approach focusing principally on hydronium adsorbing to, and desorbing from, negatively charged sulfonate groups, while conventional ionic conductivity calculations featured the use of mean square displacements that included natural atomic vibrations. The results revealed that shorter side-chains caused more hydroniums to enter the conductive state, associated with higher ion conductivity. In addition, the hydronium path tracking showed that shorter side-chains allowed hydroniums to move among host groups, facilitating chain adsorption, in agreement with a mechanism suggested in earlier studies.
Collapse
Affiliation(s)
- JinHyeok Cha
- Institute of Fundamentals and Advanced Technology, Hyundai Motor Company, 37 Cheoldobangmulgwan-ro, Uiwang-si, Gyeonggi-do, 16082, Republic of Korea.
| |
Collapse
|
13
|
Jenness GR, Bresnahan CG, Shukla MK. Adventures in DFTB: Toward an Automatic Parameterization Scheme. J Chem Theory Comput 2020; 16:6894-6903. [PMID: 33119287 DOI: 10.1021/acs.jctc.0c00842] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
As we push forward on understanding the fate of chemicals in the environment, we need a method that will allow for the simulation of the inherent heterogeneity. Density functional tight binding (DFTB) is a methodology that allows for a detailed electronic description and would be ideal for this problem. While many parameters can be derived directly from DFT, empirical parameters still exist in the confinement and repulsion potentials. In this manuscript, we examine these potentials and present solutions that will minimize the degree of empiricism. Our results show that it is possible to construct confinement potentials from examining the atomic radial wavefunctions. Moreover, we found that the heterogeneous repulsion potentials can be derived from using only homogeneous repulsion curves.
Collapse
Affiliation(s)
- Glen R Jenness
- Environmental Laboratory, US Army Engineer Research and Development Center, 3909 Halls Ferry Road, Vicksburg, Mississippi 39180, United States
| | - Caitlin G Bresnahan
- Environmental Laboratory, US Army Engineer Research and Development Center, 3909 Halls Ferry Road, Vicksburg, Mississippi 39180, United States
| | - Manoj K Shukla
- Environmental Laboratory, US Army Engineer Research and Development Center, 3909 Halls Ferry Road, Vicksburg, Mississippi 39180, United States
| |
Collapse
|
14
|
Ono J, Imai M, Nishimura Y, Nakai H. Hydroxide Ion Carrier for Proton Pumps in Bacteriorhodopsin: Primary Proton Transfer. J Phys Chem B 2020; 124:8524-8539. [DOI: 10.1021/acs.jpcb.0c05507] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Junichi Ono
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
| | - Minori Imai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| |
Collapse
|
15
|
Cuny J, Cerda Calatayud J, Ansari N, Hassanali AA, Rapacioli M, Simon A. Simulation of Liquids with the Tight-Binding Density-Functional Approach and Improved Atomic Charges. J Phys Chem B 2020; 124:7421-7432. [PMID: 32696649 DOI: 10.1021/acs.jpcb.0c04167] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Theoretical description of liquids, especially liquid water, is an ongoing subject with important implications in various domains such as homogeneous catalysis; solvation of molecular, ionic, and biomolecular species; and reactivity. Various formalisms exist to describe liquids, each one displaying its own balance between accuracy and computational cost that defines its range of applications. The present article revisits the ability of the density-functional-based tight-binding (SCC-DFTB) approach to model liquids by focusing on liquid water and liquid benzene under ambient conditions. To do so, we benchmark a recent correction for the SCC-DFTB atomic charges that allows for a drastic improvement of the pair radial distribution functions of liquid water as compared to both experimental data and density-functional theory results performed in the generalized-gradient approximation. We also report the coupling of the deMonNano and i-PI codes to perform path-integral molecular dynamics. This allows us to rationalize the impact of nuclear quantum effects on the SCC-DFTB description of liquid water. This study evidences the rather good ability of SCC-DFTB to describe liquid water and liquid benzene. As the first example of application, we also present results for a benzene molecule solvated in water with the perspectives of further studies devoted to solvent/water interfaces.
Collapse
Affiliation(s)
- Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques (LCPQ), Université de Toulouse III [UPS] and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Jesus Cerda Calatayud
- Laboratoire de Chimie et Physique Quantiques (LCPQ), Université de Toulouse III [UPS] and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Narjes Ansari
- The Abdus Salam International Center for Theoretical Physics, Condensed Matter and Statistical Physics Section, Strada Costiera 11, 34151 Trieste, Italy.,Department of Chemistry and Applied Biosciences, ETH Zurich, 3 c/o USI Campus, Via Giuseppe Buffi13, 6900 Lugano, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland
| | - Ali A Hassanali
- The Abdus Salam International Center for Theoretical Physics, Condensed Matter and Statistical Physics Section, Strada Costiera 11, 34151 Trieste, Italy
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques (LCPQ), Université de Toulouse III [UPS] and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Aude Simon
- Laboratoire de Chimie et Physique Quantiques (LCPQ), Université de Toulouse III [UPS] and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| |
Collapse
|
16
|
Stöhr M, Medrano Sandonas L, Tkatchenko A. Accurate Many-Body Repulsive Potentials for Density-Functional Tight Binding from Deep Tensor Neural Networks. J Phys Chem Lett 2020; 11:6835-6843. [PMID: 32787209 DOI: 10.1021/acs.jpclett.0c01307] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
We combine density-functional tight binding (DFTB) with deep tensor neural networks (DTNN) to maximize the strengths of both approaches in predicting structural, energetic, and vibrational molecular properties. The DTNN is used to construct a nonlinear model for the localized many-body interatomic repulsive energy, which so far has been treated in an atom-pairwise manner in DFTB. Substantially improving upon standard DFTB and DTNN, the resulting DFTB-NNrep model yields accurate predictions of atomization and isomerization energies, equilibrium geometries, vibrational frequencies, and dihedral rotation profiles for a large variety of organic molecules compared to the hybrid DFT-PBE0 functional. Our results highlight the potential of combining semiempirical electronic-structure methods with physically motivated machine learning approaches for predicting localized many-body interactions. We conclude by discussing future advancements of the DFTB-NNrep approach that could enable chemically accurate electronic-structure calculations for systems with tens of thousands of atoms.
Collapse
Affiliation(s)
- Martin Stöhr
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Leonardo Medrano Sandonas
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| |
Collapse
|
17
|
Spiegelman F, Tarrat N, Cuny J, Dontot L, Posenitskiy E, Martí C, Simon A, Rapacioli M. Density-functional tight-binding: basic concepts and applications to molecules and clusters. ADVANCES IN PHYSICS: X 2020; 5:1710252. [PMID: 33154977 PMCID: PMC7116320 DOI: 10.1080/23746149.2019.1710252] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Accepted: 12/19/2019] [Indexed: 06/10/2023] Open
Abstract
The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.
Collapse
Affiliation(s)
- Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Nathalie Tarrat
- CEMES, Université de Toulouse (UPS), CNRS, UPR8011, Toulouse, Toulouse, France
| | - Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Leo Dontot
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Evgeny Posenitskiy
- Laboratoire Collisions Agrégats et Réactivité LCAR/IRSAMC, UMR5589, Université de Toulouse (UPS) and CNRS, Toulouse, France
| | - Carles Martí
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
- Laboratoire de Chimie, UMR5182, Ecole Normale Supérieure de Lyon, Université de Lyon and CNRS, Lyon, France
| | - Aude Simon
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| |
Collapse
|
18
|
Lourenço MP, dos Santos EC, Pettersson LGM, Duarte HA. Accurate SCC-DFTB Parametrization for Bulk Water. J Chem Theory Comput 2020; 16:1768-1778. [DOI: 10.1021/acs.jctc.9b00816] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Maicon Pierre Lourenço
- Departamento de Quı́mica e Fı́sica−Centro de Ciências Exatas, Naturais e da Saúde (CCENS), Universidade Federal do Espı́rito Santo, 29500-000, Alegre, Espı́rito Santo Brasil
| | - Egon Campos dos Santos
- Departamento de Quı́mica, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, Minas Gerais Brasil
- Department of Physics, AlbaNova University Center, Stockholm University, S-10691 Stockholm, Sweden
| | - Lars G. M. Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, S-10691 Stockholm, Sweden
| | - Hélio Anderson Duarte
- Departamento de Quı́mica, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, Minas Gerais Brasil
| |
Collapse
|
19
|
Sakti AW, Nishimura Y, Nakai H. Recent advances in quantum‐mechanical molecular dynamics simulations of proton transfer mechanism in various water‐based environments. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1419] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Aditya W. Sakti
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
| | - Hiromi Nakai
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering Waseda University Tokyo Japan
| |
Collapse
|
20
|
Villa AM, Doglia SM, De Gioia L, Bertini L, Natalello A. Anomalous Intrinsic Fluorescence of HCl and NaOH Aqueous Solutions. J Phys Chem Lett 2019; 10:7230-7236. [PMID: 31689111 DOI: 10.1021/acs.jpclett.9b02163] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The unique properties of liquid water mainly arise from its hydrogen bond network. The geometry and dynamics of this network play a key role in shaping the characteristics of soft matter, from simple solutions to biosystems. Here we report an anomalous intrinsic fluorescence of HCl and NaOH aqueous solutions at room temperature that shows important differences in the excitation and emission bands between the two solutes. From ab initio time-dependent density functional theory modeling we propose that fluorescence emission could originate from hydrated ion species contained in transient cavities of the bulk solvent. These cavities, which are characterized by a stiff surface, could provide an environment that, upon trapping the excited state, suppresses the fast nonradiative decay and allows the slower radiative channel to become a possible decay pathway.
Collapse
Affiliation(s)
- Anna Maria Villa
- Department of Biotechnology and Biosciences , University of Milano-Bicocca , Piazza della Scienza 2 , 20126 Milan , Italy
| | - Silvia Maria Doglia
- Department of Biotechnology and Biosciences , University of Milano-Bicocca , Piazza della Scienza 2 , 20126 Milan , Italy
| | - Luca De Gioia
- Department of Biotechnology and Biosciences , University of Milano-Bicocca , Piazza della Scienza 2 , 20126 Milan , Italy
| | - Luca Bertini
- Department of Biotechnology and Biosciences , University of Milano-Bicocca , Piazza della Scienza 2 , 20126 Milan , Italy
| | - Antonino Natalello
- Department of Biotechnology and Biosciences , University of Milano-Bicocca , Piazza della Scienza 2 , 20126 Milan , Italy
| |
Collapse
|
21
|
Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk. J Phys Chem B 2019; 123:9846-9861. [PMID: 31647873 DOI: 10.1021/acs.jpcb.9b03258] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H2O)100H+, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.
Collapse
|
22
|
Watanabe HC, Cui Q. Quantitative Analysis of QM/MM Boundary Artifacts and Correction in Adaptive QM/MM Simulations. J Chem Theory Comput 2019; 15:3917-3928. [DOI: 10.1021/acs.jctc.9b00180] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Hiroshi C. Watanabe
- Quantum Computing Center, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
- Japan Science and Technology Agency, PRESTO, 4-18, Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
23
|
Chou CP, Sakti AW, Nishimura Y, Nakai H. Development of Divide-and-Conquer Density-Functional Tight-Binding Method for Theoretical Research on Li-Ion Battery. CHEM REC 2019; 19:746-757. [PMID: 30462370 DOI: 10.1002/tcr.201800141] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2018] [Revised: 10/25/2018] [Accepted: 10/26/2018] [Indexed: 01/24/2023]
Abstract
The density-functional tight-binding (DFTB) method is one of the useful quantum chemical methods, which provides a good balance between accuracy and computational efficiency. In this account, we reviewed the basis of the DFTB method, the linear-scaling divide-and-conquer (DC) technique, as well as the parameterization process. We also provide some refinement, modifications, and extension of the existing parameters that can be applicable for lithium-ion battery systems. The diffusion constants of common electrolyte molecules and LiTFSA salt in solution have been estimated using DC-DFTB molecular dynamics simulation with our new parameters. The resulting diffusion constants have good agreement to the experimental diffusion constants.
Collapse
Affiliation(s)
- Chien-Pin Chou
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, Tokyo, 169-8555, Japan
| | - Aditya Wibawa Sakti
- Element Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyotodaigaku-Katsura, Kyoto, 615-8520, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, Tokyo, 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, Tokyo, 169-8555, Japan.,Element Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyotodaigaku-Katsura, Kyoto, 615-8520, Japan.,Department of Chemistry and Biochemistry, School of Advanced Science and Enigineering, Waseda University, Tokyo, 169-8555, Japan
| |
Collapse
|
24
|
Simon A, Rapacioli M, Michoulier E, Zheng L, Korchagina K, Cuny J. Contribution of the density-functional-based tight-binding scheme to the description of water clusters: methods, applications and extension to bulk systems. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1554903] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- A. Simon
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - M. Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - E. Michoulier
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
- Laboratoire Collisions Agrégats et Réactivité LCAR/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - L. Zheng
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - K. Korchagina
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| | - J. Cuny
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse and CNRS, Toulouse, France
| |
Collapse
|
25
|
Ullah N, Chen S, Zhang R. Excited state dynamics study of the self-trapped exciton formation in silicon nanosheets. Phys Chem Chem Phys 2018; 20:29299-29305. [DOI: 10.1039/c8cp04806a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
After excitation to S1 (1), the exciton takes ∼450–850 femtoseconds to relax into the self-trapped (ST) state (2) with the occurrence of strong localization and a large Stokes shift, due to the significant stretching of the Si–Si bonds.
Collapse
Affiliation(s)
- Naeem Ullah
- Department of Physics, City University of Hong Kong
- Hong Kong SAR
- China
| | - Shunwei Chen
- Department of Physics, City University of Hong Kong
- Hong Kong SAR
- China
| | - Ruiqin Zhang
- Department of Physics, City University of Hong Kong
- Hong Kong SAR
- China
- Beijing Computational Science Research Center
- Beijing 100193
| |
Collapse
|
26
|
Sakti AW, Nishimura Y, Chou CP, Nakai H. Density-Functional Tight-Binding Molecular Dynamics Simulations of Excess Proton Diffusion in Ice Ih, Ice Ic, Ice III, and Melted Ice VI Phases. J Phys Chem A 2017; 122:33-40. [DOI: 10.1021/acs.jpca.7b10664] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
| | | | | | - Hiromi Nakai
- CREST, Japan Science and Technology Agency, Tokyo 102-0075, Japan
- ESICB, Kyoto University, Kyoto 615-8520, Japan
| |
Collapse
|
27
|
Molecular dynamics simulations of zinc oxide solubility: From bulk down to nanoparticles. Food Chem Toxicol 2017; 112:518-525. [PMID: 28736190 DOI: 10.1016/j.fct.2017.07.038] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2016] [Revised: 06/16/2017] [Accepted: 07/19/2017] [Indexed: 11/23/2022]
Abstract
The solubility of metal oxides is one of the key descriptors for the evaluation of their potential toxic effects, both in the bulk form and in nanoparticulated aggregates. Current work presents a new methodology for the in silico assessment of the solubility of metal oxides, which is demonstrated using a well-studied system, ZnO. The calculation of the solubility is based on statistical thermodynamics tools combined with Density Functional Tight Binding theory for the evaluation of the free energy exchange during the dissolution process. Models of small ZnO clusters are used for describing the final dissolved material, since the complete ionic dissolution of ZnO is hindered by the formation of O2- anions in solution, which are highly unstable. Results show very good agreement between the computed solubility values and experimental data for ZnO bulk, up to 0.5 mg L-1 and equivalents of 50 μg L-1 for the free Zn2+ cation in solution. However, the reference model for solid nanoparticles formed by free space nanoparticles can only give a limited quantitative solubility evaluation for ZnO nanoparticles.
Collapse
|
28
|
Sakti AW, 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. [DOI: 10.1021/acs.jpcb.6b10659] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | | | - Hiromi Nakai
- Core
Research for Evolutional Science and Technology (CREST), Japan Science and Technology Agency (JST), Chiyoda-ku, Tokyo 102-0075, Japan
- Elements
Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyoto 615-8520, Japan
| |
Collapse
|
29
|
Mazzuca JW, Schultz CP. Quantum Mechanical Enhancement of Rate Constants and Kinetic Isotope Effects for Water-Mediated Proton Transfer in a Model Biological System. J Phys Chem A 2017; 121:819-826. [DOI: 10.1021/acs.jpca.6b10337] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- James W. Mazzuca
- Chemistry Department, Alma College, Alma, Michigan 48801, United States
| | - Chase P. Schultz
- Chemistry Department, Alma College, Alma, Michigan 48801, United States
| |
Collapse
|
30
|
Watanabe HC, Kubillus M, Kubař T, Stach R, Mizaikoff B, Ishikita H. Cation solvation with quantum chemical effects modeled by a size-consistent multi-partitioning quantum mechanics/molecular mechanics method. Phys Chem Chem Phys 2017; 19:17985-17997. [DOI: 10.1039/c7cp01708a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In the condensed phase, quantum chemical properties such as many-body effects and intermolecular charge fluctuations are critical determinants of the solvation structure and dynamics.
Collapse
Affiliation(s)
- Hiroshi C. Watanabe
- Research Center for Advanced Science and Technology
- The University of Tokyo
- Tokyo 153-8904
- Japan
- Department of Applied Chemistry
| | - Maximilian Kubillus
- Institute of Physical Chemistry and Center for Functional Nanostructures
- Karlsruhe Institute of Technology
- Karlsruhe 73131
- Germany
| | - Tomáš Kubař
- Institute of Physical Chemistry and Center for Functional Nanostructures
- Karlsruhe Institute of Technology
- Karlsruhe 73131
- Germany
| | - Robert Stach
- Institute of Analytical Bioanalytical Chemistry
- Ulm University
- Ulm 89081
- Germany
| | - Boris Mizaikoff
- Institute of Analytical Bioanalytical Chemistry
- Ulm University
- Ulm 89081
- Germany
| | - Hiroshi Ishikita
- Research Center for Advanced Science and Technology
- The University of Tokyo
- Tokyo 153-8904
- Japan
- Department of Applied Chemistry
| |
Collapse
|
31
|
Christensen A, Kubař T, Cui Q, Elstner M. Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications. Chem Rev 2016; 116:5301-37. [PMID: 27074247 PMCID: PMC4867870 DOI: 10.1021/acs.chemrev.5b00584] [Citation(s) in RCA: 246] [Impact Index Per Article: 30.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2015] [Indexed: 12/28/2022]
Abstract
Semiempirical (SE) methods can be derived from either Hartree-Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.
Collapse
Affiliation(s)
- Anders
S. Christensen
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Tomáš Kubař
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Marcus Elstner
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| |
Collapse
|
32
|
John C, Spura T, Habershon S, Kühne TD. Quantum ring-polymer contraction method: Including nuclear quantum effects at no additional computational cost in comparison to ab initio molecular dynamics. Phys Rev E 2016; 93:043305. [PMID: 27176426 DOI: 10.1103/physreve.93.043305] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2016] [Indexed: 06/05/2023]
Abstract
We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.
Collapse
Affiliation(s)
- Christopher John
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Thomas Spura
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Thomas D Kühne
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany and Paderborn Center for Parallel Computing and Institute for Lightweight Design, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| |
Collapse
|
33
|
Nakai H, Sakti AW, 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-21. [DOI: 10.1021/acs.jpcb.5b12439] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hiromi Nakai
- Department
of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, Tokyo 169-8555, Japan
- Research
Institute for Science and Engineering, Waseda University, Tokyo 169-8555, Japan
- CREST, Japan Science and Technology Agency, Tokyo 102-0075, Japan
- ESICB, Kyoto University, Kyoto 615-8520, Japan
| | - Aditya Wibawa Sakti
- Department
of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| | - Yoshifumi Nishimura
- Research
Institute for Science and Engineering, Waseda University, Tokyo 169-8555, Japan
- Department
of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki 444-8585, Japan
| |
Collapse
|
34
|
Watanabe HC, Banno M, Sakurai M. An adaptive quantum mechanics/molecular mechanics method for the infrared spectrum of water: incorporation of the quantum effect between solute and solvent. Phys Chem Chem Phys 2016; 18:7318-33. [DOI: 10.1039/c5cp07136d] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Quantum effects in solute–solvent interactions, such as the many-body effect and the dipole-induced dipole, are known to be critical factors influencing the infrared spectra of species in the liquid phase.
Collapse
Affiliation(s)
- Hiroshi C. Watanabe
- Research Center for Advanced Science and Technology
- The University of Tokyo
- Meguro-ku
- Japan
| | - Misa Banno
- Center for Biological Resources and Informatics
- Tokyo Institute of Technology
- Midori-ku
- Japan
| | - Minoru Sakurai
- Center for Biological Resources and Informatics
- Tokyo Institute of Technology
- Midori-ku
- Japan
| |
Collapse
|
35
|
Cui W, Cerqueira TFT, Botti S, Marques MAL, San-Miguel A. Nanostructured water and carbon dioxide inside collapsing carbon nanotubes at high pressure. Phys Chem Chem Phys 2016; 18:19926-32. [DOI: 10.1039/c6cp03263j] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
We present simulations of the collapse under hydrostatic pressure of carbon nanotubes containing either water or carbon dioxide.
Collapse
Affiliation(s)
- Wenwen Cui
- Institut Lumière Matière
- UMR5306 Université Lyon 1-CNRS
- Université de Lyon
- F-69622 Villeurbanne Cedex
- France
| | - Tiago F. T. Cerqueira
- Institut Lumière Matière
- UMR5306 Université Lyon 1-CNRS
- Université de Lyon
- F-69622 Villeurbanne Cedex
- France
| | - Silvana Botti
- Institut Lumière Matière
- UMR5306 Université Lyon 1-CNRS
- Université de Lyon
- F-69622 Villeurbanne Cedex
- France
| | - Miguel A. L. Marques
- Institut für Physik
- Martin-Luther-Universität Halle-Wittenberg
- D-06099 Halle
- Germany
- Institut Lumière Matière
| | - Alfonso San-Miguel
- Institut Lumière Matière
- UMR5306 Université Lyon 1-CNRS
- Université de Lyon
- F-69622 Villeurbanne Cedex
- France
| |
Collapse
|
36
|
Wang S, MacKay L, Lamoureux G. Development of Semiempirical Models for Proton Transfer Reactions in Water. J Chem Theory Comput 2015; 10:2881-90. [PMID: 26588263 DOI: 10.1021/ct500164h] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
This letter presents a method for the parametrization of semiempirical models for proton transfer reactions in water clusters. Two new models are developed: AM1-W, which is a reparameterization of the classic AM1 model, and AM1PG-W, which is a modified AM1-like model including a pairwise correction to the core repulsion function. Both models show good performance on hydrogen-bonding energies and on proton transfer energy profiles, which are of great importance for proton transfer reactions in large water clusters and in proteins. The parametrization method introduced is general and can be used to develop any other system-specific semiempirical models.
Collapse
Affiliation(s)
- Shihao Wang
- Department of Chemistry and Biochemistry and Centre for Research in Molecular Modeling (CERMM) and ‡Department of Physics, Concordia University , Montréal, Canada
| | - Laurent MacKay
- Department of Chemistry and Biochemistry and Centre for Research in Molecular Modeling (CERMM) and ‡Department of Physics, Concordia University , Montréal, Canada
| | - Guillaume Lamoureux
- Department of Chemistry and Biochemistry and Centre for Research in Molecular Modeling (CERMM) and ‡Department of Physics, Concordia University , Montréal, Canada
| |
Collapse
|
37
|
Chou CP, Nishimura Y, Fan CC, Mazur G, Irle S, Witek HA. Automatized Parameterization of DFTB Using Particle Swarm Optimization. J Chem Theory Comput 2015; 12:53-64. [DOI: 10.1021/acs.jctc.5b00673] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Chien-Pin Chou
- Department
of Applied Chemistry and Institute of Molecular Science, National Chiao Tung University, Hsinchu 30010, Taiwan
| | - Yoshifumi Nishimura
- Department
of Applied Chemistry and Institute of Molecular Science, National Chiao Tung University, Hsinchu 30010, Taiwan
| | - Chin-Chai Fan
- Department
of Applied Chemistry and Institute of Molecular Science, National Chiao Tung University, Hsinchu 30010, Taiwan
| | - Grzegorz Mazur
- K. Gumiński
Department of Theoretical Chemistry, Jagiellonian University, Ingardena
3, 30-060 Cracow, Poland
| | | | - Henryk A. Witek
- Department
of Applied Chemistry and Institute of Molecular Science, National Chiao Tung University, Hsinchu 30010, Taiwan
| |
Collapse
|
38
|
Time-dependent density-functional tight-binding method with the third-order expansion of electron density. J Chem Phys 2015; 143:094108. [DOI: 10.1063/1.4929926] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
39
|
Nishimoto Y, Fedorov DG, Irle S. Third-order density-functional tight-binding combined with the fragment molecular orbital method. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.07.022] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
40
|
Wolf S, Freier E, Gerwert K. A delocalized proton-binding site within a membrane protein. Biophys J 2015; 107:174-84. [PMID: 24988352 DOI: 10.1016/j.bpj.2014.05.019] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2013] [Revised: 05/13/2014] [Accepted: 05/14/2014] [Indexed: 11/26/2022] Open
Abstract
The role of protein-bound water molecules in protein function and catalysis is an emerging topic. Here, we studied the solvation of an excess proton by protein-bound water molecules and the contribution of the surrounding amino acid residues at the proton release site of the membrane protein bacteriorhodopsin. It hosts an excess proton within a protein-bound water cluster, which is hydrogen bonded to several surrounding amino acids. Indicative of delocalization is a broad continuum absorbance experimentally observed by time-resolved Fourier transform infrared spectroscopy. In combination with site-directed mutagenesis, the involvement of several amino acids (especially Glu-194 and Glu-204) in the delocalization was elaborated. Details regarding the contributions of the glutamates and water molecules to the delocalization mode in biomolecular simulations are controversial. We carried out quantum mechanics/molecular mechanics (QM/MM) self-consistent charge density functional tight-binding simulations for all amino acids that have been experimentally shown to be involved in solvation of the excess proton, and systematically investigated the influence of the quantum box size. We compared calculated theoretical infrared spectra with experimental ones as a measure for the correct description of excess proton delocalization. A continuum absorbance can only be observed for small quantum boxes containing few amino acids and/or water molecules. Larger quantum boxes, including all experimentally shown involved amino acids, resulted in narrow absorbance bands, indicating protonation of a single binding site in contradiction to experimental results. We conclude that small quantum boxes seem to reproduce representative extreme cases of proton delocalization modes: proton delocalization only on water molecules or only between Glu-194 and Glu-204. Extending the experimental spectral region to lower wave numbers, a water-delocalized proton reproduces the observed continuum absorbance better than a glutamate-shared delocalized proton. However, a full agreement between QM simulations and experimental results on the delocalized excess proton will require a larger quantum box as well as more sophisticated QM/MM methods.
Collapse
Affiliation(s)
- Steffen Wolf
- Department of Biophysics, Chinese Academy of Sciences-Max-Planck-Gesellschaft Partner Institute for Computational Biology, Shanghai Institutes for Biological Sciences, Shanghai, P.R. China
| | - Erik Freier
- Department of Biophysics, University of Bochum, Bochum, Germany
| | - Klaus Gerwert
- Department of Biophysics, Chinese Academy of Sciences-Max-Planck-Gesellschaft Partner Institute for Computational Biology, Shanghai Institutes for Biological Sciences, Shanghai, P.R. China; Department of Biophysics, University of Bochum, Bochum, Germany.
| |
Collapse
|
41
|
Pezeshki S, Lin H. Adaptive-Partitioning QM/MM for Molecular Dynamics Simulations: 4. Proton Hopping in Bulk Water. J Chem Theory Comput 2015; 11:2398-411. [DOI: 10.1021/ct501019y] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Soroosh Pezeshki
- Chemistry Department, University of Colorado Denver, Denver, Colorado 80217-3364, United States
| | - Hai Lin
- Chemistry Department, University of Colorado Denver, Denver, Colorado 80217-3364, United States
| |
Collapse
|
42
|
Tse YLS, Knight C, Voth GA. An analysis of hydrated proton diffusion in ab initio molecular dynamics. J Chem Phys 2015; 142:014104. [DOI: 10.1063/1.4905077] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ying-Lung Steve Tse
- Department of Chemistry, James Franck Institute, and Computation Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Chris Knight
- Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, and Computation Institute, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
43
|
Goyal P, Yang S, Cui Q. Microscopic basis for kinetic gating in Cytochrome c oxidase: insights from QM/MM analysis. Chem Sci 2015; 6:826-841. [PMID: 25678950 PMCID: PMC4321873 DOI: 10.1039/c4sc01674b] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps.
Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps. For the proton pump cytochrome c oxidase, while the regulation of thermodynamic driving force for key proton transfers has been discussed in great detail, the microscopic basis for the control of proton transfer kinetics has been poorly understood. Here we carry out extensive QM/MM free energy simulations to probe the kinetics of relevant proton transfer steps and analyze the effects of local structure and hydration level. We show that protonation of the proton loading site (PLS, taken to be a propionate of heme a3) requires a concerted process in which a key glutamic acid (Glu286H) delivers the proton to the PLS while being reprotonated by an excess proton coming from the D-channel. The concerted nature of the mechanism is a crucial feature that enables the loading of the PLS before the cavity containing Glu286 is better hydrated to lower its pKa to experimentally measured range; the charged rather than dipolar nature of the process also ensures a tight coupling with heme a reduction, as emphasized by Siegbahn and Blomberg. In addition, we find that rotational flexibility of the PLS allows its protonation before that of the binuclear center (the site where oxygen gets reduced to water). Together with our recent study (P. Goyal, et al., Proc. Natl. Acad. Sci. U. S. A., 2013, 110, 18886–18891) that focused on the modulation of Glu286 pKa, the current work suggests a mechanism that builds in a natural sequence for the protonation of the PLS prior to that of the binuclear center. This provides microscopic support to the kinetic constraints revealed by kinetic network analysis as essential elements that ensure an efficient vectorial proton transport in cytochrome c oxidase.
Collapse
Affiliation(s)
- Puja Goyal
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Avenue, Madison, WI 53706
| | - Shuo Yang
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Avenue, Madison, WI 53706
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Avenue, Madison, WI 53706
| |
Collapse
|
44
|
Song R, Feng W, Jimenez-Cruz CA, Wang B, Jiang W, Wang Z, Zhou R. Water film inside graphene nanosheets: electron transfer reversal between water and graphene via tight nano-confinement. RSC Adv 2015. [DOI: 10.1039/c4ra13736a] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Electron transfer reversal between water and grapheneviatight nano-confinement.
Collapse
Affiliation(s)
- Ruixia Song
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy (Jilin University)
| | - Wei Feng
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy (Jilin University)
| | | | - Bo Wang
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy (Jilin University)
| | - Wanrun Jiang
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy (Jilin University)
| | - Zhigang Wang
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy (Jilin University)
| | - Ruhong Zhou
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun
- China
- Computational Biology Center
| |
Collapse
|
45
|
Mazzuca JW, Garashchuk S, Jakowski J. The effect of local substrate motion on quantum hydrogen transfer in soybean lipoxygenase-1 modeled with QTES-DFTB dynamics. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.08.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
46
|
Goyal P, Qian HJ, Irle S, Lu X, Roston D, Mori T, Elstner M, Cui Q. Molecular simulation of water and hydration effects in different environments: challenges and developments for DFTB based models. J Phys Chem B 2014; 118:11007-27. [PMID: 25166899 PMCID: PMC4174991 DOI: 10.1021/jp503372v] [Citation(s) in RCA: 77] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
We discuss the description of water and hydration effects that employs an approximate density functional theory, DFTB3, in either a full QM or QM/MM framework. The goal is to explore, with the current formulation of DFTB3, the performance of this method for treating water in different chemical environments, the magnitude and nature of changes required to improve its performance, and factors that dictate its applicability to reactions in the condensed phase in a QM/MM framework. A relatively minor change (on the scale of kBT) in the O-H repulsive potential is observed to substantially improve the structural properties of bulk water under ambient conditions; modest improvements are also seen in dynamic properties of bulk water. This simple change also improves the description of protonated water clusters, a solvated proton, and to a more limited degree, a solvated hydroxide. By comparing results from DFTB3 models that differ in the description of water, we confirm that proton transfer energetics are adequately described by the standard DFTB3/3OB model for meaningful mechanistic analyses. For QM/MM applications, a robust parametrization of QM-MM interactions requires an explicit consideration of condensed phase properties, for which an efficient sampling technique was developed recently and is reviewed here. The discussions help make clear the value and limitations of DFTB3 based simulations, as well as the developments needed to further improve the accuracy and transferability of the methodology.
Collapse
Affiliation(s)
- Puja Goyal
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | | | | | | | | | | | | | | |
Collapse
|
47
|
Nishimura Y, Lee YP, Irle S, Witek HA. Critical interpretation of CH– and OH– stretching regions for infrared spectra of methanol clusters (CH3OH)n (n = 2–5) using self-consistent-charge density functional tight-binding molecular dynamics simulations. J Chem Phys 2014; 141:094303. [DOI: 10.1063/1.4893952] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yoshifumi Nishimura
- Department of Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
- Department of Applied Chemistry and Institute of Molecular Science, National Chiao Tung University, 1001 Ta-Hsueh Road, Hsinchu 30010, Taiwan
| | - Yuan-Pern Lee
- Department of Applied Chemistry and Institute of Molecular Science, National Chiao Tung University, 1001 Ta-Hsueh Road, Hsinchu 30010, Taiwan
- Institute of Atomic and Molecular Sciences, Academia Sinica, P.O. Box 23-166, Taipei 10617, Taiwan
| | - Stephan Irle
- Department of Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | - Henryk A. Witek
- Department of Applied Chemistry and Institute of Molecular Science, National Chiao Tung University, 1001 Ta-Hsueh Road, Hsinchu 30010, Taiwan
| |
Collapse
|
48
|
Cui Q, Elstner M. Density functional tight binding: values of semi-empirical methods in an ab initio era. Phys Chem Chem Phys 2014; 16:14368-77. [PMID: 24850383 PMCID: PMC4836871 DOI: 10.1039/c4cp00908h] [Citation(s) in RCA: 102] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Semi-empirical (SE) methods are derived from Hartree-Fock (HF) or Density Functional Theory (DFT) by neglect and approximation of electronic integrals. Thereby, parameters are introduced which have to be determined from reference calculations and/or by fitting to available experimental data. This leads to computational methods that are about 2-3 orders of magnitude faster than the standard HF/DFT methods using medium sized basis sets while being about 3 orders of magnitude slower than empirical force field methods (Molecular Mechanics: MM). Therefore, SE methods are most appropriate for a specific range of applications. These include the study of systems that contain a large number of atoms and therefore being too large for ab initio or DFT methods and also problems where dynamic or entropic effects are particularly important. In the latter case, the errors made by considering a very limited number of molecular structures or neglecting entropic contributions can be much larger than the accuracy lost due to the use of SE methods. Another area where SE methods are attractive concerns the analysis of systems for which reliable MM models are not readily available. Therefore, even in an era when rapid progress is being made in ab initio methods, there is considerable interest in further developing SE methods. We illustrate this point by focusing on the discussion of recent development and application of the Density Functional Tight Binding method.
Collapse
|
49
|
Jahangiri S, Cai L, Peslherbe GH. Performance of density-functional tight-binding models in describing hydrogen-bonded anionic-water clusters. J Comput Chem 2014; 35:1707-15. [DOI: 10.1002/jcc.23677] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2014] [Revised: 06/19/2014] [Accepted: 06/20/2014] [Indexed: 01/04/2023]
Affiliation(s)
- Soran Jahangiri
- Department of Chemistry and Biochemistry; Centre for Research in Molecular Modeling (CERMM), Concordia University; 7141 Sherbrooke Street West Montréal Québec Canada H4B 1R6
| | - Lemin Cai
- Department of Chemistry and Biochemistry; Centre for Research in Molecular Modeling (CERMM), Concordia University; 7141 Sherbrooke Street West Montréal Québec Canada H4B 1R6
| | - Gilles H. Peslherbe
- Department of Chemistry and Biochemistry; Centre for Research in Molecular Modeling (CERMM), Concordia University; 7141 Sherbrooke Street West Montréal Québec Canada H4B 1R6
| |
Collapse
|
50
|
Giese T, Chen H, Huang M, York DM. Parametrization of an Orbital-Based Linear-Scaling Quantum Force Field for Noncovalent Interactions. J Chem Theory Comput 2014; 10:1086-1098. [PMID: 24803856 PMCID: PMC3985928 DOI: 10.1021/ct401035t] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2013] [Indexed: 01/22/2023]
Abstract
We parametrize a linear-scaling quantum mechanical force field called mDC for the accurate reproduction of nonbonded interactions. We provide a new benchmark database of accurate ab initio interactions between sulfur-containing molecules. A variety of nonbond databases are used to compare the new mDC method with other semiempirical, molecular mechanical, ab initio, and combined semiempirical quantum mechanical/molecular mechanical methods. It is shown that the molecular mechanical force field significantly and consistently reproduces the benchmark results with greater accuracy than the semiempirical models and our mDC model produces errors twice as small as the molecular mechanical force field. The comparisons between the methods are extended to the docking of drug candidates to the Cyclin-Dependent Kinase 2 protein receptor. We correlate the protein-ligand binding energies to their experimental inhibition constants and find that the mDC produces the best correlation. Condensed phase simulation of mDC water is performed and shown to produce O-O radial distribution functions similar to TIP4P-EW.
Collapse
Affiliation(s)
- Timothy
J. Giese
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| | - Haoyuan Chen
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| | - Ming Huang
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
- Scientific
Computation, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455−0431, United States
| | - Darrin M. York
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| |
Collapse
|