Nonlinear Fokker-Planck equations endowed with power-law diffusion terms have proven to be valuable tools for the study of diverse complex systems in physics, biology, and other fields. The nonlinearity appearing in these evolution equations can be interpreted as providing an effective description of a system of particles interacting via short-range forces while performing overdamped motion under the effect of an external confining potential. This point of view has been recently applied to the study of thermodynamical features of interacting vortices in type II superconductors. In the present work we explore an embedding of the nonlinear Fokker-Planck equation within a Vlasov equation, thus incorporating inertial effects to the concomitant particle dynamics. Exact time-dependent solutions of the $q$-Gaussian form (with compact support) are obtained for the Vlasov equation in the case of quadratic confining potentials.

In this paper we propose a protocol to suppress double-layer forces between two microspheres immersed in a dielectric medium, being one microsphere metallic at a controlled potential ${\psi }_{\mathrm{M}}$ and the other a charged one either metallic or dielectric. The approach is valid for a wide range of distances between them. We show that, for a given distance between the two microspheres, the double-layer force can be totally suppressed by simply tuning ${\psi }_{\mathrm{M}}$ up to values dictated by the linearized Poisson-Boltzmann equation. Our key finding is that such values can be substantially different from the ones predicted by the commonly used proximity force approximation, also known as the Derjaguin approximation, even in situations where the latter is expected to be accurate. The proposed procedure can be used to suppress the double-layer interaction in force spectroscopy experiments, thus paving the way for measurements of other surface interactions, such as Casimir dispersion forces.

Neutral grains made of the same dielectric material can attain considerable charges due to collisions and generate long-range interactions. We perform molecular dynamic simulations in three dimensions for a dilute, freely cooling granular gas of viscoelastic particles that exchange charges during collisions. As compared to the case of clustering of viscoelastic particles solely due to dissipation, we find that the electrostatic interactions due to collisional charging alter the characteristic size, morphology, and growth rate of the clusters. The average cluster size grows with time as a power law, whose exponent is relatively larger in the charged gas than the neutral case. The growth of the average cluster size is found to be independent of the ratio of characteristic Coulomb to kinetic energy, or equivalently, of the typical Bjerrum length. However, this ratio alters the crossover time of the growth. Both simulations and mean-field calculations based on Smoluchowski's equation suggest that a suppression of particle diffusion due to the electrostatic interactions helps in the aggregation process.

We introduce a powerful and easily implemented method to calculate the permeability of porous media at the pore scale using an approximation based on the Poiseulle equation to calculate permeability to fluid flow with a Laplace solver. The method consists of calculating the Euclidean distance map of the fluid phase to assign local conductivities and lends itself naturally to the treatment of multiscale problems. We compare with analytical solutions as well as experimental measurements and lattice Boltzmann calculations of permeability for Fontainebleau sandstone. The solver is significantly more stable than the lattice Boltzmann approach, uses less memory, and is significantly faster. Permeabilities are in excellent agreement over a wide range of porosities.

In the modern theory of critical phenomena, the liquid-vapor density diameter in simple fluids is generally expected to deviate from a rectilinear law approaching the critical point. However, by performing precise scannerlike optical measurements of the position of the ${\mathrm{SF}}_6$ liquid-vapor meniscus, in an approach much closer to criticality in temperature and density than earlier measurements, no deviation from a rectilinear diameter can be detected. The observed meniscus position from far ($10\mathrm{K}$) to extremely close ($1\mathrm{mK}$) to the critical temperature is analyzed using recent theoretical models to predict the complete scaling consequences of a fluid asymmetry. The temperature dependence of the meniscus position appears consistent with the law of rectilinear diameter. The apparent absence of the critical hook in ${\mathrm{SF}}_6$ therefore seemingly rules out the need for the pressure scaling field contribution in the complete scaling theoretical framework in this ${\mathrm{SF}}_6$ analysis. More generally, this work suggests a way to clarify the experimental ambiguities in the simple fluids for the near-critical singularities in the density diameter.

