Long-range hoppings in quantum disordered systems are known to yield quantum multifractality, the features of which can go beyond the characteristic properties associated with an Anderson transition. Indeed, critical dynamics of long-range quantum systems can exhibit anomalous dynamical behaviors distinct from those at the Anderson transition in finite dimensions. In this paper, we propose a phenomenological model of wave packet expansion in long-range hopping systems. We consider both their multifractal properties and the algebraic fat tails induced by the long-range hoppings. Using this model, we analytically derive the dynamics of moments and inverse participation ratios of the time-evolving wave packets, in connection with the multifractal dimension of the system. To validate our predictions, we perform numerical simulations of a Floquet model that is analogous to the power law random banded matrix ensemble. Unlike the Anderson transition in finite dimensions, the dynamics of such systems cannot be adequately described by a single parameter scaling law that solely depends on time. Instead, it becomes crucial to establish scaling laws involving both the finite size and the time. Explicit scaling laws for the observables under consideration are presented. Our findings are of considerable interest towards applications in the fields of many-body localization and Anderson localization on random graphs, where long-range effects arise due to the inherent topology of the Hilbert space.

In this work, a stochastic model of gaseous transfer in polymer–carbon-nanotube (CNT) nanocomposites is presented. The model takes into account interfacial areas, i.e., polymer depletion regions. The local regime of transport is controlled by the density of the polymer. In a dense polymer, this regime corresponds to the ordinary diffusion, while in free volume regions, it corresponds to the ballistic transport. The introduction of a free volume and/or a depleted polymer layer near to a CNT wall leads to the emergence of anomalous diffusion. We have demonstrated how the anomalous diffusion regime changes in the presence of nanotubes for different distributions of polymer density. The presented approach allows us to describe the threshold effect in the diffusion coefficient as a function of CNTs density in polymer-CNT nanocomposites.

The nonlinear supratransmission phenomenon in a Fermi-Pasta-Ulam (FPU) diatomic lattice with two forbidden bands is investigated. Using a decoupling ansatz for the motion of the two different sublattices combined with the continuum (quasidiscrete) approximation, we derived analytically the threshold amplitudes of supratransmission occurrence when a sinusoidal driving with frequency in the upper forbidden band (lower forbidden band gap between acoustic and optical modes) is applied at one end. The resulting estimate of the threshold of a lattice with a first heavy particle is different to the one obtained from a lattice with a first light particle, showing the influence of the driven particle and giving also the possibility to have two thresholds on each forbidden gap of a diatomic lattice by switching the order of light ($m$) and heavy ($M$) masses. In the lower forbidden band, the dependence of the supratransmission threshold on the mass ratio ($\mu =m/M$) has been evidenced and it appears that for large (small) values of $\mu$, that is $\mu >60\%$, the coupling between the two modes must (must not) be considered. Numerical explorations were subsequently performed with an emphasis on the dependence of the threshold on the driving frequency and also on the mass of the driven particle (light or heavy). A good agreement is found between the numerical and analytical thresholds. For the limit case where all the masses are identical, the results of the monoatomic FPU previously found in the literature are recovered.

The properties of mixed eigenstates in a generic quantum system with a classical counterpart that has mixed-type phase space, although important to understand several fundamental questions that arise in both theoretical and experimental studies, are still not clear. Here, following a recent work [Č. Lozej, D. Lukman, and M. Robnik, Phys. Rev. E 106, 054203 (2022)], we perform an analysis of the features of mixed eigenstates in a time-dependent Hamiltonian system, the celebrated kicked top model. As a paradigmatic model for studying quantum chaos, the kicked top model is known to exhibit both classical and quantum chaos. The types of eigenstates are identified by means of the phase-space overlap index, which is defined as the overlap of the Husimi function with regular and chaotic regions in classical phase space. We show that the mixed eigenstates appear due to various tunneling precesses between different phase-space structures, while the regular and chaotic eigenstates are, respectively, associated with invariant tori and chaotic components in phase space. We examine how the probability distribution of the phase-space overlap index evolves with increasing system size for different kicking strengths. In particular, we find that the relative fraction of mixed states exhibits a power-law decay as the system size increases, indicating that only purely regular and chaotic eigenstates are left in the strict semiclassical limit. We thus provide further verification of the principle of uniform semiclassical condensation of Husimi functions and confirm the correctness of the Berry-Robnik picture.

