We have studied the electrical conductance of two-dimensional (2D) random percolating networks of zero-width metallic nanowires (a mixture of rings and sticks). We took into account the nanowire resistance per unit length and the junction (nanowire-nanowire contact) resistance. Using a mean-field approximation (MFA) approach, we derived the total electrical conductance of these nanowire-based networks as a function of their geometrical and physical parameters. The MFA predictions have been confirmed by our Monte Carlo (MC) numerical simulations. The MC simulations were focused on the case when the circumferences of the rings and the lengths of the wires were equal. In this case, the electrical conductance of the network was found to be almost insensitive to the relative proportions of the rings and sticks, provided that the wire resistance and the junction resistance were equal. When the junction resistance dominated over the wire resistance, a linear dependency of the electrical conductance of the network on the proportions of the rings and sticks was observed.

We observe multiple steplike jumps in a Dy-Fe-Ga-based ferrimagnetic alloy in its magnetic hysteresis curve at 2 K. The observed jumps are found to have a stochastic character with respect to their magnitude and the field position, and the jumps do not correlate with the duration of the field. The distribution of jump size follows a power law variation indicating the scale invariance nature of the jumps. We have invoked a simple two-dimensional random bond Ising-type spin system to model the dynamics. Our computational model can qualitatively reproduce the jumps and their scale-invariant character. It also elucidates that the flipping of antiferromagnetically coupled Dy and Fe clusters is responsible for the observed jumps in the hysteresis loop. These features are described in terms of the self-organized criticality.

This work identifies a solvable (in the sense that spectral correlation functions can be expressed in terms of orthogonal polynomials), rotationally invariant random matrix ensemble with a logarithmic weakly confining potential. The ensemble, which can be interpreted as a transformed Jacobi ensemble, is in the thermodynamic limit characterized by a Lorentzian eigenvalue density. It is shown that spectral correlation functions can be expressed in terms of the nonclassical Gegenbauer polynomials $C_n^{(-1/2)}\left(x\right)$ with $n\ge 2$, which have been proven to form a complete orthogonal set with respect to the proper weight function. A procedure to sample matrices from the ensemble is outlined and used to provide a numerical verification for some of the analytical results. This ensemble is pointed out to potentially have applications in quantum many-body physics.

We study the localization properties of normal modes in harmonic chains with mass and spring weak disorder. Using a perturbative approach, an expression for the localization length $L_{\text{loc}}$ is obtained, which is valid for arbitrary correlations of the disorder (mass disorder correlations, spring disorder correlations, and mass-spring disorder correlations are allowed), and for practically the whole frequency band. In addition, we show how to generate effective mobility edges by the use of disorder with long range self-correlations and cross-correlations. The transport of phonons is also analyzed, showing effective transparent windows that can be manipulated through the disorder correlations even for relative short chain sizes. These results are connected to the problem of heat conduction in the harmonic chain; indeed, we discuss the size scaling of the thermal conductivity from the perturbative expression of $L_{\text{loc}}$. Our results may have applications in modulating thermal transport, particularly in the design of thermal filters or in manufacturing high-thermal-conductivity materials.

We reexamine results obtained with the recently proposed density functional theory framework based on forces (force-DFT) [S. M. Tschopp et al., Phys. Rev. E 106, 014115 (2022)]. We compare inhomogeneous density profiles for hard sphere fluids to results from both standard density functional theory and from computer simulations. Test situations include the equilibrium hard sphere fluid adsorbed against a planar hard wall and the dynamical relaxation of hard spheres in a switched harmonic potential. The comparison to grand canonical Monte Carlo simulation profiles shows that equilibrium force-DFT alone does not improve upon results obtained with the standard Rosenfeld functional. Similar behavior holds for the relaxation dynamics, where we use our event-driven Brownian dynamics data as benchmark. Based on an appropriate linear combination of standard and force-DFT results, we investigate a simple hybrid scheme which rectifies these deficiencies in both the equilibrium and the dynamical case. We explicitly demonstrate that although the hybrid method is based on the original Rosenfeld fundamental measure functional, its performance is comparable to that of the more advanced White Bear theory.

