Following numerous earlier studies, extensive simulations and analyses were made on the continuous interaction distribution Gaussian model and the discrete bimodal interaction distribution Ising spin glass (ISG) models in two dimensions [Lundow and Campbell, Phys. Rev. E 93, 022119 (2016)]. Here we further analyze the bimodal and Gaussian data together with data on two other continuous interaction distribution two-dimensional ISG models, the uniform and the Laplacian models, and three other discrete interaction distribution models, a diluted bimodal model, an “antidiluted” model, and a more exotic symmetric Poisson model. Comparisons between the three continuous distribution models show that not only do they share the same exponent $\eta \equiv 0$ but that to within the present numerical precision they share the same critical exponent $\nu$ also, and so lie in a single universality class. On the other hand the critical exponents of the four discrete distribution models are not the same as those of the continuous distributions, and the present data strongly indicate that they differ from one discrete distribution model to another. This is evidence that discrete distribution ISG models in two dimensions have nonzero values of the critical exponent $\eta$ and do not lie in a single universality class.

We have determined the thermal conductance of a system consisting of a two-level atom coupled to two quantum harmonic oscillators in contact with heat reservoirs at distinct temperatures. The calculation of the heat flux as well as the atomic population and the rate of entropy production are obtained by the use of a quantum Fokker-Planck-Kramers equation and by a Lindblad master equation. The calculations are performed for small values of the coupling constant. The results coming from both approaches show that the conductance is proportional to the coupling constant squared and that, at high temperatures, it is proportional to the inverse of temperature.

We present general criteria for the occurrence of infinite avalanches and critical hysteresis in the zero-temperature nonequilibrium random-field Ising model on a Bethe lattice. Drawing upon extant results as well as a result on a dilute four-coordinated ($z=4$) lattice, we show that diverging avalanches can occur if an arbitrarily small fraction of sites on a spanning cluster have connectivity $z\ge 4$.

The formalism of complex networks is extensively employed to describe the dynamics of interacting agents in several applications. The features of the connections among the nodes in a network are not always provided beforehand, hence the problem of appropriately inferring them often arises. Here, we present a method to reconstruct directed and weighted topologies of networks of heterogeneous nonlinear oscillators. We illustrate the theory on a set of representative examples.

We study the low-frequency polarization response of a surface-charged oblate spheroidal particle immersed in an electrolyte solution. Because the charged spheroid attracts counterions which form the electric double layer around the particle, using usual boundary conditions at the interface between the particle and electrolyte can be quite complicated and challenging. Hence, we generalize Fixman's boundary conditions, originally derived for spherical particles, to the case of the charged oblate spheroid. Given two different counterion distributions in the thin electric double-layer limit, we obtain analytic expressions for the polarization coefficients to the first nontrivial order in frequency. We find that the polarization response normal to the symmetry axis depends on the total amount of charge carried by the oblate spheroid while that parallel to the symmetry axis is suppressed when there is less charge on the edge of the spheroid. We further study the overall dielectric response for a dilute suspension of charged spheroids. We find that the dielectric enhancement at low frequency, which is driven by the presence of a large $\zeta$ potential (surface charge), is suppressed by high ion concentrations in the electrolyte and depends on the size of the suspended particles. In addition, spheroids with higher aspect ratios will also lead to a stronger dielectric enhancement due to the combination of the electric double layer and textural effects. The characteristic frequency associated with the dielectric enhancement scales inversely with the square of the particle size, the major radius of the spheroid, and it has a weak dependence on the shape of spheroids.

After a raindrop impacts on a granular bed, a crater is formed as both drop and target deform. After an initial, transient, phase in which the maximum crater depth is reached, the crater broadens outwards until a final steady shape is attained. By varying the impact velocity of the drop and the packing density of the bed, we find that avalanches of grains are important in the second phase and hence affect the final crater shape. In a previous paper, we introduced an estimate of the impact energy going solely into sand deformation and here we show that both the transient and final crater diameter collapse with this quantity for various packing densities. The aspect ratio of the transient crater is however altered by changes in the packing fraction.

