Inducing stratification of colloidal mixtures with a mixed binary solvent.

Molecular dynamics simulations are used to demonstrate that a binary solvent can be used to stratify colloidal mixtures when the suspension is rapidly dried. The solvent consists of two components, one more volatile than the other. When evaporated at high rates, the more volatile component becomes depleted near the evaporation front and develops a negative concentration gradient from the bulk of the mixture to the liquid-vapor interface while the less volatile solvent is enriched in the same region and exhibit a positive concentration gradient. Such gradients can be used to drive a binary mixture of colloidal particles to stratify if one is preferentially attracted to the more volatile solvent and the other to the less volatile solvent. During solvent evaporation, the fraction of colloidal particles preferentially attracted to the less volatile solvent is enhanced at the evaporation front, whereas the colloidal particles having stronger attractions with the more volatile solvent are driven away from the interfacial region. As a result, the colloidal particles show a stratified distribution after drying, even if the two colloids have the same size.

Molecular Weight Distribution of Branched Polymers: Comparison between Monte Carlo Simulation and Flory-Stockmayer Theory.

It is challenging to predict the molecular weight distribution (MWD) for a polymer with a branched architecture, though such information will significantly benefit the design and development of branched polymers with desired properties and functions. A Monte Carlo (MC) simulation method based on the Gillespie algorithm is developed to quickly compute the MWD of branched polymers formed through step-growth polymerization, with a branched polyetherimide from two backbone monomers (4,4'-bisphenol A dianhydride and m-phenylenediamine), a chain terminator (phthalic anhydride), and a branching agent (tris[4-(4-aminophenoxy)phenyl] ethane) as an example. This polymerization involves four reactions that can be all reduced to a condensation reaction between an amine group and a carboxylic anhydride group. A comparison between the MC simulation results and the predictions of the Flory-Stockmayer theory on MWD shows that the rates of the reactions are determined by the concentrations of the functional groups on the monomers involved in each reaction. It further shows that the Flory-Stockmayer theory predicts MWD well for systems below the gel point but starts to fail for systems around or above the gel point. However, for all the systems, the MC method can be used to reliably predict MWD no matter if they are below or above the gel point. Even for a macroscopic system, a converging distribution can be quickly obtained through MC simulations on a system of only a few hundred to a few thousand monomers that have the same molar ratios as in the macroscopic system.

Modeling Solution Drying by Moving a Liquid-Vapor Interface: Method and Applications.

A method of simulating the drying process of a soft matter solution with an implicit solvent model by moving the liquid-vapor interface is applied to various solution films and droplets. For a solution of a polymer and nanoparticles, we observe "polymer-on-top" stratification, similar to that found previously with an explicit solvent model. Furthermore, "polymer-on-top" is found even when the nanoparticle size is smaller than the radius of gyration of the polymer chains. For a suspension droplet of a bidisperse mixture of nanoparticles, we show that core-shell clusters of nanoparticles can be obtained via the "small-on-outside" stratification mechanism at fast evaporation rates. "Large-on-outside" stratification and uniform particle distribution are also observed when the evaporation rate is reduced. Polymeric particles with various morphologies, including Janus spheres, core-shell particles, and patchy particles, are produced from drying droplets of polymer solutions by combining fast evaporation with a controlled interaction between the polymers and the liquid-vapor interface. Our results validate the applicability of the moving interface method to a wide range of drying systems. The limitations of the method are pointed out and cautions are provided to potential practitioners on cases where the method might fail.

Long-range DNA-water interactions.

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.

Physicochemical properties of respiratory droplets and their role in COVID-19 pandemics: a critical review.

The ongoing coronavirus disease 2019 (COVID-19) pandemic is a serious challenge faced by the global community. Physical scientists can help medical workers and biomedical scientists, engineers, and practitioners, who are working on the front line, to slow down and eventually contain the spread of the COVID-19 virus. This review is focused on the physicochemical characteristics, including composition, aerodynamics, and drying behavior of respiratory droplets as a complex and multicomponent soft matter system, which are the main carrier of the virus for interpersonal transmission. The distribution and dynamics of virus particles within a droplet are also discussed. Understanding the characteristics of virus-laden respiratory droplets can lead to better design of personal protective equipment, frequently touched surfaces such as door knobs and touchscreens, and filtering equipment for indoor air circulation. Such an understanding also provides the scientific basis of public policy, including social distancing rules and public hygiene guidelines, implemented by governments around the world.

Potential interference with microtubule assembly by graphene: a tug-of-war.

With the ever-increasing demand for graphene-based materials and their promising applications in numerous nanotechnologies, the biological effects of graphene on living systems have become crucial and ought to be well understood. Previously, both the cytotoxicity of graphene towards biological cells and its potential application as a nanomedicine have been revealed experimentally and theoretically. Besides many existing anticancer drugs that target microtubules, here we investigate the possibility of using graphene as a nanomedicine, which could alter the dynamic assembly and disassembly of a microtubule. We found that when a graphene nanosheet is at the hydrophilic interface of two neighboring heterodimers (containing α and ß tubulins), it can pull one dimer away from the other through a "tug-of-war" mechanism, driven by the strong dispersive interaction exerted by the surface of the graphene nanosheet. This work demonstrates that based on the existing methods for mitigating graphene's cytotoxicity (already developed in this field), a graphene-based nanomedicine could be designed to target microtubules of cancer cells and induce cell apoptosis.

Insights into Hydration Dynamics and Cooperative Interactions in Glycerol-Water Mixtures by Terahertz Dielectric Spectroscopy.

We report relaxation dynamics of glycerol-water mixtures as probed by megahertz-to-terahertz dielectric spectroscopy in a frequency range from 50 MHz to 0.5 THz at room temperature. The dielectric relaxation spectra reveal several polarization processes at the molecular level with different time constants and dielectric strengths, providing an understanding of the hydrogen-bonding network in glycerol-water mixtures. We have determined the structure of hydration shells around glycerol molecules and the dynamics of bound water as a function of glycerol concentration in solutions using the Debye relaxation model. The experimental results show the existence of a critical glycerol concentration of ∼7.5 mol %, which is related to the number of water molecules in the hydration layer around a glycerol molecule. At higher glycerol concentrations, water molecules dispersed in a glycerol network become abundant and eventually dominate, and four distinct relaxation processes emerge in the mixtures. The relaxation dynamics and hydration structure in glycerol-water mixtures are further probed with molecular dynamics simulations, which confirm the physical picture revealed by the dielectric spectroscopy.

Stratification of drying particle suspensions: Comparison of implicit and explicit solvent simulations.

Large scale molecular dynamics simulations are used to study drying suspensions of a binary mixture of large and small particles in explicit and implicit solvents. The solvent is first modeled explicitly and then mapped to a uniform viscous medium by matching the diffusion coefficients and the pair correlation functions of the particles. "Small-on-top" stratification of the particles, with an enrichment of the smaller ones at the receding liquid-vapor interface during drying, is observed in both models under the same drying conditions. With the implicit solvent model, we are able to model much thicker films and study the effect of the initial film thickness on the final distribution of particles in the dry film. Our results show that the degree of stratification is controlled by the Péclet number defined using the initial film thickness as the characteristic length scale. When the Péclet numbers of large and small particles are much larger than 1, the degree of "small-on-top" stratification is first enhanced and then weakens as the Péclet numbers are increased.

Control of Stratification in Drying Particle Suspensions via Temperature Gradients.

A potential strategy for controlling stratification in a drying suspension of bidisperse particles is studied using molecular dynamics simulations. When the suspension is maintained at a constant temperature during fast drying, it can exhibit "small-on-top" stratification with an accumulation (depletion) of smaller (larger) particles in the top region of the drying film, consistent with the prediction of current theories based on diffusiophoresis. However, when only the region near the substrate is thermalized at a constant temperature, a negative temperature gradient develops in the suspension because of evaporative cooling at the liquid-vapor interface. Since the associated thermophoresis is stronger for larger nanoparticles, a higher fraction of larger nanoparticles migrate to the top of the drying film at fast evaporation rates. As a result, stratification is converted to "large-on-top". Very strong small-on-top stratification can be produced with a positive thermal gradient in the drying suspension. Here, we explore a way to produce a positive thermal gradient by thermalizing the vapor at a temperature higher than that of the solvent. Possible experimental approaches to realize various thermal gradients in a suspension undergoing solvent evaporation and thus to produce different stratification states in the drying film are suggested.

The meniscus on the outside of a circular cylinder: From microscopic to macroscopic scales.

We systematically study the meniscus on the outside of a small circular cylinder vertically immersed in a liquid bath in a cylindrical container that is coaxial with the cylinder. The cylinder has a radius R much smaller than the capillary length, κ-1, and the container radius, L, is varied from a small value comparable to R to ∞. In the limit of L≪κ-1, we analytically solve the general Young-Laplace equation governing the meniscus profile and show that the meniscus height, Δh, scales approximately with Rln(L/R). In the opposite limit where L≫κ-1,Δh becomes independent of L and scales with Rln(κ-1/R). We implement a numerical scheme to solve the general Young-Laplace equation for an arbitrary L and demonstrate the crossover of the meniscus profile between these two limits. The crossover region has been determined to be roughly 0.4κ-1≲L≲4κ-1. An approximate analytical expression has been found for Δh, enabling its accurate prediction at any values of L that ranges from microscopic to macroscopic scales.

Stratification in Drying Films Containing Bidisperse Mixtures of Nanoparticles.

Large scale molecular dynamics simulations for bidisperse nanoparticle suspensions with an explicit solvent are used to investigate the effects of evaporation rates and volume fractions on the nanoparticle distribution during drying. Our results show that "small-on-top" stratification can occur when Pe sϕ s ≳ c with c ∼ 1, where Pe s is the Péclet number and ϕ s is the volume fraction of the smaller particles. This threshold of Pe sϕ s for "small-on-top" is larger by a factor of ∼α2 than the prediction of the model treating solvent as an implicit viscous background, where α is the size ratio between the large and small particles. Our simulations further show that when the evaporation rate of the solvent is reduced, the "small-on-top" stratification can be enhanced, which is not predicted by existing theories. This unexpected behavior is explained with thermophoresis associated with a positive gradient of solvent density caused by evaporative cooling at the liquid/vapor interface. For ultrafast evaporation the gradient is large and drives the nanoparticles toward the liquid/vapor interface. This phoretic effect is stronger for larger nanoparticles, and consequently the "small-on-top" stratification becomes more distinct when the evaporation rate is slower (but not too slow such that a uniform distribution of nanoparticles in the drying film is produced), as thermophoresis that favors larger particles on the top is mitigated. A similar effect can lead to "large-on-top" stratification for Pe sϕ s above the threshold when Pe s is large but ϕ s is small. Our results reveal the importance of including the solvent explicitly when modeling evaporation-induced particle separation and organization and point to the important role of density gradients brought about by ultrafast evaporation.

High-Precision Megahertz-to-Terahertz Dielectric Spectroscopy of Protein Collective Motions and Hydration Dynamics.

The low-frequency collective vibrational modes in proteins as well as the protein-water interface have been suggested as dominant factors controlling the efficiency of biochemical reactions and biological energy transport. It is thus crucial to uncover the mystery of the hydration structure and dynamics as well as their coupling to collective motions of proteins in aqueous solutions. Here, we report dielectric properties of aqueous bovine serum albumin protein solutions as a model system using an extremely sensitive dielectric spectrometer with frequencies spanning from megahertz to terahertz. The dielectric relaxation spectra reveal several polarization mechanisms at the molecular level with different time constants and dielectric strengths, reflecting the complexity of protein-water interactions. Combining the effective-medium approximation and molecular dynamics simulations, we have determined collective vibrational modes at terahertz frequencies and the number of water molecules in the tightly bound and loosely bound hydration layers. High-precision measurements of the number of hydration water molecules indicate that the dynamical influence of proteins extends beyond the first solvation layer, to around 7 Å distance from the protein surface, with the largest slowdown arising from water molecules directly hydrogen-bonded to the protein. Our results reveal critical information of protein dynamics and protein-water interfaces, which determine biochemical functions and reactivity of proteins.

Ordering nanoparticles with polymer brushes.

Ordering nanoparticles into a desired super-structure is often crucial for their technological applications. We use molecular dynamics simulations to study the assembly of nanoparticles in a polymer brush randomly grafted to a planar surface as the solvent evaporates. Initially, the nanoparticles are dispersed in a solvent that wets the polymer brush. After the solvent evaporates, the nanoparticles are either inside the brush or adsorbed at the surface of the brush, depending on the strength of the nanoparticle-polymer interaction. For strong nanoparticle-polymer interactions, a 2-dimensional ordered array is only formed when the brush density is finely tuned to accommodate a single layer of nanoparticles. When the brush density is higher or lower than this optimal value, the distribution of nanoparticles shows large fluctuations in space and the packing order diminishes. For weak nanoparticle-polymer interactions, the nanoparticles order into a hexagonal array on top of the polymer brush as long as the grafting density is high enough to yield a dense brush. An interesting healing effect is observed for a low-grafting-density polymer brush that can become more uniform in the presence of weakly adsorbed nanoparticles.

Nanocapillary Adhesion between Parallel Plates.

Molecular dynamics simulations are used to study capillary adhesion from a nanometer scale liquid bridge between two parallel flat solid surfaces. The capillary force, Fcap, and the meniscus shape of the bridge are computed as the separation between the solid surfaces, h, is varied. Macroscopic theory predicts the meniscus shape and the contribution of liquid/vapor interfacial tension to Fcap quite accurately for separations as small as two or three molecular diameters (1-2 nm). However, the total capillary force differs in sign and magnitude from macroscopic theory for h ≲ 5 nm (8-10 diameters) because of molecular layering that is not included in macroscopic theory. For these small separations, the pressure tensor in the fluid becomes anisotropic. The components in the plane of the surface vary smoothly and are consistent with theory based on the macroscopic surface tension. Capillary adhesion is affected by only the perpendicular component, which has strong oscillations as the molecular layering changes.

Dispersing Nanoparticles in a Polymer Film via Solvent Evaporation.

Large-scale molecular dynamics simulations are used to study the dispersion of nanoparticles (NPs) in a polymer film during solvent evaporation. As the solvent evaporates, a dense polymer-rich skin layer forms at the liquid/vapor interface, which is either NP rich or poor depending on the strength of the NP/polymer interaction. When the NPs are strongly wet by the polymer, the NPs accumulate at the interface and form layers. However, when the NPs are only partially wet by the polymer, most NPs are uniformly distributed in the bulk of the polymer film, with the dense skin layer serving as a barrier to prevent the NPs from moving to the interface. Our results point to a possible route to employ less favorable NP/polymer interactions and fast solvent evaporation to uniformly disperse NPs in a polymer film, contrary to the common belief that strong NP/polymer attractions are needed to make NPs well dispersed in polymer nanocomposites.

Capillary adhesion at the nanometer scale.

Molecular dynamics simulations are used to study the capillary adhesion from a nonvolatile liquid meniscus between a spherical tip and a flat substrate. The atomic structure of the tip, the tip radius, the contact angles of the liquid on the two surfaces, and the volume of the liquid bridge are varied. The capillary force between the tip and substrate is calculated as a function of their separation h. The force agrees with continuum predictions based on macroscopic theory for h down to ∼5 to 10 nm. At smaller h, the force tends to be less attractive than predicted and has strong oscillations. This oscillatory component of the capillary force is completely missed in the macroscopic theory, which only includes contributions from the surface tension around the circumference of the meniscus and the pressure difference over the cross section of the meniscus. The oscillation is found to be due to molecular layering of the liquid confined in the narrow gap between the tip and substrate. This effect is most pronounced for large tip radii and/or smooth surfaces. The other two components considered by the macroscopic theory are also identified. The surface tension term, as well as the meniscus shape, is accurately described by the macroscopic theory for h down to ∼1 nm, but the capillary pressure term is always more positive than the corresponding continuum result. This shift in the capillary pressure reduces the average adhesion by a factor as large as 2 from its continuum value and is found to be due to an anisotropy in the pressure tensor. The component in the plane of the substrate is consistent with the capillary pressure predicted by the macroscopic theory (i.e., the Young-Laplace equation), but the normal pressure that determines the capillary force is always more positive than the continuum counterpart.

Supramolecular assembly of asymmetric self-neutralizing amphiphilic peptide wedges.

Mimicking the remarkable dynamic and multifunctional utility of biological nanofibers, such as microtubules, is a challenging and technologically attractive objective in synthetic supramolecular chemistry. Understanding the complex molecular interactions that govern the assembly of synthetic materials, such as peptides, is key to meeting this challenge. Using molecular dynamics simulations to guide molecular design, we explore here the self-assembly of structurally and functionally asymmetric wedge-shaped peptides. Supramolecular assembly into nanofiber gels or multilayered lamellar structures was determined by cooperative influences of hydrogen bonding, amphiphilicity (hydrophilic asymmetry), and the distribution of electrostatic charges on the aqueous self-assembly of asymmetric peptides. Molecular amphiphilicity and ß-sheet forming capacity were both identified as necessary, but not independently sufficient, to form supramolecular nanofibers. Imbalances in positive and negative charges prevented nanofiber assembly, while the asymmetric distribution of balanced charges within a peptide is believed to affect peptide conformation and subsequent self-assembly into either nanofibers or lamellar structures. Insights into cooperative molecular interactions and the effects of molecular asymmetry on assembly may aid the development of next-generation supramolecular nanomaterial assemblies.

Self-assembly of chiral tubules.

The efficient and controlled assembly of complex structures from macromolecular building blocks is a critical open question in both biological systems and nanoscience. Using molecular dynamics simulations we study the self-assembly of tubular structures from model macromolecular monomers with multiple binding sites on their surfaces [Cheng et al., Soft Matter, 2012, 8, 5666-5678]. In this work we add chirality to the model monomer and a lock-and-key interaction. The self-assembly of free monomers into tubules yields a pitch value that often does not match the chirality of the monomer (including achiral monomers). We show that this mismatch occurs because of a twist deformation that brings the lateral interaction sites into alignment when the tubule pitch differs from the monomer chirality. The energy cost for this deformation is small as the energy distributions substantially overlap for small differences in the pitch and chirality. In order to control the tubule pitch by preventing the twist deformation, the interaction between the vertical surfaces must be increased without resulting in kinetically trapped structures. For this purpose, we employ lock-and-key interactions and obtain good control of the self-assembled tubule pitch. These results explain some fundamental features of microtubules. The vertical interaction strength is larger than the lateral in microtubules because this yields a controlled assembly of tubules with the proper pitch. We also generally find that the control of the assembly into tubules is difficult, which explains the wide range of pitch values and protofilament numbers observed in microtubule assembly.

Molecular dynamics simulations of evaporation-induced nanoparticle assembly.

While evaporating solvent is a widely used technique to assemble nano-sized objects into desired superstructures, there has been limited work on how the assembled structures are affected by the physical aspects of the process. We present large scale molecular dynamics simulations of the evaporation-induced assembly of nanoparticles suspended in a liquid that evaporates in a controlled fashion. The quality of the nanoparticle crystal formed just below the liquid/vapor interface is found to be better at relatively slower evaporation rates, as less defects and grain boundaries appear. This trend is understood as the result of the competition between the accumulation and diffusion times of nanoparticles at the liquid/vapor interface. When the former is smaller, nanoparticles are deposited so fast at the interface that they do not have sufficient time to arrange through diffusion, which leads to the prevalence of defects and grain boundaries. Our results have important implications in understanding assembly of nanoparticles and colloids in non-equilibrium liquid environments.

Structure and diffusion of nanoparticle monolayers floating at liquid/vapor interfaces: a molecular dynamics study.

Large-scale molecular dynamics simulations are used to simulate a layer of nanoparticles floating on the surface of a liquid. Both a low viscosity liquid, represented by Lennard-Jones monomers, and a high viscosity liquid, represented by linear homopolymers, are studied. The organization and diffusion of the nanoparticles are analyzed as the nanoparticle density and the contact angle between the nanoparticles and liquid are varied. When the interaction between the nanoparticles and liquid is reduced the contact angle increases and the nanoparticles ride higher on the liquid surface, which enables them to diffuse faster. In this case the short-range order is also reduced as seen in the pair correlation function. For the polymeric liquids, the out-of-layer fluctuation is suppressed and the short-range order is slightly enhanced. However, the diffusion becomes much slower and the mean square displacement even shows sub-linear time dependence at large times. The relation between diffusion coefficient and viscosity is found to deviate from that in bulk diffusion. Results are compared to simulations of the identical nanoparticles in 2-dimensions.