Observation of surface layering in a nonmetallic liquid

Nonmonotonic Dynamical Correlations beneath the Surface of Glass-Forming Liquids

Collective motion over increasing length scales is a signature of the vitrification process of liquids. We demonstrate how distinct static and dynamic length scales govern the dynamics of vitrifying films. In contrast to a monotonically growing static correlation length, the dynamical correlation length that measures the extent of surface-dynamics acceleration into the bulk displays a striking nonmonotonic temperature evolution that is robust also against changes in detailed interatomic interaction. This nonmonotonic change defines a crossover temperature T_{*} that is distinct from the critical temperature T_{c} of mode-coupling theory. We connect this nonmonotonic change to a morphological change of cooperative rearrangement regions of fast particles, and to the point where the decoupling of fast-particle motion from the bulk relaxation is most sensitive to fluctuations. We propose a rigorous definition of this new crossover temperature T_{*} within a recent extension of mode-coupling theory, the stochastic β-relaxation theory.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Weakly Anisotropic Dielectric Properties of Water Droplets at the Nanoscale

The slab-confined water at the nanoscale exhibits anomalous dielectric properties compared to bulk water, for example, significantly low dielectric constant. In this work, we study the dielectric properties of nanoscale water droplets at room temperature using molecular dynamics simulations. We find that the nanoscale water droplets feature weakly anisotropic dielectric constant: the radial component of dielectric constants is distinctly smaller than the tangential component although they both decrease with reducing droplet size in a similar way. Such dielectric behavior is closely related to the orientational preference of water molecules near the convex surface. The molecular dipole prefers to slightly orientate toward the interior of droplets in contrast to the out-of-plane preference for free-standing water films and slab-confined water, which suppresses the fluctuation of dipole moments in the radial direction. Meanwhile, it facilitates the formation of the open hydrogen-bond network in the surface layer and ultimately leads to the relatively weak suppression of tangential fluctuations. The differential suppression is responsible for the anisotropic dielectric constant of water droplets. This anisotropic characteristic is also found in dielectric relaxation: both the radial and the tangential relaxation are consistently slowed down upon approaching surface but the latter is universally slower.

Comment on "Bi-layering at ionic liquid surfaces: a sum - frequency generation vibrational spectroscopy - and molecular dynamics simulation-based study" by T. Iwahashi, T. Ishiyama, Y. Sakai, A. Morita, D. Kim and Y. Ouchi, Phys. Chem. Chem. Phys., 2020, 22, 12565

This Comment raises several questions concerning the surface structure concluded in the paper referenced in the title. Specifically, that paper ignores previous experiments and simulations which demonstrate for the same ionic liquids depth-decaying, multilayered surface-normal density profiles rather than the claimed molecular mono- or bi-layers. We demonstrate that the claimed structure does not reproduce the measured X-ray reflectivity, which probes directly the surface-normal density profile. The measured reflectivities are found, however, to be well-reproduced by a multilayered density model. These results, and previous experimental and simulation results, cast severe doubt on the validity of the surface structure claimed in the paper referenced in the title.

Crystallization of highly supercooled glass-forming alloys induced by anomalous surface wetting

Crystallization in highly supercooled Cu50Zr50 films close to the glass transition is studied by using molecular dynamics simulations. Spontaneous nucleation is observed at the simulation timescale in contrast to the bulk counterpart. We find that nucleation occurs at free surfaces owing to the partial wetting of the nucleus by melt. The anomalous wetting phenomenon is closely related to strong density layering arising from the surface: the high density associated with surface layering increases surface energy of supercooled melts, resulting in that one facet of the crystalline embryo is preferentially formed on the film surface. The surface-based embryo is then developed into a stable nucleus by bridging two surfaces of thin films. The kinetics and thermodynamics analyses based on the mean first-passage time method show that the nucleation process still follows the description of the classical nucleation theory despite extremely high supercoolings. In nucleating, the slow interface dynamics becomes dominant and induces a low nucleation rate although the nucleation barrier is very low. The subsequent crystal growth is found to proceed in a quasi-two-dimensional manner with a ramified interface morphology, which is analogous to percolative crystals predicted in glass-forming liquids.

Multiscale Relaxation Dynamics in Ultrathin Metallic Glass-Forming Films

