We recently reported the experimental realization of a chiral artificial microswimmer exhibiting helical self-propulsion [T. Yamamoto and M. Sano, Soft Matter 13, 3328 (2017)]. In the experiment, cholesteric liquid crystal (CLC) droplets dispersed in surfactant solutions swam spontaneously, driven by the Marangoni flow, in helical paths whose handedness is determined by the chirality of the component molecules of CLC. To study the mechanism of the emergence of the helical self-propelled motion, we propose a phenomenological model of the self-propelled helical motion of the CLC droplets. Our model is constructed by symmetry argument in chiral systems, and it describes the dynamics of CLC droplets with coupled time-evolution equations in terms of a velocity, an angular velocity, and a tensor variable representing the symmetry of the helical director field of the droplet. We found that helical motions as well as other chiral motions appear in our model. By investigating bifurcation behaviors between each chiral motion, we found that the chiral coupling terms between the velocity and the angular velocity, the structural anisotropy of the CLC droplet, and the nonlinearity of model equations play a crucial role in the emergence of the helical motion of the CLC droplet.

The normal modes and relaxation rates of weak colloidal gels are investigated in calculations using different models of the hydrodynamic interactions between suspended particles. The relaxation spectrum is computed for freely draining, Rotne-Prager-Yamakawa, and accelerated Stokesian dynamics approximations of the hydrodynamic mobility in a normal mode analysis of a harmonic network representing several colloidal gels. We find that the density of states and spatial structure of the normal modes are fundamentally altered by long-ranged hydrodynamic coupling among the particles. Short-ranged coupling due to hydrodynamic lubrication affects only the relaxation rates of short-wavelength modes. Hydrodynamic models accounting for long-ranged coupling exhibit a microscopic relaxation rate for each normal mode, $\lambda$ that scales as $l^{-2}$, where $l$ is the spatial correlation length of the normal mode. For the freely draining approximation, which neglects long-ranged coupling, the microscopic relaxation rate scales as $l^{-\gamma }$, where $\gamma$ varies between three and two with increasing particle volume fraction. A simple phenomenological model of the internal elastic response to normal mode fluctuations is developed, which shows that long-ranged hydrodynamic interactions play a central role in the viscoelasticity of the gel network. Dynamic simulations of hard spheres that gel in response to short-ranged depletion attractions are used to test the applicability of the density of states predictions. For particle concentrations up to 30% by volume, the power law decay of the relaxation modulus in simulations accounting for long-ranged hydrodynamic interactions agrees with predictions generated by the density of states of the corresponding harmonic networks as well as experimental measurements. For higher volume fractions, excluded volume interactions dominate the stress response, and the prediction from the harmonic network density of states fails. Analogous to the Zimm model in polymer physics, our results indicate that long-ranged hydrodynamic interactions play a crucial role in determining the microscopic dynamics and macroscopic properties of weak colloidal gels.

The impact of an intruder on granular matter leads to the formation of mesoscopic force networks, which were seen particularly clearly in the recent experiments carried out with photoelastic particles [Clark et al., Phys. Rev. Lett. 114, 144502 (2015)]. These force networks are characterized by complex structure and evolve on fast time scales. While it is known that total photoelastic activity in the granular system is correlated with the acceleration of the intruder, it is not known how the structure of the force network evolves during impact, and if there are dominant features in the networks that can be used to describe the intruder's dynamics. Here, we use topological tools, in particular persistent homology, to describe these features. Persistent homology allows quantification of both structure and time evolution of the resulting force networks. We find that there is a clear correlation of the intruder's dynamics and some of the topological measures implemented. This finding allows us to discuss which properties of the force networks are most important when attempting to describe the intruder's dynamics. In particular, we find that the presence of loops in the force network, quantified by persistent homology, is strongly correlated to the deceleration of the intruder. In some cases, particularly for the impact on soft particles, the measures derived from the persistence analysis describe the deceleration of the intruder even better than the total photoelastic activity. We are also able to define an upper bound on the relevant time scale over which the force networks evolve.

We study the impact of quantum phase transitions (QPTs) on the distribution of exceptional points (EPs) of the Hamiltonian in the complex-extended parameter domain. Analyzing first- and second-order QPTs in the Lipkin-Meshkov-Glick model we find an exponentially and polynomially close approach of EPs to the respective critical point with increasing size of the system. If the critical Hamiltonian is subject to random perturbations of various kinds, the averaged distribution of EPs close to the critical point still carries decisive information on the QPT type. We therefore claim that properties of the EP distribution represent a parametrization-independent signature of criticality in quantum systems.

