Liquid sulfur consists mainly of eight-membered rings and hence can be regarded as a model of a molecular liquid. A liquid, which is built from different molecular structures, will demonstrate a wide range in relaxation processes and excitation modes. Three inelastic neutron scattering experiments have been performed to study the collective dynamics of liquid sulfur over three decades in frequencies. A wide range of wave vectors was studied to reveal the response of density fluctuations over different lengthscales. A viscoelastic model with a two-times memory function was applied to the data. The analysis revealed a slow relaxation mode, an acoustic-type excitation, and a high-frequency mode, which resembles an optic-type excitation. The wave-vector dependence of the slow relaxation mode width exhibits the signs of a de Gennes narrowing around the wave vector where the structure factor has a shoulder. This slow relaxation process could be related to diffusive particle movements. The acoustic-type modes evidence a viscoelastic reaction with a 50% enhancement of the sound velocity. This enhancement of the sound velocity and the spectral line shape is qualitatively similar to spectra of molecular liquids. The two relaxation times of the memory function are separated by about two orders of magnitude and underpin the need for a wide frequency range investigation of this complex liquid. The high-frequency response can be interpreted as optic-type modes in the liquid.

We report computer simulations on the shear deformation of CuZr metallic glasses at zero and room temperatures. Shear bands emerge in athermal alloys at strain ${\gamma }_{\text{c}}$, with a finite-size effect found. The correlation of nonaffine displacement exhibits an exponential decay even after yielding in thermal alloys, but transits to a power law at $\gamma >{\gamma }_{\mathrm{c}}$ in athermal ones. The algebraic exponent is around $-1$ for the decay inside shear bands, consistent with the theoretical prediction in random elastic media. We quantify the anisotropic correlation with harmonic projection, finding the spectrum is weak in the exponential-decay regime, while it displays a strong polar and quadrupolar symmetry in the power-law regime. The nonvanishing quadrupolar symmetry at long distance signifies the nonlocality of plastic correlation in the athermal alloys. In contrast, the plastic correlation was found to be isotropic and localized at the yielding in the thermal alloys without shear bands.

The Derjaguin-Landau-Verywey-Overbeek (DLVO) theory has been a remarkably accurate framework for the characterization of macromolecular stability in water solvent. In view of its solvent-implicit nature neglecting the electrostatics of water molecules with non-negligible charge structure and concentration, the precision of the DLVO formalism is somewhat puzzling. In order to shed light on this issue, we derive from our earlier explicit solvent formalism [S. Buyukdagli et al., Phys. Rev. E 87, 063201 (2013)] a solvent-augmented contact value theorem and assess the contribution of solvent molecules to the interaction of charged membranes. We find that in the case of hydrophobic membranes with fixed charges embedded in the membrane surface, the nearly exact cancellation of various explicit solvent effects of substantially large magnitude but opposite sign keeps the intermembrane pressure significantly close to the double-layer force of the DLVO theory. Then, in the case of hydrophilic surface charge groups within the aqueous region, due to the spatial separation of the membrane substrate from the location of the fixed charges where the nonlocal dielectric response of the structured solvent is sharply localized, the interfacial field energy and the contact charge densities remain unaffected by the explicit solvent. As a result, the hydration of the lipid head groups suppresses the signature of the solvent molecules from the membrane interaction force.

Physical interactions generally respect certain symmetries, such as reciprocity and energy conservation, which survive in coarse-grained isothermal descriptions. Active many-body systems usually break such symmetries intrinsically, on the particle level, so that their collective behavior is often more naturally interpreted as a result of information exchange. Here we study numerically how information spreads from a “leader” particle through an initially aligned flock, described by the Vicsek model without noise. In the low-speed limit of a static spin lattice, we find purely conductive spreading, reminiscent of heat transfer. Swarm motility and heterogeneity can break reciprocity and spin conservation. But what seems more consequential for the swarm response is that the dispersion relation acquires a significant convective contribution along the leader's direction of motion.

