In recent years, kirigami techniques have inspired the design of graphene-based nanodevices with exceptional stretchability and ductility. Based on an I-shaped cutting pattern, here we propose a graphene kirigami that exhibits remarkable stretchability and ductility in two independent planar directions along with negative Poisson's ratios. The deformation mechanism underlying the high stretchability of the structure is the flipping and rotation of its cutting ligaments during elongation. Molecular dynamics simulations show that the yield and fracture strains of graphene kirigami can be enhanced by factors of 6 and 10 in the two planar directions. In addition, the mechanical properties of the graphene kirigami can be tuned by altering the cutting geometric parameters as well as incorporating distinct cutting patterns in series. We develop a numerical algorithm to predict the stress-strain response of the series-connected graphene kirigami, and verify its accuracy using appropriate simulations. On this basis, the stress-strain response of the series-connected graphene kirigami can be tuned by altering its geometric parameters and the number of building blocks. This graphene kirigami could be applied to the design and development of next-generation flexible electronics such as stretchable electrodes and strain sensors.

Diffusing diffusivity models, polymers in the grand canonical ensemble and polydisperse, and continuous-time random walks all exhibit stages of non-Gaussian diffusion. Is non-Gaussian targeting more efficient than Gaussian? We address this question, central to, e.g., diffusion-limited reactions and some biological processes, through a general approach that makes use of Jensen's inequality and that encompasses all these systems. In terms of customary mean first-passage time, we show that Gaussian searches are more effective than non-Gaussian ones. A companion paper argues that non-Gaussianity becomes instead highly more efficient in applications where only a small fraction of tracers is required to reach the target.

Coexistence of segregation and interconversion or isomerization reaction among molecular species leads to fascinating structure formation in the biological and chemical worlds. Using Monte Carlo simulations of the prototype Ising model, we explore the chemical kinetics of such a system consisting of a binary mixture of isomers. Our results reveal that even though the two concerned processes are individually Arrhenius in nature, the Arrhenius behavior of the isomerization reaction gets significantly disrupted due to an interplay of the nonconserved dynamics of the reaction and the conserved diffusive dynamics of segregation. The approach used here can be potentially adapted to understand reaction kinetics of more complex reactions.

Despite the periodic kicks, a linear kicked rotor (LKR) is an integrable and exactly solvable model in which the kinetic energy term is linear in momentum. It was recently shown that spatially interacting LKRs are also integrable, and results in dynamical localization in the corresponding quantum regime. Similar localized phases exist in other nonintegrable models such as the coupled relativistic kicked rotors. This work, using a two-body LKR, demonstrates two main results; first, it is shown that chaos can be induced in the integrable linear kicked rotor through interactions between the momenta of rotors. An analytical estimate of its Lyapunov exponent is obtained. Second, the quantum dynamics of this chaotic model, upon variation of kicking and interaction strengths, is shown to exhibit a variety of phases: classically induced localization, dynamical localization, subdiffusive, and diffusive phases. We point out the signatures of these phases from the perspective of entanglement production in this system. By defining an effective Hilbert space dimension, the entanglement growth rate can be understood using appropriate random matrix averages.

Braided channel networks exhibit a complex interplay between spatial and temporal dynamics. Their behavior is characterized by both simple and multiscaling patterns, and the mechanisms underlying the stochastic processes associated with this dynamics remain incompletely understood. Leveraging Taylor's pioneering work [Nature (London) 189, 732 (1961)], which unveiled scaling relations in a plethora of natural phenomena through what is now known as the Taylor power law (TPL), we propose a physical interpretation of braided channel systems. This interpretation utilizes a specific class of transformation functions applied to the mean of fluvial geomorphic variables measured along cross sections, namely, the number of wet channels, the average width of wet channels, and the entropic braiding index. By analyzing remotely sensed data of the Brahmaputra-Jamuna River in Bangladesh we obtain valuable insight into the spatiotemporal scaling of these geomorphological variables and gather a deeper understanding of the complexity of braided channel systems. Finally, through a direct analysis employing the TPL in conjunction with a fixed-mass multifractal algorithm, we prove that braided channel networks exhibit a multiscaling behavior.

