Isomorphism between ice and silica

Both ice and silica crystallize into solid-state structures composed of tetrahedral building units that are joined together to form an infinite four-connected net. Mathematical considerations suggest that there is a vast number of such nets and thus potential crystal structures. It is therefore perhaps surprising to discover that, despite the differences in the nature of interatomic interactions in these materials, a fair number of commonly observed ice and silica phases are based on common nets. Here we use computer simulation to investigate the origin of this symmetry between the structures formed for ice and silica and to attempt to understand why it is not complete. We start from a comparison of the dense phases and then move to the relationship between the different open (zeolitic and clathratic) structures formed for both materials. We show that there is a remarkably strong correlation between the energetics of isomorphic silica and water ice structures and that this correlation arises because of the strong link between the total energy of a material and its local geometric features. Finally, we discuss a number of as yet unsynthesized low-energy structures which include a phase of ice based on quartz, a silica based on the structure of ice VI, and an ice clathrate that is isomorphic to the silicate structure nonasil.

Surface energy and surface proton order of the ice Ih basal and prism surfaces

Density-functional theory (DFT) is used to examine the basal and prism surfaces of ice Ih. Similar surface energies are obtained for the two surfaces; however, in each case a strong dependence of the surface energy on surface proton order is identified. This dependence, which can be as much as 50% of the absolute surface energy, is significantly larger than the bulk dependence (< 1%) on proton order, suggesting that the thermodynamic ground state of the ice surface will remain proton ordered well above the bulk order-disorder temperature of about 72 K. On the basal surface this suggestion is supported by Monte Carlo simulations with an empirical potential and solution of a 2D Ising model with nearest neighbor interactions taken from DFT. Order parameters that define the surface energy of each surface in terms of nearest neighbor interactions between dangling OH bonds (those which point out of the surface into vacuum) have been identified and are discussed. Overall, these results suggest that proton order-disorder effects have a profound impact on the stability of ice surfaces and will most likely have an effect on ice surface reactivity as well as ice crystal growth and morphology. S Supplementary data are available from stacks.iop.org/JPhysCM/22/074209/mmedia

On thin ice: surface order and disorder during pre-melting

The effect of temperature on the structure of the ice Ih (0001) surface is considered through a series of molecular dynamics simulations on an ice slab. At relatively low temperatures (200K) a small fraction of surface self-interstitials (i.e. admolecules) appear that are formed exclusively from molecules leaving the outermost bilayer. At higher temperatures (ca. 250 K), vacancies start to appear in the inner part of the outermost bilayer exposing the underlying bilayer and providing sites with a high concentration of the dangling hydrogen bonds. Around 250-260 K aggregates of molecules formed on top of the outermost bilayer from self-interstitials become more mobile and have diffusivities approaching that of liquid water. At similar to 270-280 K the inner bilayer of one surface noticeably destructures and it appears that at above 285 K both surfaces are melting. The observed disparity in the onset of melting between the two sides of the slab is rationalised by considering the relationship between surface energy and the spatial distribution of protons at the surface; thermodynamic stability is conferred on the surface by maximising separations between dangling protons at the crystal exterior. Local hotspots associated with a high dangling proton density are suggested to be susceptible to pre-melting and may be more efficient at trapping species at the external surface than regions with low concentrations of protons thus potentially helping ice particles to catalyse reactions. A preliminary conclusion of this work is that only about 10-20 K below the melting temperature of the particular water potential employed is major disruption of the crystalline lattice noted which could be interpreted as being "liquid", the thickness of this film being about a nanometre.

Surface Energy and Surface Proton Order of Ice Ih

Ice Ih is comprised of orientationally disordered water molecules giving rise to positional disorder of the hydrogen atoms in the hydrogen bonded network of the lattice. Here we arrive at a first principles determination of the surface energy of ice Ih and suggest that the surface of ice is significantly more proton ordered than the bulk. We predict that the proton order-disorder transition, which occurs in the bulk at similar to 72 K, will not occur at the surface at any temperature below surface melting. An order parameter which defines the surface energy of ice Ih surfaces is also identified.