High-speed particles overtaking the shock front during the propagation of compressional shocks in two-dimensional (2D) Yukawa solids are investigated using molecular dynamical simulations. When the compressional speed is lower, all particles around the shock front are almost accelerated synchronously. However, when the compressional speed is much higher, some particles penetrate the shock front to enter the preshock region. Around the shock front, it is found that the particle velocity profile at the first peak of the dispersive shock wave (DSW) is able to be described using the Gaussian distribution, so that the amplitudes of the DSW can be well characterized. As the compressional speed increases, the particle velocity corresponding to these DSW's amplitudes increase more substantially than the shock front speed. These amplitudes of the DSW are found to be able to predict the occurrence of the fast particles. Combined with the previous study of the DSW's period, it is demonstrated that the properties of the DSW are nearly not affected by the conditions of the 2D Yukawa systems, but only related to the compressional speed.

In this paper we examine the emergent structures of random networks that have undergone bond percolation an arbitrary, but finite, number of times. We define two types of sequential branching processes: a competitive branching process, in which each iteration performs bond percolation on the residual graph (RG) resulting from previous generations, and a collaborative branching process, where percolation is performed on the giant connected component (GCC) instead. We investigate the behavior of these models, including the expected size of the GCC for a given generation, the critical percolation probability, and other topological properties of the resulting graph structures using the analytically exact method of generating functions. We explore this model for Erdős-Renyi and scale-free random graphs. This model can be interpreted as a seasonal $N$-strain model of disease spreading.

Free-energy landscapes for short peptides—specifically for variants of the pH low insertion peptide (pHLIP)—in the heterogeneous environment of a lipid bilayer or cell membrane are constructed, taking into account a set of dominant interactions and the conformational preferences of the peptide backbone. Our methodology interprets broken internal H-bonds along the backbone of a polypeptide as statistically interacting quasiparticles, activated from the helix reference state. The favored conformation depends on the local environment (ranging from polar to nonpolar), specifically on the availability of external H-bonds (with ${\mathrm{H}}_2\mathrm{O}$ molecules or lipid headgroups) to replace internal H-bonds. The dominant side-chain contribution is accounted for by residue-specific transfer free energies between polar and nonpolar environments. The free-energy landscape is sensitive to the level of pH in the aqueous environment surrounding the membrane. For high pH, we identify pathways of descending free energy that suggest a coexistence of membrane-adsorbed peptides with peptides in solution. A drop in pH raises the degree of protonation of negatively charged residues and thus increases the hydrophobicity of peptide segments near the C terminus. For low pH, we identify insertion pathways between the membrane-adsorbed state and a stable trans-membrane state with the C terminus having crossed the membrane.

We consider an inertial active Ornstein-Uhlenbeck particle in an athermal bath. The particle is charged, constrained to move in a two-dimensional harmonic trap, and a magnetic field is applied perpendicular to the plane of motion. The steady-state correlations and the mean-square displacement are studied when the particle is confined as well as when it is set free from the trap. With the help of both numerical simulation and analytical calculations, we observe that inertia plays a crucial role in the dynamics in the presence of a magnetic field. In a highly viscous medium where the inertial effects are negligible, the magnetic field has no influence on the correlated behavior of position as well as velocity. In the time asymptotic limit, the overall displacement of the confined harmonic particle gets enhanced by the presence of a magnetic field and saturates for a stronger magnetic field. On the other hand, when the particle is set free, the overall displacement gets suppressed and approaches zero when the strength of the field is very high. Interestingly, it is seen that in the time asymptotic limit, the confined harmonic particle behaves like a passive particle and becomes independent of the activity, especially in the presence of a very strong magnetic field. Similarly, for a free particle the mean-square displacement in the long time limit becomes independent of activity even for a longer persistence of noise cor- relation in the dynamics.

Penetration of intruders into granular packings is well described by separately considering the dry or wet case of granular environments in previous experiments and simulations; however, the unified description of such penetration depth in these two granular media remains elusive due to lacking clear explanations about its origins. Based on three-dimensional discrete element method simulations, we introduce a power-law fitting form of the final penetration depth of a spherical intruder with low velocity vertically penetrating into dry and wet granular packings, excellently expressed on a master curve as a power-law function of a dimensionless impact number that is defined as the square root of the ratio between the inertial stress of the intruder and the linear combination of the mean gravitational stress and the cohesive stress exerted on each grain in the packings, as a remarkable extension of the inertial number in dry granular flows. This scaling robustly provides physical insights inherent in the unified description of the material properties of granular packings and the impactor penetration conditions on the final penetration depth in the impact tests, providing evidence of impact properties in different disciplines and applications in science and engineering.