A scalar Langevin-type process $X\left(t\right)$ that is driven by Ornstein-Uhlenbeck noise $\eta \left(t\right)$ is non-Markovian. However, the joint dynamics of $X$ and $\eta$ is described by a Markov process in two dimensions. But even though there exists a variety of techniques for the analysis of Markov processes, it is still a challenge to estimate the process parameters solely based on a given time series of $X$. Such a partially observed 2D process could, e.g., be analyzed in a Bayesian framework using Markov chain Monte Carlo methods. Alternatively, an embedding strategy can be applied, where first the joint dynamics of $X$ and its temporal derivative $\dot{X}$ is analyzed. Subsequently, the results can be used to determine the process parameters of $X$ and $\eta$. In this paper, we propose a more direct approach that is purely based on the moments of the increments of $X$, which can be estimated for different time-increments $\tau$ from a given time series. From a stochastic Taylor expansion of $X$, analytic expressions for these moments can be derived, which can be used to estimate the process parameters by a regression strategy.

In this work, we show that a universal quantum work relation for a quantum system driven arbitrarily far from equilibrium extends to a parity-time- ($\mathcal{PT}$-) symmetric quantum system with unbroken $\mathcal{PT}$ symmetry, which is a consequence of microscopic reversibility. The quantum Jarzynski equality, linear response theory, and Onsager reciprocal relations for the $\mathcal{PT}$-symmetric quantum system are recovered as special cases of the universal quantum work relation in a $\mathcal{PT}$-symmetric quantum system. In the regime of broken $\mathcal{PT}$ symmetry, the universal quantum work relation does not hold because the norm is not preserved during the dynamics.

The nature of social dilemmas is studied in $n$-strategy evolutionary potential games on a square lattice with nearest-neighbor interactions and the logit rule. For symmetric games with symmetric payoff matrices there are no dilemmas because of the coincidence of individual and common interests. The dilemmas are caused by the antisymmetric parts of the self- and cross-dependent payoff components if it modifies the preferred Nash equilibrium. The contentment of players and the emergence of dilemmas in the preferred Nash equilibria are illustrated on some two-dimensional cross sections of the parameter space.

In population dynamics, the concept of structural stability has been used to quantify the tolerance of a system to environmental perturbations. Yet, measuring the structural stability of nonlinear dynamical systems remains a challenging task. Focusing on the classic Lotka-Volterra dynamics, because of the linearity of the functional response, it has been possible to measure the conditions compatible with a structurally stable system. However, the functional response of biological communities is not always well approximated by deterministic linear functions. Thus, it is unclear the extent to which this linear approach can be generalized to other population dynamics models. Here, we show that the same approach used to investigate the classic Lotka-Volterra dynamics, which is called the structural approach, can be applied to a much larger class of nonlinear models. This class covers a large number of nonlinear functional responses that have been intensively investigated both theoretically and experimentally. We also investigate the applicability of the structural approach to stochastic dynamical systems and we provide a measure of structural stability for finite populations. Overall, we show that the structural approach can provide reliable and tractable information about the qualitative behavior of many nonlinear dynamical systems.

Dynamical bending, buckling, and polymorphic transformations of the flagellum are known to affect bacterial motility, but run-reverse-flick motility of monotrichous bacteria also involves the even more flexible hook connecting the flagellum to its rotary motor. Although flick initiation has been hypothesized to involve either static Euler buckling or dynamic bending of the hook, the precise mechanism of flick initiation remains unknown. Here, we find that flicks initiate via a dynamic instability requiring flexibility in both the hook and flagellum. We obtain accurate estimates of forces and torques on the hook that suggest that flicks occur for stresses below the (static) Euler buckling criterion, then provide a mechanistic model for flick initiation that requires combined bending of the hook and flagellum. We calculate the triggering torque-stiffness ratio and find that our predicted onset of dynamic instability corresponds well with experimental observations.

We investigate by means of molecular dynamics simulation a coarse-grained polymer glass model focusing on (quasistatic and dynamical) shear-stress fluctuations as a function of temperature $T$ and sampling time $\mathrm{\Delta }t$. The linear response is characterized using (ensemble-averaged) expectation values of the contributions (time averaged for each shear plane) to the stress-fluctuation relation ${\mu }_{\mathrm{sf}}$ for the shear modulus and the shear-stress relaxation modulus $G\left(t\right)$. Using 100 independent configurations, we pay attention to the respective standard deviations. While the ensemble-averaged modulus ${\mu }_{\mathrm{sf}}\left(T\right)$ decreases continuously with increasing $T$ for all $\mathrm{\Delta }t$ sampled, its standard deviation $\delta {\mu }_{\mathrm{sf}}\left(T\right)$ is nonmonotonic with a striking peak at the glass transition. The question of whether the shear modulus is continuous or has a jump singularity at the glass transition is thus ill posed. Confirming the effective time-translational invariance of our systems, the $\mathrm{\Delta }t$ dependence of ${\mu }_{\mathrm{sf}}$ and related quantities can be understood using a weighted integral over $G\left(t\right)$.

