Formation of polymorphic cluster phases for a class of models of purely repulsive soft spheres

Synergizing Steric Hindrance and Stacking Interactions To Facilitate the Controlled Assembly of Multiple 41 Metalla-Knots and Pseudo-Solomon Links

In this work, a noncoplanar terphenyl served as a building block to synthesize a novel 3,3'-substituted bipyridyl ligand (L1) which further reacted with binuclear half-sandwich units A/B, giving rise to two aesthetic 41 metalla-knots in high yields via a coordination-driven self-assembly strategy. Furthermore, given the inherent compactness of the 41 metalla-knots, it creates favorable conditions for the emergence of steric repulsion. We focused on progressively introducing nitrogen atoms featuring a lone pair of electrons (LPEs) into ligand L1 to manipulate the balance of H⋅⋅⋅H/LPEs⋅⋅⋅LPEs steric repulsion during the assembly process, ultimately achieving controlled assembly from 41 metalla-knots to the pseudo-Solomon link and then to molecular tweezer-like assembly facilitated by stacking interactions. All the assemblies were well characterized by solution-state NMR techniques, ESI-TOF/MS, and single-crystal X-ray diffraction. The evolutionary process of the topological architectures is equivalent to visualizing the synergistic effect of steric hindrance and stacking interactions on structural assembly, providing a new avenue for achieving the controlled synthesis of different topologies.

A quantitative theory and atomistic simulation study on the soft-sphere crystal-melt interfacial properties. I. Kinetic coefficients

By employing non-equilibrium molecular dynamics (NEMD) simulations and time-dependent Ginzburg-Landau (TDGL) theory for solidification kinetics [Cryst. Growth Des. 20, 7862 (2020)], we predict the kinetic coefficients of FCC(100) crystal-melt interface (CMI) of soft-spheres modeled with an inverse-sixth-power repulsive potential. The collective dynamics of the local interfacial liquid phase at the equilibrium FCC(100) CMIs are calculated based on a recently proposed algorithm [J. Chem. Phys. 157, 084 709 (2022)] and are employed as the resulting parameter that eliminates the discrepancy between the predictions of the kinetic coefficient using the NEMD simulations and the TDGL solidification theory. A speedup of the two modes of the interfacial liquid collective dynamics (at wavenumbers equal to the principal and the secondary reciprocal lattice vector of the grown crystal) is observed. With the insights provided by the quantitative predictive theory, the variation of the solidification kinetic coefficient along the crystal-melt coexistence boundary is discussed. The combined methodology (simulation and theory) presented in this study could be further applied to investigate the role of the inter-atomic potential (e.g., softness parameter s = 1/n of the inverse-power repulsive potential) in the kinetic coefficient.

Phase separation dynamics in a symmetric binary mixture of ultrasoft particles

Phase separation plays a key role in determining the self-assembly of biological and soft-matter systems. In biological systems, liquid-liquid phase separation inside a cell leads to the formation of various macromolecular aggregates. The interaction among these aggregates is soft, i.e., they can significantly overlap at a small energy cost. From a computer simulation point of view, these complex macromolecular aggregates are generally modeled by soft particles. The effective interaction between two particles is defined via the generalized exponential model of index n, with n = 4. Here, using molecular dynamics simulations, we study the phase separation dynamics of a size-symmetric binary mixture of ultrasoft particles. We find that when the mixture is quenched to a temperature below the critical temperature, the two components spontaneously start to separate. Domains of the two components form, and the equal-time order parameter reveals that the domain sizes grow with time in a power-law manner with an exponent of 1/3, which is consistent with the Lifshitz-Slyozov law for conserved systems. Furthermore, the static structure factor shows a power-law decay with an exponent of 4, consistent with the Porod law.

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.

Structural transition, orientational order, and anomalous specific heat in a two-dimensional dimer crystal of core-softened particles