A doughnut-shaped object supporting surface rotations was a hypothetical construct proposed by both Taylor and Purcell as a swimmer that would be able to propel itself in a Stokesian fluid because of the irreversibility of its stroke. Here we numerically examine the hydrodynamic interaction of pairs and trios of these free toroidal swimmers. First, we study the axisymmetric case of two toroidal swimmers placed in tandem, and show that a single torus of a corotating pair is more efficient than when it swims alone, but less efficient when paired with a counterrotating partner. Using a regularized Stokeslet framework, we study the nonaxisymmetric case of toroidal swimmers whose axes are initially parallel, but not collinear. These perturbed in tandem swimmers can exhibit qualitatively different trajectories that may, for instance, repel the swimmers or have them settle into a periodic state. We also illustrate interesting dynamics that occur for different initial configurations of three tori.

We study the problem of predictability, or “nature vs nurture,” in several disordered Ising spin systems evolving at zero temperature from a random initial state: How much does the final state depend on the information contained in the initial state, and how much depends on the detailed history of the system? Our numerical studies of the “dynamical order parameter” in Edwards-Anderson Ising spin glasses and random ferromagnets indicate that the influence of the initial state decays as dimension increases. Similarly, this same order parameter for the Sherrington-Kirkpatrick infinite-range spin glass indicates that this information decays as the number of spins increases. Based on these results, we conjecture that the influence of the initial state on the final state decays to zero in finite-dimensional random-bond spin systems as dimension goes to infinity, regardless of the presence of frustration. We also study the rate at which spins “freeze out” to a final state as a function of dimensionality and number of spins; here the results indicate that the number of “active” spins at long times increases with dimension (for short-range systems) or number of spins (for infinite-range systems). We provide theoretical arguments to support these conjectures, and also study analytically several mean-field models: the random energy model, the uniform Curie-Weiss ferromagnet, and the disordered Curie-Weiss ferromagnet. We find that for these models, the information contained in the initial state does not decay in the thermodynamic limit—in fact, it fully determines the final state. Unlike in short-range models, the presence of frustration in mean-field models dramatically alters the dynamical behavior with respect to the issue of predictability.

In the absence of external driving, a system exposed to thermal fluctuations will relax to equilibrium. However, the constant input of work makes it possible to counteract this relaxation and maintain the system in a nonequilibrium steady state. In this article, we use the stochastic thermodynamics of Markov jump processes to compute the minimum rate at which energy must be supplied and dissipated to maintain an arbitrary nonequilibrium distribution in a given energy landscape. This lower bound depends on two factors: the undriven probability current in the equilibrium state and the distance from thermal equilibrium of the target distribution. By showing the consequences of this result in a few simple examples, we suggest general implications for the required energetic costs of macromolecular repair and cytosolic protein localization.

We investigate feedback processes with measurement-induced protocols for particular tasks that drive systems in specified directions in state spaces. We focus on mutual information as a measure of correlation between system and memory, which has been known to play a crucial role for the second law of information thermodynamics. The performance of task is enhanced in the early stage of driving, along with the decrease of correlation and mutual information due to the passage from initial measurement. However, we find that the performance is suppressed if the time of driving exceeds a threshold, which we call feedback overshooting. We find that a type of correlation, anticorrelation, between system and memory is built up as a result of overshooting and gives rise to regaining mutual information. We examine the effect of overshooting in detail from two examples. We study the Szilard engine for the task of work extraction. We also study a recurrent feedback with finite time interval for the task to reduce the mean square distance of a colloid below the value by thermal fluctuation. We find that recurrent feedback is stable only for a moderate range of time intervals and the intensity of feedback protocol. We discuss the problem of divergence of mutual information for error-free measurement.

Quantum resonant activation is investigated for the archetype setup of an externally driven two-state (spin-boson) system subjected to strong dissipation by means of both analytical and extensive numerical calculations. The phenomenon of resonant activation emerges in the presence of either randomly fluctuating or deterministic periodically varying driving fields. Addressing the incoherent regime, a characteristic minimum emerges in the mean first passage time to reach an absorbing neighboring state whenever the intrinsic time scale of the modulation matches the characteristic time scale of the system dynamics. For the case of deterministic periodic driving, the first passage time probability density function (pdf) displays a complex, multipeaked behavior, which depends crucially on the details of initial phase, frequency, and strength of the driving. As an interesting feature we find that the mean first passage time enters the resonant activation regime at a critical frequency ${\nu }^{*}$ which depends very weakly on the strength of the driving. Moreover, we provide the relation between the first passage time pdf and the statistics of residence times.

