We study chaotic dynamics and anomalous transport in a Bose-Hubbard chain in the semiclassical regime (the limit when the number of particles goes to infinity). We find that the system has mixed phase space with both regular and chaotic dynamics, even for long chains with up to 100 wells. The consequence of the mixed phase space is strongly anomalous diffusion in the space of occupation numbers, with a discrete set of transport exponents. After very long times the system crosses over to the hydrodynamic regime with normal diffusion. Anomalous transport is quite universal and almost completely independent of the parameters of the model (Coulomb interaction and chemical potential): It is mainly determined by the initial distribution of particles along the chain. We corroborate our findings by analytical arguments: scaling analysis for the anomalous regime and the Langevin equation for the normal diffusion regime.

Chemical reaction fronts separate regions of reacted and unreacted substances as they propagate in liquids. These fronts may induce density gradients due to different chemical compositions and temperatures across the front. In this paper, we investigate buoyancy-induced convection driven by both types of gradients. We consider a thin front approximation where the normal front velocity depends only on the front curvature. This model applies for small curvature fronts independent of the specific type of chemical reaction. For density changes due only to heat variations near the front, we find that convection can take place for either upward or downward propagating fronts if density gradients are above a threshold. Convection can set in even if the fluid with lower density is above the higher density fluid. Our model consists of Navier-Stokes equations coupled to the front propagation equation. We carry out a linear stability analysis to determine the parameters for the onset of convection. We study the nonlinear front propagation for liquids confined in narrow two-dimensional domains. Convection leads to fronts of steady shape, propagating with constant velocities.

Real-world water wave fields exhibit significant nonlinear and nonisospectral characteristics, making it challenging to predict their evolution by relying solely on numerical simulation or exact solutions using integrable system theory. Hence, this paper introduces a fast and adaptive method of modal identification and prediction in nonisospectral water wave fields using the reduced-order nonlinear solution (RONS) scheme. Specifically, we discuss the coarse graining and mode extraction of wave field snapshots from the data-driven and physics-driven perspectives and utilize the RONS method for principle modal prediction of nonisospectral water wave fields. This is achieved by investigating the standard and nonisospectral Gardner system describing nonlinear water waves as a demonstration. Through detailed comparison and analysis, the fundamental solitary behaviors and dispersive effects in the Gardner system are discussed. Subsequently, a neighbor approximation is developed that combines the essences of symbolic precomputation and numerical computation in the RONS procedure, which exploits the locality of nonlinear interactions in water wave fields.

The existence and properties of envelope solitary waves on a periodic traveling-wave background, called traveling breathers, are investigated numerically in representative nonlocal dispersive media. Using a fixed-point computational scheme, a space-time boundary-value problem for bright traveling breather solutions is solved for the weakly nonlinear Benjamin-Bona-Mahony equation, a nonlocal, regularized shallow water wave model, and the strongly nonlinear conduit equation, a nonlocal model of viscous core-annular flows. Curves of unit-mean traveling breather solutions within a three-dimensional parameter space are obtained. Resonance due to nonconvex, rational linear dispersion leads to a nonzero oscillatory background upon which traveling breathers propagate. These solutions exhibit a topological phase jump and so act as defects within the periodic background. For small amplitudes, traveling breathers are well approximated by bright soliton solutions of the nonlinear Schrödinger equation with a negligibly small periodic background. These solutions are numerically continued into the large-amplitude regime as elevation defects on cnoidal or cnoidal-like periodic traveling-wave backgrounds. This study of bright traveling breathers provides insight into systems with nonconvex, nonlocal dispersion that occur in a variety of media such as internal oceanic waves subject to rotation and short, intense optical pulses.

