Flow-driven transport of soft particles in porous media is ubiquitous in many natural and engineering processes, such as the gel treatment for enhanced oil recovery. In many of these processes, injected deformable particles block the pores and thus increase the overall pressure drop and reduce the permeability of the particle-resided region. The change of macroscopic properties (e.g., pressure drop and permeability) is an important indicator of the system performance, yet sometimes impossible to be measured. Therefore, it is desirable to correlate these macroscopic properties with the measurable or controllable properties. In this work, we study flow-driven transport of soft particles in porous media using a generalized capillary bundle model. By modeling a homogeneous porous medium as parallel capillaries along the flow direction with periodically distributed constrictions, we first build a governing differential equation for pressure. Solving this equation gives a quantitative correlation between the total pressure drop and measurable parameters, including concentration and stiffness of particles, size ratio of particle to pore throat, and flow rate. The resultant permeability reduction is also obtained. Our results show that the total pressure drop and permeability reduction are both exponentially dependent on the particle concentration and the size ratio of particles to pore throat. With no more than two fitting parameters, our model shows excellent agreements with several reported experiments. The work not only sheds light on understanding transport of soft particles in porous media but also provides important guidance for choosing the optimal parameters in the relevant industrial processes.

Laser-driven ion acceleration has attracted global interest for its potential towards the development of a new generation of compact, low-cost accelerators. Remarkable advances have been seen in recent years with a substantial proton energy increase in experiments, when nanometer-scale ultrathin foil targets and high-contrast intense lasers are applied. However, the exact acceleration dynamics and particularly the ion energy scaling laws in this novel regime are complex and still unclear. Here, we derive a scaling law for the attainable maximum ion energy from such laser-irradiated nanometer-scale foils based on analytical theory and multidimensional particle-in-cell simulations, and further show that this scaling law can be used to accurately describe experimental data over a large range of laser and target parameters on different facilities. This provides crucial references for parameter design and experimentation of the future laser devices towards various potential applications.

Recently, asymptotic scaling behavior of the compressible randomly forced Navier-Stokes equation has been analyzed with the use of field-theoretic renormalization group near four dimensions [Phys. Rev. E 95, 033120 (2017)]. Two infrared stable nontrivial asymptotic scaling patterns have been found and their parameters determined in the one-loop approximation. Here, it is pointed out that the asymptotic scaling behavior predicted in this way may not be realized in physical fluid systems. This is a consequence of restrictions on viscosities imposed by the energy balance equation which the one-loop fixed-point value of the relative viscosity fails to meet.

We explore quantum phase transitions using two probes of quantum chaos: out-of-time-order correlators (OTOCs) and the $r$-parameter obtained from the level spacing statistics. In particular, we address $p$-spin models associated with quantum annealing or reverse annealing. Quantum annealing triggers first-order or second-order phase transitions, which is crucial for the performance of quantum devices. We find that the time-averaging OTOCs for the ground state and the average $r$-parameter change behavior around the corresponding transition points, diagnosing the phase transition. Furthermore, they can identify the order (first or second) of the phase transition by their behavior at the quantum transition point, which changes abruptly (smoothly) in the case of first-order (second-order) phase transitions.

Rayleigh-Taylor instability (RTI) has become a powerful tool for determining the mechanical properties of materials under extreme conditions. In this paper, we first present the exact and approximate linear dispersion relations for RTI in viscoelastic materials based on the Maxwell and Kelvin-Voigt models. The approximate dispersion relation produces good predictions of growth rates in comparison with the exact one. The motion of the interface in Maxwell flow is mainly controlled by viscosity and elasticity dominates this behavior in Kelvin-Voigt flow. Since elasticity plays a distinct role from viscosity, cutoff wavelengths arise only in Kelvin-Voigt flow. The variation of the maximum growth rates and their corresponding wave numbers are also carefully studied. For both types of materials, viscosity suppresses the growth of instability, while elasticity speeds it up. This is at odds with the well-known understanding that elasticity suppresses hydrodynamic instabilities. The dependence of the maximum growth rate on slab thickness is also investigated for RTI in both types of flow, since the metal slab as a pusher has been extensively employed in high-energy-density physics. The model presented here allows study of more realistic situations by considering convergent effects and shock wave interactions, for the traditional potential flow theory is not suitable. To summary, it is able to provide guidances for future experimental designs for studies of materials under high strain and high strain rate conditions, as well as allow us to study RTI theoretically in more complicated conditions.

We conducted numerical experiments to investigate the mixing of stratified suspensions containing different types of particles. We used a point-force two-way coupling method. We studied the mixing behavior of stratified suspensions and we discovered two types of mixing: microscopic (individual-particle-level) and macroscopic (vessel-scale) collective mixing. In addition, we examined the vertical mixing speed of the stratified suspension. We used a simple theoretical model to analyze the fingering settling velocity. Then we introduced a nondimensional number representing the difference in collectivities of the upper and lower suspensions while accounting for particle terminal velocities. We discovered that the proposed nondimensional parameter has a negative sign that distinguishes the mixing form of only microscopic individual-particle-level mixing and a positive value that predicts the speed of macroscopic collective mixing of stratified suspensions.

