We prove that there is no quantum speedup when using quantum Monte Carlo integration to estimate the mean (and other moments) of analytically defined log-concave probability distributions prepared as quantum states using the Grover-Rudolph method.

Volume integrals over the radial pair-distribution function, so-called Kirkwood-Buff integrals (KBIs), play a central role in the theory of solutions by linking structural with thermodynamic information. The simplest example is the compressibility equation, a fundamental relation in statistical mechanics of fluids. Until now, KBI theory could not be applied to crystals because the integrals strongly diverge when computed in the standard way. We solve the divergence problem and generalize KBI theory to crystalline matter by using the recently proposed finite-volume theory. For crystals with harmonic interaction, we derive an analytic expression for the peak shape of the pair-distribution function at finite temperature. From this we demonstrate that the compressibility equation holds exactly in harmonic crystals.

Heat current $J$ that flows through a few typical two-dimensional nonlinear lattices is systematically studied. Each lattice consists of two identical segments that are coupled by an interface with strength $k_{\mathrm{int}}$. It is found that the two-universality-class scenario that is revealed in one-dimensional systems is still valid in the two-dimensional systems. Namely, $J$ may follow $k_{\mathrm{int}}$ in two entirely different ways, depending on whether or not the interface potential energy decays to zero. We also study the dependence of $J$ on lattice width $N_Y$ and transverse interaction strength $k_Y$. Universal power-law decay or divergence is observed. Finally, we check the equipartition theorem in the systems since it is the basis of all our theoretical analyses. Surprisingly, it holds perfectly even at the interface where there is a finite temperature jump, which makes the system far from equilibrium. However, the equipartition of potential energy, which is observed in one-dimensional systems, is no longer satisfied due to the interaction between different dimensions.

Quantum quenches to or near criticality give rise to the phenomenon of aging, manifested by glassylike dynamics at short times and far from equilibrium. The recent surge of interest in the dynamics of quantum many-body systems has rejuvenated interest in this phenomenon. Motivated by the ubiquitous long-range interactions in emerging experimental platforms, it is vital to study quantum aging in such settings. In this paper, we investigate the dynamical universality and aging in the $d$-dimensional $O\left(N\right)$ model with the long-range coupling $1/x^{d+\sigma }$ and in the mean-field limit $N\to \infty$ that allows an exact treatment. An immediate consequence of long-range coupling is the emergence of nonlinear light cones. We focus on the correlation and response functions, and identify a rich scaling behavior depending on how the corresponding space-time positions are located relative to each other, via a local light cone, and to the time of the quench via a global quench light cone. We determine the initial-slip exponent that governs the short-time dependence of two-point functions. We highlight the qualitative features of aging due to the long-range coupling, in particular in the region outside the light cones. As an important consequence of long-range coupling, the correlation function decays as $1/x^{d+\sigma }$ outside the quench light cone while increasing polynomially with the total time after quench. This is while, for short-time differences, the two-time response function “equilibrates” at all distances even outside this light cone. Our analytic findings are in excellent agreement with exact numerics, and provide a useful benchmark for modern experimental platforms with long-range interactions.

Balachandran et al. [Phys. Rev. Lett. 120, 200603 (2018)] presented a segmented XXZ spin chain with zero anisotropy in one half and a large anisotropy on the other half that gave rise to a spin current rectification which is perfect in the thermodynamic limit. Here we extend the previous study to segmented chains with interacting integrable as well as nonintegrable halves, considering even cases in which no ballistic transport can emerge in either half. We demonstrate that, also in this more general case, it is possible to obtain giant rectification when the two interacting half chains are sufficiently different. We also show that the mechanism causing this effect is the emergence of an energy gap in the excitation spectrum of the out-of-equilibrium insulating steady state in one of the two biases. Finally, we demonstrate that in the thermodynamic limit there is no perfect rectification when each of the two half chains is interacting.

