The resonance-absorption condition in the laser-nanoplasma interactions has been considered to follow the wavelength dependence of the critical plasma density. We experimentally demonstrate that this assumption fails in the middle-infrared spectral range, while it is valid for visible and near-infrared wavelengths. A thorough analysis supported by molecular dynamic (MD) simulations indicates that the observed transition in the resonance condition is caused by the reduction of the electron scattering rate and the associated increase of the cluster outer-ionization contribution. An expression for the nanoplasma resonance density is derived based on experimental results and MD simulations. The findings are important for a broad range of plasma experiments and applications, since the extension of the laser-plasma interaction studies to longer wavelengths has become increasingly topical.

In this paper, we first develop a fourth-order multiple-relaxation-time lattice Boltzmann (MRT-LB) model for the one-dimensional convection-diffusion equation (CDE) with the constant velocity and diffusion coefficient, where the D1Q3 (three discrete velocities in one-dimensional space) lattice structure is used. We also perform the Chapman-Enskog analysis to recover the CDE from the MRT-LB model. Then an explicit four-level finite-difference (FLFD) scheme is derived from the developed MRT-LB model for the CDE. Through the Taylor expansion, the truncation error of the FLFD scheme is obtained, and at the diffusive scaling, the FLFD scheme can achieve the fourth-order accuracy in space. After that, we present a stability analysis and derive the same stability condition for the MRT-LB model and FLFD scheme. Finally, we perform some numerical experiments to test the MRT-LB model and FLFD scheme, and the numerical results show that they have a fourth-order convergence rate in space, which is consistent with our theoretical analysis.

In a previous work [R. Goerlich et al., Phys. Rev. E 106, 054617 (2022)] the authors studied the transition from one nonequilibrium steady state (NESS) to another NESS by changing the correlated noise that is driving a Brownian particle held in an optical trap. They found that the amount of heat that is released during the transition is directly proportional to the difference of spectral entropy between the two colored noises, in a fashion that is reminiscent of Landauer's principle. In this Comment I argue that the relation found between the released heat and the spectral entropy does not hold in general and that one can provide examples of noises where it clearly fails. I also show that, even in the case considered by the authors, the relation cannot be rigorously true and is only approximately verified experimentally.

We stress that the limitations on one of the results of our paper [R. Goerlich et al., Phys. Rev. E 106, 054617 (2022)], which are mentioned in the preceding Comment [A. Bérut, preceding Comment, Phys. Rev. E 107, 056601 (2023)], were actually already acknowledged and discussed in the original publication. Although the observed relationship between the released heat and the spectral entropy of the correlated noise is not universal (but limited to one-parameter Lorentzian spectra), the existence of such a clear relationship is a solid experimental finding. It not only gives a convincing explanation for the surprising thermodynamics observed in the transitions between nonequilibrium steady states, but also provides new tools for the analysis of nontrivial baths. In addition, by using different measures of the correlated noise information content, it may be possible to generalize these results to non-Lorentzian spectra.

We investigate analytically and computationally the dynamics of two-dimensional needle crystal growth from the melt in a narrow channel. Our analytical theory predicts that, in the low supersaturation limit, the growth velocity $V$ decreases in time $t$ as a power law $V\sim t^{-2/3}$, which we validate by phase-field and dendritic-needle-network simulations. Simulations further reveal that, above a critical channel width $\mathrm{\Lambda }\approx 5l_D$, where $l_D$ is the diffusion length, needle crystals grow with a constant $V<V_s$, where $V_s$ is the free-growth needle crystal velocity, and approaches $V_s$ in the limit $\mathrm{\Lambda }\gg l_D$.

