A series of recent publications, within the framework of network science, have focused on the coexistence of mixed attractive and repulsive (excitatory and inhibitory) interactions among the units within the same system, motivated by the analogies with spin glasses as well as to neural networks, or ecological systems. However, most of these investigations have been restricted to single layer networks, requiring further analysis of the complex dynamics and particular equilibrium states that emerge in multilayer configurations. This article investigates the synchronization properties of dynamical systems connected through multiplex architectures in the presence of attractive intralayer and repulsive interlayer connections. This setting enables the emergence of antisynchronization, i.e., intralayer synchronization coexisting with antiphase dynamics between coupled systems of different layers. We demonstrate the existence of a transition from interlayer antisynchronization to antiphase synchrony in any connected bipartite multiplex architecture when the repulsive coupling is introduced through any spanning tree of a single layer. We identify, analytically, the required graph topologies for interlayer antisynchronization and its interplay with intralayer and antiphase synchronization. Next, we analytically derive the invariance of intralayer synchronization manifold and calculate the attractor size of each oscillator exhibiting interlayer antisynchronization together with intralayer synchronization. The necessary conditions for the existence of interlayer antisynchronization along with intralayer synchronization are given and numerically validated by considering Stuart-Landau oscillators. Finally, we also analytically derive the local stability condition of the interlayer antisynchronization state using the master stability function approach.

The local microenvironment of a tumor plays an important and commonly observed role in cancer development and progression. Dynamic changes in the tissue microenvironment are thought to epigenetically disrupt healthy cellular phenotypes and drive cancer incidence. Despite the experimental work in this area there are no conceptual models to understand the interplay between the epigenetic dysregulation in the microenvironment of early tumors and the appearance of cancer driver mutations. Here, we develop a minimal model of the tissue microenvironment which considers three interacting subpopulations: healthy, phenotypically dysregulated, and mutated cancer cells. Healthy cells can epigenetically (reversibly) transition to the dysregulated phenotype, and from there to the cancer state. The epigenetic transition rates of noncancer cells can be influenced by the number of cancer cells in the microenvironment (termed microenvironment feedback). Our model delineates the regime in which microenvironment feedback accelerates the rate of cancer initiation. In addition, the model shows when and how microenvironment feedback may inhibit cancer progression. We discuss how our framework may provide resolution to some of the puzzling experimental observations of slow cancer progression.

We study the packing fraction of clusters in free-falling streams of spherical and irregularly shaped particles using flash x-ray radiography. The estimated packing fraction of clusters is low enough to correspond to coordination numbers less than 6. Such coordination numbers in numerical simulations correspond to aggregates that collide and grow without bouncing. Moreover, the streams of irregular particles evolved faster and formed clusters of larger sizes with lower packing fraction. This result on granular streams suggests that particle shape has a significant effect on the agglomeration process of granular materials.

By means of 3D particle dynamics simulations, we analyze the microstructure of granular materials subjected to isochoric (constant volume) cyclic shearing, which drives the system towards a liquefaction state characterized by loops of jamming-unjamming transition with periodic loss of strength and irreversible accumulation of shear strain. We first show that the macroscopic response obtained by these simulations agrees well with the most salient features of the well-known cyclic behavior of granular materials both before and after liquefaction. Then we investigate the evolution of particle connectivity, force transmission, and anisotropies of contact and force networks. The onset of liquefaction is marked by partial collapse of the force-bearing network with rapid drop of the coordination number and nonrattler fraction of particles, and significant broadening of the contact force probability density function, which begins in the preliquefaction period. We find that the jamming transition in each cycle occurs for a critical value of the coordination number that can be interpreted as the percolation threshold of the contact network and appears to be independent of the initial mean stress, void ratio, and cyclic shear amplitude. We show that upon unjamming in each cycle an isotropic loss of contacts occurs and is followed by the development of high contact anisotropy and a large proportion of particles with only two or three contacts. The higher mobility of the particles also involves a lower degree of frustration of particle rotations and thus lower friction mobilization and tangential force anisotropy. These findings are relevant to both undrained cyclic deformations of saturated soils and rheology of dense non-Brownian suspensions where volume change is coupled with pore liquid drainage conditions.