A new basis has been found for the theory of self-organization of transport avalanches and jet zonal flows in L-mode tokamak plasma, the so-called “plasma staircase” [Dif-Pradalier et al., Phys. Rev. E 82, 025401(R) (2010)]. The jet zonal flows are considered as a wave packet of coupled nonlinear oscillators characterized by a complex time- and wave-number-dependent wave function; in a mean-field approximation this function is argued to obey a discrete nonlinear Schrödinger equation with subquadratic power nonlinearity. It is shown that the subquadratic power leads directly to a white Lévy noise, and to a Lévy fractional Fokker-Planck equation for radial transport of test particles (via wave-particle interactions). In a self-consistent description the avalanches, which are driven by the white Lévy noise, interact with the jet zonal flows, which form a system of semipermeable barriers to radial transport. We argue that the plasma staircase saturates at a state of marginal stability, in whose vicinity the avalanches undergo an ever-pursuing localization-delocalization transition. At the transition point, the event-size distribution of the avalanches is found to be a power law $w_{\tau }\left(\mathrm{\Delta }n\right)\sim \mathrm{\Delta }{n}^{-\tau }$, with the drop-off exponent $\tau =(\sqrt{17}+1)/2\simeq 2.56$. This value is an exact result of the self-consistent model. The edge behavior bears signatures enabling to associate it with the dynamics of a self-organized critical (SOC) state. At the same time the critical exponents, pertaining to this state, are found to be inconsistent with classic models of avalanche transport based on sand piles and their generalizations, suggesting that the coupled avalanche-jet zonal flow system operates on different organizing principles. The results obtained have been validated in a numerical simulation of the plasma staircase using flux-driven gyrokinetic code for L-mode Tore-Supra plasma.

We study the response to shear deformations of packings of long spherocylindrical particles that interact via frictional forces with friction coefficient $\mu$. The packings are produced and deformed with the help of molecular dynamics simulations combined with minimization techniques performed on a GPU. We calculate the linear shear modulus $g_{\infty }$, which is orders of magnitude larger than the modulus $g_0$ in the corresponding frictionless system. The motion of the particles responsible for these large frictional forces is governed by and increases with the length $\ell$ of the spherocylinders. One consequence of this motion is that the shear modulus $g_{\infty }$ approaches a finite value in the limit $\ell \to \infty$, even though the density of the packings vanishes, $\rho \propto {\ell }^{-2}$. By way of contrast, the frictionless modulus decreases to zero, $g_0\sim {\ell }^{-2}$, in accordance with the behavior of density. Increasing the strain beyond a value ${\gamma }_c\sim \mu$, the packing strain weakens from the large frictional to the smaller frictionless modulus when contacts saturate at the Coulomb inequality and start to slide. In this regime, sliding friction contributes a “yield stress” ${\sigma }_y=g_{\infty }{\gamma }_c$ and the stress behaves as $\sigma ={\sigma }_y+g_0\gamma$. The interplay between static and sliding friction gives rise to hysteresis in oscillatory shear simulations.

We investigate the impact of composite objects. They consist of a soft layer on top of a rigid part with a hemispherical impacting end. The coefficient of restitution ($e$) of such objects is studied systematically as a function of the mass ratio and of the nature of the materials. For rather elastic materials, the coefficient of restitution is a nonmonotonic function of the mass ratio and exhibits important variations. The dynamics of the impact can be characterized by several bounces depending on the ratios between the four timescales at play. These include the duration of contact of the rigid part with the substrate and the time for the elastic waves to travel back and forth in the soft layer. In that sense, describing these projectiles requires one to take into account both the Hertzian theory of contact and the elastic waves described by Saint-Venant's approach.

In this article, the classical Rayleigh-Taylor instability is extended to situations where the fluid is completely confined, in both the vertical and horizontal directions. This article starts with the two-dimensional (2D) viscous periodic case with finite height where the effect of adding surface tension to the interface is analyzed. This problem is simulated from a dual perspective: first, the linear stability analysis obtained when the Navier-Stokes equations are linearized and regularized in terms of density and viscosity; and second, looking at the weakly compressible version of a multiphase smoothed particle hydrodynamics (WCSPH) method. The evolution and growth rates of the different fluid variables during the linear regime of the SPH simulation are compared to the computation of the eigenvalues and eigenfunctions of the viscous version of the Rayleigh-Taylor stability (VRTI) analysis with and without surface tension. The most unstable mode, which has the maximal linear growth rate obtained with both approaches, as well as other less unstable modes with more complex structures are reported. The classical horizontally periodic (VRTI) case is now adapted to the case where two additional left and right walls are included in the problem, representing the cases where a two-phase flow is confined in a accelerated tank. This 2D case where no periodic assumptions are allowed is also solved using both techniques with tanks of different sizes and a wide range of Atwood numbers. The agreement with the linear stability analysis obtained by a Lagrangian method such as multiphase WCSPH is remarkable.