Spontaneous symmetry breaking is an important phenomenon observed in various fields including physics and biology. In this connection, we here show that the trade-off between attractive and repulsive couplings can induce spontaneous symmetry breaking in a homogeneous system of coupled oscillators. With a simple model of a system of two coupled Stuart-Landau oscillators, we demonstrate how the tendency of attractive coupling in inducing in-phase synchronized (IPS) oscillations and the tendency of repulsive coupling in inducing out-of-phase synchronized oscillations compete with each other and give rise to symmetry breaking oscillatory states and interesting multistabilities. Further, we provide explicit expressions for synchronized and antisynchronized oscillatory states as well as the so called oscillation death (OD) state and study their stability. If the Hopf bifurcation parameter ($\lambda$) is greater than the natural frequency ($\omega$) of the system, the attractive coupling favors the emergence of an antisymmetric OD state via a Hopf bifurcation whereas the repulsive coupling favors the emergence of a similar state through a saddle-node bifurcation. We show that an increase in the repulsive coupling not only destabilizes the IPS state but also facilitates the reentrance of the IPS state.

One of the properties that make ecological systems so unique is the range of complex behavioral patterns that can be exhibited by even the simplest communities with only a few species. Much of this complexity is commonly attributed to stochastic factors that have very high-degrees of freedom. Orthodox study of the evolution of these simple networks has generally been limited in its ability to explain complexity, since it restricts evolutionary adaptation to an inertia-free process with few degrees of freedom in which only gradual, moderately complex behaviors are possible. We propose a model inspired by particle-mediated field phenomena in classical physics in combination with fundamental concepts in adaptation, which suggests that small but high-dimensional chaotic dynamics near to the adaptive trait optimum could help explain complex properties shared by most ecological datasets, such as aperiodicity and pink, fractal noise spectra. By examining a simple predator-prey model and appealing to real ecological data, we show that this type of complexity could be easily confused for or confounded by stochasticity, especially when spurred on or amplified by stochastic factors that share variational and spectral properties with the underlying dynamics.

Presented here are observations of the outcomes of the collapses of large single bubbles in ${\mathrm{H}}_2\mathrm{O}$ and ${\mathrm{D}}_2\mathrm{O}$ at high ambient pressures. Experiments were carried out in a high-pressure spherical resonator at ambient pressures of up to 30 MPa and acoustic pressures up to 35 MPa. Monitoring of the collapse events and their outcomes was accomplished using multiframe high-speed photography. Among the observations to be presented are the temporal and spatial evolution of light emissions produced by the collapse events, which were observed to last on the order of 30 ns and have time independent radii on the order of $30\mu \mathrm{m}$; the production of Rayleigh-Taylor jets which were observed to travel distances of up to $70\mu \mathrm{m}$ at speeds in excess of 4500 m/s; the entrainment of the light emitting regions in the jets' remnants; the production of spheroidal objects around the collapse points of the bubbles, far from any surface of the resonator; and the traversal and emergence of the Rayleigh-Taylor jets through the spherical objects. These spheroidal objects appear to behave as amorphous solids and form at locations where hydrodynamics predicts pressures in excess of the known transition pressures of water into the high-pressure crystalline ices, Ice-VI and Ice-VII.

Thorough numerical studies reveal that spatially extended dissipative systems with long-range interactions may give rise to a large-scale dynamics. This phenomenon, which generalizes mean-field chaos, can be interpreted as a form of subtle pattern formation, where a chaotic microscopic dynamics coexists with a macroscopic irregular behavior, sustained by the spontaneous emergence of long-wavelength “hydrodynamic” modes. This regime can emerge only if the coupling is sufficiently long ranged, otherwise normal space-time chaos is observed. In Stuart-Landau oscillators a further regime is found, where the amplitude of the hydrodynamic modes exhibits an anomalous intermediate scaling between that of collective and standard space-time chaos.

In the literature, different theoretical models have been proposed to describe the properties of systems which consist of magnetizable particles that are embedded into an elastomer matrix. It is well known that such magneto-sensitive elastomers display a strong magneto-mechanical coupling when subjected to an external magnetic field. Nevertheless, the predictions of available models often vary significantly since they are based on different assumptions and approximations. Up to now the actual accuracy and the limits of applicability are widely unknown. In the present work, we compare the results of a microscale continuum and a dipolar mean field approach with regard to their predictions for the magnetostrictive response of magneto-sensitive elastomers and reveal some fundamental relations between the relevant quantities in both theories. It turns out that there is a very good agreement between both modeling strategies, especially for entirely random microstructures. In contrast, a comparison of the finite-element results with a modified approach, which—similar to the continuum model—is based on calculations with discrete particle distributions, reveals clear deviations. Our systematic analysis of the differences shows to what extent the dipolar mean field approach is superior to other dipole models.

