Mixed Ice Research Articles

Thermal behavior of astrophysical amorphous molecular ices.

Ice is a major component of astrophysical environments - from interstellar molecular clouds through protoplanetary disks to evolved solar systems. Ice and complex organic matter coexist in these environments as well, and it is thought primordial ice brought the molecules of life to Earth four billion years ago, which could have kickstarted the origin of life on Earth. To understand the journey of ice and organics from their origins to becoming a part of evolved planetary systems, it is important to complement high spatial and spectral resolution telescopes such as JWST with laboratory experimental studies that provide deeper insight into the processes that occur in these astrophysical environments. Our laboratory studies are aimed at providing this knowledge. In this article we present simultaneous mass spectrometric and infrared spectroscopic investigation on how molecular ice mixtures behave at different temperatures and how this information is critical to interpret observational data from protoplanetary disks as well as comets. We find that amorphous to crystalline water ice transformation is the most critical phenomenon that differentiates between outgassing of trapped volatiles such as CO2vs. outgassing of pure molecular ice domains of the same in a mixed molecular ice. Crystalline water ice is found to trap only a small fraction of other volatiles (<5%), indicating ice grain composition in astrophysical and planetary environments must be different depending on whether the ice is in amorphous phase or transformed into crystalline phase, even if the crystalline ice undergoes radiation-induced amorphization subsequently. Crystallization of water ice is a key differentiator for many ices in astronomical environments as well as in our Solar System.

Ammonia, carbon dioxide, and the non-detection of the 2152 cm−1 CO band

Context. COa title should be crisp and newspaper-style, and direct questions are to be avoided in scientific writing. Please rephrase your title and take this into account. Your running title may be a good place to start is one of the most abundant ice components on interstellar dust grains. When it is mixed with amorphous solid water (ASW) or located on its surface, an absorption band of CO at 2152 cm−1 is always present in laboratory measurements. This spectral feature is attributed to the interaction of CO with dangling-OH bonds (dOH) in ASW. However, this band is absent in observational spectra of interstellar ices. This raises the question whether CO forms a relatively pure layer on top of ASW or is in close contact with ASW, but not via dangling bonds. Aims. We aim to determine whether the incorporation of NH3 or CO2 into ASW blocks the dOH and therefore reduces the 2152 cm−1 band. Methods. We performed laboratory experiments to simulate the layered structure of the ice mantle, that is, we grew CO ice on top of 1) pure ASW, 2) NH3:H2O=10:100 mixed ice, and 3) CO2:H2O=20:100 mixed ice. Infrared spectra were measured to quantify the strength of the 2152 cm−1 band. In addition, a second set of experiments were performed to determine how the incorporation of NH3 into ASW affects the dOH band. Results. We found that annealing the ice reduces the 2152 cm−1 band and that NH3 blocks the dOH on ASW surface and therefore reduces the 2152 cm−1 band more effectively than CO2. We suggest that this difference between NH3 and CO2 can be ascribed to the polarity of the guest molecule (NH3 is a polar species, whereas CO2 is apolar). The polarity implies that the formation of an H-bond between the N atom of ammonia and the dOH is a barrier-less reaction. We also determined the pore surface area of the ice mixtures as a function of the annealing temperature, and found that the nondetection of 2152 cm−1 band does not necessarily exclude the possibility of a porous ice mantle.

Laboratory spectroscopy of theoretical ices: Predictions for JWST and test for astrochemical models

Context. The pre-stellar core L1544 has been the subject of several observations conducted in the past years, complemented by modelling studies focused on its gas and ice-grain chemistry. The chemical composition of the ice mantles reflects the environmental physical changes along the temporal evolution, such as density and temperature. The investigation outcome hints at a layered structure of interstellar ices with abundance of H2O in the inner layers and an increasing concentration of CO near the surface. The morphology of interstellar ice analogues can be investigated experimentally assuming a composition derived from chemical models. Aims. This research presents a new approach of a three-dimensional fit where observational results are first fitted with a gas-grain chemical model predicting the exact ice composition including infrared (IR) inactive species. Then the laboratory IR spectra are recorded for interstellar ice analogues whose compositions reflect the obtained numerical results, in a layered and in a mixed morphology. These results could then be compared with the results of James Webb Space Telescope (JWST) observations. Special attention is paid to the inclusion of IR inactive species whose presence is predicted in the ice, but is typically omitted in the laboratory obtained data. This stands for N2, one of the main possible constituents of interstellar ice mantles, and O2. Methods. Ice analogue spectra were recorded at a temperature of 10 K using a Fourier transform infrared spectrometer. In the case of layered ice we deposited a H2O-CO-N2-O2 mixture on top of a H2O-CH3OH-N2 ice, while in the case of mixed ice we examined a H2O-CH3OH-N2-CO composition. The selected species are the four most abundant ice components predicted by the chemical model. Results. Following the changing composition and structure of the ice, we find differences in the absorption bands for most of the examined vibrational modes. The extent of observed changes in the IR band profiles will allow us to analyse the structure of ice mantles in L1544 from future observations by the JWST. Conclusions. Our spectroscopic measurements of interstellar ice analogues predicted by our well-received gas-grain chemical codes of pre-stellar cores will allow detailed comparison with upcoming JWST observations. This is crucial in order to put stringent constraints on the chemical and physical structure of dust icy mantles just before the formation of stars and protoplanetary disks, and to explain surface chemistry.

