The excitation of the terahertz (THz) high-intensity surface modes when the two-frequency $p$-polarized laser radiation interacts with a plasma slab is studied. It was found that the significant amplification of the laser field in the plasma slab occurs when p-polarized laser radiation is incident at the angle of total reflection. It is shown that, under the action of laser radiation ponderomotive forces, the resonant excitation of the THz mode of the plasma slab occurs if the frequency difference of the laser fields coincides with the eigenfrequency of the surface mode. It is established that the giant increase in the energy flux density of the THz mode occurs when $p$-polarized laser radiation is incident at the angle of total reflection on the near-critical plasma slab with rare electron collisions if the conditions of resonant excitation are satisfied. It is shown that in this case the energy flux density of THz mode can significantly exceed the laser intensity.

We numerically study, from first principles, the temperature $T_{\mathrm{QSS}}$ and duration $t_{\mathrm{QSS}}$ of the longstanding initial quasi-stationary state of the isolated $d$-dimensional classical inertial $\alpha \text{−}XY$ ferromagnet with two-body interactions decaying as $1/r_{ij}^{\alpha }$ $(\alpha \ge 0)$. It is shown that this temperature $T_{\mathrm{QSS}}$ (defined proportional to the kinetic energy per particle) depends, for the long-range regime $0\le \alpha /d\le 1$, on $(\alpha ,d,U,N)$ with numerically negligible changes for dimensions $d=1,2,3$, with $U$ the energy per particle and $N$ the number of particles. We verify the finite-size scaling $T_{\mathrm{QSS}}-T_{\infty }\propto 1/N^{\beta }$ where $T_{\infty }\equiv 2U-1$ for $U\lesssim U_c$, and $\beta$ appears to depend sensibly only on $\alpha /d$. Our numerical results indicate that neither the scaling with $N$ of $T_{\mathrm{QSS}}$ nor that of $t_{\mathrm{QSS}}$ depend on $U$.

We investigate density fluctuations in a liquid close to a solvophobic substrate at which a surface critical drying transition occurs. Using classical density functional theory, we determine three measures of the spatial extent and strength of the fluctuations, i.e., the local compressibility ${\chi }_{\mu }\left(z\right)$, the local thermal susceptibility ${\chi }_T\left(z\right)$, and the reduced density ${\chi }_{*}\left(z\right)$; $z$ is the distance from the substrate. While the first measure is frequently used, the second and third were introduced very recently by Eckert et al., [T. Eckert et al., Phys. Rev. Lett. 125, 268004 (2020)]. For state points in the critical drying regime, all three measures, each scaled by its bulk value, exhibit very similar forms and the ratio of ${\chi }_T\left(z\right)$ to ${\chi }_{\mu }\left(z\right)$, for fixed $z$ in the vapor-liquid transition region, is constant. Using a scaling treatment of surface thermodynamics, we show that such behavior is to be expected on general grounds.

Recent progress in studying the physics of amorphous solids has revealed that mechanical strains can be strongly screened by the formation of plastic events that are typically quadrupolar in nature. The theory stipulates that gradients in the density of the quadrupoles act as emergent dipole sources, leading to strong screening and to qualitative changes in the mechanical response, as seen, for example, in the displacement field. In this Letter we first offer direct measurements of the dipole field, independently of any theoretical assumptions, and second we demonstrate detailed agreement with the recently proposed theory. These two goals are achieved by using data from both simulations and experiments. Finally, we show how measurements of the dipole fields pinpoint the theory parameters that determine the profile of the displacement field.

The work and heat distribution densities of a classical Brownian particle immersed in a heat bath and interacting with an external field have been extensively analyzed from different approaches. In this article, a previous method based on basic principles of stochastic dynamics is extended to derive these functions. The starting point is the inertial Langevin equation and the external field is an off-center harmonic potential driven by an external protocol. Unlike previous works where the driving is arbitrary, the so-called optimal protocol that minimizes the mechanical work is used instead. The corresponding work and heat distributions are derived through a procedure based on getting a generic Fokker-Planck equation indistinctly of the variable under consideration. The work distribution is calculated for different initial conditions and values of the friction coefficient of the thermal fluid ranging from the periodic or very low-underdamped mode up to the overdamped regime. It is a Gaussian as that of previous experiments of a particle trapped in an optical tweezers moved at constant velocity. Some aspects about the heat distribution is analyzed in terms of the statistical features of the non-Gaussian noise accompanying its dynamics to give an account of experimental results. It is concluded that the easiness in calculating the work distribution cannot be applied to heat. It requires numerical calculations.

