Soft matter with soft particles

Computational investigation of the effects of polymer grafting on the effective interaction between silica nanoparticles in water

Understanding and control of the effective interaction between nanoscale building blocks (colloids or nanoparticles) dispersed in a solvent is an important prerequisite for the development of bottom-up design strategies for soft functional materials. Here, we have employed all-atom molecular dynamics simulations to investigate the impact of polymer grafting on the solvent-mediated effective interaction between the silica nanoparticles (Si-NPs) in water, and in turn, on its bulk structural and thermodynamic properties. We found that the nature of the short grafting polymers [characterized by their interaction with water (hydrophobicity or hydrophilicity) and molecular weight] has a profound effect on the range and strength of the effective interaction between the Si-NPs. The hydrophobic polymer [such as polyethylene (PE)]-grafting of Si-NP gives rise to a more attractive interaction between the Si-NPs compared to the hydrophilic polymer [such as polyethylene glycol (PEG)] and non-grafted cases. This study further provides fundamental insights into the molecular origin of the observed behavior of the effective pair interactions between the grafted Si-NPs. For PE-grafted Si-NPs, the confined water (water inside the cavity formed by a pair of Si-NPs) undergoes a partial dewetting transition on approaching below a critical inter-particle separation leading to a stronger attractive interaction. Furthermore, we report that the effective attraction between the PE-grafted Si-NPs can be reliably controlled by changing the grafting PE density. We have also investigated the bulk structural and thermodynamic behavior of the coarse-grained Si-NP system where the particles interact via effective interaction in the absence of water. We believe that the insights gained from this work are important prerequisites for formulating rational bottom-up design strategies for functional materials where nano- (or, colloidal) particles are the building blocks.

Engineering Ultrasoft Interactions in Stiff All-DNA Dendrimers by Site-Specific Control of Scaffold Flexibility

A combined experimental and theoretical study of the structural correlations in moderately concentrated suspensions of all-DNA dendrimers of the second generation (G2) with controlled scaffold rigidity is reported here. Small-angle X-ray scattering experiments in concentrated aqueous saline solutions of stiff all-DNA G2 dendritic constructs reveal a novel anomalous liquid-like phase behavior which is reflected in the calculated structure factors as a two-step increase at low scattering wave vectors. By developing a new design strategy for adjusting the particle's internal flexibility based on site-selective incorporation of single-stranded DNA linkers into the dendritic scaffold, it is shown that this unconventional type of self-organization is strongly contingent on the dendrimer's stiffness. A comprehensive computer simulation study employing dendritic models with different levels of coarse-graining, and two theoretical approaches based on effective, pair-potential interactions, remarkably confirmed the origin of this unusual liquid-like behavior. The results demonstrate that the precise control of the internal structure of the dendritic scaffold conferred by the DNA can be potentially used to engineer a rich palette of novel ultrasoft interaction potentials that could offer a route for directed self-assembly of intriguing soft matter phases and experimental realizations of a host of unusual phenomena theoretically predicted for ultrasoft interacting systems.

M13-Phage-Based Star-Shaped Particles with Internal Flexibility

We report on the construction and the dynamics of monodisperse star-shaped particles, mimicking, at the mesoscale, star polymers. Such multiarm star-like particles result from the self-assembly of gold nanoparticles, forming the core, with tip-linked filamentous viruses (M13 bacteriophages) acting as spines in a sea urchin-like structure. By combining fluorescence and dark-field microscopy with dynamic light scattering, we investigate the diffusion of these hybrid spiny particles. We reveal the internal dynamics of the star particles by probing their central metallic core, which exhibits a hindered motion that can be described as a Brownian particle trapped in a harmonic potential. We therefore show that the filamentous viruses and specifically their tip proteins behave as entropic springs, extending the relevance of the study of such hybrid mesoscopic analogues of star polymers to phage biotechnology.

The entropy-controlled strategy in self-assembling systems

Self-assembly of various building blocks has been considered as a powerful approach to generate novel materials with tailorable structures and optimal properties. Understanding physicochemical interactions and mechanisms related to structural formation and transitions is of essential importance for this approach. Although it is well-known that diverse forces and energies can significantly contribute to the structures and properties of self-assembling systems, the potential entropic contribution remains less well understood. The past few years have witnessed rapid progress in addressing the entropic effects on the structures, responses, and functions in the self-assembling systems, and many breakthroughs have been achieved. This review provides a framework regarding the entropy-controlled strategy of self-assembly, through which the structures and properties can be tailored by effectively tuning the entropic contribution and its interplay with the enthalpic counterpart. First, we focus on the fundamentals of entropy in thermodynamics and the entropy types that can be explored for self-assembly. Second, we discuss the rules of entropy in regulating the structural organization in self-assembly and delineate the entropic force and superentropic effect. Third, we introduce the basic principles, significance and approaches of the entropy-controlled strategy in self-assembly. Finally, we present the applications where this strategy has been employed in fields like colloids, macromolecular systems and nonequilibrium assembly. This review concludes with a discussion on future directions and future research opportunities for developing and applying the entropy-controlled strategy in complex self-assembling systems.

Multi-particle collision dynamics for a coarse-grained model of soft colloids applied to ring polymers

