Raman side scatter, whereby scattered light is resonant while propagating perpendicularly to a density gradient in a plasma, was identified experimentally in planar-target experiments at the National Ignition Facility at intensities orders of magnitudes below the threshold for absolute instability. We have derived a new theoretical description of convective Raman side scatter below the absolute threshold, validated by numerical simulations. We show that inertial confinement fusion experiments at full ignition scale, i.e., with mm-scale spot sizes and density scale lengths, are prone to increased coupling losses from Raman side scatter as the instability can extend from the absolute regime near the quarter-critical density to the convective regime at lower electron densities.

We investigate the propagation of electronic waves described by the Dirac equation subject to a Lévy-type disorder distribution. Our numerical calculations, based on the transfer matrix method, in a system with a distribution of potential barriers show that it presents a phase transition from anomalous to standard to anomalous localization as the incidence energy increases. In contrast, electronic waves described by the Schrödinger equation do not present such transitions. Moreover, we obtain the phase diagram delimiting anomalous and standard localization regimes, in the form of an incidence angle versus incidence energy diagram, and argue that transitions can also be characterized by the behavior of the dispersion of the transmission. We attribute this transition to an abrupt reduction in the transmittance of the system when the incidence angle is higher than a critical value which induces a decrease in the transmission fluctuations.

Many animal groups are heterogeneous and may even consist of individuals of different species, called mixed-species flocks. Mathematical and computational models of collective animal movement behavior, however, typically assume that groups and populations consist of identical individuals. In this paper, using the mathematical framework of the coagulation-fragmentation process, we develop and analyze a model of merge and split group dynamics, also called fission-fusion dynamics, for heterogeneous populations that contain two types (or species) of individuals. We assume that more heterogeneous groups experience higher split rates than homogeneous groups, forming two daughter groups whose compositions are drawn uniformly from all possible partitions. We analytically derive a master equation for group size and compositions and find mean-field steady-state solutions. We predict that there is a critical group size below which groups are more likely to be homogeneous and contain the abundant type or species. Despite the propensity of heterogeneous groups to split at higher rates, we find that groups are more likely to be heterogeneous but only above the critical group size. Monte Carlo simulation of the model show excellent agreement with these analytical model results. Thus, our model makes a testable prediction that composition of flocks are group-size-dependent and do not merely reflect the population level heterogeneity. We discuss the implications of our results to empirical studies on flocking systems.

Weak wave turbulence has been observed on a thin elastic plates since the work by Düring et al. [Phys. Rev. Lett. 97, 025503 (2006)]. Here we report theoretical, experimental, and numerical studies of wave turbulence in a forced thin elastic plate submitted to increasing tension. When increasing the tension (or decreasing the bending stiffness of the plate) the plate evolves progressively from a plate into an elastic membrane as in drums. We first consider a thin plate and increase the tension in experiments and numerical simulations. We observe that the system remains in a state of weak turbulence of weakly dispersive waves. This observation is in contrast with what has been observed in water waves when decreasing the water depth, which also changes the waves from dispersive to weakly dispersive. The weak turbulence observed in the deep water case evolves into a solitonic regime. Here no such transition is observed for the stretched plate. We then apply the weak turbulence theory to the membrane case and show with numerical simulations that indeed the weak turbulence framework remains valid for the membrane and no formation of singular structures (shocks) should be expected in contrast with acoustic wave turbulence.

Quantum heat transfer is analyzed in nonequilibrium two-qubits systems by applying the nonequilibrium polaron-transformed Redfield equation combined with full counting statistics. Steady-state heat currents with weak and strong qubit-bath couplings are clearly unified. Within the two-terminal setup, the negative differential thermal conductance is unraveled with strong qubit-bath coupling and finite qubit splitting energy. The partially strong spin-boson interaction is sufficient to show the negative differential thermal conductance. Based on the three-terminal setup, in which two qubits are asymmetrically coupled to three thermal baths, a giant heat amplification factor is observed with strong qubit-bath coupling. Moreover, the strong interaction of either the left or right spin-boson coupling is able to exhibit the apparent heat amplification effect.

