Wall turbulence has a sublayer where one-point statistics, e.g., the mean velocity and the variances of some velocity fluctuations, vary logarithmically with the distance from the wall. This logarithmic scaling is found here for two-point statistics or specifically two-point cumulants of those fluctuations by means of experiments in a wind tunnel. As for corresponding statistics of the rate of the energy dissipation, the scaling is found to be not logarithmic. We reproduce these scaling laws with some mathematics and also with a model of energy-containing eddies that are attached to the wall.

We investigate the energy-loss characteristics of an ion in warm dense matter (WDM) and dense plasmas concentrating on the influence of electronic correlations. The basis for our analysis is a recently developed ab initio quantum Monte Carlo– (QMC) based machine learning representation of the static local field correction (LFC) [Dornheim et al., J. Chem. Phys. 151, 194104 (2019)], which provides an accurate description of the dynamical density response function of the electron gas at the considered parameters. We focus on the polarization-induced stopping power due to free electrons, the friction function, and the straggling rate. In addition, we compute the friction coefficient which constitutes a key quantity for the adequate Langevin dynamics simulation of ions. Considering typical experimental WDM parameters with partially degenerate electrons, we find that the friction coefficient is of the order of $\gamma /{\omega }_{pi}=0.01$, where ${\omega }_{pi}$ is the ionic plasma frequency. This analysis is performed by comparing QMC-based data to results from the random-phase approximation (RPA), the Mermin dielectric function, and the Singwi-Tosi-Land-Sjölander (STLS) approximation. It is revealed that the widely used relaxation time approximation (Mermin dielectric function) has severe limitations regarding the description of the energy loss of ions in a correlated partially degenerate electrons gas. Moreover, by comparing QMC-based data with the results obtained using STLS, we find that the ion energy-loss properties are not sensitive to the inaccuracy of the static local field correction (LFC) at large wave numbers, $k/k_F>2$ (with $k_F$ being the Fermi wave number), but that a correct description of the static LFC at $k/k_F\lesssim 1.5$ is important.

Cotranslational folding is one of the most important features of protein folding in vivo. Although many studies have shown that the folding pathways of cotranslational folding are different from free folding in vitro, the detailed mechanism of folding dynamics is lacking. Here we combine all-atom molecular simulations with an ideal ribosome tunnel model to investigate the cotranslational folding of villin headpiece subdomain HP35. By comparing the folding dynamics between cotranslational folding and free folding, we found that cotranslational folding tends to fold along the pathway that is easier to fold into native state in the latter. In addition, the roles of the ribosome tunnel and sequential folding are analyzed separately. Our results show that the ribosome can prevent the untimely folding of the $C$ segment of HP35 to reduce the non-native interactions, while the translation speed can regulate the amounts of native and non-native interactions and the balance between them. Overall, these results give insights into the general mechanisms of cotranslational protein folding.

We perform systematic particle dynamics simulations of granular flows composed of breakable particles in a 2D rotating drum to investigate the evolution of the mean particle size and specific surface as a function of system parameters such as drum size, rotation speed, filling degree, and particle shape and size. The specific surface increases at a nearly constant rate up to a point where particle breakage begins to slow down. The rates of particle breakage for all values of system parameters are found to collapse on a master curve when the times are scaled by the characteristic time defined in the linear regime. We determine the characteristic time as a function of all system parameters, and we show that the rate of particle breakage can be expressed as a linear function of a general scaling parameter that incorporates all our system parameters. This scaling behavior provides a general framework for the upscaling of drum grinding process from laboratory to industrial scale.

An integral expression for the translational velocity of a perfectly slipping spherical particle under a time-dependent applied force in unsteady Stokes flow is derived. For example, when the ratio of particle density to fluid density is small, our analysis pertains to an inviscid bubble in a viscous fluid. We determine an explicit form of the particle velocity under an impulsive force, wherefrom the velocity autocorrelation function and mean-squared displacement of a perfectly slipping sphere undergoing Brownian motion are obtained. The above results are contrasted against the time-dependent diffusion of a rigid sphere with no hydrodynamic slip. Finally, the thermal force power spectral density is analytically calculated for a diffusing sphere with arbitrary slip length. We suggest this quantity to be suitable to infer slip length from the measurement of nondiffusive Brownian motion.