The boundary treatment is fundamental for modeling fluid flows especially in the lattice Boltzmann method; the curved boundary conditions effectively improve the accuracy of single-phase simulations with complex-geometry boundaries. However, the conventional curved boundary conditions usually cause dramatic mass leakage or increase when they are directly used for multiphase flow simulations. We find that the principal reason for this is the absence of a nonideal effect in the curved boundary conditions, followed by a calculation error. In this paper, incorporating the nonideal effect into the linear interpolation scheme and compensating for the interpolating error, we propose a multiphase curved boundary condition to treat the wetting boundaries with complex geometries. A series of static and dynamic multiphase simulations with large density ratio verify that the present scheme is accurate and ensures mass conservation.

We show that application of boundary constraints generates unusual folding behaviors in responsive (swellable) helical bilayer strips. Unlike the smooth folding trajectories typical of free helical bilayers, the boundary-constrained bilayers exhibit intermittent folding behaviors characterized by rapid, steplike movements. We experimentally study bilayer strips as they swell and fold, and we propose a simple model to explain the emergence of ratchetlike behavior. Experiments and model predictions are then compared to simulations, which enable calculation of elastic energy during swelling. We investigate the dependence of this steplike behavior as a function of elastic boundary condition strength, strip length, and strip shape; interestingly, “V-shape” strips with the same boundary conditions fold smoothly.

Brilliantly colored butterflies and weevils are known to utilize photonic crystals for their coloration. Interestingly, the morphology of such crystals made of cuticle is based on triply periodic minimal surfaces such as gyroid and diamond surfaces. Recently, a different minimal-surface-based photonic crystal, the I-WP surface, was discovered inside the scale of a longhorn beetle. The letter I is derived from expressing the body center symmetry and WP is derived from a wrapped package. It was reported that the brilliant green color is produced by the photonic band gap existing along the [110] direction of this crystal. In this study, the polarization dependence of the reflection from this photonic crystal was investigated. A peculiar reflectance spectrum with two peaks was observed under the crossed polarizers. This characteristic is theoretically reproduced by calculating the reflectance from a finite-sized photonic crystal, and the spectral shape is explained based on the symmetry of the electromagnetic modes. In addition, inspired by this longhorn beetle, a photonic crystal structure consisting of colloidal particles is proposed, which has a similar polarization effect.

We know that a one-dimensional (1D) modulated system can simulate 2D topological states by expanding the dimension. This scenario provides a justifiable avenue to test the dilatation of the dynamical quantum phase transition (DQPT). Through a generalized Su-Schrieffer-Heeger model, we have shown how the Loschmidt echo, Fisher zero, and Dynamical topological order parameter (DTOP) transit from one to two dimensions. Owing to the introduced pseudomomentum, the derivative of the return rate does not always capture the DQPT well, but the Fisher zero and the DTOP can be treated as faithful indicators. A topology-independent parameter will also affect the occurrence of the DQPTs for quenches inside a given phase. Moreover, a comparison with the Haldane model owning the same phase diagram implies that a pair of fixed points will lead to different critical momentum distributions, thus different robustness, further reminding us that the correspondences between the equilibrium and dynamical phases transitions are multifarious.

Discontinuous phase transitions are particularly interesting from a social point of view because of their relationship to social hysteresis and critical mass. In this paper, we show that the replacement of a time-varying (annealed, situation-based) disorder by a static (quenched, personality-based) one can lead to a change from a continuous to a discontinuous phase transition. This is a result beyond the state of the art, because so far numerous studies on various complex systems (physical, biological, and social) have indicated that the quenched disorder can round or destroy the existence of a discontinuous phase transition. To show the possibility of the opposite behavior, we study a multistate $q$-voter model, with two types of disorder related to random competing interactions (conformity and anticonformity). We confirm, both analytically and through Monte Carlo simulations, that indeed discontinuous phase transitions can be induced by a static disorder.