Information storage, reflecting the capability of a dynamical system to keep predictable information during its evolution over time, is a key element of intrinsic distributed computation, useful for the description of the dynamical complexity of several physical and biological processes. Here we introduce a parametric approach which allows one to compute information storage across multiple timescales in stochastic processes displaying both short-term dynamics and long-range correlations (LRC). Our analysis is performed in the popular framework of multiscale entropy, whereby a time series is first “coarse grained” at the chosen timescale through low-pass filtering and downsampling, and then its complexity is evaluated in terms of conditional entropy. Within this framework, our approach makes use of linear fractionally integrated autoregressive (ARFI) models to derive analytical expressions for the information storage computed at multiple timescales. Specifically, we exploit state space models to provide the representation of lowpass filtered and downsampled ARFI processes, from which information storage is computed at any given timescale relating the process variance to the prediction error variance. This enhances the practical usability of multiscale information storage, as it enables a computationally reliable quantification of a complexity measure which incorporates the effects of LRC together with that of short-term dynamics. The proposed measure is first assessed in simulated ARFI processes reproducing different types of autoregressive dynamics and different degrees of LRC, studying both the theoretical values and the finite sample performance. We find that LRC alter substantially the complexity of ARFI processes even at short timescales, and that reliable estimation of complexity can be achieved at longer timescales only when LRC are properly modeled. Then, we assess multiscale information storage in physiological time series measured in humans during resting state and postural stress, revealing unprecedented responses to stress of the complexity of heart period and systolic arterial pressure variability, which are related to the different role played by LRC in the two conditions.

This paper addresses a keystone problem for the progress of phononics: the proposal of efficient thermal diodes. Aiming the disclosure of an easy itinerary for the building of a heat rectifier, I investigate unsophisticated systems linked to simple thermal baths, precisely, asymmetric quantum Ising models, i.e., simple quadratic models, involving only one spin component. I analytically show the occurrence of thermal rectification for the case of a chain with interactions long enough to connect the first to the last site. Moreover, I describe cases of a perfect rectification, i.e., finite heat flow in one direction and zero current in the opposite direction. I argue to indicate that the ingredients for the rectification are just given by the quantum nature of the baths and dynamics, and by the structural asymmetry of the system, here in the intersite interactions. I believe that the description of a perfect thermal rectification in a simple many-body quantum model, that is, the presentation of a simple itinerary for the building of a diode shall stimulate theoretical and experimental research on the theme.

We present a simple thermodynamically consistent method for solving time-dependent Fokker-Planck equations (FPE) for overdamped stochastic processes, also known as Smoluchowski equations. It yields both transition and steady-state behavior and allows for computations of moment-generating and large-deviation functions of observables defined along stochastic trajectories, such as the fluctuating particle current, heat, and work. The key strategy is to approximate the FPE by a master equation with transition rates in configuration space that obey a local detailed balance condition for arbitrary discretization. Its time-dependent solution is obtained by a direct computation of the time-ordered exponential, representing the propagator of the FPE, by summing over all possible paths in the discretized space. The method thus not only preserves positivity and normalization of the solutions but also yields a physically reasonable total entropy production, regardless of the discretization. To demonstrate the validity of the method and to exemplify its potential for applications, we compare it against Brownian-dynamics simulations of a heat engine based on an active Brownian particle trapped in a time-dependent quartic potential.

Identifying the qualitative changes in time-series data provides insights into the dynamics associated with such data. Such qualitative changes can be detected through topological approaches, which first embed the data into a high-dimensional space using a time-delay parameter and subsequently extract topological features describing the shape of the data from the embedded points. However, the essential topological features that are extracted using a single time delay are considered to be insufficient for evaluating the aforementioned qualitative changes, even when a well-selected time delay is used. We therefore propose a delay-variant embedding method that constructs the extended topological features by considering the time delay as a variable parameter instead of considering it as a single fixed value. This delay-variant embedding method reveals multiple-timescale patterns in a time series by allowing the observation of the variations in topological features, with the time delay serving as an additional dimension in the topological feature space. We theoretically prove that the constructed topological features are robust when the time series is perturbed by noise. Furthermore, we combine these features with the kernel technique in machine learning algorithms to classify the general time-series data. We demonstrate the effectiveness of our method for classifying the synthetic noisy biological and real time-series data. Our method outperforms a method that is based on a single time delay and, surprisingly, achieves the highest classification accuracy on an average among the standard time-series analysis techniques.

The magnetization dynamics of nano-oscillators may be excited by both magnetic fields and spin-polarized currents. While the dynamics of single oscillators has been well characterized, the synchronization of several ones is not fully understood yet. An analytical and numerical study of the nonlinear dynamics of two magnetostatically coupled spin valves driven by spin-transfer torques is presented under the macrospin approximation. The oscillators interact via magnetostatic fields and exhibit a robust synchronized magnetization motion. We describe the magnetization dynamics of the system using the Landau-Lifshitz-Gilbert-Slonczewski equation. Using a modal decomposition technique, we describe the dynamics, synchronization, and competition of oscillatory modes as a function of the current density, and the geometrical parameters of the setup. Simulations of the Landau-Lifshitz-Gilbert-Slonczewski equation show good agreement with an approximate analytic solution.