The statistics of the diffusive motion of particles often serve as an experimental proxy for their interaction with the environment. However, inferring the physical properties from the observed trajectories is challenging. Inspired by a recent experiment, here we analyze the problem of particles undergoing two-dimensional Brownian motion with transient tethering to the surface. We model the problem as a hidden Markov model where the physical position is observed and the tethering state is hidden. We develop an alternating maximization algorithm to infer the hidden state of the particle and estimate the physical parameters of the system. The crux of our method is a saddle-point-like approximation, which involves finding the most likely sequence of hidden states and estimating the physical parameters from it. Extensive numerical tests demonstrate that our algorithm reliably finds the model parameters and is insensitive to the initial guess. We discuss the different regimes of physical parameters and the algorithm's performance in these regimes. We also provide a free software implementation of our algorithm.

Liquid state entropy formulas based on configurational probability distributions are examined for Lennard-Jones fluids across a range temperatures and densities. These formulas are based on expansions of the entropy in a series of $n$-body distribution functions. We focus on two special cases. One, which we term the “perfect gas” series, starts with the entropy of an ideal gas; the other, which we term the “dense liquid” series, removes a many-body contribution from the ideal gas entropy and reallocates it among the subsequent $n$-body terms. We show that the perfect gas series is most accurate at low density, while the dense liquid series is most accurate at high density. We propose empirical interpolation methods that are capable of connecting the two series and giving consistent predictions in most situations.

We study a nonequilibrium ferromagnetic mean-field spin model exhibiting a phase with spontaneous temporal oscillations of the magnetization, on top of the usual paramagnetic and ferromagnetic phases. This behavior is obtained by introducing dynamic field variables coupled to the spins through nonreciprocal couplings. We determine a nonequilibrium generalization of the Landau free energy in terms of the large deviation function of the magnetization and of an appropriately defined smoothed stochastic time derivative of the magnetization. While the transition between paramagnetic and oscillating phase is continuous, the transition between ferromagnetic and oscillating phases is found to be discontinuous, with a coexistence of both phases, one being stable and the other one metastable. Depending on parameter values, the ferromagnetic points may either be inside or outside the limit cycle, leading to different transition scenarios. The stability of these steady states is determined from the large deviation function. We also show that in the coexistence region, the entropy production has a pronounced maximum as a function of system size.

The motivation for the proposed work is drawn from the attachment-detachment observed in biological and physical transport processes that entail finite resources. We investigate the influence of limited particle availability on particle dynamics within two parallel totally asymmetric simple exclusion lanes, with one lane incorporating only particle detachment and the other considering particle attachment. We establish a theoretical framework by employing vertical mean-field theory in conjunction with singular perturbation technique. The analytical findings are supported by numerical and stochastic validation using a finite-difference scheme and the Gillespie algorithm. By utilizing these approaches, we scrutinize various stationary properties, including particle densities, phase boundaries, and particle currents for both lanes. Our analysis reveals that the complexity of the phase diagram exhibits a nonmonotonic trend in the number of stationary phases as the particle count increases. Each phase diagram is constructed with respect to the intrinsic boundary parameters, illustrating both bulk and surface transitions occurring within the lanes. The interplay between finite resources and coupling mechanisms gives rise to two phases involving upward shock in one of the lanes, while two phases exhibit synchronized downward shock in both lanes. Finally, we delve into shock dynamics to comprehend critical phase transitions occurring in the system.

Systems featuring hard-core–soft-shell repulsive pair potentials can form ordered phases, where particles organize themselves in aggregates with nontrivial geometries. The dimer crystal formed by one such potential, namely, the hard-core plus generalized exponential model of order 4, has been recently investigated, revealing a low-temperature structural phase transition, with the onset of nematic ordering of the dimers. In the present work, we aim to characterize this phase transition via a mean-field theory, by which a detailed analysis of the low-temperature properties of the system is carried out under quadrupole approximation. We determine the transition temperature and identify its order parameter, highlighting the link between the structural transition and the nematic ordering of the system. The first-order character of the transition is established and supported by the Landau expansion of the free energy in powers of the order parameter. The theory is subsequently generalized to take into account lattice vibrations and dimer length fluctuations. Finally, we provide an explanation for the anomalous behavior displayed by the specific heat in the vanishing-temperature limit, which is also supported by Monte Carlo simulations.