We study the ordering dynamics of nonlinear voter models with multiple states, also providing a discussion of the two-state model. The rate with which an individual adopts an opinion scales as the $q\mathrm{th}$ power of the number of the individual's neighbors in that state. For $q>1$ the dynamics favor the opinion held by the most agents. The ordering to consensus is driven by deterministic drift, and noise plays only a minor role. For $q<1$ the dynamics favors minority opinions, and for multistate models the ordering proceeds through a noise-driven succession of metastable states. Unlike linear multistate systems, the nonlinear model cannot be reduced to an effective two-state model. We find that the average density of active interfaces in the model with multiple opinion states does not show a single exponential decay in time for $q<1$, again at variance with the linear model. This highlights the special character of the conventional (linear) voter model, in which deterministic drift is absent. As part of our analysis, we develop a pair approximation for the multistate model on graphs, valid for any positive real value of $q$, improving on previous approximations for nonlinear two-state voter models.

A harmonically trapped active Brownian particle exhibits two types of positional distributions—one has a single peak and the other has a single well—that signify steady-state dynamics with low and high activity, respectively. Adding inertia to the translational motion preserves this strict classification of either single-peak or single-well densities but shifts the dividing boundary between the states in the parameter space. We characterize this shift for the dynamics in one spatial dimension using the static Fokker-Planck equation for the full joint distribution of the state space. We derive local results analytically with a perturbation method for a small rotational velocity and then extend them globally with a numerical approach.

The National Ignition Facility has recently achieved successful burning plasma and ignition using the inertial confinement fusion (ICF) approach. However, there are still many fundamental physics phenomena that are not well understood, including the kinetic processes in the hohlraum. Shan et al. [Phys. Rev. Lett. 120, 195001 (2018)] utilized the energy spectra of neutrons to investigate the kinetic colliding plasma in a hohlraum of indirect drive ICF. However, due to the typical large spatial-temporal scales, this experiment could not be well simulated by using available codes at that time. Utilizing our advanced high-order implicit PIC code, LAPINS, we were able to successfully reproduce the experiment on a large scale of both spatial and temporal dimensions, in which the original computational scale was increased by approximately seven to eight orders of magnitude. Not only is the validity of the explanation of the experiment confirmed by our simulations, i.e., the abnormally large width of neutron spectra comes from beam-target nuclear fusions, but also a different physical insight into the source of energetic deuterium ions is provided. The acceleration of deuterium ions can be categorized into two components: one is propelled by a sheath electric field created by the charge separation at the onset, while the other is a result of the reflection of the potential of the shock wave. The robustness of the acceleration mechanism is analyzed with varying initial conditions, e.g., temperatures, drifting velocity, and ion components. This paper might serve as a reference for benchmark simulations of upcoming simulation codes and may be relevant for future research on mixtures and entropy increments at plasma interfaces.

Two quantum systems, each described as a random-matrix ensemble, are coupled to each other via a number of transition states. Each system is strongly coupled to a large number of channels. The average transmission probability is the product of three factors describing, respectively, formation of the first system from the entrance channel, decay of the second system through the exit channel, and transport through the transition states. Each of the transition states contributes a Breit-Wigner resonance. In general, the resonances overlap.

It is expected that the energy-diffusion propagator in a one-dimensional nonlinear lattice with three conserved quantities: energy, momentum, and stretch, consists of a central heat mode and two sound modes. The heat mode follows a Lévy distribution. Consequently, the heat diffusion is super, i.e., the second moment of the diffusion propagator diverges as $t^{\beta }$ with $\beta >1$; and the heat conduction is anomalous, i.e., the heat conductivity is size dependent and diverges with size $N$ by $N^{\alpha }$, with $\alpha >0$. In this paper, we study a one-dimensional lattice with two-dimensional transverse motions, in which the total angular momentum also conserves. More importantly, the diffusion of this conserved quantity is ballistic. Surprisingly, the above pictures and the values of the mentioned power exponents keep unchanged. The universality of the scalings is then further extended. On the other hand, the detailed strengths of heat transports are largely enhanced. Such a counterintuitive finding can be explained by the change of the phonon mean-free path of the lattices.

A quantum wave function with localization on classical periodic orbits in a mesoscopic elliptic billiard has been achieved by appropriately superposing nearly degenerate eigenstates expressed as products of Mathieu functions. We analyze and discuss the rotational and librational regimes of motion in the elliptic billiard. Simplified line equations corresponding to the classical trajectories can be extracted from the quantum state as an integral equation involving angular Mathieu functions. The phase factors appearing in the integrals are connected to the classical initial positions and velocity components. We analyze the probability current density, phase maps, and vortex distributions of the periodic orbit quantum states for both rotational and librational motions; furthermore, they may represent traveling and standing trajectories inside the elliptic billiard.

