DOCSLIB.ORG

Melting Dynamics of Ice in the Mesoscopic Regime

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Melting Dynamics of Ice in the Mesoscopic Regime Melting dynamics of ice in the mesoscopic regime Margherita Citronia,b,1, Samuele Fanettia,c, Naomi Falsinia, Paolo Foggia,d,e, and Roberto Binia,b aLENS–European Laboratory for Non-Linear Spectroscopy, 50019 Sesto, Florence, Italy; bDipartimento di Chimica Ugo Schiff, University of Florence, 50019 Sesto, Florence, Italy; cIstituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, 50019 Sesto, Florence, Italy; dIstituto Nazionale di Ottica, Consiglio Nazionale delle Ricerche, 50125 Florence, Italy; and eDipartimento di Chimica, University of Perugia, I-06123 Perugia, Italy Edited by Russell J. Hemley, The George Washington University, Washington, DC, and approved May 2, 2017 (received for review December 6, 2016) How does a crystal melt? How long does it take for melt nuclei The temperature increases by the redistribution of the absorbed to grow? The melting mechanisms have been addressed by sev- energy via the vibrational relaxation over the lattice modes. This eral theoretical and experimental works, covering a subnanosec- thermalization process takes ∼20 ps in water ice (13) and can ond time window with sample sizes of tens of nanometers and be slower in non–H-bonded molecular crystals. The structural thus suitable to determine the onset of the process but unable to changes following the heating pulse have been monitored in alu- unveil the following dynamics. On the other hand, macroscopic minum (14) and gold (15) by electron diffraction with a sample observations of phase transitions, with millisecond or longer time thickness of 20 nm, where the main signatures of melting com- resolution, account for processes occurring at surfaces and time plete in a few picoseconds. Water ice melting after an ultrafast limited by thermal contact with the environment. Here, we fill the T jump has been studied by time-resolved infrared spectra (13, gap between these two extremes, investigating the melting of 16), with a sample thickness of 1.6 µm and over a time window ice in the entire mesoscopic regime. A bulk ice Ih or ice VI sample of 250 ps. Here, the spectral changes attributed by the authors to is homogeneously heated by a picosecond infrared pulse, which melting mostly occur in the first 150 ps after the pulse (τ = 37 ps) delivers all of the energy necessary for complete melting. The evo- (16). At 250 ps an incomplete melting was still observed, consis- lution of melt/ice interfaces thereafter is monitored by Mie scat- tent with the amount of energy delivered. tering with nanosecond resolution, for all of the time needed for The knowledge at this picosecond–nanometer scale barely the sample to reequilibrate. The growth of the liquid domains, connects to the experiences in the macroscopic world. On the over distances of micrometers, takes hundreds of nanoseconds, macroscopic side, photographic investigations have been per- a time orders of magnitude larger than expected from simple formed with millisecond resolution to investigate premelting H-bond dynamics. structures in colloidal crystals (17), ice nucleation events upon cooling (18), or the rate and topology of ice and hydrate forma- Mie scattering j temperature jump j superheating j laser heating j tion under dynamic compression (19, 20). anvil cell Here, we use a picosecond T -jump technique and time- resolved Mie scattering to explore the melting and the successive elting of ice is one of the most common occurrences on our refreezing on bulk samples of ice Ih and ice VI. The main pecu- Mplanet, from polar ices to glaciers, to the morning frost and liarities of this experiment are the bulk sample thickness (50 µm) the preparation of our drinks. Ices are present as well in extrater- and the time sampling: We access a time window of tens of mil- restrial environments, planetary and intersellar, going through liseconds, which is the time needed to completely restore the extremely different pressure and temperature conditions, where initial conditions, with nanosecond resolution. We observe the many phase boundaries are encountered. The molecular mecha- scattering at ice/melt interfaces and the propagation of melting nisms of melting and the characteristic timescales of the process over distances of micrometers, accessing the characteristic time are not well elucidated, as for any other phase transition. In fact, of the process. As will be shown, the P, T, and V conditions in this great effort is presently being made with state of the art exper- experiment are quasi-static over the timescale of melting evolu- tion. The kinetics of melting propagation in a bulk crystal are imental techniques to understand the intrinsic kinetics of struc- CHEMISTRY tural transformations, particularly under dynamic compression, thus explored over the entire mesoscopic regime, connecting the by which phase-transition boundaries up to extreme pressures microscopic and the macroscopic realms. (P) and temperatures (T) can be accessed (1–6). Melting is the least hindered of phase transitions, requiring the disruption of Significance the crystalline order and the achievement of the local structure of the liquid phase. The complete unfolding of the melting process, from the The dynamics of melting have been studied up to the present nucleation onset to the achievement of a bulk liquid structure, in the picosecond timescale, observing structural changes on a is not known to date for the timescales and the sizes of the length scale of few nanometers, in the attempt to explain the heterogeneities involved. Here, we observe in ice the archety- fundamental mechanisms of the transition onset. In both simu- pal H-bonded molecular crystal, the evolution of solid/liquid lations and experiments, micrometric- or submicrometric-sized interfaces generated by the absorption of a picosecond infra- crystals can be rapidly heated above the melting temperature red pulse. The pulse delivers all of the energy needed for com- (Tm ), into a superheated state where melting occurs. Metals and plete melting; therefore, the lattice is homogenously heated water ice have been the test systems. Simulations have revealed from within, above its melting temperature. The growth of a nucleation and growth mechanism (7–11), and nucleation has liquid domains over distances of micrometers takes hundreds been especially addressed, in trying to determine the minimum of nanoseconds and is orders of magnitude slower than vibra- number of atoms/molecules needed to form a critical nucleus tional dynamics. for the transitions, the solid/melt interfacial energy, the tran- sient local structures achieved during the transformation, and the Author contributions: M.C., S.F., P.F., and R.B. designed research; M.C., S.F., and N.F. per- role of defects and surfaces. Experimentally, the sample can be formed research; M.C., S.F., and N.F. analyzed data; and M.C. and R.B. wrote the paper. rapidly heated by a pulse of infrared light at a wavelength that The authors declare no conflict of interest. the sample can absorb (T jump). With this “direct laser heat- This article is a PNAS Direct Submission. ing” the sample is heated from within, and transition nuclei can 1To whom correspondence should be addressed. Email: [email protected]fi.it. be formed in the bulk, overcoming the dominance of heteroge- This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. neous nucleation, induced by surfaces or preexisting defects (12). 1073/pnas.1620039114/-/DCSupplemental. www.pnas.org/cgi/doi/10.1073/pnas.1620039114 PNAS j June 6, 2017 j vol. 114 j no. 23 j 5935–5940 Downloaded by guest on September 25, 2021 The experimental insight into the melting dynamics of ices with which we measure the transmission through the sample with has paramount relevance in understanding natural phenomena a fast Si photodiode (Fig. 1C). The focus of the probe (∼150 µm and is a fundamental benchmark for theories at any level of in diameter) lies inside the area irradiated by the pump. When a accuracy. Ices, ubiquitous in terrestrial and extraterrestrial envi- pump pulse is shot, the probe transmission drops to a minimum ronments, are characterized by hydrogen-bond networks that in tens or hundreds of nanoseconds, depending on the amount determine the great complexity of their physical and chemical of deposited energy, and then returns to a static value in a time behavior. The melting line of water borders at least five crys- (micro- to milliseconds) strongly dependent on T0. The trans- talline phases (Ih , III, V, VI, and VII) with different structures mission change is due to scattering at ice/melt interfaces, form- and H-bond arrangements. In this work we investigate the melt- ing with the T jump and evolving in time during melting (due ing dynamics in ice Ih and VI, accessed respectively by the use to superheating) and refreezing (due to the thermal contact with of low temperature and high static pressure, providing an unique the environment that decreases T to the thermal bath tempera- experimental study in the mesoscopic domain and opening the ture T0). When the transmitted beam is blocked (Fig. 1D), a pos- way to investigations in different P and T conditions. itive signal due to off-axis scattered light is consistently observed, with a similar time dependence to that of the transmission change. Results In both configurations the signal amplitude increases as λprobe Temperature Jump and Transient Transmission Loss. Ice Ih and ice decreases, as expected for scattering. As also expected, the trans- VI samples 50 µm thick, at fixed temperature (T0) and pres- mission change is not observed when liquid samples are heated sure (P0), are “instantaneously” heated by a 15-ps pulse of by the T -jump pulse, except at the highest E values (E ≥ 4.5 mJ) light at 1,930 nm (FWHM = 100 nm) (21) with a focus diame- where a strong oscillation of the transmission is observed, likely ter of 400 µm.
Load more