Using the two-temperature model for ultrafast matter (UFM), we compare the equation of state, pair-distribution functions $g\left(r\right)$, and phonons using the neutral pseudoatom (NPA) model with results from density functional theory (DFT) codes and molecular dynamics (MD) simulations for Al, Li, and Na. The NPA approach uses state-dependent first-principles pseudopotentials from an “all-electron” DFT calculation with finite-$T$ exchange-correlation functional (XCF). It provides pair potentials, structure factors, the “bound” and “free” states, as well as a mean ionization $\overline{Z}$ unambiguously. These are not easily accessible via DFT+MD calculations which become prohibitive for $T/T_F$ exceeding $\sim 0.6$, where $T_F$ is the Fermi temperature. Hence, both DFT+MD and NPA methods can be compared up to $\sim 8\mathrm{eV}$, while higher $T$ can be addressed via the NPA. The high-$T_e$ phonon calculations raise the question of UFM lattice stability and surface ablation in thin UFM samples. The ablation forces in a UFM slab are used to define an “ablation time” competing with phonon formation times in thin UFM samples. Excellent agreement for all properties is found between NPA and standard DFT codes, even for Li where a strongly nonlocal pseudopotential is used in DFT codes. The need to use pseudopotentials appropriate to the ionization state $\overline{Z}$ is emphasized. The effect of finite-$T$ XCF is illustrated via its effect on the pressure and the electron-density distribution at a nucleus.

In this paper, different force treatments are analyzed in detail for a pseudopotential lattice Boltzmann equation (LBE), and the contribution of third-order error terms to pressure tensor with a force scheme is analyzed by a higher-order Chapman-Enskog expansion technique. From the theoretical analysis, the performance of the original force treatment of Shan-Chen (SC), Ladd, Guo et al., and the exact difference method (EDM) are ${\varepsilon }_{\mathrm{Ladd}}<{\varepsilon }_{\mathrm{Guo}}<{\varepsilon }_{\mathrm{EDM}}\le {\varepsilon }_{\mathrm{SC}}$ with the relaxation time $\tau \ge 1$, while ${\varepsilon }_{\mathrm{Ladd}}<{\varepsilon }_{\mathrm{Guo}}<{\varepsilon }_{\mathrm{SC}}<{\varepsilon }_{\mathrm{EDM}}$ with $\tau <1$; here $\varepsilon$ is a parameter related to the mechanical stability and the subscripts are the corresponding force scheme. To be consistent with the thermodynamic theory, a force term is introduced to modify the coefficients in the pressure tensor. Some numerical simulations are conducted to show that the predictions of modified force treatment of the pseudopotential LBE are all in good agreement with the analytical solution and other predictions.

We show that chimera states, where differentiated subsets of synchronized and desynchronized dynamical elements coexist, can emerge in networks of hyperbolic chaotic oscillators subject to global interactions. As local dynamics we employ Lozi maps, which possess hyperbolic chaotic attractors. We consider a globally coupled system of these maps and use two statistical quantities to describe its collective behavior: the average fraction of elements belonging to clusters and the average standard deviation of state variables. Chimera states, clusters, complete synchronization, and incoherence are thus characterized on the space of parameters of the system. We find that chimera states are related to the formation of clusters in the system. In addition, we show that chimera states arise for a sufficiently long range of interactions in nonlocally coupled networks of these maps. Our results reveal that, under some circumstances, hyperbolicity does not impede the formation of chimera states in networks of coupled chaotic systems, as it had been previously hypothesized.

Suspensions of cornstarch in water exhibit strong dynamic shear thickening. We show that partly replacing water with ethanol strongly alters the suspension rheology. We perform steady and nonsteady rheology measurements combined with atomic force microscopy to investigate the role of fluid chemistry on the macroscopic rheology of the suspensions and its link with the interactions between cornstarch grains. Upon increasing the ethanol content, the suspension goes through a yield-stress fluid state and ultimately becomes a shear-thinning fluid. On the cornstarch grain scale, atomic force microscopy measurements reveal the presence of polymers on the cornstarch surface, which exhibit a cosolvency effect. At intermediate ethanol content, a maximum of polymer solubility induces high microscopic adhesion which we relate to the macroscopic yield stress.