We use numerical linked cluster expansions (NLCEs) to study the site-diluted transverse-field Ising model on the square lattice at $T=0$. NLCE with a self-consistent mean field on the boundary of the clusters is used to obtain the ground-state magnetization, susceptibility, and structure factor as a function of transverse field $h$ and exchange constant $J$. Adding site dilution to the model turns NLCE into a series expansion in the dilution parameter $p$. Studying the divergence of the structure factor allows us to establish the phase diagram in the $h/J$ and $p$ plane. By studying the magnetization of the system in a longitudinal field, we investigate the Griffiths-McCoy singularities. We find that the magnetization develops nonlinearities in the Griffiths phase with exponents that vary continuously with $h$. Additionally, the probability distribution of the local susceptibility develops long tails in the Griffiths phase, which is studied in terms of its moments.

The Ornstein-Zernike (OZ) integral equation theory is a powerful approach to simple liquids due to its low computational cost and the fact that, when combined with an appropriate closure equation, the theory is thermodynamically complete. However, approximate closures proposed to date exhibit pressure or free energy inconsistencies that produce inaccurate or ambiguous results, limiting the usefulness of the Ornstein-Zernike approach. To address this problem, we combine methods to enforce both pressure and free energy consistency to create a new closure approximation and test it for a single-component Lennard-Jones fluid. The closure is a simple power series in the direct and total correlation functions for which we have derived analytical formulas for the excess Helmholtz free energy and chemical potential. These expressions contain a partial molar volumelike term, similar to excess chemical potential correction terms recently developed. Using our bridge approximation, we have calculated the pressure, Helmholtz free energy, and chemical potential for the Lennard-Jones fluid using the Kirkwood charging, thermodynamic integration techniques, and analytic expressions. These results are compared with those from the hypernetted chain equation and the Verlet-modified closure against Monte Carlo and equations-of-state data for reduced densities of ${\rho }^{*}<1$ and temperatures of $T^{*}=1.5$, 2.74, and 5. Our closure shows consistency among all thermodynamic paths, except for one expression of the Gibbs-Duhem relation, whereas the hypernetted chain equation and the Verlet-modified closure exhibit consistency between only a few relations. Accuracy of the closure is comparable to the Verlet-modified closure and a significant improvement to results obtained from the hypernetted chain equation.

We present a mechanism to generate a sequence of shearless curves or attractors to form a band of transport barriers. We consider the labyrinthic nontwist standard map to prepare a scenario with three shearless curves. Dissipation is introduced and three shearless attractors coexist, very close to each other. In both cases a collective transport barrier is formed.

Quantum chaotic kicked top model is implemented experimentally in a two-qubit system comprising of a pair of spin-1/2 nuclei using nuclear magnetic resonance techniques. The essential nonlinear interaction was realized using indirect spin-spin coupling, while the linear kicks were realized using radio-frequency pulses. After a variable number of kicks, quantum state tomography was employed to reconstruct the single-qubit reduced density matrices, using which we could extract von Neumann entropies and Husimi distributions. These measures enabled the study of correspondence with classical phase space as well as probing distinct features of quantum chaos, such as symmetries and temporal periodicity in the two-qubit kicked top.

We extend the one-dimensional cochlear model to include the effect of viscosity and compressibility of the cochlear fluid. The resulting boundary-value problem is solved exactly using numerical techniques and semianalytically using the WKB approximation. Our results indicate the general trend of basilar membrane response increasing with an increase in compressibility or a decrease in viscosity. However, in the physiologically relevant range of these parameters, the change in the response is insignificant, justifying the assumption made in one-dimensional cochlear models that these effects are unimportant. Using the semianalytical WKB algorithm, we also demonstrate the simultaneous existence of forward- and backward-traveling waves on the basilar membrane.