The simulation of polymer solutions often requires the development of methods that accurately include hydrodynamic interactions. Resolution on the atomistic scale is too computationally expensive to cover mesoscopic time and length scales on which the interesting polymer phenomena are observed. Therefore, coarse-graining methods have to be applied. In this work, the solvent is simulated using the well-established multi-particle collision dynamics scheme, and for the polymer, different coarse-graining methods are employed and compared against the monomer resolved Kremer-Grest model by their resulting diffusion coefficients. This research builds on previous work [Ruiz-Franco et al., J. Chem. Phys. 151, 074902 (2019)], in which star polymers and linear chains in a solvent were simulated and two different coarse-graining methods were developed, in order to increase computational efficiency. The present work extends this approach to ring polymers and seeks to refine one of the authors' proposed model: the penetrable soft colloid model. It was found that both proposed models are not well suited to ring polymers; however, the introduction of a factor to the PSC model delivers satisfying results for the diffusion behavior by regulating the interaction intensity with the solvent.

Electron hopping heat transport in molecules

The realization of single-molecule thermal conductance measurements has driven the need for theoretical tools to describe conduction processes that occur over atomistic length scales. In macroscale systems, the principle that is typically used to understand thermal conductivity is Fourier's law. At molecular length scales, however, deviations from Fourier's law are common in part because microscale thermal transport properties typically depend on the complex interplay between multiple heat conduction mechanisms. Here, the thermal transport properties that arise from electron transfer across a thermal gradient in a molecular conduction junction are examined theoretically. We illustrate how transport in a model junction is affected by varying the electronic structure and length of the molecular bridge in the junction as well as the strength of the coupling between the bridge and its surrounding environment. Three findings are of note: First, the transport properties can vary significantly depending on the characteristics of the molecular bridge and its environment; second, the system's thermal conductance commonly deviates from Fourier's law; and third, in properly engineered systems, the magnitude of electron hopping thermal conductance is similar to what has been measured in single-molecule devices.

Association of two polyelectrolytes in salt solutions

The association of polyelectrolytes (PEs) in solution affects a wealth of structural and dynamic behaviors, and is also fundamentally important for an understanding of protein association and aggregation. Here, we theoretically study the association of two PE chains by addressing the stability and morphology of the non-spherical associates. Our theory predicts that an elongated pearl-necklace (PN) associate can be stable at high salt concentrations due to the screened electrostatic repulsion. This contradicts the implication of scaling theory. In addition, there is no one-to-one correspondence between the morphology of the associate and its constituting unimers, which is demonstrated by the existence of different association modes.

Impact of free energy of polymers on polymorphism of polymer-grafted nanoparticles

Colloidal crystals have gathered wide attention as a model material for optical applications because of their feasibility in controlling the propagation of light by their crystal structure and lattice spacing as well as the simplicity of their fabrication. However, due to the simple interaction between colloids, the colloidal crystal structures that can be formed are limited. It is also difficult to adjust the lattice spacing. Furthermore, colloidal crystals are fragile compared to other crystals. In this study, we focused on polymer-grafted nanoparticles (PGNP) as a possible solution to these unresolved issues. We expected that PGNPs, composed of two distinct layers (the hard core of a nanoparticle and the soft corona of grafted polymers on the surface), will demonstrate similar behaviors as star polymers and hard spheres. We also predicted that PGNPs may exhibit polymorphism because the interaction between PGNPs strongly depends upon their grafting density and the length of the grafted polymer chains. Moreover, we expected that crystals made from PGNPs will be structurally tough due to the entanglement of grafted polymers. From exploration of crystal polymorphs of PGNPs by molecular dynamics simulations, we found face-centered cubic (FCC)/hexagonal close-packed (HCP) and body-centered cubic (BCC) crystals, depending on the length of the grafted polymer chains. When the chains were short, PGNPs behaved like hard spheres and crystals were arranged in FCC/HCP structure, much like the phase transition observed in an Alder transition. When the chains were long enough, the increase in the free energy of grafted polymers was no longer negligible and crystals were arranged in BCC structure, which has a lower density than FCC/HCP. When the chains were not too short or long, FCC/HCP structures were first observed when the volume fraction of system was small, but a phase transition occurred when the system was further compressed and the crystals arranged themselves in a BCC structure. These results most likely have laid strong foundations for future simulations and experimental studies of PGNP crystals.

Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations

A generic model was developed for studies of the polymerization process of regular branched macromolecules. Monte Carlo simulations were performed employing the Dynamic Lattice Liquid algorithm to study this process. A core-first methodology was used in a living polymerization of stars with up to 32 arms, and dendrimers consisted of 4-functional segments. The kinetics of the synthesis process for stars with different numbers of branches and dendrimers was compared. The size and structure of star-branched polymers and dendrimers during the synthesis were studied. The influence of the functionality of well-defined cores on the structure and on the dispersity of the system was also examined. The differences in the kinetics in the formation of both architectures, as well as changes to their structures, were described and discussed.

Bond Orientational Order Parameters for Classifying Solid-like Clusters in a Lennard-Jones System near Liquid-Solid Transition and at Solid States