Proton ordering energetics in ice phases

Results from first-principles calculations on the subtle energetics of proton ordering in ice phases are shown only to depend on the electrostatic components of the total energy. Proton ordered ice phases can therefore be predicted using electronic structure methods or a tailored potential model. However, analysis of the electron density reveals that high order multipole components, up to hexadecapole, are needed to adequately capture total energy differences between proton ordered and disordered phases. This suggests that current potential models may be unable to reproduce the position of proton ordered ice phases in the phase diagram without extensions to describe high order electrostatics. (c) 2006 Elsevier B.V. All rights reserved.

Lorentz TEM investigation of the switching behavior in artificial spin ice building blocks

Frustration – the inability to simultaneously satisfy all interactions – occurs in a wide range of systems including neural networks, water ice and magnetic systems. An example of the latter is the so called spin-ice in pyrochlore materials [1] which have attracted a lot of interest not least due to the emergence of magnetic monopole defects when the ‘ice rules’ governing the local ordering breaks down [2]. However it is not possible to directly measure the frustrated property – the direction of the magnetic moments – in such spin ice systems with current experimental techniques. This problem can be solved by instead studying artificial spin-ice systems where the molecular magnetic moments are replaced by nanoscale ferromagnetic islands [3-8]. Two different arrangements of the ferromagnetic islands have been shown to exhibit spin ice behaviour: a square lattice maintaining four moments at each vertex [3,8] and the Kagome lattice which has only three moments per vertex but equivalent interactions between them [4-7]. Magnetic monopole defects have been observed in both types of lattices [7-8]. One of the challenges when studying these artificial spin-ice systems is that it is difficult to arrive at the fully demagnetised ground-state [6-8].
Here we present a study of the switching behaviour of building blocks of the Kagome lattice influenced by the termination of the lattice. Ferromagnetic islands of nominal size 1000 nm by 100 nm were fabricated in five island blocks using electron-beam lithography and lift-off techniques of evaporated 18 nm Permalloy (Ni80Fe20) films. Each block consists of a central island with four arms terminated by a different number and placement of ‘injection pads’, see Figure 1. The islands are single domain and magnetised along their long axis. The structures were grown on a 50 nm thick electron transparent silicon nitride membrane to allow TEM observation, which was back-coated with a 5 nm film of Au to prevent charge build-up during the TEM experiments.
To study the switching behaviour the sample was subjected to a magnetic field strong enough to magnetise all the blocks in one direction, see Figure 1. Each block obeys the Kagome lattice ‘ice-rules’ of “2-in, 1-out” or “1-in, 2-out” in this fully magnetised state. Fresnel mode Lorentz TEM images of the sample were then recorded as a magnetic field of increasing magnitude was applied in the opposite direction. While the Fresnel mode is normally used to image magnetic domain structures [9] for these types of samples it is possible to deduce the direction of the magnetisation from the Lorentz contrast [5]. All images were recorded at the same over-focus judged to give good Lorentz contrast.
The magnetisation was found to switch at different magnitudes of the applied field for nominally identical blocks. However, trends could still be identified: all the blocks with any injection pads, regardless of placement and number, switched the direction of the magnetisation of their central island at significantly smaller magnitudes of the applied magnetic field than the blocks without injection pads. It can therefore be concluded that the addition of an injection pad lowers the energy barrier to switching the connected island, acting as a nucleation site for monopole defects. In these five island blocks the defects immediately propagate through to the other side, but in a larger lattice the monopoles could potentially become trapped at a vertex and observed [10].
References

