1
|
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
|
2
|
Lourenço MP, Hostaš J, Herrera LB, Calaminici P, Köster AM, Tchagang A, Salahub DR. GAMaterial-A genetic-algorithm software for material design and discovery. J Comput Chem 2023; 44:814-823. [PMID: 36444916 DOI: 10.1002/jcc.27043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 09/26/2022] [Accepted: 11/06/2022] [Indexed: 11/30/2022]
Abstract
Genetic algorithms (GAs) are stochastic global search methods inspired by biological evolution. They have been used extensively in chemistry and materials science coupled with theoretical methods, ranging from force-fields to high-throughput first-principles methods. The methodology allows an accurate and automated structural determination for molecules, atomic clusters, nanoparticles, and solid surfaces, fundamental to understanding chemical processes in catalysis and environmental sciences, for instance. In this work, we propose a new genetic algorithm software, GAMaterial, implemented in Python3.x, that performs global searches to elucidate the structures of atomic clusters, doped clusters or materials and atomic clusters on surfaces. For all these applications, it is possible to accelerate the GA search by using machine learning (ML), the ML@GA method, to build subsequent populations. Results for ML@GA applied for the dopant distributions in atomic clusters are presented. The GAMaterial software was applied for the automatic structural search for the Ti6 O12 cluster, doping Al in Si11 (4Al@Si11 ) and Na10 supported on graphene (Na10 @graphene), where DFTB calculations were used to sample the complex search surfaces with reasonably low computational cost. Finally, the global search by GA of the Mo8 C4 cluster was considered, where DFT calculations were made with the deMon2k code, which is interfaced with GAMaterial.
Collapse
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, Espírito Santo, Brazil
| | - Jiří Hostaš
- Department of Chemistry, Department of Physics and Astronomy, CMS Centre for Molecular Simulation, IQST Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, Calgary, Alberta, Canada
| | - Lizandra Barrios Herrera
- Department of Chemistry, Department of Physics and Astronomy, CMS Centre for Molecular Simulation, IQST Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, Calgary, Alberta, Canada
| | | | | | - Alain Tchagang
- Digital Technologies Research Centre, National Research Council of Canada, Ottawa, Ontario, Canada
| | - Dennis R Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS Centre for Molecular Simulation, IQST Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, Calgary, Alberta, Canada
| |
Collapse
|
3
|
Medvedev N, Voronkov R, Volkov AE. Metallic water: Transient state under ultrafast electronic excitation. J Chem Phys 2023; 158:074501. [PMID: 36813717 DOI: 10.1063/5.0139802] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
The modern means of controlled irradiation by femtosecond lasers or swift heavy ion beams can transiently produce such energy densities in samples that reach collective electronic excitation levels of the warm dense matter state, where the potential energy of interaction of the particles is comparable to their kinetic energies (temperatures of a few eV). Such massive electronic excitation severely alters the interatomic potentials, producing unusual nonequilibrium states of matter and different chemistry. We employ density functional theory and tight binding molecular dynamics formalisms to study the response of bulk water to ultrafast excitation of its electrons. After a certain threshold electronic temperature, the water becomes electronically conducting via the collapse of its bandgap. At high doses, it is accompanied by nonthermal acceleration of ions to a temperature of a few thousand Kelvins within sub-100 fs timescales. We identify the interplay of this nonthermal mechanism with the electron-ion coupling, enhancing the electron-to-ions energy transfer. Various chemically active fragments are formed from the disintegrating water molecules, depending on the deposited dose.
Collapse
Affiliation(s)
- Nikita Medvedev
- Department of Radiation and Chemical Physics, Institute of Physics, Czech Academy of Sciences, Na Slovance 1999/2, 182 21 Prague 8, Czech Republic
| | - Roman Voronkov
- P.N. Lebedev Physical Institute of the Russian Academy of Sciences, Leninskij pr., 53, 119991 Moscow, Russia
| | - Alexander E Volkov
- P.N. Lebedev Physical Institute of the Russian Academy of Sciences, Leninskij pr., 53, 119991 Moscow, Russia
| |
Collapse
|
4
|
Souza TG, Olusegun SJ, Galvao BR, Da Silva JL, Mohallem ND, Ciminelli VS. Mechanism of amoxicillin adsorption by ferrihydrites: experimental and computational approaches. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
5
|
Wróbel P, Kubisiak P, Eilmes A. Hydrogen Bonding and Infrared Spectra of Ethyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide/Water Mixtures: A View from Molecular Dynamics Simulations. J Phys Chem B 2022; 126:10922-10932. [PMID: 36516319 PMCID: PMC9806834 DOI: 10.1021/acs.jpcb.2c06947] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Simulations of ab initio molecular dynamics have been performed for mixtures of ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI) ionic liquid and water. Statistics of donors and acceptors of hydrogen bonds has revealed that with increasing water content, hydrogen bonds between EMIM cations and TFSI anions are replaced by bonds to water molecules. In the mixture of liquids, the total number of bonds (from EMIM cations or water molecules) formed by TFSI acceptors increases. IR spectra obtained from ab initio molecular dynamics trajectories are in good agreement with literature data for ionic liquid/water systems. Analysis of oscillations of individual C-H and O-H bonds has shown correlations between vibrational frequencies and hydrogen bonds formed by an EMIM cation or water molecule and has indicated that the changes in the IR spectrum result from the decreased number of water-water hydrogen bonds in the mixture. The tests of DFTB methodology with tailored parameterizations have yielded reasonably good description of the IR spectrum of bulk water, whereas available parameterizations have failed in satisfactory reproduction of the IR spectrum of EMIM-TFSI/water mixtures in the region above 3000 cm-1.
Collapse
|
6
|
Gregory KP, Elliott GR, Wanless EJ, Webber GB, Page AJ. A quantum chemical molecular dynamics repository of solvated ions. Sci Data 2022; 9:430. [PMID: 35864118 PMCID: PMC9304403 DOI: 10.1038/s41597-022-01527-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2022] [Accepted: 06/30/2022] [Indexed: 12/29/2022] Open
Abstract
The importance of ion-solvent interactions in predicting specific ion effects in contexts ranging from viral activity through to electrolyte viscosity cannot be underestimated. Moreover, investigations of specific ion effects in nonaqueous systems, highly relevant to battery technologies, biochemical systems and colloid science, are severely limited by data deficiency. Here, we report IonSolvR – a collection of more than 3,000 distinct nanosecond-scale ab initio molecular dynamics simulations of ions in aqueous and non-aqueous solvent environments at varying effective concentrations. Density functional tight binding (DFTB) is used to detail the solvation structure of up to 55 solutes in 28 different protic and aprotic solvents. DFTB is a fast quantum chemical method, and as such enables us to bridge the gap between efficient computational scaling and maintaining accuracy, while using an internally-consistent simulation technique. We validate the database against experimental data and provide guidance for accessing individual IonSolvR records. Measurement(s) | solvation structure | Technology Type(s) | quantum chemistry computational method • Molecular Dynamics |
Collapse
Affiliation(s)
- Kasimir P Gregory
- Discipline of Chemistry, School of Environmental & Life Sciences, University of Newcastle, Callaghan, NSW, 2308, Australia.,Department of Materials Physics, Research School of Physics, Australian National University, Canberra, ACT, 0200, Australia
| | - Gareth R Elliott
- Discipline of Chemistry, School of Environmental & Life Sciences, University of Newcastle, Callaghan, NSW, 2308, Australia
| | - Erica J Wanless
- Discipline of Chemistry, School of Environmental & Life Sciences, University of Newcastle, Callaghan, NSW, 2308, Australia
| | - Grant B Webber
- Discipline of Chemical Engineering, School of Engineering, University of Newcastle, Callaghan, NSW, 2308, Australia
| | - Alister J Page
- Discipline of Chemistry, School of Environmental & Life Sciences, University of Newcastle, Callaghan, NSW, 2308, Australia.
| |
Collapse
|
7
|
Lourenço MP, Herrera LB, Hostaš J, Calaminici P, Köster AM, Tchagang A, Salahub DR. A new active learning approach for adsorbate-substrate structural elucidation in silico. J Mol Model 2022; 28:178. [PMID: 35654918 DOI: 10.1007/s00894-022-05173-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Accepted: 05/19/2022] [Indexed: 11/29/2022]
Abstract
Adsorbate interactions with substrates (e.g. surfaces and nanoparticles) are fundamental for several technologies, such as functional materials, supramolecular chemistry, and solvent interactions. However, modeling these kinds of systems in silico, such as finding the optimum adsorption geometry and energy, is challenging, due to the huge number of possibilities of assembling the adsorbate on the surface. In the current work, we have developed an artificial intelligence (AI) approach based on an active learning (AL) method for adsorption optimization on the surface of materials. AL uses machine learning (ML) regression algorithms and their uncertainties to make a decision (based on a policy) for the next unexplored structures to be computed, increasing, though, the probability of finding the global minimum with a small number of calculations. The methodology allows an accurate and automated structural elucidation of the adsorbate on the surface, based on the minimization of the total electronic energy. The new AL method for adsorption optimization was developed and implemented in the quantum machine learning software/agent for material design and discovery (QMLMaterial) program and was applied for C60@TiO2 anatase (101). It marks another software extension with a new feature in addition to the automatic structural elucidation of defects in materials and of nanoparticles as well. SCC-DFTB calculations were used to build the complex search surfaces with a reasonably low computational cost. An artificial neural network (NN) was employed in the AL framework evaluated together with two uncertainty quantification methods: K-fold cross-validation and non-parametric bootstrap (BS) resampling. Also, two different acquisition functions for decision-making were used: expected improvement (EI) and the lower confidence bound (LCB).
Collapse
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.
| | - Lizandra Barrios Herrera
- Department of Chemistry, Department of Physics and Astronomy, Quantum Alberta, CMS Centre for Molecular Simulation, IQST Institute for Quantum Science and Technology, University of Calgary, 2500 University Drive NW, Calgary, AB, T2N 1N4, Canada
| | - Jiří Hostaš
- Department of Chemistry, Department of Physics and Astronomy, Quantum Alberta, CMS Centre for Molecular Simulation, IQST Institute for Quantum Science and Technology, University of Calgary, 2500 University Drive NW, Calgary, AB, T2N 1N4, Canada
| | - Patrizia Calaminici
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional 2508, AP 14-740, México City, D.F., 07000, México
| | - Andreas M Köster
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional 2508, AP 14-740, México City, D.F., 07000, México
| | - Alain Tchagang
- Digital Technologies Research Centre, National Research Council of Canada, 1200 Montréal Road, Ottawa, ON, K1A 0R6, Canada
| | - Dennis R Salahub
- Department of Chemistry, Department of Physics and Astronomy, Quantum Alberta, CMS Centre for Molecular Simulation, IQST Institute for Quantum Science and Technology, University of Calgary, 2500 University Drive NW, Calgary, AB, T2N 1N4, Canada
| |
Collapse
|
8
|
Lindsey RK, Huy Pham C, Goldman N, Bastea S, Fried LE. Machine‐Learning a Solution for Reactive Atomistic Simulations of Energetic Materials. PROPELLANTS EXPLOSIVES PYROTECHNICS 2022. [DOI: 10.1002/prep.202200001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Rebecca K. Lindsey
- Physical and Life Sciences Directorate Lawrence Livermore National Laboratory Livermore California 94550 USA
| | - Cong Huy Pham
- Physical and Life Sciences Directorate Lawrence Livermore National Laboratory Livermore California 94550 USA
| | - Nir Goldman
- Physical and Life Sciences Directorate Lawrence Livermore National Laboratory Livermore California 94550 USA
- Department of Chemical Engineering University of California, Davis Davis California 95616 USA
| | - Sorin Bastea
- Physical and Life Sciences Directorate Lawrence Livermore National Laboratory Livermore California 94550 USA
| | - Laurence E. Fried
- Physical and Life Sciences Directorate Lawrence Livermore National Laboratory Livermore California 94550 USA
| |
Collapse
|
9
|
Goldman N, Kweon KE, Sadigh B, Heo TW, Lindsey RK, Pham CH, Fried LE, Aradi B, Holliday K, Jeffries JR, Wood BC. Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields. J Chem Theory Comput 2021; 17:4435-4448. [PMID: 34128678 DOI: 10.1021/acs.jctc.1c00172] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. We have focused our workflow on TiH2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.
Collapse
Affiliation(s)
- Nir Goldman
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States.,Department of Chemical Engineering, University of California, Davis, Davis, California 95616, United States
| | - Kyoung E Kweon
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Babak Sadigh
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Tae Wook Heo
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Rebecca K Lindsey
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - C Huy Pham
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Laurence E Fried
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, Universität Bremen, P.O.B. 330440, Bremen D-28334, Germany
| | - Kiel Holliday
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Jason R Jeffries
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Brandon C Wood
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| |
Collapse
|
10
|
Galvão BRL, Viegas LP, Salahub DR, Lourenço MP. Reliability of semiempirical and DFTB methods for the global optimization of the structures of nanoclusters. J Mol Model 2020; 26:303. [PMID: 33064203 DOI: 10.1007/s00894-020-04484-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2020] [Accepted: 07/21/2020] [Indexed: 11/26/2022]
Abstract
In this work, we explore the possibility of using computationally inexpensive electronic structure methods, such as semiempirical and DFTB calculations, for the search of the global minimum (GM) structure of chemical systems. The basic prerequisite that these inexpensive methods will need to fulfill is that their lowest energy structures can be used as starting point for a subsequent local optimization at a benchmark level that will yield its GM. If this is possible, one could bypass the global optimization at the expensive method, which is currently impossible except for very small molecules. Specifically, we test our methods with clusters of second row elements including systems of several bonding types, such as alkali, metal, and covalent clusters. The results reveal that the DFTB3 method yields reasonable results and is a potential candidate for this type of applications. Even though the DFTB2 approach using standard parameters is proven to yield poor results, we show that a re-parametrization of only its repulsive part is enough to achieve excellent results, even when applied to larger systems outside the training set.
Collapse
Affiliation(s)
- Breno R L Galvão
- Centro Federal de Educação Tecnológica de Minas Gerais, CEFET-MG, Av. Amazonas 5253, Belo Horizonte, Minas Gerais, 30421-169, Brazil.
| | - Luís P Viegas
- Coimbra Chemistry Center and Chemistry Department, University of Coimbra, 3004-535, Coimbra, Portugal
| | - Dennis R Salahub
- Department of Chemistry, CMS - Centre for Molecular Simulation, IQST - Institute for Quantum Science and Technology, University of Calgary, 2500 University Drive NW, Calgary, AB, T2N 1N4, Canada
| | - Maicon P 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, Alegre, Espírito Santo, 29500-000, Brazil
| |
Collapse
|
11
|
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
|
12
|
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
|