Systems featuring hard-core-soft-shell repulsive pair potentials can form ordered phases, where particles organize themselves in aggregates with nontrivial geometries. The dimer crystal formed by one such potential, namely, the hard-core plus generalized exponential model of order 4, has been recently investigated, revealing a low-temperature structural phase transition, with the onset of nematic ordering of the dimers. In the present work, we aim to characterize this phase transition via a mean-field theory, by which a detailed analysis of the low-temperature properties of the system is carried out under quadrupole approximation. We determine the transition temperature and identify its order parameter, highlighting the link between the structural transition and the nematic ordering of the system. The first-order character of the transition is established and supported by the Landau expansion of the free energy in powers of the order parameter. The theory is subsequently generalized to take into account lattice vibrations and dimer length fluctuations. Finally, we provide an explanation for the anomalous behavior displayed by the specific heat in the vanishing-temperature limit, which is also supported by Monte Carlo simulations.

Cluster Formation in Solutions of Polyelectrolyte Rings

We use molecular dynamics simulations to explore concentrated solutions of semiflexible polyelectrolyte ring polymers, akin to the DNA mini-circles, with counterions of different valences. We find that the assembly of rings into nanoscopic cylindrical stacks is a generic feature of the systems, but the morphology and dynamics of such a cluster can be steered by the counterion conditions. In general, a small addition of trivalent ions can stabilize the emergence of clusters due to the counterion condensation, which mitigates the repulsion between the like-charged rings. Stoichiometric addition of trivalent ions can even lead to phase separation of the polyelectrolyte ring phase due to the ion-bridging effects promoting otherwise entropically driven clustering. On the other hand, monovalent counterions cause the formation of stacks to be re-entrant with density. The clusters are stable within a certain window of concentration, while above the window the polyelectrolytes undergo an osmotic collapse, disfavoring ordering. The cluster phase exhibits characteristic cluster glass dynamics with arrest of collective degrees of freedom but not the self-ones. On the other hand, the collapsed phase shows arrest on both the collective and single level, suggesting an incipient glass-to-glass transition, from a cluster glass of ring clusters to a simple glass of rings.

Mixing Linear Polymers with Rings and Catenanes: Bulk and Interfacial Behavior

We derive and parameterize effective interaction potentials between a multitude of different types of ring polymers and linear chains, varying the bending rigidity and solvent quality for the former species. We further develop and apply a density functional treatment for mixtures of both disconnected (chain-ring) and connected (chain-polycatenane) mixtures of the same, drawing coexistence binodals and exploring the ensuing response functions as well as the interface and wetting behavior of the mixtures. We show that worsening of the solvent quality for the rings brings about a stronger propensity for macroscopic phase separation in the linear-polycatenane mixtures, which is predominantly of the demixing type between phases of similar overall particle density. We formulate a simple criterion based on the effective interactions, allowing us to determine whether any specific linear-ring mixture will undergo a demixing phase separation.

A computational framework for testing hypotheses of the minimal mechanical requirements for cell aggregation using early annual killifish embryogenesis as a model

Introduction: Deciphering the biological and physical requirements for the outset of multicellularity is limited to few experimental models. The early embryonic development of annual killifish represents an almost unique opportunity to investigate de novo cellular aggregation in a vertebrate model. As an adaptation to seasonal drought, annual killifish employs a unique developmental pattern in which embryogenesis occurs only after undifferentiated embryonic cells have completed epiboly and dispersed in low density on the egg surface. Therefore, the first stage of embryogenesis requires the congregation of embryonic cells at one pole of the egg to form a single aggregate that later gives rise to the embryo proper. This unique process presents an opportunity to dissect the self-organizing principles involved in early organization of embryonic stem cells. Indeed, the physical and biological processes required to form the aggregate of embryonic cells are currently unknown. Methods: Here, we developed an in silico, agent-based biophysical model that allows testing how cell-specific and environmental properties could determine the aggregation dynamics of early Killifish embryogenesis. In a forward engineering approach, we then proceeded to test two hypotheses for cell aggregation (cell-autonomous and a simple taxis model) as a proof of concept of modeling feasibility. In a first approach (cell autonomous system), we considered how intrinsic biophysical properties of the cells such as motility, polarity, density, and the interplay between cell adhesion and contact inhibition of locomotion drive cell aggregation into self-organized clusters. Second, we included guidance of cell migration through a simple taxis mechanism to resemble the activity of an organizing center found in several developmental models. Results: Our numerical simulations showed that random migration combined with low cell-cell adhesion is sufficient to maintain cells in dispersion and that aggregation can indeed arise spontaneously under a limited set of conditions, but, without environmental guidance, the dynamics and resulting structures do not recapitulate in vivo observations. Discussion: Thus, an environmental guidance cue seems to be required for correct execution of early aggregation in early killifish development. However, the nature of this cue (e.g., chemical or mechanical) can only be determined experimentally. Our model provides a predictive tool that could be used to better characterize the process and, importantly, to design informed experimental strategies.