Granular materials are used in several fields and in a wide variety of processes. An important feature of these materials is the diversity of grain sizes, commonly referred to as polydispersity. When granular materials are sheared, they exhibit a predominant small elastic range. Then, the material yields, with or without a peak shear strength depending on the initial density. Finally, the material reaches a stationary state, in which it deforms at a constant shear stress, which can be linked to the residual friction angle ${\varphi }_{\mathrm{r}}$. However, the role of polydispersity on the shear strength of granular materials is still a matter of debate. In particular, a series of investigations have proved, using numerical simulations, that ${\varphi }_{\mathrm{r}}$ is independent of polydispersity. This counterintuitive observation remains elusive to experimentalists, and especially for some technical communities that use ${\varphi }_{\mathrm{r}}$ as a design parameter (e.g., the soil mechanics community). In this Letter, we studied experimentally the effects of polydispersity on ${\varphi }_{\mathrm{r}}$. In order to do so, we built samples of ceramic beads and then sheared these samples in a triaxial apparatus. We varied polydispersity, building monodisperse, bidisperse, and polydisperse granular samples; this allowed us to study the effects of grain size, size span, and grain size distribution on ${\varphi }_{\mathrm{r}}$. We find that ${\varphi }_{\mathrm{r}}$ is indeed independent of polydispersity, confirming the previous findings achieved through numerical simulations. Our work fairly closes the gap of knowledge between experiments and simulations.

By using trajectory averaging to analyze the statistics of energy dissipation in the nonequilibrium energy-state transitions of a driven two-state system, we show that the average energy dissipation induced by external driving is connected to its fluctuations about equilibrium through the simple relation $2k_{\mathrm{B}}T\langle Q\rangle =\langle \delta Q^2\rangle$, which is preserved by an adiabatic approximation scheme. We use this scheme to obtain the heat statistics of a single-electron box with a superconducting lead in the slow-driving regime, where the dissipated heat becomes normally distributed with a relatively high probability to be extracted from the environment rather than dissipated. We also discuss the validity of heat fluctuation relations beyond driven two-state transitions and the slow-driving regime.

This work focuses on the study of time-periodic solutions, including breathers, in a nonlinear lattice consisting of elements whose contacts alternate between strain hardening and strain softening. The existence, stability, and bifurcation structure of such solutions, as well as the system dynamics in the presence of damping and driving, are studied systematically. It is found that the linear resonant peaks in the system bend toward the frequency gap in the presence of nonlinearity. The time-periodic solutions that lie within the frequency gap compare well to Hamiltonian breathers if the damping and driving are small. In the Hamiltonian limit of the problem, we use a multiple scale analysis to derive a nonlinear Schrödinger equation to construct both acoustic and optical breathers. The latter compare very well with the numerically obtained breathers in the Hamiltonian limit.

The stiffness of biological membranes determines the work required by cellular machinery to form and dismantle vesicles and other lipidic shapes. Model membrane stiffness can be determined from the equilibrium distribution of giant unilamellar vesicle surface undulations observable by phase contrast microscopy. With two or more components, lateral fluctuations of composition will couple to surface undulations depending on the curvature sensitivity of the constituent lipids. The result is a broader distribution of undulations whose complete relaxation is partially determined by lipid diffusion. In this work, kinetic analysis of the undulations of giant unilamellar vesicles made of phosphatidylcholine-phosphatidylethanolamine mixtures validates the molecular mechanism by which the membrane is made 25% softer than a single-component one. The mechanism is relevant to biological membranes, which have diverse and curvature-sensitive lipids.

Self-replicability is a unique attribute observed in all living organisms, and the question of how the life was physically initiated could be equivalent to the question of how self-replicating informative polymers were formed in the abiotic material world. It has been suggested that the present DNA and proteins world was preceded by an RNA world in which genetic information of RNA molecules was replicated by the mutual catalytic function of RNA molecules. However, the important question of how the transition occurred from a material world to the very early pre-RNA world remains unsolved both experimentally and theoretically. We present an onset model of mutually catalytic self-replicative systems formed in an assembly of polynucleotides. A quantitative expression of the critical condition for the onset of growing fluctuation towards self-replication in this model is obtained by analytical and numerical calculations.

