Condensed Matter

Classical dynamical density functional theory (DDFT) is one of the cornerstones of modern statistical mechanics. It is an extension of the highly successful method of classical density functional theory (DFT) to nonequilibrium systems. Originally developed for the treatment of simple and complex fluids, DDFT is now applied in fields as diverse as hydrodynamics, materials science, chemistry, biology, and plasma physics. In this review, we give a broad overview over classical DDFT. We explain its theoretical foundations and the ways in which it can be derived. The relations between the different forms of deterministic and stochastic DDFT as well as between DDFT and related theories, such as quantum-mechanical time-dependent DFT, mode coupling theory, and phase field crystal models, are clarified. Moreover, we discuss the wide spectrum of extensions of DDFT, which covers methods with additional order parameters (like extended DDFT), exact approaches (like power functional theory), and systems with more complex dynamics (like active matter). Finally, the large variety of applications, ranging from fluid mechanics and polymer physics to solidification, pattern formation, biophysics, and electrochemistry, is presented.

The colonisation of a soft passive material by motile cells such as bacteria is common in biology. The resulting colonies of the invading cells are often observed to exhibit intricate patterns whose morphology and dynamics can depend on a number of factors, particularly the mechanical properties of the substrate and the motility of the individual cells. We use simulations of a minimal 2D model of self-propelled rods moving through with a passive compliant medium consisting of particles that offer elastic resistance before being plastically displaced from their equilibrium positions. It is observed that the motility-induced clustering of active (self-propelled) particles is crucial for understanding the morphodynamics of colonisation. Clustering enables motile colonies to spread faster than they would have as isolated particles. The colonisation rate depends non-monotonically on substrate stiffness with a distinct maximum at a non-zero value of substrate stiffness. This is observed to be due to a change in the morphology of clusters. Furrow networks created by the active particles have a fractal-like structure whose dimension varies systematically with substrate stiffness but is less sensitive to particle activity. The power-law growth exponent of the furrowed area is smaller than unity, suggesting that, to sustain such extensive furrow networks, colonies must regulate their overall growth rate.

We investigate coarse-grained models of suspended self-thermophoretic microswimmers. Upon heating, the Janus spheres, with hemispheres made of different materials, induce a heterogeneous local solvent temperature that causes the self-phoretic particle propulsion. Starting from atomistic molecular dynamics simulations, we verify the coarse-grained description of the fluid in terms of a local molecular temperature field, and its role for the particle's thermophoretic self-propulsion and hot Brownian motion. The latter is governed by effective nonequilibrium temperatures, which are measured from simulations by confining the particle position and orientation. They are theoretically shown to remain relevant for any further spatial coarse-graining towards a hydrodynamic description of the entire suspension as a homogeneous complex fluid.

Coarse-grained (CG) conformational surface hopping (SH) adapts the concept of multisurface dynamics, initially developed to describe electronic transitions in chemical reactions, to accurately describe classical molecular dynamics at a reduced level. The SH scheme couples distinct conformational basins (states), each described by its own force field (surface), resulting in a significant improvement of the approximation to the many-body potential of mean force [Phys. Rev. Lett. 121, 256002 (2018)]. The present study first describes CG SH in more detail, through both a toy model and a three-bead model of hexane. We further extend the methodology to non-bonded interactions and report its impact on liquid properties. Finally, we investigate the transferability of the surfaces to distinct systems and thermodynamic state points, through a simple tuning of the state probabilities. In particular, applications to variations in temperature and chemical composition show good agreement with reference atomistic calculations, introducing a promising "weak-transferability regime," where CG force fields can be shared across thermodynamic and chemical neighborhoods.

The conformation of a long linear polymer dissolved in fluid and exposed to an extensional flow is well-known to exhibit a "coil-stretch" transition, which for sufficiently long chains can lead to bistability. The present work reports computations indicating that an analogous "compact-stretched" transition arises in the dynamics of a thin elastic sheet. Sheets of nominally circular, square or rectangular shape are simulated in planar and biaxial flows using a finite element method for the sheet conformations and a regularized Stokeslet method for the fluid flow. If a neo-Hookean constitutive model is used for the sheet elasticity, the sheets will stretch without bound once a critical extension rate, as characterized nondimensionally by a capillary number, is exceeded. Nonlinear elasticity, represented with the Yeoh model, arrests the stretching, leading to a highly-stretched steady state once the critical capillary number is exceeded. For all shapes and in both planar and biaxial extension, a parameter regime exists in which both weakly stretched (compact) and strongly stretched states can be found, depending on initial conditions. I.e. this parameter regime displays bistability. As in the long-chain polymer case, the bistable behavior arises from the hydrodynamic interaction between distant elements of the sheet, and vanishes if these interactions are artificially screened by use of a Brinkman model for the fluid motion. While the sheets can transiently display wrinkled shapes, all final shapes in planar and biaxial extension are planar.