In this paper, we introduced an order parameter, named the local structure similarity (LSS), to measure the resemblance of a cluster structure in a liquid with respect to a perfect crystal. The LSS is based on a dot product of two bond orientational order complex vectors, with one vector associated with a particle in a liquid and the other vector with a particle in a crystal. The calculation of the LSS should scan the entire space of the Euler angles determined by the two coordinate frames describing individually the liquid and the crystal. The effectiveness of the LSS was examined by solid-like clusters in a Lennard-Jones (LJ) system near its liquid-solid phase transition and at solid states below its melting point, where the thermodynamic states of the LJ system were obtained by simulation annealing. The LSS measure was utilized to scrutinize the fcc-like, hcp-like, and bcc-like clusters classified by criteria based on W₄ and W₆ order parameters. As indicated by our results, the two ways of classification are consistent for fcc-like and hcp-like clusters, which are in a close resemblance to their crystalline counterparts. However, the classification with positive W₆ for bcc-like clusters is inconsistent with the results of the LSS measure, which was confirmed by clusters in a LJ system confined between two parallel slabs of particles in the bcc structure arrangement.

Validity of Effective Potentials in Crowded Solutions of Linear and Ring Polymers with Reversible Bonds

We perform simulations to compute the effective potential between the centers-of-mass of two polymers with reversible bonds. We investigate the influence of the topology on the potential by employing linear and ring backbones for the precursor (unbonded) polymer, finding that it leads to qualitatively different effective potentials. In the linear and ring cases the potentials can be described by Gaussians and generalized exponentials, respectively. The interactions are more repulsive for the ring topology, in analogy with known results in the absence of bonding. We also investigate the effect of the specific sequence of the reactive groups along the backbone (periodic or with different degrees of randomness), establishing that it has a significant impact on the effective potentials. When the reactive sites of both polymers are chemically orthogonal so that only intramolecular bonds are possible, the interactions become more repulsive the closer to periodic the sequence is. The opposite effect is found if both polymers have the same types of reactive sites and intermolecular bonds can be formed. We test the validity of the effective potentials in solution, in a broad range of concentrations from high dilution to far above the overlap concentration. For this purpose, we compare simulations of the effective fluid and test particle route calculations with simulations of the real all-monomer system. Very good agreement is found for the reversible linear polymers, indicating that unlike in their nonbonding counterparts many-body effects are minor even far above the overlap concentration. The agreement for the reversible rings is less satisfactory, and at high concentration the real system does not show the clustering behavior predicted by the effective potential. Results similar to the former ones are found for the partial self-correlations in ring/linear mixtures. Finally, we investigate the possibility of creating, at high concentrations, a gel of two interpenetrated reversible networks. For this purpose we simulate a 50/50 two-component mixture of reversible polymers with orthogonal chemistry for the reactive sites, so that intermolecular bonds are only formed between polymers of the same component. As predicted by both the theoretical phase diagram and the simulations of the effective fluid, the two networks in the all-monomer mixture do not interpenetrate, and phase separation (demixing) is observed instead.

Rheology Applied to Microgels: Brief (Revision of the) State of the Art

The ability of polymer microgels to rapidly respond to external stimuli is of great interest in sensors, lubricants, and biomedical applications, among others. In most of their uses, microgels are subjected to shear, deformation, and compression forces or a combination of them, leading to variations in their rheological properties. This review article mainly refers to the rheology of microgels, from the hard sphere versus soft particles' model. It clearly describes the scaling theories and fractal structure formation, in particular, the Shih et al. and Wu and Morbidelli models as a tool to determine the interactions among microgel particles and, thus, the viscoelastic properties. Additionally, the most recent advances on the characterization of microgels' single-particle interactions are also described. The review starts with the definition of microgels, and a brief introduction addresses the preparation and applications of microgels and hybrid microgels.

Glass quantization of the Gaussian core model

We use the replica method to study the dynamical glass transition of the Gaussian core model, a system of ultrasoft repulsive spheres interacting via a Gaussian potential, focusing on low temperatures and low-to-moderate densities. At constant temperature, an amorphous glassy state is entered upon a first compression but this glass melts as the density is further increased. In addition to this reentrant transition, a second, smooth transition is discovered between a continuous and a discretized glass. The properties of the former are continuous functions of temperatures, whereas the latter exhibits a succession of stripes, characterized by discontinuous jumps of the glassiness parameters. The glass physics of ultrasoft particles is hence richer than that of impenetrable particles for reasons that can be attributed to the ability of the former to create and break out-of-equilibrium clusters of overlapping particles.

Phase transitions of liposomes: when light meets heat

Phase transitions of liposomes are normally studied by differential scanning calorimetry. A suspension of liposomes is subjected to an increase (decrease) of temperature and when heat is absorbed (released), the liposomes transit from a gel (liquid) to a liquid (gel) phase. This endothermic (exothermic) process takes place at a temperature called the melting temperatureT_m, which is distinctive of the type of lipids forming the vesicles. The vesicles, though, also modify their size in the transition. Indeed, the thickness of the membranes decreases (increases) because carbon tails misalign (align). Concomitant with the modifications in the membrane thickness, the diameter (D) of the liposomes changes too. Therefore, when they are inspected by light, the scattered signal carries information from such dilatation (contraction) process. We performed careful experiments using dynamic light scattering as a function of temperature to detect the size changes of different liposomes. Gaussian fits of the derivatives of theDvsTcurves coincide within 1% with thermograms, which hints to the possibility of performing thermodynamic studies of lipid systems employing light.

Multiple 2D crystal structures in bilayered lamellae from direct self-assembly of 3D systems of soft Janus particles