We analytically examine fluctuations of vorticity excited by an external random force in two-dimensional fluid. We develop the perturbation theory enabling one to calculate nonlinear corrections to correlation functions of the flow fluctuations found in the linear approximation. We calculate the correction to the pair correlation function and the triple correlation function. It enables us to establish the criterion of validity of the perturbation theory for different ratios of viscosity and bottom friction. We find that the corrections to the second moment are anomalously weak in the cases of small bottom friction and small viscosity and relate the weakness to the energy and enstrophy balances. We demonstrate that at small bottom friction the triple correlation function is characterized by universal scaling behavior in some region of lengths. The developed perturbation method was verified and confirmed by direct numerical simulations.

We characterize universal features of the sample-to-sample fluctuations of global geometrical observables, such as the area, width, length, and center-of-mass position, in random growing planar clusters. Our examples are taken from simulations of both continuous and discrete models of kinetically rough interfaces, including several universality classes, such as Kardar-Parisi-Zhang. We mostly focus on the scaling behavior with time of the sample-to-sample deviation for those global magnitudes, but we have also characterized their histograms and correlations.

We analyze the stationary current of Bose particles across the Bose-Hubbard chain connected to a battery, focusing on the effect of interparticle interactions. It is shown that the current magnitude drastically decreases as the strength of interparticle interactions exceeds the critical value which marks the transition to quantum chaos in the Bose-Hubbard Hamiltonian. We found that this transition is well reflected in the nonequilibrium many-body density matrix of the system. Namely, the level-spacing distribution for eigenvalues of the density matrix changes from Poisson to Wigner-Dyson distributions. With the further increase of the interaction strength, the Wigner-Dyson spectrum statistics change back to the Poisson statistics which now marks fermionization of the Bose particles. With respect to the stationary current, this leads to the counter-intuitive dependence of the current magnitude on the particle number.

Besides its original spin representation, the Ising model is known to have the Fortuin-Kasteleyn (FK) bond and loop representations, of which the former was recently shown to exhibit two upper critical dimensions $(d_c=4,d_p=6)$. Using a lifted worm algorithm, we determine the critical coupling as $K_c=0.07770891\left(4\right)$ for $d=7$, which significantly improves over the previous results, and then study critical geometric properties of the loop Ising clusters on tori for spatial dimensions $d=5$ to 7. We show that as the spin representation, the loop Ising model has only one upper critical dimension at $d_c=4$. However, sophisticated finite-size scaling (FSS) behaviors, such as two length scales, two configuration sectors, and two scaling windows, still exist as the interplay effect of the Gaussian fixed point and complete-graph asymptotics. Moreover, using the loop-cluster algorithm, we provide an intuitive understanding of the emergence of the percolation-like upper critical dimension $d_p=6$ in the FK-Ising model. As a consequence, a unified physical picture is established for the FSS behaviors in all three representations of the Ising model above $d_c=4$.

The recently proposed Ehrenfest $M$-urn model with interactions on a ring is considered as a paradigm model which can exhibit a variety of distinct nonequilibrium steady states. Unlike the previous three-urn model on a ring which consists of a uniform steady state and a nonuniform nonequilibrium steady state, it is found that for even $M\ge 4$, an additional nonequilibrium steady state can coexist with the original ones. Detailed analysis reveals that this additional nonequilibrium steady state emerged via a pitchfork bifurcation which cannot occur if $M$ is odd. Properties of this nonequilibrium steady state, such as stability, and steady-state flux are derived analytically for the four-urn case. The full phase diagram with the phase boundaries is also derived explicitly. The associated thermodynamic stability is also analyzed, confirming its stability. These theoretical results are also explicitly verified by direct Monte Carlo simulations for the three-urn and four-urn ring models.

War and peace, spanning history, deeply affect society, economy, and individuals. Grasping their dynamics is vital to lessen harm and foster global peace. Yet, quantifying them remains hard. Our goal is to create a simple qualitative model using landscape and flux theory, exploring war and peace mechanisms. In this symmetric network, they appear as separate attraction basins, dynamically shifting. Analyzing landscape shape gives insights into global stability. Near critical points, indicators like cross correlations, autocorrelation times, and flickering frequency surge, as warnings. We also calculate the irreversible path between war and peace due to rotational flux. Global sensitivity analysis identifies history's role in system stability. In summary, our research unveils a way to understand war and peace complexities, enhancing knowledge of key elements that lead to conflict, aiding resolution.