We numerically investigate the dynamics of strongly disordered 1D lattices under single-particle displacements, using both the Hertzian model, describing a granular chain, and the $\alpha +\beta$ Fermi-Pasta-Ulam-Tsingou model (FPUT). The most profound difference between the two systems is the discontinuous nonlinearity of the granular chain appearing whenever neighboring particles are detached. We therefore sought to unravel the role of these discontinuities in the destruction of Anderson localization and their influence on the system's chaotic dynamics. Our results show that the dynamics of both models can be characterized by: (i) localization with no chaos; (ii) localization and chaos; (iii) spreading of energy, chaos, and equipartition. The discontinuous nonlinearity of the Hertzian model is found to trigger energy spreading at lower energies. More importantly, a transition from Anderson localization to energy equipartition is found for the Hertzian chain and is associated with the “propagation” of the discontinuous nonlinearity in the chain. On the contrary, the FPUT chain exhibits an alternate behavior between localized and delocalized chaotic behavior which is strongly dependent on the initial energy excitation.

Homoclinic and heteroclinic orbits provide a skeleton of the full dynamics of a chaotic dynamical system and are the foundation of semiclassical sums for quantum wave packets, coherent states, and transport quantities. Here, the homoclinic orbits are organized according to the complexity of their phase-space excursions, and exact relations are derived expressing the relative classical actions of complicated orbits as linear combinations of those with simpler excursions plus phase-space cell areas bounded by stable and unstable manifolds. The total number of homoclinic orbits increases exponentially with excursion complexity, and the corresponding cell areas decrease exponentially in size as well. With the specification of a desired precision, the exponentially proliferating set of homoclinic orbit actions is expressible by a slower-than-exponentially increasing set of cell areas, which may present a means for developing greatly simplified semiclassical formulas.

Several experiments have shown that the stimulation of a single neuron in the cortex can influence the local network activity and even the behavior of an animal. From the theoretical point of view, it is not clear how stimulating a single cell in a cortical network can evoke a statistically significant change in the activity of a large population. Our previous study considered a random network of integrate-and-fire neurons and proposed a way of detecting the stimulation of a single neuron in the activity of a local network: a threshold detector biased toward a specific subset of neurons. Here, we revisit this model and extend it by introducing a second network acting as a readout. In the simplest scenario, the readout consists of a collection of integrate-and-fire neurons with no recurrent connections. In this case, the ability to detect the stimulus does not improve. However, a readout network with both feed-forward and local recurrent inhibition permits detection with a very small bias, if compared to the readout scheme introduced previously. The crucial role of inhibition is to reduce global input cross correlations, the main factor limiting detectability. Finally, we show that this result is robust if recurrent excitatory connections are included or if a different kind of readout bias (in the synaptic amplitudes instead of connection probability) is used.

We calculate the probabilities that a trajectory of a fractional Brownian motion with arbitrary fractal dimension $d_f$ visits the same spot $n\ge 3$ times, at given moments $t_1,...,t_n$, and obtain a determinant expression for these probabilities in terms of a displacement-displacement covariance matrix. Except for the standard Brownian trajectories with $d_f=2$, the resulting many-body contact probabilities cannot be factorized into a product of single-loop contributions. Within a Gaussian network model of a self-interacting polymer chain, which we suggested recently [K. Polovnikov et al., Soft Matter 14, 6561 (2018)], the probabilities we calculate here can be interpreted as probabilities of multibody contacts in a fractal polymer conformation with the same fractal dimension $d_f$. This Gaussian approach, which implies a mapping from fractional Brownian motion trajectories to polymer conformations, can be used as a semiquantitative model of polymer chains in topologically stabilized conformations, e.g., in melts of unconcatenated rings or in the chromatin fiber, which is the material medium containing genetic information. The model presented here can be used, therefore, as a benchmark for interpretation of the data of many-body contacts in genomes, which we expect to be available soon in, e.g., Hi-C experiments.

