A scheme for enhanced quantum electrodynamics (QED) production of electron-positron-pair plasmas is proposed that uses two ultraintense lasers irradiating a thin solid foil from opposite sides. In the scheme, under a proper matching condition, in addition to the skin-depth emission of $\gamma$-ray photons and Breit-Wheeler creation of pairs on each side of the foil, a large number of high-energy electrons and photons from one side can propagate through it and interact with the laser on the other side, leading to much enhanced $\gamma$-ray emission and pair production. More importantly, the created pairs can be collected later and confined to the center by opposite laser radiation pressures when the foil becomes transparent, resulting in the formation of unprecedentedly overdense and high-energy pair plasmas. Two-dimensional QED particle-in-cell simulations show that electron-positron-pair plasmas with overcritical density ${10}^{22}{\mathrm{cm}}^{-3}$ and a high energy of 100s of MeV are obtained with 10 PW lasers at intensities ${10}^{23}\mathrm{W}/{\mathrm{cm}}^2$, which are of key significance for laboratory astrophysics studies.

We present a self-consistent semianalytical model of the relativistic plasma mirror, based on the exact computation of the laser-driven electron surface oscillations within the cold-fluid approximation. Valid for arbitrary solid densities, laser incidence angle, and a large set of laser intensities (${10}^{18}–{10}^{21}\text{W}/{\text{cm}}^2$), the model unravels different regimes of harmonic generation. In particular, it is found that efficient conversion of $p$-polarized laser pulses into high-order harmonics well above the plasma frequency requires either high laser intensities, low plasma densities, or incidence angles larger than a threshold value. This critical angle corresponds to a transition between a regime where the electron surface dynamics is mostly governed by the laser $J\times B$ force and a “cyclotron Brunel” regime, where electrons perform many cyclotron gyrations when moving into the vacuum. Under conditions relevant to current laser experiments, the latter regime gives rise to nonmonotonic variations of the harmonic yield with the laser field. Our predictions are supported by an extensive parametric study performed with highly resolved one-dimensional particle-in-cell simulations.

We report here energetic $(>100\mathrm{keV})$ negative hydrogen ions $\left({\mathrm{H}}^{-}\right)$ generated in the interaction of moderately intense $\left({10}^{18}\mathrm{W}\mathrm{c}{\mathrm{m}}^{-2}\right)$ ultrashort laser pulses (45 fs) with transparent hydrogen containing solid targets. An unambiguous and consistent detection of negative hydrogen ions, with a flux of $8\times {10}^{11}{\mathrm{H}}^{-}$ ions/sr, has been observed in every single laser shot, using a Thomson parabola ion spectrograph. Simple estimates based on charge transfer cross sections match well with experimental observations. Our method offers the implementation of an intense, ultrashort laser based negative-ion source at a higher repetition rate, which can be important for various applications.

In this work, a physical system described by the Hamiltonian ${\mathbf{H}}_{\omega }={\mathbf{H}}_0+{\mathbf{V}}_{\omega }\left(t\right)$ consisting of a solvable model ${\mathbf{H}}_0$ and external random and time-dependent potential ${\mathbf{V}}_{\omega }\left(t\right)$ is investigated. Under the conditions in which, for each realization, the potential changes smoothly so that the evolution of the system follows the Schrödinger dynamics, and that the average external potential with respect to all realizations is constant in time, an adjusted equilibrium state can be defined as a reference state and the mean dynamics can be derived from taking the average of the equation with respect to the configuration parameter $\omega$. It provides extra contributions from the deviations of the Hamiltonian and evolves the state along the time by the Heisenberg and Liouville–von Neumann equations. Consequently, the Kubo formula and the fluctuation-dissipation relation (FDR) are modified in the sense that the contribution from the information of randomness and memory effects from the time dependence is also present. The modified Kubo formula now has a contribution from two terms. The first term is an antisymmetric cross correlation between two observables measured by a probe as expected, and the latter term is an accumulation of the propagation of the effects from the randomness. When the considered system is in the adjusted equilibrium state at the time the measurement probe interacts, the latter contribution vanishes, and the standard FDR is recovered.