Numerous crystals and Frank-Kasper phases in two-dimensional (2D) systems of soft particles have been presented by theoretical investigations. How to realize 2D crystals or Frank-kasper phases by direct self-assembly of...

Structural transition in two-dimensional Hertzian spheres in the presence of random pinning

Using molecular dynamics simulation we have investigated the influence of random pinning on the phase diagram and melting scenarios of a two-dimensional system with the Hertz potential for α=5/2. It has been shown that random pinning can cardinally change the mechanism of first-order transition between the different crystalline phases (triangular and square) by virtue of generating hexatic and tetratic phases: a triangular crystal to hexatic transition is of the continuous Berezinskii-Kosterlitz-Thouless (BKT) type, a hexatic to tetratic transition is of first order, and finally, there is a continuous BKT-type transition from tetratic to the square crystal.

Self-assembly of binary solutions to complex structures

Self-assembly in natural and synthetic molecular systems can create complex aggregates or materials whose properties and functionalities rise from their internal structure and molecular arrangement. The key microscopic features that control such assemblies remain poorly understood, nevertheless. Using classical density functional theory, we demonstrate how the intrinsic length scales and their interplay in terms of interspecies molecular interactions can be used to tune soft matter self-assembly. We apply our strategy to two different soft binary mixtures to create guidelines for tuning intermolecular interactions that lead to transitions from a fully miscible, liquid-like uniform state to formation of simple and core-shell aggregates and mixed aggregate structures. Furthermore, we demonstrate how the interspecies interactions and system composition can be used to control concentration gradients of component species within these assemblies. The insight generated by this work contributes toward understanding and controlling soft multi-component self-assembly systems. Additionally, our results aid in understanding complex biological assemblies and their function and provide tools to engineer molecular interactions in order to control polymeric and protein-based materials, pharmaceutical formulations, and nanoparticle assemblies.

General Synthesis of Ultrafine Monodispersed Hybrid Nanoparticles from Highly Stable Monomicelles

Ultrafine nanoparticles with organic-inorganic hybridization have essential roles in myriad applications. Over the past three decades, although various efforts on the formation of organic or inorganic ultrasmall nanoparticles have been made, ultrafine organic-inorganic hybrid nanoparticles have scarcely been achieved. Herein, a family of ultrasmall hybrid nanoparticles with a monodisperse, uniform size is synthesized by a facile thermo-kinetics-mediated copolymer monomicelle approach. These thermo-kinetics-mediated monomicelles with amphiphilic ABC triblock copolymers are structurally robust due to their solidified polystyrene core, endowing them with ultrahigh thermodynamic stability, which is difficult to achieve using Pluronic surfactant-based micelles (e.g., F127). This great stability combined with a core-shell-corona structure makes the monodispersed monomicelles a robust template for the precise synthesis of ultrasmall hybrid nanoparticles with a highly uniform size. As a demonstration, the obtained micelles/SiO₂ hybrid nanoparticles display ultrafine sizes, excellent uniformity, monodispersity, and tunable structural parameters (diameters: 24-47 nm and thin shell thickness: 2.0-4.0 nm). Notably, this approach is universal for creating a variety of multifunctional ultrasmall hybrid nanostructures, involving organic/organic micelle/polymers (polydopamine) nanoparticles, organic/inorganic micelle/metal oxides (ZnO, TiO₂ , Fe₂ O₃ ), micelle/hydroxides (Co(OH)₂ ), micelle/noble metals (Ag), and micelle/TiO₂ /SiO₂ hybrid composites. As a proof of concept, the ultrasmall micelle/SiO₂ hybrid nanoparticles demonstrate superior toughness as biomimetic materials.

From predictive modelling to machine learning and reverse engineering of colloidal self-assembly

An overwhelming diversity of colloidal building blocks with distinct sizes, materials and tunable interaction potentials are now available for colloidal self-assembly. The application space for materials composed of these building blocks is vast. To make progress in the rational design of new self-assembled materials, it is desirable to guide the experimental synthesis efforts by computational modelling. Here, we discuss computer simulation methods and strategies used for the design of soft materials created through bottom-up self-assembly of colloids and nanoparticles. We describe simulation techniques for investigating the self-assembly behaviour of colloidal suspensions, including crystal structure prediction methods, phase diagram calculations and enhanced sampling techniques, as well as their limitations. We also discuss the recent surge of interest in machine learning and reverse-engineering methods. Although their implementation in the colloidal realm is still in its infancy, we anticipate that these data-science tools offer new paradigms in understanding, predicting and (inverse) design of novel colloidal materials.

Volume fraction determination of microgel composed of interpenetrating polymer networks of PNIPAM and polyacrylic acid

Interpenetrated polymer network microgels, composed of crosslinked networks of poly(N-isopropylacrylamide) and polyacrylic acid (PAAc), have been investigated through rheological measurements at four different amounts of PAAc. Both PAAc content and crosslinking degree modify particle dimensions, mass and softness, thereby strongly affecting the volume fraction and the system viscosity. Here the volume fraction is derived from the flow curves at low concentrations by fitting the zero-shear viscosity with the Einstein-Batchelor equation which provides a parameterkto shift weight concentration to volume fraction. We find that particles with higher PAAc content and crosslinker are characterized by a greater value ofkand therefore by larger volume fractions when compared to softer particles. The packing fractions obtained from rheological measurements are compared with those from static light scattering for two PAAc contents revealing a good agreement. Moreover, the behaviour of the viscosity as a function of packing fraction, at room temperature, has highlighted an Arrhenius dependence for microgels synthesized with low PAAc content and a Vogel-Fulcher-Tammann dependence for the highest investigated PAAc concentration. A comparison with the hard spheres behaviour indicates a steepest increase of the viscosity with decreasing particles softness. Finally, the volume fraction dependence of the viscosity at a fixed PAAc and at two different temperatures, below and above the volume phase transition, shows a quantitative agreement with the structural relaxation time measured through dynamic light scattering indicating that interpenetrated polymer network microgels softness can be tuned with PAAc and temperature and that, depending on particle softness, two different routes are followed.

