1
|
Dey S, Folkestad SD, Paul AC, Koch H, Krylov AI. Core-ionization spectrum of liquid water. Phys Chem Chem Phys 2024; 26:1845-1859. [PMID: 38174659 DOI: 10.1039/d3cp02499g] [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/2024]
Abstract
We present state-of-the-art calculations of the core-ionization spectrum of water. Despite significant progress in procedures developed to mitigate various experimental complications and uncertainties, the experimental determination of ionization energies of solvated species involves several non-trivial steps such as assessing the effect of the surface potential, electrolytes, and finite escape depths of photoelectrons. This provides a motivation to obtain robust theoretical values of the intrinsic bulk ionization energy and the corresponding solvent-induced shift. Here we develop theoretical protocols based on coupled-cluster theory and electrostatic embedding. Our value of the intrinsic solvent-induced shift of the 1sO ionization energy of water is -1.79 eV. The computed absolute position and the width of the 1sO peak in photoelectron spectrum of water are 538.47 eV and 1.44 eV, respectively, agreeing well with the best experimental values.
Collapse
Affiliation(s)
- Sourav Dey
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA.
| | - Sarai Dery Folkestad
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA.
| | - Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA.
| |
Collapse
|
2
|
Mazlee MTF, Heidelberg T, Ariffin A, Zain SM. Cation-stimulated drug delivery via lipid assembly comprising macrocyclized disaccharides - A DFT study. Carbohydr Res 2023; 532:108923. [PMID: 37598565 DOI: 10.1016/j.carres.2023.108923] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2023] [Revised: 08/07/2023] [Accepted: 08/09/2023] [Indexed: 08/22/2023]
Abstract
In the attempt to create a delivery system for an alkali-cation stimulated drug release, a computational study was conducted, aiming for the evaluation of synthetic access towards glycolipid crown ethers analogs and their potential for coordination-induced changes of packing constraints for molecular assemblies. The results disfavor amide-linkages for the creation of macrocycles around the inter-glycosidic bond of a disaccharide. Conformational changes upon cation coordination of the macrocycle decrease the intersection area for easily accessible macrocycles based on lactose. This leads to shrinking intersection areas upon alkali complexation. Maltose-based analogs, on the other hand, exhibited the targeted increase of the glycolipid intersection area and, hence, may be considered as a promising resource.
Collapse
Affiliation(s)
| | - Thorsten Heidelberg
- Chemistry Department, Faculty of Science, University of Malaya, 50603, Kuala Lumpur, Malaysia.
| | - Azhar Ariffin
- Chemistry Department, Faculty of Science, University of Malaya, 50603, Kuala Lumpur, Malaysia
| | - Sharifuddin Md Zain
- Chemistry Department, Faculty of Science, University of Malaya, 50603, Kuala Lumpur, Malaysia
| |
Collapse
|
3
|
Huo J, Chen J, Liu P, Hong B, Zhang J, Dong H, Li S. Microscopic Mechanism of Proton Transfer in Pure Water under Ambient Conditions. J Chem Theory Comput 2023. [PMID: 37365994 DOI: 10.1021/acs.jctc.3c00244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2023]
Abstract
Water molecules and the associated proton transfer (PT) are prevalent in chemical and biological systems and have been a hot research topic. Spectroscopic characterization and ab initio molecular dynamics (AIMD) simulations have previously revealed insights into acidic and basic liquids. Presumably, the situation in the acidic/basic solution is not necessarily the same as in pure water; in addition, the autoionization constant for water is only 10-14 under ambient conditions, making the study of PT in pure water challenging. To overcome this issue, we modeled periodic water box systems containing 1000 molecules for tens of nanoseconds based on a neural network potential (NNP) with quantum mechanical accuracy. The NNP was generated by training a dataset containing the energies and atomic forces of 17 075 configurations of periodic water box systems, and these data points were calculated at the MP2 level that considers electron correlation effects. We found that the size of the system and the duration of the simulation have a significant impact on the convergence of the results. With these factors considered, our simulations showed that hydronium (H3O+) and hydroxide (OH-) ions in water have distinct hydration structures, thermodynamic and kinetic properties, e.g., the longer-lasting and more stable hydrated structure of OH- ions than that of H3O+, as well as a significantly higher free energy barrier for the OH--associated PT than that of H3O+, leading the two to exhibit completely different PT behaviors. Given these characteristics, we further found that PT via OH- ions tends not to occur multiple times or between many molecules. In contrast, PT via H3O+ can synergistically occur among multiple molecules and prefers to adopt a cyclic pattern among three water molecules, while it occurs mostly in a chain pattern when more water molecules are involved. Therefore, our studies provide a detailed and solid microscopic explanation for the PT process in pure water.
Collapse
Affiliation(s)
- Jun Huo
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China
| | - Jianghao Chen
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China
- School of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Pei Liu
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| | - Benkun Hong
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| | - Jian Zhang
- School of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Hao Dong
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China
- State Key Laboratory of Analytical Chemistry for Life Science, Nanjing University, Nanjing 210023, China
- Institute for Brain Sciences, Nanjing University, Nanjing 210023, China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| |
Collapse
|
4
|
Garcia Carcamo RA, Zhang X, Estejab A, Zhou J, Hare BJ, Sievers C, Sarupria S, Getman RB. Differences in solvation thermodynamics of oxygenates at Pt/Al 2O 3 perimeter versus Pt(111) terrace sites. iScience 2023; 26:105980. [PMID: 36756373 PMCID: PMC9900392 DOI: 10.1016/j.isci.2023.105980] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2022] [Revised: 12/26/2022] [Accepted: 01/11/2023] [Indexed: 01/19/2023] Open
Abstract
A prominent role of water in aqueous-phase heterogeneous catalysis is to modify free energies; however, intuition about how is based largely on pure metal surfaces or even homogeneous solutions. Using multiscale modeling with explicit liquid water molecules, we show that the influence of water on the free energies of adsorbates at metal/support interfaces is different than that on pure metal surfaces. We specifically compute free energies of solvation for methanol and its constituents on a Pt/Al2O3 catalyst and compare the results to analogous values calculated on a pure Pt catalyst. We find that the more hydrophilic Pt/Al2O3 interface leads to smaller (more positive) free energies of solvation due to an increased entropy penalty resulting from the additional work necessary to disrupt the interfacial water structure and accommodate the interfacial species. The results will be of interest in other fields, including adsorption and proteins.
Collapse
Affiliation(s)
| | - Xiaohong Zhang
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC 29634, USA
| | - Ali Estejab
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC 29634, USA
| | - Jiarun Zhou
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC 29634, USA
| | - Bryan J. Hare
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - Carsten Sievers
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - Sapna Sarupria
- Department of Chemistry and Chemical Theory Center, University of Minnesota, Minneapolis, MN 55455, USA
| | - Rachel B. Getman
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC 29634, USA,Corresponding author
| |
Collapse
|
5
|
Liu J, He X. Recent advances in quantum fragmentation approaches to complex molecular and condensed‐phase systems. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy China Pharmaceutical University Nanjing China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering East China Normal University Shanghai China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering East China Normal University Shanghai China
- New York University‐East China Normal University Center for Computational Chemistry New York University Shanghai Shanghai China
| |
Collapse
|
6
|
Liu J, Lan J, He X. Toward High-level Machine Learning Potential for Water Based on Quantum Fragmentation and Neural Networks. J Phys Chem A 2022; 126:3926-3936. [PMID: 35679610 DOI: 10.1021/acs.jpca.2c00601] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Accurate and efficient simulation of liquids, such as water and salt solutions, using high-level wave function theories is still a formidable task for computational chemists owing to the high computational costs. In this study, we develop a deep machine learning potential based on fragment-based second-order Møller-Plesset perturbation theory (DP-MP2) for water through neural networks. We show that the DP-MP2 potential predicts the structural, dynamical, and thermodynamic properties of liquid water in better agreement with the experimental data than previous studies based on density functional theory (DFT). The nuclear quantum effects (NQEs) on the properties of liquid water are also examined, which are noticeable in affecting the structural and dynamical properties of liquid water under ambient conditions. This work provides a general framework for quantitative predictions of the properties of condensed-phase systems with the accuracy of high-level wave function theory while achieving significant computational savings compared to ab initio simulations.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China.,Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Jinggang Lan
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,New York University-East China Normal University Center for Computational Chemistry, NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
7
|
Gao Y, Fang H, Ni K, Feng Y. Water clusters and density fluctuations in liquid water based on extended hierarchical clustering methods. Sci Rep 2022; 12:8036. [PMID: 35577839 PMCID: PMC9110331 DOI: 10.1038/s41598-022-11947-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2022] [Accepted: 04/25/2022] [Indexed: 11/29/2022] Open
Abstract
The microscopic structures of liquid water at ambient temperatures remain a hot debate, which relates with structural and density fluctuations in the hydrogen bond network. Here, we use molecular dynamics simulations of liquid water to study the properties of three-dimensional cage-like water clusters, which we investigate using extended graph-based hierarchical clustering methods. The water clusters can cover over 95% of hydrogen bond network, among which some clusters maximally encompass thousands of molecules extending beyond 3.0 nm. The clusters imply fractal behaviors forming percolating networks and the morphologies of small and large clusters show different scaling rules. The local favored clusters and the preferred connections between adjacent clusters correspond to lower energy and conformational entropy depending on cluster topologies. Temperature can destroy large clusters into small ones. We show further that the interior of clusters favors high-density patches. The water molecules in the small clusters, inside which are the void regarded as hydrophobic objects, have a preference for being more tetrahedral. Our results highlight the properties and changes of water clusters as the fundamental building blocks of hydrogen bond networks. In addition, the water clusters can elucidate structural and density fluctuations on different length scales in liquid water.
Collapse
|
8
|
Kurskaya EA, Podorozhko EA, Afanasyev ES, Kononova EG, Askadskii AA. Trends in Cryotropic Gelation of Semidilute Aqueous Solutions of Poly(vinyl alcohol) with Different Thermal History. POLYMER SCIENCE SERIES A 2022. [DOI: 10.1134/s0965545x22010060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
9
|
Liu J, He X. Ab initio molecular dynamics simulation of liquid water with fragment-based quantum mechanical approach under periodic boundary conditions. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- New York University-East China Normal University Center for Computational Chemistry at New York University Shanghai, Shanghai 200062, China
| |
Collapse
|
10
|
Heindel JP, Xantheas SS. Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD). J Chem Theory Comput 2021; 17:7341-7352. [PMID: 34723531 DOI: 10.1021/acs.jctc.1c00780] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-molecular dynamics (MD) classical and nuclear quantum dynamical simulations, in which the energies and forces of the full system are approximated by the two-, three-, and four-body terms of the MBE, and compare the average potential and the vibrational density of states with the full simulation, i.e., the one for which no MBE is used. Our results indicate that the thermally averaged potential energy from the MBE up to the four-body term converges with near-identical behavior to the one from the full simulation. The three-body makes a substantial contribution (∼20%) to the energy, whereas the four-body is necessary for obtaining quantitatively accurate energetics and forces, albeit making a small contribution to each (∼2%). We further show that the harmonic frequencies are reproduced to within a few wavenumbers (cm-1) at the four-body level and that the slowest modes to converge with the MBE rank are those involving the strongest hydrogen bonds. Anharmonicity exacerbates this effect, so that a four-body description of the energies and forces is needed to achieve accurate anharmonic vibrational frequencies in the hydrogen-bonded OH-stretching region. We also discuss the asymptotic scaling of the MBE-MD protocol with respect to the cost of the underlying potential energy evaluation, suggesting that electronic structure methods that scale at least as N4, N being the size of the system, are needed to result in savings over the traditional full MD simulation. We anticipate that the MBE-MD protocol can evolve into a powerful and practical method, which will allow for highly accurate ab initio MD simulations on a much broader range of molecular systems than can be currently handled.
Collapse
Affiliation(s)
- Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States.,Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| |
Collapse
|
11
|
Towards complete assignment of the infrared spectrum of the protonated water cluster H +(H 2O) 21. Nat Commun 2021; 12:6141. [PMID: 34686665 PMCID: PMC8536673 DOI: 10.1038/s41467-021-26284-x] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2021] [Accepted: 09/22/2021] [Indexed: 11/08/2022] Open
Abstract
The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H+(H2O)n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H+(H2O)21 cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster. Protonated water species have been the subject of numerous experimental and computational studies. Here the authors provide a nearly complete assignment of the experimental IR spectrum of the H+(H2O)21 water cluster based on high-level wavefunction theory and anharmonic vibrational quasi-degenerate perturbation theory.
Collapse
|
12
|
Yao Y, Kanai Y. Nuclear Quantum Effect and Its Temperature Dependence in Liquid Water from Random Phase Approximation via Artificial Neural Network. J Phys Chem Lett 2021; 12:6354-6362. [PMID: 34231366 DOI: 10.1021/acs.jpclett.1c01566] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We report structural and dynamical properties of liquid water described by the random phase approximation (RPA) correlation together with the exact exchange energy (EXX) within density functional theory. By utilizing thermostated ring polymer molecular dynamics, we examine the nuclear quantum effects and their temperature dependence. We circumvent the computational limitation of performing direct first-principles molecular dynamics simulation at this high level of electronic structure theory by adapting an artificial neural network model. We show that the EXX+RPA level of theory accurately describes liquid water in terms of both dynamical and structural properties.
Collapse
Affiliation(s)
- Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Durham, North Carolina 27599, United States
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Durham, North Carolina 27599, United States
| |
Collapse
|
13
|
Eriksen JJ. Decomposed Mean-Field Simulations of Local Properties in Condensed Phases. J Phys Chem Lett 2021; 12:6048-6055. [PMID: 34165982 DOI: 10.1021/acs.jpclett.1c01375] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The present work demonstrates a robust protocol for probing localized electronic structure in condensed-phase systems, operating in terms of a recently proposed theory for decomposing the results of Kohn-Sham density functional theory in a basis of spatially localized molecular orbitals. In an initial application to liquid, ambient water and the assessment of the solvation energy and the embedded dipole moment of H2O in solution, we find that both properties are amplified on average-in accordance with expectation-and that correlations are indeed observed to exist between them. However, the simulated solvent-induced shift to the dipole moment of water is found to be significantly dampened with respect to typical literature values. The local nature of our methodology has further allowed us to evaluate the convergence of bulk properties with respect to the extent of the underlying one-electron basis set, ranging from single-ζ to full (augmented) quadruple-ζ quality. Albeit a pilot example, our work paves the way toward future studies of local effects and defects in more complex phases, e.g., liquid mixtures and even solid-state crystals.
Collapse
Affiliation(s)
- Janus J Eriksen
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg. 206, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
14
|
Liu J, Liu Y, Yang J, Zeng XC, He X. Directional Proton Transfer in the Reaction of the Simplest Criegee Intermediate with Water Involving the Formation of Transient H 3O . J Phys Chem Lett 2021; 12:3379-3386. [PMID: 33784110 DOI: 10.1021/acs.jpclett.1c00448] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The reaction of Criegee intermediates with water vapor has been widely known as a key Criegee reaction in the troposphere. Herein, we investigated the reaction of the smallest Criegee intermediate, CH2OO, with a water cluster through fragment-based ab initio molecular dynamics simulations at the MP2/aug-cc-pVDZ level. Our results show that the CH2OO-water reaction could occur not only at the air/water interface but also inside the water cluster. Moreover, more than one reactive water molecules are required for the CH2OO-water reaction, which is always initiated from the Criegee carbon atom and ends at the terminal Criegee oxygen atom via a directional proton transfer process. The observed reaction pathways include the loop-structure-mediated and stepwise mechanisms, and the latter involves the formation of transient H3O+. The lifetime of transient H3O+ is on the order of a few picoseconds, which may impact the atmospheric budget of the other trace gases in the actual atmosphere.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Yanqing Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - Jinrong Yang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Xiao Cheng Zeng
- Department of Chemistry, University of Nebraska, Lincoln, Nebraska 68588, United States
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU, Shanghai, Shanghai, 200062, China
| |
Collapse
|
15
|
Hao X, Liu J, Ali I, Luo H, Han Y, Hu W, Liu J, He X, Li J. Ab initio determination of crystal stability of di-p-tolyl disulfide. Sci Rep 2021; 11:7076. [PMID: 33782489 PMCID: PMC8007795 DOI: 10.1038/s41598-021-86519-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2020] [Accepted: 03/16/2021] [Indexed: 11/09/2022] Open
Abstract
With the rapid growth of energy demand and the depletion of existing energy resources, the new materials with superior performances, low costs and environmental friendliness for energy production and storage are explored. Di-p-tolyl disulfide (p-Tol2S2) is a typical lubricating material, which has been applied in the field of energy storage. The conformational properties and phase transformations of p-Tol2S2 have been studied by pioneers, but their polymorphs and the polymorphism induced crystal structure changes require further analysis. In this study, we perform the crystal structural screening, prediction and optimization of p-Tol2S2 crystal with quantum mechanical calculations, i.e., density functional theory (DFT) and second-order Møller–Plesset perturbation (MP2) methods. A series of crystal structures with different molecular arrangements are generated based on the crystal structure screening. As compared to long-established lattice energy calculation, we take an advantage of using more accurate technique, which is Gibbs free energy calculation. It considers the effects of entropy and temperature to predict the crystal structures and energy landscape. By comparing the Gibbs free energies between predicted and experimental structures, we found that phase α is the most stable structure for p-Tol2S2 crystal at ambient temperature and standard atmospheric pressure. Furthermore, we provide an efficient method to discriminate different polymorphs that are otherwise difficult to be identified based on the Raman/IR spectra. The proposed work enable us to evaluate the quality of various crystal polymorphs rapidly.
Collapse
Affiliation(s)
- Xuan Hao
- School of Chemistry and Molecular Engineering, Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, East China Normal University, Shanghai, 200062, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jinfeng Liu
- State Key Laboratory of Natural Medicines, Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - Imran Ali
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Hongyuan Luo
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Yanqiang Han
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Wenxin Hu
- School of Computer Science and Software Engineering, The Computer Center, East China Normal University, Shanghai, 200062, China
| | - Jinyun Liu
- Anhui Provincial Engineering Laboratory for New-Energy Vehicle Battery Energy-Storage Materials, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu, 241002, Anhui, China.
| | - Xiao He
- School of Chemistry and Molecular Engineering, Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, East China Normal University, Shanghai, 200062, China. .,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China.
| | - Jinjin Li
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China.
| |
Collapse
|
16
|
Zhu T, Van Voorhis T. Understanding the Dipole Moment of Liquid Water from a Self-Attractive Hartree Decomposition. J Phys Chem Lett 2021; 12:6-12. [PMID: 33296211 DOI: 10.1021/acs.jpclett.0c03300] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The dipole moment of a single water molecule in liquid water has been a critical concept for understanding water's dielectric properties. In this work, we investigate the dipole moment of liquid water through a self-attractive Hartree (SAH) decomposition of total electron density computed by density functional theory, on water clusters sampled from ab initio molecular dynamics simulation of bulk water. By adjusting one parameter that controls the degree of density localization, we reveal two distinct pictures of water dipoles that are consistent with bulk dielectric properties: a localized picture with smaller and less polarizable monomer dipoles and a delocalized picture with larger and more polarizable monomer dipoles. We further uncover that the collective dipole-dipole correlation is stronger in the localized picture and is key to connecting individual dipoles with bulk dielectric properties. On the basis of these findings, we suggest considering both individual and collective dipole behaviors when studying the dipole moment of liquid water and propose new design strategies for developing water models.
Collapse
Affiliation(s)
- Tianyu Zhu
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
17
|
de Kock MB, Azim S, Kassier GH, Miller RJD. Determining the radial distribution function of water using electron scattering: A key to solution phase chemistry. J Chem Phys 2020; 153:194504. [DOI: 10.1063/5.0024127] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- M. B. de Kock
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, Bldg. 99 (CFEL), 22761 Hamburg, Germany
| | - S. Azim
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, Bldg. 99 (CFEL), 22761 Hamburg, Germany
| | - G. H. Kassier
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, Bldg. 99 (CFEL), 22761 Hamburg, Germany
| | - R. J. D. Miller
- Departments of Chemistry and Physics, University of Toronto, 80 St. George Street, Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
18
|
|
19
|
Jin X, Glover WJ, He X. Fragment Quantum Mechanical Method for Excited States of Proteins: Development and Application to the Green Fluorescent Protein. J Chem Theory Comput 2020; 16:5174-5188. [PMID: 32551640 DOI: 10.1021/acs.jctc.9b00980] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Understanding the excited-state properties of luminescent biomolecules is of central importance to their biophysical applications. In this study, we develop the Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps (EE-GMFCC) method for quantitatively characterizing properties of covalently bonded systems with localized excitations (i.e., involving a single chromophore), such as fluorescent proteins. The excitation energy, transition dipole moment, and oscillator strength of wild-type Green Fluorescent Protein (wt-GFP) calculated by EE-GMFCC are found to be in excellent agreement with full system time-dependent density functional theory results. We also applied the Polarized Protein-Specific Charge model to wt-GFP, and found that electronic polarization of the protein is critical in stabilizing hydrogen bonding interactions in wt-GFP, which influences its absorption spectrum. The predicted absorption spectra of wt-GFP in the A and B states qualitatively agree with experiment. The fragmentation approach further allows a straightforward per residue decomposition of the excitation which reveals the influence of the protein environment on the absorption spectra of wt-GFP A and B states. Our results demonstrate that the EE-GMFCC method is both accurate and efficient for excited-state property calculations on proteins.
Collapse
Affiliation(s)
- Xinsheng Jin
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - William J Glover
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University, New York, New York 10003, United States
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
20
|
Ricard TC, Iyengar SS. Efficient and Accurate Approach To Estimate Hybrid Functional and Large Basis-Set Contributions to Condensed-Phase Systems and Molecule–Surface Interactions. J Chem Theory Comput 2020; 16:4790-4812. [DOI: 10.1021/acs.jctc.9b01089] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Timothy C. Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
21
|
Liu J, He X. Fragment-based quantum mechanical approach to biomolecules, molecular clusters, molecular crystals and liquids. Phys Chem Chem Phys 2020; 22:12341-12367. [PMID: 32459230 DOI: 10.1039/d0cp01095b] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
To study large molecular systems beyond the system size that the current state-of-the-art ab initio electronic structure methods could handle, fragment-based quantum mechanical (QM) approaches have been developed over the past years, and proved to be efficient in dealing with large molecular systems at various ab initio levels. According to the fragmentation approach, a large molecular system can be divided into subsystems (fragments), and subsequently the property of the whole system can be approximately obtained by taking a proper combination of the corresponding terms of individual fragments. Therefore, the standard QM calculation of a large system could be circumvented by carrying out a series of calculations on small fragments, which significantly promotes computational efficiency. The electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method is one of the fragment-based QM approaches which has been developed by our research group in recent years. This Perspective presents the theoretical framework of this fragmentation method and its applications in biomolecules, molecular clusters, molecular crystals and liquids, including total energy calculation, protein-ligand/protein binding affinity prediction, geometry optimization, vibrational spectrum simulation, ab initio molecular dynamics simulation, and prediction of excited-state properties.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | | |
Collapse
|
22
|
Huang L, Han Y, Liu J, He X, Li J. Ab Initio Prediction of the Phase Transition for Solid Ammonia at High Pressures. Sci Rep 2020; 10:7546. [PMID: 32372007 PMCID: PMC7200730 DOI: 10.1038/s41598-020-64030-3] [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/22/2019] [Accepted: 04/08/2020] [Indexed: 11/11/2022] Open
Abstract
Ammonia is one of the most basic components on the planet and its high-pressure characteristics play an important role in planetary science. Solid ammonia crystals frequently adopt multiple distinct polymorphs exhibiting different properties. Predicting the crystal structure of these polymorphs and under what thermodynamic conditions these polymorphs are stable would be of great value to environmental industry and other fields. Theoretical calculations based on the classical force fields and density-functional theory (DFT) are versatile methods but lack of accurate description of weak intermolecular interactions for molecular crystals. In this study, we employ an ab initio computational study on the solid ammonia at high pressures, using the second-order Møller-Plesset perturbation (MP2) theory and the coupled cluster singles, doubles, and perturbative triples (CCSD(T)) theory along with the embedded fragmentation method. The proposed algorithm is capable of performing large-scale calculations using high-level wavefunction theories, and accurately describing covalent, ionic, hydrogen bonding, and dispersion interactions within molecular crystals, and therefore can predict the crystal structures, Raman spectra and phase transition of solid ammonia phases I and IV accurately. We confirm the crystal structures of solid ammonia phases I and IV that have been controversial for a long time and predict their phase transition that occurs at 1.17 GPa and 210 K with small temperature dependence, which is in line with experiment.
Collapse
Affiliation(s)
- Lei Huang
- Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Yanqiang Han
- Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jinyun Liu
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu, Anhui, 241000, China.
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China.
| | - Jinjin Li
- Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai, 200240, China.
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu, Anhui, 241000, China.
| |
Collapse
|
23
|
Zhang W, Liu J, Jin X, Gu X, Zeng XC, He X, Li H. Quantitative Prediction of Aggregation‐Induced Emission: A Full Quantum Mechanical Approach to the Optical Spectra. Angew Chem Int Ed Engl 2020; 59:11550-11555. [DOI: 10.1002/anie.202003326] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Indexed: 12/23/2022]
Affiliation(s)
- Wei Zhang
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
- Department of Chemistry University of Nebraska Lincoln NE 68588 USA
| | - Jinfeng Liu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
- Department of Basic Medicine and Clinical Pharmacy China Pharmaceutical University Nanjing 210009 China
| | - Xinsheng Jin
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
| | - Xinggui Gu
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
| | - Xiao Cheng Zeng
- Department of Chemistry University of Nebraska Lincoln NE 68588 USA
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai Shanghai 200062 China
| | - Hui Li
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
| |
Collapse
|
24
|
Zhang W, Liu J, Jin X, Gu X, Zeng XC, He X, Li H. Quantitative Prediction of Aggregation‐Induced Emission: A Full Quantum Mechanical Approach to the Optical Spectra. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202003326] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Wei Zhang
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
- Department of Chemistry University of Nebraska Lincoln NE 68588 USA
| | - Jinfeng Liu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
- Department of Basic Medicine and Clinical Pharmacy China Pharmaceutical University Nanjing 210009 China
| | - Xinsheng Jin
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
| | - Xinggui Gu
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
| | - Xiao Cheng Zeng
- Department of Chemistry University of Nebraska Lincoln NE 68588 USA
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai Shanghai 200062 China
| | - Hui Li
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
| |
Collapse
|
25
|
Wang Z, Han Y, Li J, He X. Combining the Fragmentation Approach and Neural Network Potential Energy Surfaces of Fragments for Accurate Calculation of Protein Energy. J Phys Chem B 2020; 124:3027-3035. [PMID: 32208716 DOI: 10.1021/acs.jpcb.0c01370] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Accurate and efficient all-atom quantum mechanical (QM) calculations for biomolecules still present a challenge to computational physicists and chemists. In this study, an extensible generalized molecular fractionation with a conjugate caps method combined with neural networks (NN-GMFCC) is developed for efficient QM calculation of protein energy. In the NN-GMFCC scheme, the total energy of a given protein is calculated by taking a proper combination of the high-precision neural network potential energies of all capped residues and overlapping conjugate caps. In addition, the two-body interaction energies of residue pairs are calculated by molecular mechanics (MM). With reference to the GMFCC/MM calculation at the ωB97XD/6-31G* level, the overall mean unsigned errors of the energy deviations and atomic force root-mean-squared errors calculated by NN-GMFCC are only 2.01 kcal/mol and 0.68 kcal/mol/Å, respectively, for 14 proteins (containing up to 13,728 atoms). Meanwhile, the NN-GMFCC approach is about 4 orders of magnitude faster than the GMFCC/MM method. The NN-GMFCC method could be systematically improved by inclusion of two-body QM interaction and multibody electronic polarization effect. Moreover, the NN-GMFCC approach can also be applied to other macromolecular systems such as DNA/RNA, and it is capable of providing a powerful and efficient approach for exploration of structures and functions of proteins with QM accuracy.
Collapse
Affiliation(s)
- Zhilong Wang
- Key Laboratory of Thin Film and Micro Fabrication, Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Yanqiang Han
- Key Laboratory of Thin Film and Micro Fabrication, Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Jinjin Li
- Key Laboratory of Thin Film and Micro Fabrication, Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
26
|
Self-assembly Processes in Hydrated Montmorillonite by FTIR Investigations. MATERIALS 2020; 13:ma13051100. [PMID: 32121630 PMCID: PMC7084999 DOI: 10.3390/ma13051100] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/04/2020] [Revised: 02/21/2020] [Accepted: 02/27/2020] [Indexed: 11/17/2022]
Abstract
Experimental findings obtained by FTIR and Raman spectroscopies on montmorillonite-water mixtures at three concentration values are presented. To get some insight into the hydrogen bond network of water within the montmorillonite network, FTIR and Raman spectra have been collected as a function of time and then analyzed following two complementary approaches: An analysis of the intramolecular OH stretching mode in the spectral range of 2700–3900 cm−1 in terms of two Gaussian components, and an analysis of the same OH stretching mode by wavelet cross-correlation. The FTIR and Raman investigations have been carried as a function of time for a montmorillonite-water weight composition (wt%) of 20–80%, 25–75%, and 35–65%, until the dehydrated state where the samples appear as a homogeneous rigid layer of clay. In particular, for both the FTIR and Raman spectra, the decomposition of the OH stretching band into a “closed” and an “open” contribution and the spectral wavelet analysis allow us to extract quantitative information on the time behavior of the system water content. It emerges that, the total water contribution inside the montmorillonite structure decreases as a function of time. However, the relative weight of the ordered water contribution diminishes more rapidly while the relative weight of the disordered water contribution increases, indicating that a residual water content, characterized by a highly structural disorder, rests entrapped in the montmorillonite layer structure for a longer time. From the present study, it can be inferred that the montmorillonite dehydration process promotes the layer self-assembly.
Collapse
|
27
|
QM Implementation in Drug Design: Does It Really Help? Methods Mol Biol 2020. [PMID: 32016884 DOI: 10.1007/978-1-0716-0282-9_2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/13/2023]
Abstract
Computational chemistry allows one to characterize the structure, dynamics, and energetics of protein-ligand interactions, which makes it a valuable tool in drug discovery in both academic research and pharmaceutical industry. Molecular mechanics (MM)-based approaches are widely utilized to assist the discovery of new drug candidates. However, the complexity of protein-ligand interactions challenges the accuracy and efficiency of the commonly used empirical methods. Aiming to provide better accuracy in the description of protein-ligand interactions, quantum mechanics (QM)-based approaches are becoming increasingly explored. In principle, QM calculation includes all contributions to the energy, accounting for terms usually missing in empirical force fields, and provides a greater degree of transferability. The usefulness of QM in drug design cannot be overemphasized. In this chapter, we present recent developments and applications of fragment-based QM method in studying the protein-ligand and protein-protein interactions. We critically discuss the performance of the fragment-based QM method at different ab initio levels while trying to answer a critical question: do QM-based methods really help in drug design?
Collapse
|
28
|
Xu J, Liu J, Liu J, Hu W, He X, Li J. Phase Transition of Ice at High Pressures and Low Temperatures. Molecules 2020; 25:molecules25030486. [PMID: 31979295 PMCID: PMC7037513 DOI: 10.3390/molecules25030486] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Revised: 01/08/2020] [Accepted: 01/21/2020] [Indexed: 11/21/2022] Open
Abstract
The behavior of ice under extreme conditions undergoes the change of intermolecular binding patterns and leads to the structural phase transitions, which are needed for modeling the convection and internal structure of the giant planets and moons of the solar system as well as H2O-rich exoplanets. Such extreme conditions limit the structural explorations in laboratory but open a door for the theoretical study. The ice phases IX and XIII are located in the high pressure and low temperature region of the phase diagram. However, to the best of our knowledge, the phase transition boundary between these two phases is still not clear. In this work, based on the second-order Møller–Plesset perturbation (MP2) theory, we theoretically investigate the ice phases IX and XIII and predict their structures, vibrational spectra and Gibbs free energies at various extreme conditions, and for the first time confirm that the phase transition from ice IX to XIII can occur around 0.30 GPa and 154 K. The proposed work, taking into account the many-body electrostatic effect and the dispersion interactions from the first principles, opens up the possibility of completing the ice phase diagram and provides an efficient method to explore new phases of molecular crystals.
Collapse
Affiliation(s)
- Jinjin Xu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China;
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Jinfeng Liu
- State Key Laboratory of Natural Medicines, Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China;
| | - Jinyun Liu
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu 241000, China
- Correspondence: (J.L.); (X.H.); (J.L.)
| | - Wenxin Hu
- The Computer Center, School of Computer Science and Software Engineering, East China Normal University, Shanghai 200062, China;
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China;
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Correspondence: (J.L.); (X.H.); (J.L.)
| | - Jinjin Li
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu 241000, China
- Correspondence: (J.L.); (X.H.); (J.L.)
| |
Collapse
|
29
|
Affiliation(s)
- Jongbaik Ree
- Department of Chemistry EducationChonnam National University Gwangju 61186 South Korea
| | - Yoo Hang Kim
- Department of ChemistryInha University Incheon 22212 South Korea
| | - Hyung Kyu Shin
- Department of ChemistryUniversity of Nevada Reno 89557 USA
| |
Collapse
|
30
|
Huang L, Han Y, He X, Li J. Ab initio-enabled phase transition prediction of solid carbon dioxide at ultra-high temperatures. RSC Adv 2020; 10:236-243. [PMID: 35492555 PMCID: PMC9049158 DOI: 10.1039/c9ra06478h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2019] [Accepted: 11/12/2019] [Indexed: 12/27/2022] Open
Abstract
Carbon dioxide is one of the fundamental chemical species on Earth, but its solid-phase behavior at high pressures is still far from well understood and some phases remain uncertain or unknown, which increases the challenge to predict its structures.
Collapse
Affiliation(s)
- Lei Huang
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education
- Department of Micro/Nano-Electronics
- Shanghai Jiao Tong University
- Shanghai 200240
- China
| | - Yanqiang Han
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education
- Department of Micro/Nano-Electronics
- Shanghai Jiao Tong University
- Shanghai 200240
- China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development
- School of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
- China
| | - Jinjin Li
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education
- Department of Micro/Nano-Electronics
- Shanghai Jiao Tong University
- Shanghai 200240
- China
| |
Collapse
|
31
|
Orlando F, Artiglia L, Yang H, Kong X, Roy K, Waldner A, Chen S, Bartels-Rausch T, Ammann M. Disordered Adsorbed Water Layers on TiO 2 Nanoparticles under Subsaturated Humidity Conditions at 235 K. J Phys Chem Lett 2019; 10:7433-7438. [PMID: 31725306 DOI: 10.1021/acs.jpclett.9b02779] [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 interaction of water with TiO2 is of substantial scientific and technological interest as it determines the activity of TiO2 in photocatalytic and environmental applications in nanoparticle suspensions in water, in complex appliances, or in pure form interacting with water vapor. The influence of TiO2 nanoparticles on the hydrogen bonding structure of water molecules is an important factor that controls hydration of other species, reactions, or nucleation processes. We use a combination of ambient-pressure X-ray photoelectron spectroscopy and electron yield near-edge X-ray absorption fine structure (NEXAFS) spectroscopy at the oxygen K-edge to investigate the hydrogen bonding structure of adsorbed water on titania nanoparticles in equilibrium with nearly saturated water vapor at 235 K. The results clearly show that the net NEXAFS spectrum of adsorbed water resembles that of liquid, disordered water at 235 K, a temperature at which both homogeneous and heterogeneous freezing of bulk water is anticipated.
Collapse
Affiliation(s)
- Fabrizio Orlando
- Laboratory of Environmental Chemistry , Paul Scherrer Institut , 5232 Villigen , Switzerland
| | - Luca Artiglia
- Laboratory of Environmental Chemistry , Paul Scherrer Institut , 5232 Villigen , Switzerland
| | - Huanyu Yang
- Laboratory of Environmental Chemistry , Paul Scherrer Institut , 5232 Villigen , Switzerland
- Institute of Atmospheric and Climate Science , ETH Zürich , 8092 Zürich , Switzerland
| | - Xiangrui Kong
- Laboratory of Environmental Chemistry , Paul Scherrer Institut , 5232 Villigen , Switzerland
- Department of Chemistry and Molecular Biology , University of Gothenburg , 41296 Gothenburg , Sweden
| | - Kanak Roy
- Institute for Chemical and Bioengineering , ETH Zürich , 8093 Zürich , Switzerland
| | - Astrid Waldner
- Laboratory of Environmental Chemistry , Paul Scherrer Institut , 5232 Villigen , Switzerland
| | - Shuzhen Chen
- Laboratory of Environmental Chemistry , Paul Scherrer Institut , 5232 Villigen , Switzerland
| | - Thorsten Bartels-Rausch
- Laboratory of Environmental Chemistry , Paul Scherrer Institut , 5232 Villigen , Switzerland
| | - Markus Ammann
- Laboratory of Environmental Chemistry , Paul Scherrer Institut , 5232 Villigen , Switzerland
| |
Collapse
|
32
|
Acharya S, Nandi UK, Bhattacharyya SM. Comparative Study of Anomalous Size Dependence of Charged and Neutral Solute Diffusion in Water. J Phys Chem B 2019; 123:10275-10285. [PMID: 31697084 DOI: 10.1021/acs.jpcb.9b08023] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a comparative study of size dependence of diffusion for charged and neutral solutes in water. Although both show nonmonotonicity of the size dependence of diffusion, their nature and origin are quite different. For neutral solutes, the peak position and the value of diffusion at the maximum are both independent of the solute-water interaction. Interestingly, for charged solutes, with an increase in solute-water interaction strength, the peak position shifts to lower solute sizes and with an increase in charge, it shifts to higher solute sizes. The diffusion value at the peak reduces with an increase in both solute-water interaction and solute charge. We show that all these features observed for charged solutes can be understood in terms of the interplay between ionic and nonionic interactions which is definitely absent for neutral solutes. Some of the earlier studies addressing the nonmonotonicity in diffusion did suggest the interplay between the two interactions to be the cause. However, this is the first time we show that such an interplay gives rise to the nonmonotonicity in the potential energy which is a prerequisite for obtaining the nonmonotonicity in the diffusion. Such nonmonotonicity in the potential energy is absent for neutral solutes.
Collapse
Affiliation(s)
- Sayantan Acharya
- Polymer Science and Engineering Department , National Chemical Laboratory , Pune 411008 , India.,Academy of Scientific and Innovative Research (AcSIR) , Ghaziabad 201002 , India
| | - Ujjwal K Nandi
- Polymer Science and Engineering Department , National Chemical Laboratory , Pune 411008 , India.,Academy of Scientific and Innovative Research (AcSIR) , Ghaziabad 201002 , India
| | - Sarika Maitra Bhattacharyya
- Polymer Science and Engineering Department , National Chemical Laboratory , Pune 411008 , India.,Academy of Scientific and Innovative Research (AcSIR) , Ghaziabad 201002 , India
| |
Collapse
|
33
|
Abstract
Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
34
|
Liu J, Rana B, Liu KY, Herbert JM. Variational Formulation of the Generalized Many-Body Expansion with Self-Consistent Charge Embedding: Simple and Correct Analytic Energy Gradient for Fragment-Based ab Initio Molecular Dynamics. J Phys Chem Lett 2019; 10:3877-3886. [PMID: 31251619 DOI: 10.1021/acs.jpclett.9b01214] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The many-body expansion (MBE) and its extension to overlapping fragments, the generalized (G)MBE, constitute the theoretical basis for most fragment-based approaches for large-scale quantum chemistry. We reformulate the GMBE for use with embedding charges determined self-consistently from the fragment wave functions, in a manner that preserves the variational nature of the underlying self-consistent field method. As a result, the analytic gradient retains the simple "sum of fragment gradients" form that is often assumed in practice, sometimes incorrectly. This obviates (without approximation) the need to solve coupled-perturbed equations, and we demonstrate stable, fragment-based ab initio molecular dynamics simulations using this technique. Energy conservation fails when charge-response contributions to the Fock matrix are neglected, even while geometry optimizations and vibrational frequency calculations may yet be accurate. Stable simulations can be recovered by means of straightforward modifications introduced here, providing a general paradigm for fragment-based ab initio molecular dynamics.
Collapse
Affiliation(s)
- Jie Liu
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| | - Bhaskar Rana
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| | - Kuan-Yu Liu
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| | - John M Herbert
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| |
Collapse
|
35
|
Liu J, Sun H, Glover WJ, He X. Prediction of Excited-State Properties of Oligoacene Crystals Using Fragment-Based Quantum Mechanical Method. J Phys Chem A 2019; 123:5407-5417. [PMID: 31187994 DOI: 10.1021/acs.jpca.8b12552] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
A fundamental understanding of the excited-state properties of molecular crystals is of central importance for their optoelectronics applications. In this study, we developed the electrostatically embedded generalized molecular fractionation (EE-GMF) method for the quantitative characterization of the excited-state properties of locally excited molecular clusters. The accuracy of the EE-GMF method is systematically assessed for oligoacene crystals. Our result demonstrates that the EE-GMF method is capable of providing the lowest vertical singlet (S1) and triplet excitation energies (T1), in excellent agreement with the full-system quantum mechanical calculations. Using this method, we also investigated the performance of different density functionals in predicting the excited-state properties of the oligoacene crystals. Among the 13 tested functionals, B3LYP and MN15 give the two lowest overall mean unsigned errors with reference to the experimental S1 and T1 excitation energies. The EE-GMF approach can be readily utilized for studying the excited-state properties of large-scale organic solids at diverse ab initio levels.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy , China Pharmaceutical University , Nanjing 210009 , China
| | | | - William J Glover
- NYU Shanghai , Shanghai 200122 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062 , China.,Department of Chemistry , New York University , New York , New York 10003 , United States
| | - Xiao He
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062 , China
| |
Collapse
|
36
|
Progress in molecular-simulation-based research on the effects of interface-induced fluid microstructures on flow resistance. Chin J Chem Eng 2019. [DOI: 10.1016/j.cjche.2019.02.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
37
|
Crystal Structure Optimization and Gibbs Free Energy Comparison of Five Sulfathiazole Polymorphs by the Embedded Fragment QM Method at the DFT Level. CRYSTALS 2019. [DOI: 10.3390/cryst9050256] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Molecular crystal plays an important role in many fields of science and technology, but it often crystallizes in different polymorphs with different physical properties. To guide the experimental synthesis of candidate materials, the atomic-scale model is frequently used to predict the most stable polymorph and its structural properties. Here, we show how an ab initio method can be used to achieve a rapid and accurate prediction of sulfathiazole crystal polymorphs (an antibiotic drug), based on the Gibbs free energy calculation and Raman spectra analysis. At the atmospheric pressure and the temperature of 300 K, we demonstrate that form III (FIII) is the most stable structure of sulfathiazole. The agreement between the predicted and experimental crystal structures corresponds to the order of stability for five sulfathiazole polymorphs as FI < FV < FIV < FII < FIII, which is achieved by employing the density functional theory (DFT) calculations.
Collapse
|
38
|
Ansari N, Dandekar R, Caravati S, Sosso GC, Hassanali A. High and low density patches in simulated liquid water. J Chem Phys 2018; 149:204507. [PMID: 30501251 DOI: 10.1063/1.5053559] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
We present insights into the nature of structural heterogeneities in liquid water by characterizing the empty space within the hydrogen bond network. Using molecular dynamics simulations, we show that density fluctuations create regions of empty space characterized by a diverse morphology - from spherical to fractal-like voids. These voids allow for the identification of low and high density patches of the liquid, encompassing short (0.3-0.5 nm) as well as long (1-2 nm) length-scales. In addition, we show that the formation of these patches is coupled to collective fluctuations involving the topology of hydrogen-bonded rings of water molecules. In particular, water molecules in the high density patches tend to be slightly more tetrahedral - which is consistent with the predictions of the hydrophobic effect.
Collapse
Affiliation(s)
- N Ansari
- The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| | - R Dandekar
- The Institute of Mathematical Sciences-HBNI, 4th Cross Street, CIT Campus, Tharamani, Chennai, India
| | - S Caravati
- Department of Chemistry, University of Zurich, Winterhurerstrasse 190, Zurich CH-8057, Switzerland
| | - G C Sosso
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Gibbet Hill, Coventry CV4 7AL, United Kingdom
| | - A Hassanali
- The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
39
|
Liu J, Zhang JZH, He X. Probing the Ion-Specific Effects at the Water/Air Interface and Water-Mediated Ion Pairing in Sodium Halide Solution with Ab Initio Molecular Dynamics. J Phys Chem B 2018; 122:10202-10209. [PMID: 30351119 DOI: 10.1021/acs.jpcb.8b09513] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The ion-specific effects at the water/air interface represent a fundamentally essential topic of research, and high-level ab initio simulations are still demanding to reveal the microscopic picture of the interactions between ions and water at the solvation interface. In this work, we present a fragment-based ab initio molecular dynamics (AIMD) simulation of sodium halide solution droplet (in a neutral mixture of Na+, F-, Cl-, and Br- ions) at the MP2/aug-cc-pVDZ level. We show that the studied halide ions exhibit surface preference in the order (F- < Cl- < Br-) which is in accordance with the experimental observation. The resulting potential of mean force (PMF) for Br- produces a distinct minimum at the water/air interface, while the minimum of the PMF for F- appears in the bulk region. The ion-pairing interactions between halide anions and Na+ cations are characterized, and it reveals that the specific solvent-separated ion pairs (SIPs) are more preferred than the direct contact ion pairs (CIPs). The transition between different types of SIPs is observed. Other structural and dynamical properties of ions and ion-hydration shells are investigated. These results provide broader and new physical insights for understanding the ion-specific behavior in interfacial solvation at the atomistic level.
Collapse
Affiliation(s)
- Jinfeng Liu
- State Key Laboratory of Natural Medicines, Department of Basic Medicine and Clinical Pharmacy , China Pharmaceutical University , Nanjing , 210009 , China.,Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China
| | - John Z H Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai , 200062 , China.,Department of Chemistry , New York University , New York , New York 10003 , United States
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai , 200062 , China
| |
Collapse
|
40
|
Liu J, He X, Zhang JZH, Qi LW. Hydrogen-bond structure dynamics in bulk water: insights from ab initio simulations with coupled cluster theory. Chem Sci 2018; 9:2065-2073. [PMID: 29675248 PMCID: PMC5885775 DOI: 10.1039/c7sc04205a] [Citation(s) in RCA: 68] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2017] [Accepted: 12/04/2017] [Indexed: 02/04/2023] Open
Abstract
An accurate and efficient ab initio molecular dynamics (AIMD) simulation of liquid water was made possible using the fragment-based approach (J. F. Liu, X. He and J. Z. H. Zhang, Phys. Chem. Chem. Phys., 2017, 19, 11931-11936). In this study, we advance the AIMD simulations using the fragment-based coupled cluster (CC) theory, more accurately revealing the structural and dynamical properties of liquid water under ambient conditions. The results show that the double-donor hydrogen-bond configurations in liquid water are nearly in balance with the single-donor configurations, with a slight bias towards the former. Our observation is in contrast to the traditional tetrahedral water structure. The hydrogen-bond switching dynamics in liquid water are very fast, with a hydrogen-bond life time of around 0.78 picoseconds, determined using AIMD simulation at the CCD/aug-cc-pVDZ level. This time scale is remarkably shorter than the ∼3.0 picoseconds that is commonly obtained from traditional nonpolarized force fields and density functional theory (DFT) based first-principles simulations. Additionally, the obtained radial distribution functions, triplet oxygen angular distribution, diffusion coefficient, and the dipole moment of the water molecule are uniformly in good agreement with the experimental observations. The current high-level AIMD simulation sheds light on the understanding of the structural and dynamical properties of liquid water.
Collapse
Affiliation(s)
- Jinfeng Liu
- State Key Laboratory of Natural Medicines , Department of Basic Medicine and Clinical Pharmacy , China Pharmaceutical University , Nanjing , 210009 , China .
| | - Xiao He
- School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China
- NYU-ECNU Center for Computational Chemistry , NYU Shanghai , Shanghai , 200062 , China
| | - John Z H Zhang
- School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China
- NYU-ECNU Center for Computational Chemistry , NYU Shanghai , Shanghai , 200062 , China
- Department of Chemistry , New York University , New York , NY 10003 , USA
| | - Lian-Wen Qi
- State Key Laboratory of Natural Medicines , Department of Basic Medicine and Clinical Pharmacy , China Pharmaceutical University , Nanjing , 210009 , China .
| |
Collapse
|
41
|
Liu J, He X. Accurate prediction of energetic properties of ionic liquid clusters using a fragment-based quantum mechanical method. Phys Chem Chem Phys 2017; 19:20657-20666. [DOI: 10.1039/c7cp03356g] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Accurate prediction of physicochemical properties of ionic liquids (ILs) is of great significance to understand and design novel ILs with unique properties.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy
- China Pharmaceutical University
- Nanjing
- China
| | - Xiao He
- School of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
- China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai
| |
Collapse
|