Atmospheric icing measurement and online ice type recognition for aircraft utilizing optical fiber sensor and machine learning algorithms

For the detection of atmospheric icing process, an ice type recognition and thickness measurement method based on optical fiber sensor is proposed and evaluated. Finite element models based on different structures are built for the four corresponding typical ice types in the atmospheric icing environment. Finite element analysis results show the characteristics of sensor response for different ice types and the potential of real-time ice type recognition using optical fiber sensor. In the laboratory test, under the atmospheric icing conditions, 1655 sets of valid data including glazed ice, mixed ice, rime ice and frost are collected. Cascade models based on machine learning algorithms are built and evaluated using a round-robin leave-one-out method. The model realizes instant ice type recognition and thickness measurement without relying on time-domain curve information. The evaluation results show that the proposed method achieves 92.8 % ice type recognition accuracy and 0.27 mm RMSE for ice thickness measurement.

On modelling cosmic ray sputtering of interstellar grain ices

ABSTRACT In the interstellar medium (ISM), the formation of complex organic molecules (COMs) is largely facilitated by surface reactions. However, in cold dark clouds, thermal desorption of COMs is inefficient because of the lack of thermal energy to overcome binding energies to the grain surface. Non-thermal desorption methods are therefore important explanations for the gas-phase detection of many COMs that are primarily formed on grains. Here, we present a new non-thermal desorption process: cosmic ray sputtering of grain ice surfaces based on water, carbon dioxide, and a simple mixed ice. Our model applies estimated rates of sputtering to the three-phase rate equation model nautilus-1.1, where this inclusion results in enhanced gas-phase abundances for molecules produced by grain reactions such as methanol (CH3OH) and methyl formate (HCOOCH3). Notably, species with efficient gas-phase destruction pathways exhibit less of an increase in models with sputtering compared to other molecules. These model results suggest that sputtering is an efficient, non-specific method of non-thermal desorption that should be considered as an important factor in future chemical models.

Accurate crystal structure of ice VI from X-ray diffraction with Hirshfeld atom refinement.

Water is an essential chemical compound for living organisms, and twenty of its different crystal solid forms (ices) are known. Still, there are many fundamental problems with these structures such as establishing the correct positions and thermal motions of hydrogen atoms. The list of ice structures is not yet complete as DFT calculations have suggested the existence of additional and - to date - unknown phases. In many ice structures, neither neutron diffraction nor DFT calculations nor X-ray diffraction methods can easily solve the problem of hydrogen atom disorder or accurately determine their anisotropic displacement parameters (ADPs). Here, accurate crystal structures of H2O, D2O and mixed (50%H2O/50%D2O) ice VI obtained by Hirshfeld atom refinement (HAR) of high-pressure single-crystal synchrotron and laboratory X-ray diffraction data are presented. It was possible to obtain O-H/D bond lengths and ADPs for disordered hydrogen atoms which are in good agreement with the corresponding single-crystal neutron diffraction data. These results show that HAR combined with X-ray diffraction can compete with neutron diffraction in detailed studies of polymorphic forms of ice and crystals of other hydrogen-rich compounds. As neutron diffraction is relatively expensive, requires larger crystals which can be difficult to obtain and access to neutron facilities is restricted, cheaper and more accessible X-ray measurements combined with HAR can facilitate the verification of the existing ice polymorphs and the quest for new ones.

Design of Capacitance and Impedance Dual-Parameters Planar Electrode Sensor for Thin Ice Detection of Aircraft Wings