The density layering phenomenon originating from a free surface gives rise to the layerlike dynamics and stress heterogeneity in ultrathin Cu-Zr glassy films, which facilitates the occurrence of multistep relaxations in the timescale of computer simulations. Taking advantage of this condition, we trace the relaxation decoupling and evolution with temperature simply via the intermediate scattering function. We show that the β relaxation hierarchically follows fast and slow modes in films, and there is a β-relaxation transition as the film is cooled close to the glass transition. We provide the direct observation of particle motions responsible for the β relaxation and reveal the dominant mechanism varying from the thermal activated to the cooperative jumps across the transition.

Effect of Concentration on the Interfacial and Bulk Structure of Ionic Liquids in Aqueous Solution

Bio and aqueous applications of ionic liquids (IL) such as catalysis in micelles formed in aqueous IL solutions or extraction of chemicals from biologic materials rely on surface-active and self-assembly properties of ILs. Here, we discuss qualitative relations of the interfacial and bulk structuring of a water-soluble surface-active IL ([C₈MIm][Cl]) on chemically controlled surfaces over a wide range of water concentrations using both force probe and X-ray scattering experiments. Our data indicate that IL structuring evolves from surfactant-like surface adsorption at low IL concentrations, to micellar bulk structure adsorption above the critical micelle concentration, to planar bilayer formation in ILs with <1 wt % of water and at high charging of the surface. Interfacial structuring is controlled by mesoscopic bulk structuring at high water concentrations. Surface chemistry and surface charges decisively steer interfacial ordering of ions if the water concentration is low and/or the surface charge is high. We also demonstrate that controlling the interfacial forces by using self-assembled monolayer chemistry allows tuning of interfacial structures. Both the ratio of the head group size to the hydrophobic tail volume as well as the surface charging trigger the bulk structure and offer a tool for predicting interfacial structures. Based on the applied techniques and analyses, a qualitative prediction of molecular layering of ILs in aqueous systems is possible.

Surface structure evolution in a homologous series of ionic liquids

Interfaces of room temperature ionic liquids (RTILs) are important for both applications and basic science and are therefore intensely studied. However, the evolution of their interface structure with the cation's alkyl chain length [Formula: see text] from Coulomb to van der Waals interaction domination has not yet been studied for even a single broad homologous RTIL series. We present here such a study of the liquid-air interface for [Formula: see text], using angstrom-resolution X-ray methods. For [Formula: see text], a typical "simple liquid" monotonic surface-normal electron density profile [Formula: see text] is obtained, like those of water and organic solvents. For [Formula: see text], increasingly more pronounced nanoscale self-segregation of the molecules' charged moieties and apolar chains yields surface layering with alternating regions of headgroups and chains. The layering decays into the bulk over a few, to a few tens, of nanometers. The layering periods and decay lengths, their linear [Formula: see text] dependence, and slopes are discussed within two models, one with partial-chain interdigitation and the other with liquid-like chains. No surface-parallel long-range order is found within the surface layer. For [Formula: see text], a different surface phase is observed above melting. Our results also impact general liquid-phase issues like supramolecular self-aggregation and bulk-surface structure relations.

Self-segregated nanostructure in room temperature ionic liquids

The nanosegregated bulk structure, and its evolution with the cation's alkyl length n, are studied by X-ray scattering for an unprecedentedly broad homologous series of a model room-temperature ionic liquid, [C_nMIM][NTf₂] (n = 4-22). A tri-periodic local structure is found, with the lateral periodicities, d_II and d_III independent of n, and a longitudinal one, d_I, linearly increasing with n. The results are consistent with a local structure comprising alternating layers of polar headgroups and apolar, interdigitated, partly overlapping, cations' alkyl tails, of an average macroscopic mass density close to that of liquid alkanes. A slope decrease in the linear d_I(n) suggests a change from a lower to a higher rate of increase with n of chain overlap for n ≥ 12. The order decay lengths of the layering, and of the lateral chain packing, increase with n, as expected from the increasing van der Waals interaction's domination of the structure. The headgroups' lateral packing decay length decreases with n, due to increasing frustration between the longer lateral periodicity preferred by the headgroups, and the shorter lateral periodicity preferred by the chains. A comparison of the bulk and surface structures highlights the surface's ordering effect, which, however, does not induce here a surface phase different from the bulk, as it does in liquid crystals and liquid alkanes.