Multistability is one of the most important phenomena in dynamical systems, e.g., bistability enables the implementation of logic gates and therefore computation. Recently multistability has attracted a greatly renewed interest related to memristors and graphene structures, to name only a few. We investigate tristability in velocity dynamics of a Brownian particle subjected to a tilted periodic potential. It is demonstrated that the origin of this effect is attributed to the arcsine law for the velocity dynamics at the zero temperature limit. We analyze the impact of thermal fluctuations and construct the phase diagram for the stability of the velocity dynamics. It suggests an efficient strategy to control the multistability by changing solely the force acting on the particle or temperature of the system. Our findings for the paradigmatic model of nonequilibrium statistical physics apply to, inter alia, Brownian motors, Josephson junctions, cold atoms dwelling in optical lattices, and colloidal systems.

Periodic environmental changes are commonly observed in nature from the amount of daylight to seasonal temperature. These changes usually affect individuals' death or birth rates, dragging the system from its previous stable states. When the fluctuation of abundance is amplified due to such changes, extinction of species may be accelerated. To see this effect, we examine how the abundance and the mean time to extinction respond to the periodic environmental changes. We consider a population wherein two species coexist together implemented by three rules—birth, spontaneous death, and death from competitions. As the interspecific interaction strength is varied, we observe the resonance behavior in both fluctuations of abundances and the mean time to extinction. Our result suggests that neither too high nor too low competition rates make the system more susceptible to environmental changes.

We study the effects of higher-order antinematic interactions on the critical behavior of the antiferromagnetic (AFM) $XY$ model on a triangular lattice using Monte Carlo simulations. The parameter $q$ of the generalized antinematic (ANq) interaction is found to have a pronounced effect on the phase diagram topology by inducing new quasi-long-range ordered phases due to competition with the conventional AFM interaction as well as geometrical frustration. For values of $q$ divisible by 3, the ground-state competition between the two interactions results in a frustrated canted AFM phase appearing at low temperatures wedged between the AFM and ANq phases. For $q$ nondivisible by 3, with the increase of $q$ one can observe the evolution of the phase diagram topology featuring two $(q=2)$, three $(q=4,5)$, and four $(q\ge 7)$ ordered phases. In addition to the two phases previously found for $q=2$, the first new phase with solely AFM ordering arises for $q=4$ in the limit of strong AFM coupling and higher temperatures by separating from the phase with the coexisting AFM and ANq orderings. For $q=7$, another phase with AFM ordering but multimodal spin distribution in each sublattice appears at intermediate temperatures. All these algebraic phases also display standard and generalized chiral long-range orderings, which decouple at higher temperatures in the regime of dominant ANq (AFM) interaction for $q\ge 4$ $(q\ge 7)$ preserving only the generalized (standard) chiral ordering.

We consider conditions under which an isolated quantum system approaches a microcanonical equilibrium state. A key component is the eigenstate thermalization hypothesis, which proposes that all energy eigenstates appear thermal. We introduce a weaker version of this requirement, applying only to the average distinguishability of eigenstates from the thermal state, and investigate its necessity and sufficiency for thermalization.

We model $Q$-switched pulses in passively mode-locked lasers using the cubic-quintic complex Ginzburg-Landau equation (CGLE). We show that a wide set of parameters of this equation leads to $Q$-switched pulses of triangular shape that consist of a periodic sequence of evolving dissipative solitons. Bifurcation diagrams demonstrating various transformations of these pulses are calculated in terms of five major parameters of the CGLE. The diagrams show period doubling transformations as well as the transition to a chaotic evolution of the $Q$-switched pulses.

The dynamics of matter waves in the atomic to molecular condensate transition with a time-modulated atomic scattering length is investigated. Both the cases of rapid and slow modulations are studied. In the case of rapid modulations, the average over oscillations for the system is derived. The corresponding conditions for dynamical suppression of the association of atoms into the molecular field, or of second-harmonic generation in nonlinear optical systems, are obtained. For the case of slow modulations, we find resonant enhancement in the molecular field. We then illustrate chaos in the atomic-molecular BEC system. We suggest a sequential application of the two types of modulations, slow and rapid, when producing molecules.

The spatial correlations of entangled polymer dynamics are examined by molecular dynamics simulations and neutron spin-echo spectroscopy. Due to the soft nature of topological constraints, the initial spatial decays of intermediate scattering functions of entangled chains are, to the first approximation, surprisingly similar to those of an unentangled system in the functional forms. However, entanglements reveal themselves as a long tail in the reciprocal-space correlations, implying a weak but persistent dynamic localization in real space. Comparison with a number of existing theoretical models of entangled polymers suggests that they cannot fully describe the spatial correlations revealed by simulations and experiments. In particular, the strict one-dimensional diffusion idea of the original tube model is shown to be flawed. The dynamic spatial correlation analysis demonstrated in this work provides a useful tool for interrogating the dynamics of entangled polymers. Lastly, the failure of the investigated models to even qualitatively predict the spatial correlations of collective single-chain density fluctuations points to a possible critical role of incompressibility in polymer melt dynamics.

