Jianguo Zhang

An Atomistic Simulation for 4-Cyano-4'-pentylbiphenyl and Its Homologue with a Reoptimized Force Field

Liquid crystals playing a crucial role in material sciences show increasing potential applications in nanotechnology and industry. Generally, thermodynamic and dynamic properties of liquid crystals strongly depend on the corresponding force fields (FF); thus, it is necessary and urgent for us to establish a reliable force field for a given liquid crystal system. In this paper, we develop a new set of FF parameters for the 5CB (4-cyano-4-pentylbiphenyl) molecule by reoptimizing some parameters of TraPPE-UA in order to reproduce the bulk density. This strategy for the construction of 5CB FF is rather advisable as it not only provides reliable values for the Lennard-Jones parameters but also reduces the computational cost and maintains FF transferability. Indeed, our simulation results show that the phase behavior, the order parameter, conformational features, neighboring molecular pair arrangements, and diffusion properties of 5CB can be reproduced very well. We further validate the transferability of this 5CB FF by extending it to the 8CB (4-cyano-4-octylbiphenyl) system. As a result, both the nematic and the partial bilayer smectic phases (Sm-A(d)) and the nematic-isotropic and the smectic-nematic transition temperatures as well as the diffusion properties of 8CB are successfully reproduced. Therefore, this set of FF parameters originally designed for the 5CB molecule is reliable and transferable. Its effectiveness to model nCB series and molecules with similar chemical structures is expected.

Evaporation of nanodroplets on heated substrates: a molecular dynamics simulation study.

Molecular dynamics simulations of Lennard-Jones particles have been performed to study the evaporation behavior of nanodroplets on heated substrates. The influence of the liquid-substrate interaction strength on the evaporation properties was addressed. Our results show that, during the temperature-raising evaporation, the gas is always hotter than the droplet. In contrast to the usual experimental conditions, the droplet sizes in our simulations are in the nanometer scale range and the substrates are ideally smooth and chemically homogeneous. As a result, no pinning was observed in our simulations for substrates denoted either hydrophilic (contact angle θ < 90°) or hydrophobic (contact angle θ > 90°). The evaporative mass flux is stronger with increasing hydrophilicity of the substrate because the heat transfer from the substrate to the droplet is more efficient for stronger attraction between the solid and the droplet. Evaporation and heat transfer to the gas phase occur preferentially in the vicinity of the three-phase contact line in the hydrophilic system. However, in the case of a hydrophobic substrate, there is no preferential location for mass and heat fluxes. During the whole evaporation process, no pure behavior according to either the constant-angle or the constant-radius model was found; both the contact angle and contact radius decrease for the droplets on hydrophilic and hydrophobic substrates alike.

Coarse-grained molecular dynamics simulations of the phase behavior of the 4-cyano-4'-pentylbiphenyl liquid crystal system.

In this paper, with the aim to establish a rational coarse-grained (CG) model for the 4-cyano-4-pentylbiphenyl (5CB) molecule, we construct three possible CG models (5P, 6P, and 7P) and then determine the bonded and nonbonded interaction parameters separately. For the intramolecular bonded interactions, the bond and angle distributions of the 5CB bulk phase are used as the target properties. For the nonbonded interactions between CG particles, we combine the structure-based and thermodynamic quantities-based methods for the parametrization of CG interaction potentials and attempt to use several fragment molecular systems to derive the CG nonbonded interaction parameters in order to maintain the transferability of our CG models to some extent. Finally, we fix the optimal nonbonded LJ parameters between CG bead pairs such that the results from CG simulations not only correctly reproduce the experimental density and the nematic LC state at 300 K and 1 atm but also reasonably approximate the local structural properties calculated from the underlying atomistic model. Through comparison of the resulting CG data with target properties, the 6P model is found to be the best one among the three, and then we use this model to investigate the phase behavior and dynamic properties. Our results show that the phase transition temperature from nematic to isotropic phase and the diffusion coefficients are reproduced very well, demonstrating the rationality of the 6P model. Our coarse-grained process should have implications for constructing CG models for nCB series or molecules with similar architectures.

Pinning of the Contact Line during Evaporation on Heterogeneous Surfaces: Slowdown or Temporary Immobilization? Insights from a Nanoscale Study

The question of the effect of surface heterogeneities on the evaporation of liquid droplets from solid surfaces is addressed through nonequilibrium molecular dynamics simulations. The mechanism behind contact line pinning which is still unclear is discussed in detail on the nanoscale. Model systems with the Lennard-Jones interaction potential were employed to study the evaporation of nanometer-sized cylindrical droplets from a flat surface. The heterogeneity of the surface was modeled through alternating stripes of equal width but two chemical types. The first type leads to a contact angle of 67°, and the other leads to a contact angle of 115°. The stripe width was varied between 2 and 20 liquid-particle diameters. On the surface with the narrowest stripes, evaporation occurred at constant contact angle as if the surface was homogeneous, with a value of the contact angle as predicted by the regular Cassie-Baxter equation. When the width was increased, the contact angle oscillated during evaporation between two boundaries whose values depend on the stripe width. The evaporation behavior was thus found to be a direct signature of the typical size of the surface heterogeneity domains. The contact angle both at equilibrium and during evaporation could be predicted from a local Cassie-Baxter equation in which the surface composition within a distance of seven fluid-particle diameters around the contact line was considered, confirming the local nature of the interactions that drive the wetting behavior of droplets. More importantly, we propose a nanoscale explanation of pinning during evaporation. Pinning should be interpreted as a drastic slowdown of the contact line dynamics rather than a complete immobilization of it during a transition between two contact angle boundaries.