Using the path integral representation of the nonequilibrium dynamics, we compute the most probable path between arbitrary starting and final points that is followed by an active particle driven by persistent noise. We focus our attention on the case of active particles immersed in harmonic potentials, where the trajectory can be computed analytically. Once we consider the extended Markovian dynamics where the self-propulsive drive evolves according to an Ornstein-Uhlenbeck process, we can compute the trajectory analytically with arbitrary conditions on position and self-propulsion velocity. We test the analytical predictions against numerical simulations and we compare the analytical results with those obtained within approximated equilibriumlike dynamics.

We modeled the change in the temperature of maximum density (TMD) of a waterlike solvent when small amounts of solute are added to the mixture. The solvent is modeled as a two length scales potential, which is known to exhibit waterlike characteristic anomalies, while the solute is chosen to have an attractive interaction with the solvent which is tuned from small to large attractive potential. We show that if the solute exhibits high attraction with the solvent it behaves as a structure maker and the TMD increases with the addition of solute, while if the solute shows a low attraction with the solvent the TMD decreases, with the solute behaving as a structure breaker.

The fracture behavior of brittle and ductile materials can be strongly influenced by thermal fluctuations, especially in micro- and nanodevices as well as in rubberlike and biological materials. However, temperature effects, in particular on the brittle-to-ductile transition, still require a deeper theoretical investigation. As a step in this direction we propose a theory, based on equilibrium statistical mechanics, able to describe the temperature-dependent brittle fracture and brittle-to-ductile transition in prototypical discrete systems consisting in a lattice with breakable elements. Concerning the brittle behavior, we obtain closed form expressions for the temperature-dependent fracture stress and strain, representing a generalized Griffith criterion, ultimately describing the fracture as a genuine phase transition. With regard to the brittle-to-ductile transition, we obtain a complex critical scenario characterized by a threshold temperature between the two fracture regimes (brittle and ductile), an upper and a lower yield strength, and a critical temperature corresponding to the complete breakdown. To show the effectiveness of the proposed models in describing thermal fracture behaviors at small scales, we successfully compare our theoretical results with molecular dynamics simulations of Si and GaN nanowires.

We study ultrasonic resonance in a coflow system comprising a pair of immiscible liquids in a microchannel exposed to bulk acoustic waves. We show using an analytical model that there are two resonating frequencies corresponding to each of the coflowing liquids, which depend on the speed of sound and stream width of the liquid. We perform a frequency domain analysis using numerical simulations to reveal that resonance can be achieved by actuating both liquids at a single resonating frequency that depends on the speeds of sound, densities, and widths of the liquids. In a coflow system with equal speeds of sound and densities of the pair of fluids, the resonating frequency is found to be independent of the relative width of the two streams. In coflow systems with unequal speeds of sound or densities, even with matching characteristic acoustic impedances, the resonating frequency depends on the stream width ratio, and the value increases with an increase in the stream width of the liquid with a higher speed of sound. We show that a pressure nodal plane can be realized at the channel center by operating at a half-wave resonating frequency when the speeds of sound and densities are equal. However, the pressure nodal plane is found to shift away from the center of the microchannel when the speeds of sound and densities of the two liquids are unequal. The results of the model and simulations are verified experimentally via acoustic focusing of microparticles suggesting the formation of a pressure nodal plane and hence a resonance condition. Our study will find relevance in acoustomicrofluidics involving immiscible coflow systems.

Giant half-cycle attosecond pulse generation in the coherent bremsstrahlung emission regime is proposed for laser pulses with normal incidence on a double-foil target, where the first foil is transparent and the second foil is opaque. The presence of the second opaque target contributes to the formation of a relativistic flying electron sheet (RFES) from the first foil target. After the RFES has passed through the second opaque target, it is decelerated sharply, and bremsstrahlung emission occurs, which results in the generation of an isolated half-cycle attosecond pulse having an intensity of $\sim 1.4\times {10}^{22}\mathrm{W}/{\mathrm{cm}}^2$ and a duration of 3.6 $\mathrm{as}$. The generation mechanism does not require extra filters and may open a regime of nonlinear attosecond science.