We develop a real-space renormalization-group (RSRG) scheme by appropriately inserting the long-range hopping $t\sim r^{-\alpha }$ with nearest-neighbor interaction to study the entanglement entropy and maximum block size for the many-body localization (MBL) transition. We show that for $\alpha <2$ there exists a localization transition with renormalized disorder that depends logarithmically on the finite size of the system. The transition observed for $\alpha >2$ does not need a rescaling in disorder strength. Most important, we find that even though the MBL transition for $\alpha >2$ falls in the same universality class as that of the short-range models, the transition for $\alpha <2$ belongs to a different universality class. Because of the intrinsic nature of the RSRG flow toward delocalization, MBL phase for $\alpha >2$ might suffer an instability in the thermodynamic limit while the underlying systems support algebraic localization. Moreover, we verify these findings by inserting microscopic details to the RSRG scheme where we additionally find a more appropriate rescaling function for disorder strength; we indeed uncover a power-law scaling with a logarithmic correction and a distinctly different stretched exponential scaling for $\alpha <2$ and $\alpha >2$, respectively, by analyzing system with finite size. This finding further suggests that microscopic RSRG scheme is able to give a hint of instability of the MBL phase for $\alpha >2$ even considering systems of finite size.

The physics of supercritical states is understood to a much lesser degree compared to subcritical liquids. Carbon dioxide, in particular, has been intensely studied, yet little is known about the supercritical part of its phase diagram. Here, we combine neutron scattering experiments and molecular dynamics simulations and demonstrate the structural crossover at the Frenkel line. The crossover is seen at pressures as high as 14 times the critical pressure and is evidenced by changes of the main features of the structure factor and pair distribution functions.

In the past few years, the exponential expansion analysis of time autocorrelation functions has provided profound insight into the leading microscopic processes driving the atomic-scale dynamics and has made it possible to highlight the presence and the role of various relaxation channels through which the fundamental correlation functions decay with time. Here we apply this method to the determination of the full time dependence of a correlation function $c\left(t\right)$ in a quantum system at nonzero temperature, by making explicit its relationship with its Kubo transform $c_{\mathrm{K}}\left(t\right)$, which in some cases can be approximately computed with the presently available quantum simulation techniques. We obtain an exact expression for $c\left(t\right)$ in terms of the exponential modes that describe the time behavior of $c_{\mathrm{K}}\left(t\right)$. The relative importance of the various modes in determining the overall shape of $c\left(t\right)$ can then be studied in detail. This work extends to the full time domain the results of a previous paper [Guarini et al., Phys. Rev. Lett. 123, 135301 (2019)], in which we employed the same method to calculate the zero time value of the velocity autocorrelation function, to obtain a microscopic description of the quantum mean kinetic energy in a fluid. In particular, we show that the decay constants and the frequencies of the dominant microscopic modes of $c\left(t\right)$ are the same as those of $c_{\mathrm{K}}\left(t\right)$, but the dynamics of the quantum system also contains an additional term decaying on a time scale determined solely by temperature of the system.

We consider a one-dimensional Brownian motion of fixed duration $T$. Using a path-integral technique, we compute exactly the probability distribution of the difference $\tau =t_{min}-t_{max}$ between the time $t_{min}$ of the global minimum and the time $t_{max}$ of the global maximum. We extend this result to a Brownian bridge, i.e., a periodic Brownian motion of period $T$. In both cases, we compute analytically the first few moments of $\tau$, as well as the covariance of $t_{max}$ and $t_{min}$, showing that these times are anticorrelated. We demonstrate that the distribution of $\tau$ for Brownian motion is valid for discrete-time random walks with $n$ steps and with a finite jump variance, in the limit $n\to \infty$. In the case of Lévy flights, which have a divergent jump variance, we numerically verify that the distribution of $\tau$ differs from the Brownian case. For random walks with continuous and symmetric jumps we numerically verify that the probability of the event “$\tau =n$” is exactly $1/\left(2n\right)$ for any finite $n$, independently of the jump distribution. Our results can be also applied to describe the distance between the maximal and minimal height of $(1+1)$-dimensional stationary-state Kardar-Parisi-Zhang interfaces growing over a substrate of finite size $L$. Our findings are confirmed by numerical simulations. Some of these results have been announced in a recent Letter [Phys. Rev. Lett. 123, 200201 (2019)].