Recommended publications
  • Ice Ic” Werner F
    Extent and relevance of stacking disorder in “ice Ic” Werner F. Kuhsa,1, Christian Sippela,b, Andrzej Falentya, and Thomas C. Hansenb aGeoZentrumGöttingen Abteilung Kristallographie (GZG Abt. Kristallographie), Universität Göttingen, 37077 Göttingen, Germany; and bInstitut Laue-Langevin, 38000 Grenoble, France Edited by Russell J. Hemley, Carnegie Institution of Washington, Washington, DC, and approved November 15, 2012 (received for review June 16, 2012) “ ” “ ” A solid water phase commonly known as cubic ice or ice Ic is perfectly cubic ice Ic, as manifested in the diffraction pattern, in frequently encountered in various transitions between the solid, terms of stacking faults. Other authors took up the idea and liquid, and gaseous phases of the water substance. It may form, attempted to quantify the stacking disorder (7, 8). The most e.g., by water freezing or vapor deposition in the Earth’s atmo- general approach to stacking disorder so far has been proposed by sphere or in extraterrestrial environments, and plays a central role Hansen et al. (9, 10), who defined hexagonal (H) and cubic in various cryopreservation techniques; its formation is observed stacking (K) and considered interactions beyond next-nearest over a wide temperature range from about 120 K up to the melt- H-orK sequences. We shall discuss which interaction range ing point of ice. There was multiple and compelling evidence in the needs to be considered for a proper description of the various past that this phase is not truly cubic but composed of disordered forms of “ice Ic” encountered. cubic and hexagonal stacking sequences. The complexity of the König identified what he called cubic ice 70 y ago (11) by stacking disorder, however, appears to have been largely over- condensing water vapor to a cold support in the electron mi- looked in most of the literature.
    [Show full text]
  • A Primer on Ice
    A Primer on Ice L. Ridgway Scott University of Chicago Release 0.3 DO NOT DISTRIBUTE February 22, 2012 Contents 1 Introduction to ice 1 1.1 Lattices in R3 ....................................... 2 1.2 Crystals in R3 ....................................... 3 1.3 Comparingcrystals ............................... ..... 4 1.3.1 Quotientgraph ................................. 4 1.3.2 Radialdistributionfunction . ....... 5 1.3.3 Localgraphstructure. .... 6 2 Ice I structures 9 2.1 IceIh........................................... 9 2.2 IceIc........................................... 12 2.3 SecondviewoftheIccrystalstructure . .......... 14 2.4 AlternatingIh/Iclayeredstructures . ........... 16 3 Ice II structure 17 Draft: February 22, 2012, do not distribute i CONTENTS CONTENTS Draft: February 22, 2012, do not distribute ii Chapter 1 Introduction to ice Water forms many different crystal structures in its solid form. These provide insight into the potential structures of ice even in its liquid phase, and they can be used to calibrate pair potentials used for simulation of water [9, 14, 15]. In crowded biological environments, water may behave more like ice that bulk water. The different ice structures have different dielectric properties [16]. There are many crystal structures of ice that are topologically tetrahedral [1], that is, each water molecule makes four hydrogen bonds with other water molecules, even though the basic structure of water is trigonal [3]. Two of these crystal structures (Ih and Ic) are based on the same exact local tetrahedral structure, as shown in Figure 1.1. Thus a subtle understanding of structure is required to differentiate them. We refer to the tetrahedral structure depicted in Figure 1.1 as an exact tetrahedral structure. In this case, one water molecule is in the center of a square cube (of side length two), and it is hydrogen bonded to four water molecules at four corners of the cube.
    [Show full text]
  • Arxiv:2004.08465V2 [Cond-Mat.Stat-Mech] 11 May 2020
    Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations Pablo M. Piaggi1 and Roberto Car2 1)Department of Chemistry, Princeton University, Princeton, NJ 08544, USA a) 2)Department of Chemistry and Department of Physics, Princeton University, Princeton, NJ 08544, USA (Dated: 13 May 2020) We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. We analyze the effect of the system size on the ice Ih-liquid free energy differences and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih. I. INTRODUCTION ture forms in an orientation compatible with the simulation box9. The study of phase equilibria using computer simulations is of central importance to understand the behavior of a given model. However, finding the thermodynamic condition at II. CRYSTAL STRUCTURE OF ICE Ih which two or more phases coexist is particularly hard in the presence of first order phase transitions.
    [Show full text]
  • Modeling the Ice VI to VII Phase Transition
    Modeling the Ice VI to VII Phase Transition Dawn M. King 2009 NSF/REU PROJECT Physics Department University of Notre Dame Advisor: Dr. Kathie E. Newman July 31, 2009 Abstract Ice (solid water) is found in a number of different structures as a function of temperature and pressure. This project focuses on two forms: Ice VI (space group P 42=nmc) and Ice VII (space group Pn3m). An interesting feature of the structural phase transition from VI to VII is that both structures are \self clathrate," which means that each structure has two sublattices which interpenetrate each other but do not directly bond with each other. The goal is to understand the mechanism behind the phase transition; that is, is there a way these structures distort to become the other, or does the transition occur through the breaking of bonds followed by a migration of the water molecules to the new positions? In this project we model the transition first utilizing three dimensional visualization of each structure, then we mathematically develop a common coordinate system for the two structures. The last step will be to create a phenomenological Ising-like spin model of the system to capture the energetics of the transition. It is hoped the spin model can eventually be studied using either molecular dynamics or Monte Carlo simulations. 1 Overview of Ice The known existence of many solid states of water provides insight into the complexity of condensed matter in the universe. The familiarity of ice and the existence of many structures deem ice to be interesting in the development of techniques to understand phase transitions.
    [Show full text]
  • 11Th International Conference on the Physics and Chemistry of Ice, PCI
    11th International Conference on the Physics and Chemistry of Ice (PCI-2006) Bremerhaven, Germany, 23-28 July 2006 Abstracts _______________________________________________ Edited by Frank Wilhelms and Werner F. Kuhs Ber. Polarforsch. Meeresforsch. 549 (2007) ISSN 1618-3193 Frank Wilhelms, Alfred-Wegener-Institut für Polar- und Meeresforschung, Columbusstrasse, D-27568 Bremerhaven, Germany Werner F. Kuhs, Universität Göttingen, GZG, Abt. Kristallographie Goldschmidtstr. 1, D-37077 Göttingen, Germany Preface The 11th International Conference on the Physics and Chemistry of Ice (PCI- 2006) took place in Bremerhaven, Germany, 23-28 July 2006. It was jointly organized by the University of Göttingen and the Alfred-Wegener-Institute (AWI), the main German institution for polar research. The attendance was higher than ever with 157 scientists from 20 nations highlighting the ever increasing interest in the various frozen forms of water. As the preceding conferences PCI-2006 was organized under the auspices of an International Scientific Committee. This committee was led for many years by John W. Glen and is chaired since 2002 by Stephen H. Kirby. Professor John W. Glen was honoured during PCI-2006 for his seminal contributions to the field of ice physics and his four decades of dedicated leadership of the International Conferences on the Physics and Chemistry of Ice. The members of the International Scientific Committee preparing PCI-2006 were J.Paul Devlin, John W. Glen, Takeo Hondoh, Stephen H. Kirby, Werner F. Kuhs, Norikazu Maeno, Victor F. Petrenko, Patricia L.M. Plummer, and John S. Tse; the final program was the responsibility of Werner F. Kuhs. The oral presentations were given in the premises of the Deutsches Schiffahrtsmuseum (DSM) a few meters away from the Alfred-Wegener-Institute.
    [Show full text]
  • 2Growth, Structure and Properties of Sea
    Growth, Structure and Properties 2 of Sea Ice Chris Petrich and Hajo Eicken 2.1 Introduction The substantial reduction in summer Arctic sea ice extent observed in 2007 and 2008 and its potential ecological and geopolitical impacts generated a lot of attention by the media and the general public. The remote-sensing data documenting such recent changes in ice coverage are collected at coarse spatial scales (Chapter 6) and typically cannot resolve details fi ner than about 10 km in lateral extent. However, many of the processes that make sea ice such an important aspect of the polar oceans occur at much smaller scales, ranging from the submillimetre to the metre scale. An understanding of how large-scale behaviour of sea ice monitored by satellite relates to and depends on the processes driving ice growth and decay requires an understanding of the evolution of ice structure and properties at these fi ner scales, and is the subject of this chapter. As demonstrated by many chapters in this book, the macroscopic properties of sea ice are often of most interest in studies of the interaction between sea ice and its environment. They are defi ned as the continuum properties averaged over a specifi c volume (Representative Elementary Volume) or mass of sea ice. The macroscopic properties are determined by the microscopic structure of the ice, i.e. the distribution, size and morphology of ice crystals and inclusions. The challenge is to see both the forest, i.e. the role of sea ice in the environment, and the trees, i.e. the way in which the constituents of sea ice control key properties and processes.
    [Show full text]
  • Mechanical Characterization Via Full Atomistic Simulation: Applications to Nanocrystallized Ice
    Mechanical Characterization via Full Atomistic Simulation: Applications to Nanocrystallized Ice A thesis presented By Arvand M.H. Navabi to The Department of Civil and Environmental Engineering in partial fulfillment of the requirements for the degree of Master of Science in the field of Civil Engineering Northeastern University Boston, Massachusetts August, 2016 Submitted to Prof. Steven W. Cranford Acknowledgements I acknowledge generous support from my thesis advisor Dr. Steven W. Cranford whose encouragement and availability was crucial to this thesis and also my parents for allowing me to realize my own potential. The simulations were made possible by LAMMPS open source program. Visualization has been carried out using the VMD visualization package. 3 Abstract This work employs molecular dynamic (MD) approaches to characterize the mechanical properties of nanocrystalline materials via a full atomistic simulation using the ab initio derived ReaxFF potential. Herein, we demonstrate methods to efficiently simulate key mechanical properties (ultimate strength, stiffness, etc.) in a timely and computationally inexpensive manner. As an illustrative example, the work implements the described methodology to perform full atomistic simulation on ice as a material platform, which — due to its complex behavior and phase transitions upon pressure, heat exchange, energy transfer etc. — has long been avoided or it has been unsuccessful to ascertain its mechanical properties from a molecular perspective. This study will in detail explain full atomistic MD methods and the particulars required to correctly simulate crystalline material systems. Tools such as the ReaxFF potential and open-source software package LAMMPS will be described alongside their fundamental theories and suggested input methods to simulate further materials, encompassing both periodic and finite crystalline models.
    [Show full text]
  • Ice-Seven (Ice VII)
    Ice-seven (Ice VII) Ice-seven (ice VII) [1226] is formed from liquid water above 3 GPa by lowering its temperature to ambient temperatures (see Phase Diagram). It can be obtained at low temperature and ambient pressure by decompressing (D2O) ice-six below 95 K and is metastable over a wide range of pressure, transforming into LDA above 120 K [948]. Note that in this structural diagram the hydrogen bonding is ordered whereas in reality it is random (obeying the 'ice rules': two hydrogen atoms near each oxygen, one hydrogen atom on each O····O bond). As the H-O-H angle does not vary much from that of the isolated molecule, the hydrogen bonds are not straight (although shown so in the figures). The Ice VII unit cell, which forms cubic crystal ( , 224; Laue class symmetry m-3m) consists of two interpenetrating cubic ice lattices with hydrogen bonds passing through the center of the water hexamers and no connecting hydrogen-bonds between lattices. It has a density of about 1.65 g cm- 3 (at 2.5 GPa and 25 °C [8]), which is less than twice the cubic ice density as the intra-network O····O distances are 8% longer (at 0.1 MPa) to allow for the interpenetration. The cubic crystal (shown opposite) has cell dimensions 3.3501 Å (a, b, c, 90º, 90º, 90º; D2O, at 2.6 GPa and 22 °C [361]) and contains two water molecules. All molecules experience identical molecular environments. The hydrogen bonding is disordered and constantly changing as in hexagonal ice but ice-seven undergoes a proton disorder-order transition to ice-eight at about 5 °C; ice-seven and ice-eight having identical structures apart from the proton ordering.
    [Show full text]
  • © Cambridge University Press Cambridge
    Cambridge University Press 978-0-521-80620-6 - Creep and Fracture of Ice Erland M. Schulson and Paul Duval Index More information Index 100-year wave force, 336 friction and fracture, 289, 376 60° dislocations, 17, 82, 88 indentation failure, 345 microstructure, 45, 70, 237, 255, 273 abrasion, 337 multiscale fracture and frictional accommodation processes of basal slip, 165 sliding, 386 acoustic emission, 78, 90, 108 nested envelopes, 377 across-column cleavage cracks, 278 pressure–area relationship, 349, 352 across-column confinement, 282 S2 growth texture, 246, 273 across-column cracks, 282, 306, SHEBA faults, 371 across-column loading, 275 SHEBA stress states, 377 across-column strength, 246, 249, 275 Arctic Ocean, 1, 45, 190, 361 activation energy, 71, 84, 95, 111, 118, 131 aspect ratio, 344 activation volume, 114, 182 atmospheric ice, 219, 221, 241, 243 activity of pyramidal slip systems, 158 atmospheric icing, 31 activity of slip systems, 168 atmospheric impurities, 113 adiabatic heating, 291, 348 atomic packing factor, 9 adiabatic softening, 291 audible report, 240 affine self-consistent model, 160 avalanches, 206 air bubbles, 38 air-hydrate crystals, 37 bands, 89 albedo, 363 basal activity, 162 aligned first-year sea ice, 246 basal dislocations, 77 along-column confinement, 282 basal planes, 214, along-column confining stress, 282, basal screw dislocations, 77, 87 along-column strength, 244, 275 basal shear bands, 163 ammonia dihydrate, 181 basal slip, 18, 77, 127, 228 ammonia–water system, 186 basal slip lines, 77 amorphous forms
    [Show full text]
  • Hexagonal Ice (Ice Ih)
    Hexagonal Ice (Ice Ih) Some physical properties Ice nucleation and growth Is ice slippery? Hexagonal ice (ice Ih) [1969]. is the form of all natural snow and ice on Earth (see Phase Diagram), as evidenced in the six-fold symmetry in ice crystals grown from water vapor (that is, snow flakes). Hexagonal ice (Space group P63/mmc, 194; symmetry D6h, Laue class symmetry 6/mmm; analogous to β-tridymite silica or lonsdaleite) possesses a fairly open low-density structure, where the packing efficiency is low (~1/3) compared with simple cubic (~1/2) or face centered cubic (~3/4) structuresa (and in contrast to face centered cubic close packed solid hydrogen sulfide). The crystals may be thought of as consisting of sheets lying on top of each other. The basic structure consists of a hexameric box where planes consist of chair-form hexamers (the two horizontal planes, opposite) or boat-form hexamers (the three vertical planes, opposite). In this diagram the hydrogen bonding is shown ordered whereas in reality it is random,b as protons can move between (ice) water molecules at temperatures above about 130 K [1504]. The water molecules have a staggered arrangement of hydrogen bonding with respect to three of their neighbors, in the plane of the chair-form hexamers. The fourth neighbor (shown as vertical links opposite) has an eclipsed arrangement of hydrogen bonding. There is a small deviation from ideal hexagonal symmetry as the unit cellc is 0.3 % shorter in the c- direction (in the direction of the eclipsed hydrogen bonding, shown as vertical links in the figures).
    [Show full text]
  • The Fracture of Water Ice Ih: a Short Overview
    Meteoritics & Planetary Science 41, Nr 10, 1497–1508 (2006) Abstract available online at http://meteoritics.org The fracture of water ice Ih: A short overview Erland M. SCHULSON Thayer School of Engineering, Dartmouth College, Hanover, New Hampshire 03755, USA E-mail: [email protected] (Received 25 October 2005; revision accepted 09 April 2006) Abstract–This paper presents a short overview of the fracture of water ice Ih. Topics include the ductile-to-brittle transition, tensile and compressive strength, compressive failure under multiaxial loading, compressive failure modes, and brittle failure on the geophysical scale (Arctic sea ice cover, Europa’s icy crust). Emphasis is placed on the underlying physical mechanisms. Where appropriate, comment is made on the formation of high-latitude impact craters on Mars. INTRODUCTION (Weeks and Gow 1980). (The terms S2 and S3 are taken from the classification of texture by Michel and Ramseier 1971). Water ice can adopt one of 13 crystalline forms or one of Subsurface Martian ice may be a textured polycrystal as well, least two amorphous states (Petrenko and Whitworth 1999). should it have formed from the directional solidification of Under low pressure and at temperatures above about −120 °C, water. water ice (or simply ice) adopts a hexagonal crystal structure, Whether in the form of a single crystal or a polycrystal, ice denoted Ih, reflected in the shape of snowflakes. Each exhibits two kinds of inelastic behavior. When slowly loaded it molecule within the Ih unit cell is connected to four others flows plastically (i.e., creeps). This is manifested, for instance, (owing to the 104.5° H-O-H bond angle), and each unit cell by strains in excess of unity within alpine glaciers where the contains four molecules.
    [Show full text]
  • Permeation of Light Gases Through Hexagonal Ice
    Materials 2012, 5, 1593-1601; doi:10.3390/ma5091593 OPEN ACCESS materials ISSN 1996-1944 www.mdpi.com/journal/materials Article Permeation of Light Gases through Hexagonal Ice Joana Durão 1 and Luis Gales 1,2,* 1 IBMC—Institute for Molecular and Cell Biology, University of Porto Porto, Rua do Campo Alegre 823, Porto 4150-180, Portugal; E-Mail: [email protected] 2 ICBAS—Institute of Biomedical Sciences Abel Salazar, Porto 4099-003, Portugal * Author to whom correspondence should be addressed; E-Mail: [email protected]; Tel.: +351-222-062-200. Fax: +351-222-062-232. Received: 16 July 2012; in revised form: 21 August 2012 / Accepted: 27 August 2012 / Published: 5 September 2012 Abstract: Gas separation using porous solids have attracted great attention due to their energetic applications. There is an enormous economic and environmental interest in the development of improved technologies for relevant processes, such as H2 production, CO2 separation or O2 and N2 purification from air. New materials are needed for achieving major improvements. Crystalline materials, displaying unidirectional and single-sized pores, preferentially with low pore tortuosity and high pore density, are promising candidates for membrane synthesis. Herein, we study hexagonal ice crystals as an example of this class of materials. By slowly growing ice crystals inside capillary tubes we were able to measure the permeation of several gas species through ice crystals and investigate its relation with both the size of the guest molecules and temperature of the crystal. Keywords: ice; light gases; diffusion 1. Introduction Separation or purification of light gases is crucial in many industrial activities.
    [Show full text]