[1] M J Harris et al, Phys Rev Lett 79 (1997) p.2554.
[2] C Castelnovo, R Moessner and S L Sondhi, Nature 451 (2008) p. 42.
[3] R F Wang et al, Nature 439 (2006) 303.
[4] M Tanaka et al, Phys Rev B 73 (2006) 052411.
[5] Y Qi, T Brintlinger and J Cumings, Phys Rev B 77 (2008) 094418.
[6] E Mengotti et al, Phys Rev B 78 (2008) 144402.
[7] S Ladak et al, Nature Phys 6 (2010) 359.
[8] C Phatak et al, Phys Rev B 83 (2011) 174431.
[9] J N Chapman, J Phys D 17 (1984) 623.
[10] The authors gratefully acknowledge funding from the EPSRC under grant number EP/D063329/1.

Coleoptera and microbe biomass in Antarctic Dry Valley paleosols adjacent to the Inland Ice: Implications for Mars

Bulk paleosol samples collected from a Middle to Early Miocene moraine in the New Mountain area of the Dry Valleys, Antarctica, yielded Coleoptera exoskeletons and occasional endoskeletons showing considerable diagenetic effects along with several species of bacteria, all lodged in a dry-frozen but salt-rich horizon at shallow depth to the land surface. The till is at the older end of a chronologic sequence of glacial deposits, thought to have been deposited before the transition from wet-based to cold-based ice (similar to 15 Ma), and hence, entirely weathered in contact with the subaerial atmosphere. It is possible, though not absolutely verifiable, that the skeletons date from this early stage of emplacement having undergone modifications whenever light snowmelt occurred or salt concentrations lowered the freezing temperature to maintain water as liquid. Correlation of the Coleoptera species with cultured bacteria in the sample and the likelihood of co-habitation with Beauveria bassiani found in two adjacent, although younger paleosols, leads to new questions about the antiquity of the Coleoptera and the source of N and glucose from chitinase derived from the insects. The skeletons in the 831 section may date close to the oldest preserved chitin (Oligocene) yet found on Earth. While harsh Martian conditions make it seemingly intolerable for complex, multicellular organisms such as insects to exist in the near-surface and subaerially, life within similar cold, dry paleosol microenvironments (Cryosols) of Antarctica point to life potential for the Red Planet, especially when considering the relatively diverse microbe (bacteria and fungi) population. (C) 2011 Elsevier Ltd. All rights reserved.

Ash from Changbaishan Millennium eruption recorded in Greenland ice: Implications for determining the eruption's timing and impact

Major volcanic eruptions can impact on global climate by injecting large quantities of aerosols and ash into the atmosphere that alter the radiative balance and chemical equilibrium of the stratosphere. The Millennium eruption of Tianchi (Paektu), China/North Korea, was one of the largest Late Holocene eruptions. Uncertainty about the precise timing of the eruption has hindered the recognition of its climate impact in palaeoclimate and historical records. Here we report the compelling identification of the eruption's volcanic signal in Greenland ice cores through the association of geochemically-characterized volcanic glass, represented in by bimodal populations that compare with proximal material from the source eruption. The eruption most probably occurred in the AD 940?s, seven years after the Eldgjá eruption on Iceland. We examine the eruption's potential for climate forcing using the sulfate records from the ice-cores and conclude that it was unlikely to have had a global or extra-regional impact.

Asynchrony between the Greenland ice core chronology and the radiocarbon calibration curve? Evidence from Irish bogs

A significant cold event, deduced from the Greenland ice cores, took place between 8200 and 8000 cal. BP. Modeling of the event suggests that higher northern latitudes would have also experienced considerable decreases in precipitation and that Ireland would have witnessed one of the greatest depressions. However, no well-dated proxy record exists from the British Isles to test the model results. Here we present independent evidence for a phase of major pine recruitment on Irish bogs at around 8150 cal. BP. Dendrochronological dating of subfossil trees from three sites reveal synchronicity in germination across the region, indicative of a regional forcing, and allows for high-precision radiocarbon based dating. The inner-rings of 40% of all samples from the north of Ireland dating to the period 8500-7500 cal. BP fall within a 25-yr window. The concurrent colonization of pine on peatland is interpreted as drier conditions in the region and provides the first substantive proxy data in support of a significant hydrological change in the north of Ireland accompanying the 8.2 ka event. The dating uncertainties associated with the Irish pine record and the Greenland Ice Core Chronology 2005 (GICC05) do not allow for any overlap between the two. Our results indicate that the discrepancy could be an artifact of dating inaccuracy, and support a similar claim by Lohne et al. (2013) for the Younger Dryas boundaries. If real, this asynchrony will most likely have affected interpretations of previous proxy alignments.