Continuum mechanics of nonideal crystals: Microscopic approach based on projection-operator formalism

We present a microscopic derivation of the laws of continuum mechanics of nonideal ordered solids including dissipation, defect diffusion, and heat transport. The starting point is the classical many-body Hamiltonian. The approach relies on the Zwanzig-Mori projection operator formalism to connect microscopic fluctuations to thermodynamic derivatives and transport coefficients. Conservation laws and spontaneous symmetry breaking, implemented via Bogoliubov's inequality, determine the selection of the slow variables. Density fluctuations in reciprocal space encode the displacement field and the defect concentration. Isothermal and adiabatic elastic constants are obtained from equilibrium correlations, while transport coefficients are given as Green-Kubo formulas, providing the basis for their measurement in atomistic simulations or colloidal experiments. The approach to the linearized continuum mechanics and results are compared to others from the literature.

Ultralong organic phosphorescence from isolated molecules with repulsive interactions for multifunctional applications

Intermolecular interactions, including attractive and repulsive interactions, play a vital role in manipulating functionalization of the materials from micro to macro dimensions. Despite great success in generation of ultralong organic phosphorescence (UOP) by suppressing non-radiative transitions through attractive interactions recently, there is still no consideration of repulsive interactions on UOP. Herein, we proposed a feasible approach by introducing carboxyl groups into organic phosphors, enabling formation of the intense repulsive interactions between the isolated molecules and the matrix in rigid environment. Our experimental results show a phosphor with a record lifetime and quantum efficiency up to 3.16 s and 50.0% simultaneously in film under ambient conditions. Considering the multiple functions of the flexible films, the potential applications in anti-counterfeiting, afterglow display and visual frequency indicators were demonstrated. This finding not only outlines a fundamental principle to achieve bright organic phosphorescence in film, but also expands the potential applications of UOP materials.

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.

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 diagram of SALR fluids on spherical surfaces

We investigate the phase diagram of a fluid of hard-core disks confined to the surface of a sphere and whose interaction potential contains a short-range attraction followed by a long-range repulsive tail (SALR). Based on previous work in the bulk we derive a stability criterion for the homogeneous phase of the fluid, and locate a region of instability linked to the presence of a negative minimum in the spherical harmonics expansion of the interaction potential. The inhomogeneous phases contained within this region are characterized using a mean-field density functional theory. We find several inhomogeneous patterns that can be separated into three broad classes: cluster crystals, stripes, and bubble crystals, each containing topological defects. Interestingly, while the periodicity of inhomogeneous phases at large densities is mainly determined by the position of the negative minimum of the potential expansion, the finite size of the system induces a richer behavior at smaller densities.

Self assembling cluster crystals from DNA based dendritic nanostructures