Transferability of coarse-grained force field for nCB liquid crystal systems.

In this paper, the transferability of the coarse-grained (CG) force field originally developed for the liquid crystal (LC) molecule 5CB ( Zhang et al. J. Phys. Chem. B 2012 , 116 , 2075 - 2089 ) was investigated by its homologues 6CB and 8CB molecules. Note that, to construct the 5CB CG force field, we combined the structure-based and thermodynamic quantities-based methods and at the same time attempted to use several fragment molecular systems to derive the CG nonbonded interaction parameters. The resultant 5CB CG force field exhibits a good transferability to some extent. For example, not only the experimental densities, the local packing of atom groups, and the antiparallel arrangements of nearest neighboring molecules, but also the unique LC mesophases as well as the nematic-isotropic phase transition temperatures of 6CB and 8CB were reproduced. Meanwhile, the limitations of this 5CB CG force field were also observed, such as the phase transition from nematic to smectic was postponed to the lower temperature and the resulting smectic phase structure is single-layer-like instead of partially interdigitated bilayer-like as observed in underlying atomistic model. Apparently, more attention should be paid when applying a CG force field to the state point which is quite different from which the force field is explicitly parametrized for. The origin of the above limitations can be potentially traced back to the inherent simplifications and some approximations often adopted in the creation process of CG force field, for example, choosing symmetric CG potentials which do not explicitly include electrostatic interactions and are parametrized by reproducing the target properties of the specific nematic 5CB phase at 300 K and 1 atm, as well as using soft nonbonded potential and excluding torsion barriers. Moreover, although by construction this CG force field could inevitably incorporate both thermodynamic and local structural information on the nematic 5CB phase, the anisotropic diffusion coefficient ratios for different LC phases in both 6CB and 8CB systems are reproduced well. All these findings suggest that the multiproperty parametrization route together with fragment-based method provides a new approach to maximize the possibility to simultaneously reproduce multiple physical properties of a given molecule or related molecules with similar chemical structures at other state points.

A steady-state non-equilibrium molecular dynamics approach for the study of evaporation processes

Two non-equilibrium methods (called bubble method and splitting method, respectively) have been developed and tested to study the steady state evaporation of a droplet surrounded by its vapor, where the evaporation continuously occurs at the vapor-liquid interface while the droplet size remains constant. In the bubble method, gas molecules are continuously reinserted into a free volume (represented by a bubble) located at the centre of mass of the droplet to keep the droplet size constant. In the splitting method, a molecule close to the centre of mass of the droplet is split into two: In this way, the droplet size is also maintained during the evaporation. By additional local thermostats confined to the area of insertion, the effect of frequent insertions on properties such as density and temperature can be limited to the immediate insertion area. Perturbations are not observed in other parts of the droplet. In the end, both the bubble method and the splitting method achieve steady-state droplet evaporation. Although these methods have been developed using an isolated droplet, we anticipate that they will find a wide range of applications in the study of the evaporation of isolated films and droplets or thin films on heated substrates or under confinement. They can in principle also be used to study the steady-state of other physical processes, such as the diffusion or permeation of gas molecules or ions in a pressure gradient or a concentration gradient.

Formation of coffee-stain patterns at the nanoscale: The role of nanoparticle solubility and solvent evaporation rate

When droplets of nanoparticle suspension evaporate from surfaces, they leave behind a deposit of nanoparticles. The mechanism of evaporation-induced pattern formation in the deposit is studied by molecular dynamics simulations for sessile nanodroplets. The influence of the interaction between nanoparticles and liquid molecules and the influence of the evaporation rate on the final deposition pattern are addressed. When the nanoparticle-liquid interaction is weaker than the liquid-liquid interaction, an interaction-driven or evaporation-induced layer of nanoparticles appears at the liquid-vapor interface and eventually collapses onto the solid surface to form a uniform deposit independently of the evaporation rate. When the nanoparticle-liquid and liquid-liquid interactions are comparable, the nanoparticles are dispersed inside the droplet and evaporation takes place with the contact line pinned at a surface defect. In such a case, a pattern with an approximate ring-like shape is found with fast evaporation, while a more uniform distribution is observed with slower evaporation. When the liquid-nanoparticle interaction is stronger than the liquid-liquid interaction, evaporation always occurs with receding contact line. The final deposition pattern changes from volcano-like to pancake-like with decreasing evaporation rate. These findings might help to design nanoscale structures like nanopatterns or nanowires on surface through controlled solvent evaporation.

Extending reverse nonequilibrium molecular dynamics to the calculation of mutual diffusion coefficients in molecular fluid mixtures

Abstract In a recent publication, a reverse nonequilibrium molecular dynamics method was introduced to calculate the mutual diffusion (binary diffusion, interdiffusion) coefficient of binary Lennard-Jones mixture directly. In this work, modification is made to this method for calculation the mutual diffusion coefficient of complicated binary mixture of molecules containing the same number of particle in calculation. As test, this method is applied to calculate the mutual diffusion coefficient of water/deuterium oxide (H2O/D2O) and benzene/cyclohexane (C6H6/C6H12) systems. Steady mass flows are induced by suitably exchanging molecule positions and velocities in different regions. The analysis of the resulting steady state concentration profiles allows the calculation of mutual diffusion coefficient. The obtained results are qualitative agreement with the references. The excess energy introduced by swapping molecules becomes strong when the studied mixture contains large and/or complicated molecules. For the systems studied, however, a standard thermostat is still able to remove the excess energy efficiently.