Theoretical and Experimental Design of Heavy Metal-Mopping Magnetic Nanoparticles

Herein, we show a comprehensive experimental, theoretical, and computational study aimed at designing macromolecules able to adsorb a cargo at the nanoscale. Specifically, we focus on the adsorption properties of star diblock copolymers, i.e., macromolecules made by a number f of H-T diblock copolymer arms tethered on a central core; the H monomeric heads, which are closer to the tethering point, are attractive toward a specific target, while the T monomeric tails are neutral to the cargo. Experimentally, we exploited the adaptability of poly(2-oxazoline)s (POxs) to realize block copolymer-coated nanoparticles with a proper functionalization able to interact with heavy metals and show or exhibit a thermoresponsive behavior in aqueous solution. We here present the synthesis and analysis of the properties of a high molecular mass block copolymer featured by (i) a polar side chain, capable of exploiting electrostatic and hydrophilic interaction with a predetermined cargo, and (ii) a thermoresponsive scaffold, able to change the interaction with the media by tuning the temperature. Afterward, the obtained polymers were grafted onto iron oxide nanoparticles and the thermoresponsive properties were investigated. Through isothermal titration calorimetry, we then analyzed the adsorption properties of the synthesized superparamagnetic nanoparticles for heavy metal ions in aqueous solution. Additionally, we use a combination of scaling theories and simulations to link equilibrium properties of the system to a prediction of the loading properties as a function of size ratio and effective interactions between the considered species. The comparison between experimental results on adsorption and theoretical prediction validates the whole design process.

Structure of the simple harmonic-repulsive system in liquid and glassy states studied by the triple correlation function

An efficient description of the structures of liquids and, in particular, the structural changes that happen with liquids on supercooling remains to be a challenge. The systems composed of soft particles are especially interesting in this context because they often demonstrate non-trivial local orders that do not allow to introduce the concept of the nearest-neighbor shell. For this reason, the use of some methods, developed for the structure analysis of atomic liquids, is questionable for the soft-particle systems. Here we report about our investigations of the structure of the simple harmonic-repulsive liquid in 3D using the triple correlation function (TCF), i.e., the method that does not rely on the nearest neighbor concept. The liquid is considered at reduced pressure (P = 1.8) at which it exhibits remarkable stability against crystallization on cooling. It is demonstrated that the TCF allows addressing the development of the orientational correlations in the structures that do not allow drawing definite conclusions from the studies of the bond-orientational order parameters. Our results demonstrate that the orientational correlations, if measured by the heights of the peaks in the TCF, significantly increase on cooling. This rise in the orientational ordering is not captured properly by the Kirkwood's superposition approximation. Detailed considerations of the peaks' shapes in the TCF suggest the existence of a link between the orientational ordering and the slowdown of the system's dynamics. Our findings support the view that the development of the orientational correlations in liquids may play a significant role in the liquids' dynamics and that the considerations of the pair distribution function may not be sufficient to understand intuitively all the structural changes that happen with liquids on supercooling. In general, our results demonstrate that the considerations of the TCF are useful in the discussions of the liquid's structures beyond the pair density function and interpreting the results obtained with the bond-orientational order parameters.

Reentrant melting and multiple occupancy crystals of bounded potentials: Simple theory and direct observation by molecular dynamics simulations

Aspects of the phase coexistence behavior of the generalized exponential model (GEM-m) and bounded versions of inverse power potentials based on theory and molecular dynamics (MD) simulation data are reported. The GEM-m potential is ϕ(r)=exp(-r^{m}), where r is the pair separation and m is an adjustable exponent. A simple analytic formula for the fluid-solid envelope of the Gaussian core model which takes account of the known low- and high-density limiting forms is proposed and shown to represent the simulation data well. The bounded inverse power (BIP) potential is ϕ(r)=1/(a^{q}+r^{q})^{n/q}, where a, n, and q are positive constants. The BIP potential multiple occupancy crystal or cluster crystals are predicted to form when q>2 and a>0, for n>3, which compares with the corresponding GEM-m condition of m>2. Reentrant melting should occur for the BIP potential when q≤2 and a>0. MD simulations in which the system was gradually compressed at constant temperature using the BIP potential produced cluster states in the parameter domain expected but it was not possible to establish conclusively whether a multiply occupied crystal or a cluster fluid had formed owing to assembly structural fluctuations. The random phase approximation reproduces very well the BIP MD energy per particle without any discontinuities at the phase boundaries. The Lindemann melting rule is shown analytically to give a more rapidly decaying reentrant melting curve boundary than the so-called melting indicator (MI) empirical melting criterion which has also been investigated in this study. The MI model gives a better match to the high-density phase boundary for small m and q values.