Strong nonlinearity of a self-resonant radio-frequency (rf) superconducting-quantum-interference-device (SQUID) meta-atom is explored via intermodulation (IM) measurements. Previous work in zero dc magnetic flux showed a sharp onset of IM response as the frequency sweeps through the resonance. A second onset at higher frequency was also observed, creating a prominent gap in the IM response. By extending those measurements to nonzero dc flux, different dynamics are revealed, including dc flux tunability of the aforementioned gaps and enhanced IM response near geometric resonance of the rf SQUID. These features observed experimentally are understood and analyzed theoretically through a combination of a steady-state analytical modeling and a full numerical treatment of the rf SQUID dynamics. The latter in addition predicts the presence of chaos in narrow parameter regimes. The understanding of intermodulation in rf SQUID metamaterials is important for producing low-noise amplification of microwave signals and tunable filters.

In this paper, we study the dynamics of an interesting class of vector solitons in the long-wave–short-wave resonance interaction (LSRI) system. The model that we consider here describes the nonlinear interaction of long wave and two short waves and it generically appears in several physical settings. To derive this class of nondegenerate vector soliton solutions we adopt the Hirota bilinear method with the more general form of admissible seed solutions with nonidentical distinct propagation constants. We express the resultant fundamental as well as multisoliton solutions in a compact way using Gram-determinants. The general fundamental vector soliton solution possesses several interesting properties. For instance, the double-hump or a single-hump profile structure including a special flattop profile form results in when the soliton propagates in all the components with identical velocities. Interestingly, in the case of nonidentical velocities, the soliton number is increased to two in the long-wave component, while a single-humped soliton propagates in the two short-wave components. We establish through a detailed analysis that the nondegenerate multisolitons in contrast to the already known vector solitons (with identical wave numbers) can undergo three types of elastic collision scenarios: (i) shape-preserving, (ii) shape-altering, and (iii) a shape-changing collision, depending on the choice of the soliton parameters. Here, by shape-altering we mean that the structure of the nondegenerate soliton gets modified slightly during the collision process, whereas if the changes occur appreciably then we call such a collision as shape-changing collision. We distinguish each of the collision scenarios, by deriving a zero phase shift criterion with the help of phase constants. Very importantly, the shape-changing behavior of the nondegenerate vector solitons is observed in the long-wave mode also, along with corresponding changes in the short-wave modes, and this nonlinear phenomenon has not been observed in the already known vector solitons. In addition, we point out the coexistence of nondegenerate and degenerate solitons simultaneously along with the associated physical consequences. We also indicate the physical realizations of these general vector solitons in nonlinear optics, hydrodynamics, and Bose-Einstein condensates. Our results are generic and they will be useful in these physical systems and other closely related systems including plasma physics when the long-wave–short-wave resonance interaction is taken into account.

The present work reports an investigation into the collisional dynamics of particles in the vicinity of the outlet of a two-dimensional silo using molecular dynamics simulations. Most studies on this granular system focus in the bulk of the medium. In this region, contacts are permanent or long-lived, so continuous approximations are able to yield results for velocity distributions or mass flow. Close to the exit, however, the density of the medium decreases and contacts are instantaneous. Thus, the collisional nature of the dynamics becomes significant, warranting a dedicated investigation as carried out in this work. More interesting, the vicinity of the outlet is the region where the arches that block the flow for small apertures are formed. It is found that the transition from the clogging regime (at small apertures) to the continuous flow regime is smooth in collisional variables. Furthermore, the dynamics of particles as reflected by the distributions of the velocities is as well unaffected. This result implies that there is no critical outlet size that separates both regimes, as had been proposed in the literature. Instead, the results achieved support the alternative picture in which a clog is possible for any outlet size.