In this work we study the dynamical buckling process of a thin filament immersed in a highly viscous medium. We perform an experimental study to track the shape evolution of the filament during a constant velocity compression. Numerical simulations reproduce the dynamical features observed from the experimental data and allow quantifying the filament's load. We observe that both the filament's load and the wave number evolve in a stepwise manner. In order to achieve a physical insight of the process, we apply a theoretical model to describe the buckling of a filament in a viscous medium. We solve a hydrodynamic equation in terms of normal modes for clamped-clamped boundary conditions and constant applied load. We find a good agreement between experimental data and simulations, suggesting that the proposed mechanistic model captures the essential features underlying the dynamical buckling process.

In this paper, we consider the Boltzmann equation with respect to orthonormal vielbein fields in conservative form. This formalism allows the use of arbitrary coordinate systems to describe the space geometry, as well as of an adapted coordinate system in the momentum space, which is linked to the physical space through the use of vielbeins. Taking advantage of the conservative form, we derive the macroscopic equations in a covariant tensor notation, and show that the hydrodynamic limit can be obtained via the Chapman-Enskog expansion in the Bhatnaghar-Gross-Krook approximation for the collision term. We highlight that in this formalism, the component of the momentum which is perpendicular to some curved boundary can be isolated as a separate momentum coordinate, for which the half-range Gauss-Hermite quadrature can be applied. We illustrate the capabilities of this formalism by considering two applications. The first one is the circular Couette flow between rotating coaxial cylinders, for which benchmarking data are available for all degrees of rarefaction, from the hydrodynamic to the ballistic regime. The second application concerns the flow in a gradually expanding channel. We employ finite-difference lattice Boltzmann models based on half-range Gauss-Hermite quadratures for the implementation of diffuse reflection, together with the fifth-order weighted essentially nonoscillatory and third-order total variation diminishing Runge-Kutta numerical methods for the advection and time stepping, respectively.

This Rapid Communication reports on an experimental study of granular craters formed by a mechanism, namely, optical energy, via a pulsed laser focused onto the surface of a granular bed. This represents an insight into granular cratering for two reasons; first, there is no physical contact between the initiation mechanism and the granular media (as typical for impact or explosion craters). Second, the resulting craters are millimetric in scale, which facilitates a test of energy scalings down to a previously unobserved lengthscale. Indeed, we observe a range of energy scalings conforming to $D_c\sim E^{\beta }$ with $\beta \approx 0.31\text{–}0.43$ depending on the characteristics of the granular media.

We address the calculation of transition probabilities in multiplicative noise stochastic differential equations using a path integral approach. We show the equivalence between the conditional probability and the propagator of a quantum particle with variable mass. Introducing a time reparametrization, we are able to transform the problem of multiplicative noise fluctuations into an equivalent additive one. We illustrate the method by showing the explicit analytic computation of the conditional probability of a harmonic oscillator in a nonlinear multiplicative environment.

In this paper, after a brief recall of the derivation of the unidirectional pulse propagation equation generalized to structured media, a fast modal transform linking the spatiotemporal representation of the field and its modal distribution is presented. This transform is used for solving the propagation equation by using a split-step algorithm in an efficient way. As an example, we present, to the best of our knowledge, the first numerical evidence of the generation of conical waves in highly multimodes waveguides.

Environments with immobile obstacles or void regions that inhibit and alter the motion of individuals within that environment are ubiquitous. Correlation in the location of individuals within such environments arises as a combination of the mechanisms governing individual behavior and the heterogeneous structure of the environment. Measures of spatial structure and correlation have been successfully implemented to elucidate the roles of the mechanisms underpinning the behavior of individuals. In particular, the pair correlation function has been used across biology, ecology, and physics to obtain quantitative insight into a variety of processes. However, naively applying standard pair correlation functions in the presence of obstacles may fail to detect correlation, or suggest false correlations, due to a reliance on a distance metric that does not account for obstacles. To overcome this problem, here we present an analytic expression for calculating a corrected pair correlation function for lattice-based domains containing obstacles. We demonstrate that this obstacle pair correlation function is necessary for isolating the correlation associated with the behavior of individuals, rather than the structure of the environment. Using simulations that mimic cell migration and proliferation we demonstrate that the obstacle pair correlation function recovers the short-range correlation known to be present in this process, independent of the heterogeneous structure of the environment. Further, we show that the analytic calculation of the obstacle pair correlation function derived here is significantly faster to implement than the corresponding numerical approach.