Preparation of nanoemulsions of triglyceride oils in water usually requires high mechanical energy and sophisticated equipment. Recently, we showed that alpha-to-beta (viz. gel-to-crystal) phase transition, observed with most lipid substances (triglycerides, diglycerides, phospholipids, alkanes, etc.), may cause spontaneous disintegration of micro-particles of these lipids, dispersed in aqueous solutions of appropriate surfactants, into nanometer particles/drops using a simple cooling/heating cycle of the lipid dispersion (Cholakova et al. ACS Nano 14 (2020) 8594). In the current study we show that this "cold-burst process" is observed also with natural oils of high practical interest, incl. coconut oil, palm kernel oil and cocoa butter. Mean drop diameters of ca. 50 to 100 nm were achieved with some of the studied oils. From the results of dedicated model experiments we conclude that intensive nano-fragmentation is observed when the following requirements are met: (1) The three phase contact angle at the air-water-solid lipid interface is below ca. 30 degrees; (2) The equilibrium surface tension of the surfactant solution is below ca. 30 mN/m and the dynamic surface tension decreases rapidly. (3) The surfactant solution contains non-spherical surfactant micelles. e.g. ellipsoidal micelles or bigger supramolecular aggregates; (4) The three phase contact angle measured at the contact line (frozen oil-melted oil-surfactant solution) is also relatively low. The mechanism(s) of the particle bursting process is revealed and, on this basis, the role of all these factors is clarified and discussed. We explain all main effects observed experimentally and define guiding principles for optimization of the cold-burst process in various, practically relevant lipid-surfactant systems.

The phenomenon of solidification of a substance from its liquid phase is of the greatest practical and theoretical importance, and atomistic simulations can provide precious information towards its understanding and control. Unfortunately, the time scale for crystallization is much larger than what can be explored in standard simulations. Enhanced sampling methods can overcome this time scale hurdle. Here we employ the on-the-fly probability enhanced sampling method that is a recent evolution of metadynamics. This method, like many others, relies on the definition of appropriate collective variables able to capture the slow degrees of freedom. To this effect we introduce collective coordinates of general applicability to crystallization simulations. They are based on the peaks of the three-dimensional structure factor that are combined non-linearly via the Deep Linear Discriminant Analysis machine learning method. We apply the method to the study of crystallization of a multicomponent system, Sodium Chloride and a molecular system, Carbon Dioxide.

García-Aguilar et al. [Phys. Rev. Lett 126, 038001 (2021)] have shown that the deformations of "shape-shifting droplets" are consistent with an elastic model, that, unlike previous models, includes the intrinsic curvature of the frozen surfactant layer. In this Comment, we show that the interplay between surface tension and intrinsic curvature in their model is in fact mathematically equivalent to a physically very different phase-transition mechanism of the same process that we developed previously [Phys. Rev. Lett. 118, 088001 (2017); Phys. Rev. Res. 1, 023017 (2019)]. The mathematical models cannot therefore distinguish between the two mechanisms, and hence it is not possible to claim that one mechanism underlies all observed shape-shifting phenomena without a much more detailed comparison of experiment and theory.

Huerta et al. [Phys. Rev. Research 2, 033351 (2020)] report a power-law decay of positional order in numerical simulations of hard disks confined within hard parallel walls, which they interpret as a Kosterlitz-Thouless-type caging-uncaging transition. The proposed existence of such a transition in a quasi-one-dimensional (q1D) system, however, contradicts long-held physical expectations. To clarify if the proposed ordering persists in the thermodynamic limit, we introduce an exact transfer matrix approach to expeditiously generate equilibrium configurations for systems of arbitrary size. The power-law decay of positional order is found to extend only over finite distances. We conclude that the numerical simulation results reported are associated with a crossover, and not a proper thermodynamic phase transition.

We study the buckling of a one fiber composite whose matrix stiffness is slightly dependent on the compressive force. We show that the equilibrium curves of the system exhibit a limit load when the induced stiffness parameter gets bigger than a threshold. This limit load increases when increasing the stiffness parameter and is related to a possible localized path in the post-buckling domain. Such a change in the maximum load may be very desirable from a structural stand point.