We report on a fundamental role of a non-normalized formal steady state, i.e., an infinite invariant density, in a semi-Markov process where the state is determined by the interevent time of successive renewals. The state describes certain observables found in models of anomalous diffusion, e.g., the velocity in the generalized Lévy walk model and the energy of a particle in the trap model. In our model, the interevent-time distribution follows a fat-tailed distribution, which makes the state value more likely to be zero because long interevent times imply small state values. We find two scaling laws describing the density for the state value, which accumulates in the vicinity of zero in the long-time limit. These laws provide universal behaviors in the accumulation process and give the exact expression of the infinite invariant density. Moreover, we provide two distributional limit theorems for time-averaged observables in these nonstationary processes. We show that the infinite invariant density plays an important role in determining the distribution of time averages.

Nonequilibrium systems exhibit particle-type solutions. Oscillons are one of the best-known localized states of systems with time-dependent forcing or parametrically driven systems. We investigate the transition from nonradiative to radiative oscillons in the parametrically driven sine-Gordon model in two spatial dimensions. The bifurcation takes place when the strength of the forcing (frequency) increases (decreases) above a certain threshold. As a result of this transition, the oscillon emits radially symmetric evanescent waves. Numerically, we provide the phase diagram and show the supercritical nature of this transition. For small oscillations, based on the amplitude equation approach, the sine-Gordon equation with time-dependent forcing is transformed into the parametrically driven damped nonlinear Schrödinger model in two spatial dimensions. This amplitude equation exhibits a transition between nonradiative to radiative localized structures, consistently. Both models show quite good agreement.

Ring vaccination is a mitigation strategy that consists in seeking and vaccinating the contacts of a sick patient, in order to provide immunization and halt the spread of disease. We study an extension of the susceptible-infected-recovered (SIR) epidemic model with ring vaccination in complex and spatial networks. Previously, a correspondence between this model and a link percolation process has been established, however, this is only valid in complex networks. Here, we propose that the SIR model with ring vaccination is equivalent to a mixed percolation process of links and nodes, which offers a more complete description of the process. We verify that this approach is valid in both complex and spatial networks, the latter being built according to the Waxman model. This model establishes a distance-dependent cost of connection between individuals arranged in a square lattice. We determine the epidemic-free regions in a phase diagram based on the wiring cost and the parameters of the epidemic model (vaccination and infection probabilities and recovery time). In addition, we find that for long recovery times this model maps into a pure node percolation process, in contrast to the SIR model without ring vaccination, which maps into link percolation.

We present a study of the landscape structure of hard and soft spheres as a function of the packing fraction and of the energy. We find that, on approaching the jamming transition, the number of different configurations available to the system increases steeply and a hierarchical organization of the landscape emerges. We use the knowledge of the structure of the landscape to predict the values of thermodynamic observables on the edge of the transition.

We investigate the inviscid compressible flow (Euler) equations constrained by an “isentropic” equation of state (EOS), whose functional form in pressure is an arbitrary function of density alone. Under the aforementioned condition, we interrogate using symmetry methods the scale-invariance of the homentropic inviscid Euler equations. We find that under general conditions, we can reduce the inviscid Euler equations into a system of two coupled ordinary differential equations. To exemplify the utility of these results, we formulate two example scale-invariant, self-similar solutions. The first example includes a shock-free expanding bubble scenario, featuring a modified Tait EOS. The second example features the classical Noh problem, coupled to an arbitrary isentropic EOS. In this case, to satisfy the conditions set forth in the classical Noh problem, we find that the solution for the flow is given by a transcendental algebraic equation in the shocked density.

The nonequilibrium thermodynamics of an open (classical or quantum) system in strong contact with a single heat bath can be conveniently described in terms of the Hamiltonian of mean force. However, the conventional formulation is limited by the necessity to measure differences in equilibrium properties of the system-bath composite. We make use of the freedom involved in defining thermodynamic quantities, which leaves the thermodynamics unchanged, to show that the Hamiltonian of mean force can be inferred from measurements on the system alone, up to that irrelevant freedom. In doing so, we refute a key criticism expressed in the works by P. Talkner and P. Hänggi [Phys. Rev. E 94, 022143 (2016) and arXiv:1911.11660]. We also discuss the remaining part of the criticism.