Symmetries in a network regulate its organization into functional clustered states. Given a generic ensemble of nodes and a desirable cluster (or group of clusters), we exploit the direct connection between the elements of the eigenvector centrality and the graph symmetries to generate a network equipped with the desired cluster(s), with such a synthetical structure being furthermore perfectly reflected in the modular organization of the network's functioning. Our results solve a relevant problem of designing a desired set of clusters and are of generic application in all cases where a desired parallel functioning needs to be blueprinted.

Light sources that switch periodically over time have a wide range of application value in life and engineering, and generally require additional controller to periodically switch circuits to achieve periodic lighting. In this paper, a self-oscillating spring oscillator based on optically responsive liquid crystal elastomer (LCE) fiber is constructed, which consists of a embedded light source and a LCE fiber. The spring oscillator can oscillate autonomously to achieve periodic switching of the light source. On the basis of the well-established dynamic LCE model, a nonlinear dynamic model is proposed and its dynamic behavior is studied. Numerical calculations demonstrate that the spring oscillator presents two motion regimes, namely the self-oscillation regime and the static regime. The self-oscillation of spring oscillator is maintained by the energy competition between light energy and damping dissipation. Furthermore, the critical conditions for triggering self-oscillation are also investigated in detail, as well as the key system parameters that affect its frequency and amplitude. Different from the existing abundant self-oscillating systems, this self-oscillating structure with simple structure and convenient fabrication does not require complex controller to obtain periodic lighting, and it is expected to provide more diversified design ideas for soft robots and sensors.

The theoretical background of the nonperturbative method of spectral response to stochastic processes (SRSP) for measuring the nonreciprocal interparticle effective interactions in strongly coupled underdamped systems is described. Analytical expressions for vibrational spectral density of confined Brownian particles with a nonreciprocal effective interaction are presented. The changes in the vibrational spectral density with varying different parameters of the system (nonreciprocity, viscosity, ratios of particle sizes, and intensities of random processes acting on each particle) are discussed using the example of a pair of nonidentical particles in a harmonic trap. The SRSP method is compared to three other nonperturbative methods. The SRSP method demonstrates an undeniable advantage when processing particle trajectories with errors in particle tracking.

This paper proposes a simple and accurate lattice Boltzmann model for simulating thermocapillary flows, which can deal with the contrast between thermodynamic parameters. In this model, two lattice Boltzmann equations are utilized to solve the conservative Allen-Cahn equation and the incompressible Navier-Stokes equations, while another lattice Boltzmann equation is used for solving the temperature field, where the collision term is delicately designed such that the influence of the contrast between thermodynamic parameters is incorporated. In contrast to the previous lattice Boltzmann models for thermocapillary flows, the most distinct feature of the current model is that the forcing term used in the present thermal lattice Boltzmann equation is not needed to calculate space derivatives of the heat capacitance or the order parameter, making the scheme much more straightforward and able to retain the main merits of the lattice Boltzmann method. The developed model is first validated by considering the thermocapillary flows in a heated microchannel with two superimposed planar fluids. It is then used to simulate the thermocapillary migration of a two-dimensional deformable droplet, and its accuracy is consistent with the theoretical prediction when the Marangoni number approaches zero. Finally, we numerically study the motion of two recalcitrant bubbles in a two-dimensional channel where the relationship between surface tension and temperature is assumed to be a parabolic function. It is observed that due to the competition between the inertia and thermal effects, the bubbles can move against the liquid's bulk motion and towards areas with low surface tension.

We present an explicit finite-difference method to simulate the nonideal multiphase fluid flow. The local density and momentum transport are modeled by the Navier-Stokes equations and the pressure is computed by the van der Waals equation of the state. The static droplet and the dynamics of liquid-vapor separation simulations are performed as validations of this numerical scheme. In particular, to maintain the thermodynamic consistency, we propose a general wetting energy boundary condition at the contact line between fluids and the solid boundary. We conduct a series of comparisons between the current boundary condition and the constant contact angle boundary condition as well as the stress-balanced boundary condition. This boundary condition alleviates the instability induced by the constant contact angle boundary condition at $\theta \approx 0$ and $\theta \approx \pi$. Using this boundary condition, the equilibrium contact angle is correctly recovered and the contact line dynamics are consistent with the simulation by applying a stress-balanced boundary condition. Nevertheless, unlike the stress-balanced boundary condition for which we need to further introduce the interface thickness parameter, the current boundary condition implicitly incorporates the interface thickness information into the wetting energy.