We investigate the dynamics of an ensemble of smooth inelastic hard spheres confined between two horizontal plates separated by a distance smaller than twice the diameter of the particles, in such a way that the system is quasi-two-dimensional. The bottom wall is vibrating and, therefore, it injects energy into the system in the vertical direction and a stationary state is reached. It is found that if the size of the plates is small enough, the stationary state is homogeneous. Otherwise, a cluster of particles is developed. The instability is understood by using some effective hydrodynamic equations in the horizontal plane. Moreover, the theoretical prediction for the size of the system above which it becomes unstable agrees very well with molecular dynamics simulation results without any fitting parameter.

Available macroscopic theories—such as the Johnson-Kendall-Roberts (JKR) model—predict spherical particles to stick to each other at small collision velocities $v$; above the bouncing velocity, $v_b$, they bounce. We study the details of the bouncing threshold using molecular dynamics simulation for crystalline nanoparticles where atoms interact via the Lennard-Jones potential. We show that the bouncing velocity strongly depends on the nanoparticle orientation during collision; for some orientations, nanoparticles stick at all velocities. The dependence of bouncing on orientation is caused by energy dissipation during dislocation activity. The bouncing velocity decreases with increasing nanoparticle radius in reasonable agreement with JKR theory. For orientations for which bouncing exists, nanoparticles stick again at a higher velocity, the fusion velocity, $v_f$, such that bouncing only occurs in a finite range of velocities—the bouncing window. The fusion velocity is rather independent of the nanoparticle radius.

We develop a numerical model for a two-phase flow of a pair of immiscible liquids within a capillary tube. We assume that a capillary is initially saturated with one liquid and the other liquid is injected via one of the capillary's ends. The governing equations are solved in the velocity-pressure formulation, so the pressure levels are imposed at the capillary inlet and outlet ends. We model the structure of the flow and the shape of the interface. We are able to reproduce the flow for a wide range of capillary numbers, when the meniscus preserves its shape moving together with the flow, and when the meniscus constantly stretches resembling the transport of a passive impurity. We demonstrate that the phase-field approach is capable of reproducing all features of the liquid-liquid displacement, including the motion of a contact line, the dynamic changes of the capillary pressure, and the dynamic changes of the apparent contact angle.

We consider a Brownian particle diffusing in a one-dimensional interval with absorbing end points. We study the ramifications when such motion is interrupted and restarted from the same initial configuration. We provide a comprehensive study of the first-passage properties of this trapping phenomena. We compute the mean first-passage time and derive the criterion on which restart always expedites the underlying completion. We show how this set-up is a manifestation of a success-failure problem. We obtain the success and failure rates and relate them with the splitting probabilities, namely the probability that the particle will eventually be trapped on either of the boundaries without hitting the other one. Numerical studies are presented to support our analytic results.

We study a metacommunity model of consumer-resource populations coupled via dispersal under an environment-dependent framework, and we explore the occurrence of multistability and critical transition. By emphasizing two magnitudes acting on a dynamic environment at temporal and spatial scales, the coupled system with simple diffusive coupling and the nonlinear environmental coupling enables various interesting complex dynamics such as bistability, multistability, and critical transitions. Using the basin stability measure, we find the probability of attaining each alternative state in a multistable region. In addition, critical transitions (one from a high to a low species density and the other from a low to a high species density) are identified at different magnitudes in the presence of stochastic fluctuations. We also explore the robustness of critical slowing-down indicators, e.g., lag-1 autocorrelation and variance, to forewarn the critical transition in the metacommunity model. Further, a network structure also identifies synchronization and multiclustering for a different choice of initial conditions. In contrast with the earlier studies on dynamic environmental coupling, our results based on the defined magnitudes provide important insights into environmental heterogeneity, which determines the set of environmental conditions to predict metacommunity stability and persistence.