Self-propelled Brownian particles show rich out-of-equilibrium physics, for instance, the motility-induced phase separation (MIPS). While decades of studying the structure of liquids have established a deep understanding of passive systems, not much is known about correlations in active suspensions. In this work we derive an approximate analytic theory for three-body correlations and forces in systems of active Brownian disks starting from the many-body Smoluchowski equation. We use our theory to predict the conditional forces that act on a tagged particle and their dependence on the propulsion speed of self-propelled disks. We identify preferred directions of these forces in relation to the direction of propulsion and the positions of the surrounding particles. We further relate our theory to the effective swimming speed of the active disks, which is relevant for the physics of MIPS. To test and validate our theory, we additionally run particle-resolved computer simulations, for which we explicitly calculate the three-body forces. In this context, we discuss the modeling of active Brownian swimmers with nearly hard interaction potentials. We find very good agreement between our simulations and numerical solutions of our theory, especially for the nonequilibrium pair-distribution function. For our analytical results, we carefully discuss their range of validity in the context of the different levels of approximation we applied. This discussion allows us to study the individual contribution of particles to three-body forces and to the emerging structure. Thus, our work sheds light on the collective behavior, provides the basis for further studies of correlations in active suspensions, and makes a step towards an emerging liquid state theory.

We probe the effects of particle shape on the global and local behavior of a two-dimensional granular pillar, acting as a proxy for a disordered solid, under uniaxial compression. This geometry allows for direct measurement of global material response, as well as tracking of all individual particle trajectories. In general, drawing connections between local structure and local dynamics can be challenging in amorphous materials due to lower precision of atomic positions, so this study aims to elucidate such connections. We vary local interactions by using three different particle shapes: discrete circular grains (monomers), pairs of grains bonded together (dimers), and groups of three bonded in a triangle (trimers). We find that dimers substantially strengthen the pillar and the degree of this effect is determined by orientational order in the initial condition. In addition, while the three particle shapes form void regions at distinct rates, we find that anisotropies in the local amorphous structure remain robust through the definition of a metric that quantifies packing anisotropy. Finally, we highlight connections between local deformation rates and local structure.

We simulate numerically the compression-driven jamming of athermal, frictionless, soft-core spherocylinders in two dimensions, for a range of particle aspect ratios $\alpha$. We find the critical packing fraction ${\varphi }_J\left(\alpha \right)$ for the jamming transition and the average number of contacts per particle $z_J\left(\alpha \right)$ at jamming. We find that both are nonmonotonic, with a peak at $\alpha \approx 1$. We find that configurations at the compression-driven jamming point are always hypostatic for all $\alpha$, with $z_J<z_{\mathrm{iso}}=2d_f=6$ the isostatic value. We show that, for moderately elongated spherocylinders, there is no orientational ordering upon athermal compression through jamming. We analyze in detail the eigenmodes of the dynamical matrix close to the jamming point for a few different values of the aspect ratio, from nearly circular to moderately elongated. We find that there are low frequency bands containing $N(z_{\mathrm{iso}}-z_J)/2$ modes, such that the frequencies of these modes vanish as $\varphi \to {\varphi }_J$. We consider the extended versus localized nature of these low frequency modes, and the extent to which they involve translational or rotational motion, and find many low frequency sliding modes where particles can move with little rotation. We highlight the importance of treating side-to-side contacts, along flat sides of the spherocylinder, properly for the correct determination of $z_J$. We note the singular nature of taking the $\alpha \to 0$ limit. We discuss the similarities and differences with previous work on jammed ellipses and ellipsoids, to illustrate the effects that different particle shapes have on configurations at jamming.

Thermal convection is numerically computed in three-dimensional (3D) fluid saturated isotropically fractured porous media. Fractures are randomly inserted as two-dimensional (2D) convex polygons. Flow is governed by Darcy's 2D and 3D laws in the fractures and in the porous medium, respectively; exchanges take place between these two structures. Results for unfractured porous media are in agreement with known theoretical predictions. The influence of parameters such as the fracture aperture (or fracture transmissivity) and the fracture density on the heat released by the whole system is studied for Rayleigh numbers up to 150 in cubic boxes with closed-top conditions. Then, fractured media are compared to homogeneous porous media with the same macroscopic properties. Three major results could be derived from this study. The behavior of the system, in terms of heat release, is determined as a function of fracture density and fracture transmissivity. First, the increase in the output flux with fracture density is linear over the range of fracture density tested. Second, the increase in output flux as a function of fracture transmissivity shows the importance of percolation. Third, results show that the effective approach is not always valid, and that the mismatch between the full calculations and the effective medium approach depends on the fracture density in a crucial way.