Cluster crystals are periodic structures with lattice sites occupied by several, overlapping building blocks, featuring fluctuating site occupancy, whose expectation value depends on thermodynamic conditions. Their assembly from atomic or mesoscopic units is long-sought-after, but its experimental realization still remains elusive. Here, we show the existence of well-controlled soft matter cluster crystals. We fabricate dendritic-linear-dendritic triblock composed of a thermosensitive water-soluble polymer and nanometer-scale all-DNA dendrons of the first and second generation. Conclusive small-angle X-ray scattering (SAXS) evidence reveals that solutions of these triblock at sufficiently high concentrations undergo a reversible phase transition from a cluster fluid to a body-centered cubic (BCC) cluster crystal with density-independent lattice spacing, through alteration of temperature. Moreover, a rich concentration-temperature phase diagram demonstrates the emergence of various ordered nanostructures, including BCC cluster crystals, birefringent cluster crystals, as well as hexagonal phases and cluster glass-like kinetically arrested states at high densities.

Experimental realization of cluster crystals- periodic structures with lattice sites occupied by several, overlapping building blocks, has been elusive. Here, the authors show the existence of well-controlled soft matter cluster crystals composed of a thermosensitive water-soluble polymer and nanometer-scale all-DNA dendrons.

Low-temperature ordering of the dimer phase of a two-dimensional model of core-softened particles

Purely pairwise interactions of the core-softened type, i.e., featuring a soft repulsion followed by a hard-core interaction at shorter distance, give rise to nontrivial equilibrium structures entirely different from the standard close packing of spheres. In particular, in a suitable low-temperature region of their phase diagram, such interactions are well known to favor a transition from a fluid to a cluster crystal. The residual mutual interaction between individual clusters can lead to the formation of patterns of their reciprocal orientations. In this work, we investigate two examples of such models in two dimensions, at the density most appropriate to the dimer phase, whereby clusters consist of just two particles, studying them with optimization techniques and Monte Carlo simulations. We focus on the dimer crystal, and unveil a second phase transition at extremely low temperature. This transition leads from a triangular dimer lattice with randomly disordered dimer orientations at high temperature to a reduced-symmetry ground state with nematic orientational order and a slightly distorted structure characterized by a centered-rectangular lattice at low temperature.

On the yielding of a point-defect-rich model crystal under shear: insights from molecular dynamics simulations

In real crystals and at finite temperatures point defects are inevitable. Under shear their dynamics severely influence the mechanical properties of these crystals, giving rise to non-linear effects, such as ductility. In an effort to elucidate the complex behavior of crystals under plastic deformation it is crucial to explore and to understand the interplay between the timescale related to the equilibrium point-defect diffusion and the shear-induced timescale. Based on extensive non-equilibrium molecular dynamics simulations we present a detailed investigation on the yielding behavior of cluster crystals, an archetypical model for a defect-rich crystal: in such a system clusters of overlapping particles occupy the lattice sites of a regular (FCC) structure. In equilibrium particles diffuse via site-to-site hopping while maintaining the crystalline structure intact. We investigate these cluster crystals at a fixed density and at different temperatures where the system remains in the FCC structure: temperature allows us to vary the diffusion timescale appropriately. We then expose the crystal to shear, thereby choosing shear rates which cover timescales that are both higher and lower than the equilibrium diffusion timescales. We investigate the macroscopic and microscopic response of our cluster crystal to shear and find that the yielding scenario of such a system does not rely on the diffusion of the particles - it is rather related to the plastic deformation of the underlying crystalline structure. The local bond order parameters and the measurement of local angles between neighboring clusters confirm the cooperative movement of the clusters close to the yield point. Performing complementary, related simulations for an FCC crystal formed by harshly repulsive particles reveals similarities in the yielding behavior between both systems. Still we find that the diffusion of particles does influence characteristic features in the cluster crystal, such as a less prominent increase of order parameters close to the yield point. Our simulations provide for the first time an insight into the role of the diffusion of defects in the yielding behavior of a defect-rich crystal under shear. These observations will thus be helpful in the development of theories for the plastic deformation of defect-rich crystals.

Intercluster Exchange-Stabilized Novel Complex Colloidal χc Phase