Softness- and Size-Dependent Packing Symmetries of Polymer-Grafted Nanoparticles

Achieving ordered arrays of nanoparticles (NPs) with controlled packing symmetry and interparticle spacing is of great importance to design complex metamaterials. Herein, we report softness- and size-dependent self-assembly behavior of polystyrene-grafted Au NPs (Au@PS NPs). We varied the core size of Au NPs from 1.9 to 9.6 nm and the number-average molecular weight (M_n) of thiol-terminated polystyrene from 1.8 to 7.9 kg mol^-1. The optimal packing model based on an "effective softness" parameter λ_eff that accounts for close-packed and semidilute brush regimes could predict the effective radius of Au@PS NPs (within ±9%) for a wide range of PS M_n, grafting density, and Au core size. With increasing λ_eff, the self-assembled Au@PS NP superlattices undergo a symmetry transition from hexagonal close packed (hcp) to body-centered tetragonal (bct) to body-centered cubic (bcc). This work demonstrates the effective softness model as a simple but robust tool for the design of NP superlattices with precisely controlled interparticle distance and packing symmetry, both of which are critical for the development of sophisticated materials through control of nanoscale structure.

Stress driven fractionalization of vacancies in regular packings of elastic particles

Elucidating the interplay of defects and stress at the microscopic level is a fundamental physical problem that has a strong connection with materials science. Here, based on the two-dimensional crystal model, we show that the instability mode of vacancies with varying size and morphology conforms to a common scenario. A vacancy under compression is fissioned into a pair of dislocations that glide and vanish at the boundary. This neat process is triggered by the local shear stress around the vacancy. The remarkable fractionalization of vacancies creates rich modes of interaction between vacancies and other topological defects, and provides a new dimension for mechanical engineering of defects in extensive crystalline structures.

Stress-induced ordering of two-dimensional packings of elastic spheres

Packing of particles in confined environments is a common problem in multiple fields. Here, based on the two-dimensional Hertzian particle model, we study the packing of deformable spherical particles under compression and reveal the crucial role of stress as an ordering field in regulating particle arrangement. Specifically, under increasing compression, the squeezed particles spontaneously organize into regular packings in the sequence of triangular and square lattices, pentagonal tessellation, and the reentrant triangular lattice. The rich ordered patterns and complex structures revealed in this work suggest a fruitful organizational strategy based on the interplay of external stress and intrinsic elastic instability of particle arrays.

Melting scenarios of two-dimensional Hertzian spheres with a single triangular lattice

We present a molecular dynamics simulation study of the phase diagram and melting scenarios of two-dimensional Hertzian spheres with exponent 7/2. We have found multiple re-entrant melting of a single crystal with a triangular lattice in a wide range of densities from 0.5 to 10.0. Depending on the position on the phase diagram, the triangular crystal has been shown to melt through both two-stage melting with a first-order hexatic-isotropic liquid transition and a continuous solid-hexatic transition as well as in accordance with the Berezinskii-Kosterlitz-Thouless-Halperin-Nelson-Young (BKTHNY) scenario (two continuous transitions with an intermediate hexatic phase). We studied the behavior of heat capacity and have shown that despite two-stage melting, the heat capacity has one peak which seems to correspond to a solid-hexatic transition.

Understanding, predicting, and tuning the fragility of vitrimeric polymers

Fragility is an empirical property that describes how abruptly a glass-forming material solidifies upon supercooling. The degree of fragility carries important implications for the functionality and processability of a material, as well as for our fundamental understanding of the glass transition. However, the microstructural properties underlying fragility still remain poorly understood. Here, we explain the microstructure-fragility link in vitrimeric networks, a novel type of high-performance polymers with unique bond-swapping functionality and unusual glass-forming behavior. Our results are gained from coarse-grained computer simulations and first-principles mode-coupling theory (MCT) of star-polymer vitrimers. We first demonstrate that the vitrimer fragility can be tuned over an unprecedentedly broad range, from fragile to strong and even superstrong behavior, by decreasing the bulk density. Remarkably, this entire phenomenology can be reproduced by microscopic MCT, thus challenging the conventional belief that existing first-principles theories cannot account for nonfragile behaviors. Our MCT analysis allows us to rationally identify the microstructural origin of the fragile-to-superstrong crossover, which is rooted in the sensitivity of the static structure factor to temperature variations. On the molecular scale, this behavior stems from a change in dominant length scales, switching from repulsive excluded-volume interactions to intrachain attractions as the vitrimer density decreases. Finally, we develop a simplified schematic MCT model which corroborates our microscopically founded conclusions and which unites our findings with earlier MCT studies. Our work sheds additional light on the elusive structure-fragility link in glass-forming matter and provides a first-principles-based platform for designing amorphous materials with an on-demand dynamic response.

Anomalous behavior and structure of a liquid of particles interacting through the harmonic-repulsive pair potential near the crystallization transition

