We present a family of graphs with remarkable properties. They are obtained by connecting the points of a random walk when their distance is smaller than a given scale. Their degree (number of neighbors) does not depend on the graph's size but only on the considered scale. It follows a gamma distribution and thus presents an exponential decay. Levy flights are particular random walks with some power-law increments of infinite variance. When building the geometric graphs from them, we show from dimensional arguments that the number of connected components (clusters) follows an inverse power of the scale. The distribution of the size of their components, properly normalized, is scale invariant, which reflects the self-similar nature of the underlying process. This allows to test if a graph (including nonspatial ones) could possibly result from an underlying Levy process. When the scale increases, these graphs never tend towards a single cluster, the giant component. In other words, while the autocorrelation of the process scales as a power of the distance, they never undergo a phase transition of percolation type. The Levy graphs may find applications in community detection and in the analysis of collective behaviors as in face-to-face interaction networks.

Avalanches are often defined as signals higher than some detection level in bursty systems. The choice of the detection threshold affects the number of avalanches, but it can also affect their temporal correlations. We simulated the depinning of a long-range elastic interface and applied different thresholds including a zero one on the data to see how the sizes and durations of events change and how this affects temporal avalanche clustering. Higher thresholds result in steeper size and duration distributions and cause the avalanches to cluster temporally. Using methods from seismology, the frequency of the events in the clusters was found to decrease as a power-law of time, and the size of an event in a cluster was found to help predict how many events it is followed by. The results bring closer theoretical studies of this class of models to real experiments, but also highlight how different phenomena can be obtained from the same set of data.

Rotors of reaction and diffusion are phase singularities that give rise to spiral waves of chemical activity, which are very similar to spatiotemporal patterns observed across several excitable media. Here we carry out experiments with the Belousov-Zhabotinsky reaction system and numerical simulations based on a reaction-diffusion model to show the possible interactions of multiple spiral rotors. When the cores of two spirals come very close to each other, they could either repel, attract, or remain stationary, depending on their relative chirality, phase, and distance separating them. Multiple pairs of spiral waves, in proximity to each other, could alter the paths of the individual rotors. A spiral core will be influenced most by the rotor that is closest to it, depending on the nature of the corresponding spiral wave arm. We observed rotors lying within a limiting distance of each other attract and annihilate. Otherwise, they push each other away until they reach a critical distance, beyond which all interactions seem to cease. We have established a relationship of this critical distance to the properties of the spiral wave. We also observed spontaneous symmetry-breaking instability for a system of up to eight rotors. Our experiments with the Belousov-Zhabotinsky reaction have successfully demonstrated the validity of the numerical results. A thorough understanding of the dynamics of several spiral rotors within a small area could help us perceive the nature of such excitation waves in cardiac tissue and cell membranes.

The system of a transversely pumped Bose-Einstein condensate (BEC) coupled to a lossy ring cavity can favor a supersolid steady state. Here we find the existence of supersolid gap soliton in such a driven-dissipative system. By numerically solving the mean-field atom-cavity field coupling equations, gap solitons of a few different families have been identified. Their dynamical properties, including stability, propagation, and soliton collision, are also studied. Due to the feedback atom-intracavity field interaction, these supersolid gap solitons show numerous new features compared with the usual BEC gap solitons in static optical lattices.

We study an asymmetric version of the threshold model of binary decision making with anticonformity under an asynchronous update mode that mimics continuous time. We analyze this model on a complete graph using three different approaches: the mean-field approximation, Monte Carlo simulation, and the Markov chain approach. The latter approach yields analytical results for arbitrarily small systems, in contrast to the mean-field approach, which is strictly correct only for an infinite system. We show that, for sufficiently large systems, all three approaches produce the same results, as expected. We consider two cases: (1) homogeneous, in which all agents have the same tolerance threshold, and (2) heterogeneous, in which thresholds are given by a beta distribution parametrized by two positive shape parameters $\alpha$ and $\beta$. The heterogeneous case can be treated as a generalized model that reduces to a homogeneous model in special cases. We show that particularly interesting behaviors, including social hysteresis and critical mass reported in innovation diffusion, arise only for values of $\alpha$ and $\beta$ that yield the shape of the distribution observed in reality.

We propose a method to modulate the drifting motion of overdamped circle swimmers in steady fluid flows by means of static sinusoidal potentials. Using Langevin formalism, we study drift velocity as a function of potential strength and wavelength with and without diffusional motion. Drift velocity is essentially quantized without diffusion, but in the presence of noise, the displacement per cycle has a continuous range. As a function of dimensionless potential wave number, domains of damped oscillatory and plateau regimes are observed in the drift velocity diagram. At weak potential and fluid velocity less than powered velocity, there is also a regime where drift velocity exceeds the fluid velocity. Methods based on these results can be used to separate biological and artificial circle swimmers based on their dynamical properties.