Fractional Brownian motion, a stochastic process with long-time correlations between its increments, is a prototypical model for anomalous diffusion. We analyze fractional Brownian motion in the presence of a reflecting wall by means of Monte Carlo simulations. Whereas the mean-square displacement of the particle shows the expected anomalous diffusion behavior $\langle x^2\rangle \sim t^{\alpha }$, the interplay between the geometric confinement and the long-time memory leads to a highly non-Gaussian probability density function with a power-law singularity at the barrier. In the superdiffusive case $\alpha >1$, the particles accumulate at the barrier leading to a divergence of the probability density. For subdiffusion $\alpha <1$, in contrast, the probability density is depleted close to the barrier. We discuss implications of these findings, in particular, for applications that are dominated by rare events.

We address the issue of whether the golden-mean $[\psi =(\sqrt{5}+1)/2\simeq 1.618]$ universality class, as predicted by several theoretical models, can be observed in the dynamical scaling of thermal transport. Remarkably, we show strong evidence that $\psi$ appears to be the scaling exponent of heat mode correlation in a purely quartic anharmonic chain. This observation seems to somewhat deviate from the previous expectation and we explain it by the unusual slow decay of the cross correlation between heat and sound modes. Whenever the cubic anharmonicity is included, this cross correlation gradually dies out and another universality class with scaling exponent $\gamma =5/3$, as commonly predicted by theories, seems recovered. However, this recovery is accompanied by two interesting phase transition processes characterized by a change of symmetry of the potential and a clear variation of the dynamic structure factor, respectively. Due to these transitions, an additional exponent close to $\gamma \simeq 1.580$ emerges. All this evidence suggests that, to gain a full prediction of the scaling of thermal transport, more ingredients should be taken into account.

We address systematically an apparent nonphysical behavior of the free-energy moment generating function for several instances of the logarithmically correlated models: the fractional Brownian motion with Hurst index $H=0$ (fBm0) (and its bridge version), a one-dimensional model appearing in decaying Burgers turbulence with log-correlated initial conditions and, finally, the two-dimensional log-correlated random-energy model (logREM) introduced in Cao et al. [Phys. Rev. Lett. 118, 090601 (2017)] based on the two-dimensional Gaussian free field with background charges and directly related to the Liouville field theory. All these models share anomalously large fluctuations of the associated free energy, with a variance proportional to the log of the system size. We argue that a seemingly nonphysical vanishing of the moment generating function for some values of parameters is related to the termination point transition (i.e., prefreezing). We study the associated universal log corrections in the frozen phase, both for logREMs and for the standard REM, filling a gap in the literature. For the above mentioned integrable instances of logREMs, we predict the nontrivial free-energy cumulants describing non-Gaussian fluctuations on the top of the Gaussian with extensive variance. Some of the predictions are tested numerically.

We characterize the soft modes of the dynamical matrix at the depinning transition, and compare the matrix with the properties of the Anderson model (and long-range generalizations). The density of states at the edge of the spectrum displays a universal linear tail, different from the Lifshitz tails. The eigenvectors are instead very similar in the two matrix ensembles. We focus on the ground state (soft mode), which represents the epicenter of avalanche instabilities. We expect it to be localized in all finite dimensions, and make a clear connection between its localization length and the Larkin length of the depinning model. In the fully connected model, we show that the weak-strong pinning transition coincides with a peculiar localization transition of the ground state.

The influence of the strength of repulsion between particles on the thermodynamic curvature scalar $R$ for the fluid and solid states is investigated for particles interacting with the inverse power ($r^{-n}$) potential, where $r$ is the pair separation and $1/n$ is the softness. Exact results are obtained for $R$ in certain limiting cases, and the $R$ behavior determined for the systems in the fluid and solid phases. It is found that in such systems the thermodynamic curvature can be positive for very soft particles, negative for steeply repulsive (or large $n$) particles across almost the entire density range, and can change sign between negative and positive at a certain density. The relationship between $R$ and the form of the interaction potential is more complex than previously suggested, and it may be that $R$ is an indicator of the relative importance of energy and entropy contributions to the thermodynamic properties of the system.