Double-layer targets (DLTs), made of a low-density foam on top of a solid substrate, can efficiently convert the energy of a high-intensity laser to provide sources of photons and protons. We investigate a 30-fs pulse with a peak intensity of $I\sim 8.7\times {10}^{20}\mathrm{W}/{\mathrm{cm}}^2$ and a peak power of $\sim 120$ TW interacting with a DLT using three-dimensional (3D) particle-in-cell simulations. We focus on providing quantitative results in full 3D geometry on the foam thickness dependence; on the competition between two photon-generating processes in DLTs, i.e., nonlinear inverse Compton scattering (NICS) and bremsstrahlung (BS); and on the acceleration of protons via enhanced target-normal sheath acceleration. We discuss conversion efficiency, average energy, and angular distributions of such multiradiation sources. We find that NICS can prevail over BS if the DLT's substrate is thin enough $(\sim \mu \mathrm{m})$ and that the optimal foam thickness that maximizes the conversion efficiency in NICS and BS photons and the proton cutoff energy, among those considered, is the same $\left(15\mu \mathrm{m}\right)$. These results show that DLTs constitute an excellent tool for developing relatively compact and optimized laser-driven multicomponent radiation sources.

We consider scattering processes where a quantum system comprises an inner subsystem and a boundary and is subject to Haar-averaged random unitaries acting on the boundary-environment Hilbert space only. We show that, regardless of the initial state, a single scattering event will disentangle the unconditional state (i.e., the scattered state when no information about the applied unitary is available) across the inner subsystem-boundary partition. Also, we apply Lévy's lemma to constrain the trace norm fluctuations around the unconditional state. Finally, we derive analytical formulas for the mean scattered purity for initial globally pure states and provide one with numerical evidence of the reduction of fluctuations around such mean values with increasing environmental dimension.

Molecular motors fulfill critical functions within all living beings. Understanding their underlying working principles is therefore of great interest. Here we develop a simple model inspired by the two-component biomolecular motor ${\text{F}}_{\text{o}}\text{−}{\text{F}}_1$ ATP synthase. We analyze its energetics and characterize information flows between the machine's components. At maximum output power we find that information transduction plays a minor role for free-energy transduction. However, when the two components are coupled to different environments (e.g., when in contact with heat baths at different temperatures), we show that information flow becomes a resource worth exploiting to maximize free-energy transduction. Our findings suggest that real-world powerful and efficient information engines could be found in machines whose components are subjected to fluctuations of different strength, since in this situation the benefit gained from using information for work extraction can outweigh the costs of information generation.

We address the question of the time needed by $N$ particles, initially located on the first sites of a finite one-dimensional lattice of size $L$, to exit that lattice when they move according to a TASEP transport model. Using analytical calculations and numerical simulations, we show that when $N\ll L$, the mean exit time of the particles is asymptotically given by $T_N\left(L\right)\sim L+{\beta }_N\sqrt{L}$ for large lattices. Building upon exact results obtained for two particles, we devise an approximate continuous space and time description of the random motion of the particles that provides an analytical recursive relation for the coefficients ${\beta }_N$. The results are shown to be in very good agreement with numerical results. This approach sheds some light on the exit dynamics of $N$ particles in the regime where $N$ is finite while the lattice size $L\to \infty$. This complements previous asymptotic results obtained by Johansson [Commun. Math. Phys. 209, 437 (2000)] in the limit where both $N$ and $L$ tend to infinity while keeping the particle density $N/L$ finite.

Many technical and biomedical applications of magnetic nanoparticles rely on their response to time-varying magnetic fields. While well-established models exist for either immobile or thermally blocked nanoparticles, the intermediate regime where Brownian as well as Néel relaxation occur at the same time is less well explored. Here, we use an efficient model that allows us to study the nonlinear dynamics of individual magnetic nanoparticles in response to different time-varying magnetic fields over a broad range of field parameters, taking into account both relaxation mechanisms. We provide quasiexact solutions for the longitudinal dynamics as well as approximate formulas from dynamic mean-field theory. Our results are relevant, e.g., for magnetorelaxometry, magnetic fluid hyperthermia, and magnetic particle imaging. For these example applications, we show that the ratio of characteristic Brownian to Néel relaxation time can have a profound impact on characteristic response quantities, especially at large field strengths.