One of the simplest models that is used to study the isotropic-nematic tranition in liquid crystalline systems is the Lebwohl-Lasher model. Several extensions of this model further enhanced its applicability. We combine two of these extensions (a generalization and the inclusion of a chiral term) and study the phase behavior and the nature of phase transitions of the resulting generalized chiral Lebwohl-Lasher model using Monte Carlo simulations. We find that the type (and even the existence) of the transition depends on the combination of the width of the interaction potential, the strength of the chiral part of the interaction, and the geometry of the system. As well, the pitch of the cholesteric bulk phase changes on approaching the phase transition if the interaction width differs from the one of the original Lebwohl-Lasher model. Thus, we identify parameter combinations that allow one to tune the properties of the ordered phase and the nature of the order-disorder transition.

We present a mathematical description of amorphous solid deformation and plasticity by extending the concept of instantaneous normal modes (INMs) to deformed systems, which allows us to retain the effect of strain on the vibrational density of states (VDOS). Starting from the nonaffine lattice dynamics (NALD) description of elasticity and viscoelasticity of glasses, we formulate the linear response theory up to large deformations by considering the strain-dependent tangent modulus at finite values of shear strain. The (nonaffine) tangent shear modulus is computed from the VDOS of affinely strained configurations at varying strain values. The affine strain, found analytically on the static (undeformed) snapshot of the glass, leads to configurations that are rich with soft low-energy modes as well as unstable modes (negative eigenvalues) that are otherwise completely “washed out” and lost if one lets the system fully relax after strain. This procedure is consistent with the structure of NALD. The INM spectrum of deformed states allows for the analytical prediction of the stress-strain curve of a model glass. Good parameter-free quantitative agreement is shown between the prediction and simulations of athermal quasistatic shear of a coarse-grained polymer glass.

Accurate calculations of the spectral density in a strongly correlated quantum many-body system are of fundamental importance to study its dynamics in the linear response regime. Typical examples are the calculation of inclusive and semiexclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter. Integral transform techniques play an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned. In the present work we extend a recently proposed quantum algorithm which circumvents this problem. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigorous error estimates. An appropriate expansion in Chebyshev polynomials allows for efficient simulations also on classical computers. We apply our idea to obtain the local density of states for graphene in a magnetic field as a proof of principle. This paves the way for future applications in nuclear and condensed matter physics.

We develop a unified Hamiltonian approach to the diffusion of a particle coupled to a dissipative environment, an archetypal model widely invoked to interpret condensed phase phenomena, such as polymerization and cold-atom diffusion in optical lattices. By appropriate choices of the coupling functions, we reformulate phenomenological diffusion models by adding otherwise ignored space-momentum terms. We thus numerically predict a variety of diffusion regimes, from diffusion saturation to superballistic diffusion. With reference to ultracold atoms in optical lattices, we also show that time correlated external noises prevent superdiffusion from exceeding Richardson's law. Some of these results are unexpected and call for experimental validation.

Scattering of an acoustic wave by particles gives rise to microstreaming, as well as to acoustic radiation and interaction forces on the particles. We numerically study these steady, nonlinear phenomena for a case of two elastic spheres in a standing wave. We show that if one or both spheres are smaller or comparable to the viscous boundary layer, the microstreaming close to the pressure node can lead to an interparticle attraction along the direction of the pressure gradient of the wave. Similar behavior is observed when, instead of size, density of one of the spheres is sufficiently larger relative to the other sphere. These findings could promote the acoustic manipulation of nanoparticles and bacteria.

Higher implosion efficiency is of great significance in direct-drive fusion research. We demonstrated the critical role played by the intensity gradient of the main drive laser pulse in improving efficiency of direct-drive implosions, using a double-gradient nanosecond pulse. Compared with a square pulse, the burn-through depth was increased by over 200%, and the shell velocity was increased by $\sim 2.1$ times with an optimized double-gradient pulse. As the result, the implosion efficiency was enhanced by $\sim$ six times. It was found that by limiting the intensity gradient of the main drive pulse to no more than $\sim 2.5\times {10}^{15}\mathrm{W}/\left(\mathrm{c}{\mathrm{m}}^2\mathrm{ns}\right)$, heat flux inhibition by nonlocal electron thermal transport effects could be eliminated, and ultimately an efficient mass ablation process was achieved. These results have relevance for pulse designs in ignition-scale direct-drive implosions.