We consider a reaction-diffusion system with linear, stochastic activator-inhibitor kinetics where the time evolution of concentration of a species at any spatial location depends on the relative average concentration of its neighbors. This self-regulating nature of kinetics brings in spatial correlation between the activator and the inhibitor. An interplay of this correlation in kinetics and disparity of diffusivities of the two species leads to symmetry breaking non-equilibrium transition resulting in stationary pattern formation. The role of initial noise strength and the linear reaction terms has been analyzed for pattern selection.

We study the collective dynamics of self-propelled rods in an inhomogeneous motility field. At the interface between two regions of constant but different motility, a smectic rod layer is spontaneously created through aligning interactions between the active rods, reminiscent of an artificial, semipermeable membrane. This “active membrane” engulfes rods which are locally trapped in low-motility regions and thereby further enhances the trapping efficiency by self-organization, an effect which we call “self-encapsulation.” Our results are gained by computer simulations of self-propelled rod models confined on a two-dimensional planar or spherical surface with a stepwise constant motility field, but the phenomenon should be observable in any geometry with sufficiently large spatial inhomogeneity. We also discuss possibilities to verify our predictions of active-membrane formation in experiments of self-propelled colloidal rods and vibrated granular matter.

The low-temperature thermal properties of glasses are anomalous with respect to those of crystals. These thermal anomalies indicate that the low-frequency vibrational properties of glasses differ from those of crystals. Recent studies revealed that, in the simplest model of glasses, i.e., the harmonic potential system, phonon modes coexist with soft localized modes in the low-frequency (continuum) limit. However, the nature of low-frequency vibrational modes of more realistic models is still controversial. In the present work, we study the Lennard-Jones (LJ) system using large-scale molecular-dynamics (MD) simulation and establish that the vibrational property of the LJ glass converges to coexistence of the phonon modes and the soft localized modes in the continuum limit as in the case of the harmonic potential system. Importantly, we find that the low-frequency vibrations are rather sensitive to the numerical scheme of potential truncation, which is usually implemented in the MD simulation, and this is the reason why contradictory arguments have been reported by previous works. We also discuss the physical origin of this sensitiveness by means of a linear stability analysis.

We systematically investigate the role of different swimming patterns on the concentration distribution of bacterial suspensions confined between two flat walls, by considering wild-type motility Escherichia coli and Pseudomonas aeruginosa, which perform Run and Tumble and Run and Reverse patterns, respectively. The experiments count motile bacteria at different distances from the bottom wall. In agreement with previous studies, an accumulation of motile bacteria close to the walls is observed. Different wall separations, ranging from 100 to $250\mu \mathrm{m}$, are tested. The concentration profiles result to be independent on the motility pattern and on the walls' separation. These results are confirmed by numerical simulations, based on a collection of self-propelled dumbbells-like particles interacting only through steric interactions. The good agreement with the simulations suggests that the behavior of the investigated bacterial suspensions is determined mainly by steric collisions and self-propulsion, as well as hydrodynamic interactions.

This paper is about the kinetic equation for gas-particle flows, in particular its well-posedness and realizability and its relationship to the generalized Langevin model (GLM) probability density function (PDF) equation. Previous analyses, e.g. [J.-P. Minier and C. Profeta, Phys. Rev. E 92, 053020 (2015)], have concluded that this kinetic equation is ill posed, that in particular it has the properties of a backward heat equation, and as a consequence, its solution will in the course of time exhibit finite-time singularities. We show that this conclusion is fundamentally flawed because it ignores the coupling between the phase space variables in the kinetic equation and the time and particle inertia dependence of the phase space diffusion tensor. This contributes an extra $\mathit{positive}$ diffusion that always outweighs the $\mathit{negative}$ diffusion associated with the dispersion along one of the principal axes of the phase space diffusion tensor. This is confirmed by a numerical evaluation of analytic solutions of these $\mathit{positive}$ and $\mathit{negative}$ contributions to the particle diffusion coefficient along this principal axis. We also examine other erroneous claims and assumptions made in previous studies that demonstrate the apparent superiority of the GLM PDF approach over the kinetic approach. In so doing, we have drawn attention to the limitations of the GLM approach, which these studies have ignored or not properly considered, to give a more balanced appraisal of the benefits of both PDF approaches.