In this paper, we present experimental observations and phenomenological modeling of the intermittent dynamics that emerge while mitigating thermoacoustic instability by rotating the otherwise static swirler in a lean premixed, laboratory-scale combustor. Starting with a self-excited thermoacoustically unstable combustor, here we find that a progressive increase in swirler rotation rate does not uniformly decrease amplitudes of coherent, sinusoidal pressure or heat-release-rate oscillations. Instead, these oscillations emerge as high-amplitude bursts separated by low-amplitude noise in the signal. At increased rotational speeds, the high-amplitude coherent oscillations become scarce and their duration in the signal reduces. The velocity field from high-speed particle image velocimetry and simultaneous pressure and chemiluminescence data support these observations. Such an intermittent route to instability mitigation is reminiscent of the opposite transition implemented by changing the Reynolds number from a fully chaotic state to a fully unstable state. To model such dynamics phenomenologically, we discretize the swirling turbulent premixed flame into an ensemble of flamelet oscillators arranged circumferentially around the center body of the swirler. The Kuramoto model is proposed for these flamelet oscillators which is subsequently used to analyze their synchronization dynamics. The order parameter $r$, which is a measure of the synchronization between the oscillator phases, provides critical insights on the transition from the thermoacoustically unstable to stable states via intermittency. Finally, it is shown that the Kuramoto model for flamelet oscillator can qualitatively reproduce the time-averaged and intermittent dynamics while transitioning from the state of thermoacoustic instability to a state of incoherent noisy oscillations.

An adaptive coupling based on a low-pass filter (LPF) is proposed to manipulate dynamic activity of diffusively coupled dynamical systems. A theoretical analysis shows that tracking either the external or internal signal in the coupling via a LPF gives rise to distinctly different ways of regulating the rhythmicity of the coupled systems. When the external signals of the coupling are attenuated by a LPF, the macroscopic oscillations of the coupled system are quenched due to the emergence of amplitude or oscillation death. If the internal signals of the coupling are further filtered by a LPF, amplitude and oscillation deaths are effectively revoked to restore dynamic behaviors. The applicability of this approach is demonstrated in laboratory experiments of coupled oscillatory electrochemical reactions by inducing coupling through LPFs. Our study provides additional insight into (ar)rhythmogenesis in diffusively coupled systems.

The simplest model of DNA mechanics describes the double helix as a continuous rod with twist and bend elasticity. Recent work has discussed the relevance of a little-studied coupling $G$ between twisting and bending, known to arise from the groove asymmetry of the DNA double helix. Here the effect of $G$ on the statistical mechanics of long DNA molecules subject to applied forces and torques is investigated. We present a perturbative calculation of the effective torsional stiffness $C_{\text{eff}}$ for small twist-bend coupling. We find that the “bare” $G$ is “screened” by thermal fluctuations, in the sense that the low-force, long-molecule effective free energy is that of a model with $G=0$ but with long-wavelength bending and twisting rigidities that are shifted by $G$-dependent amounts. Using results for torsional and bending rigidities for freely fluctuating DNA, we show how our perturbative results can be extended to a nonperturbative regime. These results are in excellent agreement with numerical calculations for Monte Carlo “triad” and molecular dynamics “oxDNA” models, characterized by different degrees of coarse graining, validating the perturbative and nonperturbative analyses. While our theory is in generally good quantitative agreement with experiment, the predicted torsional stiffness does systematically deviate from experimental data, suggesting that there are as-yet-uncharacterized aspects of DNA twisting-stretching mechanics relevant to low-force, long-molecule mechanical response, which are not captured by widely used coarse-grained models.