We present the inertial active dynamics of an Ornstein-Uhlenbeck particle in a piecewise sawtooth ratchet potential. Using the Langevin simulation and matrix continued fraction method (MCFM), the particle transport, steady-state diffusion, and coherence in transport are investigated in different parameter regimes of the model. Spatial asymmetry is found to be a key criterion for the possibility of directed transport in the ratchet. The MCFM results for net particle current of overdamped dynamics of the particle agree well with the simulation results. The simulated particle trajectories for the inertial dynamics and the corresponding position and velocity distribution functions reveal that the system passes through an activity-induced transition in the transport from the running phase to the locked phase of the dynamics. This is further corroborated by the mean square displacement (MSD) calculations, where the MSD gets suppressed with increase in the persistent duration of activity or self-propulsion in the medium and finally approaches zero for a very large value of self propulsion time. The nonmonotonic behavior of the particle current and Péclet number with self-propulsion time confirms that the particle transport and its coherence can be enhanced or reduced by fine tuning the persistent duration of activity. Moreover, for intermediate ranges of self-propulsion time as well as mass of the particle, even though the particle current shows a pronounced unusual maximum with mass, there is no enhancement in the Péclet number, instead the Péclet number decreases with mass, confirming the degradation of coherence in transport.

A ferrofluid droplet confined in a Hele-Shaw cell can be deformed into a stably spinning “gear,” using crossed magnetic fields. Previously, fully nonlinear simulation revealed that the spinning gear emerges as a stable traveling wave along the droplet's interface bifurcates from the trivial (equilibrium) shape. In this work, a center manifold reduction is applied to show the geometrical equivalence between a two-harmonic-mode coupled system of ordinary differential equations arising from a weakly nonlinear analysis of the interface shape and a Hopf bifurcation. The rotating complex amplitude of the fundamental mode saturates to a limit cycle as the periodic traveling wave solution is obtained. An amplitude equation is derived from a multiple-time-scale expansion as a reduced model of the dynamics. Then, inspired by the well-known delay behavior of time-dependent Hopf bifurcations, we design a slowly time-varying magnetic field such that the timing and emergence of the interfacial traveling wave can be controlled. The proposed theory allows us to determine the time-dependent saturated state resulting from the dynamic bifurcation and delayed onset of instability. The amplitude equation also reveals hysteresislike behavior upon time reversal of the magnetic field. The state obtained upon time reversal differs from the state obtained during the initial (forward-time) period, yet it can still be predicted by the proposed reduced-order theory.

The centrifugal force in flow through a curved channel initiates a hydrodynamic instability that results in the development of Dean vortices, a pair of counter-rotating roll cells across the channel that deflect the high velocity fluid in the center toward the outer (concave) wall. When this secondary flow toward the concave (outer) wall is too strong to be dissipated by viscous effects, an additional pair of vortices emerges near the outer wall. Combining numerical simulation and dimensional analysis, we find that the critical condition for the onset of the second vortex pair depends on ${\gamma }^{1/2}\mathrm{Dn}$ ($\gamma$: channel aspect ratio; $\mathrm{Dn}$: Dean number). We also investigate the development length for the additional vortex pair in channels with different aspect ratios and curvatures. The larger centrifugal force at higher Dean numbers produces the additional vortices further upstream, with the required development length being inversely proportional to the Reynolds number and increasing linearly with the radius of curvature of the channel.

The effect of helicity in magnetohydrodynamic turbulence on the effective turbulent magnetic diffusion is considered here. The helical correction to turbulent diffusivity is analytically calculated with the use of the renormalization group approach. In agreement with previous numerical findings, this correction is shown to be negative and proportional to the second power of the magnetic Reynolds number, when the latter is small. In addition, the helical correction to turbulent diffusivity is found to obey a power-law-type dependence on the wave number of the most energetic turbulent eddies, $k_{\ell }$, of the form $k_{\ell }^{-10/3}$.

We present a general and numerically efficient method for calculation of collision integrals for interacting quantum gases on a discrete momentum lattice. Here we employ the original analytical approach based on Fourier transform covering a wide spectrum of solid-state problems with various particle statistics and arbitrary interaction models, including the case of momentum-dependent interaction. The comprehensive set of the transformation principles is given in detail and realized as a computer Fortran 90 library FLBE (Fast Library for Boltzmann Equation).