While bouncing walking silicone oil droplets (walkers) do show many quantumlike phenomena, the original, most intriguing, double-slit experiment with walkers has been shown to be an overinterpretation of data. Several experiments and numerical simulations have proven that for at least some parameter region there is no randomness. Still, true randomness was claimed to be observed in an experiment on chaotically bouncing walkers. Also, most of the available phase space has not been investigated. The main goal of this paper is therefore to look for true interference and chaos in the entire phase space. Recently, we made an extensive investigation of drops interacting with slits, but still in a limited range. However, the outcome was always deterministic and only incidentally mimicked the statistics of the corresponding quantum case. We also showed that the extra interference, already seen by others, in the double-slit case was caused by reflection of waves from the outlet of the unused slit after passage and thus was not a true double-slit effect. A new theoretical treatment of bouncing drop dynamics has since given analytic relations for the associated wave field, leading to a proposal for criteria for the possible occurrence of true interference in the double-slit experiment. Satisfying these criteria, requires working close to the onset of the Faraday instability, with two spatial conditions favoring slow walkers, and a temporal condition favoring fast walkers. The regions of high velocity, where the walkers bounce periodically, and of very low velocity, with chaotically bouncing walkers, have not been comprehensively investigated so far. We have therefore looked at these regions, probing the limits for interaction with slits. Furthermore, noting that a short transit time is essential to fulfill the criteria, experiments were done using double-slit barriers only 0.5 and 2 mm broad. Nowhere was true interference or a chaotic response found. As the theory has implications for many of the observed quantumlike phenomena involving walkers as, e.g., tunneling and interaction between drops, we have measured the spatial and temporal decay of the wave field. A comparison with the theory shows very good agreement but leads to a reformulation of the above-mentioned criteria.

The elastic-plastic Richtmyer-Meshkov instability of multiple interfaces is investigated by numerical simulation using a multimaterial solid mechanics algorithm based on an Eulerian framework. This Richtmyer-Meshkov instability problem is realized by a copper layer that is flanked by vacuum and a copper block of different material strength. The research efforts are directed to reveal the influence of the layer thickness and material strength on the deformation of the perturbed solid-vacuum interface impacted by an initial shock. By varying the initial thickness ($x_{\mathrm{I}}$) of the copper layer and the yield stress (${\sigma }_{\mathrm{Y}2}$) of the copper block, two deformation modes, which have been identified as the broken mode and the stable mode, are closely scrutinized. For a fixed $x_{\mathrm{I}}$ and a decreasing ${\sigma }_{\mathrm{Y}2}$, the reflected rarefaction waves (RRWs), developing after the initial shock impacts the perturbed interface 1 (I1) between vacuum and the copper layer, become stronger after traveling across the interface 2 (I2). Subsequently, the velocity of I2 becomes larger, causing the width of I1 to grow larger. This width growth of I1 leads to a final separation of the spike from I1 and, consequently, the deformation mode changes from the stable mode to the broken mode. For a fixed ${\sigma }_{\mathrm{Y}2}$ and a decreasing $x_{\mathrm{I}}$, the RRWs impact I2 at an earlier moment with a greater strength and thus the deformation mode changes from the stable mode to the broken mode. Meanwhile, the comparison of the spike width of cases whose deformation mode is the broken mode shows that there exists a maximum value of rescaled spike width, at which the deformation mode changes from the stable mode to the broken mode.

The double-cone ignition (DCI) scheme has been proposed as one of the alternative approaches to inertial confinement fusion, based on direct-drive and fast-ignition, in order to reduce the requirement for the driver energy. To evaluate the conical implosion energetics from the laser beams to the plasma flows, a series of experiments have been systematically conducted. The results indicate that $89\%–96\%$ of the laser energy was absorbed by the target, with moderate stimulated Raman scatterings. Here $2\%–6\%$ of the laser energy was coupled into the plasma jets ejected from the cone tips, which was mainly restricted by the mass reductions during the implosions inside the cones. The supersonic dense jets with a Mach number of 4 were obtained, which is favorable for forming a high-density, nondegenerated plasma core after the head-on collisions. These findings show encouraging results in terms of energy transport of the conical implosions in the DCI scheme.

We study the dynamical large deviations of the classical stochastic symmetric simple exclusion process (SSEP) by means of numerical matrix product states. We show that for half filling, long-time trajectories with a large enough imbalance between the number hops in even and odd bonds of the lattice belong to distinct symmetry-protected topological (SPT) phases. Using tensor network techniques, we obtain the large deviation (LD) phase diagram in terms of counting fields conjugate to the dynamical activity and the total hop imbalance. We show the existence of high activity trivial and nontrivial SPT phases (classified according to string order parameters) separated by either a critical phase or a critical point. Using the leading eigenstate of the tilted generator, obtained from infinite-system density-matrix renormalization group simulations, we construct a near-optimal dynamics for sampling the LDs, and show that the SPT phases manifest at the level of rare stochastic trajectories. We also show how to extend these results to other filling fractions, and discuss generalizations to asymmetric SEPs.