The definition of $boundary$ at the nanoscale has been a matter of dispute for years. Addressing this issue, the nonequilibrium molecular dynamics (NEMD) simulations in this work investigate the flow characteristics of a simple liquid in a single-walled carbon nanotube (SWCNT), and equilibrium molecular dynamics simulations support the range of the NEMD results. The inconsistencies in defining the flow boundary at the nanoscale are understood through the first law of thermodynamics: Local thermodynamic properties (the effects of the density distribution, pressure, viscosity, and temperature) define the boundary. We have selected different boundary positions in the CNT to demonstrate the probability of density distribution that also indicates the coexistence of multiple thermodynamic states. Altering the interaction parameters, we produce convergence between the NEMD result and the no-slip Hagen-Poiseuille assumptions. Meanwhile, the results indicate that the boundary position varies between the innermost solid wall and peak density position of the CNT as a function of the input energy or work done in the system. Finally, we reveal that the ratio between the potential energy barrier and the kinetic energy is proportional to the shift of the boundary position away from the innermost solid wall.

When the nature of a magnetosonic pulse propagating in a bounded magnetized plasma slab is successively transformed from compression to rarefaction and vice versa upon reflection at a plasma-vacuum interface, both the energy and the longitudinal electromagnetic (EM) momentum of the plasma-pulse system are found to oscillate between two states. Simple analytical models and particle-in-cell simulations show that these oscillations are associated with EM radiation to and from the surrounding magnetized vacuum. For partial reflection supplemental losses in total pulse energy and mechanical momentum are identified and shown to follow, respectively, Fresnel's transmission coefficients for the energy and the magnetic perturbation. This mechanical momentum loss upon partial reflection is traced to the momentarily nonzero volume-integrated Lorentz force, which in turn supports that mechanical and EM momentum transfers are, respectively, associated with the magnetic and electric parts of the momentum flux density.

A computer simulation technique has been applied to the modeling of radiation redistribution functions in low- and moderate-density magnetized hydrogen plasmas. The radiating dipole is described within the Heisenberg picture, and perturbations by the plasma microfield are accounted for through a time-dependent Stark effect term in the Hamiltonian. Numerical applications are presented for the first Lyman and Balmer lines at plasma conditions relevant to tokamak divertors and magnetized white dwarf atmospheres. In both cases, the collisional redistribution of the radiation frequency is shown to be incomplete. Comparisons with a previously developed impact model are performed, and results are discussed.

A hybrid mechanism of ion acceleration is investigated which demonstrates the higher spectral density of protons at high energies. The interaction of few-cycle terrawatt laser pulses with near-critical density gas target is studied with the help of two-dimensional particle-in-cell simulation. The generation of few MeV protons with high spectral concentration near cutoff is attributed to the propagation of solitary waves in the decaying density profile of the gas jet. Plasma dynamics at longer time scale is explained by semianalytical modeling and conditions for solitary wave breaking are presented.

An optimal finite-time process drives a given initial distribution to a given final one in a given time at the lowest cost as quantified by total entropy production. We prove that for a system with discrete states this optimal process involves nonconservative driving, i.e., a genuine driving affinity, in contrast to the case of a system with continuous states. In a multicyclic network, the optimal driving affinity is bounded by the number of states within each cycle. If the driving affects forward and backwards rates nonsymmetrically, the bound additionally depends on a structural parameter characterizing this asymmetry.

We implement parallel versions of the generalized atmospheric Rosenbluth methods and Wang-Landau algorithms for stars and for acyclic uniform branched networks in the square lattice. These are models of monodispersed branched polymers, and we estimate the star vertex exponents ${\sigma }_f$ for $f$ stars, and the entropic exponent ${\gamma }_{\mathcal{G}}$ for networks with comb and brush connectivity in two dimensions. Our results verify the predicted (but not rigorously proven) exact values of the vertex exponents and we test the scaling relation [B. Duplantier, J. Stat. Phys. 54, 581 (1989)] $${\gamma }_{\mathcal{G}}-1=\sum _{f\ge 1}{m}_f{\sigma }_f$$for several acyclic branched networks in two dimensions.

In thermodynamics, entropy production and work quantify irreversibility and the consumption of useful energy, respectively, when a system is driven out of equilibrium. For quantum systems, these quantities can be identified at the stochastic level by unravelling the system's evolution in terms of quantum jump trajectories. We here derive a general formula for computing the joint statistics of work and entropy production in Markovian driven quantum systems, whose instantaneous steady states are of Gibbs form. If the driven system remains close to the instantaneous Gibbs state at all times, then we show that the corresponding two-variable cumulant generating function implies a joint detailed fluctuation theorem so long as detailed balance is satisfied. As a corollary, we derive a modified fluctuation-dissipation relation (FDR) for the entropy production alone, applicable to transitions between arbitrary steady states, and for systems that violate detailed balance. This FDR contains a term arising from genuinely quantum fluctuations, and extends an analogous relation from classical thermodynamics to the quantum regime.