Influence of van der Waals interactions on morphology and dynamics in ultrathin liquid films at silicon oxide interfaces

Single molecule tracer diffusion studies of evaporating (thinning) ultrathin tetrakis-2-ethyl-hexoxysilane (TEHOS) films on silicon with 100 nm thermal oxide reveal a considerable slowdown of the molecular mobility within less than 4 nm above the substrate (corresponding to a few molecular TEHOS layers). This is related to restricted mobility and structure formation of the liquid in this region, in agreement with information obtained from a long-time ellipsometric study of thinning TEHOS films on silicon substrates with 100 nm thermal or 2 nm native oxide. Both show evidence for the formation of up to four layers. Additionally, on thermal oxide, a lateral flow of the liquid is observed, while the film on the native oxide forms an almost flat surface and shows negligible flow. Thus, on the 2 nm native oxide the liquid mobility is even more restricted in close vicinity to the substrate as compared to the 100 nm thermal oxide. In addition, we found a significantly smaller initial film thickness in case of the native oxide under similar dipcoating conditions. We ascribe these differences to van der Waals interactions with the underlying silicon in case of the native oxide, whereas the thermal oxide suffices to shield those interactions.

Melting of monatomic glass with free surfaces

Melting of monatomic glass with free surfaces has been studied by molecular dynamics simulations in models with Lennard-Jones-Gauss interatomic potential. Models have been heated up from a glassy state toward a normal liquid state. Atomic mechanism of melting has been analyzed via monitoring spatio-temporal arrangements of liquid-like atoms occurred during heating process. Liquid-like atoms are detected via the Lindemann criterion of melting. It is clear that the transition from glass into supercooled liquid of our "ordinary" glass with free surfaces exhibits a non-heterogeneous behavior, i.e., although liquid-like atoms initiate/grow mainly in the surface shell, significant amount of liquid-like atoms also initiates/grows simultaneously in the interior during heating process. We found three characteristic temperatures of melting of glass with a free surface. Temperature dependence of structure and various thermodynamic quantities of the system upon heating is also presented and discussed.

Probing the protic ionic liquid surface using X-ray reflectivity

The structure of the free liquid surface of three protic ionic liquids, ethylammonium nitrate (EAN), propylammonium nitrate (PAN), and ethylammonium formate (EAF), has been elucidated using X-ray reflectivity. The results show all three liquids have an extended interfacial region, spanning at least five ion pairs, which can be divided into two parts. Adjacent to the gas phase are aggregates consisting of multiple cations and anions. Below this are layers oriented parallel to the macroscopic surface that are alternately enriched and depleted in cation alkyl chains and polar domains of cation ammonium groups and their anions, gradually decaying to the isotropic sponge-like bulk structure. The most pronounced layering is observed for PAN, driven by strong solvophobic interactions, while reduced hydrogen bonding in EAF results in the least structured and least extensive interfacial region.

Temperature dependence of multilayering at the free surface of ionic liquids probed by X-ray reflectivity measurements

The effect of the temperature on the surface layering of ionic liquids has been studied for two ionic liquids, trioctylmethylammonium bis(nonafluorobutanesulfonyl)amide([TOMA(+)][C(4)C(4)N(-)]) and trihexyltetradecylphosphonium bis(nonafluorobutanesulfonyl)amide ([THTDP(+)][C(4)C(4)N(-)]), using X-ray reflectivity measurements at 285, 300, and 315 K. Both [TOMA(+)][C(4)C(4)N(-)] and [THTDP(+)][C(4)C(4)N(-)] develop multilayers at the surface. The structure of the multilayers at the [TOMA(+)][C(4)C(4)N(-)] surface shows little temperature-dependent change, whereas that at the [THTDP(+)][C(4)C(4)N(-)] surface clearly becomes diffused with increasing temperature. The different temperature dependence seems to be related to the difference in the recently reported ultraslow dynamics of the interfacial structure of [TOMA(+)][C(4)C(4)N(-)] and [THTDP(+)][C(4)C(4)N(-)] at the ionic liquid|water interface.

Structural signal of a dynamic glass transition