Phase transitions of the $J_1\text{−}{J}_2$ Ising model on a square lattice are studied using the higher-order tensor renormalization group (HOTRG) method. This system involves a competition between the ferromagnetic interaction $J_1$ and antiferromagnetic interaction $J_2$, and in previous studies, weak first-order and second-order transitions were observed near the ratio $g=J_2/\left|J_1\right|=1/2$. It has also been suggested that the universality class of the second-order phase transition connected to the first-order transition line for $g>1/2$ belongs to the Ashkin-Teller class, which is characterized by a continuously varying critical exponent with $g$, as predicted by field-theoretical and other studies. Our results, based on the HOTRG calculations for significantly larger sizes, indicate that the region of the first-order transition is marginally narrower than that in previous studies. Furthermore, it is suggested that the region where the critical exponent changes does not necessarily coincide with the Ashkin-Teller region.

In quasi-two-dimensional experiments with photoelastic particles confined to an annular region, an intruder constrained to move in a circular path halfway between the annular walls experiences stick-slip dynamics. We discuss the response of the granular medium to the driven intruder, focusing on the evolution of the force network during sticking periods. Because the available experimental data do not include precise information about individual contact forces, we use an approach developed in our previous work [Basak et al., J. Eng. Mech. 147, 04021100 (2021)] based on networks constructed from measurements of the integrated strain magnitude on each particle. These networks are analyzed using topological measures based on persistence diagrams, revealing that force networks evolve smoothly but in a nontrivial manner throughout each sticking period, even though the intruder and granular particles are stationary. Characteristic features of persistence diagrams show identifiable slip precursors. In particular, the number of generators describing the structure and complexity of force networks increases consistently before slips. Key features of the dynamics are similar for granular materials composed of disks or pentagons, but some details are consistently different. In particular, we find significantly larger fluctuations of the measures computed based on persistence diagrams and, therefore, of the underlying networks, for systems of pentagonal particles.

We present a scheme for calculating the shape of two well-known conical motifs: the d-Cone and the e-Cone. Each begins as a thin, flat disk, before buckling during loading into a deformed shape with distinctive, asymmetrical conical features and a localised apex. Various deformed equilibrium models rightly assume a developable shape, with a particular focus on determining how much of the disk detaches from how it is supported during buckling; they are, nevertheless, extensively curated analytically, and must confront (some, ingeniously) the question of singular, viz., infinite properties at the conical apex. In this study, we find an approximate description of shape that reveals the extent of detachment, from an analogous mobile vertex that packages optimally according to its constraints. To this end, we further develop the usage of Gauss's Mapping and the associated spherical image, which has been used previously, but only to confirm known properties of deformed shape. Despite the simplicity of our approach, remarkably good predictions are availed, perhaps because such problems of extreme deformation are geometrically (rather than equilibrium) dominated.

Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k\text{−}2L\text{−}a\text{−}\mathcal{C}$ and $k\text{−}\varphi \text{−}L\text{−}a\text{−}\mathcal{C}$ models by Morgan [Phys. Rev. E 103, 053108 (2021); Phys. Rev. E 105, 045104 (2022)]. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment [J. Fluids Eng. 136, 091212 (2014)] and in simulation of a shock-accelerated localized patch of turbulence. It is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.

We provide an algorithm, i-SPin 2, for evolving general spin-$s$ Gross-Pitaevskii or nonlinear Schrödinger systems carrying a variety of interactions, where the $2s+1$ components of the “spinor” field represent the different spin-multiplicity states. We consider many nonrelativistic interactions up to quartic order in the Schrödinger field (both short and long range, and spin-dependent and spin-independent interactions), including explicit spin-orbit couplings. The algorithm allows for spatially varying external and/or self-generated vector potentials that couple to the spin density of the field. Our work can be used for scenarios ranging from laboratory systems such as spinor Bose-Einstein condensates (BECs), to cosmological or astrophysical systems such as self-interacting bosonic dark matter. As examples, we provide results for two different setups of spin-1 BECs that employ a varying magnetic field and spin-orbit coupling, respectively, and also collisions of spin-1 solitons in dark matter. Our symplectic algorithm is second-order accurate in time, and is extensible to the known higher-order-accurate methods.