The generation of hot, directional electrons via laser-driven stimulated Raman scattering (SRS) is a topic of great importance in inertial confinement fusion (ICF) schemes. Little recent research has been dedicated to this process at high laser intensity, in which back, side, and forward scatter simultaneously occur in high energy density plasmas, of relevance to, for example, shock ignition ICF. We present an experimental and particle-in-cell (PIC) investigation of hot electron production from SRS in the forward and near-forward directions from a single speckle laser of wavelength ${\lambda }_0=1.053\mu \mathrm{m}$, peak laser intensities in the range $I_0=0.2–1.0\times {10}^{17}{\mathrm{Wcm}}^{-2}$ and target electron densities between $n_e=0.3–1.6\%n_c$, where $n_c$ is the plasma critical density. As the intensity and density are increased, the hot electron spectrum changes from a sharp cutoff to an extended spectrum with a slope temperature $T=34\pm 1\mathrm{keV}$ and maximum measured energy of 350 keV experimentally. Multidimensional PIC simulations indicate that the high energy electrons are primarily generated from SRS-driven electron plasma wave phase fronts with k vectors angled $\sim {50}^{\circ }$ with respect to the laser axis. These results are consistent with analytical arguments that the spatial gain is maximized at an angle which balances the tendency for the growth rate to be larger for larger scattered light wave angles until the kinetic damping of the plasma wave becomes important. The efficiency of generated high energy electrons drops significantly with a reduction in either laser intensity or target electron density, which is a result of the rapid drop in growth rate of Raman scattering at angles in the forward direction.

Understanding complex systems with their reduced model is one of the central roles in scientific activities. Although physics has greatly been developed with the physical insights of physicists, it is sometimes challenging to build a reduced model of such complex systems on the basis of insight alone. We propose a framework that can infer hidden conservation laws of a complex system from deep neural networks (DNNs) that have been trained with physical data of the system. The purpose of the proposed framework is not to analyze physical data with deep learning but to extract interpretable physical information from trained DNNs. With Noether's theorem and by an efficient sampling method, the proposed framework infers conservation laws by extracting the symmetries of dynamics from trained DNNs. The proposed framework is developed by deriving the relationship between a manifold structure of a time-series data set and the necessary conditions for Noether's theorem. The feasibility of the proposed framework has been verified in some primitive cases in which the conservation law is well known. We also apply the proposed framework to conservation law estimation for a more practical case, that is, a large-scale collective motion system in the metastable state, and we obtain a result consistent with that of a previous study.

The irreversible thermodynamics of a multicomponent fluid is reviewed. This includes a discussion of the role of individual component fluxes. It is argued that their differences vanish on the same time scale as that which establishes local thermodynamic equilibrium and thus do not play an independent role in fluid dynamics, but only arise in response to gradients in conserved thermodynamic variables. The contributions to the energy flux are examined and it is argued that there should be explicit contributions associated with the various component fluxes, which are not mentioned in standard kinetic theory presentations. Three different thermodynamic perspectives are discussed as to their form, with the respective equations for the entropy flux and production described and contrasted. The Onsager reciprocal relations are considered to be a consequence of the single-valuedness of the entropy production with the chemical potential gradients as the driving forces for diffusion. These are specialized to ideal gas mixtures using the component density gradients associated with Fick's laws and to using the mole fraction gradients that are standardly used in gas kinetic theory. The ideal gas Onsager relations are identical to those deduced from the Boltzmann equation. Irving and Kirkwood's statistical mechanics treatment of the evolution equations of a one-component fluid [J. Chem. Phys. 18, 817 (1950)] is generalized to multicomponent fluids and agrees with the thermodynamic perspective that treats the energy transfers as reversible.

We study relaxation of long-wavelength density perturbations in a one-dimensional conserved Manna sandpile. Far from criticality where correlation length $\xi$ is finite, relaxation of density profiles having wave numbers $k\to 0$ is diffusive, with relaxation time ${\tau }_R\sim k^{-2}/D$ with $D$ being the density-dependent bulk-diffusion coefficient. Near criticality with $k\xi \gtrsim 1$, the bulk diffusivity diverges and the transport becomes anomalous; accordingly, the relaxation time varies as ${\tau }_R\sim k^{-z}$, with the dynamical exponent $z=2-(1-\beta )/{\nu }_{\perp }<2$, where $\beta$ is the critical order-parameter exponent and ${\nu }_{\perp }$ is the critical correlation-length exponent. Relaxation of initially localized density profiles on an infinite critical background exhibits a self-similar structure. In this case, the asymptotic scaling form of the time-dependent density profile is analytically calculated: we find that, at long times $t$, the width $\sigma$ of the density perturbation grows anomalously, $\sigma \sim t^w$, with the growth exponent $\omega =1/(1+\beta )>1/2$. In all cases, theoretical predictions are in reasonably good agreement with simulations.