We develop and study the hydrodynamic theory of flocking with autochemotaxis. This describes large collections of self-propelled entities all spontaneously moving in the same direction, each emitting a substance which attracts the others (e.g., ants). The theory combines features of the Keller-Segel model for autochemotaxis with the Toner-Tu theory of flocking. We find that sufficiently strong autochemotaxis leads to an instability of the uniformly moving state (the “flock”), in which bands of different density form moving parallel to the mean flock velocity with different speeds. This instability is, therefore, completely different from the well-known “banding instability,” in which bands form perpendicular to the mean flock velocity. The bands we find, which are reminiscent of ant trails, coarsen over time to reach a phase-separated state, in which one high-density and one low-density band fill the entire system. The same instability, described by the same hydrodynamic theory, can occur in flocks phase separating due to any microscopic mechanism (e.g., sufficiently strong attractive interactions). Although in many ways analogous to equilibrium phase separation via spinodal decomposition, the two steady-state densities here are determined not by a common tangent construction, as in equilibrium, but by an $\mathit{un}$common tangent construction very similar to that found for motility-induced phase separation of disordered active particles. Our analytic theory agrees well with our numerical simulations of our equations of motion.

We use a standard minimal active Brownian model to investigate the osmotic-like effective pressure generated by active fluids on fixed hollow inclusions. These inclusions are enclosed by a permeable (albeit nonflexible) membrane, and the interior and exterior regions of the inclusions have different particle motility strengths. We consider both rectangular and disklike inclusions and analyze the effects of various system parameters, such as excluded volume interaction between active particles, hardness of membrane, and active particle density, on the effective pressure produced on the enclosing membrane. We focus on the range of intermediate to high motility strengths and analyze the effective pressure in the steady state. Our findings for the active pressure produced in the interior and exterior regions of the inclusion indicate that the pressure is higher in the region with lower motility due to the relatively stronger accumulation of active particles.

The contribution of cross- and self-correlations to the dielectric and light-scattering spectra of supercooled polar glass formers has recently become a most challenging problem. Herein, we employ dielectric spectroscopy, depolarized dynamic light scattering (DDLS), and rheology to thoroughly examine the dynamics of van der Waals liquid 1,2-Diphenylvinylene. Carbonate (DVC), which is a polar counterpart of canonical glass former ortho-Terphenyl (OTP). We show that the light-scattering data correspond well with the dielectric permittivity function over a wide $T$ range. This pattern is very different from the peaks’ separation ${\omega }_{\max }^{\mathrm{DDLS}}/{\omega }_{\max }^{\mathrm{BDS}}=3.7$ reported recently for tributyl phosphate (TBP), despite the same dielectric characteristics of these two glass formers (${\beta }_{\mathrm{KWW}}=0.75, \mathrm{\Delta }\varepsilon =20$ for both TBP and DVC; KWW stands for Kohlrausch-Williams-Watts). This indicates different influence of orientational correlations in both methods for these two systems. We also show the results of the computer simulations of the model, polar molecules, which clearly indicate that the contribution of the cross-term to the correlation function probed in the DDLS experiment can be significant.

The energy transfer between bouncing particles and rigid boundaries during impacts is crucially influenced not only by restitution coefficients of the material but also by particle shapes. This is particularly important when such particles are mechanically agitated with vibrating plates. Inertial measurement units are able to measure all acceleration and rotational velocity components of an object and store these data for subsequent analysis. We employ them to measure the dynamics of cubes and icosahedra on vibrating plates to study the efficiency of energy transfer into the individual degrees of freedom (DOFs) of the excited object. The rotational DOFs turn out to be much less excited than the vertical translational motion. Most remarkably, there is only little difference between the two Platonic solids in both the absolute energies and the energy partition ratios.