An ice types recognition and ice thickness detection system was designed based on capacitance and impedance double parameters. By studying the influence of electrode structure on sensor performance based on finite element simulation, the finger electrode was selected with better performance, and a planar electrode sensor was designed suitable for wing thin ice detection. After the thin ice (2.5 mm) making experiment, the capacitance and impedance spectra of three ice types (glaze ice, rime ice and mixed ice) were obtained. Among them, capacitance can reflect the growth of ice well, and the impedance spectra of different ice types are significantly different. Therefore, a capacitance-based regression model and an impedance-based classification model were built. The results of regression model show that about 90% of the calculation error of thin ice thickness is concentrated in the range of [−0.2 mm, 0.2 mm]. Finally, a field test was conducted in the ice wind tunnel to preliminarily determine the system effectiveness.

Titik Kritis Keamanan Pangan pada Tahap Pengolahan dan Penyajian Beberapa Jenis Minuman Es

In small vendors, iced beverages potentially contain microbiological and chemical hazards. The objective of this study was to identify critical points of hazard during processing and serving steps various types of most popular iced beverages. Based on processing step, there were identified nine types of iced beverages i.e. coconut, grass jelly, sliced fruit, orange, mixed, cendol, juice, blended and flavouring ice. The study was conducted by observation 90 vendors addressed to processing step, utilities, and environ-mental condition at low, middle and high economic scale of vendors. The critical points were determined according to the principles of HACCP. The results showed that more complex of processing step, lower economic scale and poor of facilities caused more critical points. Mixed ice (containing fruits, grass jelly and jelly) that sold at lower economic scale vendor has the most critical points number (11 out of 14), i.e. ice cube, water mixture, syrup, ice washing, ice downsizing, ice storage, fruits and ingredients washing process, fruit and jelly slicing process and storage, and mixing process. On the medium economic scale vendor, there were 1-3 critical points number and on the highest economic scale there was only one critical point number which was in the mixing process.

Investigation on one-step preparation and anti-icing experiments of robust super-hydrophobic surface on wind turbine blades

Wind turbines are often built in areas prone to icing in winter, such as high altitudes and humidity. When wind turbine blades rotate, they collide and capture super-cooled water droplets to form ice, which causes significant hidden dangers to the safe operation of wind turbines. The super-hydrophobic surface is an attractive topic in the field of anti-icing research. However, the low weather resistance of the existing super-hydrophobic surface limits its large-scale application. This paper proposes a one-step preparation method for a robust super-hydrophobic surface. Its hydrophobicity, wear resistance, and other properties are tested. We have experimented with super-hydrophobic blades and control blades in a natural icing environment. The results show that the prepared super-hydrophobic blade has better wear resistance, acid, and alkali resistance. The super-hydrophobic surface has a retardation effect on glaze and mixed ice only in the initial icing stage. The robust super-hydrophobicity limits the frozen area. Throughout the icing period, the effectiveness of super-hydrophobic surfaces is limited. Although the prepared surface effectively improves the robustness, the super-hydrophobic surface is still very challenging in anti-icing applications.

Mapping liquid water content in snow at the millimeter scale: an intercomparison of mixed-phase optical property models using hyperspectral imaging and in situ measurements

Abstract. It is well understood that the distribution and quantity of liquid water in snow is relevant for snow hydrology and avalanche forecasting, yet detecting and quantifying liquid water in snow remains a challenge from the micro- to the macro-scale. Using near-infrared (NIR) spectral reflectance measurements, previous case studies have demonstrated the capability to retrieve surface liquid water content (LWC) of wet snow by leveraging shifts in the complex refractive index between ice and water. However, different models to represent mixed-phase optical properties have been proposed, including (1) internally mixed ice and water spheres, (2) internally mixed water-coated ice spheres, and (3) externally mixed interstitial ice and water spheres. Here, from within a controlled laboratory environment, we determined the optimal mixed-phase optical property model for simulating wet snow reflectance using a combination of NIR hyperspectral imaging, radiative transfer simulations (Discrete Ordinate Radiative Transfer model, DISORT), and an independent dielectric LWC measurement (SLF Snow Sensor). Maps of LWC were produced by finding the lowest residual between measured reflectance and simulated reflectance in spectral libraries, generated for each model with varying LWC and grain size, and assessed against the in situ LWC sensor. Our results show that the externally mixed model performed the best, retrieving LWC with an uncertainty of ∼1 %, while the simultaneously retrieved grain size better represented wet snow relative to the established scaled band area method. Furthermore, the LWC retrieval method was demonstrated in the field by imaging a snowpit sidewall during melt conditions and mapping LWC distribution in unprecedented detail, allowing for visualization of pooling water and flow features.

