Multi-scale study of condensation in water jets using ellipsoidal-statistical Bhatnagar-Gross-Krook and molecular dynamics modeling

Properties of water and argon clusters developed in supersonic expansions

Using adiabatic molecular dynamics coupled with the fluid dynamics equations, we model nucleation in an expanding beam of water vapor and argon on a microsecond scale. The size distribution of clusters, their temperature, and pickup cross sections in dependence on velocity are investigated and compared to the geometric cross sections and the experiment. The clusters are warmer than the expanding gas because of the time scale of relaxation processes. We also suggest that their translational and rotational kinetic energies are modified due to evaporative cooling. The pickup cross sections determined for the final clusters using molecules of the same kind increase with decreasing velocity, still obeying the (a+bN1/3)2 law.

Molecular Dynamics of Heterogeneous Systems on GPUs and Their Application to Nucleation in Gas Expanding to a Vacuum

Expansion of water vapor through a small orifice to a vacuum produces liquid or frozen clusters which in the experiment serve as model particles for atmospheric aerosols. Yet, there are controversies about the shape of these clusters, suggesting that the nucleation process is not fully understood. Such questions can be answered by molecular dynamics simulations; however, they require microsecond-scale runs with thousands of molecules and accurate energy conservation. The available highly parallel codes typically utilize domain decomposition and are inefficient for heterogeneous systems as clusters in a dilute gas. In this work, we present an implementation of molecular dynamics on graphics processing units based on the Verlet list and apply it to several systems for which experimental data are available. We reproduce sufficiently sized clusters but not the experimentally observed clusters of irregular shape.

Large cluster ions: soft local probes and tools for organic and bio surfaces

Ionised cluster beams have been produced and employed for thin film deposition and surface processing for half a century. In the last two decades, kiloelectronvolt cluster ions have also proved to be outstanding for surface characterisation by secondary ion mass spectrometry (SIMS), because their sputter and ion yields are enhanced in a non-linear fashion with respect to monoatomic projectiles, with a resulting step change of sensitivity for analysis and imaging. In particular, large gas cluster ion beams, or GCIB, have now become a reference in organic surface and thin film analysis using SIMS and X-ray photoelectron spectroscopy (XPS). The reason is that they induce soft molecular desorption and offer the opportunity to conduct damageless depth-profiling and 3D molecular imaging of the most sensitive organic electronics and biological samples, with a nanoscale depth resolution. In line with these recent developments, the present review focuses on rather weakly-bound, light-element cluster ions, such as noble or other gas clusters, and water or alcohol nanodroplets (excluding clusters made of metals, inorganic salts or ionic liquids) and their interaction with surfaces (essentially, but not exclusively, organic). The scope of this article encompasses three aspects. The first one is the fundamentals of large cluster impacts with surfaces, using the wealth of information provided by molecular dynamics simulations and experimental observations. The second focus is on recent applications of large cluster ion beams in surface characterisation, including mass spectrometric analysis and 2D localisation of large molecules, molecular depth-profiling and 3D molecular imaging. Finally, the perspective explores cutting edge developments, involving (i) new types of clusters with a chemistry designed to enhance performance for mass spectrometry imaging, (ii) the use of cluster fragment ion backscattering to locally retrieve physical surface properties and (iii) the fabrication of new biosurface and thin film architectures, where large cluster ion beams are used as tools to transfer biomolecules in vacuo from a target reservoir to any collector substrate.

Multiscale modeling of damaged surface topology in a hypersonic boundary

In this work, we used molecular dynamics (MD) to perform trajectory simulations of ice-like argon and amorphous silica aggregates on atomically smooth highly ordered pyrolytic graphite (HOPG) and a comparatively rougher quartz surface. It was found that at all incidence velocities, the quartz surface was stickier than the HOPG surface. The sticking probabilities and elastic moduli obtained from MD were then used to model surface evolution at a micron length scale using kinetic Monte Carlo (kMC) simulations. Rules were derived to control the number of sites available for the process execution in kMC to accurately model erosion of HOPG by atomic oxygen (AO) attack and ice-nucleation on surfaces. It was observed that the effect of defects was to increase the material erosion rate, while that of aggregate nucleation was to lower it. Similarly, simulations were performed to study the effects of AO attack and N₂ adsorption-desorption on surface evolution and it was found that N₂ adsorption-desorption limits the surface available for erosion by AO attack.

Study of thermal conductivity of ice clusters after impact deposition on the silica surfaces using the ReaxFF reactive force field

During aircraft or spacecraft missions, ice accumulates on different parts of their surface elements. An important parameter affecting the ability to remove this ice from the surface is the heat transfer characteristics of the accumulated ice. The ice heat transfer is related to the process of ice formation and its density and internal structure. In this study we investigate the effects of the ice and silica structure and the ice cluster attachment mechanism to the silica surface on the thermal conductivity (TC) of the attached ice cluster using the ReaxFF reactive force field. The purpose of this study is to investigate the thermal transport in amorphous and crystalline ice after high-velocity deposition on the silica surfaces. A dual thermostat method has been applied for the calculation of TC values. The validity of this method has been verified by comparing the calculated values of TC for crystal and amorphous ice with available experimental values. Our calculations show that the TC values of both crystal and amorphous ice drop after deposition on the silica surfaces. This decrease in the TC is more significant for the ice deposition on suboxide silica surfaces. Furthermore, crystal ice shows higher TC values than amorphous ice after accumulation. However, when crystal ice impacts on the silica surface at 1 km s(-1) impact speed, the crystalline shape of the ice cluster is lost to a considerable level and the TC values obtained for the ice clusters in such cases are closer to amorphous ice TC values. We observed a decrease in the TC values when ionic species are added inside the ice clusters.