Collective motion and its connection to the energy landscape in 2D soft crystals.

We examine the collective motion in computational models of a two-dimensional dusty plasma crystal and a charged colloidal suspension as they approach their respective melting transitions. To unambiguously identify rearrangement events in the crystal, we map the trajectory of configurations from an equilibrium molecular dynamics simulation to the corresponding sequence of configurations of local potential energy minima ("inherent structures"). This inherent structure (IS) trajectory eliminates the ambiguity that arises from localized vibrational motion.

Computationally efficient machine-learned model for GST phase change materials via direct and indirect learning.

Phase change materials such as Ge2Sb2Te5 (GST) are ideal candidates for next-generation, non-volatile, solid-state memory due to the ability to retain binary data in the amorphous and crystal phases and rapidly transition between these phases to write/erase information. Thus, there is wide interest in using molecular modeling to study GST. Recently, a Gaussian Approximation Potential (GAP) was trained for GST to reproduce Density Functional Theory (DFT) energies and forces at a fraction of the computational cost [Zhou et al.

Hierarchy of topological transitions in a network liquid.

The representation of complex systems as networks has become a critical tool across many fields of science. In the context of physical networks, such as biological neural networks, vascular networks, or network liquids where the nodes and edges occupy volume in three-dimensional space, the question of how they become densely packed is of special importance. Here, we investigate a model network liquid, which is known to densify via two successive liquid-liquid phase transitions (LLPTs).

Solution properties of spherical gold nanoparticles with grafted DNA chains from simulation and theory.

There has been a rapidly growing interest in the use of functionalized Au nanoparticles (NPs) as platforms in multiple applications in medicine and manufacturing. The sensing and targeting characteristics of these NPs, and the realization of precisely organized structures in manufacturing applications using such NPs, depend on the control of their surface functionalization. NP functionalization typically takes the form of polymer grafted layers, and a detailed knowledge of the chemical and structural properties of these layers is required to molecularly engineer the particle characteristics for specific applications.

Relating dynamic free volume to cooperative relaxation in a glass-forming polymer composite.

There are a variety of complementary descriptions of the temperature dependence of the structural relaxation time τ in glass-forming materials, which we interpret positively as suggesting an underlying unified description. We examine the inter-relation between the string model, an outgrowth of the Adam and Gibbs approach that emphasizes collective particle exchange motion, and the localization model, which emphasizes the volume explored by particles in their caged states, a kind of dynamic "free volume." Each model of liquid dynamics is described by a limited set of parameters that must be interrelated if both descriptions simultaneously describe the relaxation behavior.

The effect of nanoparticle softness on the interfacial dynamics of a model polymer nanocomposite.

The introduction of soft organic nanoparticles (NPs) into polymer melts has recently expanded the material design space for polymer nanocomposites, compared to traditional nanocomposites that utilize rigid NPs, such as silica, metallic NPs, and other inorganic NPs. Despite advances in the fabrication and characterization of this new class of materials, the effect of NP stiffness on the polymer structure and dynamics has not been systematically investigated. Here, we use molecular dynamics to investigate the segmental dynamics of the polymer interfacial region of isolated NPs of variable stiffness in a polymer matrix.

Reactive Molecular Dynamics Simulations of the Depolymerization of Polyethylene Using Graphene-Oxide-Supported Platinum Nanoparticles.

While plastic materials offer many benefits to society, the slow degradation and difficulty in recycling plastics raise important environmental and sustainability concerns. Traditional recycling efforts often lead to materials with inferior properties and correspondingly lower value, making them uneconomical to recycle. Recent efforts have shown promising chemical pathways for converting plastic materials into a wide range of value-added products, feedstocks or monomers.

Explaining the Sensitivity of Polymer Segmental Relaxation to Additive Size Based on the Localization Model.

We use molecular simulations to examine how the dynamics of a coarse-grained polymer melt are altered by additives of variable size and interaction strength with the polymer matrix. The effect of diluent size σ on polymer dynamics changes significantly when its size is comparable to the polymer segment size. For each σ, we show that the localization model (LM) quantitatively describes the dependence of the segmental relaxation time τ on temperature T in terms of dynamic free volume, quantified by the Debye-Waller factor ⟨u^{2}⟩.

Activation free energy gradient controls interfacial mobility gradient in thin polymer films.

We examine the mobility gradient in the interfacial region of substrate-supported polymer films using molecular dynamics simulations and interpret these gradients within the string model of glass-formation. No large gradients in the extent of collective motion exist in these simulated films, and an analysis of the mobility gradient on a layer-by-layer basis indicates that the string model provides a quantitative description of the relaxation time gradient. Consequently, the string model indicates that the interfacial mobility gradient derives mainly from a gradient in the high-temperature activation enthalpy ΔH and entropy ΔS as a function of depth z, an effect that exists even in the high-temperature Arrhenius relaxation regime far above the glass transition temperature.

Detecting bound polymer layers in attractive polymer-nanoparticle hybrids.

When polymer-nanoparticle (NP) attractions are sufficiently strong, a bound polymer layer with a distinct dynamic signature spontaneously forms at the NP interface. A similar phenomenon occurs near a fixed attractive substrate for thin polymer films. While our previous simulations fixed the NPs to examine the dilute limit, here, we allow the NP to move.

Predictive relation for the α-relaxation time of a coarse-grained polymer melt under steady shear.

We examine the influence of steady shear on structural relaxation in a simulated coarse-grained unentangled polymer melt over a wide range of temperature and shear rates. Shear is found to progressively suppress the α-relaxation process observed in the intermediate scattering function, leading ultimately to a purely inertially dominated β-relaxation at high shear rates, a trend similar to increasing temperature. On the basis of a scaling argument emphasizing dynamic heterogeneity in cooled liquids and its alteration under material deformation, we deduce and validate a parameter-free scaling relation for both the structural relaxation time τ from the intermediate scattering function and the "stretching exponent" β quantifying the extent of dynamic heterogeneity over the entire range of temperatures and shear rates that we can simulate.

Reconciling computational and experimental trends in the temperature dependence of the interfacial mobility of polymer films.

Many measurements have indicated that thin polymer films in their glass state exhibit a mobile interfacial layer that grows in thickness upon heating, while some measurements indicate the opposite trend. Moreover, simulations and limited measurements on glass-forming liquids at temperatures above the glass transition temperature T exhibit a growing interfacial mobility scale ξ upon cooling. To better understand these seemingly contradictory trends, we perform molecular dynamics simulations over a temperature regime for which our simulated polymer film enters a non-equilibrium glassy state and find that the relaxation time τ within the film interior, relative to the polymer-air interfacial layer, exhibits a maximum near the computational T.

Dynamic heterogeneity and collective motion in star polymer melts.

While glass formation of linear chain polymer melts has widely been explored, comparatively little is known about glass formation in star polymer melts. We study the segmental dynamics of star polymer melts via molecular dynamics simulations and examine the cooperative nature of segmental motion in star melts. In particular, we quantify how the molecular architecture of star polymers, i.

What does the instantaneous normal mode spectrum tell us about dynamical heterogeneity in glass-forming fluids?

We examine the instantaneous normal mode spectrum of model metallic and polymeric glass-forming liquids. We focus on the localized modes in the unstable part of the spectrum [unstable localized (UL) modes] and find that the particles making the dominant contribution to the participation ratio form clusters that grow upon cooling in a fashion similar to the dynamical heterogeneity in glass-forming fluids, i.e.

The interfacial zone in thin polymer films and around nanoparticles in polymer nanocomposites.

We perform coarse-grained simulations of model unentangled polymer materials to quantify the range over which interfaces alter the structure and dynamics in the vicinity of the interface. We study the interfacial zone around nanoparticles (NPs) in model polymer-NP composites with variable NP diameter, as well as the interfacial zone at the solid substrate and free surface of thin supported polymer films. These interfaces alter both the segmental packing and mobility in an interfacial zone.

State variables for glasses: The case of amorphous ice.

Glasses are out-of-equilibrium systems whose state cannot be uniquely defined by the usual set of equilibrium state variables. Here, we seek to identify an expanded set of variables that uniquely define the state of a glass. The potential energy landscape (PEL) formalism is a useful approach within statistical mechanics to describe supercooled liquids and glasses.

Collective Motion in the Interfacial and Interior Regions of Supported Polymer Films and Its Relation to Relaxation.

To understand the role of collective motion in the often large changes in interfacial molecular mobility observed in polymer films, we investigate the extent of collective motion in the interfacial regions of a thin supported polymer film and within the film interior by molecular dynamics simulation. Contrary to commonly stated expectations, we find that the extent of collective motion, as quantified by string-like molecular exchange motion, is similar in magnitude in the polymer-air interfacial layer as the film interior and distinct from the bulk material. This finding is consistent with Adam-Gibbs description of the segmental dynamics within mesoscopic film regions, where the extent of collective motion is related to the configurational entropy of the film as a rather than a locally defined extent of collective motion or configurational entropy.

The Stability of a Nanoparticle Diamond Lattice Linked by DNA.

The functionalization of nanoparticles (NPs) with DNA has proven to be an effective strategy for self-assembly of NPs into superlattices with a broad range of lattice symmetries. By combining this strategy with the DNA origami approach, the possible lattice structures have been expanded to include the cubic diamond lattice. This symmetry is of particular interest, both due to the inherent synthesis challenges, as well as the potential valuable optical properties, including a complete band-gap.

Diminishing Interfacial Effects with Decreasing Nanoparticle Size in Polymer-Nanoparticle Composites.

Using molecular simulations on model polymer nanocomposites at fixed filler loading, we show that interfacial polymer dynamics are affected less with decreasing nanoparticle (NP) size. However, the glass transition temperature T_{g} changes substantially more for an extremely small NP. The reason for this apparent contradiction is that the mean NP spacing decreases with decreasing particle size.

Why we need to look beyond the glass transition temperature to characterize the dynamics of thin supported polymer films.

There is significant variation in the reported magnitude and even the sign of [Formula: see text] shifts in thin polymer films with nominally the same chemistry, film thickness, and supporting substrate. The implicit assumption is that methods used to estimate [Formula: see text] in bulk materials are relevant for inferring dynamic changes in thin films. To test the validity of this assumption, we perform molecular simulations of a coarse-grained polymer melt supported on an attractive substrate.

String-like collective motion in the α- and β-relaxation of a coarse-grained polymer melt.

Relaxation in glass-forming liquids occurs as a multi-stage hierarchical process involving cooperative molecular motion. First, there is a "fast" relaxation process dominated by the inertial motion of the molecules whose amplitude grows upon heating, followed by a longer time α-relaxation process involving both large-scale diffusive molecular motion and momentum diffusion. Our molecular dynamics simulations of a coarse-grained glass-forming polymer melt indicate that the fast, collective motion becomes progressively suppressed upon cooling, necessitating large-scale collective motion by molecular diffusion for the material to relax approaching the glass-transition.

Valence, loop formation and universality in self-assembling patchy particles.

Patchy particles have emerged as an attractive model to mimic phase separation and self-assembly of globular proteins solutions, colloidal patchy particles, and molecular fluids where directional interactions are operative. In our previous work, we extensively explored the coupling of directional and isotropic interactions on both the phase separation and self-assembly in a system of patchy particles with five spots. Here, we extend this work to consider different patch valences and isotropic interaction strengths with an emphasis on self-assembly.

Molecular rigidity and enthalpy-entropy compensation in DNA melting.

Enthalpy-entropy compensation (EEC) is observed in diverse molecular binding processes of importance to living systems and manufacturing applications, but this widely occurring phenomenon is not sufficiently understood from a molecular physics standpoint. To gain insight into this fundamental problem, we focus on the melting of double-stranded DNA (dsDNA) since measurements exhibiting EEC are extensive for nucleic acid complexes and existing coarse-grained models of DNA allow us to explore the influence of changes in molecular parameters on the energetic parameters by using molecular dynamics simulations. Previous experimental and computational studies have indicated a correlation between EEC and changes in molecular rigidity in certain binding-unbinding processes, and, correspondingly, we estimate measures of DNA molecular rigidity under a wide range of conditions, along with resultant changes in the enthalpy and entropy of binding.

Influence of sample preparation on the transformation of low-density to high-density amorphous ice: An explanation based on the potential energy landscape.

Experiments and computer simulations of the transformations of amorphous ices display different behaviors depending on sample preparation methods and on the rates of change of temperature and pressure to which samples are subjected. In addition to these factors, simulation results also depend strongly on the chosen water model. Using computer simulations of the ST2 water model, we study how the sharpness of the compression-induced transition from low-density amorphous ice (LDA) to high-density amorphous ice (HDA) is influenced by the preparation of LDA.

Effects of a "bound" substrate layer on the dynamics of supported polymer films.

It is widely appreciated that an attractive polymer-substrate interaction can slow relaxation in thin supported polymer films and polymer nanocomposites. Recent measurements and simulations on nancomposites have indicated that this slowing of polymer dynamics occurs more strongly near a highly attractive particle surface where a "bound" layer having a much lower mobility can form, strongly influencing the thermodynamics and dynamics of the film. Here we use molecular simulations to show that a bound interfacial layer having a very similar nature arises in thin supported polymer films when the polymer-polymer attraction is stronger than the polymer-polymer interaction strength.