Systems that can be effectively described as a localized spin-$s$ particle subject to time-dependent fields have attracted a great deal of interest due to, among other things, their relevance for quantum technologies. Establishing analytical relationships between the topological features of the applied fields and certain time-averaged quantities of the spin can provide important information for the theoretical understanding of these systems. Here, we address this question in the case of a localized spin-$s$ particle subject to a static magnetic field coplanar to a coexisting elliptically rotating magnetic field. The total field periodically traces out an ellipse which encloses the origin of the coordinate system or not, depending on the values taken on by the static and the rotating components. As a result, two regimes with different topological properties characterized by the winding number of the total field emerge: the winding number is 1 if the origin lies inside the ellipse, and 0 if it lies outside. We show that the time average of the energy associated with the rotating component of the magnetic field is always proportional to the time average of the out-of-plane component of the expectation value of the spin. Moreover, the product of the signs of these two time-averaged quantities is uniquely determined by the topology of the total field and, consequently, provides a measurable indicator of this topology. We also propose an implementation of these theoretical results in a trapped-ion quantum system. Remarkably, our findings are valid in the totality of the parameter space and regardless of the initial state of the spin. In particular, when the system is prepared in a Floquet state, we demonstrate that the quasienergies, as a function of the driving amplitude at constant eccentricity, have stationary points at the topological transition boundary. The ability of the topological indicator proposed here to accurately locate the abrupt topological transition can have practical applications for the determination of unknown parameters appearing in the Hamiltonian. In addition, our predictions about the quasienergies can assist in the interpretation of conductance measurements in transport experiments with spin carriers in mesoscopic rings.

Machine learning (ML) has been well applied to studying equilibrium phase transition models by accurately predicating critical thresholds and some critical exponents. Difficulty will be raised, however, for integrating ML into nonequilibrium phase transitions. The extra dimension in a given nonequilibrium system, namely time, can greatly slow down the procedure toward the steady state. In this paper we find that by using some simple techniques of ML, non-steady-state configurations of directed percolation (DP) suffice to capture its essential critical behaviors in both ($1+1$) and ($2+1$) dimensions. With the supervised learning method, the framework of our binary classification neural networks can identify the phase transition threshold, as well as the spatial and temporal correlation exponents. The characteristic time $t_c$, specifying the transition from active phases to absorbing ones, is also a major product of the learning. Moreover, we employ the convolutional autoencoder, an unsupervised learning technique, to extract dimensionality reduction representations and cluster configurations of ($1+1$) bond DP. It is quite appealing that such a method can yield a reasonable estimation of the critical point.

We investigate the dynamics of a charged particle confined to move on a toroidal helix while being driven by an external time-dependent electric field. The underlying phase space is analyzed for linearly and circularly polarized fields. For small driving amplitudes and a linearly polarized field, we find a split up of the chaotic part of the phase space, which prevents the particle from inverting its direction of motion. This allows for a nonzero average velocity of chaotic trajectories without breaking the well-known symmetries commonly responsible for directed transport. Within our chosen normalized units, the resulting average transport velocity is constant and does not change significantly with the driving amplitude. A very similar effect is found in case of the circularly polarized field and low driving amplitudes. Furthermore, when driving with a circularly polarized field, we unravel a second mechanism of the split up of the chaotic phase space region for very large driving amplitudes. There exists a wide range of parameter values for which trajectories may travel between the two chaotic regions by crossing a permeable cantorus. The limitations of these phenomena, as well as their implication on manipulating directed transport in helical geometries are discussed.

Affinity maturation (AM) is the process through which the immune system is able to develop potent antibodies against new pathogens it encounters, and is at the base of the efficacy of vaccines. At its core AM is analogous to a Darwinian evolutionary process, where B cells mutate and are selected on the base of their affinity for an antigen (Ag), and Ag availability tunes the selective pressure. In cases when this selective pressure is high, the number of B cells might quickly decrease and the population might risk extinction in what is known as a population bottleneck. Here we study the probability for a B-cell lineage to survive this bottleneck scenario as a function of the progenitor affinity for the Ag. Using recursive relations and probability generating functions we derive expressions for the average extinction time and progeny size for lineages that go extinct. We then extend our results to the full population, both in the absence and presence of competition for T-cell help, and quantify the population survival probability as a function of Ag concentration and initial population size. Our study suggests the population bottleneck phenomenology might represent a limit case in the space of biologically plausible maturation scenarios, whose characterization could help guide the process of vaccine development.