A characteristic property of many soft matter systems is an ultrasoft effective interaction between their structural units. This softness often leads to complex behavior. In particular, ultrasoft systems under pressure demonstrate polymorphism of complex crystal and quasicrystal structures. Therefore, it is of interest to investigate how different can be the structure of the fluid state in such systems at different pressures. Here we address this issue for a model liquid composed of particles interacting through the harmonic-repulsive pair potential. This system can form different crystal structures as the liquid is cooled. We find that, at certain pressures, the liquid exhibits unusual properties, such as a negative thermal expansion coefficient. Besides, the volume and the potential energy of the system can increase during crystallization. At certain pressures, the system demonstrates high stability against crystallization and it is hardly possible to crystallize it on the timescales of the simulations. To address the liquid's structure at high pressures, we consider the scaled pair distribution function (PDF) and the bond-orientational order (BOO) parameters. The marked change happening with the PDF, as pressure increases, is the splitting of the first peak which is caused by the appearance of non-negligible interactions with the second neighbors and the following rearrangement of the structure. Our findings suggest that non-trivial effects, usually explained by different interactions at different spatial scales, can also be observed in one-component systems with simple one-length-scale ultrasoft repulsive interactions.

Monte Carlo simulation and integral equation study of Hertzian spheres in the low-temperature regime

We investigate the behavior of Hertzian spheres in the fluid phase and in proximity of the freezing threshold by using Monte Carlo (MC) simulations and integral equation theories, based on the Ornstein-Zernike (OZ) approach. The study is motivated by the importance of the Hertzian model in representing a large class of systems interacting via soft interactions, such as star polymers or microgels. Radial distribution functions, structure factors, and excess entropy clearly show the reentrant behavior typical of the Hertzian fluid, well caught by both MC simulations and OZ theory. Then, we make use of some phenomenological one-phase criteria for testing their reliability in detecting the freezing threshold. All criteria provide evidence of the fluid-solid transition with different degrees of accuracy. This suggests the possibility to adopt these empirical rules to provide a quick and reasonable estimate of the freezing transition in model potentials of wide interest for soft matter systems.

Multi-particle collision dynamics for a coarse-grained model of soft colloids

The growing interest in the dynamical properties of colloidal suspensions, both in equilibrium and under an external drive such as shear or pressure flow, requires the development of accurate methods to correctly include hydrodynamic effects due to the suspension in a solvent. In the present work, we generalize Multiparticle Collision Dynamics (MPCD) to be able to deal with soft, polymeric colloids. Our methods build on the knowledge of the monomer density profile that can be obtained from monomer-resolved simulations without hydrodynamics or from theoretical arguments. We hereby propose two different approaches. The first one simply extends the MPCD method by including in the simulations effective monomers with a given density profile, thus neglecting monomer-monomer interactions. The second one considers the macromolecule as a single penetrable soft colloid (PSC), which is permeated by an inhomogeneous distribution of solvent particles. By defining an appropriate set of rules to control the collision events between the solvent and the soft colloid, both linear and angular momenta are exchanged. We apply these methods to the case of linear chains and star polymers for varying monomer lengths and arm number, respectively, and compare the results for the dynamical properties with those obtained within monomer-resolved simulations. We find that the effective monomer method works well for linear chains, while the PSC method provides very good results for stars. These methods pave the way to extend MPCD treatments to complex macromolecular objects such as microgels or dendrimers and to work with soft colloids at finite concentrations.

Deriving phase field crystal theory from dynamical density functional theory: Consequences of the approximations

Phase field crystal (PFC) theory is extensively used for modeling the phase behavior, structure, thermodynamics, and other related properties of solids. PFC theory can be derived from dynamical density functional theory (DDFT) via a sequence of approximations. Here, we carefully identify all of these approximations and explain the consequences of each. One approximation that is made in standard derivations is to neglect a term of form ∇·[n∇Ln], where n is the scaled density profile and L is a linear operator. We show that this term makes a significant contribution to the stability of the crystal, and that dropping this term from the theory forces another approximation, that of replacing the logarithmic term from the ideal gas contribution to the free energy with its truncated Taylor expansion, to yield a polynomial in n. However, the consequences of doing this are (i) the presence of an additional spinodal in the phase diagram, so the liquid is predicted first to freeze and then to melt again as the density is increased; and (ii) other periodic structures, such as stripes, are erroneously predicted to be thermodynamic equilibrium structures. In general, L consists of a nonlocal convolution involving the pair direct correlation function. A second approximation sometimes made in deriving PFC theory is to replace L with a gradient expansion involving derivatives. We show that this leads to the possibility of the density going to zero, with its logarithm going to -∞ while being balanced by the fourth derivative of the density going to +∞. This subtle singularity leads to solutions failing to exist above a certain value of the average density. We illustrate all of these conclusions with results for a particularly simple model two-dimensional fluid, the generalized exponential model of index 4 (GEM-4), chosen because a DDFT is known to be accurate for this model. The consequences of the subsequent PFC approximations can then be examined. These include the phase diagram being both qualitatively incorrect, in that it has a stripe phase, and quantitatively incorrect (by orders of magnitude) regarding the properties of the crystal phase. Thus, although PFC models are very successful as phenomenological models of crystallization, we find it impossible to derive the PFC as a theory for the (scaled) density distribution when starting from an accurate DDFT, without introducing spurious artifacts. However, we find that making a simple one-mode approximation for the logarithm of the density distribution lnρ(x) rather than for ρ(x) is surprisingly accurate. This approach gives a tantalizing hint that accurate PFC-type theories may instead be derived as theories for the field lnρ(x), rather than for the density profile itself.

Controlled self-aggregation of polymer-based nanoparticles employing shear flow and magnetic fields

