Recovery of ground-state degeneracy in two-dimensional square ice is a significant challenge in the field of geometric frustration with far-reaching fundamental implications, such as realization of vertex models and understanding the effect of dimensionality reduction. We combine experiments, theory, and numerical simulations to demonstrate that sheared square colloidal ice partially recovers the ground-state degeneracy for a wide range of field strengths and lattice shear angles. Our method could inspire engineering a novel class of frustrated microstructures and nanostructures based on sheared magnetic lattices in a wide range of soft- and condensed-matter systems.

We perform friction experiments between a compliant gel and a rigid cylinder at sliding velocities comparable to the Rayleigh wave or secondary wave velocity of the gel. We find that, when the sliding velocity exceeds the wave velocities, the contact state transitions from Hertzian like to flat punch like, resulting in the breakdown of the lubricating oil film and the abrupt increase in the friction coefficient. We succeed in deriving theoretical solutions for the contact pressure distributions and the deformation profiles in the presence of friction, which are consistent with our experimental observations.

In situ measurements of the free energy difference between the open and closed states of ion channels are challenging due to hysteresis effects and inactivation. Exploiting recent developments in statistical physics, we present a general formalism to extract the free energy difference $\mathrm{\Delta }F$ between the closed and open states of mechanosensitive ion channels from nonequilibrium work distributions associated with the opening and closing of the channels (gating) in response to ramp stimulation protocols recorded in native patches. We show that the work distributions obtained from the gating of MscS channels in E. coli membrane satisfy the strong symmetry relation predicted by the Crooks fluctuation theorem. Our approach enables the determination of $\mathrm{\Delta }F$ using patch-clamp experiments, which are often inherently restricted to the nonequilibrium regime.

Computer simulations of the fluid-to-solid phase transition in the hard sphere system were instrumental for our understanding of crystallization processes. But while colloid experiments and theory have been predicting the stability of several binary hard sphere crystals for many years, simulations were not successful to confirm this phenomenon. Here, we report the growth of binary hard sphere crystals isostructural to Laves phases, ${\mathrm{AlB}}_2$, and ${\mathrm{NaZn}}_{13}$ in simulation directly from the fluid. We analyze particle kinetics during Laves phase growth using event-driven molecular dynamics simulations with and without swap moves that speed up diffusion. The crystallization process transitions from nucleation and growth to spinodal decomposition already deep within the fluid-solid coexistence regime. Finally, we present packing fraction–size ratio state diagrams in the vicinity of the stability regions of three binary crystals.

Transport properties of dense fluids are fundamentally challenging, because the powerful approaches of equilibrium statistical physics cannot be applied. Polar fluids compound this problem, because the long-range interactions preclude the use of a simple effective diameter approach based solely on hard spheres. Here, we develop a kinetic theory for dipolar hard-sphere fluids that is valid up to high density. We derive a mathematical approximation for the radial distribution function at contact directly from the equation of state, and use it to obtain the shear viscosity. We also perform molecular-dynamics simulations of this system and extract the shear viscosity numerically. The theoretical results compare favorably to the simulations.

The higher-order interactions of complex systems, such as the brain, are captured by their simplicial complex structure and have a significant effect on dynamics. However, the existing dynamical models defined on simplicial complexes make the strong assumption that the dynamics resides exclusively on the nodes. Here we formulate the higher-order Kuramoto model which describes the interactions between oscillators placed not only on nodes but also on links, triangles, and so on. We show that higher-order Kuramoto dynamics can lead to an explosive synchronization transition by using an adaptive coupling dependent on the solenoidal and the irrotational component of the dynamics.

Two-dimensional systems may admit a hexatic phase and hexatic-liquid transitions of different natures. The determination of their phase diagrams proved challenging, and indeed, those of hard disks, hard regular polygons, and inverse power-law potentials have only recently been clarified. In this context, the role of attractive forces is currently speculative, despite their prevalence at both the molecular and colloidal scale. Here, we demonstrate, via numerical simulations, that attraction promotes a discontinuous melting scenario with no hexatic phase. At high-temperature, Lennard-Jones particles and attractive polygons follow the shape-dominated melting scenario observed in hard disks and hard polygons, respectively. Conversely, all systems melt via a first-order transition with no hexatic phase at low temperature, where attractive forces dominate. The intermediate temperature melting scenario is shape dependent. Our results suggest that, in colloidal experiments, the tunability of the strength of the attractive forces allows for the observation of different melting scenarios in the same system.