In strongly magnetized neutral plasmas, electron motion is reduced perpendicular to the magnetic field direction. This changes dynamical plasma properties such as temperature equilibration, spatial density evolution, electron pressure, and thermal and electrical conductivity. In this paper we report measurements of free plasma expansion in the presence of a strong magnetic field. We image laser-induced fluorescence from an ultracold neutral ${\mathrm{Ca}}^{+}$ plasma to map the plasma size as a function of time for a range of magnetic field strengths. The asymptotic expansion velocity perpendicular to the magnetic field direction falls rapidly with increasing magnetic field strength. We observe that the initially Gaussian spatial distribution remains Gaussian throughout the expansion in both the parallel and perpendicular directions. We compare these observations with a diffusion model and with a self-similar expansion model and show that neither of these models reproduces the observed behavior over the entire range of magnetic fields used in this study. Modeling the expansion of a magnetized ultracold plasma poses a nontrivial theoretical challenge.

We propose an unsupervised learning methodology with descriptors based on topological data analysis (TDA) concepts to describe the local structural properties of materials at the atomic scale. Based only on atomic positions and without a priori knowledge, our method allows for an autonomous identification of clusters of atomic structures through a Gaussian mixture model. We apply successfully this approach to the analysis of elemental Zr in the crystalline and liquid states as well as homogeneous nucleation events under deep undercooling conditions. This opens the way to deeper and autonomous study of complex phenomena in materials at the atomic scale.

A nominally two-dimensional spin model wrapped onto a cylinder can profitably be viewed, especially for long cylinders, as a one-dimensional chain. Each site of such a chain is a ring of spins with a complex state space. Traditional correlation functions are inadequate for the study of correlations in such a system and need to be replaced with something like mutual information. Being induced purely by frustration, the disorder of a cylindrical zero-temperature triangular Ising antiferromagnet (TIAFM) and attendant correlations have a chance of evading the consequences of the Perron-Frobenius theorem which describes and constrains correlations in thermally disordered one-dimensional systems. Correlations in such TIAFM systems and the aforementioned evasion are studied here through a fermionic representation. For cylindrical TIAFM models with open boundary conditions, we explain and derive the following characteristics of end-to-end mutual information: period-three oscillation of the decay length, halving of the decay length compared to what Perron-Frobenius predicts on the basis of transfer matrix eigenvalues, and subexponential decay—inverse square in the length—for certain systems.

The canonical ensemble plays a crucial role in statistical mechanics in and out of equilibrium. For example, the standard derivation of the fluctuation theorem relies on the assumption that the initial state of the heat bath is the canonical ensemble. On the other hand, the recent progress in the foundation of statistical mechanics has revealed that a thermal equilibrium state is not necessarily described by the canonical ensemble but can be a quantum pure state or even a single energy eigenstate, as formulated by the eigenstate thermalization hypothesis (ETH). Then a question raised is how these two pictures, the canonical ensemble and a single energy eigenstate as a thermal equilibrium state, are compatible in the fluctuation theorem. In this paper, we theoretically and numerically show that the fluctuation theorem holds in both of the long- and short-time regimes, even when the initial state of the bath is a single energy eigenstate of a many-body system. Our proof of the fluctuation theorem in the long-time regime is based on the ETH, while it was previously shown in the short-time regime on the basis of the Lieb-Robinson bound and the ETH [Phys. Rev. Lett. 119, 100601 (2017)]. The proofs for these time regimes are theoretically independent and complementary, implying the fluctuation theorem in the entire time domain. We also perform a systematic numerical simulation of hard-core bosons by exact diagonalization and verify the fluctuation theorem in both of the time regimes by focusing on the finite-size scaling. Our results contribute to the understanding of the mechanism that the fluctuation theorem emerges from unitary dynamics of quantum many-body systems and can be tested by experiments with, e.g., ultracold atoms.

The appealing theoretical measure of irreversibility in a stochastic process, as the ratio of the probabilities of a trajectory and its time reversal, cannot be accessed directly in experiments since the probability of a single trajectory is zero. We regularize this definition by considering, instead, the limiting ratio of probabilities for trajectories to remain in the tubular neighborhood of a smooth path and its time reversal. The resulting pathwise medium entropy production agrees with the formal expression from stochastic thermodynamics and can be obtained from measurable tube probabilities. Estimating the latter from numerically sampled trajectories for Langevin dynamics yields excellent agreement with theory. By combining our measurement of pathwise entropy production with a Markov chain Monte Carlo algorithm, we infer the entropy-production distribution for a transition path ensemble directly from short recorded trajectories. Our work enables the measurement of irreversibility along individual paths and path ensembles in a model-free manner.