We study theoretically the rotation induced on an uncharged metal nanocylinder immersed in an electrolyte by AC electric fields. We consider the rotation of the cylinder when subjected to a rotating electric field (electrorotation) and the orientation of the cylinder in an AC field with constant direction (electro-orientation). The cylinder rotation is due to two mechanisms: the electric field interaction with the induced dipole on the particle and the hydrodynamic stress on the particle originated by the induced-charge electro-osmotic (ICEO) flow around the particle. The cylinder rotation induced by the ICEO mechanism can be calculated by using the Lorentz reciprocal theorem, while the rotation due to the induced dipole is calculated from the cylinder polarizability. We employ 3D numerical computations using finite elements for the general case as well as analytical methods for slender cylinders. Both calculations use the thin-double-layer approximation. We compare the results for slender cylinders of both methods showing good agreement. The electro-orientation (EOr) due to dipole torque aligns the axis of slender cylinders with the applied field, but aligns the axis of short cylinders perpendicularly to the field. The EOr due to ICEO torque always aligns the axis of cylinders with the field. The rotation induced by ICEO torque tends to disappear for frequencies of the applied field much greater than the characteristic frequency for charging the double-layer capacitance of the metal-electrolyte interface.

We study the low-frequency dielectric response of highly charged spheres arranged in a cubic lattice and immersed in an electrolyte solution. We focus on the influence of the out-of-phase current in the regime where the ionic charge is neutral. We consider the case where the charged spheres have no surface conductance and no frequency-dependent surface capacitance. Hence, the frequency dispersion of the dielectric constant is dominated by the effect of neutral currents outside the electric double layer. In the thin double-layer limit, we use Fixman's boundary condition at the outer surface of the double layer to capture the interaction between the electric field and the flow of the ions. For periodic conditions, we combine the methods developed by Lord Rayleigh for understanding the electric conduction across rectangularly arranged obstacles and by Korringa, Kohn, and Rostoker for the electronic band structure computation. When the charged spheres occupy a very small volume fraction, smaller than 1%, our solution becomes consistent with the Maxwell Garnett mixing formula together with the single-particle polarization response, as expected, because interparticle interactions become less prominent in the dilute limit. By contrast, the interparticle interaction greatly alters the dielectric response even when charged spheres occupy only 2% of the volume. We found that the characteristic frequency shifts to a higher value compared to that derived from the single-particle polarization response. At the same time, the low-frequency dielectric enhancement, a signature of charged spheres immersed in an electrolyte, becomes less prominent for the periodic array of charged spheres. Our results imply that the signature of the dielectric response of a system consisting of densely packed charged spheres immersed in an electrolyte can differ drastically from a dilute suspension.

The theoretical description of the reorientational dynamics in the case of a hybrid-oriented two-dimensional (2D) liquid crystal (LC) cell, under the influence of a temperature gradient $\nabla T$, caused by a heat flux $\mathbf{\text{q}}$ directed at an angle $\alpha$ across the lower bounding surface with the orientational defect, has been presented. Our calculations, based on the appropriate nonlinear extension of the classical Ericksen-Leslie theory, show that due to interaction between $\nabla T$ and the gradient of the director field $\nabla \widehat{\mathbf{\text{n}}}$ in the LC sample a thermally excited vortical fluid flow is maintained in the bulk of the nematic volume. In order to elucidate the role of both the orientational defect and the heat flux $\mathbf{\text{q}}$ in formation of the vortical flow in the microsized LC volume, we have analyzed the response of the LC material confined in the microsized 2D volume on the effect of the laser beam focused on the bounding surface.

We propose a surface energy for conically degenerate anchorings of uniaxial liquid crystal mesogens by modifying tensorial Nobili-Durand surface energy that is usually employed for fixed anchoring orientations with preferred polar angles. By minimizing Landau–de Gennes free energy and the proposed surface energy, we obtain the equilibrium director configuration around a spherical colloid in the uniform nematic liquid crystal. Our calculations show that the proposed surface energy can cause boojum or/and Saturn-ring defect textures depending on the equilibrium conic angle. We also study the interactions between two spherical colloids with the equilibrium conic angle ${45}^{\circ }$, where the surface energy provides both boojum and Saturn-ring defects on the surface of particles. We compare the calculated anisotropic colloidal interactions with experimental observations [B. Senyuk et al., Nat. Commun. 7, 10659 (2016)]. In agreement with experiment, our results show two stable angular assemblies in the close particle-particle separations. Also, the long-range elastic interactions are almost consistent with the hexadecapolar elastic distortion.

In this paper we use simulations of the magnetohydrodynamic equations coupled with heat to show the generation of magnetic field by the dynamical interaction of a pair of vortices in a fluid electrically conducting within a cylindrical domain nonhomogeneously heated from below, setting in a rotation frame. For large enough rotation rates we show that the formation of a pair of vortices inside the primary whirl gives rise to a magnetic field. The magnetic field has a strong horizontal component, and the magnetic lines go from one vortex to the other.