Many phenomena in nature are described by excitable systems driven by colored noise. The temporal correlations in the fluctuations hinder an analytical treatment. We here present a general method of reduction to a white-noise system, capturing the color of the noise by effective and time-dependent boundary conditions. We apply the formalism to a model of the excitability of neuronal membranes, the leaky integrate-and-fire neuron model, revealing an analytical expression for the linear response of the system valid up to moderate frequencies. The closed form analytical expression enables the characterization of the response properties of such excitable units and the assessment of oscillations emerging in networks thereof.

A long-standing goal of nonequilibrium statistical mechanics has been to extend the conceptual power of the Boltzmann distribution to driven systems. We report some new progress towards this goal. Instead of writing the nonequilibrium steady-state distribution in terms of perturbations around thermal equilibrium, we start from the linearized driven dynamics of observables about their stable fixed point, and expand in the strength of the nonlinearities encountered during typical fluctuations away from the fixed point. The first terms in this expansion retain the simplicity of known expansions about equilibrium, but can correctly describe the statistics of a certain class of systems even under strong driving. We illustrate this approach by comparison with a numerical simulation of a sheared Brownian colloid, where we find that the first two terms in our expansion are sufficient to account for the shear thinning behavior at high shear rates.

In this paper we have presented the heat exchange between the two fermionic thermal reservoirs which are connected by a fermionic system. We have calculated the heat flux using solution of the $c$-number Langevin equation for the system. Assuming small temperature difference between the baths we have defined the thermal conductivity for the process. It first increases as a nonlinear function of average temperature of the baths to a critical value then decreases to a very low value such that the heat flux almost becomes zero. There is a critical temperature for the fermionic case at which the thermal conductivity is maximum for the given coupling strength and the width of the frequency distribution of bath modes. The critical temperature grows if these quantities become larger. It is a sharp contrast to the Bosonic case where the thermal conductivity monotonically increases to the limiting value. The change of the conductivity with increase in width of the frequency distribution of the bath modes is significant at the low temperature regime for the fermionic case. It is highly contrasting to the Bosonic case where the signature of the enhancement is very prominent at high temperature limit. We have also observed that thermal conductivity monotonically increases as a function of damping strength to the limiting value at the asymptotic limit. There is a crossover between the high and the low temperature results in the variation of the thermal conductivity as a function of the damping strength for the fermionic case. Thus it is apparent here that even at relatively high temperature, the fermionic bath may be an effective one for the strong coupling between system and reservoir. Another interesting observation is that at the low temperature limit, the temperature dependence of the heat flux is the same as the Stefan-Boltzmann law. This is similar to the bosonic case.

We present a linear stability analysis of a planar metal electrode during steady electrodeposition. We extend the previous work of Sundstrom and Bark by accounting for the extended space-charge density, which develops at the cathode once the applied voltage exceeds a few thermal voltages. In accordance with Chazalviel's conjecture, the extended space-charge region is found to greatly affect the morphological stability of the electrode. To supplement the numerical solution of the stability problem, we have derived analytical expressions valid in the limit of low and high voltage, respectively.

We describe experiments and computer simulations of molecular deposition on a substrate in which the molecules (substituted adenine derivatives) self-assemble into ordered structures. The resulting structures depend strongly on the deposition rate (flux). In particular, there are two competing surface morphologies $(\alpha$ and $\beta )$, which differ by their topology (interdigitated vs lamellar structure). Experimentally, the $\alpha$ phase dominates at both low and high flux, with the $\beta$ phase being most important in the intermediate regime. A similar nonmonotonic behavior is observed on varying the substrate temperature. To understand these effects from a theoretical perspective, a lattice model is devised which reproduces qualitatively the topological features of both phases. Via extensive Monte Carlo studies we can, on the one hand, reproduce the experimental results and, on the other hand, obtain a microscopic understanding of the mechanisms behind this anomalous behavior. The results are discussed in terms of an interplay between kinetic trapping and temporal exploration of configuration space.