Designing complex cluster crystals with a specific function using simple colloidal building blocks remains a challenge in materials science. Herein, we propose a conceptually new design strategy for constructing complex cluster crystals via hierarchical self-assembly of simple soft Janus colloids. A novel and previously unreported colloidal cluster-χ (χ_c) phase, which resembles the essential structural features of α-manganese but at a larger length scale, is obtained through molecular dynamics simulations. The formation of the χ_c phase undergoes a remarkable two-step self-assembly process, that is, the self-assembly of clusters with specific size dispersity from Janus colloids, followed by the highly ordered organization of these clusters. More importantly, the dynamic exchange of particles between these clusters plays a critical role in stabilizing the χ_c phase. Such a conceptual design framework based on intercluster exchange has the potential to effectively construct novel complex cluster crystals by hierarchical self-assembly of colloidal building blocks.

Unraveling the Role of Architecture in Polymer-Based Glues for Hydrogels

We investigate polymers of different architectures as potential candidates for the development of glues for hydrogels. Using a combination of coarse-grained modeling and molecular dynamics simulations, we systematically characterize the link between experimentally tunable parameters and adhesion energy. We find that, for a broad set of parameters, adhesion is controlled almost exclusively by the total amount of glue at the interface and by the glue-hydrogel affinity. Instead, it is largely independent of changes in polymer architecture and size, a conclusion that shines new light on previously observed experimental trends. Additionally, we show that the scaling behavior of the properties we measure can be explained by modeling the glue as an ensemble of ideal, noninteracting, and linear polymer segments. We expect that the fundamental insights herein provided will aid the design of new polymer-based adhesives for hydrogels.

Direct Correlation Function of a Crystalline Solid

Direct correlation functions (DCFs), linked to the second functional derivative of the free energy with respect to the one-particle density, play a fundamental role in a statistical mechanics description of matter. This holds, in particular, for the ordered phases: DCFs contain information about the local structure including defects and encode the thermodynamic properties of crystalline solids; they open a route to the elastic constants beyond low temperature expansions. Via a demanding numerical approach, we have explicitly calculated for the first time the DCF of a solid: based on the fundamental measure concept, we provide results for the DCF of a hard sphere crystal. We demonstrate that this function differs at coexistence significantly from its liquid counterpart-both in shape as well as in its order of magnitude-because it is dominated by vacancies. We provide evidence that the traditional use of liquid DCFs in functional Taylor expansions of the free energy is conceptually wrong and show that the emergent elastic constants are in good agreement with simulation-based results.

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.

Density Distribution in Soft Matter Crystals and Quasicrystals

The density distribution in solids is often represented as a sum of Gaussian peaks (or similar functions) centered on lattice sites or via a Fourier sum. Here, we argue that representing instead the logarithm of the density distribution via a Fourier sum is better. We show that truncating such a representation after only a few terms can be highly accurate for soft matter crystals. For quasicrystals, this sum does not truncate so easily, nonetheless, representing the density profile in this way is still of great use, enabling us to calculate the phase diagram for a three-dimensional quasicrystal-forming system using an accurate nonlocal density functional theory.

Cluster self-assembly condition for arbitrary interaction potentials

We present a sufficient criterion for the emergence of cluster phases in an ensemble of interacting classical particles with repulsive two-body interactions. Through a zero-temperature analysis in the low density region we determine the relevant characteristics of the interaction potential that make the energy of a two-particle cluster-crystal become smaller than that of a simple triangular lattice in two dimensions. The method leads to a mathematical condition for the emergence of cluster crystals in terms of the sum of Fourier components of a regularized interaction potential, which can be in principle applied to any arbitrary shape of interactions. We apply the formalism to several examples of bounded and unbounded potentials with and without cluster-forming ability. In all cases, the emergence of self-assembled cluster crystals is well captured by the presented analytic criterion and verified with known results from molecular dynamics simulations at vanishingly temperatures. Our work generalises known results for bounded potentials to repulsive potentials of arbitrary shape.

Fluctuation Profiles in Inhomogeneous Fluids