We deduce a thermodynamically consistent diffuse interface model to study the line tension phenomenon of sessile droplets. By extending the standard Cahn-Hilliard model via modifying the free energy functional due to the spatial reflection asymmetry at the substrate, we provide an alternative interpretation for the wall energy. In particular, we find the connection of the line tension effect with the droplet-matrix-substrate triple interactions. This finding reveals that the apparent contact angle deviating from Young's law is contributed by the wall energy reduction as well as the line energy minimization. Besides, the intrinsic negative line tension resulting from the curvature effect is observed in our simulations and shows good accordance with recent experiments [Tan et al. Phys. Rev. Lett. 130, 064003 (2023)]. Moreover, our model sheds light upon the understanding of the wetting edge formation which results from the vying effect of wall energy and line tension.

We propose and study a one-dimensional (1D) model consisting of two lanes with open boundaries. One of the lanes executes diffusive and the other lane driven unidirectional or asymmetric exclusion dynamics, which are mutually coupled through particle exchanges in the bulk. We elucidate the generic nonuniform steady states in this model. We show that in a parameter regime, where hopping along the TASEP lane, diffusion along the SEP lane, and the exchange of particles between the TASEP and SEP lanes compete, the SEP diffusivity $D$ appears as a tuning parameter for both the SEP and TASEP densities for a given exchange rate in the nonequilibrium steady states of this model. Indeed, $D$ can be tuned to achieve phase coexistence in the asymmetric exclusion dynamics together with spatially smoothly varying density in the diffusive dynamics in the steady state. We obtain phase diagrams of the model using mean field theories, and corroborate and complement the results with stochastic Monte Carlo simulations. This model reduces to an isolated open totally asymmetric exclusion process (TASEP) and an open TASEP with bulk particle nonconserving Langmuir kinetics (LK), respectively, in the limits of vanishing and diverging particle diffusivity in the lane executing diffusive dynamics. Thus, this model works as an overarching general model, connecting both pure TASEPs and TASEPs with LK in different asymptotic limits. We further define phases in the SEP and obtain phase diagrams and show their correspondence with the TASEP phases. In addition to its significance as a 1D driven, diffusive model, this model also serves as a simple reduced model for cell biological transport by molecular motors undergoing diffusive and directed motion inside eukaryotic cells.

Using the theory of large deviations, macroscopic fluctuation theory provides a framework to understand the behavior of nonequilibrium dynamics and steady states in diffusive systems. We extend this framework to a minimal model of a nonequilibrium nondiffusive system, specifically an open linear network on a finite graph. We explicitly calculate the dissipative bulk and boundary forces that drive the system towards the steady state, and the nondissipative bulk and boundary forces that drive the system in orbits around the steady state. Using the fact that these forces are orthogonal in a certain sense, we provide a decomposition of the large-deviation cost into dissipative and nondissipative terms. We establish that the purely nondissipative force turns the dynamics into a Hamiltonian system. These theoretical findings are illustrated by numerical examples.

We consider the nonlinear Schrödinger equation with nonlocal derivatives in a two-dimensional periodic domain. For certain orders of derivatives, we find a type of quasi-breather solution dominating the field evolution at low nonlinearity levels. With the increase of nonlinearity, the structures break down, giving way to Rayleigh-Jeans (or wave turbulence) spectra. Phase-space trajectories associated with the quasibreather solutions are found to be close to that of the linear system and almost periodic. We employ two methods to search for nearby periodic solutions (e.g., exact breathers), yet none are found. Given these distinguishing behaviors, we interpret this structure in the context of Kolmogorov-Arnold-Moser (KAM) theory.

Complex biological processes involve collective behavior of entities (bacteria, cells, animals) over many length and time scales and can be described by discrete models that track individuals or by continuum models involving densities and fields. We consider hybrid stochastic agent-based models of branching morphogenesis and angiogenesis (new blood vessel creation from preexisting vasculature), which treat cells as individuals that are guided by underlying continuous chemical and/or mechanical fields. In these descriptions, leader (tip) cells emerge from existing branches and follower (stalk) cells build the new sprout in their wake. Vessel branching and fusion (anastomosis) occur as a result of tip and stalk cell dynamics. Coarse graining these hybrid models in appropriate limits produces continuum partial differential equations (PDEs) for endothelial cell densities that are more analytically tractable. While these models differ in nonlinearity, they produce similar equations at leading order when chemotaxis is dominant. We analyze this leading order system in a simple quasi-one-dimensional geometry and show that the numerical solution of the leading order PDE is well described by a soliton wave that evolves from vessel to source. This wave is an attractor for intermediate times until it arrives at the hypoxic region releasing the growth factor. The mathematical techniques used here thus identify common features of discrete and continuum approaches and provide insight into general biological mechanisms governing their collective dynamics.