1
|
Avallone N, Huppert S, Depondt P, Andriambariarijaona L, Datchi F, Ninet S, Plé T, Spezia R, Finocchi F. Orientational Disorder Drives Site Disorder in Plastic Ammonia Hemihydrate. PHYSICAL REVIEW LETTERS 2024; 133:106102. [PMID: 39303235 DOI: 10.1103/physrevlett.133.106102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 07/24/2024] [Indexed: 09/22/2024]
Abstract
In the 2-10 GPa pressure range, ammonia hemihydrate H_{2}O:(NH_{3})_{2} (AHH) is a molecular solid in which intermolecular interactions are ruled by distinct types of hydrogen bonds. Upon heating, the low-temperature ordered P2_{1}/c crystal (AHH-II) transits to a bcc phase (AHH-pbcc) where each site is randomly occupied by water or ammonia. In addition to the site disorder, experiments suggest that AHH-pbcc is a plastic solid, but the physical origin and mechanisms at play for the rotational and site disordering remain unknown. Using large-scale (∼10^{5} atoms) and long-time (>10 ns) simulations, we show that, as temperature rises above the transition line, orientational disorder sets in, breaking the strongest hydrogen bonds that provide the largest contribution to the cohesion of the ordered AHH-II phase and enabling the molecules to migrate from a crystal site to a neighboring one. This generates a plastic molecular alloy with site disorder while the solid state is overall maintained until melting at a higher temperature. The case of high (P,T) plastic ammonia hemihydrate can be extended to other water-ammonia alloys where a similar interplay between distinct hydrogen bonds occurs.
Collapse
|
2
|
Creazzo F, Luber S. Water-air interface revisited by means of path-integral ab initio molecular dynamics. Phys Chem Chem Phys 2024; 26:21290-21302. [PMID: 39078670 PMCID: PMC11305098 DOI: 10.1039/d4cp02500h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2024] [Accepted: 07/21/2024] [Indexed: 07/31/2024]
Abstract
Although nuclear quantum effects (NQEs) have been considered on bulk liquid water, the impact of these latter on the air-water interface has not yet been reported. Herein, by performing and comparing ab initio molecular dynamics (AIMD) and path integral AIMD (PI-AIMD) simulations, we reveal the impact of NQEs on structural, dynamical and electronic properties as well as IR spectra of the air-water interface at room temperature. NQEs, being able to describe a more accurate proton delocalization in H-bonded system than AIMD, reveal a different structural arrangement and dynamical behaviour of both bulk and interfacial water molecules in comparison to AIMD results. A more de-structured and de-bound water arrangement and coordination are identified when the quantum nature of nuclei are considered for both bulk and interfacial water molecules. Structural properties, such as inter-/intra-molecular bond lengths, coordination numbers and H-bonding angles of bulk and interfacial water molecules here calculated, are affected by NQEs mitigating the overstructured description given by AIMD. Further evidences of an AIMD overstructured description of bulk water are in the computed IR spectra, where an increased absorption peak intensity and an increased strength of the hydrogen-bond network are alleviated by NQEs. In addition, NQEs show a valuable impact on the electronic structure of the air-water interface, reducing the total valence bandwidth and the electronic energy band-gap when passing from bulk to interfacial water. This work proves how NQEs significantly affect properties and features of the air-water interface, that are essential to accurately describe H-bonded interfacial systems.
Collapse
Affiliation(s)
- Fabrizio Creazzo
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| | - Sandra Luber
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| |
Collapse
|
3
|
Li M, Wang P, Yu X, Su Y, Zhao J. Impact of Nuclear Quantum Effects on the Structural Properties of Protonated Water Clusters. J Phys Chem A 2024; 128:5954-5962. [PMID: 39007820 DOI: 10.1021/acs.jpca.4c03340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/16/2024]
Abstract
Nuclear quantum effects (NQEs) play a crucial role in hydrogen-bonded systems due to quantum tunneling and proton fluctuation. Our understanding of how NQEs affect microstructures mainly focuses on bulk phases of liquids and solids but remains deficient for water clusters, including their hydrogen nuclei, hydrogen-bonded configurations, and temperature dependence. Here, we conducted ab initio molecular dynamics (MD) and path integral MD simulations to investigate the influence of NQEs on the structural properties of protonated water clusters H+(H2O)n (n = 3, 6, 9, 12). The results reveal that the NQEs become less evident as the cluster size increases due to the competition between NQEs and electrostatic interactions. Simulations of several H+(H2O)6 isomers at different temperatures indicate that the effect of elevated temperature on proton transfer is related to the initial structure. Interestingly, the process of proton transfer also involves the interconversion between Zundel-type and Eigen-type isomers. These findings significantly deepen our understanding of ion-water and water-water interactions, opening new avenues for the study of hydrated ion clusters and related systems.
Collapse
Affiliation(s)
- Mengxu Li
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | | | - Xueke Yu
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | - Yan Su
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | - Jijun Zhao
- Guangdong Basic Research Center of Excellence for Structure and Fundamental Interactions of Matter, Guangdong Provincial Key Laboratory of Quantum Engineering and Quantum Materials, School of Physics, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
4
|
Escamilla-Herrera LF, López-Picón JL, Torres-Arenas J, Gil-Villegas A. Semiclassical thermodynamic geometry. Phys Rev E 2024; 109:064145. [PMID: 39020900 DOI: 10.1103/physreve.109.064145] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Accepted: 05/23/2024] [Indexed: 07/20/2024]
Abstract
In this work the thermodynamic geometry (TG) of semiclassical fluids is analyzed. We present results for two models. The first one is a semiclassical hard-sphere (SCHS) fluid whose Helmholtz free energy is obtained from path-integral Monte Carlo simulations. It is found that, due to quantum contributions in the thermodynamic potential, the anomaly found in TG for the classical hard-sphere fluid related to the sign of the scalar curvature is now avoided in a considerable region of the thermodynamic space. The second model is a semiclassical square-well fluid, described by a SCHS repulsive interaction coupled with a classical attractive square-well contribution. The behavior of the semiclassical curvature scalar as a function of the thermal de Broglie wavelength λ_{B} is analyzed for several attractive-potential ranges. A description of the semiclassical R Widom lines, defined by the maxima of the curvature scalar, is also obtained and results are compared with the corresponding classical systems for different square-well ranges.
Collapse
Affiliation(s)
- L F Escamilla-Herrera
- División de Ciencias e Ingenierías Campus León, Universidad de Guanajuato, AP E-143, CP 37150, León, Guanajuato, México
| | - J L López-Picón
- División de Ciencias e Ingenierías Campus León, Universidad de Guanajuato, AP E-143, CP 37150, León, Guanajuato, México
| | - José Torres-Arenas
- División de Ciencias e Ingenierías Campus León, Universidad de Guanajuato, AP E-143, CP 37150, León, Guanajuato, México
| | - Alejandro Gil-Villegas
- División de Ciencias e Ingenierías Campus León, Universidad de Guanajuato, AP E-143, CP 37150, León, Guanajuato, México
| |
Collapse
|
5
|
Sionkowski P, Kruszewska N, Kreitschitz A, Gorb SN, Domino K. Application of Recurrence Plot Analysis to Examine Dynamics of Biological Molecules on the Example of Aggregation of Seed Mucilage Components. ENTROPY (BASEL, SWITZERLAND) 2024; 26:380. [PMID: 38785629 PMCID: PMC11119629 DOI: 10.3390/e26050380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/25/2024] [Revised: 04/24/2024] [Accepted: 04/26/2024] [Indexed: 05/25/2024]
Abstract
The goal of the research is to describe the aggregation process inside the mucilage produced by plant seeds using molecular dynamics (MD) combined with time series algorithmic analysis based on the recurrence plots. The studied biological molecules model is seed mucilage composed of three main polysaccharides, i.e. pectins, hemicellulose, and cellulose. The modeling of biological molecules is based on the assumption that a classical-quantum passage underlies the aggregation process in the mucilage, resulting from non-covalent interactions, as they affect the macroscopic properties of the system. The applied recurrence plot approach is an important tool for time series analysis and data mining dedicated to analyzing time series data originating from complex, chaotic systems. In the current research, we demonstrated that advanced algorithmic analysis of seed mucilage data can reveal some features of the dynamics of the system, namely temperature-dependent regions with different dynamics of increments of a number of hydrogen bonds and regions of stable oscillation of increments of a number of hydrophobic-polar interactions. Henceforth, we pave the path for automatic data-mining methods for the analysis of biological molecules with the intermediate step of the application of recurrence plot analysis, as the generalization of recurrence plot applications to other (biological molecules) datasets is straightforward.
Collapse
Affiliation(s)
- Piotr Sionkowski
- Institute of Theoretical and Applied Informatics, Polish Academy of Sciences, ul. Bałtycka 5, 44-100 Gliwice, Poland; (P.S.); (K.D.)
| | - Natalia Kruszewska
- Group of Modeling of Physicochemical Processes, Faculty of Chemical Technology and Engineering, Bydgoszcz University of Science and Technology, 85-796 Bydgoszcz, Poland
| | - Agnieszka Kreitschitz
- Department of Plant Developmental Biology, University of Wrocław, ul. Kanonia 6/8, 50-328 Wrocław, Poland;
| | - Stanislav N. Gorb
- Department of Functional Morphology and Biomechanics, Kiel University, Am Botanischen Garten 1-9, D-24098 Kiel, Germany;
| | - Krzysztof Domino
- Institute of Theoretical and Applied Informatics, Polish Academy of Sciences, ul. Bałtycka 5, 44-100 Gliwice, Poland; (P.S.); (K.D.)
| |
Collapse
|
6
|
Zhao L, Xiao PP, Wang Y, Lu Y, Karim TM, Gies H, Yokoi T. Modulation of Al Distribution in High-Silica ZSM-5 Zeolites for Enhancing Catalytic Performance. ACS APPLIED MATERIALS & INTERFACES 2024; 16:17701-17714. [PMID: 38546502 DOI: 10.1021/acsami.4c02311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Abstract
The spatial distribution of framework Al (AlF) has been one of the important factors that affect the catalytic properties of zeolites in diverse chemical reactions; however, the synthesis of high-silica zeolites with special AlF distribution remains a challenge. In this study, we successfully synthesized high-silica ZSM-5 zeolites with a unique AlF distribution by employing pentaerythritol (PET) as an additive in the presence of a few tetrapropylammonium hydroxide (TPAOH). The results demonstrated that the introduction of PET led to a higher proportion of Al atoms located at the sinusoidal and/or straight channels. It was observed that the addition of PET prevented the interaction between TPA+ and tetrahedral [AlO4]- during the crystallization process, resulting in enhanced availability of TPA species in the form of ion-paired TPA+. This effect leads to AlF atoms dominantly distributed away from the intersection and located in narrow channels, where acidic sites more effectively inhibit hydrogen transfer and coke formation. In the reaction of dimethyl ether (DME) to olefins, the catalyst with a unique Al distribution exhibited a significant prolonged catalytic lifetime, surpassing traditional TPA-ZSM-5 by more than 2-fold and maintaining DME conversion above 90% for a maximum of 148 h. The results of multiple pulse experiments also showed that these PET-assisted ZSM-5 zeolites significantly enhanced the selectivity of propene and butene. This approach provides an effective strategy to regulate AlF distribution in high-silica ZSM-5 catalysts with the assistance of neutral alcohol. It holds great potential for application in the synthesis of other high-silica zeolites, thereby enriching the diversity of zeolite catalysis.
Collapse
Affiliation(s)
- Liang Zhao
- Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta, Midori-ku, Yokohama 226-8503, Japan
| | - Pei-Pei Xiao
- Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta, Midori-ku, Yokohama 226-8503, Japan
| | - Yong Wang
- Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta, Midori-ku, Yokohama 226-8503, Japan
- iPEACE223 Inc., Konwa Building, Tsukiji, Chuo-ku, Tokyo 1-12-22, Japan
| | - Yao Lu
- Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta, Midori-ku, Yokohama 226-8503, Japan
| | - Tahta Muslim Karim
- Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta, Midori-ku, Yokohama 226-8503, Japan
| | - Hermann Gies
- Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta, Midori-ku, Yokohama 226-8503, Japan
| | - Toshiyuki Yokoi
- Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta, Midori-ku, Yokohama 226-8503, Japan
- iPEACE223 Inc., Konwa Building, Tsukiji, Chuo-ku, Tokyo 1-12-22, Japan
| |
Collapse
|
7
|
Udagawa T, Tanaka H, Kuwahata K, Tachikawa M. Location of the Shared Proton in Proton-Bound Dimer Compound of Hydrogen Sulfate and Formate: Path Integral Molecular Dynamics Study. J Phys Chem A 2024; 128:2103-2110. [PMID: 38414141 DOI: 10.1021/acs.jpca.4c00239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024]
Abstract
The structure of the proton-bound dimer compound of hydrogen sulfate and formate has been studied by considering nuclear quantum effects (NQEs) using the path integral molecular dynamics method. This study unveiled the location of the shared proton and answered the following question: "Is the shared proton localized on either an anion or located around the center of two anions?" We have elucidated that the shared proton is distributed in the region beyond the transition state due to the NQEs, even though the shared proton did not completely overcome the transition state for the proton shuttle.
Collapse
Affiliation(s)
- Taro Udagawa
- Department of Chemistry and Biomolecular Science, Gifu University, Yanagido 1-1, Gifu 501-1193, Japan
| | - Hikaru Tanaka
- Department of Chemistry and Biomolecular Science, Gifu University, Yanagido 1-1, Gifu 501-1193, Japan
| | - Kazuaki Kuwahata
- Graduate School of NanobioScience, Yokohama City University, 22-2 Seto, Kanazawa-ku, Yokohama 236-0027, Japan
| | - Masanori Tachikawa
- Graduate School of NanobioScience, Yokohama City University, 22-2 Seto, Kanazawa-ku, Yokohama 236-0027, Japan
| |
Collapse
|
8
|
Zhang YY, Zheng H, Wang T, Jiang S, Yan W, Wang C, Zhao Y, Lu JB, Hu HS, Yang J, Zhang W, Wu G, Xie H, Li G, Jiang L, Yang X, Li J. Spectroscopic and Theoretical Identifications of Two Structural Motifs of (H 2O) 10 Cluster. J Phys Chem Lett 2024; 15:3055-3060. [PMID: 38466221 DOI: 10.1021/acs.jpclett.4c00210] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
Precise characterization of archetypal systems of aqueous hydrogen-bonding networks is essential for developing accurate potential functions and universal models of water. The structures of water clusters (H2O)n (n = 2-9) have been verified recently through size-specific infrared spectroscopy with a vacuum ultraviolet free electron laser (VUV-FEL) and quantum chemical studies. For (H2O)10, the pentagonal prism and butterfly motifs were proposed to be important building blocks and were observed in previous experiments. Here we report the size-specific infrared spectra of (H2O)10 via a joint experimental and theoretical study. Well-resolved spectra provide a unique signature for the coexistence of pentagonal prism and butterfly motifs. These (H2O)10 motifs develop from the dominant structures of (H2O)n (n = 8, 9) clusters. This work provides an intriguing prelude to the diverse structure of liquid water and opens avenues for size-dependent measurement of larger systems to understand the stepwise formation mechanism of hydrogen-bonding networks.
Collapse
Affiliation(s)
- Yang-Yang Zhang
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
- Department of Chemistry and Guangdong Provincial Key Laboratory of Catalytic Chemistry, Southern University of Science and Technology, Shenzhen 518055, China
| | - Huijun Zheng
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, 19A Yuquan Road, Beijing 100049, China
| | - Tiantong Wang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, 19A Yuquan Road, Beijing 100049, China
| | - Shuai Jiang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, 19A Yuquan Road, Beijing 100049, China
| | - Wenhui Yan
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, 19A Yuquan Road, Beijing 100049, China
| | - Chong Wang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, 19A Yuquan Road, Beijing 100049, China
| | - Ya Zhao
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Jun-Bo Lu
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
- Department of Chemistry and Guangdong Provincial Key Laboratory of Catalytic Chemistry, Southern University of Science and Technology, Shenzhen 518055, China
| | - Han-Shi Hu
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Jiayue Yang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Weiqing Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Guorong Wu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Hua Xie
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Gang Li
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Ling Jiang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- Hefei National Laboratory, Hefei 230088, China
| | - Xueming Yang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- Department of Chemistry and Guangdong Provincial Key Laboratory of Catalytic Chemistry, Southern University of Science and Technology, Shenzhen 518055, China
- Hefei National Laboratory, Hefei 230088, China
| | - Jun Li
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
- Department of Chemistry and Guangdong Provincial Key Laboratory of Catalytic Chemistry, Southern University of Science and Technology, Shenzhen 518055, China
| |
Collapse
|
9
|
Wang J, Long R. Nuclear Quantum Effects Accelerate Charge Recombination but Boost the Stability of Inorganic Perovskites in Mild Humidity. NANO LETTERS 2024; 24:3476-3483. [PMID: 38445608 DOI: 10.1021/acs.nanolett.4c00245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Experiments have demonstrated that mild humidity can enhance the stability of the CsPbBr3 perovskite, though the underlying mechanism remains unclear. Utilizing ab initio molecular dynamics, ring polymer molecular dynamics, and non-adiabatic molecular dynamics, our study reveals that nuclear quantum effects (NQEs) play a crucial role in stabilizing the lattice rigidity of the perovskite while simultaneously shortening the charge carrier lifetime. NQEs reduce the extent of geometric disorder and the number of atomic fluctuations, diminish the extent of hole localization, and thereby improve the electron-hole overlap and non-adiabatic coupling. Concurrently, these effects significantly suppress phonon modes and slow decoherence. As a result, these factors collectively accelerate charge recombination by a factor of 1.42 compared to that in scenarios excluding NQEs. The resulting sub-10 ns recombination time scales align remarkably well with experimental findings. This research offers novel insight into how moisture resistance impacts the stability and charge carrier lifetime in all-inorganic perovskites.
Collapse
Affiliation(s)
- Jiao Wang
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Run Long
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People's Republic of China
| |
Collapse
|
10
|
Arandhara M, Ramesh SG. Nuclear Quantum Effects in Hydroxide Hydrate Along the H-Bond Bifurcation Pathway. J Phys Chem A 2024; 128:1600-1610. [PMID: 38393819 DOI: 10.1021/acs.jpca.3c08027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2024]
Abstract
Path integral (PI) simulations are used to explore nuclear quantum effects (NQEs) in hydroxide hydrate and its perdeuterated isotopomer along the H-bond bifurcation pathway. Toward this, a new potential energy surface using the symmetric gradient domain machine learning method with ab initio data at the CCSD(T)/aug-cc-pVTZ level is built. From PI umbrella sampling (US) simulations, free energy profiles along the bifurcation coordinate are explored as a function of temperature. At ambient temperature, the bifurcation barrier is increased upon inclusion of NQEs. At low temperatures in the deep tunneling regime, the barrier is strongly decreased and flattened. These trends are examined, and the role of the O-O distance is also investigated through two-dimensional US simulations.
Collapse
Affiliation(s)
- Mrinal Arandhara
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| | - Sai G Ramesh
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
11
|
Byju S, Hassan A, Whitford PC. The energy landscape of the ribosome. Biopolymers 2024; 115:e23570. [PMID: 38051695 DOI: 10.1002/bip.23570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2023] [Revised: 10/17/2023] [Accepted: 11/08/2023] [Indexed: 12/07/2023]
Abstract
The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly.
Collapse
Affiliation(s)
- Sandra Byju
- Center for Theoretical Biological Physics, Northeastern University, Boston, Massachusetts, USA
- Department of Physics, Northeastern University, Boston, Massachusetts, USA
| | - Asem Hassan
- Department of Chemistry, The University of Texas at Austin, Austin, Texas, United States
| | - Paul C Whitford
- Center for Theoretical Biological Physics, Northeastern University, Boston, Massachusetts, USA
- Department of Physics, Northeastern University, Boston, Massachusetts, USA
| |
Collapse
|
12
|
Arandhara M, Ramesh SG. Nuclear quantum effects in gas-phase 2-fluoroethanol. Phys Chem Chem Phys 2024; 26:6885-6902. [PMID: 38333949 DOI: 10.1039/d3cp05657k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Torsional motions along the FCCO and HOCC dihedrals lead to the five unique conformations of 2-fluoroethanol, of which the conformer that is gauche along both dihedrals has the lowest energy. In this work, we explore how nuclear quantum effects (NQEs) manifest in the structural parameters of the lowest energy conformer, in the intramolecular free energy landscape along the FCCO and HOCC dihedrals, and also in the infrared spectrum of the title molecule, through the use of path integral simulations. We have first developed a full dimensional potential energy surface using the reaction surface Hamiltonian framework. On this potential, we have carried out path integral molecular dynamics simulations at several temperatures starting from the minimum energy well to explore structural influences of NQEs including geometrical markers of the interaction between the OH and F groups. From the computed free energy landscapes, significant reduction of the torsional barrier is found at low temperature near the cis region of the dihedrals, which can be understood through the trends in the radii of gyration of the atomic ring polymers. We find that the inclusion of NQEs in the computation of infrared spectrum is important to obtain good agreement with the experimental band positions.
Collapse
Affiliation(s)
- Mrinal Arandhara
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India.
| | - Sai G Ramesh
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India.
| |
Collapse
|
13
|
DeLuca M, Sensale S, Lin PA, Arya G. Prediction and Control in DNA Nanotechnology. ACS APPLIED BIO MATERIALS 2024; 7:626-645. [PMID: 36880799 DOI: 10.1021/acsabm.2c01045] [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] [Indexed: 03/08/2023]
Abstract
DNA nanotechnology is a rapidly developing field that uses DNA as a building material for nanoscale structures. Key to the field's development has been the ability to accurately describe the behavior of DNA nanostructures using simulations and other modeling techniques. In this Review, we present various aspects of prediction and control in DNA nanotechnology, including the various scales of molecular simulation, statistical mechanics, kinetic modeling, continuum mechanics, and other prediction methods. We also address the current uses of artificial intelligence and machine learning in DNA nanotechnology. We discuss how experiments and modeling are synergistically combined to provide control over device behavior, allowing scientists to design molecular structures and dynamic devices with confidence that they will function as intended. Finally, we identify processes and scenarios where DNA nanotechnology lacks sufficient prediction ability and suggest possible solutions to these weak areas.
Collapse
Affiliation(s)
- Marcello DeLuca
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Sebastian Sensale
- Department of Physics, Cleveland State University, Cleveland, Ohio 44115, United States
| | - Po-An Lin
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Gaurav Arya
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
14
|
Dufils T, Schran C, Chen J, Geim AK, Fumagalli L, Michaelides A. Origin of dielectric polarization suppression in confined water from first principles. Chem Sci 2024; 15:516-527. [PMID: 38179530 PMCID: PMC10763014 DOI: 10.1039/d3sc04740g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Accepted: 11/23/2023] [Indexed: 01/06/2024] Open
Abstract
It has long been known that the dielectric constant of confined water should be different from that in bulk. Recent experiments have shown that it is vanishingly small, however the origin of the phenomenon remains unclear. Here we used ab initio molecular dynamics simulations (AIMD) and AIMD-trained machine-learning potentials to understand water's structure and electronic properties underpinning this effect. For the graphene and hexagonal boron-nitride substrates considered, we find that it originates in the spontaneous anti-parallel alignment of the water dipoles in the first two water layers near the solid interface. The interfacial layers exhibit net ferroelectric ordering, resulting in an overall anti-ferroelectric arrangement of confined water. Together with constrained hydrogen-bonding orientations, this leads to much reduced out-of-plane polarization. Furthermore, we directly contrast AIMD and simple classical force-field simulations, revealing important differences. This work offers insight into a property of water that is critical in modulating surface forces, the electric-double-layer formation and molecular solvation, and shows a way to compute it.
Collapse
Affiliation(s)
- T Dufils
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - C Schran
- Cavendish Laboratory, Department of Physics, University of Cambridge Cambridge CB3 0HE UK
- Lennard-Jones Centre, University of Cambridge Trinity Ln Cambridge CB2 1TN UK
| | - J Chen
- School of Physics, Peking University Beijing 100871 China
| | - A K Geim
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - L Fumagalli
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - A Michaelides
- Lennard-Jones Centre, University of Cambridge Trinity Ln Cambridge CB2 1TN UK
- Yusuf Hamied Department of Chemistry, University of Cambridge Lensfield Road Cambridge CB2 1EW UK
| |
Collapse
|
15
|
Georgiev DD. Evolution of Consciousness. Life (Basel) 2023; 14:48. [PMID: 38255663 PMCID: PMC10817314 DOI: 10.3390/life14010048] [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: 09/06/2023] [Revised: 12/01/2023] [Accepted: 12/25/2023] [Indexed: 01/24/2024] Open
Abstract
The natural evolution of consciousness in different animal species mandates that conscious experiences are causally potent in order to confer any advantage in the struggle for survival. Any endeavor to construct a physical theory of consciousness based on emergence within the framework of classical physics, however, leads to causally impotent conscious experiences in direct contradiction to evolutionary theory since epiphenomenal consciousness cannot evolve through natural selection. Here, we review recent theoretical advances in describing sentience and free will as fundamental aspects of reality granted by quantum physical laws. Modern quantum information theory considers quantum states as a physical resource that endows quantum systems with the capacity to perform physical tasks that are classically impossible. Reductive identification of conscious experiences with the quantum information comprised in quantum brain states allows for causally potent consciousness that is capable of performing genuine choices for future courses of physical action. The consequent evolution of brain cortical networks contributes to increased computational power, memory capacity, and cognitive intelligence of the living organisms.
Collapse
Affiliation(s)
- Danko D Georgiev
- Institute for Advanced Study, 30 Vasilaki Papadopulu Str., 9010 Varna, Bulgaria
| |
Collapse
|
16
|
Plé T, Lagardère L, Piquemal JP. Force-field-enhanced neural network interactions: from local equivariant embedding to atom-in-molecule properties and long-range effects. Chem Sci 2023; 14:12554-12569. [PMID: 38020379 PMCID: PMC10646944 DOI: 10.1039/d3sc02581k] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2023] [Accepted: 10/03/2023] [Indexed: 12/01/2023] Open
Abstract
We introduce FENNIX (Force-Field-Enhanced Neural Network InteraXions), a hybrid approach between machine-learning and force-fields. We leverage state-of-the-art equivariant neural networks to predict local energy contributions and multiple atom-in-molecule properties that are then used as geometry-dependent parameters for physically-motivated energy terms which account for long-range electrostatics and dispersion. Using high-accuracy ab initio data (small organic molecules/dimers), we trained a first version of the model. Exhibiting accurate gas-phase energy predictions, FENNIX is transferable to the condensed phase. It is able to produce stable Molecular Dynamics simulations, including nuclear quantum effects, for water predicting accurate liquid properties. The extrapolating power of the hybrid physically-driven machine learning FENNIX approach is exemplified by computing: (i) the solvated alanine dipeptide free energy landscape; (ii) the reactive dissociation of small molecules.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, LCT, UMR 7616 CNRS F-75005 Paris France thomas.ple@sorbonne-université louis.lagardere@sorbonne-université jean-philip.piquemal@sorbonne-université
| | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS F-75005 Paris France thomas.ple@sorbonne-université louis.lagardere@sorbonne-université jean-philip.piquemal@sorbonne-université
| | - Jean-Philip Piquemal
- Sorbonne Université, LCT, UMR 7616 CNRS F-75005 Paris France thomas.ple@sorbonne-université louis.lagardere@sorbonne-université jean-philip.piquemal@sorbonne-université
| |
Collapse
|
17
|
Fung YKC, Perkin S. Structure and anomalous underscreening in ethylammonium nitrate solutions confined between two mica surfaces. Faraday Discuss 2023; 246:370-386. [PMID: 37458200 PMCID: PMC10568257 DOI: 10.1039/d3fd00042g] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 03/02/2023] [Indexed: 10/13/2023]
Abstract
The observation of long-range interactions across ionic liquids and highly concentrated electrolytes, extending far beyond the Debye-Hückel prediction and beyond the range predicted in liquid state theory, has been called 'anomalous underscreening'. A number of theoretical and experimental works have explored this phenomenon over recent years, although its origin is not yet fully understood. Most of the experimental studies of anomalous underscreening until now involved aprotic ionic liquids, and so it is of interest to explore interactions in protic ionic liquids where the distribution of charge in the fluid is different in nature. Here we present direct measurements of the interaction force as a function of separation distance, measured using a surface force balance, across solutions of a protic ionic liquid ethylammonium nitrate (EAN) and its mixtures with water over a range of volume fractions from 10 vol% to 100 vol% EAN. The results reveal intricate details about near-surface ordering and dynamics at the EAN-mica interface as well as anomalous underscreening consistent with that observed in the past with aprotic ionic liquids.
Collapse
Affiliation(s)
- Y K Catherine Fung
- Physical & Theoretical Chemistry Laboratory, University of Oxford, Oxford, UK.
| | - Susan Perkin
- Physical & Theoretical Chemistry Laboratory, University of Oxford, Oxford, UK.
| |
Collapse
|
18
|
Dvores MP, Çarçabal P, Gerber RB. Selective reactivity of glycosyl cation stereoisomers: the role of intramolecular hydrogen bonding. Phys Chem Chem Phys 2023; 25:26737-26747. [PMID: 37779496 DOI: 10.1039/d3cp03326k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/03/2023]
Abstract
The impact of the stereochemistry of the glycosyl cation species upon its dynamic properties is examined together with their vibrational spectra in order to gain insight into the effects of configurational isomerism on conformer dynamics and proton mobility. Ab initio molecular dynamics (AIMD) simulations and infrared multiple photon dissociation (IRMPD) spectroscopy explore the conformational and reactive dynamics of two pairs of glycosyl cation isomers: (1) protonated α- and β- anomers of methyl-D-galactopyranoside and (2) the oxocarbenium ions of the D-aldohexose C2 epimers galactose and talose. Analysis of these simulations together with experimental spectroscopy, interpreted by anharmonic calculations, points to the key role played by the intramolecular hydrogen bonds which are present in a unique pattern and extent in each isomer. We find that the reactivity of galactoside stereoisomers toward acid-catalyzed nucleophilic substitution, as gauged by the ability to form free oxocarbenium ions, differs markedly in a way that agrees with experimental measurements in the condensed phase. Other properties such as conformer stability and vibrational transitions were also found to reflect the characteristic hydrogen bonding interactions present in each isomer. In both systems, the stereochemistry is shown to determine the strength of intramolecular hydrogen bonding as well as between which substituents proton transfer is possible. We expect that the critical impact of non-covalent interactions on stereoisomer selectivity may be a widely found phenomenon whose effects should be further investigated.
Collapse
Affiliation(s)
- M P Dvores
- Fritz Haber Centre for Molecular Dynamics, The Institute of Chemistry, The Hebrew University, Jerusalem 91904, Israel.
| | - P Çarçabal
- Insitut des Sciences Moléculaires d'Orsay, ISMO, Univ Paris-Sud, CNRS, bat 210, Univ Paris-Sud, 91405 Orsay Cedex, France
| | - R B Gerber
- Fritz Haber Centre for Molecular Dynamics, The Institute of Chemistry, The Hebrew University, Jerusalem 91904, Israel.
- Department of Chemistry, University of California Irvine, CA, 92697, USA
| |
Collapse
|
19
|
Jiang Y, Hu Z, Zhong C, Yang Y, Wang XB, Sun Z, Sun H, Liu Z, Peng P. Locking water molecules via ternary O-H⋯O intramolecular hydrogen bonds in perhydroxylated closo-dodecaborate. Phys Chem Chem Phys 2023; 25:25810-25817. [PMID: 37724455 DOI: 10.1039/d3cp03555g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/20/2023]
Abstract
A multitude of applications related to perhydroxylated closo-dodecaborate B12(OH)122- in the condensed phase are inseparable from the fundamental mechanisms underlying the high water orientation selectivity based on the base B12(OH)122-. Herein, we directly compare the structural evolution of water clusters, ranging from monomer to hexamer, oriented by functional groups in the bases B12H122-, B12H11OH2- and B12(OH)122- using multiple theoretical methods. A significant revelation is made regarding B12(OH)122-: each additional water molecule is locked into the intramolecular hydrogen bond B-O-H ternary ring in an embedded form. This new pattern of water cluster growth suggests that B-(H-O)⋯H-O interactions prevail over the competition from water-hydrogen bonds (O⋯H-O), distinguishing it from the behavior observed in B12H122- and B12H11OH2- bases, in which competition arises from a mixed competing model involving dihydrogen bonds (B-H⋯H-O), conventional hydrogen bonds (B-(H-O)⋯H-O) and water hydrogen bonds (O⋯H-O). Through aqueous solvation and ab initio molecular dynamics analysis, we further demonstrate the largest water clusters in the first hydrated shell with exceptional thermodynamic stability around B12(OH)122-. These findings provide a solid scientific foundation for the design of boron cluster chemistry incorporating hydroxyl-group-modified borate salts with potential implications for various applications.
Collapse
Affiliation(s)
- Yanrong Jiang
- Center for Transformative Science, ShanghaiTech University, Shanghai 201210, China.
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China.
| | - Zhubin Hu
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China.
| | - Cheng Zhong
- College of Chemistry & Molecular Sciences, Wuhan University, Wuhan, Hubei 430072, China
| | - Yan Yang
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China.
| | - Xue-Bin Wang
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| | - Zhenrong Sun
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China.
| | - Haitao Sun
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China.
| | - Zhi Liu
- Center for Transformative Science, ShanghaiTech University, Shanghai 201210, China.
| | - Peng Peng
- Center for Transformative Science, ShanghaiTech University, Shanghai 201210, China.
| |
Collapse
|
20
|
Schaack S, Mangaud E, Fallacara E, Huppert S, Depondt P, Finocchi F. When Quantum Fluctuations Meet Structural Instabilities: The Isotope- and Pressure-Induced Phase Transition in the Quantum Paraelectric NaOH. PHYSICAL REVIEW LETTERS 2023; 131:126101. [PMID: 37802932 DOI: 10.1103/physrevlett.131.126101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Revised: 05/17/2023] [Accepted: 08/15/2023] [Indexed: 10/08/2023]
Abstract
Anhydrous sodium hydroxide, a common and structurally simple compound, shows spectacular isotope effects: NaOD undergoes a first-order transition, which is absent in NaOH. By combining ab initio electronic structure calculations with Feynman path integrals, we show that NaOH is an unusual example of a quantum paraelectric: zero-point quantum fluctuations stretch the weak hydrogen bonds (HBs) into a region where they are unstable and break. By strengthening the HBs via isotope substitution or applied pressure, the system can be driven to a broken-symmetry antiferroelectric phase. In passing, we provide a simple quantitative criterion for HB breaking in layered crystals and show that nuclear quantum effects are crucial in paraelectric to ferroelectric transitions in hydrogen-bonded hydroxides.
Collapse
Affiliation(s)
- Sofiane Schaack
- Sorbonne Université, CNRS UMR 7588, Institut des NanoSciences de Paris, INSP, 75005 Paris, France
| | - Etienne Mangaud
- Sorbonne Université, CNRS UMR 7588, Institut des NanoSciences de Paris, INSP, 75005 Paris, France
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Erika Fallacara
- Sorbonne Université, CNRS UMR 7588, Institut des NanoSciences de Paris, INSP, 75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, CNRS UMR 7588, Institut des NanoSciences de Paris, INSP, 75005 Paris, France
| | - Philippe Depondt
- Sorbonne Université, CNRS UMR 7588, Institut des NanoSciences de Paris, INSP, 75005 Paris, France
| | - Fabio Finocchi
- Sorbonne Université, CNRS UMR 7588, Institut des NanoSciences de Paris, INSP, 75005 Paris, France
| |
Collapse
|
21
|
Zhou S, Liu Y, Wang S, Wang L. Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds. Sci Rep 2023; 13:13741. [PMID: 37612311 PMCID: PMC10447522 DOI: 10.1038/s41598-023-40614-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2023] [Accepted: 08/14/2023] [Indexed: 08/25/2023] Open
Abstract
There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.
Collapse
Affiliation(s)
| | - Yuanhao Liu
- Department of Statistics, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA
| | - Sijian Wang
- Department of Statistics, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA.
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA.
| |
Collapse
|
22
|
Liu R, Chen M. Characterization of the Hydrogen-Bond Network in High-Pressure Water by Deep Potential Molecular Dynamics. J Chem Theory Comput 2023; 19:5602-5608. [PMID: 37535904 DOI: 10.1021/acs.jctc.3c00445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/05/2023]
Abstract
The hydrogen-bond (H-bond) network of high-pressure water is investigated by neural-network-based molecular dynamics (MD) simulations with first-principles accuracy. The static structure factors (SSFs) of water at three densities, i.e., 1, 1.115, and 1.24 g/cm3, are directly evaluated from 512 water MD trajectories, which are in quantitative agreement with the experiments. We propose a new method to decompose the computed SSF and identify the changes in the SSF with respect to the changes in H-bond structures. We find that a larger water density results in a higher probability for one or two non-H-bonded water molecules to be inserted into the inner shell, explaining the changes in the tetrahedrality of water under pressure. We predict that the structure of the accepting end of water molecules is more easily influenced by the pressure than by the donating end. Our work sheds new light on explaining the SSF and H-bond properties in related fields.
Collapse
Affiliation(s)
- Renxi Liu
- HEDPS, CAPT, College of Engineering, Peking University, Beijing 100871, P. R. China
- Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 90871, P. R. China
| | - Mohan Chen
- HEDPS, CAPT, College of Engineering, Peking University, Beijing 100871, P. R. China
- Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 90871, P. R. China
- AI for Science Institute, Beijing 100080, P. R. China
| |
Collapse
|
23
|
Atsango AO, Morawietz T, Marsalek O, Markland TE. Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations. J Chem Phys 2023; 159:074101. [PMID: 37581418 DOI: 10.1063/5.0162066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Accepted: 07/31/2023] [Indexed: 08/16/2023] Open
Abstract
The transport of excess protons and hydroxide ions in water underlies numerous important chemical and biological processes. Accurately simulating the associated transport mechanisms ideally requires utilizing ab initio molecular dynamics simulations to model the bond breaking and formation involved in proton transfer and path-integral simulations to model the nuclear quantum effects relevant to light hydrogen atoms. These requirements result in a prohibitive computational cost, especially at the time and length scales needed to converge proton transport properties. Here, we present machine-learned potentials (MLPs) that can model both excess protons and hydroxide ions at the generalized gradient approximation and hybrid density functional theory levels of accuracy and use them to perform multiple nanoseconds of both classical and path-integral proton defect simulations at a fraction of the cost of the corresponding ab initio simulations. We show that the MLPs are able to reproduce ab initio trends and converge properties such as the diffusion coefficients of both excess protons and hydroxide ions. We use our multi-nanosecond simulations, which allow us to monitor large numbers of proton transfer events, to analyze the role of hypercoordination in the transport mechanism of the hydroxide ion and provide further evidence for the asymmetry in diffusion between excess protons and hydroxide ions.
Collapse
Affiliation(s)
- Austin O Atsango
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Tobias Morawietz
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Ondrej Marsalek
- Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
24
|
Aarabi M, Sarka J, Pandey A, Nieman R, Aquino AJA, Eckert J, Poirier B. Quantum Dynamical Investigation of Dihydrogen-Hydride Exchange in a Transition-Metal Polyhydride Complex. J Phys Chem A 2023; 127:6385-6399. [PMID: 37494557 DOI: 10.1021/acs.jpca.3c01863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
The ongoing shift toward clean, sustainable energy is a primary driving force behind hydrogen fuel research. Safe and effective storage of hydrogen is a major challenge (particularly for mobile applications) and requires a detailed understanding of the atomic level interactions of hydrogen with its host materials. The light mass of hydrogen, however, implies that quantum effects are important, so a quantum dynamical treatment is required to properly account for these effects in computational simulations. As one such example, we describe herein the hydrogen exchange dynamics between a hydride and a dihydrogen ligand in the [FeH(H2)(PH3)4]+ model complex. A global three-dimensional (3D) potential energy surface (PES) was constructed by fitting to and interpolating from a discrete set of grid points computed using density functional theory; exact quantum dynamical calculations were then carried out on the 3D PES using discrete variable representation basis sets. Energy levels and their quantum tunneling splittings were computed up to 3000 cm-1 above the ground state. Within that energy range, all three fundamentals have been identified using wave function plots, as well as the first three overtones of the exchange (reaction coordinate) motion and several of its combination bands. From the tunneling splittings, the Boltzmann-averaged tunneling rates were computed. The Arrhenius plot of the total exchange rate shows a clear transition around 150 K, below which the activation energy is essentially zero and above which it is less than half of the electronic structure barrier. This indicates that exchange rates are governed by quantum tunneling throughout the relevant temperature range with the low-temperature regime dominated by a single quantum (ground) state. This work is the first-ever fully quantum dynamical study to investigate the hydrogen exchange dynamics between hydride and dihydrogen ligands coordinated to a transition-metal complex.
Collapse
Affiliation(s)
- Mahdi Aarabi
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
| | - János Sarka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
- Institute of Chemistry, Eötvös Loránd University, H-1117 Budapest, Hungary
| | - Ankit Pandey
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
| | - Reed Nieman
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
| | - Adelia J A Aquino
- Department of Mechanical Engineering, Texas Tech University, Lubbock, Texas 79409-1061, United States
| | - Juergen Eckert
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
| |
Collapse
|
25
|
Bhatt H, Verma AK, Modak P. Polyether phases of formic acid revealed under high pressure. Chem Commun (Camb) 2023; 59:9888-9891. [PMID: 37493352 DOI: 10.1039/d3cc02129g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2023]
Abstract
Evolutionary structure searches reveal that at 20 GPa, formic acid (FA) transforms to an orthorhombic ether-chain polymer phase, ruling out hydrogen bond symmetrization, followed by a novel crowned cyclic-ether tetragonal phase above 60 GPa. Emergence of characteristic polyether and new OH stretching modes in infrared experiments validate the findings. Resemblance of polymer chain with a cosmopolymer polyoxymethylene, shows an engrossing multifaceted evolution of FA.
Collapse
Affiliation(s)
- Himal Bhatt
- High Pressure & Synchrotron Radiation Physics Division, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India.
- Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India
| | - Ashok K Verma
- High Pressure & Synchrotron Radiation Physics Division, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India.
- Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India
| | - P Modak
- High Pressure & Synchrotron Radiation Physics Division, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India.
- Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India
| |
Collapse
|
26
|
Xu X. Constrained Nuclear-Electronic Orbital Density Functional Theory with a Dielectric Continuum Solvent Model. J Phys Chem A 2023. [PMID: 37470267 DOI: 10.1021/acs.jpca.3c02507] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/21/2023]
Abstract
Solvent effects are crucial for simulating chemical and biological processes in solutions. The continuum solvation model is widely used for incorporating solvent effects with different levels of theoretical descriptions of solutes. For solutes and solutions containing hydrogen atoms, nuclear quantum effects can also be nonnegligible for reliable simulations. In this work, we couple our recently developed constrained nuclear-electronic orbital density functional theory with a dielectric continuum solvation model to cover nuclear quantum effects and solvent effects simultaneously. This approach is applied to the formate ion, where an anomalous solvatochromic shift in C-H stretch frequency was reported in experiments. By using this new approach to account for nuclear quantum effects and solvent effects, we show that the vibrational frequency of the C-H stretch and the solvatochromic shift are accurately described.
Collapse
Affiliation(s)
- Xi Xu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| |
Collapse
|
27
|
Chow M, Lambros E, Li X, Hammes-Schiffer S. Nuclear-Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics. J Chem Theory Comput 2023. [PMID: 37329317 DOI: 10.1021/acs.jctc.3c00361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear-electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. These developments and initial examples provide the foundation for future studies of nuclear-electronic quantum dynamics in complex chemical and biological environments.
Collapse
Affiliation(s)
- Mathew Chow
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Eleftherios Lambros
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | | |
Collapse
|
28
|
Zhu YC, Yang S, Zeng JX, Fang W, Jiang L, Zhang DH, Li XZ. Accurate calculation of tunneling splittings in water clusters using path-integral based methods. J Chem Phys 2023; 158:2895223. [PMID: 37290067 DOI: 10.1063/5.0146562] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Accepted: 04/03/2023] [Indexed: 06/10/2023] Open
Abstract
Tunneling splittings observed in molecular rovibrational spectra are significant evidence for tunneling motion of hydrogen nuclei in water clusters. Accurate calculations of the splitting sizes from first principles require a combination of high-quality inter-atomic interactions and rigorous methods to treat the nuclei with quantum mechanics. Many theoretical efforts have been made in recent decades. This Perspective focuses on two path-integral based tunneling splitting methods whose computational cost scales well with the system size, namely, the ring-polymer instanton method and the path-integral molecular dynamics (PIMD) method. From a simple derivation, we show that the former is a semiclassical approximation to the latter, despite that the two methods are derived very differently. Currently, the PIMD method is considered to be an ideal route to rigorously compute the ground-state tunneling splitting, while the instanton method sacrifices some accuracy for a significantly smaller computational cost. An application scenario of such a quantitatively rigorous calculation is to test and calibrate the potential energy surfaces of molecular systems by spectroscopic accuracy. Recent progress in water clusters is reviewed, and the current challenges are discussed.
Collapse
Affiliation(s)
- Yu-Cheng Zhu
- State Key Laboratory for Artificial Microstructure and Mesoscopic Physics, Frontier Science Center for Nano-optoelectronics and School of Physics, Peking University, Beijing 100871, People's Republic of China
- Interdisciplinary Institute of Light-Element Quantum Materials, Research Center for Light-Element Advanced Materials, and Collaborative Innovation Center of Quantum Materials, Peking University, Beijing 100871, People's Republic of China
| | - Shuo Yang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
- University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
| | - Jia-Xi Zeng
- State Key Laboratory for Artificial Microstructure and Mesoscopic Physics, Frontier Science Center for Nano-optoelectronics and School of Physics, Peking University, Beijing 100871, People's Republic of China
- Interdisciplinary Institute of Light-Element Quantum Materials, Research Center for Light-Element Advanced Materials, and Collaborative Innovation Center of Quantum Materials, Peking University, Beijing 100871, People's Republic of China
| | - Wei Fang
- Department of Chemistry, Fudan University, Shanghai 200438, People's Republic of China
| | - Ling Jiang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Xin-Zheng Li
- State Key Laboratory for Artificial Microstructure and Mesoscopic Physics, Frontier Science Center for Nano-optoelectronics and School of Physics, Peking University, Beijing 100871, People's Republic of China
- Interdisciplinary Institute of Light-Element Quantum Materials, Research Center for Light-Element Advanced Materials, and Collaborative Innovation Center of Quantum Materials, Peking University, Beijing 100871, People's Republic of China
- Peking University Yangtze Delta Institute of Optoelectronics, Nantong, Jiangsu 226010, People's Republic of China
| |
Collapse
|
29
|
Bore SL, Paesani F. Realistic phase diagram of water from "first principles" data-driven quantum simulations. Nat Commun 2023; 14:3349. [PMID: 37291095 PMCID: PMC10250386 DOI: 10.1038/s41467-023-38855-1] [Citation(s) in RCA: 26] [Impact Index Per Article: 26.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2023] [Accepted: 05/12/2023] [Indexed: 06/10/2023] Open
Abstract
Since the experimental characterization of the low-pressure region of water's phase diagram in the early 1900s, scientists have been on a quest to understand the thermodynamic stability of ice polymorphs on the molecular level. In this study, we demonstrate that combining the MB-pol data-driven many-body potential for water, which was rigorously derived from "first principles" and exhibits chemical accuracy, with advanced enhanced-sampling algorithms, which correctly describe the quantum nature of molecular motion and thermodynamic equilibria, enables computer simulations of water's phase diagram with an unprecedented level of realism. Besides providing fundamental insights into how enthalpic, entropic, and nuclear quantum effects shape the free-energy landscape of water, we demonstrate that recent progress in "first principles" data-driven simulations, which rigorously encode many-body molecular interactions, has opened the door to realistic computational studies of complex molecular systems, bridging the gap between experiments and simulations.
Collapse
Affiliation(s)
- Sigbjørn Løland Bore
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA, 92093, USA.
- Materials Science and Engineering, University of California San Diego, La Jolla, CA, 92093, USA.
- Halicioğlu Data Science Institute, University of California San Diego, La Jolla, CA, 92093, USA.
- San Diego Supercomputer Center, University of California San Diego, La Jolla, CA, 92093, USA.
| |
Collapse
|
30
|
Plé T, Mauger N, Adjoua O, Inizan TJ, Lagardère L, Huppert S, Piquemal JP. Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: From Force Fields to Machine Learning Potentials. J Chem Theory Comput 2023; 19:1432-1445. [PMID: 36856658 DOI: 10.1021/acs.jctc.2c01233] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2023]
Abstract
We report the implementation of a multi-CPU and multi-GPU massively parallel platform dedicated to the explicit inclusion of nuclear quantum effects (NQEs) in the Tinker-HP molecular dynamics (MD) package. The platform, denoted Quantum-HP, exploits two simulation strategies: the Ring-Polymer Molecular Dynamics (RPMD) that provides exact structural properties at the cost of a MD simulation in an extended space of multiple replicas and the adaptive Quantum Thermal Bath (adQTB) that imposes the quantum distribution of energy on a classical system via a generalized Langevin thermostat and provides computationally affordable and accurate (though approximate) NQEs. We discuss some implementation details, efficient numerical schemes, and parallelization strategies and quickly review the GPU acceleration of our code. Our implementation allows an efficient inclusion of NQEs in MD simulations for very large systems, as demonstrated by scaling tests on water boxes with more than 200,000 atoms (simulated using the AMOEBA polarizable force field). We test the compatibility of the approach with Tinker-HP's recently introduced Deep-HP machine learning potentials module by computing water properties using the DeePMD potential with adQTB thermostatting. Finally, we show that the platform is also compatible with the alchemical free energy estimation capabilities of Tinker-HP and fast enough to perform simulations. Therefore, we study how NQEs affect the hydration free energy of small molecules solvated with the recently developed Q-AMOEBA water force field. Overall, the Quantum-HP platform allows users to perform routine quantum MD simulations of large condensed-phase systems and will help to shed new light on the quantum nature of important interactions in biological matter.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Nastasia Mauger
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | | | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Simon Huppert
- Institut des Nanosciences de Paris (INSP), CNRS UMR 7588, and Sorbonne Université, F-75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France.,Institut Universitaire de France, 75005 Paris, France.,Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
31
|
Méndez E, Videla PE, Laria D. Collective Proton Transfers in Cyclic Water-Ammonia Tetramers: A Path Integral Machine-Learning Study. J Phys Chem A 2023; 127:1839-1848. [PMID: 36794937 DOI: 10.1021/acs.jpca.2c07994] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
We present results from machine-learning-based path integral molecular dynamics simulations that describe isomerization paths articulated via collective proton transfers along mixed, cyclic tetramers combining water and ammonia at cryogenic conditions. The net result of such isomerizations is a reverse of the chirality of the global hydrogen-bonding architecture along the different cyclic moieties. In monocomponent tetramers, the classical free energy profiles associated with these isomerizations present the usual symmetric double-well characteristics whereas the reactive paths exhibit full concertedness among the different intermolecular transfer processes. Contrastingly, in mixed water/ammonia tetramers, the incorporation of a second component introduces imbalances in the strengths of the different hydrogen bonds leading to a partial loss of concertedness, most notably at the vicinity of the transition state. As such, the highest and lowest degrees of progression are registered along OH···N and O···HN coordinations, respectively. These characteristics lead to polarized transition state scenarios akin to solvent-separated ion-pair configurations. The explicit incorporation of nuclear quantum effects promotes drastic depletions in the activation free energies and modifications in the overall shape of the profiles which include central plateau-like stages, indicating the prevalence of deep tunneling regimes. On the other hand, the quantum treatment of the nuclei partially restores the degree of concertedness among the evolutions of the individual transfers.
Collapse
Affiliation(s)
- Emilio Méndez
- Departamento de Química Inorgánica, Analítica y Química-Física and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E Videla
- Department of Chemistry and Energy Sciences Institute, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Daniel Laria
- Departamento de Química Inorgánica, Analítica y Química-Física and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina.,Departamento de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
32
|
Bocus M, Goeminne R, Lamaire A, Cools-Ceuppens M, Verstraelen T, Van Speybroeck V. Nuclear quantum effects on zeolite proton hopping kinetics explored with machine learning potentials and path integral molecular dynamics. Nat Commun 2023; 14:1008. [PMID: 36823162 PMCID: PMC9950054 DOI: 10.1038/s41467-023-36666-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Accepted: 02/10/2023] [Indexed: 02/25/2023] Open
Abstract
Proton hopping is a key reactive process within zeolite catalysis. However, the accurate determination of its kinetics poses major challenges both for theoreticians and experimentalists. Nuclear quantum effects (NQEs) are known to influence the structure and dynamics of protons, but their rigorous inclusion through the path integral molecular dynamics (PIMD) formalism was so far beyond reach for zeolite catalyzed processes due to the excessive computational cost of evaluating all forces and energies at the Density Functional Theory (DFT) level. Herein, we overcome this limitation by training first a reactive machine learning potential (MLP) that can reproduce with high fidelity the DFT potential energy surface of proton hopping around the first Al coordination sphere in the H-CHA zeolite. The MLP offers an immense computational speedup, enabling us to derive accurate reaction kinetics beyond standard transition state theory for the proton hopping reaction. Overall, more than 0.6 μs of simulation time was needed, which is far beyond reach of any standard DFT approach. NQEs are found to significantly impact the proton hopping kinetics up to ~473 K. Moreover, PIMD simulations with deuterium can be performed without any additional training to compute kinetic isotope effects over a broad range of temperatures.
Collapse
Affiliation(s)
- Massimo Bocus
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Ruben Goeminne
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Aran Lamaire
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Maarten Cools-Ceuppens
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Toon Verstraelen
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | | |
Collapse
|
33
|
Kurapothula PJ, Shepherd S, Wilkins DM. Competing Nuclear Quantum Effects and Hydrogen-Bond Jumps in Hydrated Kaolinite. J Phys Chem Lett 2023; 14:1542-1547. [PMID: 36745462 PMCID: PMC9940297 DOI: 10.1021/acs.jpclett.2c03896] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 02/01/2023] [Indexed: 06/18/2023]
Abstract
Recent work has shown that the dynamics of hydrogen bonds in pure clays are affected by nuclear quantum fluctuations, with different effects for the hydrogen bonds holding different layers of the clay together and for those within the same layer. At the clay-water interface there is an even wider range of types of hydrogen bond, suggesting that the quantum effects may be yet more varied. We apply classical and thermostated ring polymer molecular dynamics simulations to show that nuclear quantum effects accelerate hydrogen-bond dynamics to varying degrees. By interpreting the results in terms of the extended jump model of hydrogen-bond switching, we can understand the origins of these effects in terms of changes in the quantum kinetic energy of hydrogen atoms during an exchange. We also show that the extended jump mechanism is applicable not only to the hydrogen bonds involving water, but also those internal to the clay.
Collapse
|
34
|
Liu J, Liu R, Cao Y, Chen M. Solvation structures of calcium and magnesium ions in water with the presence of hydroxide: a study by deep potential molecular dynamics. Phys Chem Chem Phys 2023; 25:983-993. [PMID: 36519362 DOI: 10.1039/d2cp04105g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
The solvation structures of calcium (Ca2+) and magnesium (Mg2+) ions with the presence of hydroxide (OH-) ion in water are essential for understanding their roles in biological and chemical processes but have not been fully explored. Ab initio molecular dynamics (AIMD) is an important tool to address this issue, but two challenges exist. First, an accurate description of OH- from AIMD needs an appropriate exchange-correlation functional. Second, a long trajectory is needed to reach an equilibrium state for the Ca2+-OH- and Mg2+-OH- ion pairs in aqueous solutions. Herein, we adopt a deep potential molecular dynamics (DPMD) method to simulate 1 ns trajectories for the Ca2+-OH- and Mg2+-OH- ion pairs in water; the DPMD method provides efficient machine-learning-based models that have the accuracy of the SCAN exchange-correlation functional within the framework of density functional theory. The solvation structures of the cations and the OH- in terms of three different species have been systematically investigated. On the one hand, we find that OH- have more significant effects on the solvation structure of Ca2+ than that of Mg2+. We observe that the OH- substantially affects the orientation angles of water molecules surrounding the cation. Through the time correlation functions, we conclude that the water molecules in the first solvation shell of Ca2+ change their preferred orientation faster than those of Mg2+. On the other hand, with the presence of the cation in the first solvation shell of OH-, we find that the hydrogen bonds of OH- are severely altered, and the adjacent water molecules of OH- are squeezed. The two cations have substantially different effects on the solvation structure of OH-. Our work provides new insight into the solvation structures of Ca2+ and Mg2+ in water with the presence of OH-.
Collapse
Affiliation(s)
- Jianchuan Liu
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing, 100871, China.
| | - Renxi Liu
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing, 100871, China. .,Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China
| | - Yu Cao
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing, 100871, China.
| | - Mohan Chen
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing, 100871, China. .,Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China
| |
Collapse
|
35
|
Patkar D, Bharati Ahirwar M, Deshmukh MM. A Tug of War between the Self- and Cross-associating Hydrogen Bonds in Neutral Ammonia-Water Clusters: Energetic Insights by Molecular Tailoring Approach. Chemphyschem 2022; 23:e202200476. [PMID: 36127809 DOI: 10.1002/cphc.202200476] [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: 07/05/2022] [Revised: 09/09/2022] [Indexed: 01/05/2023]
Abstract
In the present work, the energies of various types of individual HBs observed in neutral (NH3 )m (H2 O)n , (m+n=2 to 7) clusters were estimated using the molecular tailoring approach (MTA)-based method. The calculated individual HB energies suggest that the O-H…N HBs are the strongest (1.21 to 12.49 kcal mol-1 ). The next ones are the O-H…O (3.97 to 9.30 kcal mol-1 ) HBs. The strengths of N-H…N (1.09 to 5.29 kcal mol-1 ) and N-H…O (2.85 to 5.56 kcal mol-1 ) HBs are the weakest. The HB energies in dimers also follow this rank ordering. However, the HB energies in dimers are much smaller than those obtained by the MTA-based method due to the loss in cooperativity contribution in the dimers. Thus, the calculated cooperativity contributions, for different types of HBs, fall in the range 0.64 to 5.73 kcal mol-1 . We wish to emphasize based on the energetic rank ordering obtained by the MTA-based method that the O-H of water is a better HB donor than the N-H of ammonia. The reasons for the observed energetic rank ordering are two folds: (i) intrinsically stronger O-H…N HBs than the O-H…O ones as revealed by dimer energies and (ii) the higher cooperativity contribution in the former than the later ones. Indeed, the MTA-based method is useful in providing the missing energetic rank ordering of various type of HBs in neutral (NH3 )m (H2 O)n clusters, in the literature.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| |
Collapse
|
36
|
Complexes of carbon dioxide with methanol and its monohalogen-substituted: Beyond the tetrel bond. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.140158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
37
|
Mauger N, Plé T, Lagardère L, Huppert S, Piquemal JP. Improving Condensed-Phase Water Dynamics with Explicit Nuclear Quantum Effects: The Polarizable Q-AMOEBA Force Field. J Phys Chem B 2022; 126:8813-8826. [PMID: 36270033 DOI: 10.1021/acs.jpcb.2c04454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We introduce a new parametrization of the AMOEBA polarizable force field for water denoted Q-AMOEBA, for use in simulations that explicitly account for nuclear quantum effects (NQEs). This study is made possible thanks to the recently introduced adaptive Quantum Thermal Bath (adQTB) simulation technique which computational cost is comparable to classical molecular dynamics. The flexible Q-AMOEBA model conserves the initial AMOEBA functional form, with an intermolecular potential including an atomic multipole description of electrostatic interactions (up to quadrupole), a polarization contribution based on the Thole interaction model and a buffered 14-7 potential to model van der Waals interactions. It has been obtained by using a ForceBalance fitting strategy including high-level quantum chemistry reference energies and selected condensed-phase properties targets. The final Q-AMOEBA model is shown to accurately reproduce both gas-phase and condensed-phase properties, notably improving the original AMOEBA water model. This development allows the fine study of NQEs on water liquid phase properties such as the average H-O-H angle compared to its gas-phase equilibrium value, isotope effects, and so on. Q-AMOEBA also provides improved infrared spectroscopy prediction capabilities compared to AMOEBA03. Overall, we show that the impact of NQEs depends on the underlying model functional form and on the associated strength of hydrogen bonds. Since adQTB simulations can be performed at near classical computational cost using the Tinker-HP package, Q-AMOEBA can be extended to organic molecules, proteins, and nucleic acids opening the possibility for the large-scale study of the importance of NQEs in biophysics.
Collapse
Affiliation(s)
- Nastasia Mauger
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Thomas Plé
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| |
Collapse
|
38
|
Dupertuis N, Tarun OB, Lütgebaucks C, Roke S. Three-Dimensional Confinement of Water: H 2O Exhibits Long-Range (>50 nm) Structure while D 2O Does Not. NANO LETTERS 2022; 22:7394-7400. [PMID: 36067223 DOI: 10.1021/acs.nanolett.2c02206] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Water is the liquid of life thanks to its three-dimensional adaptive hydrogen (H)-bond network. Confinement of this network may lead to dramatic structural changes influencing chemical and physical transformations. Although confinement effects occur on a <1 nm length scale, the upper length scale limit is unknown. Here, we investigate 3D-confinement over lengths scales ranging from 58-140 nm. By confining water in zwitterionic liposomes of different sizes and measuring the change in H-bond network conformation using second harmonic scattering (SHS), we determined long-range confinement effects in light and heavy water. D2O displays no detectable 3D-confinement effects <58 nm (<3 × 106 D2O molecules). H2O is distinctly different. The vesicle enclosed inner H-bond network has a different conformation compared to the outside network and the SHS response scales with the volume of the confining space. H2O displays confinement effects over distances >100 nm (>2 × 107 H2O molecules).
Collapse
Affiliation(s)
- Nathan Dupertuis
- Laboratory for Fundamental BioPhotonics (LBP), Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Orly B Tarun
- Laboratory for Fundamental BioPhotonics (LBP), Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Cornelis Lütgebaucks
- Laboratory for Fundamental BioPhotonics (LBP), Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Sylvie Roke
- Laboratory for Fundamental BioPhotonics (LBP), Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
- Institute of Materials Science and Engineering (IMX), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
- Lausanne Centre for Ultrafast Science, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
39
|
Jiang Y, Yuan Q, Cao W, Hu Z, Yang Y, Zhong C, Yang T, Sun H, Wang XB, Sun Z. Unraveling hydridic-to-protonic dihydrogen bond predominance in monohydrated dodecaborate clusters. Chem Sci 2022; 13:9855-9860. [PMID: 36128244 PMCID: PMC9430482 DOI: 10.1039/d2sc03986a] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2022] [Accepted: 08/02/2022] [Indexed: 12/01/2022] Open
Abstract
Hydridic-to-protonic dihydrogen bonds (DHBs) are involved in comprehensive structural and energetic evolution, and significantly affect reactivity and selectivity in solution and solid states. Grand challenges exist in understanding DHBs' bonding nature and strength, and how to harness DHBs. Herein we launched a combined photoelectron spectroscopy and multiscale theoretical investigation using monohydrated closo-dodecaborate clusters B12X12 2-·H2O (X = H, F, I) to address such challenges. For the first time, a consistent and unambiguous picture is unraveled demonstrating that B-H⋯H-O DHBs are superior to the conventional B-X⋯H-O HBs, being 1.15 and 4.61 kcal mol-1 stronger than those with X = F and I, respectively. Energy decomposition analyses reveal that induction and dispersion terms make pronounced contributions resulting in a stronger B-H⋯H-O DHB. These findings call out more attention to the prominent roles of DHBs in water environments and pave the way for efficient and eco-friendly catalytic dihydrogen production based on optimized hydridic-to-protonic interactions.
Collapse
Affiliation(s)
- Yanrong Jiang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University Shanghai 200241 China
| | - Qinqin Yuan
- Physical Sciences Division, Pacific Northwest National Laboratory 902 Battelle Boulevard, P. O. Box 999, MS K8-88 Richland Washington 99352 USA
- Department of Chemistry, Anhui University Hefei Anhui 230601 China
| | - Wenjin Cao
- Physical Sciences Division, Pacific Northwest National Laboratory 902 Battelle Boulevard, P. O. Box 999, MS K8-88 Richland Washington 99352 USA
| | - Zhubin Hu
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University Shanghai 200241 China
| | - Yan Yang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University Shanghai 200241 China
| | - Cheng Zhong
- College of Chemistry & Molecular Sciences, Wuhan University Wuhan Hubei 430072 China
| | - Tao Yang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University Shanghai 200241 China
| | - Haitao Sun
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University Shanghai 200241 China
- Collaborative Innovation Center of Extreme Optics, Shanxi University Taiyuan Shanxi 030006 China
| | - Xue-Bin Wang
- Physical Sciences Division, Pacific Northwest National Laboratory 902 Battelle Boulevard, P. O. Box 999, MS K8-88 Richland Washington 99352 USA
| | - Zhenrong Sun
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University Shanghai 200241 China
- Collaborative Innovation Center of Extreme Optics, Shanxi University Taiyuan Shanxi 030006 China
| |
Collapse
|
40
|
Non-Covalent Interaction on the Self-Healing of Mechanical Properties in Supramolecular Polymers. Int J Mol Sci 2022; 23:ijms23136902. [PMID: 35805906 PMCID: PMC9266855 DOI: 10.3390/ijms23136902] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2022] [Revised: 06/10/2022] [Accepted: 06/14/2022] [Indexed: 01/27/2023] Open
Abstract
Supramolecular polymers are widely utilized and applied in self-assembly or self-healing materials, which can be repaired when damaged. Normally, the healing process is classified into two types, including extrinsic and intrinsic self-healable materials. Therefore, the aim of this work is to review the intrinsic self-healing strategy based on supramolecular interaction or non-covalent interaction and molecular recognition to obtain the improvement of mechanical properties. In this review, we introduce the main background of non-covalent interaction, which consists of the metal–ligand coordination, hydrogen bonding, π–π interaction, electrostatic interaction, dipole–dipole interaction, and host–guest interactions, respectively. From the perspective of mechanical properties, these interactions act as transient crosslinking points to both prevent and repair the broken polymer chains. For material utilization in terms of self-healing products, this knowledge can be applied and developed to increase the lifetime of the products, causing rapid healing and reducing accidents and maintenance costs. Therefore, the self-healing materials using supramolecular polymers or non-covalent interaction provides a novel strategy to enhance the mechanical properties of materials causing the extended cycling lifetime of products before replacement with a new one.
Collapse
|
41
|
Liu R, Zhang C, Liang X, Liu J, Wu X, Chen M. Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling. J Chem Phys 2022; 157:024503. [DOI: 10.1063/5.0094944] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water has been a challenging task from ab initio molecular dynamics (AIMD). The difficulty mainly comes from a lack of accurate and efficient exchange-correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly-constrained and appropriately normed (SCAN) meta-GGA functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange-correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice in describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. Our AIMD simulations indicate hydroxide diffuses slower than hydronium in water, which is consistent with the experiments.
Collapse
Affiliation(s)
| | | | | | | | - Xifan Wu
- Physics, Temple University, United States of America
| | - Mohan Chen
- College of Engineering, Peking University, China
| |
Collapse
|
42
|
Chu W, Tan S, Zheng Q, Fang W, Feng Y, Prezhdo OV, Wang B, Li XZ, Zhao J. Ultrafast charge transfer coupled to quantum proton motion at molecule/metal oxide interface. SCIENCE ADVANCES 2022; 8:eabo2675. [PMID: 35714193 DOI: 10.1126/sciadv.abo2675] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Understanding how the nuclear quantum effects (NQEs) in the hydrogen bond (H-bond) network influence the photoexcited charge transfer at semiconductor/molecule interface is a challenging problem. By combining two kinds of emerging molecular dynamics methods at the ab initio level, the path integral-based molecular dynamics and time-dependent nonadiabatic molecular dynamics, and choosing CH3OH/TiO2 as a prototypical system to study, we find that the quantum proton motion in the H-bond network is strongly coupled with the ultrafast photoexcited charge dynamics at the interface. The hole trapping ability of the adsorbed methanol molecule is notably enhanced by the NQEs, and thus, it behaves as a hole scavenger on titanium dioxide. The critical role of the H-bond network is confirmed by in situ scanning tunneling microscope measurements with ultraviolet light illumination. It is concluded the quantum proton motion in the H-bond network plays a critical role in influencing the energy conversion efficiency based on photoexcitation.
Collapse
Affiliation(s)
- Weibin Chu
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, CA 90089, USA
- Key Laboratory of Computational Physical Sciences (Ministry of Education), Institute of Computational Physical Sciences, Fudan University, Shanghai 200433, People's Republic of China
| | - Shijing Tan
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Qijing Zheng
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Wei Fang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
- Department of Chemistry, Fudan University, Shanghai 200438, People's Republic of China
- Laboratory of Physical Chemistry, ETH Zurich, CH-8093 Zürich, Switzerland
| | - Yexin Feng
- School of Physics and Electronics, Hunan University, Changsha 410082, People's Republic of China
| | - Oleg V Prezhdo
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, CA 90089, USA
| | - Bing Wang
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Xin-Zheng Li
- Interdisciplinary Institute of Light-Element Quantum Materials, Research Center for Light-Element Advanced Materials, State Key Laboratory for Artificial Microstructure and Mesoscopic Physics, Frontier Science Center for Nano-optoelectronics and School of Physics, Peking University, Beijing 100871, People's Republic of China
- Peking University Yangtze Delta Institute of Optoelectronics, Nantong, Jiangsu 226010, People's Republic of China
| | - Jin Zhao
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
- Hefei National Laboratory, Hefei 230088, People's Republic of China
| |
Collapse
|
43
|
Guo J, Jiang Y. Submolecular Insights into Interfacial Water by Hydrogen-Sensitive Scanning Probe Microscopy. Acc Chem Res 2022; 55:1680-1692. [PMID: 35678704 DOI: 10.1021/acs.accounts.2c00111] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
ConspectusWater-solid interfaces have attracted extensive attention because of their crucial roles in a wide range of chemical and physical processes, such as ice nucleation and growth, dissolution, corrosion, heterogeneous catalysis, and electrochemistry. To understand these processes, enormous efforts have been made to obtain a molecular-level understanding of the structure and dynamics of water on various solid surfaces. By the use of scanning probe microscopy (SPM), many remarkable structures of H-bonding networks have been directly visualized, significantly advancing our understanding of the delicate competition between water-water and water-solid interactions. Moreover, the detailed dynamics of water molecules, such as diffusion, clustering, dissociation, and intermolecular and intramolecular proton transfer, have been investigated in a well-controlled manner by tip manipulation. However, resolving the submolecular structure of surface water has remained a great challenge for a long time because of the small size and light mass of protons. Discerning the position of hydrogen in water is not only crucial for the accurate determination of the structure of H-bonding networks but also indispensable in probing the proton transfer dynamics and the quantum nature of protons.In this Account, we focus on the recent advances in the H-sensitive SPM technique and its applications in probing the structures, dynamics, and nuclear quantum effects (NQEs) of surface water and ion hydrates at the submolecular level. First, we introduce the development of high-resolution scanning tunneling microscopy/spectroscopy (STM/S) and qPlus-based atomic force microscopy (qPlus-AFM), which allow access to the degrees of freedom of protons in both real and energy space. qPlus-AFM even allows imaging of interfacial water in a weakly perturbative manner by measuring the high-order electrostatic force between the CO-terminated tip and the polar water molecule, which enables the subtle difference of OH directionality to be discerned. Next we showcase the applications of H-sensitive STM/AFM in addressing several key issues related to water-solid interfaces. The surface wetting behavior and the H-bonding structure of low-dimensional ice on various hydrophilic and hydrophobic solid surfaces are characterized at the atomic scale. Then we discuss the quantitative assessment of NQEs of surface water, including proton tunneling and quantum delocalization. Moreover, the weakly perturbative and H-sensitive SPM technique can be also extended to investigations of water-ion interactions on solid surfaces, revealing the effect of hydration structure on the interfacial ion transport. Finally, we provide an outlook on the further directions and challenges for SPM studies of water-solid interfaces.
Collapse
Affiliation(s)
- Jing Guo
- College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Ying Jiang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China.,Collaborative Innovation Center of Quantum Matter, Beijing 100871, People's Republic of China.,Interdisciplinary Institute of Light-Element Quantum Materials and Research Center for Light-Element Advanced Materials, Peking University, Beijing 100871, People's Republic of China.,CAS Center for Excellence in Topological Quantum Computation, University of Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| |
Collapse
|
44
|
Patkar D, Ahirwar MB, Deshmukh MM. Energetic Ordering of Hydrogen Bond Strengths in Methanol-Water Clusters: Insights via Molecular Tailoring Approach. Chemphyschem 2022; 23:e202200143. [PMID: 35302702 DOI: 10.1002/cphc.202200143] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Revised: 03/17/2022] [Indexed: 12/31/2022]
Abstract
In this work, we examine the strength of various types of individual hydrogen bond (HB) in mixed methanol-water Mn Wm , (n+m=2 to 7) clusters, with an aim to understand the relative order of their strength, using our recently proposed molecular tailoring-based approach (MTA). Among all the types of HB, it is observed that the OM -H…OW HBs are the strongest (6.9 to 12.4 kcal mol-1 ). The next ones are OM -H…OM HBs (6.5 to 11.6 kcal mol-1 ). The OW -H…OW (0.2 to 10.9 kcal mol-1 ) and OW -H…OM HBs (0.3 to 10.3 kcal mol-1 ) are the weakest ones. This energetic ordering of HBs is seen to be different from the respective HB energies in the dimer i. e., OM -H…OM (5.0 to 6.0 kcal mol-1 )>OW -H…OM (1.5 to 6.0 kcal mol-1 )>OM -H…OW (3.8 to 5.6 kcal mol-1 )>OW -H…OW (1.2 to 5.0 kcal mol-1 ). The plausible reason for the difference in the HB energy ordering may be attributed to the increase or decrease in HB strengths due to the formation of cooperative or anti-cooperative HB networks. For instance, the cooperativity contribution towards the different types of HB follows: OM -H…OW (2.4 to 8.6 kcal mol-1 )>OM -H…OM (1.3 to 6.3 kcal mol-1 )>OW -H…OW (-1.0 to 6.5 kcal mol-1 )>OW -H…OM (-1.2 to 5.3 kcal mol-1 ). This ordering of cooperativity contribution is similar to the HB energy ordering obtained by the MTA-based method. It is emphasized here that, the interplay between the cooperative and anti-cooperative contributions are indispensable for the correct energetic ordering of these HBs.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, A Central University), Sagar, 470003, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, A Central University), Sagar, 470003, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, A Central University), Sagar, 470003, India
| |
Collapse
|
45
|
Dahl PJ, Yi SM, Gu Y, Acharya A, Shipps C, Neu J, O’Brien JP, Morzan UN, Chaudhuri S, Guberman-Pfeffer MJ, Vu D, Yalcin SE, Batista VS, Malvankar NS. A 300-fold conductivity increase in microbial cytochrome nanowires due to temperature-induced restructuring of hydrogen bonding networks. SCIENCE ADVANCES 2022; 8:eabm7193. [PMID: 35544567 PMCID: PMC9094664 DOI: 10.1126/sciadv.abm7193] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Accepted: 03/28/2022] [Indexed: 06/10/2023]
Abstract
Although proteins are considered as nonconductors that transfer electrons only up to 1 to 2 nanometers via tunneling, Geobacter sulfurreducens transports respiratory electrons over micrometers, to insoluble acceptors or syntrophic partner cells, via nanowires composed of polymerized cytochrome OmcS. However, the mechanism enabling this long-range conduction is unclear. Here, we demonstrate that individual nanowires exhibit theoretically predicted hopping conductance, at rate (>1010 s-1) comparable to synthetic molecular wires, with negligible carrier loss over micrometers. Unexpectedly, nanowires show a 300-fold increase in their intrinsic conductance upon cooling, which vanishes upon deuteration. Computations show that cooling causes a massive rearrangement of hydrogen bonding networks in nanowires. Cooling makes hemes more planar, as revealed by Raman spectroscopy and simulations, and lowers their reduction potential. We find that the protein surrounding the hemes acts as a temperature-sensitive switch that controls charge transport by sensing environmental perturbations. Rational engineering of heme environments could enable systematic tuning of extracellular respiration.
Collapse
Affiliation(s)
- Peter J. Dahl
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| | - Sophia M. Yi
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| | - Yangqi Gu
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
- Department of Molecular, Cellular and Developmental Biology, Yale University, New Haven, CT, USA
| | - Atanu Acharya
- Department of Chemistry, Yale University, New Haven, CT, USA
| | - Catharine Shipps
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| | - Jens Neu
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| | - J. Patrick O’Brien
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| | - Uriel N. Morzan
- Department of Chemistry, Yale University, New Haven, CT, USA
| | | | - Matthew J. Guberman-Pfeffer
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| | - Dennis Vu
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| | - Sibel Ebru Yalcin
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| | | | - Nikhil S. Malvankar
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
- Microbial Sciences Institute, Yale University, New Haven, CT, USA
| |
Collapse
|
46
|
Simulation of Nuclear Quantum Effects in Condensed Matter Systems via Quantum Baths. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12094756] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/10/2022]
Abstract
This paper reviews methods that aim at simulating nuclear quantum effects (NQEs) using generalized thermal baths. Generalized (or quantum) baths simulate statistical quantum features, and in particular zero-point energy effects, through non-Markovian stochastic dynamics. They make use of generalized Langevin Equations (GLEs), in which the quantum Bose–Einstein energy distribution is enforced by tuning the random and friction forces, while the system degrees of freedom remain classical. Although these baths have been formally justified only for harmonic oscillators, they perform well for several systems, while keeping the cost of the simulations comparable to the classical ones. We review the formal properties and main characteristics of classical and quantum GLEs, in relation with the fluctuation–dissipation theorems. Then, we describe the quantum thermostat and quantum thermal bath, the two generalized baths currently most used, providing several examples of applications for condensed matter systems, including the calculation of vibrational spectra. The most important drawback of these methods, zero-point energy leakage, is discussed in detail with the help of model systems, and a recently proposed scheme to monitor and mitigate or eliminate it—the adaptive quantum thermal bath—is summarised. This approach considerably extends the domain of application of generalized baths, leading, for instance, to the successful simulation of liquid water, where a subtle interplay of NQEs is at play. The paper concludes by overviewing further development opportunities and open challenges of generalized baths.
Collapse
|
47
|
Štoček JR, Socha O, Císařová I, Slanina T, Dračínský M. Importance of Nuclear Quantum Effects for Molecular Cocrystals with Short Hydrogen Bonds. J Am Chem Soc 2022; 144:7111-7116. [PMID: 35394771 DOI: 10.1021/jacs.1c10885] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Many efforts have been recently devoted to the design and investigation of multicomponent pharmaceutical solids, such as salts and cocrystals. The experimental distinction between these solid forms is often challenging. Here, we show that the transformation of a salt into a cocrystal with a short hydrogen bond does not occur as a sharp phase transition but rather a smooth shift of the positional probability of the hydrogen atoms. A combination of solid-state NMR spectroscopy, X-ray diffraction, and diffuse reflectance measurements with density functional theory calculations that include nuclear quantum effects (NQEs) provides evidence of temperature-induced hydrogen atom shift in cocrystals with short hydrogen bonds. We demonstrate that for the predictions of the salt/cocrystal solid forms with short H-bonds, the computations have to include NQEs (particularly hydrogen nuclei delocalization) and temperature effects.
Collapse
Affiliation(s)
- Jakub Radek Štoček
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 160 00, Czech Republic.,Department of Organic Chemistry, Faculty of Science, Charles University in Prague, Hlavova 2030, Prague 2 12840, Czech Republic
| | - Ondřej Socha
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 160 00, Czech Republic
| | - Ivana Císařová
- Department of Inorganic Chemistry, Faculty of Science, Charles University in Prague, Hlavova 2030, Prague 2 12840, Czech Republic
| | - Tomáš Slanina
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 160 00, Czech Republic
| | - Martin Dračínský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 160 00, Czech Republic
| |
Collapse
|
48
|
Zhang M, Li Z, Jia Y, Wang F, Tian J, Zhang C, Han T, Xing R, Ye W, Wang C. Observing Mesoscopic Nucleic Acid Capacitance Effect and Mismatch Impact via Graphene Transistors. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2105890. [PMID: 35072345 DOI: 10.1002/smll.202105890] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/26/2021] [Revised: 12/05/2021] [Indexed: 06/14/2023]
Abstract
This work reports a molecular-scale capacitance effect of the double helical nucleic acid duplex structure for the first time. By quantitatively conducting large sample measurements of the electrostatic field effect using a type of high-accuracy graphene transistor biosensor, an unusual charge-transport behavior is observed in which the end-immobilized nucleic acid duplexes can store a part of ionization electrons like molecular capacitors, other than electric conductors. To elucidate this discovery, a cascaded capacitive network model is proposed as a novel equivalent circuit of nucleic acid duplexes, expanding the point-charge approximation model, by which the partial charge-transport observation is reasonably attributed to an electron-redistribution behavior within the capacitive network. Furthermore, it is experimentally confirmed that base-pair mismatches hinder the charge transport in double helical duplexes, and lead to directly identifiable alterations in electrostatic field effects. The bioelectronic principle of mismatch impact is also self-consistently explained by the newly proposed capacitive network model. The mesoscopic nucleic acid capacitance effect may enable a new kind of label-free nucleic acid analysis tool based on electronic transistor devices. The in situ and real-time nucleic acid detections for virus biomarkers, somatic mutations, and genome editing off-target may thus be predictable.
Collapse
Affiliation(s)
- Mingfeng Zhang
- Tianjin Key Laboratory of Wireless Mobile Communications and Power Transmission, College of Electronic and Communication Engineering, Tianjin Normal University, Tianjin, 300387, China
| | - Zhibo Li
- Tianjin Key Laboratory of Wireless Mobile Communications and Power Transmission, College of Electronic and Communication Engineering, Tianjin Normal University, Tianjin, 300387, China
| | - Yuan Jia
- Industrialization Center of Micro & Nano ICs and Devices Sino-German College of Intelligent Manufacturing, Shenzhen Technology University, Shenzhen, 518118, China
| | - Fuquan Wang
- Tianjin Key Laboratory of Wireless Mobile Communications and Power Transmission, College of Electronic and Communication Engineering, Tianjin Normal University, Tianjin, 300387, China
| | - Jinpeng Tian
- Industrialization Center of Micro & Nano ICs and Devices Sino-German College of Intelligent Manufacturing, Shenzhen Technology University, Shenzhen, 518118, China
| | - Cuiping Zhang
- Tianjin Key Laboratory of Wireless Mobile Communications and Power Transmission, College of Electronic and Communication Engineering, Tianjin Normal University, Tianjin, 300387, China
| | - Tingting Han
- Tianjin Key Laboratory of Wireless Mobile Communications and Power Transmission, College of Electronic and Communication Engineering, Tianjin Normal University, Tianjin, 300387, China
- Department of Intelligence Science and Technology, College of Artificial Intelligence, Tianjin Normal University, Tianjin, 300387, China
| | - Ruiqing Xing
- Tianjin Key Laboratory of Wireless Mobile Communications and Power Transmission, College of Electronic and Communication Engineering, Tianjin Normal University, Tianjin, 300387, China
- Department of Intelligence Science and Technology, College of Artificial Intelligence, Tianjin Normal University, Tianjin, 300387, China
| | - Weixiang Ye
- Department of Physics, School of Science, Hainan University, Haikou, 570228, China
- Key Laboratory of Engineering Modeling and Statistical Computation of Hainan Province, School of Science, Hainan University, Haikou, 570228, China
| | - Cheng Wang
- Tianjin Key Laboratory of Wireless Mobile Communications and Power Transmission, College of Electronic and Communication Engineering, Tianjin Normal University, Tianjin, 300387, China
- Department of Intelligence Science and Technology, College of Artificial Intelligence, Tianjin Normal University, Tianjin, 300387, China
| |
Collapse
|
49
|
Competitive nuclear quantum effect and H/D isotope effect on torsional motion of H2O2: An ab initio path integral molecular dynamics study. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2021.113542] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
50
|
Zhou S, Liu Y, Wang S, Wang L. Effective prediction of short hydrogen bonds in proteins via machine learning method. Sci Rep 2022; 12:469. [PMID: 35013487 PMCID: PMC8748993 DOI: 10.1038/s41598-021-04306-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Accepted: 12/20/2021] [Indexed: 12/11/2022] Open
Abstract
Short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms lie within 2.7 Å, exhibit prominent quantum mechanical characters and are connected to a wide range of essential biomolecular processes. However, exact determination of the geometry and functional roles of SHBs requires a protein to be at atomic resolution. In this work, we analyze 1260 high-resolution peptide and protein structures from the Protein Data Bank and develop a boosting based machine learning model to predict the formation of SHBs between amino acids. This model, which we name as machine learning assisted prediction of short hydrogen bonds (MAPSHB), takes into account 21 structural, chemical and sequence features and their interaction effects and effectively categorizes each hydrogen bond in a protein to a short or normal hydrogen bond. The MAPSHB model reveals that the type of the donor amino acid plays a major role in determining the class of a hydrogen bond and that the side chain Tyr-Asp pair demonstrates a significant probability of forming a SHB. Combining electronic structure calculations and energy decomposition analysis, we elucidate how the interplay of competing intermolecular interactions stabilizes the Tyr-Asp SHBs more than other commonly observed combinations of amino acid side chains. The MAPSHB model, which is freely available on our web server, allows one to accurately and efficiently predict the presence of SHBs given a protein structure with moderate or low resolution and will facilitate the experimental and computational refinement of protein structures.
Collapse
Affiliation(s)
- Shengmin Zhou
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA
| | - Yuanhao Liu
- Department of Statistics, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA
| | - Sijian Wang
- Department of Statistics, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA.
| |
Collapse
|