We present a study of the stopping power of plasmas using two main approaches: the collisional (scattering theory) and the dielectric formalisms. In the former case, we use a semiclassical method based on quantum scattering theory. In the latter case, we use the full description given by the extension of the Lindhard dielectric function for plasmas of all degeneracies. We compare these two theories and show that the dielectric formalism has limitations when it is used for slow heavy ions or atoms in dense plasmas. We present a study of these limitations and show the regimes where the dielectric formalism can be used, with appropriate corrections to include the usual quantum and classical limits. On the other hand, the semiclassical method shows the correct behavior for all plasma conditions and projectile velocity and charge. We consider different models for the ion charge distributions, including bare and dressed ions as well as neutral atoms.

High-density, collisionally pumped plasma-based soft-x-ray lasers have recently delivered hundreds of femtosecond pulses, breaking the longstanding barrier of one picosecond. To pump these amplifiers an intense infrared pulse must propagate focused throughout all the length of the amplifier, which spans several Rayleigh lengths. However, strong nonlinear effects hinder the propagation of the laser beam. The use of a plasma waveguide allows us to overcome these drawbacks provided the hydrodynamic processes that dominate the creation and posterior evolution of the waveguide are controlled and optimized. In this paper we present experimental measurements of the radial density profile and transmittance of such waveguide, and we compare them with numerical calculations using hydrodynamic and particle-in-cell codes. Controlling the properties (electron density value and radial gradient) of the waveguide with the help of numerical codes promises the delivery of ultrashort (tens of femtoseconds), coherent soft-x-ray pulses.

Nonequilibrium systems commonly exhibit Lévy noise. This means that the distribution for the size of the Brownian fluctuations has a “fat” power-law tail. Large Brownian kicks are then more common as compared to the ordinary Gaussian distribution. We consider a two-state system, i.e., two wells and a barrier in between. The barrier is sufficiently high for a barrier crossing to be a rare event. When the noise is Lévy, we do not get a Boltzmann distribution between the two wells. Instead we get a situation where the distribution between the two wells also depends on the height of the barrier that is in between. Ordinarily, a catalyst, by lowering the barrier between two states, speeds up the relaxation to an equilibrium, but does not change the equilibrium distribution. In an environment with Lévy noise, on the other hand, we have the possibility of epicatalysis, i.e., a catalyst effectively altering the distribution between two states through the changing of the barrier height. After deriving formulas to quantitatively describe this effect, we discuss how this idea may apply in nuclear reactors and in the biochemistry of a living cell.

The thermoelectric properties of a three-terminal quantum spin Hall (QSH) sample are examined. The inherent helicity of the QSH sample helps to generate a large charge power efficiently. Along with charge the system can be designed to work as a highly efficient spin heat engine too. The advantage of a helical over a chiral sample is that, while a multiterminal quantum Hall sample can only work as a quantum heat engine due to broken time reversal (TR) symmetry, a multiterminal QSH system can work effectively as both a charge or spin heat engine and as a charge or spin refrigerator as the TR symmetry is preserved.

We address the detailed study of the energy current and its components, heat and work, in the boundary-driven one-dimensional $XXZ$ quantum model. We carry out the investigation by considering two different approaches present in the literature. First, we take the repeated interaction scheme and derive the expressions for the currents of heat and work, exchanged between system and baths. Then we perform the derivation of the energy current by means of a Lindblad master equation together with a continuity equation, another approach which is recurrently used. A comparison between the obtained expressions allows us to show the consistency of both approaches, and, in the latter expression derived from the Lindblad equation, it allows us to split the energy, which comes from the baths to the system, into heat and work. The recognition of work in the process, which is recurrently ignored in studies of transport, enables us to understand thermodynamical aspects and to solve some imbroglios in the physics behind the energy current in the $XXZ$ spin chain.