Pentaphenyl trimethyl trisiloxane is an isotropic liquid at room temperature with a dynamic glass transition at 224 K. Using x-ray reflectivity, we see surface density oscillations (layers) develop below 285 K, similar to those seen in other metallic and dielectric liquids and in computer simulations. The layering threshold is approximately 0.23 times the liquid-gas critical temperature. Upon cooling further, there is a sharp increase at 224 K in the persistence of the surface layers into the bulk material, i.e., an apparently discontinuous change in static structure at the glass transition.

An inherent structure view of liquid-vapor interfaces

The formal statistical mechanical theory describing liquid-vapor interfaces at thermal equilibrium has been incomplete, owing partly to discrepancies between two primary views, the local free energy approach originated by van der Waals and the capillary wave approach initiated by Mandelstam. The former provides detailed prescriptions for interfacial density profiles and surface tensions, and has recently been tailored to conform to nonclassical critical phenomena. The latter postulates a crude discontinuous density profile that is geometrically delocalized by surface excitations, but which qualitatively incorporates basic gravitational effects that are missing in the van der Waals method. The present analysis provides a formalism within which both approaches can be reconciled. It draws upon the inherent structure mapping procedure to define an intrinsic liquid-vapor density profile which invokes no arbitrary parameters that are not already present in the many-particle potential energy function or its thermodynamics. By construction, this intrinsic profile plays a role conventionally given to a van der Waals interface profile, and although it is free of capillary waves it can serve as a starting point for evaluating the gravitational implications of those interfacial fluctuation effects.

Observation of a liquid-to-layered transition in thin liquid films when surface and interface regions overlap

We have used x-ray reflectivity to study the coupling of surface and interface layering in a molecularly thin normal liquid [tetrakis(2-ethylhexoxy)silane (TEHOS)], as a function of temperature and film thickness. The best fits to the data were obtained with an electron density model that consists of a uniform density component superimposed upon molecular-scale density oscillations (layers). The two types of layer profiles were observed to vary with temperature from 187-286 K . The amount of material in the molecular layers increases as that in the uniform density layer decreases, with the onset of liquid-to-layered transition occurring at a total film thickness of approximately 40 A (about twice the bulk correlation length of TEHOS).

Combined atomic force microscopy and fluorescence correlation spectroscopy measurements to study the dynamical structure of interfacial fluids

We have studied the dynamic structure of thin (approximately a few nanometers) liquid films of a nearly spherical, nonpolar molecule tetrakis(2-ethylhexoxy)silane (TEHOS) by using a combination of atomic force microscopy (AFM) and fluorescence correlation spectroscopy (FCS). Ultra-sensitive interferometer-based AFM was used to determine the stiffness (force gradient) and the damping coefficient of the liquid film. The experiments show oscillations in the damping coefficient with a period of approximately 1 nm, which is consistent with the molecular dimension of TEHOS as well as previous X-ray reflectivity measurements. Additionally, we performed FCS experiments for direct determination of the molecular dynamics within the liquid film. From the fluctuation autocorrelation curve, we measured the translational diffusion of the probe molecule embedded within the fluid film formed on a solid substrate. The autocorrelation function was best fitted with two components, which indicate that the dynamics are heterogeneous in nature. However, the heterogeneity is not as pronounced as had been previously observed for molecularly thin liquid films sandwiched between two solid substrates.

The intrinsic structure of the water surface

An operational procedure to obtain the intrinsic structure of liquid surfaces is applied here to a molecular dynamics simulation of water, with a model of point charges for the molecular interactions. The method, which had been recently proposed and used for simple fluids, is successfully extended to a molecular liquid with the complex bond structure of water. The elimination of the capillary wave fluctuations, in the intrinsic density and orientation profiles, gives a new overall view of the water surface, at the sharpest molecular level, and without the size-dependent broadening observed in the mean profiles. The molecules belonging to the outer liquid layer are clearly identified, and we find that only these molecules exhibit a clear preferential orientation to lie flat on the surface. Moreover, there is a strong correlation between the dipolar structure and the local curvatures of the intrinsic surface, so that at the extrusions of the intrinsic surface the molecular dipoles point preferentially toward the vapor side of the interface. Finally, we have found an intrinsic density layering structure, although the inner structure is strongly damped beyond the second layer.