The aim of this paper is to investigate the response of a homogeneously aligned nematic nanosized hollow cavity (HANNHC) confined between two charged horizontal coaxial cylinders and subjected to both a radially applied electrostatic field $\mathbf{\text{E}}$, arising from the surface charge density $\kappa$ and the temperature gradient $\mathbf{\nabla }T$ set between these cylinders. This was done within the framework of an extension of the classical Ericksen-Leslie theory, supplemented by thermomechanical correction of the shear stress and Rayleigh dissipation function, as well as taking into account the entropy balance equation. The physical mechanism responsible for the excitation of the hydrodynamic flow in the HANNHC is based on the interaction of the director and temperature gradients and the static electric field. Calculations show that under the influence of both the $\mathbf{\nabla }T$ and $\mathbf{\text{E}}$, a stationary flow $u^{\mathrm{st}}$ is excited in the HANNHC in the horizontal direction. It is shown that the electric force enforced by the flexoelectric polarization plays a crucial role in the excitation of $u^{\mathrm{st}}$ between these cylinders.

Outcrops of columnar joints made of solidified lava flows are often covered by semiordered polygonal cracks. The polygon diameters are fairly uniform at each outcrop, but their shapes largely vary in the number of sides and internal angles. Herein, we unveil that the statistical variation in the polygon shape follows an extreme value distribution class: the Gumbel distribution. The Gumbel law was found to hold for different columnar joints, regardless of the locality, lithologic composition, and typical diameter. A common distribution for columnar joints implies a universal class that may integrate the polygonal crack networks observed on the surface of various fractured brittle materials.

In this paper, direct numerical simulations of particle-laden homogeneous isotropic turbulence are performed using lattice Boltzmann method incorporating interpolated bounce-back scheme. Four different particle-fluid density ratios are considered to explore how particles with different particle-fluid density ratios respond to the turbulence. Overall particle dynamics in the homogeneous isotropic turbulence such as the Lagrangian statistics of single particle and the preferential concentration of particles are investigated. Results show that particle acceleration and angular acceleration are more intermittent than velocity and angular velocity for finite-size particles with different particle-fluid density ratios. The preferential concentration of particles is investigated using radial distribution function and Voronoï tessellation, and the preferential concentration is more profound for particles with two intermediate particle-fluid density ratios. The Voronoï analysis indicates that the distribution of Voronoï cells satisfy the log-normal distribution better than the gamma distribution. The mechanism of preferential concentration is analyzed using the sweep-stick mechanism and drift mechanism. Results show that although a higher probability of having particles located near the sticky points is found, the sticky mechanism is very weak for large density ratios. The particle clustering is then found to be better qualitatively described by the drift mechanism.

We present a simple approach to high-accuracy calculations of critical properties for the three-dimensional Ising model, without prior knowledge of the critical temperature. The iterative method uses a modified block-spin transformation with a tunable parameter to improve convergence in the Monte Carlo renormalization group trajectory. We found experimentally that the iterative method enables the calculation of the critical temperature simultaneously with a critical exponent.

Wigner crystallization in macroscopic systems provides a simple avenue for studying classical many-body physics. The Coulomb repulsion of millimeter-sized conducting spheres rolling on a high potential surface enables the self-assembly of ordered structures. Prior attempts at formulating pairwise force models for Wigner islands have resulted in different functional forms with no apparent consensus. In this study, we explore fundamental particle interactions using the finite element method to solve the electrostatic problem for one-, two-, and many-sphere configurations. Analysis of charge and force shows three interaction regimes parametrized by the packing fraction, which includes many-body force behavior for compact geometries. We find that the force relationship between spheres is accurately represented by an aligned dipole model over all configurations. Both the finite element method and aligned dipole models were validated in absolute terms using an experimental setup featuring a gravity potential. Future dynamic studies, which rely on accurate force relationships, can now be treated.

The specific-heat capacity $c_v$ of glass formers undergoes a hysteresis when subjected to a cooling-heating cycle, with a larger $c_v$ and a more pronounced hysteresis for fragile glasses than for strong ones. Here we show that these experimental features, including the unusually large magnitude of $c_v$ of fragile glasses, are well reproduced by kinetic Monte Carlo and equilibrium study of a distinguishable particle lattice model incorporating a two-state picture of particle interactions. The large $c_v$ in fragile glasses is caused by a dramatic transfer of probabilistic weight from high-energy particle interactions to low-energy ones as temperature decreases.