Results and method of a three-loop renormalization-group calculation in the model of a Bose gas with a local density-density interaction in the formalism of time-dependent Green functions at finite temperature are presented. The results provide support to the recent conjecture [J. Honkonen, M. V. Komarova, Y. G. Molotkov, and M. Y. Nalimov, Nucl. Phys. B 939, 105 (2019); Y. A. Zhavoronkov, M. V. Komarova, Y. G. Molotkov, M. Y. Nalimov, and J. Honkonen, Theor. Math. Phys. 200, 1237 (2019)] that the dynamics of the superfluid phase transition is described by a model which belongs to the same universality class as the stochastic model A.

We examine how disordering joint position influences the linear elastic behavior of lattice materials via numerical simulations in two-dimensional beam networks. Three distinct initial crystalline geometries are selected as representative of mechanically isotropic materials with low connectivity, mechanically isotropic materials with high connectivity, and mechanically anisotropic materials with intermediate connectivity. Introducing disorder generates spatial fluctuations in the elasticity tensor at the local (joint) scale. Proper coarse-graining reveals a well-defined continuum-level scale elasticity tensor. Increasing disorder aids in making initially anisotropic materials more isotropic. The disorder impact on the material stiffness depends on the lattice connectivity: Increasing the disorder softens lattices with high connectivity and stiffens those with low connectivity, without modifying the scaling between elastic modulus and density (linear scaling for high connectivity and cubic scaling for low connectivity). Introducing disorder in lattices with intermediate fixed connectivity reveals both scaling: the linear scaling occurs for low density, the cubic one at high density, and the crossover density increases with disorder. Contrary to classical formulations, this work demonstrates that connectivity is not the sole parameter governing elastic modulus scaling. It offers a promising route to access novel mechanical properties in lattice materials via disordering the architectures.

We investigate the distributions of residence time for in-line chaotic mixers; in particular, we consider the Kenics, the F-mixer, and the multilevel laminating mixer and also a synthetic model that mimics their behavior and allows exact mathematical calculations. We show that whatever the number of elements of mixer involved, the distribution possesses a $t^{-3}$ tail, so that its shape is always far from Gaussian. This $t^{-3}$ tail also invalidates the use of second-order moment and variance. As a measure for the width of the distribution, we consider the mean absolute deviation and show that, unlike the standard deviation, it converges in the limit of large sample size. Finally, we analyze the performances of the different in-line mixers from the residence-time point of view when varying the number of elements and the shape of the cross section.

The interaction of light beams with helical defects in optical materials generates optical vortices. Understanding and manipulating the dynamics of helical defects allows for the creation of versatile sources of optical vortex beams. Using a magnetic ring on a nematic liquid crystal cell, we trapped helical defects identified as matter vortices. We observe oscillatory rotating and beating matter vortices by applying a low-frequency voltage. Experimentally, we determine the region of parameters where these vortices are observed. The amplitude of oscillatory rotating vortices decays with the inverse of the voltage frequency. We propose an adequate amplitude equation, which allows us to describe the vortex dynamics; theoretical findings have a qualitative agreement with the experimental observations.

Low-dimensional representations of underdamped systems often provide useful insights and analytical tractability. Here, we build such representations via information projections, obtaining an optimal model that captures the most information on observed spatial trajectories. We show that, in paradigmatic systems, the minimization of the information loss drives the appearance of a discontinuous transition in the optimal model parameters. Our results raise serious warnings for general inference approaches, and they unravel fundamental properties of effective dynamical representations impacting several fields, from biophysics to dimensionality reduction.

A Maxwell demon can reduce the entropy of a quantum system by performing measurements on its environment. The nonsignaling theorem prevents the demon from affecting the average state of the system. We study the preparations of quantum correlations from a system qubit and an auxiliary qubit, assisted by a demon who obtains information of the system qubit from measurements on its environment. The demon can affect the postmeasured states of system by choosing different measurements, which establishes the relationships between quantum steering and other correlations in the thermodynamic framework. We present the optimal protocols for creating mutual information, entanglement, and Bell-nonlocality. These maximal correlations are found to relate exactly to the steerable boundary of the system-environment state with maximally mixed marginals. We also present upper bounds of the prepared correlations by utilizing classical environment-system correlation, which can be regarded as steering-type inequalities bounding the correlations created with the aid of classical demons.