The FADE mass-stat: A technique for inserting or deleting particles in molecular dynamics simulations

Multiscale simulation of fluids: coupling molecular and continuum

Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the same simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between the two main types of coupling, state and flux, and discuss the varying assumptions in their use. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines.

Shape optimization of a meniscus-adherent nanotip

A soluble tip can dissolve into a tip with curvature when partially immersed in a liquid. This process has been used in the manufacture of sophisticated tips. However, it is difficult to observe the dissolution process in the laboratory, and the dissolution mechanisms at the nanoscale still need to be better understood. Here we utilize molecular dynamics simulations to study the dissolution process of a meniscus-adherent nanotip. The tip apex curvature radius reaches its minimum in the intermediate state. The shape of this state is defined as the optimized shape, which can be used as the termination criterion in applications. In addition, the shape of one optimized tip can be well-fitted to a double-Boltzmann function. The upper Boltzmann curve of this function forms via the competition between the chemical potential influence and the intermolecular forces, while the formation of the lower Boltzmann curve is controlled by the chemical potential influence. The parameters of the double-Boltzmann function are strongly correlated with the nanotip's initial configuration and dissolubility. A shape factor ξ is proposed to characterize the sharpness of optimized tips. Theory and simulations show that optimized tips possess a greater ability to shield the capillary effect than common tips. Our findings elucidate the meniscus-adherent nanotip's dissolution process and provide theoretical support for nano-instrument manufacture.

Moving Contact Lines: Linking Molecular Dynamics and Continuum-Scale Modeling

Despite decades of research, the modeling of moving contact lines has remained a formidable challenge in fluid dynamics whose resolution will impact numerous industrial, biological, and daily life applications. On the one hand, molecular dynamics (MD) simulation has the ability to provide unique insight into the microscopic details that determine the dynamic behavior of the contact line, which is not possible with either continuum-scale simulations or experiments. On the other hand, continuum-based models provide a link to the macroscopic description of the system. In this Feature Article, we explore the complex range of physical factors, including the presence of surfactants, which governs the contact line motion through MD simulations. We also discuss links between continuum- and molecular-scale modeling and highlight the opportunities for future developments in this area.

Electric fields can control the transport of water in carbon nanotubes

The properties of water confined inside nanotubes are of considerable scientific and technological interest. We use molecular dynamics to investigate the structure and average orientation of water flowing within a carbon nanotube. We find that water exhibits biaxial paranematic liquid crystal ordering both within the nanotube and close to its ends. This preferred molecular ordering is enhanced when an axial electric field is applied, affecting the water flow rate through the nanotube. A spatially patterned electric field can minimize nanotube entrance effects and significantly increase the flow rate.

Wetting and evaporation of salt-water nanodroplets: A molecular dynamics investigation

We employ molecular dynamics simulations to study the wetting and evaporation of salt-water nanodroplets on platinum surfaces. Our results show that the contact angle of the droplets increases with the salt concentration. To verify this, a second simulation system of a thin salt-water film on a platinum surface is used to calculate the various surface tensions. We find that both the solid-liquid and liquid-vapor surface tensions increase with salt concentration and as a result these cause an increase in the contact angle. However, the evaporation rate of salt-water droplets decreases as the salt concentration increases, due to the hydration of salt ions. When the water molecules have all evaporated from the droplet, two forms of salt crystals are deposited, clump and ringlike, depending on the solid-liquid interaction strength and the evaporation rate. To form salt crystals in a ring, it is crucial that there is a pinned stage in the evaporation process, during which salt ions can move from the center to the rim of the droplets. With a stronger solid-liquid interaction strength, a slower evaporation rate, and a higher salt concentration, a complete salt crystal ring can be deposited on the surface.