Three one-body profiles that correspond to local fluctuations in energy, in entropy, and in particle number are used to describe the equilibrium properties of inhomogeneous classical many-body systems. Local fluctuations are obtained from thermodynamic differentiation of the density profile or equivalently from average microscopic covariances. The fluctuation profiles follow from functional generators and they satisfy Ornstein-Zernike relations. Computer simulations reveal markedly different fluctuations in confined fluids with Lennard-Jones, hard sphere, and Gaussian core interactions.

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.

Defying the Gibbs Phase Rule: Evidence for an Entropy-Driven Quintuple Point in Colloid-Polymer Mixtures

Using a minimal algebraic model for the thermodynamics of binary rod-polymer mixtures, we provide evidence for a quintuple phase equilibrium; an observation that seems to be at odds with the Gibbs phase rule for two-component systems. Our model is based on equations of state for the relevant liquid crystal phases that are in quantitative agreement with computer simulations. We argue that the appearance of a quintuple equilibrium, involving an isotropic fluid, a nematic and smectic liquid crystal, and two solid phases, can be reconciled with a generalized Gibbs phase rule in which the two intrinsic length scales of the athermal colloid-polymer mixture act as additional field variables.

Freezing of a soft-core fluid in a one-dimensional potential: Predictions based on a pressure-balance equation

Using concepts from classical density functional theory, we investigate the freezing of a two-dimensional system of ultrasoft particles in a one-dimensional external potential, a phenomenon often called laser-induced freezing (LIF). In the first part of the paper, we present numerical results from free minimization of a mean-field density functional for a system of particles interacting via the GEM-4 potential. We show that the system does indeed display a LIF transition, although the interaction potential is markedly different from the cases studied before. We also introduce a suitably defined effective density within the potential wells ρ[over ¯]_{eff} as a control parameter of LIF, rather than the amplitude of the external potential as in the common LIF scenario. In the second part, we suggest a theoretical description of the onset of LIF which is based on the pressure-balance equation relating the pressure tensor and the external potential. Evaluating this equation for the modulated liquid phase at effective density ρ[over ¯]_{eff} and combining it with the (known) stability threshold of the corresponding bulk fluid, we can predict the critical effective density or, equivalently, the potential amplitude related to the onset of LIF. Our approach yields very good results for the model at hand and it is transferable, in principle, to other model systems.

A comparative study between two models of active cluster crystals

We study a system of active particles with soft repulsive interactions that lead to an active cluster-crystal phase in two dimensions. We use two different modelizations of the active force - Active Brownian particles (ABP) and Ornstein-Uhlenbeck particles (AOUP) - and focus on analogies and differences between them. We study the different phases appearing in the system, in particular, the formation of ordered patterns drifting in space without being altered. We develop an effective description which captures some properties of the stable clusters for both ABP and AOUP. As an additional point, we confine such a system in a large channel, in order to study the interplay between the cluster crystal phase and the well-known accumulation near the walls, a phenomenology typical of active particles. For small activities, we find clusters attached to the walls and deformed, while for large values of the active force they collapse in stripes parallel to the walls.

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.

Supersolid Stripe Crystal from Finite-Range Interactions on a Lattice

Strong, long-range interactions present a unique challenge for the theoretical investigation of quantum many-body lattice models, due to the generation of large numbers of competing states at low energy. Here, we investigate a class of extended bosonic Hubbard models with off-site terms interpolating between short and infinite range, thus allowing for an exact numerical solution for all interaction strengths. We predict a novel type of stripe crystal at strong coupling. Most interestingly, for intermediate interaction strengths we demonstrate that the stripes can turn superfluid, thus leading to a self-assembled array of quasi-one-dimensional superfluids. These bosonic superstripes turn into an isotropic supersolid with decreasing the interaction strength. The mechanism for stripe formation is based on cluster self-assembling in the corresponding classical ground state, reminiscent of classical soft-matter models of polymers, different from recently proposed mechanisms for cold gases of alkali or dipolar magnetic atoms.

Phase diagram for two-dimensional layer of soft particles