Long-term mass balance of perennial firn and ice in an Alpine cave (Austria):Constraints from radiocarbon-dated wood fragments

Hundsalm ice cave located at 1520 m altitude in a karst region of western Austria contains up to 7-m-thick deposits of snow, firn and congelation ice. Wood fragments exposed in the lower parts of an ice and firn wall were radiocarbon accelerator mass spectrometry (AMS) dated. Although the local stratigraphy is complex, the 19 individual dates - the largest currently available radiocarbon dataset for an Alpine ice cave - allow to place constraints on the accumulation and ablation history of the cave ice. Most of the cave was either ice free or contained only a small firn and ice body during the 'Roman Warm Period'; dates of three wood fragments mark the onset of firn and ice build-up in the 6th and 7th century ad. In the central part of the cave, the oldest samples date back to the 13th century and record ice growth coeval with the onset of the 'Little Ice Age'. The majority of the ice and firn deposit, albeit compromised by a disturbed stratigraphy, appears to have been formed during the subsequent centuries, supported by wood samples from the 15th to the 17th century. The oldest wood remains found so far inside the ice is from the end of the Bronze Age and implies that local relics of prehistoric ice may be preserved in this cave. The wood record from Hundsalm ice cave shows parallels to the Alpine glacier history of the last three millennia, for example, the lack of preserved wood remains during periods of known glacier minima, and underscores the potential of firn and ice in karst cavities as a long-term palaeoclimate archive, which has been degrading at an alarming rate in recent years. © The Author(s) 2013.

Irradiation of water ice by C+ ions in the cosmic environment

We present a first principles molecular dynamics (FPMD) study of the interaction of low energy, positively charged, carbon (C+) projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon ion at an initial energy of 11 eV and 1.7 eV with 30-molecule clusters have been investigated. Simulations indicate that the neutral isoformyl radical, COH, and carbon monoxide, CO, are the dominant products of these reactions. All these reactions are accompanied by the transfer of a proton from the reacting water molecule to the ice, where it forms a hydronium ion. We find that COH is formed either via a direct, "knock-out", mechanism following the impact of the C+ projectile upon a water molecule or by creation of a COH_2^+ intermediate. The direct mechanism is more prominent at higher energies. CO is generally produced following the dissociation of COH. More frequent production of the formyl radical, HCO, is observed here than in gas phase calculations. A less commonly occurring product is the dihydroxymethyl, CH(OH)_2, radical. Although a minor result, its existence gives an indication of the increasing chemical complexity which is possible in such heterogeneous environments.

Tephrochronology and the extended intimate (integration of ice‐core, marine and terrestrial records) event stratigraphy 8-128kab2k

The comparison of palaeoclimate records on their own independent timescales is central to the work of the INTIMATE (INTegrating Ice core, MArine and TErrestrial records) network. For the North Atlantic region, an event stratigraphy has been established from the high-precision Greenland ice-core records and the integrated GICC05 chronology. This stratotype provides a palaeoclimate signal to which the timing and nature of palaeoenvironmental change recorded in marine and terrestrial archives can be compared. To facilitate this wider comparison, without assuming synchroneity of climatic change/proxy response, INTIMATE has also focussed on the development of tools to achieve this. In particular the use of time-parallel marker horizons e.g. tephra layers (volcanic ash). Coupled with the recent temporal extension of the Greenland stratotype, as part of this special issue, we present an updated INTIMATE event stratigraphy highlighting key tephra horizons used for correlation across Europe and the North Atlantic. We discuss the advantages of such an approach, and the key challenges for the further integration of terrestrial palaeoenvironmental records with those from ice cores and the marine realm.