We introduce a model of a quantum walk on a graph in which a particle jumps between neighboring nodes and interacts with independent spins sitting on the edges. Entanglement propagates with the walker. We apply this model to the case of a one-dimensional lattice to investigate its magnetic and entanglement properties. In the continuum limit, we recover a Landau-Lifshitz equation that describes the precession of spins. A rich dynamics is observed, with regimes of particle propagation and localization, together with spin oscillations and relaxation. Entanglement of the asymptotic states follows a volume law for most parameters (the coin rotation angle and the particle-spin coupling).

Based on mean-field theory (MFT) arguments, a general description for discontinuous phase transitions in the presence of temporal disorder is considered. Our analysis extends the recent findings [C. E. Fiore et al., Phys. Rev. E 98, 032129 (2018)] by considering discontinuous phase transitions beyond those with a single absorbing state. The theory is exemplified in one of the simplest (nonequilibrium) order-disorder (discontinuous) phase transitions with “up-down” $Z_2$ symmetry: the inertial majority vote model for two kinds of temporal disorder. As for absorbing phase transitions, the temporal disorder does not suppress the occurrence of discontinuous phase transitions, but remarkable differences emerge when compared with the pure (disorderless) case. A comparison between the distinct kinds of temporal disorder is also performed beyond the MFT for random-regular complex topologies. Our work paves the way for the study of a generic discontinuous phase transition under the influence of an arbitrary kind of temporal disorder.

We study a recently proposed spin-1 model with competing antiferromagnetic first-neighbor interaction and a third-neighbor coupling mediated by nonmagnetic states, which reproduces topological features of the phase diagrams of high-$T_c$ superconductors [S. A. Cannas and D. A. Stariolo, Phys. Rev. E 99, 042137 (2019)]. We employ a cluster mean-field approach to investigate effects of crystal field anisotropy on the phase transitions hosted by this model. At low temperatures, the temperature-crystal field phase diagram exhibits superantiferromagnetic (SAF), antiferromagnetic (AF), and paramagnetic (PM) phases. In addition, we found a thermally driven state between SAF and PM phases. This thermally driven state and the SAF phase appears in the phase diagram as a domelike structure. Our calculations indicate that only second-order phase transitions occur in the PM-AF phase boundary, as suggested by previous Monte Carlo simulations.

Carreras, Dobson, and colleagues have studied empirical data on the sizes of the blackouts in real grids and modeled them with computer simulations using the direct current approximation. They have found that the resulting blackout sizes are distributed as a power law and suggested that this is because the grids are driven to the self-organized critical state. In contrast, more recent studies found that the distribution of cascades is bimodal resulting in either a very small blackout or a very large blackout, engulfing a finite fraction of the system. Here we reconcile the two approaches and investigate how the distribution of the blackouts changes with model parameters, including the tolerance criteria and the dynamic rules of failure of the overloaded lines during the cascade. In addition, we study the same problem for the Motter and Lai model and find similar results, suggesting that the physical laws of flow on the network are not as important as network topology, overload conditions, and dynamic rules of failure.

Phase-field models for strongly anisotropic surface energy need to be regularized to remove the ill posedness of the dynamic equations. Regularization introduces a new length scale, the corner size, also called the bending length. For large corner size, with respect to interface thickness, the phase-field method is known to converge asymptotically toward the sharp-interface theory when the appropriate approximation of the Willmore energy is used. In this work we study the opposite limit, i.e., corner size smaller than the interface width, and show that the shape of corners, at equilibrium, differs from the sharp-interface picture. However, we find that the phase transition at the interface is preserved and presents the same properties as the classical problem.

Critical properties of frictionless spherical particles below jamming are studied using extensive numerical simulations, paying particular attention to the nonaffine part of the displacements during the athermal quasistatic compression. It is shown that the squared norm of the nonaffine displacement exhibits a power-law divergence toward the jamming transition point. A possible connection between this critical exponent and that of the shear viscosity is discussed. The participation ratio of the displacements vanishes in the thermodynamic limit at the transition point, meaning that the nonaffine displacements are localized marginally with a fractal dimension. Furthermore, the distribution of the displacement is shown to have a power-law tail, the exponent of which is related to the fractal dimension.