3D Nanoprinting by Electron-Beam with an Ice Resist

Following the general trend in the miniaturization of electronic devices, techniques that enable 3D printing at the nanometer scale are gaining momentum. As a widely used planar processing method, electron-beam lithography (EBL) can be employed to create 3D nanostructures in a layer-by-layer fashion. However, compared with other 3D printing techniques, EBL is limited by the stringent requirement of a range of fabrication equipment and complex fabrication processes. Here, we have demonstrated that EBL can be developed to a controllable 3D nanoprinting technology with the aid of ice resists. With carefully selected accelerating voltage, electron dose, and ice thickness, 3D objects can be efficiently printed in a single vacuum system through an iterative process of ice deposition and e-beam exposure. Mixed ice resists containing solid anisole and water are also introduced into the printing process, which offer a flexible control of the thickness of printed layers. Apart from carbonaceous objects obtained with our method, 3D printing of metals is also promising by employing organometallic compounds as ice resists. This study provides a fresh perspective in EBL-based nanofabrication and expands the spectrum of modern additive manufacturing.

New insights in glaciers characterization by differential diagnosis integrating GPR and remote sensing techniques: A case study for the Eastern Gran Zebrù glacier (Central Alps)

In GPR profiles, ice is usually imaged as a mostly electromagnetic transparent facies. However, diffraction events, as well as internal layering, can be also observed. In some cases, the bedrock below glaciers is masked by dense diffractions usually interpreted as the effect of liquid water pockets inside the so-called warm ice. However, the actual physical meaning of such GPR facies is not always obvious, because it can be related also to mixed debris and ice deposits.We adopted a strategy well known in medical sciences and referred as “differential diagnosis” in order to infer which is the actual meaning of a high scattering facies imaged within the Eastern Gran Zebrù glacier (Central Italian Alps) and, more generally, of all the internal glacier features. In fact, in many cases, there is no direct information to limit the subjectivity of geophysical interpretation; therefore, we provide all the discriminative hypotheses based on both independent and integrated criteria including GPR attribute analysis, imaging effects, reflection analysis, GPR frequency evaluations combined with geomorphological and remote sensing data obtained by two photogrammetric UAV and thermal infrared surveys.On the basis of the differential diagnosis, we concluded that the high scattering zone embedded within the studied glacier is most likely related to a mixture of ice and debris probably formed during a past shrinking phase. Beside this case study, this approach could be helpful in other GPR glaciological surveys, in which the target is related not only to the bedrock detection, but also to a detailed analysis of the internal facies of a glacier.

A Wind Tunnel Experimental Study of Icing on NACA0012 Aircraft Airfoil with Silicon Compounds Modified Polyurethane Coatings.

Ice formation on the aerodynamic surfaces of an aircraft is regarded as a major problem in the aerospace industry. Ice accumulation may damage parts, sensors and controllers and alter the aerodynamics of the airplane, leading to a range of undesired consequences, including flight delays, emergency landings, damaged parts and increased energy consumption. There are various approaches to reducing ice accretion, one of them being the application of icephobic coatings. In this work, commercially available polyurethane-based coatings were modified and deposited on NACA 0012 aircraft airfoils. A hybrid modification of polyurethane (PUR) topcoats was adopted by the addition of nanosilica and three-functional spherosilicates (a variety of silsesqioxane compound), which owe their unique properties to the presence of three different groups. The ice accretion on the manufactured nanocomposites was determined in an icing wind tunnel. The tests were performed under three different icing conditions: glaze ice, rime ice and mixed ice. Furthermore, the surface topography and wetting behavior (static contact angle and contact angle hysteresis) were investigated. It was found that the anti-icing properties of polyurethane nanocomposite coatings strongly depend on the icing conditions under which they are tested. Moreover, the addition of nanosilica and spherosilicates enabled the reduction of accreted ice by 65% in comparison to the neat topcoat.

Spectroscopic measurements of CH3OH in layered and mixed interstellar ice analogues