The phase diagram of a monolayer of soft particles described by the Daoud-Cotton model for star polymers is presented. Ground state calculations and grand canonical Monte Carlo simulations are used to determine the phase behavior as a function of the number of arms on the star and the areal coverage of the soft particles. The phase diagram exhibits rich behavior including reentrant melting and freezing and solid-solid transitions with triangular, stripe, honeycomb and kagome phases. These structures in 2D are analogous to the structures observed in 3D. The evolution of the structure factor with density is qualitatively similar to that measured in experiments for polymer grafted nanocrystals [Chen et al., Macromolecules, 2017, 50, 9636].

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.

Formation of cluster crystals in an ultra-soft potential model on a spherical surface

We investigate the formation of cluster crystals with multiply occupied lattice sites on a spherical surface in systems of ultra-soft particles interacting via repulsive, bounded pair potentials. Not all interactions of this kind lead to clustering: we generalize the criterion devised in C. N. Likos et al., Phys. Rev. E, 2001, 63, 031206 to spherical systems in order to distinguish between cluster-forming systems and fluids which display reentrant melting. We use both DFT and Monte Carlo simulations to characterize the behavior of the system, and obtain semi-quantitative agreement between the two. We find that the number of clusters is determined by the ratio between the size σ of the ultra-soft particles and the radius R of the sphere in such a way that each stable configuration spans a certain interval of σ/R. Furthermore, we study the effect of topological frustration on the system due to the sphere curvature by comparing the properties of disclinations, i.e., clusters with fewer than six neighbors, and non-defective clusters. Disclinations are shown to be less stable, contain fewer particles, and be closer to their neighbors than other lattice points: these properties are explained on the basis of geometric and energetic considerations.

Some general features of mesophase formation in hard-core plus tail potentials

The formation of mesophases in fluids with hard-core plus tail interactions is investigated and compared with the occurrence of cluster crystals in ultra-soft repulsive potentials by using a simple variational expression for the Helmholtz free energy. The purpose of this study is mostly qualitative, i.e., to explain the origin of the different behavior of these systems, and the reason why, in the hard-core case, interactions which are apparently quite different display a common pattern for the phase diagram, featuring spheres, cylinders, lamellae, inverted cylinders, and inverted spheres as the density is increased. In the limit of zero temperature, our approach also yields some simple predictions for the densities at which the transitions between different mesophases are expected to take place, as well as for the size of their clusters at the transitions. We find that these results compare favorably with those obtained in a former study of a model fluid with competing attractive and repulsive interactions by density-functional theory with numerical minimization.

Dynamical density-functional-theory-based modeling of tissue dynamics: Application to tumor growth

We present a theoretical framework based on an extension of dynamical density-functional theory (DDFT) for describing the structure and dynamics of cells in living tissues and tumors. DDFT is a microscopic statistical mechanical theory for the time evolution of the density distribution of interacting many-particle systems. The theory accounts for cell-pair interactions, different cell types, phenotypes, and cell birth and death processes (including cell division), to provide a biophysically consistent description of processes bridging across the scales, including describing the tissue structure down to the level of the individual cells. Analysis of the model is presented for single-species and two-species cases, the latter aimed at describing competition between tumor and healthy cells. In suitable parameter regimes, model results are consistent with biological observations. Of particular note, divergent tumor growth behavior, mirroring metastatic and benign growth characteristics, are shown to be dependent on the cell-pair-interaction parameters.

Hunting mermaids in real space: known knowns, known unknowns and unknown unknowns

We review efforts to realise so-called mermaid (or short-ranged attraction/long ranged repulsion) interactions in 3d real space. The repulsive and attractive contributions to these interactions in charged colloids and colloid-polymer mixtures, may be accurately realised, by comparing particle-resolved studies with colloids to computer simulation. However, when we review work where these interactions have been combined, despite early indications of behaviour consistent with predictions, closer analysis reveals that in the non-aqueous systems used for particle-resolved studies, the idea of summing the attractive and repulsive components leads to wild deviations with experiment. We suggest that the origin lies in the weak ion dissociation in these systems with low dielectric constant solvents. Ultimately this leads even to non-centro-symmetric interactions and a new level of complexity in these systems.