Cyclic dominance models, such as the classic rock-paper-scissors (RPS) game, have found real-world applications in biology, ecology, and sociology. A key quantity of interest in such models is the coexistence time, i.e., the time until at least one population type goes extinct. Much recent research has considered conditions that lengthen coexistence times in an RPS model. A general finding is that coexistence is promoted by localized spatial interactions (low mobility), while extinction is fostered by global interactions (high mobility). That is, there exists a mobility threshold which separates a regime of long coexistence from a regime of rapid collapse of coexistence. The key finding of our paper is that if zealots (i.e., nodes able to defeat others while themselves being immune to defeat) of even a single type exist, then system coexistence time can be significantly prolonged, even in the presence of global interactions. This work thus highlights a crucial determinant of system survival time in cyclic dominance models.

We report on the multistability of chaotic networks of silicon neurons and demonstrate how spatiotemporal sequences of voltage oscillations are selected with timed current stimuli. A three neuron central pattern generator was built by interconnecting Hodgkin-Huxley neurons with mutually inhibitory links mimicking gap junctions. By systematically varying the timing of current stimuli applied to individual neurons, we generate the phase lag maps of neuronal oscillators and study their dependence on the network connectivity. We identify up to six attractors consisting of triphasic sequences of unevenly spaced pulses propagating clockwise and anticlockwise. While confirming theoretical predictions, our experiments reveal more complex oscillatory patterns shaped by the ratio of the pulse width to the oscillation period. Our work contributes to validating the command neuron hypothesis.

The field-theoretic renormalization group and the operator product expansion are applied to the model of passive vector (magnetic) field advected by a random turbulent velocity field. The latter is governed by the Navier-Stokes equation for compressible fluid, subject to external random force with the covariance $\propto \delta (t-t^{\prime })k^{4-d-y}$, where $d$ is the dimension of space and $y$ is an arbitrary exponent. From physics viewpoints, the model describes magnetohydrodynamic turbulence in the so-called kinematic approximation, where the effects of the magnetic field on the dynamics of the fluid are neglected. The original stochastic problem is reformulated as a multiplicatively renormalizable field-theoretic model; the corresponding renormalization group equations possess an infrared attractive fixed point. It is shown that various correlation functions of the magnetic field and its powers demonstrate anomalous scaling behavior in the inertial-convective range already for small values of $y$. The corresponding anomalous exponents, identified with scaling (critical) dimensions of certain composite fields (“operators” in the quantum-field terminology), can be systematically calculated as series in $y$. The practical calculation is performed in the leading one-loop approximation, including exponents in anisotropic contributions. It should be emphasized that, in contrast to Gaussian ensembles with finite correlation time, the model and the perturbation theory presented here are manifestly Galilean covariant.

It is shown that the filamentation instability of relativistically intense laser pulses in plasmas can be mitigated in the case where the laser beam has an elliptically distributed beam profile. A high-power elliptical Gaussian laser beam would break up into a regular filamentation pattern—in contrast to the randomly distributed filaments of a circularly distributed laser beam—and much more laser power would be concentrated in the central region. A highly elliptically distributed laser beam experiences anisotropic self-focusing and diffraction processes in the plasma channel ensuring that the unstable diffractive rings of the circular case cannot be produced. The azimuthal modulational instability is thereby suppressed. These findings are verified by three-dimensional particle-in-cell simulations.

Measurements of the structure of the electrostatic fields produced by the expansion of a laser-produced plasma into a background magnetized plasma are presented. The three-dimensional measurements of the electrostatic field are made using an emissive probe that measures the time-varying plasma potential on two orthogonal planes, one across and one containing the background magnetic field. The inductive electric field is also calculated from probe measurements of the time-varying magnetic fields. Deviations from local charge neutrality at the level of ${10}^{-4}$ generate a radial electrostatic field with peak strength an order of magnitude larger than the corresponding inductive field. The electrostatic energy density near full expansion is over an order of magnitude larger than that of the induced azimuthal electric field. These measurements show that electrostatic fields must be included in theoretical and computational models of collisionless coupling in magnetized point explosions of laser-produced plasmas and their relation to similar phenomena such as magnetospheric chemical releases.

