2024-08-15 | | Total: 5

Distribution functions of collisionless systems are known to show non-thermal power law tails. Interestingly, collisionless plasmas in various physical scenarios, (e.g., the ion population of the solar wind) feature a $v^{-5}$ tail in the velocity ($v$) distribution, whose origin has been a long-standing mystery. We show this power law tail to be a natural outcome of the self-consistent collisionless relaxation of driven electrostatic plasmas. We perform a quasilinear analysis of the perturbed Vlasov-Poisson equations to show that the coarse-grained mean distribution function (DF), $f_0$, follows a quasilinear diffusion equation with a diffusion coefficient $D(v)$ that depends on $v$ through the plasma dielectric constant. If the plasma is isotropically forced on scales much larger than the Debye length with a white noise-like electric field, then $D(v)\sim v^4$ for $\sigma<v<\omega_{\mathrm{P}}/k$, with $\sigma$ the thermal velocity, $\omega_{\mathrm{P}}$ the plasma frequency and $k$ the maximum wavenumber of the perturbation; the corresponding $f_0$, in the quasi-steady state, develops a $v^{-\left(d+2\right)}$ tail in $d$ dimensions ($v^{-5}$ tail in 3D), while the energy ($E$) distribution develops an $E^{-2}$ tail irrespective of the dimensionality of space. Any redness of the noise only alters the scaling in the high $v$ end. Non-resonant particles moving slower than the phase-velocity of the plasma waves ($\omega_{\mathrm{P}}/k$) experience a Debye-screened electric field, and significantly less (power law suppressed) acceleration than the near-resonant particles. Thus, a Maxwellian DF develops a power law tail. The Maxwellian core ($v<\sigma$) eventually also heats up, but over a much longer timescale than that over which the tail forms. We definitively show that self-consistency (ignored in test-particle treatments) is crucial for the development of the universal $v^{-5}$ tail.

Correlations of the displacement field in a two-dimensional model colloidal liquid is studied numerically and analytically. By calculating the displacement correlations and the shear strain correlations from the numerical data of particle simulations, the displacement field is shown to have nontrivial correlations, even in liquids that are only slightly glassy with the area fraction as low as 0.5. It is suggested analytically and demonstrated numerically that the displacement correlations are more informative than the shear correlations: the former behaves logarithmically with regard to the spatial distance at shorter scales, while the corresponding information is missing from the shear correlations. The logarithmic behavior of the displacement correlations is interpreted as manifesting a long-lived aspect of the cage effect.

The hydrodynamics of crystals with vacancies is developed on the basis of local-equilibrium thermodynamics, where the chemical potential of vacancies plays a key role together with a constraint relating the concentration of vacancies to the density of mass and the strain tensor. The microscopic foundations are established, leading to Green-Kubo and Einstein-Helfand formulas for the transport coefficients, including the vacancy conductivities and the coefficients of vacancy thermodiffusion. As a consequence of having introduced the chemical potential of vacancies, a relationship is obtained between the conductivities and the Fickian diffusion coefficients for the vacancies. The macroscopic equations are linearized around equilibrium to deduce the dispersion relations of the eight hydrodynamic modes. The theoretical predictions are confirmed by numerical simulations of the hard-sphere crystal with vacancies. The study explicitly shows that the eighth hydrodynamic mode of nonperfect monatomic crystals is indeed a mode of vacancy diffusion.

In the microcanonical thermal pure quantum (mTPQ) method, the canonical ensemble is derived using Taylor series expansions. We prove that the truncation error decreases exponentially with system size when the effective temperature of the mTPQ state is smaller than the target temperature, and otherwise, the error remains constant. We also show the discipline to set the mTPQ parameter by considering the trade-off between the error and the numerical cost.

Hydrodynamic turbulence exhibits nonequilibrium behaviour with $k^{-5/3}$ energy spectrum, and equilibrium behaviour with $k^{d-1}$ energy spectrum and zero viscosity, where $d$ is the space dimension. Using recursive renormalization group {in Craya-Herring basis}, we show that the nonequilibrium solution is valid only for $d < 6$, whereas equilibrium solution with zero viscosity is the only solution for $d>6$. Thus, $d=6$ is the critical dimension for hydrodynamic turbulence. In addition, we show that the energy flux changes sign from positive to negative near $d=2.15$. We also compute the energy flux and Kolmogorov's constants for various $d$'s, and observe that our results are in good agreement with past numerical results.