Star polymers with magnetically functionalized end groups are presented as a novel polymeric system whose morphology, self-aggregation, and orientation can easily be tuned by exposing these macromolecules simultaneously to an external magnetic field and to shear forces. Our investigations are based on a specialized simulation technique which faithfully takes into account the hydrodynamic interactions of the surrounding, Newtonian solvent. We find that the combination of magnetic field (including both strength and direction) and shear rate controls the mean number of magnetic clusters, which in turn is largely responsible for the static and dynamic behavior. While some properties are similar to comparable non-magnetic star polymers, others exhibit novel phenomena; examples of the latter include the breakup and reorganization of the clusters beyond a critical shear rate, and a strong dependence of the efficiency with which shear rate is translated into whole-body rotations on the direction of the magnetic field.

Influence of Branching on the Configurational and Dynamical Properties of Entangled Polymer Melts

We probe the influence of branching on the configurational, packing, and density correlation function properties of polymer melts of linear and star polymers, with emphasis on molecular masses larger than the entanglement molecular mass of linear chains. In particular, we calculate the conformational properties of these polymers, such as the hydrodynamic radius R h , packing length p, pair correlation function g ( r ) , and polymer center of mass self-diffusion coefficient, D, with the use of coarse-grained molecular dynamics simulations. Our simulation results reproduce the phenomenology of simulated linear and branched polymers, and we attempt to understand our observations based on a combination of hydrodynamic and thermodynamic modeling. We introduce a model of "entanglement" phenomenon in high molecular mass polymers that assumes polymers can viewed in a coarse-grained sense as "soft" particles and, correspondingly, we model the emergence of heterogeneous dynamics in polymeric glass-forming liquids to occur in a fashion similar to glass-forming liquids in which the molecules have soft repulsive interactions. Based on this novel perspective of polymer melt dynamics, we propose a functional form for D that can describe our simulation results for both star and linear polymers, covering both the unentangled to entangled polymer melt regimes.

Microgels Adsorbed at Liquid-Liquid Interfaces: A Joint Numerical and Experimental Study

Soft particles display highly versatile properties with respect to hard colloids and even more so at fluid-fluid interfaces. In particular, microgels, consisting of a cross-linked polymer network, are able to deform and flatten upon adsorption at the interface due to the balance between surface tension and internal elasticity. Despite the existence of experimental results, a detailed theoretical understanding of this phenomenon is still lacking due to the absence of appropriate microscopic models. In this work, we propose an advanced modeling of microgels at a flat water/oil interface. The model builds on a realistic description of the internal polymeric architecture and single-particle properties of the microgel and is able to reproduce its experimentally observed shape at the interface. Complementing molecular dynamics simulations with in situ cryo-electron microscopy experiments and atomic force microscopy imaging after Langmuir-Blodgett deposition, we compare the morphology of the microgels for different values of the cross-linking ratios. Our model allows for a systematic microscopic investigation of soft particles at fluid interfaces, which is essential to develop predictive power for the use of microgels in a broad range of applications, including the stabilization of smart emulsions and the versatile patterning of surfaces.

Slow dynamics coupled with cluster formation in ultrasoft-potential glasses

We numerically investigate the slow dynamics of a binary mixture of ultrasoft particles interacting with the generalized Hertzian potential. If the softness parameter, α, is small, the particles at high densities start penetrating each other, form clusters, and eventually undergo the glass transition. We find multiple cluster-glass phases characterized by a different number of particles per cluster, whose boundary lines are sharply separated by the cluster size. Anomalous logarithmic slow relaxation of the density correlation functions is observed in the vicinity of these glass-glass phase boundaries, which hints the existence of the higher-order dynamical singularities predicted by the mode-coupling theory. Deeply in the cluster glass phases, it is found that the dynamics of a single particle is decoupled from that of the collective fluctuations.

Structure and stimuli-responsiveness of all-DNA dendrimers: theory and experiment

We present a comprehensive theoretical and experimental study of the solution phase properties of a DNA-based family of nanoparticles - dendrimer-like DNA molecules (DL-DNA). These charged DNA dendrimers are novel macromolecular aggregates, which hold high promise in targeted self-assembly of soft matter systems in the bulk and at interfaces. To describe the behaviour of this family of dendrimers (with generations ranging from G1 to G7), we use a theoretical model in which base-pairs of a single DL-DNA molecule are modeled by charged monomers, whose interactions are chosen to mimic the equilibrium properties of DNA correctly. Experimental results on the sizes and conformations of DL-DNA are based on static and dynamic light scattering; and molecular dynamics simulations are employed to model the equilibrium properties of DL-DNA, which compare favorably to the findings from experiments while at the same time providing a host of additional information and insight into the molecular structure of the nanostructures. We also examine the salt-responsiveness of these macromolecules, finding that despite the strong screening of electrostatic interactions brought about by the added salt, the macromolecules shrink only slightly, their size robustness stemming from the high bending rigidity of the DNA-segments. The study of these charged dendrimer systems is an important field of research in the area of soft matter due to their potential role for various interdisciplinary applications, ranging from molecular cages and carriers for drug delivery in a living organism to the development of dendrimer- and dendron-based ultra-thin films in the area of nanotechnology. These findings are essential to determine if DL-DNA is a viable candidate for the experimental realization of cluster crystals in the bulk, a novel form of solid with multiple site occupancy.