In this work, an efficient gas kinetic scheme is presented for simulation of two-dimensional incompressible thermal flows. In the scheme, the macroscopic governing equations for mass, momentum, and energy conservation are discretized by the finite volume method and the numerical fluxes at the cell interface are reconstructed by the local solution of the Boltzmann equation. To compute these fluxes, two distribution functions are involved. One is the circular function, which is used to calculate the numerical fluxes of mass and momentum equations. Due to the incompressible limit, the circle at the cell interface can be approximately considered to be symmetric so that the expressions for the conservative variables and numerical fluxes at the cell interface can be given explicitly and concisely. Another one is the D2Q4 model, which is utilized to compute the numerical flux of the energy equation. By following the process for derivation of numerical fluxes of mass and momentum equations, the numerical flux of the energy equation can also be given explicitly. The accuracy, efficiency, and stability of the present scheme are validated by simulating several thermal flow problems. Numerical results showed that the present scheme can provide accurate numerical results for incompressible thermal flows at a wide range of Rayleigh numbers with less computational cost than that needed by the thermal lattice Boltzmann flux solver (TLBFS), which has been proven to be more efficient than the thermal lattice Boltzmann method (TLBM).

Disordered multihyperuniform many-particle systems are exotic amorphous states that allow exquisite color sensing capabilities due to their anomalous suppression of density fluctuations for distinct subsets of particles, as recently evidenced in photoreceptor mosaics in avian retina. Motivated by this biological finding, we present a statistical-mechanical model that rigorously achieves disordered multihyperuniform many-body systems by tuning interactions in binary mixtures of nonadditive hard-disk plasmas. We demonstrate that multihyperuniformity competes with phase separation and stabilizes a clustered phase. Our work provides a systematic means to generate disordered multihyperuniform solids, and hence lays the groundwork to explore their potentially unique photonic, phononic, electronic, and transport properties.

We report measurements of the equations of state of a fully relaxed magnetohydrodynamic (MHD) laboratory plasma. Parcels of magnetized plasma, called Taylor states, are formed in a coaxial magnetized plasma gun, and are allowed to relax and drift into a closed flux conserving volume. Density, ion temperature, and magnetic field are measured as a function of time as the Taylor states compress and heat. The theoretically predicted MHD and double adiabatic equations of state are compared to experimental measurements. We find that the MHD equation of state is inconsistent with our data.

We present extensive numerical studies of the crossover from three-dimensional to two-dimensional systems in the nonequilibrium zero-temperature random-field Ising model with metastable dynamics. Bivariate finite-size scaling hypotheses are presented for systems with sizes $L\times L\times l$ which explain the size-driven critical crossover from two dimensions ($l=\text{const}$, $L\to \infty$) to three dimensions ($l\propto L\to \infty$). A model of effective critical disorder $R_c^{\mathrm{eff}}(l,L)$ with a unique fitting parameter and no free parameters in the $R_c^{\mathrm{eff}}(l,L\to \infty )$ limit is proposed, together with expressions for the scaling of avalanche distributions bringing important implications for related experimental data analysis, especially in the case of thin three-dimensional systems.

A well-known stochastic model for intermittent fluctuations in physical systems is investigated. The model is given by a superposition of uncorrelated exponential pulses, and the degree of pulse overlap is interpreted as an intermittency parameter. Expressions for excess time statistics, that is, the rate of level crossings above a given threshold and the average time spent above the threshold, are derived from the joint distribution of the process and its derivative. Limits of both high and low intermittency are investigated and compared to previously known results. In the case of a strongly intermittent process, the distribution of times spent above threshold is obtained analytically. This expression is verified numerically, and the distribution of times above threshold is explored for other intermittency regimes. The numerical simulations compare favorably to known results for the distribution of times above the mean threshold for an Ornstein-Uhlenbeck process. This contribution generalizes the excess time statistics for the stochastic model, which find applications in a wide diversity of natural and technological systems.

We consider the contact process near an extended surface defect, where the local control parameter deviates from the bulk one by an amount of $\lambda \left(l\right)-\lambda \left(\infty \right)=A{l}^{-s}$, with $l$ being the distance from the surface. We concentrate on the marginal situation $s=1/{\nu }_{\perp }$, where ${\nu }_{\perp }$ is the critical exponent of the spatial correlation length, and study the local critical properties of the one-dimensional model by Monte Carlo simulations. The system exhibits a rich surface critical behavior. For weaker local activation rates $A<A_c$, the phase transition is continuous, having an order-parameter critical exponent, which varies continuously with $A$. For stronger local activation rates $A>A_c$, the phase transition is of mixed order: the surface order parameter is discontinuous; at the same time the temporal correlation length diverges algebraically as the critical point is approached, but with different exponents on the two sides of the transition. The mixed-order transition regime is analogous to that observed recently at a multiple junction and can be explained by the same type of scaling theory.