A discrete-velocity Boltzmann model for a nine-velocity lattice is considered. Compared to the conventional lattice Boltzmann (LB) schemes the collisions for the model are defined explicitly. Space and time discretization of the model is based on the collide and stream method; in addition, the regularization of the collision term is proposed. It is demonstrated that the regularized model can be represented as a two-relaxation-time LB model of a special type. The scheme is compared to the Onsager regularized (a specific filtered Galilean invariant model) and recursively regularized LB equations in terms of stability and dissipation properties, and linear stability analysis is performed. Several numerical experiments are carried out: double shear layer, lid-driven cavity flow, and propagation of acoustic and shear waves are considered for different grid resolutions, Mach and Reynolds numbers. It is shown that free parameters in the model corresponding to collision cross sections can be adjusted in such a way that the dissipation and stability properties can be controlled.

In quantum Monte Carlo (QMC) methods, energy estimators are calculated as (functions of) statistical averages of quantities sampled during a calculation. Associated statistical errors of these averages are often estimated. This error estimation is not straightforward and there are several choices of the error estimation methods. We evaluate the performance of three methods (the Straatsma method, an autoregressive model, and a blocking analysis based on von Neumann's ratio test for randomness) for the energy time series given by three QMC methods [diffusion Monte Carlo, full configuration interaction Quantum Monte Carlo (FCIQMC), and coupled cluster Monte Carlo (CCMC)]. From these analyses, we describe a hybrid analysis method which provides reliable error estimates for a series of various lengths of FCIQMC and CCMC's time series. Equally important is the estimation of the appropriate start point of the equilibrated phase. We establish that a simple mean squared error rule method as described by White [K. P. White, Jr., Simulation 69(6), 323 (1997)] can provide reasonable estimations.

We explore the charging of a quantum battery based on spin systems through periodic modulation of a transverse-field-like Ising Hamiltonian. In the integrable limit, we find that resonance tunneling can lead to a higher transfer of energy to the battery and better stability of the stored energy at specific drive frequencies. When the integrability is broken in the presence of an additional longitudinal field, we find that the effective Floquet Hamiltonian contains terms which may lead to a global charging of the battery. However, we do not find any quantum advantage in the charging power, thus demonstrating that global charging is only a necessary and not sufficient condition for achieving quantum advantage.

Localized perturbations in a real-world network have the potential to trigger cascade failures at the whole system level, hindering its operations and functions. Standard approaches analytically tackling this problem are mostly based either on static descriptions, such as percolation, or on models where the failure evolves through first-neighbor connections, crucially failing to capture the nonlocal behavior typical of real cascades. We introduce a dynamical model that maps the failure propagation across the network to a self-avoiding random walk that, at each step, has a probability to perform nonlocal jumps toward operational systems' units. Despite the inherent non-Markovian nature of the process, we are able to characterize the critical behavior of the system out of equilibrium, as well as the stopping time distribution of the cascades. Our numerical experiments on synthetic and empirical biological and transportation networks are in excellent agreement with theoretical expectation, demonstrating the ability of our framework to quantify the vulnerability to nonlocal cascade failures of complex systems with interconnected structure.

Interaction equations governing slow time evolution of the coordinates and phases of two interacting temporal cavity solitons in a delay differential equation model of a nonlinear mirror mode-locked laser are derived and analyzed. It is shown that long-range soliton interaction due to gain depletion and recovery can lead either to a development of a harmonic mode-locking regime or to a formation of closely packed incoherent soliton bound state with weakly oscillating intersoliton time separation. Short-range soliton interaction via electric field tails can result in an antiphase or in-phase stationary and breathing harmonic mode-locking regimes.

Convection is a well-studied topic in fluid dynamics, yet it is less understood in the context of network flows. Here, we incorporate techniques from topological data analysis (namely, persistent homology) to automate the detection and characterization of convective flows (also called cyclic or chiral flows) over networks, particularly those that arise for irreversible Markov chains. As two applications, we study convection cycles arising under the PageRank algorithm and we investigate chiral edge flows for a stochastic model of a bimonomer's configuration dynamics. Our experiments highlight how system parameters—e.g., the teleportation rate for PageRank and the transition rates of external and internal state changes for a monomer—can act as homology regularizers of convection, which we summarize with persistence barcodes and homological bifurcation diagrams. Our approach establishes a connection between the study of convection cycles and homology, the branch of mathematics that formally studies cycles, which has diverse potential applications throughout the sciences and engineering.

We show that two active particles can form a bound state by coupling to a driven nonequilibrium environment. We specifically investigate the case of two mutually noninteracting run-and-tumble probes moving on a ring, each in short-range interaction with driven colloids. Under conditions of timescale separation, these active probes become trapped in bound states. In fact, the bound state appears at high enough persistence (low effective temperature). From the perspective of a comoving frame, where colloids are in thermal equilibrium and the probes are active and driven, an appealing analogy appears with Cooper pairing, as electrons can be viewed as run-and-tumble particles in a pilot-wave picture.

Electrical image potentials can be important in small spaces, such as nanoscale pores in porous electrodes, which are used in capacitive desalination and in supercapacitors, as argued by Bazant's group at Massachusetts Institute of Technology. It will be shown here that inside pores in porous metallic materials the image potentials can be considerably larger than near flat walls, as a result of the fact that the dielectric constant for an electric field perpendicular to a wall is much smaller than the bulk dielectric constant of water. Calculations will be presented for the image potential in spherical and cylindrically shaped pores. The calculations for cylindrical pores can also be applied to nanotubes. It was believed for a long time, on the basis of molecular dynamics simulations, that in order to push a salt solution through a small radius nanotube, work must be done against the solvation energy of the ions, which is larger inside a narrow nanotube than it is in the bulk. The relatively large image charge potential energy in narrow nanotubes, however, tends to oppose this increase in the solvation energy. The degree to which the image potential facilitates the flow of the salt ions into nanotubes will be discussed.

Intermittent motion, called stick-slip, is a friction instability that commonly occurs during relative sliding of two elastic solids. In adhesive polymer contacts, where elasticity and interface adhesion are strongly coupled, stick-slip arises due to the propagation of slow detachment waves at the interface. Here we analyze two distinct detachment waves moving parallel (Schallamach wave) and antiparallel (separation wave) to applied remote sliding. Both waves cause slip in the same direction, travel at speeds much lesser than any elastic wave speed, and are therefore describable using the same perturbative elastodynamic framework with identical boundary conditions. A numerical scheme is used to obtain interface stresses and velocities for arbitrary Poisson ratio, along with closed-form solutions for incompressible solids. Our calculations reveal a close correspondence between moving detachment waves and bimaterial interface cracks, including the nature of the singularity and the functional forms of the stresses. Based on this correspondence, and coupled with a fracture analogy for dynamic friction, we develop a phase diagram showing domains of possible occurrence of stick-slip via detachment waves vis-á-vis steady interface sliding. Our results have interesting implications for sliding and stick-slip phenomena at soft interfaces.

It is fascinating that fish groups spontaneously form different formations. The collective locomotions of two and multiple undulatory self-propelled foils swimming in a fluid are numerically studied and the deep reinforcement learning (DRL) is applied to control the locomotion. We explored whether typical patterns emerge spontaneously under the driven two DRL strategies. One strategy is that only the following fish gets hydrodynamic advantages. The other is that all individuals in the group take advantage of the interaction. In the DRL strategy, we use swimming efficiency as the reward function, and the visual information is included. We also investigated the effect of involving hydrodynamic force information, which is an analogy to that detected by the lateral line of fish. Each fish can adjust its undulatory phase to achieve the goal. Under the two strategies, collective patterns with different characteristics, i.e., the staggered-following, tandem-following phalanx and compact modes emerge. They are consistent with the results in the literature. The hydrodynamic mechanism of the above high-efficiency collective traveling modes is analyzed by the vortex-body interaction and thrust. We also found that the time sequence feature and hydrodynamic information in the DRL are essential to improve the performance of collective swimming. Our research can reasonably explain the controversial issue observed in the relevant experiments. The paper may be helpful for the design of bionic fish.

Resonant four-photon scattering could nearly double frequencies of intense laser pulses in plasma. However, transverse slippage between pulses presents a technological challenge, while collinear four-photon scattering is forbidden for classical light dispersion in plasma. Nonlinear renormalization of intense laser pulses can enable collinear four-photon resonance. However, such a very intensity-sensitive resonance is difficult to maintain for evolving pulses. Remarkably, there is a lower-dimensionality submanifold of the resonant four-photon manifold where the evolving pulses stay in resonance. This could enable an all-optical frequency doubling of mildly relativistic-intense laser pulses in collinear geometry, advantageously free of the transverse slippage challenges.

The main motivation of this paper is to introduce the permutation Jensen-Shannon distance, a symbolic tool able to quantify the degree of similarity between two arbitrary time series. This quantifier results from the fusion of two concepts, the Jensen-Shannon divergence and the encoding scheme based on the sequential ordering of the elements in the data series. The versatility and robustness of this ordinal symbolic distance for characterizing and discriminating different dynamics are illustrated through several numerical and experimental applications. Results obtained allow us to be optimistic about its usefulness in the field of complex time-series analysis. Moreover, thanks to its simplicity, low computational cost, wide applicability, and less susceptibility to outliers and artifacts, this ordinal measure can efficiently handle large amounts of data and help to tackle the current big data challenges.

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. We perform a variety of tests, and in all cases find rapid convergence of the chemical potential—inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

We study numerically the critical behavior and marginal stability of the shear jamming transition for frictionless soft spheres, observed to occur over a finite range of densities, associated with isotropic jamming for densities above the minimum jamming ($J$-point) density. Several quantities are shown to scale near the shear jamming point in the same way as the isotropic jamming point. We compute the exponents associated with the small force distribution and the interparticle gap distribution and show that the corresponding exponents are consistent with the marginal stability condition observed for isotropic jamming and with predictions of the mean-field theory of jamming in hard spheres.

We study in detail an open quantum generalization of a classical kinetically constrained model—the East model—known to exhibit slow glassy dynamics stemming from a complex hierarchy of metastable states with distinct lifetimes. Using the recently introduced theory of classical metastability for open quantum systems, we show that the driven open quantum East model features a hierarchy of classical metastabilities at low temperature and weak driving field. We find that the effective long-time description of its dynamics not only is classical, but shares many properties with the classical East model, such as obeying an effective detailed balance condition and lacking static interactions between excitations, but with this occurring within a modified set of metastable phases which are coherent, and with an effective temperature that is dependent on the coherent drive.

We present a Monte Carlo study of the backbone and the shortest-path exponents of the two-dimensional $Q$-state Potts model in the Fortuin-Kasteleyn bond representation. We first use cluster algorithms to simulate the critical Potts model on the square lattice and obtain the backbone exponents $d_B=1.7320\left(3\right)$ and 1.794(2) for $Q=2,3$, respectively. However, for large $Q$, the study suffers from serious critical slowing down and slowly converging finite-size corrections. To overcome these difficulties, we consider the $\mathrm{O}\left(n\right)$ loop model on the honeycomb lattice in the densely packed phase, which is regarded to correspond to the critical Potts model with $Q=n^2$. With a highly efficient cluster algorithm, we determine from domains enclosed by the loops $d_B=1.64339\left(5\right),1.73227\left(8\right),1.7938\left(3\right),1.8384\left(5\right),1.8753\left(6\right)$ for $Q=1,2,3,2+\sqrt{3},4$, respectively, and $d_{\min }=1.0945\left(2\right),1.0675\left(3\right),1.0475\left(3\right),1.0322\left(4\right)$ for $Q=2,3,2+\sqrt{3},4$, respectively. Our estimates significantly improve over the existing results for both $d_B$ and $d_{\min }$. Finally, by studying finite-size corrections in backbone-related quantities, we conjecture an exact formula as a function of $n$ for the leading correction exponent.

We study the underdamped motion of a passive particle in an active environment. Using the phase space path integral method we find the probability distribution function of position and velocity for a free and a harmonically bound particle. The environment is characterized by an active noise which is described as the Ornstein-Uhlenbeck process (OUP). Taking two similar, yet slightly different OUP models, it is shown how inertia along with other relevant parameters affect the dynamics of the particle. Further we investigate the work fluctuations of a harmonically trapped particle by considering the trap center being pulled at a constant speed. Finally, the fluctuation theorem of work is validated with an effective temperature in the steady-state limit.

We study a system of self-propelled particles (SPPs) in which individual particles are allowed to switch between a fast aligning and a slow nonaligning state depending upon the degree of the alignment in the neighborhood. The switching is modeled using a threshold for the local order parameter. This additional attribute gives rise to a mixed phase, in contrast to the ordered phases found in clean SPP systems. As the threshold is increased from zero, we find the sudden appearance of clusters of nonaligners. Clusters of nonaligners coexist with moving clusters of aligners with continual coalescence and fragmentation. The behavior of the system with respect to the clustering of nonaligners appears to be very different for values of low and high global densities. In the low density regime, for an optimal value of the threshold, the largest cluster of nonaligners grows in size up to a maximum that varies logarithmically with the total number of particles. However, on further increasing the threshold the size decreases. In contrast, for the high density regime, an initial abrupt rise is followed by the appearance of a giant cluster of nonaligners. The latter growth can be characterized as a continuous percolation transition. In addition, we find that the speed differences between aligners and nonaligners is necessary for the segregation of aligners and nonaligners.