Monte Carlo simulations on the structural change of hard sphere-platelet mixtures were performed to investigate the effect of particle size. We quantitatively analyzed local equilibrium structures of sphere-platelet mixtures with varying size ratios under various sphere and platelet density conditions. Based on the simulation results, we investigated the structural transitions such as isotropic to anisotropic, clustering, and so on. When a small amount of small-sized sphere is added to a large-sized platelet system, the mixture structure transitions from isotropic to nematic ones as the platelet number density increases. On the other hand, the platelet forms clusters with the addition of a large number of spheres. In a small platelet–large sphere system, the spheres form aggregates by increasing platelet density instead. The platelet and spherical particles exhibit different structural transitions depending on the size and density. In the limit of small and large size ratios, the structures of the platelet-sphere mixture obtained from the Monte Carlo simulation are close to those shown by previous theoretical and experimental studies, respectively. Because the primary actor shifts from sphere to platelet as the size ratio changes, the transition boundary shifts continuously. When the size ratio is close to unity, the most complicated behavior is observed, with both the platelet and sphere simultaneously acting the leading part.

An analysis of the IR absorbance for the segmented functional groups of liquid crystal dimers mesogen and linker enabled the orientation order to be determined and information about the dipole interactions in the nematic and twist-bend nematic phases to be obtained. The long axis orientational order increases as the temperature decreases in the nematic phase, although much more slowly than for the classical nematics, and then reverses this trend in the twist-bend nematic phase due to the tilt of the molecules. In the nematic phase, the short axis of the molecule performs an isotropic uniform rotation and has a uniaxial alignment. In the twist-bend nematic phase, however, biaxial ordering occurs and grows significantly in accordance with the helical deformation of the director. Changes in the mean absorbance in the twist-bend nematic phase were observed: a decrease for the longitudinal dipole at the nematic-twist-bend nematic phase transition, thus emphasizing the antiparallel axial interaction of the dipoles, while the absorbance of the transverse dipoles remains unchanged up to 340 K, and then the latter become parallelly correlated.

We report the discovery and elucidation of giant spatiotemporal orientational fluctuations in nematic liquid crystal drops with radial orientation of the nematic anisotropy axis producing a central “hedgehog” defect. We study the spatial and temporal properties of the fluctuations experimentally using polarized optical microscopy, and theoretically, by calculating the eigenspectrum of the Frank elastic free energy of a nematic drop of radius $R_2$, containing a spherical central core of radius $R_1$ and constrained by perpendicular boundary conditions on all surfaces. We find that the hedgehog defect with radial orientation has a complex excitation spectrum with a single critical mode whose energy vanishes at a critical value ${\mu }_c$ of the ratio $\mu =R_2/R_1$. When $\mu <{\mu }_c$, the mode has positive energy, indicating that the radial hedgehog state is stable; when $\mu >{\mu }_c$, it has negative energy indicating that the radial state is unstable to the formation of a lower-energy state. This mode gives rise to the large-amplitude director fluctuations we observe near the core, for $\mu$ near ${\mu }_c$. A collapse of the experimental data corroborates model predictions for $\mu <{\mu }_c$ and provides an estimate of the defect core size.

We present calculations of eigenmode energies and wave functions of both azimuthal and polar distortions of the nematic director relative to a radial hedgehog trapped in a spherical drop with a smaller concentric spherical droplet at its core. All surfaces interior to the drop have perpendicular (homeotropic) boundary conditions. We also calculate director correlation functions and their relaxation times. Of particular interest is a critical mode whose energy, with fixed Frank constants, vanishes as the ratio $\mu =R_2/R_1$ increases toward a critical value ${\mu }_c$, where $R_2$ is the radius of the drop and $R_1$ that of the inner droplet, and then becomes negative for $\mu >{\mu }_c$. Our calculations form a basis for interpreting experimental measurements of director fluctuations relative to a radial hedgehog state in a spherical drop. We compare results with those obtained by previous investigations, which use a calculational approach different from ours, and with our experimental observations.