We determine the caloric curves of classical self-gravitating systems at statistical equilibrium in general relativity. In the classical limit, the caloric curves of a self-gravitating gas depend on a unique parameter $\nu =GNm/R{c}^2$, called the compactness parameter, where $N$ is the particle number and $R$ the system's size. Typically, the caloric curves have the form of a double spiral. The “cold spiral,” corresponding to weakly relativistic configurations, is a generalization of the caloric curve of nonrelativistic classical self-gravitating systems. The “hot spiral,” corresponding to strongly relativistic configurations, is similar (but not identical) to the caloric curve of the ultrarelativistic self-gravitating black-body radiation. We introduce two types of normalization of energy and temperature to obtain asymptotic caloric curves describing, respectively, the cold and the hot spirals in the limit $\nu \to 0$. As the number of particles increases, the cold and the hot spirals approach each other, merge at ${\nu }_S^{\prime }=0.128$, form a loop above ${\nu }_S=0.1415$, reduce to a point at ${\nu }_{\max }=0.1764$, and finally disappear. Therefore, the double spiral shrinks when the compactness parameter $\nu$ increases, implying that general relativistic effects render the system more unstable. We discuss the nature of the gravitational collapse at low and high energies with respect to a dynamical (fast) or a thermodynamical (slow) instability. We also provide an historical account of the developments of the statistical mechanics of classical self-gravitating systems in Newtonian gravity and general relativity.

Under certain conditions, a counterintuitive behavior—an initially hotter sample freezes faster when quenched to a cold bath than an identical system initialled at a lower temperature—is known as the Mpemba effect (ME). Here we identify the existence of the ME in mean-field systems (MFS). Specifically, the thermal contact between MFS and a large thermal reservoir is built up using the microcanonical Monte Carlo algorithm. The simulation results unambiguously demonstrate that an initial hotter system undergoes the paramagnetic-ferromagnetic phase transition faster than the initial cooler one. The ME here originates from the back-reaction of the MFS system on the reservoir, which is thus an embodiment of non-Markovianness in relaxation. In addition, we confirm that the ME survives in the thermodynamic limit. And the significance of reservoir size is also explored: A smaller heat reservoir facilitates the overall relaxation process. In general, this work establishes a theoretical non-Markovian ME framework, which may shed light on widening the understanding of the mechanism behind the ME in other substances, including water.

We study the influence that collective transition phenomena have on quantum metrological protocols. The single spherical quantum spin (SQS) serves as stereotypical toy model that allows analytical insights on a mean-field level. First, we focus on equilibrium quantum criticality in the SQS and obtain the quantum Fisher information analytically, which is associated with the minimum lower bound for the precision of estimation of the parameter driving the phase transition. We compare it with the Fisher information for a specific experimental scenario where photon-counting-like measurements are employed. We find how quantum criticality and squeezing are useful resources in the metrological scenario. Second, we obtain the quantum Fisher information for the out-of-equilibrium transition in the dissipative nonequilibrium steady state and investigate how the presence of dissipation affects the parameter estimation. In this scenario, it is known that the critical point is shifted by an amount which depends on the dissipation rate. This is used here to design high precision protocols for a whole range of the transition-driving parameter in the ordered phase. In fact, for certain values of the parameter being estimated, dissipation may be used to obtain higher precision when compared to the equilibrium scenario.

Low-dimensional dynamics of large networks is the focus of many theoretical works, but controlled laboratory experiments are comparatively very few. Here, we discuss experimental observations on a mean-field coupled network of hundreds of semiconductor lasers, which collectively display effectively low-dimensional mixed mode oscillations and chaotic spiking typical of slow-fast systems. We demonstrate that such a reduced dimensionality originates from the slow-fast nature of the system and of the existence of a critical manifold of the network where most of the dynamics takes place. Experimental measurement of the bifurcation parameter for different network sizes corroborates the theory.