Context. The molecular composition of interstellar ice mantles is defined by gas-grain processes in molecular clouds, with the main components being H2O, CO, and CO2. Methanol (CH3OH) ice is detected towards the denser pre-stellar cores and star-forming regions, where large amounts of CO molecules freeze out and get hydrogenated on top of the icy grains. The thermal heating from nearby protostars can further change the ice structure and composition. Despite the several observations of icy features carried out towards molecular clouds and along the line of site of protostars, it is not yet clear if interstellar ices are mixed or if they have a layered structure. Aims. We aim to examine the effect of mixed and layered ice growth in dust grain mantle analogues, with specific focus on the position and shape of methanol infrared bands, so dedicated future observations could shed light on the structure of interstellar ices in different environments. Methods. Mixed and layered ice samples were deposited on a cold substrate kept at a temperature of 10 K using a closed-cycle cryostat placed in a vacuum chamber. The spectroscopic features were analysed by Fourier transform infrared spectroscopy. Different proportions of the most abundant four molecular species in ice mantles, namely H2O, CO, CO2, and CH3OH, were investigated, with a special attention placed on the analysis of the CH3OH bands. Results. We measure changes in the position and shape of the CH and CO stretching bands of CH3OH depending on the mixed or layered nature of the ice sample. Spectroscopic features of methanol are also found to change due to heating. Conclusions. A layered ice structure best reproduces the CH3OH band position recently observed towards a pre-stellar core and in star-forming regions. Based on our experimental results, we conclude that observations of CH3OH ice features in space can provide information about the structure of interstellar ices, and we expect the James Webb Space Telescope to put stringent constraints on the layered or mixed nature of ices in different interstellar environments, from molecular clouds to pre-stellar cores to protostars and protoplanetary discs.

PENERAPAN METODE COST REDUCTION TERHADAP BAHAN BAKU PEDAGANG KULINER MINUMAN ES TELER

Es teler is an authentic Indonesian fruit cocktail. Avocado, young coconut, grass jelly, jackfruit and other fruits are served with coconut milk, sweetened condensed milk, Pandanus amaryllifolius leaves (usually in the form of coconut and pandan flavored syrup), sugar, and a pinch of salt. Es teler is an original drink from Central Java, Indonesia. The origin or history turns out to be from Sukoharjo Regency. In 1957, this drink appeared because it was discovered by Hj. Samijem Darmo Putro (Mother of Tukiman Darmowijono), a resident of Juron village, Nguter District, Sukoharjo Regency, Central Java. Ice telers are culinary foods that can be found in almost all regions of Indonesia. This study discusses minimizing expenses in making ice telers by using a cost reduction strategy, namely by reducing unnecessary or important costs incurred for making teler ice. Es teler sales business includes micro, small and medium enterprises (MSMEs) which can be carried out by anyone who is able to develop it. This business can be started at medium to lower costs without using high capital. High sales and small expenses can keep the business afloat. The conclusion obtained from this research is that the price of the raw materials used in Es Teler durian bang dens which is located in Karawang and Es Campur 88 which is located in Purwakarta has differences in terms of price and quality of these raw materials. Es Teler durian bang dens does not have salaried employees, while Es Campur 88 has employees who need to be paid so that in mixed ice 88 they spend more.

UAV icing: the influence of airspeed and chord length on performance degradation

PurposeThe main purpose of this paper is to investigate the effects of icing on unmanned aerial vehicles (UAVs) at low Reynolds numbers and to highlight the differences to icing on manned aircraft at high Reynolds numbers. This paper follows existing research on low Reynolds number effects on ice accretion. This study extends the focus to how variations of airspeed and chord length affect the ice accretions, and aerodynamic performance degradation is investigated.Design/methodology/approachA parametric study with independent variations of airspeed and chord lengths was conducted on a typical UAV airfoil (RG-15) using icing computational fluid dynamic methods. FENSAP-ICE was used to simulate ice shapes and aerodynamic performance penalties. Validation was performed with two experimental ice shapes obtained from a low-speed icing wind tunnel. Three meteorological conditions were chosen to represent the icing typologies of rime, glaze and mixed ice. A parameter study with different chord lengths and airspeeds was then conducted for rime, glaze and mixed icing conditions.FindingsThe simulation results showed that the effect of airspeed variation depended on the ice accretion regime. For rime, it led to a minor increase in ice accretion. For mixed and glaze, the impact on ice geometry and penalties was substantially larger. The variation of chord length had a substantial impact on relative ice thicknesses, ice area, ice limits and performance degradation, independent from the icing regime.Research limitations/implicationsThe implications of this manuscript are relevant for highlighting the differences between icing on manned and unmanned aircraft. Unmanned aircraft are typically smaller and fly slower than manned aircraft. Although previous research has documented the influence of this on the ice accretions, this paper investigates the effect on aerodynamic performance degradation. The findings in this work show that UAVs are more sensitive to icing conditions compared to larger and faster manned aircraft. By consequence, icing conditions are more severe for UAVs.Practical implicationsAtmospheric in-flight icing is a severe risk for fixed-wing UAVs and significantly limits their operational envelope. As UAVs are typically smaller and operate at lower airspeeds compared to manned aircraft, it is important to understand how the differences in airspeed and size affect ice accretion and aerodynamic performance penalties.Originality/valueEarlier work has described the effect of Reynolds number variations on the ice accretion characteristics for UAVs. This work is expanding on those findings by investigating the effect of airspeed and chord length on ice accretion shapes separately. In addition, this study also investigates how these parameters affect aerodynamic performance penalties (lift, drag and stall).