The evolution of high-dimensional phenotypes is investigated using a statistical physics model consisting of interacting spins, in which phenotypes, genotypes, and environments are represented by spin configurations, interaction matrices, and external fields, respectively. We found that phenotypic changes upon diverse environmental change and genetic variation are highly correlated across all spins, consistent with recent experimental observations of biological systems. The dimension reduction in phenotypic changes is shown to be a result of the evolution of the robustness to thermal noise, achieved at the replica symmetric phase.

The link between local structure and dynamical slowdown in glassy fluids has been the focus of intense debate for the better part of a century. Nonetheless, a simple method to predict the dynamical behavior of a fluid purely from its local structural features is still missing. Here, we demonstrate that the diffusivity of perhaps the most fundamental family of glass formers—hard sphere mixtures—can be accurately predicted based on just the packing fraction and a simple order parameter measuring the tetrahedrality of the local structure. Essentially, we show that the number of tetrahedral clusters in a hard sphere mixture is directly linked to its global diffusivity. Moreover, the same order parameter is capable of locally pinpointing particles in the system with high and low mobility. We attribute the power of the local tetrahedrality for predicting local and global dynamics to the high stability of tetrahedral clusters, the most fundamental building and densest-packing building blocks for a disordered fluid.

We investigate the aggregation and phase separation of thin, living T. tubifex worms that behave as active polymers. Randomly dispersed active worms spontaneously aggregate to form compact, highly entangled blobs, a process similar to polymer phase separation, and for which we observe power-law growth kinetics. We find that the phase separation of active polymerlike worms does not occur through Ostwald ripening, but through active motion and coalescence of the phase domains. Interestingly, the growth mechanism differs from conventional growth by droplet coalescence: the diffusion constant characterizing the random motion of a worm blob is independent of its size, a phenomenon that can be explained from the fact that the active random motion arises from the worms at the surface of the blob. This leads to a fundamentally different phase-separation mechanism that may be unique to active polymers.

Network flows often exhibit a hierarchical treelike structure that can be attributed to the minimization of dissipation. The common feature of such systems is a single source and multiple sinks (or vice versa). In contrast, here we study networks with only a single source and sink. These systems can arise from secondary purposes of the networks, such as blood sugar regulation through insulin production. Minimization of dissipation in these systems leads to vascular shunting, a single vessel connecting the inlet and outlet. We show instead how optimizing the transport time yields network topologies that match those observed in the insulin-producing pancreatic islets. These are patterns of periphery-to-center and center-to-periphery flows. The obtained flow networks are broadly independent of how the flow velocity depends on the flow flux, but continuous and discontinuous phase transitions appear at extreme flux dependencies. Lastly, we show how constraints on flows can lead to buckling of the branches of the network, a feature that is also observed in pancreatic islets.

A choreographic time crystal is a dynamic lattice structure in which the points comprising the lattice move in a coordinated fashion. These structures were initially proposed for understanding the motion of synchronized satellite swarms. Using simulations, we examine colloids interacting with a choreographic crystal consisting of traps that could be created optically. As a function of the trap strength, speed, and colloidal filling fraction, we identify a series of phases including states where the colloids organize into a dynamic chiral loop lattice as well as a frustrated induced liquid state and a choreographic lattice state. We show that transitions between these states can be understood in terms of vertex frustration effects that occur during a certain portion of the choreographic cycle. Our results can be generalized to a broader class of systems of particles coupled to choreographic structures, such as vortices, ions, cold atoms, and soft matter systems.

We introduce a new model to study the effect of surface roughness on the jamming transition. By performing numerical simulations, we show that for a smooth surface, the jamming transition density and the contact number at the transition point both increase upon increasing asphericity, as for ellipsoids and spherocylinders. Conversely, for a rough surface, both quantities decrease, in quantitative agreement with the behavior of frictional particles. Furthermore, in the limit corresponding to the Coulomb friction law, the model satisfies a generalized isostaticity criterion proposed in previous studies. We introduce a counting argument that justifies this criterion and interprets it geometrically. Finally, we propose a simple theory to predict the contact number at finite friction from the knowledge of the force distribution in the infinite friction limit.

Telechelic polymers contain two chain ends that are able to promote self-assembly into “flowerlike” or interconnected micellar structures. Here, we investigate the molecular exchange kinetics of such micelles using time-resolved small-angle neutron scattering. We show that the activation energies of monofunctional and telechelic chain exchange are identical. This demonstrates that the two chain ends are not simultaneously released in a single event. Instead, the results show that, contrary to regular micelles, the kinetics occurs in a multistep process involving a collision-induced single-molecule exchange mechanism where the exchange rate is directly proportional to the polymer concentration. We show that this novel mechanism can be quantitatively explained by a simple kinetic model.