Sliding states of a soft-colloid cluster crystal: Cluster versus single-particle hopping

We study a two-dimensional model for interacting colloidal particles which displays spontaneous clustering. Within this model we investigate the competition between the pinning to a periodic corrugation potential and a sideways constant pulling force which would promote a sliding state. For a few sample particle densities and amplitudes of the periodic corrugation potential we investigate the depinning from the statically pinned to the dynamically sliding regime. This sliding state exhibits the competition between a dynamics where entire clusters are pulled from a minimum to the next and a dynamics where single colloids or smaller groups leave a cluster and move across the corrugation energy barrier to join the next cluster downstream in the force direction. Both kinds of sliding states can occur either coherently across the entire sample or asynchronously: the two regimes result in different average mobilities. Finite temperature tends to destroy separate sliding regimes, generating a smoother dependence of the mobility on the driving force.

Pattern formation and self-assembly driven by competing interactions

Colloidal fluids interacting via effective potentials which are attractive at the short range and repulsive at the long range have long been raising considerable attention because such an instance provides a simple mechanism leading to pattern formation even for isotropic interactions. If the competition between attraction and repulsion is strong enough, the gas-liquid phase transition is suppressed, and replaced by the formation of mesophases, i.e., inhomogeneous phases displaying periodic density modulations whose length, although being larger than the particle size, cannot nevertheless be considered macroscopic. We describe a fully numerical implementation of density-functional theory in three dimensions, tailored to periodic phases. The results for the equilibrium phase diagram of the model are compared with those already obtained in previous investigations for the present system as well as for other systems which form mesophases. The phase diagram which we find shows a strong similarity with that of block copolymer melts, in which self-assembly also results from frustration of a macroscopic phase separation. As the inhomogeneous region is swept by increasing the density from the low-density side, one encounters clusters, bars, lamellae, inverted bars, and inverted clusters. Moreover, a bicontinuous gyroid phase consisting of two intertwined percolating networks is predicted in a narrow domain between the bar and lamellar phases.

Dynamical density functional theory analysis of the laning instability in sheared soft matter

Using dynamical density functional theory (DDFT) methods we investigate the laning instability of a sheared colloidal suspension. The nonequilibrium ordering at the laning transition is driven by nonaffine particle motion arising from interparticle interactions. Starting from a DDFT which incorporates the nonaffine motion, we perform a linear stability analysis that enables identification of the regions of parameter space where lanes form. We illustrate our general approach by applying it to a simple one-component fluid of soft penetrable particles.

Quantum Critical Behavior of One-Dimensional Soft Bosons in the Continuum

We consider a zero-temperature one-dimensional system of bosons interacting via the soft-shoulder potential in the continuum, typical of dressed Rydberg gases. We employ quantum Monte Carlo simulations, which allow for the exact calculation of imaginary-time correlations, and a stochastic analytic continuation method, to extract the dynamical structure factor. At finite densities, in the weakly interacting homogeneous regime, a rotonic spectrum marks the tendency to clustering. With strong interactions, we indeed observe cluster liquid phases emerging, characterized by the spectrum of a composite harmonic chain. Luttinger theory has to be adapted by changing the reference lattice density field. In both the liquid and cluster liquid phases, we find convincing evidence of a secondary mode, which becomes gapless only at the transition. In that region, we also measure the central charge and observe its increase towards c=3/2, as recently evaluated in a related extended Bose-Hubbard model, and we note a fast reduction of the Luttinger parameter. For two-particle clusters, we then interpret such observations in terms of the compresence of a Luttinger liquid and a critical transverse Ising model, related to the instability of the reference lattice density field towards coalescence of sites, typical of potentials which are flat at short distances. Even in the absence of a true lattice, we are able to evaluate the spatial correlation function of a suitable pseudospin operator, which manifests ferromagnetic order in the cluster liquid phase, exponential decay in the liquid phase, and algebraic order at criticality.