The nature of thermonuclear explosions of white-dwarf stars is a fundamental astrophysical issue, the first principle interpretation of which is still commonly regarded as an unresolved problem. There is a general consensus that stellar explosions are a manifestation of the deflagration-to-detonation transition of an outward propagating self-accelerating thermonuclear flame subjected to instability-induced corrugations. A similar problem arises in unconfined terrestrial flames where a positive feedback mechanism leading to the pressure runaway has been identified. The present study is an application of this finding to the stellar environment. Notwithstanding a substantial modification of the equation of state the runaway effect endures. Approaching the runaway point the pretransition flame may stay perfectly subsonic, which challenges the view that to ensure the transition the flame should cross the threshold of Chapman-Jouguet deflagration.

Polymer molecules in a flow undergo a coil-stretch phase transition on an increase of the velocity gradients. Model-independent identification and characterization of the transition in a random flow has been lacking so far. Here we suggest to use the entropy of the extension statistics as a proper measure due to strong fluctuations around the transition. We measure experimentally the entropy as a function of the local Weisenberg number and show that it has a maximum, which identifies and quantifies the transition. We compare the new approach with the traditional one based on the theory using either linear Oldroyd-B or nonlinear finite extensible nonlinear elastic polymer models.

The random Lorentz gas (RLG) is a minimal model of transport in heterogeneous media that exhibits a continuous localization transition controlled by void space percolation. The RLG also provides a toy model of particle caging, which is known to be relevant for describing the discontinuous dynamical transition of glasses. In order to clarify the interplay between the seemingly incompatible percolation and caging descriptions of the RLG, we consider its exact mean-field solution in the infinite-dimensional $d\to \infty$ limit and perform numerics in $d=2...20$. We find that for sufficiently high $d$ the mean-field caging transition precedes and prevents the percolation transition, which only happens on timescales diverging with $d$. We further show that activated processes related to rare cage escapes destroy the glass transition in finite dimensions, leading to a rich interplay between glassiness and percolation physics. This advance suggests that the RLG can be used as a toy model to develop a first-principle description of particle hopping in structural glasses.

The emergent photoactive materials obtained through photochemistry make it possible to directly convert photon energy to mechanical work. There has been much recent work in developing appropriate materials, and a promising system is semicrystalline polymers of the photoactive molecule azobenzene. We develop a phase field model with two order parameters for the crystal-melt transition and the trans-cis photoisomerization to understand such materials, and the model describes the rich phenomenology. We find that the photoreaction rate depends sensitively on temperature: At temperatures below the crystal-melt transition temperature, photoreaction is collective, requires a critical light intensity, and shows an abrupt first-order phase transition manifesting nucleation and growth; at temperatures above the transition temperature, photoreaction is independent and follows first-order kinetics. Further, the phase transition depends significantly on the exact forms of spontaneous strain during the crystal-melt and trans-cis transitions. A nonmonotonic change of photopersistent cis ratio with increasing temperature is observed accompanied by a reentrant crystallization of trans below the melting temperature. A pseudo phase diagram is subsequently presented with varying temperature and light intensity along with the resulting actuation strain. These insights can assist the further development of these materials.

We consider a thin fluid film flowing down an inclined substrate subjected to localized external sources of momentum and heat flux that induce deformations of the fluid's free surface. This scenario is encountered in several industrial processes and of particular interest is the case where these deformations are undesirable. When the substrate is thin and the temperature along its underside is freely imposed by an active cooling mechanism, temperature gradients are generated at the fluid surface which drive a thermocapillary flow and influence the deformations. This naturally leads us to pose the optimal control problem of choosing the temperature profile that minimizes the unwanted free-surface deformations. Numerical computations reveal that the external forces generate deflections in a region near their peak beyond which all deflections are suppressed by the optimal control. Where nonzero deflections occur, the control is of bang-bang type (taking either its upper or lower bound), while the control is obtained in closed form for regions where the deflections are suppressed. Strikingly, in switching between these regions the optimal control chatters, that is, it switches infinitely many times over a finite interval. By appealing to Pontryagin's maximum principle and leveraging a symmetry embedded in the adjoint problem we uncover the underlying fractal structure of the chattering. Finally, we present practical approaches to avoid the infinite switching while retaining significantly reduced free-surface deformations.