1
|
Lv G, Tu Y, Zhang JH, Chen G. Photomolecular effect: Visible light interaction with air-water interface. Proc Natl Acad Sci U S A 2024; 121:e2320844121. [PMID: 38652751 PMCID: PMC11067046 DOI: 10.1073/pnas.2320844121] [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: 11/27/2023] [Accepted: 03/08/2024] [Indexed: 04/25/2024] Open
Abstract
Although water is almost transparent to visible light, we demonstrate that the air-water interface interacts strongly with visible light via what we hypothesize as the photomolecular effect. In this effect, transverse-magnetic polarized photons cleave off water clusters from the air-water interface. We use 14 different experiments to demonstrate the existence of this effect and its dependence on the wavelength, incident angle, and polarization of visible light. We further demonstrate that visible light heats up thin fogs, suggesting that this process can impact weather, climate, and the earth's water cycle and that it provides a mechanism to resolve the long-standing puzzle of larger measured clouds absorption to solar radiation than theory could predict based on bulk water optical constants. Our study suggests that the photomolecular effect should happen widely in nature, from clouds to fogs, ocean to soil surfaces, and plant transpiration and can also lead to applications in energy and clean water.
Collapse
Affiliation(s)
- Guangxin Lv
- Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA02139
| | - Yaodong Tu
- Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA02139
| | - James H. Zhang
- Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA02139
| | - Gang Chen
- Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA02139
| |
Collapse
|
2
|
Santis GD, Herman KM, Heindel JP, Xantheas SS. Descriptors of water aggregation. J Chem Phys 2024; 160:054306. [PMID: 38341703 DOI: 10.1063/5.0179815] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2023] [Accepted: 01/05/2024] [Indexed: 02/13/2024] Open
Abstract
We rely on a total of 23 (cluster size, 8 structural, and 14 connectivity) descriptors to investigate structural patterns and connectivity motifs associated with water cluster aggregation. In addition to the cluster size n (number of molecules), the 8 structural descriptors can be further categorized into (i) one-body (intramolecular): covalent OH bond length (rOH) and HOH bond angle (θHOH), (ii) two-body: OO distance (rOO), OHO angle (θOHO), and HOOX dihedral angle (ϕHOOX), where X lies on the bisector of the HOH angle, (iii) three-body: OOO angle (θOOO), and (iv) many-body: modified tetrahedral order parameter (q) to account for two-, three-, four-, five-coordinated molecules (qm, m = 2, 3, 4, 5) and radius of gyration (Rg). The 14 connectivity descriptors are all many-body in nature and consist of the AD, AAD, ADD, AADD, AAAD, AAADD adjacencies [number of hydrogen bonds accepted (A) and donated (D) by each water molecule], Wiener index, Average Shortest Path Length, hydrogen bond saturation (% HB), and number of non-short-circuited three-membered cycles, four-membered cycles, five-membered cycles, six-membered cycles, and seven-membered cycles. We mined a previously reported database of 4 948 959 water cluster minima for (H2O)n, n = 3-25 to analyze the evolution and correlation of these descriptors for the clusters within 5 kcal/mol of the putative minima. It was found that rOH and % HB correlated strongly with cluster size n, which was identified as the strongest predictor of energetic stability. Marked changes in the adjacencies and cycle count were observed, lending insight into changes in the hydrogen bond network upon aggregation. A Principal Component Analysis (PCA) was employed to identify descriptor dependencies and group clusters into specific structural patterns across different cluster sizes. The results of this study inform our understanding of how water clusters evolve in size and what appropriate descriptors of their structural and connectivity patterns are with respect to system size, stability, and similarity. The approach described in this study is general and can be easily extended to other hydrogen-bonded systems.
Collapse
Affiliation(s)
- Garrett D Santis
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MSIN J7-10, Richland, Washington 99352, USA
- Computational and Theoretical Chemistry Institute (CTCI), Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
3
|
Herman KM, Xantheas SS. An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2-25. Phys Chem Chem Phys 2023; 25:7120-7143. [PMID: 36853239 DOI: 10.1039/d2cp03241d] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2-11, 16-17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (De's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2-25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83-119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (De's) throughout the entire cluster range (n = 2-25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ∼4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6-3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.
Collapse
Affiliation(s)
- Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA.
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA. .,Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, WA, 99352, USA.
| |
Collapse
|
4
|
Tokmachev AM. Statistical Equivalence of Quantum Chemical Methods for Energy Distribution in Water Clusters. Chemphyschem 2023; 24:e202300019. [PMID: 36787108 DOI: 10.1002/cphc.202300019] [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: 01/09/2023] [Revised: 02/13/2023] [Accepted: 02/14/2023] [Indexed: 02/15/2023]
Abstract
Methods of quantum chemistry are instrumental in understanding molecular structures and properties. However, the results demonstrate significant variability, which is difficult to predict and rationalize. The fundamental question is whether some molecular systems exhibit properties invariant with respect to the computational method. The idea explored here is that collective properties of statistical ensembles should be more robust than characteristics of individual molecules and their arbitrary sets. This effect is demonstrated for the complete set of hydrogen-bond topologies of the dodecahedral water cluster (H2 O)20 . Non-Gaussian energy distributions produced by various methods have the same functional form despite strong differences in mean values and standard deviations. The conclusion is tested on methods of different complexity and origin employing a number of criteria. A linear mapping between the energies produced by different methods is discussed. The significance of the results is in establishing a collective equivalence property of quantum chemical methods.
Collapse
Affiliation(s)
- Andrey M Tokmachev
- National Research Center "Kurchatov Institute", Kurchatov Sq. 1, 123182, Moscow, Russia
| |
Collapse
|
5
|
Becker D, Dierking CW, Suchan J, Zurheide F, Lengyel J, Fárník M, Slavíček P, Buck U, Zeuch T. Temperature evolution in IR action spectroscopy experiments with sodium doped water clusters. Phys Chem Chem Phys 2021; 23:7682-7695. [PMID: 33496289 DOI: 10.1039/d0cp05390b] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The combination of supersonic expansions with IR action spectroscopy techniques is the basis of many successful approaches to study cluster structure and dynamics. The effects of temperature and temperature evolution are important with regard to both the cluster synthesis and the cluster dynamics upon IR excitation. In the past the combination of the sodium doping technique with IR excitation enhanced near threshold photoionization has been successfully applied to study neutral, especially water clusters. In this work we follow an overall examination approach for inspecting the interplay of cluster temperature and cluster structure in the initial cooling process and in the IR excitation induced heating of the clusters. In molecular simulations, we study the temperature dependent photoionization spectra of the sodium doped clusters and the evaporative cooling process by water molecule ejection at the cluster surface. We present a comprehensive analysis that provides constraints for the temperature evolution from the nozzle to cluster detection in the mass spectrometer. We attribute the IR action effect to the strong temperature dependence of sodium solvation in the IR excited clusters and we discuss the effects of geometry changes during the IR multi-photon absorption process with regard to application prospects of the method.
Collapse
Affiliation(s)
- Daniel Becker
- Universität Göttingen, Institut für Physikalische Chemie, Tammannstr. 6, 37077 Göttingen, Germany.
| | | | | | | | | | | | | | | | | |
Collapse
|
6
|
Abstract
Many efforts undertaken to study the solvation process have led to general theories that may describe mean properties, but are unable to provide a detailed understanding at the molecular level. Remarkably, the basic question of how many solvent molecules are necessary to solvate one solute molecule is still open. By exploring several water aggregates of increasing complexity, in this contribution we employ semiclassical spectroscopy to determine on quantum dynamical grounds the minimal network of surrounding water molecules to make the central one display the same vibrational features of liquid water. We find out that double-acceptor double-donor tetrahedral coordination constituting the standard picture is necessary but not sufficient, and that particular care must be reserved for the quantum description of the combination band due to the coupling of the central monomer bending mode with network librations. It is actually our ability to investigate the combination band with a quantum-derived approach that allows us to answer the titular question. The minimal structure eventually responsible for proper solvation is made of a total of 21 water molecules and includes two complete solvation shells, of which the whole first one is tetrahedrally coordinated to the central molecule. How quantum spectroscopic simulations can explain water solvation by comparison with experimental spectra.![]()
Collapse
Affiliation(s)
- Alessandro Rognoni
- Dipartimento di Chimica, Università degli Studi di Milano Via Golgi 19 20133 Milano Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano Via Golgi 19 20133 Milano Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano Via Golgi 19 20133 Milano Italy
| |
Collapse
|
7
|
Rakshit A, Bandyopadhyay P, Heindel JP, Xantheas SS. Atlas of putative minima and low-lying energy networks of water clusters n = 3-25. J Chem Phys 2019; 151:214307. [PMID: 31822087 DOI: 10.1063/1.5128378] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report a database consisting of the putative minima and ∼3.2 × 106 local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3-25 using an improved version of the Monte Carlo temperature basin paving (MCTBP) global optimization procedure in conjunction with the ab initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as low-lying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation (MP2)/aug-cc-pVTZ level of theory. In total, we have identified 3 138 303 networks corresponding to local minima of the clusters n = 3-25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1-3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but not limited to neural networks and machine and deep learning) to understand the properties of water, nature's most important substance.
Collapse
Affiliation(s)
- Avijit Rakshit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
8
|
Abstract
The appearance of ice I in the smallest possible clusters and the nature of its phase coexistence with liquid water could not thus far be unraveled. The experimental and theoretical infrared spectroscopic and free-energy results of this work show the emergence of the characteristic hydrogen-bonding pattern of ice I in clusters containing only around 90 water molecules. The onset of crystallization is accompanied by an increase of surface oscillator intensity with decreasing surface-to-volume ratio, a spectral indicator of nanoscale crystallinity of water. In the size range from 90 to 150 water molecules, we observe mixtures of largely crystalline and purely amorphous clusters. Our analysis suggests that the liquid-ice I transition in clusters loses its sharp 1st-order character at the end of the crystalline-size regime and occurs over a range of temperatures through heterophasic oscillations in time, a process without analog in bulk water.
Collapse
|
9
|
Zhao W, Huang H, Bi Q, Xu Y, Lü Y. One-dimensional water nanowires induced by electric fields. Phys Chem Chem Phys 2019; 21:19414-19422. [PMID: 31460524 DOI: 10.1039/c9cp02788b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Self-aggregation of water vapour molecules under external electric fields is systemically investigated by using molecular dynamics simulations. It is found that small water clusters aggregate into one-dimensional water nanowires along the electric field direction. The electric field strength plays a crucial role in tuning the nanowire structure. Under relatively weak electric fields such as E = 0.1 V Å-1, square and pentagonal prism-like structures are preferred; when intermediate strength electric fields are applied (E = 1.0 V Å-1), water nanowires featuring a disordered mixture of four-, five- and six-membered rings are formed; and an open ordered structure which is reminiscent of two-dimensional (2D) ice is observed when the field strength becomes very high (E > 3.0 V Å-1). Bond parameter analysis based on density-functional theory calculations shows that the electric field affects anisotropically the conformation of water molecules as well as the hydrogen-bond properties. Along the electric field, the H-O bond is stretched and the hydrogen bond shrinks with field strength in contrast to the changes perpendicular to the electric field. As a result, the hydrogen bonding is enhanced along the electric field. Under very high electric fields, the anisotropic hydrogen-bond network opens up via breaking of the bonds perpendicular to the electric field and ultimately relaxes into a loose quasi-2D ordered network.
Collapse
Affiliation(s)
- Wan Zhao
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| | - Haishen Huang
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| | - Qingling Bi
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| | - Yujia Xu
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| | - Yongjun Lü
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.
| |
Collapse
|
10
|
Lippe M, Szczepaniak U, Hou GL, Chakrabarty S, Ferreiro JJ, Chasovskikh E, Signorell R. Infrared Spectroscopy and Mass Spectrometry of CO2 Clusters during Nucleation and Growth. J Phys Chem A 2019; 123:2426-2437. [DOI: 10.1021/acs.jpca.9b01030] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Martina Lippe
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog Weg 2, CH-8093 Zürich, Switzerland
| | - Urszula Szczepaniak
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog Weg 2, CH-8093 Zürich, Switzerland
| | - Gao-Lei Hou
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog Weg 2, CH-8093 Zürich, Switzerland
| | - Satrajit Chakrabarty
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog Weg 2, CH-8093 Zürich, Switzerland
| | - Jorge J. Ferreiro
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog Weg 2, CH-8093 Zürich, Switzerland
| | - Egor Chasovskikh
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog Weg 2, CH-8093 Zürich, Switzerland
| | - Ruth Signorell
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog Weg 2, CH-8093 Zürich, Switzerland
| |
Collapse
|
11
|
Dierking CW, Zurheide F, Zeuch T, Med J, Parez S, Slavíček P. Revealing isomerism in sodium-water clusters: Photoionization spectra of Na(H2O)n (n = 2–90). J Chem Phys 2017; 146:244303. [DOI: 10.1063/1.4986520] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Affiliation(s)
- Christoph W. Dierking
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, D-37077 Göttingen, Germany
| | - Florian Zurheide
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, D-37077 Göttingen, Germany
| | - Thomas Zeuch
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, D-37077 Göttingen, Germany
| | - Jakub Med
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, Prague 6, Czech Republic
| | - Stanislav Parez
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, Prague 6, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, Prague 6, Czech Republic
| |
Collapse
|
12
|
Douberly GE, Miller RE, Xantheas SS. Formation of Exotic Networks of Water Clusters in Helium Droplets Facilitated by the Presence of Neon Atoms. J Am Chem Soc 2017; 139:4152-4156. [DOI: 10.1021/jacs.7b00510] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gary E. Douberly
- Department
of Chemistry, University of Georgia, Athens, Georgia 30602-2556, United States
| | - Roger E. Miller
- Department
of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - Sotiris S. Xantheas
- Physical
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
13
|
Wang B, Jiang W, Gao Y, Zhang Z, Sun C, Liu F, Wang Z. Energetics competition in centrally four-coordinated water clusters and Raman spectroscopic signature for hydrogen bonding. RSC Adv 2017. [DOI: 10.1039/c6ra28335g] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Viaseparating the H-bonded neighbour molecules of centrally four-coordinated water molecules from other molecules in outer cages, the calculations discover these two regions interact competitively with the central molecule.
Collapse
Affiliation(s)
- Bo Wang
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun 130012
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy
| | - Wanrun Jiang
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun 130012
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy
| | - Yang Gao
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun 130012
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy
| | - Zhiyuan Zhang
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun 130012
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy
| | - Changqing Sun
- School of Electrical and Electronic Engineering
- Nanyang Technological University
- Singapore 639798
- Singapore
| | - Fang Liu
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun 130012
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy
| | - Zhigang Wang
- Institute of Atomic and Molecular Physics
- Jilin University
- Changchun 130012
- China
- Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy
| |
Collapse
|
14
|
DiTucci MJ, Williams ER. Nanometer patterning of water by tetraanionic ferrocyanide stabilized in aqueous nanodrops. Chem Sci 2016; 8:1391-1399. [PMID: 28451280 PMCID: PMC5361863 DOI: 10.1039/c6sc03722d] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2016] [Accepted: 10/16/2016] [Indexed: 11/21/2022] Open
Abstract
Formation of the small, highly charged tetraanion ferrocyanide, Fe(CN)64–, stabilized in aqueous nanodrops and its influence to the surrounding hydrogen-bonding network of water is reported.
Formation of the small, highly charged tetraanion ferrocyanide, Fe(CN)64–, stabilized in aqueous nanodrops is reported. Ion–water interactions inside these nanodrops are probed using blackbody infrared radiative dissociation, infrared photodissociation (IRPD) spectroscopy, and molecular modeling in order to determine how water molecules stabilize this highly charged anion and the extent to which the tetraanion patterns the hydrogen-bonding network of water at long distance. Fe(CN)64–(H2O)38 is the smallest cluster formed directly by nanoelectrospray ionization. Ejection of an electron from this ion to form Fe(CN)63–(H2O)38 occurs with low-energy activation, but loss of a water molecule is favored at higher energy indicating that water molecule loss is entropically favored over loss of an electron. The second solvation shell is almost complete at this cluster size indicating that nearly two solvent shells are required to stabilize this highly charged anion. The extent of solvation necessary to stabilize these clusters with respect to electron loss is substantially lower through ion pairing with either H+ or K+ (n = 17 and 18, respectively). IRPD spectra of Fe(CN)64–(H2O)n show the emergence of a free O–H water molecule stretch between n = 142 and 162 indicating that this ion patterns the structure of water molecules within these nanodrops to a distance of at least ∼1.05 nm from the ion. These results provide new insights into how water stabilizes highly charged ions and demonstrate that highly charged anions can have a significant effect on the hydrogen-bonding network of water molecules well beyond the second and even third solvation shells.
Collapse
Affiliation(s)
- Matthew J DiTucci
- Department of Chemistry , University of California , B42 Hildebrand Hall , Berkeley , CA 94270 , USA .
| | - Evan R Williams
- Department of Chemistry , University of California , B42 Hildebrand Hall , Berkeley , CA 94270 , USA .
| |
Collapse
|
15
|
Miliordos E, Aprà E, Xantheas SS. A New, Dispersion-Driven Intermolecular Arrangement for the Benzene–Water Octamer Complex: Isomers and Analysis of their Vibrational Spectra. J Chem Theory Comput 2016; 12:4004-14. [DOI: 10.1021/acs.jctc.6b00668] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Evangelos Miliordos
- Physical
Sciences Division, Pacific Northwest National Laboratory, 902 Battelle
Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| | - Edoardo Aprà
- Environmental
Molecular Sciences Laboratory, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Sotiris S. Xantheas
- Physical
Sciences Division, Pacific Northwest National Laboratory, 902 Battelle
Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| |
Collapse
|
16
|
Sahu N, Gadre SR. Accurate vibrational spectra via molecular tailoring approach: a case study of water clusters at MP2 level. J Chem Phys 2016; 142:014107. [PMID: 25573553 DOI: 10.1063/1.4905004] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In spite of the recent advents in parallel algorithms and computer hardware, high-level calculation of vibrational spectra of large molecules is still an uphill task. To overcome this, significant effort has been devoted to the development of new algorithms based on fragmentation methods. The present work provides the details of an efficient and accurate procedure for computing the vibrational spectra of large clusters employing molecular tailoring approach (MTA). The errors in the Hessian matrix elements and dipole derivatives arising due to the approximation nature of MTA are reduced by grafting the corrections from a smaller basis set. The algorithm has been tested out for obtaining vibrational spectra of neutral and charged water clusters at Møller-Plesset second order level of theory, and benchmarking them against the respective full calculation (FC) and/or experimental results. For (H2O)16 clusters, the estimated vibrational frequencies are found to differ by a maximum of 2 cm(-1) with reference to the corresponding FC values. Unlike the FC, the MTA-based calculations including grafting procedure can be performed on a limited hardware, yet take a fraction of the FC time. The present methodology, thus, opens a possibility of the accurate estimation of the vibrational spectra of large molecular systems, which is otherwise impossible or formidable.
Collapse
Affiliation(s)
- Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| |
Collapse
|
17
|
|
18
|
Sahu N, Gadre SR, Rakshit A, Bandyopadhyay P, Miliordos E, Xantheas SS. Low energy isomers of (H2O)25 from a hierarchical method based on Monte Carlo temperature basin paving and molecular tailoring approaches benchmarked by MP2 calculations. J Chem Phys 2015; 141:164304. [PMID: 25362296 DOI: 10.1063/1.4897535] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We report new global minimum candidate structures for the (H2O)25 cluster that are lower in energy than the ones reported previously and correspond to hydrogen bonded networks with 42 hydrogen bonds and an interior, fully coordinated water molecule. These were obtained as a result of a hierarchical approach based on initial Monte Carlo Temperature Basin Paving sampling of the cluster's Potential Energy Surface with the Effective Fragment Potential, subsequent geometry optimization using the Molecular Tailoring Approach with the fragments treated at the second order Møller-Plesset (MP2) perturbation (MTA-MP2) and final refinement of the entire cluster at the MP2 level of theory. The MTA-MP2 optimized cluster geometries, constructed from the fragments, were found to be within <0.5 kcal/mol from the minimum geometries obtained from the MP2 optimization of the entire (H2O)25 cluster. In addition, the grafting of the MTA-MP2 energies yields electronic energies that are within <0.3 kcal/mol from the MP2 energies of the entire cluster while preserving their energy rank order. Finally, the MTA-MP2 approach was found to reproduce the MP2 harmonic vibrational frequencies, constructed from the fragments, quite accurately when compared to the MP2 ones of the entire cluster in both the HOH bending and the OH stretching regions of the spectra.
Collapse
Affiliation(s)
- Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Avijit Rakshit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Evangelos Miliordos
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| | - Sotiris S Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| |
Collapse
|
19
|
Tabor DP, Kusaka R, Walsh PS, Zwier TS, Sibert EL. Local Mode Approach to OH Stretch Spectra of Benzene–(H2O)n Clusters, n = 2–7. J Phys Chem A 2015; 119:9917-30. [PMID: 26340135 DOI: 10.1021/acs.jpca.5b06954] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Daniel P. Tabor
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706, United States
| | - Ryoji Kusaka
- Department
of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Patrick S. Walsh
- Department
of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Timothy S. Zwier
- Department
of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Edwin L. Sibert
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706, United States
| |
Collapse
|
20
|
Sahu N, Khire SS, Gadre SR. Structures, energetics and vibrational spectra of (H2O)32clusters: a journey from model potentials to correlated theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1062150] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
21
|
Gimelshein N, Gimelshein S, Pradzynski CC, Zeuch T, Buck U. The temperature and size distribution of large water clusters from a non-equilibrium model. J Chem Phys 2015; 142:244305. [DOI: 10.1063/1.4922312] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - S. Gimelshein
- University of Southern California, Los Angeles, California 90089, USA
| | - C. C. Pradzynski
- Institut für Physikalische Chemie, Universität Göttingen, Tammanstr. 6, D-37077 Göttingen, Germany
| | - T. Zeuch
- Institut für Physikalische Chemie, Universität Göttingen, Tammanstr. 6, D-37077 Göttingen, Germany
| | - U. Buck
- Max-Planck-Institut für Dynamik und Selbstorganisation, Am Faßberg 17, D-37077 Göttingen, Germany
| |
Collapse
|
22
|
Miliordos E, Xantheas SS. An accurate and efficient computational protocol for obtaining the complete basis set limits of the binding energies of water clusters at the MP2 and CCSD(T) levels of theory: Application to (H2O)m, m = 2-6, 8, 11, 16, and 17. J Chem Phys 2015; 142:234303. [DOI: 10.1063/1.4922262] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Evangelos Miliordos
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| | - Sotiris S. Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| |
Collapse
|
23
|
Dieterich JM, Hartke B. Observable-targeting global cluster structure optimization. Phys Chem Chem Phys 2015; 17:11958-61. [DOI: 10.1039/c5cp01910a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Global cluster structure optimization not towards lowest energy but towards best fit of cluster properties to experimental data provides theoretical support for cluster experiments under non-equilibrium conditions.
Collapse
Affiliation(s)
| | - Bernd Hartke
- Institute for Physical Chemistry
- Christian-Albrechts-University
- 24098 Kiel
- Germany
| |
Collapse
|
24
|
Shishido R, Li YC, Tsai CW, Bing D, Fujii A, Kuo JL. An infrared spectroscopic and theoretical study on (CH3)3N–H+–(H2O)n, n = 1–22: highly polarized hydrogen bond networks of hydrated clusters. Phys Chem Chem Phys 2015; 17:25863-76. [DOI: 10.1039/c5cp01487e] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Highly polarized water networks are found in the micro hydaration of protonated trimethylamine.
Collapse
Affiliation(s)
- Ryunosuke Shishido
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| | - Ying-Cheng Li
- Department of Physics
- National Taiwan University
- Taipei 10617
- Taiwan
- Institute of Atomic and Molecular Sciences
| | - Chen-Wei Tsai
- Department of Physics
- National Taiwan University
- Taipei 10617
- Taiwan
- Institute of Atomic and Molecular Sciences
| | - Dan Bing
- Pujiang Institute
- Nanjing Tech University
- Nanjing
- China
| | - Asuka Fujii
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences
- Academia Sinica
- Taipei
- Taiwan
| |
Collapse
|
25
|
Zurheide F, Dierking CW, Pradzynski CC, Forck RM, Flüggen F, Buck U, Zeuch T. Size-Resolved Infrared Spectroscopic Study of Structural Transitions in Sodium-Doped (H2O)n Clusters Containing 10–100 Water Molecules. J Phys Chem A 2014; 119:2709-20. [DOI: 10.1021/jp509883m] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Florian Zurheide
- Institut für Physikalische Chemie, Tammannstaße
6, Georg-August-Universität Göttingen, D-37077 Göttingen, Germany
| | - Christoph W. Dierking
- Institut für Physikalische Chemie, Tammannstaße
6, Georg-August-Universität Göttingen, D-37077 Göttingen, Germany
| | - Christoph C. Pradzynski
- Institut für Physikalische Chemie, Tammannstaße
6, Georg-August-Universität Göttingen, D-37077 Göttingen, Germany
| | - Richard M. Forck
- Institut für Physikalische Chemie, Tammannstaße
6, Georg-August-Universität Göttingen, D-37077 Göttingen, Germany
| | - Florian Flüggen
- Institut für Physikalische Chemie, Tammannstaße
6, Georg-August-Universität Göttingen, D-37077 Göttingen, Germany
| | - Udo Buck
- Max-Planck-Institut für Dynamik und Selbstorganisation, Am Faßberg 17, D-37077 Göttingen, Germany
| | - Thomas Zeuch
- Institut für Physikalische Chemie, Tammannstaße
6, Georg-August-Universität Göttingen, D-37077 Göttingen, Germany
| |
Collapse
|