We propose a model for a chain of particles coupled by nonlinear springs in which each mass has an internal mass and all interactions are assumed to be nonlinear. We show how to construct an asymptotic solution of this system using multiple timescales, the systematic solution of coupled equations by repeated application of a consistency condition. Our results show that for some combinations of nonlinearity the dynamics are governed by the NLS as in the more usual mass-in-mass chains with linear interactions between inner and outer masses. However, when both nonlinearities have quadratic components, we show that the asymptotic reduction results in a Ginzburg-Landau equation instead of NLS.

Alcohols are organic compounds characterized by one or more hydroxyl groups attached to a carbon atom of an alkyl group. They can be considered as organic derivatives of water in which one of the hydrogen atoms is replaced by an alkyl group. In this work, the Mercedes-Benz model of water is used to design simple two-dimensional (2D) models of lower alcohols. The structural and thermodynamic properties of the constructed simple models are studied by conducting Monte Carlo simulations in the isothermal-isobaric ensemble. We show that 2D models display similar trends in structuring and thermodynamics as in experiments. The present work on the smallest amphiphilc organic solutes provides a simple testing ground to study the competition between polar and non-polar effects within the molecule and physical properties.

Glycerol is one of the glass-forming liquids selected by Robert H. Cole in 1950 to start his study of molecular dynamics by dielectric spectroscopy. Seventy-one years have gone by and remarkably no consensus has been reached on the nature and identity of the relaxation processes observed in the dielectric spectra. The macroscopic dielectric relaxation data allow different interpretations to yield contrasting results, and it is not possible to determine which one is most plausible. Coming to the rescue is the application of the nuclear $\gamma$-resonance time-domain interferometry (TDI) to glycerol by Saito et al. [Phys. Rev. E 105, L012605 (2022)]. Their microscopic TDI data potentially can decide which interpretation of the dielectric spectra of glycerol is most plausible. The attempt was made by Saito et al., but there is a problem in their analysis of the dielectric data of glycerol and hence their conclusion is untenable. In this paper, we critically compare four major interpretations with the TDI data in an effort to identify the most plausible interpretation of the relaxation processes constituting the dielectric spectra of glycerol.

We derive an exact law for compressible pressure-anisotropic magnetohydrodynamic turbulence. For a gyrotropic pressure tensor, we study the double-adiabatic case and show the presence of new flux and source terms in the exact law, reminiscent of the plasma instability conditions due to pressure anisotropy. The Hall term is shown to bring ion-scale corrections to the exact law without affecting explicitly the pressure terms. In the pressure isotropy limit we recover all known results obtained for isothermal and polytropic closures. The incompressible limit of the gyrotropic system leads to a generalization of the Politano and Pouquet's law where a new incompressible source term is revealed and reflects exchanges of the magnetic and kinetic energies with the no-longer-conserved internal energy. We highlight the possibilities offered by the new laws to investigate potential links between turbulence cascade and instabilities widely observed in laboratory and astrophysical plasmas.

We consider a simple model for multidimensional conewise linear dynamics around cusplike equilibria. We assume that the local linear evolution is either ${\mathbf{v}}^{\prime }=\mathbb{A}\mathbf{v}$ or $\mathbb{B}\mathbf{v}$ (with $\mathbb{A}, \mathbb{B}$ independently drawn from a rotationally invariant ensemble of symmetric $N\times N$ matrices) depending on the sign of the first component of $\mathbf{v}$. We establish strong connections with the random diffusion persistence problem. When $N\to \infty$, we find that the Lyapunov exponent is non-self-averaging, i.e., one can observe apparent stability and apparent instability for the same system, depending on time and initial conditions. Finite $N$ effects are also discussed and lead to cone trapping phenomena.

In the thermodynamics of nanoscopic systems the relation between classical and quantum mechanical description is of particular importance. To scrutinize this correspondence we have chosen two two-dimensional billiard systems. Both systems are studied in the classical and the quantum mechanical settings. The classical conditional probability density $p(E,L|E_0,L_0)$ as well as the quantum mechanical transition probability $P(n,l|n_0,l_0)$ are calculated, which build the basis for the statistical analysis. We calculate the work distribution for one particle. The results in the quantum case in particular are of special interest since a suitable definition of mechanical work in small quantum systems is already controversial. Furthermore, we analyze the probability of both zero work and zero angular momentum difference. Using connections to an exactly solvable system analytical formulas are given for both systems. In the quantum case we get numerical results with some interesting relations to the classical case.