We reexamine the well-studied one-dimensional spin-$1/2 XY$ model to reveal its nontrivial energy spectrum, in particular the energy gap between the ground state and the first excited state. In the case of the isotropic $XY$ model, the $XX$ model, the gap behaves very irregularly as a function of the system size at a second order transition point. This is in stark contrast to the usual power-law decay of the gap and is reminiscent of the similar behavior at the first order phase transition in the infinite-range quantum $XY$ model. The gap also shows nontrivial oscillatory behavior for the phase transitions in the anisotropic model in the incommensurate phase. We observe a close relation between this anomalous behavior of the gap and the correlation functions. These results, those for the isotropic case in particular, are important from the viewpoint of quantum annealing where the efficiency of computation is strongly affected by the size dependence of the energy gap.

Purely dissipative evolution equations are often cast as gradient flow structures, $\dot{\mathbf{z}}=K\left(\mathbf{z}\right)DS\left(\mathbf{z}\right)$, where the variable $\mathbf{z}$ of interest evolves towards the maximum of a functional $S$ according to a metric defined by an operator $K$. While the functional often follows immediately from physical considerations (e.g., the thermodynamic entropy), the operator $K$ and the associated geometry does not necessarily do so (e.g., Wasserstein geometry for diffusion). In this paper, we present a variational statement in the sense of maximum entropy production that directly delivers a relationship between the operator $K$ and the constraints of the system. In particular, the Wasserstein metric naturally arises here from the conservation of mass or energy, and depends on the Onsager resistivity tensor, which, itself, may be understood as another metric, as in the steepest entropy ascent formalism. This variational principle is exemplified here for the simultaneous evolution of conserved and nonconserved quantities in open systems. It thus extends the classical Onsager flux-force relationships and the associated variational statement to variables that do not have a flux associated to them. We further show that the metric structure $K$ is intimately linked to the celebrated Freidlin-Wentzell theory of stochastically perturbed gradient flows, and that the proposed variational principle encloses an infinite-dimensional fluctuation-dissipation statement.

We study droplet formation in granular suspensions by systematically varying the volume fractions ($\phi$) and particle diameters ($d$). For suspensions with water as the suspending liquid, we find three different regimes. For dilute suspensions ($\phi \le 45\%$), drop formation follows the predictions for inertial breakup and exhibits identical dynamics to that of pure water. The breakup is strongly asymmetrical in this case. Only for more concentrated suspensions ($\phi >45\%$) does the presence of particles change the dynamics and two other regimes, a symmetrical inertial regime and a Bagnoldian regime, are uncovered. We construct and discuss a phase diagram that allows us to understand and predict the breakup behavior in granular suspensions.

To elucidate mechanisms of mesoscopic turbulence exhibited by active particles, we experimentally study turbulent states of nonliving self-propelled particles. We realize an experimental system with dense suspensions of asymmetrical colloidal particles (Janus particles) self-propelling on a two-dimensional surface under an ac electric field. Velocity fields of the Janus particles in the crowded situation can be regarded as a sort of turbulence because it contains many vortices and their velocities change abruptly. Correlation functions of their velocity field reveal the coexistence of polar alignment and antiparallel alignment interactions, which is considered to trigger mesoscopic turbulence. Probability distributions of local order parameters for polar and nematic orders indicate the formation of local clusters with particles moving in the same direction. A broad peak in the energy spectrum of the velocity field appears at the spatial scales where the polar alignment and the cluster formation are observed. Energy is injected at the particle scale and conserved quantities such as energy could be cascading toward the larger clusters.

It has been hypothesized that in the era just before the last universal common ancestor emerged, life on earth was fundamentally collective. Ancient life forms shared their genetic material freely through massive horizontal gene transfer (HGT). At a certain point, however, life made a transition to the modern era of individuality and vertical descent. Here we present a minimal model for stochastic processes potentially contributing to this hypothesized “Darwinian transition.” The model suggests that HGT-dominated dynamics may have been intermittently interrupted by selection-driven processes during which genotypes became fitter and decreased their inclination toward HGT. Stochastic switching in the population dynamics with three-point (hypernetwork) interactions may have destabilized the HGT-dominated collective state and essentially contributed to the emergence of vertical descent and the first well-defined species in early evolution. A systematic nonlinear analysis of the stochastic model dynamics covering key features of evolutionary processes (such as selection, mutation, drift and HGT) supports this view. Our findings thus suggest a viable direction out of early collective evolution, potentially enabling the start of individuality and vertical Darwinian evolution.