Motivated by the interplay of multiple species in several real world transport processes, we propose a bidirectional totally asymmetric simple exclusion process with two finite particle reservoirs regulating the inflow of oppositely directed particles corresponding to two different species. The system's stationary characteristics, such as densities, currents, etc., are investigated using a theoretical framework based on mean-field approximation and are supported by extensive Monte Carlo simulations. The impact of individual species populations, quantified by filling factor, has been comprehensively analyzed considering both equal and unequal conditions. For the equal case, the system exhibits the spontaneous symmetry-breaking phenomena and admits both symmetric as well as asymmetric phases. Moreover, the phase diagram exhibits a different asymmetric phase and displays a nonmonotonic variation in the number of phases with respect to the filling factor. For unequal filling factors, the phase schema can display at most five phases including a phase that shows maximal current for one of the species.

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate $p$, and their autocatalytic creation at unoccupied sites with $n\ge 2$ occupied neighbors at rate $k_n$. Kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the $p$ value for equistability of coexisting populated and vacuum states, $p_{\mathrm{eq}}\left(S\right)$, depends on the orientation or slope, $S$, of a planar interface separating those phases. The vacuum state displaces the populated state for $p>p_{\mathrm{eq}}\left(S\right)$, and the opposite applies for $p<p_{\mathrm{eq}}\left(S\right)$ for $0<S<\infty$. The special “combinatorial” rate choice $k_n=n\left(n-1\right)/12$ facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that $p_{\mathrm{eq}}\left(\max \right)=p_{\mathrm{eq}}(S=1)=0.09645$ and $p_{\mathrm{eq}}\left(\min \right)=p_{\mathrm{eq}}(S\to \infty )=0.08827$, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all $p<p_{\mathrm{eq}}(S=\infty )=0.08907$, a value exceeding $p_{\mathrm{eq}}(S\to \infty )$. One can regard an interface for large $S\to \infty$ as a vertical interface decorated with isolated kinks. For $p<p_{\mathrm{eq}}\left(S=\infty \right)$, the kink can move in either direction along this otherwise stationary interface depending upon $p$, but for $p=p_{\mathrm{eq}}\left(\min \right)$ the kink is also stationary.

A density oscillator is a fluid system in which oscillatory flow occurs between different density fluids through the pore connecting them. We investigate the synchronization in coupled density oscillators using two-dimensional hydrodynamic simulation and analyze the stability of the synchronous state based on the phase reduction theory. Our results show that the antiphase, three-phase, and 2-2 partial-in-phase synchronization modes spontaneously appear as stable states in two, three, and four coupled oscillators, respectively. The phase dynamics of coupled density oscillators is interpreted with their sufficiently large first Fourier components of the phase coupling function.

We study the wrinkling dynamics of three-dimensional vesicles in a time-dependent elongation flow by utilizing an immersed boundary method. For a quasispherical vesicle, our numerical results well match the predictions of perturbation analysis, where similar exponential relationships between wrinkles' characteristic wavelength and the flow strength are observed. Using the same parameters as in the experiments by Kantsler et al. [V. Kantsler et al., Phys. Rev. Lett. 99, 178102 (2007)], our simulations of an elongated vesicle are in good agreement with their results. In addition, we get rich three-dimensional morphological details, which are favorable to comprehend the two-dimensional snapshots. This morphological information helps identify wrinkle patterns. We analyze the morphological evolution of wrinkles using spherical harmonics. We find discrepancies in elongated vesicle dynamics between simulations and perturbation analysis, highlighting the importance of the nonlinear effects. Finally, we investigate the unevenly distributed local surface tension, which largely determines the position of wrinkles excited on the vesicle membrane.