A New Ice Accretion Model for Aircraft Icing Based on Phase-Field Method

Aircraft icing presents a serious threat to the aerodynamic performance and safety of aircraft. The numerical simulation method for the accurate prediction of icing shape is an important method to evaluate icing hazards and develop aircraft icing protection systems. Referring to the phase-field method, a new ice accretion mathematical model is developed to predict the ice shape. The mass fraction of ice in the mixture is selected as the phase parameter, and the phase equation is established with a freezing coefficient. Meanwhile, the mixture thickness and temperature are determined by combining mass conservation and energy balance. Ice accretions are simulated under typical ice conditions, including rime ice, glaze ice and mixed ice, and the ice shape and its characteristics are analyzed and compared with those provided by experiments and LEWICE. The results show that the phase-field ice accretion model can predict the ice shape under different icing conditions, especially reflecting some main characteristics of glaze ice.

The relationship between water‐holding capacities of soybean–whey mixed protein and ice crystal size for ice cream

AbstractThis study investigated the water‐holding properties of a mixture of whey protein (WP) and soy protein (SP) (WP–SP) and their effect on ice crystal size using low‐field nuclear magnetic resonance (LF‐NMR) and freezing microscopy techniques. The size of ice crystals determines the perceived quality of ice cream. Furthermore, the addition of plant protein facilitates the study of the water‐holding properties of mixed proteins and limits the size of ice crystals. WP and SP proteins were mixed in different ratios (10:0–0:10), the water‐holding capacity was measured by LF‐NMR, and the size of the ice crystals observed by freezing microscopy. The results showed that the water‐holding capacity of the mixed proteins increased with increasing SP content. A mixed protein ratio of 3:7 had the highest water‐holding capacity but also displayed a strong inhibitory effect on the size of the ice crystals formed. Hence, the ice crystal area of SP solution alone was 128.96 μm2 with an average diameter of 11.74 μm, and, in contrast, a mixed protein ratio of 3:7 had an ice crystal area of 79.49 μm2, and an average diameter of 8.61 μm.Practical ApplicationsThe study provides insight into the design of an ice cream formula to assist with the industrial production of ice cream products of improved quality and consumer appeal.

Liquid Water and Interfacial, Cubic, and Hexagonal Ice Classification through Eclipsed and Staggered Conformation Template Matching.

We propose a novel method based on template matching for the recognition of liquid water, cubic ice (ice Ic), hexagonal ice (ice Ih), clathrate hydrates, and different interfacial structures in atomistic and coarse-grained simulations of water and ice. The two template matrices represent staggered and eclipsed conformations, which are the building blocks of hexagonal and cubic ice and clathrate crystals. The algorithm is rotationally invariant and highly robust against imperfections in the ice structure, and its sensitivity for recognizing ice-like structures can be tuned for different applications. Unlike most other algorithms, it can discriminate between cubic, hexagonal, clathrate, mixed, and other interfacial ice types and is therefore well suited to study complex systems and heterogeneous ice nucleation.

Dwarf planet (1) Ceres surface bluing due to high porosity resulting from sublimation

The Dawn mission found that the dominant colour variation on the surface of dwarf planet Ceres is a change of the visible spectral slope, where fresh impact craters are surrounded by blue (negative spectral-sloped) ejecta. The origin of this colour variation is still a mystery. Here we investigate a scenario in which an impact mixes the phyllosilicates present on the surface of Ceres with the water ice just below. In our experiment, Ceres analogue material is suspended in liquid water to create intimately mixed ice particles, which are sublimated under conditions approximating those on Ceres. The sublimation residue has a highly porous, foam-like structure made of phyllosilicates that scattered light in similar blue fashion as the Ceres surface. Our experiment provides a mechanism for the blue colour of fresh craters that can naturally emerge from the Ceres environment.