We consider a sheared granular system experiencing intermittent dynamics of stick-slip type via discrete element simulations. The considered setup consists of a two-dimensional system of soft frictional particles sandwiched between solid walls, one of which is exposed to a shearing force. The slip events are detected using stochastic state space models applied to various measures describing the system. The amplitudes of the events spread over more than four decades and present two distinctive peaks, one for the microslips and the other for the slips. We show that the measures describing the forces between the particles provide earlier detection of an upcoming slip event than the measures based solely on the wall movement. By comparing the detection times obtained from the considered measures, we observe that a typical slip event starts with a local change in the force network. However, some local changes do not spread globally over the force network. For the changes that become global, we find that their size strongly influences the further behavior of the system. If the size of a global change is large enough, then it triggers a slip event; if it is not, then a much weaker microslip follows. Quantification of the changes in the force network is made possible by formulating clear and precise measures describing their static and dynamic properties.

Using the Jacobian matrix, we obtain a theoretical expression of rigidity and the density of states of two-dimensional amorphous solids consisting of frictional grains in the linear response to an infinitesimal strain, in which we ignore the dynamical friction caused by the slip processes of contact points. The theoretical rigidity agrees with that obtained by molecular dynamics simulations. We confirm that the rigidity is smoothly connected to the value in the frictionless limit. We find that there are two modes in the density of states for sufficiently small $k_T/k_N$, which is the ratio of the tangential to normal stiffness. Rotational modes exist at low frequencies or small eigenvalues, whereas translational modes exist at high frequencies or large eigenvalues. The location of the rotational band shifts to the high-frequency region with an increase in $k_T/k_N$ and becomes indistinguishable from the translational band for large $k_T/k_N$.

A mesoscopic simulation model to study the phase separation in a binary fluid mixture in three dimensions (3D) is presented here by augmenting the existing particle-based multiparticle collision dynamics (MPCD) algorithm. The approach describes the nonideal equation of the fluid state by incorporating the excluded-volume interaction between the two components within the framework of stochastic collision, which depends on the local fluid composition and velocity. Calculating the nonideal contribution to the pressure both from simulation and analytics shows the model to be thermodynamically consistent. A phase diagram to explore the range of parameters that give rise to phase separation in the model is investigated. The interfacial width and phase growth obtained from the model agree with the literature for a wide range of temperatures and parameters.

Linear diffusions are used to model a large number of stochastic processes in physics, including small mechanical and electrical systems perturbed by thermal noise, as well as Brownian particles controlled by electrical and optical forces. Here we use techniques from large deviation theory to study the statistics of time-integrated functionals of linear diffusions, considering three classes of functionals or observables relevant for nonequilibrium systems which involve linear or quadratic integrals of the state in time. For these, we derive exact results for the scaled cumulant generating function and the rate function, characterizing the fluctuations of observables in the long-time limit, and study in an exact way the set of paths or effective process that underlies these fluctuations. The results give a complete description of how fluctuations arise in linear diffusions in terms of effective forces that remain linear in the state or, alternatively, in terms of fluctuating densities and currents that solve Riccati-type equations. We illustrate these results using two common nonequilibrium models, namely, transverse diffusions in two dimensions involving a nonconservative rotating force, and two interacting particles in contact with heat baths at different temperatures.

The zero-temperature Ising model is known to reach a fully ordered ground state in sufficiently dense random graphs. In sparse random graphs, the dynamics gets absorbed in disordered local minima at magnetization close to zero. Here, we find that the nonequilibrium transition between the ordered and the disordered regime occurs at an average degree that slowly grows with the graph size. The system shows bistability: The distribution of the absolute magnetization in the reached absorbing state is bimodal, with peaks only at zero and unity. For a fixed system size, the average time to absorption behaves nonmonotonically as a function of average degree. The peak value of the average absorption time grows as a power law of the system size. These findings have relevance for community detection, opinion dynamics, and games on networks.