This paper extends the partially saturated method (PSM), used for curved or complex walls, to the lattice Boltzmann (LB) pseudopotential multicomponent model and adapts the wetting boundary condition to model the contact angle. The pseudopotential model is widely used for various complex flow simulations due to its simplicity. To simulate the wetting phenomenon within this model, the mesoscopic interaction force between the boundary fluid and solid nodes is used to mimic the microscopic adhesive force between the fluid and the solid wall, and the bounce-back (BB) method is normally adopted to achieve the no-slip boundary condition. In this paper, the pseudopotential interaction forces are computed with eighth-order isotropy since fourth-order isotropy leads to the condensation of the dissolved component on curved walls. Due to the staircase approximation of curved walls in the BB method, the contact angle is sensitive to the shape of corners on curved walls. Furthermore, the staircase approximation makes the movement of the wetting droplet on curved walls not smooth. To solve this problem, the curved boundary method may be used, but due to the interpolation or extrapolation process, most curved boundary conditions suffer from massive mass leakage when applied to the LB pseudopotential model. Through three test cases, it is found that the improved PSM scheme is mass conservative, that nearly identical static contact angles are observed on flat and curved walls under the same wetting condition, and that the movement of a wetting droplet on curved and inclined walls is smoother compared to the usual BB method. The present method is expected to be a promising tool for modeling flows in porous media and in microfluidic channels.

We introduce a family of generalized continuous Maxwell demons (GCMDs) operating on idealized single-bit equilibrium devices that combine the single-measurement Szilard and the repeated measurements of the continuous Maxwell demon protocols. We derive the cycle distributions for extracted work, information content, and time and compute the power and information-to-work efficiency fluctuations for the different models. We show that the efficiency at maximum power is maximal for an opportunistic protocol of continuous type in the dynamical regime dominated by rare events. We also extend the analysis to finite-time work extracting protocols by mapping them to a three-state GCMD. We show that dynamical finite-time correlations in this model increase the information-to-work conversion efficiency, underlining the role of temporal correlations in optimizing information-to-energy conversion. The effect of finite-time work extraction and demon memory resetting is also analyzed. We conclude that GCMD models are thermodynamically more efficient than the single-measurement Szilard and preferred for describing biological processes in an information-redundant world.

An exact expression for the average velocity of cold atoms in a driven, dissipative optical lattice in terms of the amplitudes of atomic density waves is derived from semiclassical equations for the phase space densities of the Zeeman ground-state sublevels. The calculations are for a $J_g=1/2\to J_e=3/2$ transition, as it is customary in theoretical studies of Sisyphus cooling. While the driver, an additional beam of small amplitude, sets the atoms into directed motion, the new expression permits the quantification of the contribution to the atomic motion of a specific atomic wave, revealing unexpected counterpropagating contributions from many modes. Additionally, the method is shown to provide the generic threshold for the transition into the regime of infinite density, regardless of the details, or even the presence, of driving.

This work is dedicated to the topological analysis of complex transitional networks for dynamic state detection. Transitional networks are formed from time series data and they leverage graph theory tools to reveal information about the underlying dynamic system. However, traditional tools can fail to summarize the complex topology present in such graphs. In this work, we leverage persistent homology from topological data analysis to study the structure of these networks. We contrast dynamic state detection from time series using a coarse-grained state-space network (CGSSN) and topological data analysis (TDA) to two state of the art approaches: ordinal partition networks (OPNs) combined with TDA and the standard application of persistent homology to the time-delay embedding of the signal. We show that the CGSSN captures rich information about the dynamic state of the underlying dynamical system as evidenced by a significant improvement in dynamic state detection and noise robustness in comparison to OPNs. We also show that because the computational time of CGSSN is not linearly dependent on the signal's length, it is more computationally efficient than applying TDA to the time-delay embedding of the time series.

We study the transport properties of diffusing particles restricted to confined regions on curved surfaces. We relate particle mobility to the curvature of the surface where they diffuse and the constraint due to confinement. Applying the Fick-Jacobs procedure to diffusion in curved manifolds shows that the local diffusion coefficient is related to average geometric quantities such as constriction and tortuosity. Macroscopic experiments can record such quantities through an average surface diffusion coefficient. We test the accuracy of our theoretical predictions of the effective diffusion coefficient through finite-element numerical solutions of the Laplace-Beltrami diffusion equation. We discuss how this work contributes to understanding the link between particle trajectories and the mean-square displacement.