We investigated the properties of cholesteric liquid crystals (CLCs) being in an external static magnetic field directed along the helix axis. We considered a dichroic CLC, that is, CLC with parameters $\mathrm{Re}\mathrm{\Delta }=\frac{\mathrm{Re}{\varepsilon }_1-\mathrm{Re}{\varepsilon }_2}{2}=0$ and $\mathrm{Im}\mathrm{\Delta }=\frac{\mathrm{Im}{\varepsilon }_1-\mathrm{Im}{\varepsilon }_2}{2}\ne 0$, where ${\varepsilon }_{1,2}$ are the principal values of the local dielectric permittivity tensor. We have shown that in the case of the wavelength dependence of the magneto-optic activity parameter, new features appear in the optics of dichroic CLCs, in particular, in this case new Dirac points appear. Dirac points are points where there is an intersection of any two wave vector curves (they degenerate) and a linear law of the wave vector dependence on the frequency near these points. And moreover, at some Dirac points photonic band gaps (PBGs) appear; at others, lines of magnetically induced transparency (MIT), that is, a full transmission band appears, in an absorbing medium. In this case a polarization-sensitive transmission band appears too. At certain values of the helix pitch of the CLC and of the magnitude of the external magnetic field, three PBGs of different nature appear: a transmittance band, two narrow lines of MIT, and at others a broadband MIT, etc. This system is nonreciprocal, and the nonreciprocity changes over a wide range. It is observed both for reflection and transmittance and for absorption. The soft-matter nature of CLCs and their response to external influences lead to easily tunable multifunctional devices that can find a variety of applications. They can be applied as tunable narrow-band or broadband filters and mirrors, a highly tunable broad/narrow-band coherent perfect absorber, transmitter, ideal optical diode, and in other devices.

The boundary layer near a cooled inclined plate, which is immersed in a stably stratified fluid rotating about an axis parallel to the direction of gravity, is a model for katabatic flows at high latitudes. In this paper the base flow of such an inclined buoyancy layer is solved analytically for arbitrary Prandtl numbers. By applying linear stability analyses, five unstable modes are identified for both the fixed temperature and the isoflux boundary conditions, i.e., the stationary longitudinal roll (LR) mode, the oblique roll with low streamwise wave-number (OR-1) and high streamwise wave-number (OR-2) modes, and the Tolmien-Schlichting (TS) wave with low streamwise wave-number (TS-1) and high streamwise wave-number (TS-2) modes. It is indicated that the Coriolis effect induced by the rotation leads the critical modes to be three dimensional, and a larger tilt angle of the plate and stronger Coriolis effect cause both TS wave modes to be more unstable for both thermal boundary conditions. When the Coriolis effect is considered, the OR-1 and OR-2 modes are the most unstable mode at low and high tilt angles, respectively, but the TS-1 wave mode may be the most unstable one when the plate is nearly vertical. In addition, the spanwise phase velocities of the TS wave modes change directions as the tilt angle passes some threshold values for both thermal boundary conditions except for the TS-1 wave mode with a fixed temperature boundary condition, which propagates in the same spanwise direction for all explored tilt angles.

We study the self-organization of spherical particles in an oscillating flow through experiments inside an oscillating box. The interactions between the particles and the time-averaged (steady streaming) flow lead to the formation of either one-particle-thick chains or multiple-particle-wide bands, depending on the oscillatory conditions. Both the chains and the bands are oriented perpendicular to the direction of oscillation with a regular spacing between them. For all our experiments, this spacing is only a function of the relative particle-fluid excursion length normalized by the particle diameter, $A_r/D$, implying that it is an intrinsic quantity that is established only by the hydrodynamics. In contrast, the width of the bands depends on both $A_r/D$ and the confinement, characterized by the particle coverage fraction $\varphi$. Using the relation for the chain spacing, we accurately predict the transition from one-particle-thick chains to wider bands as a function of $\varphi$ and $A_r/D$. Our experimental results are complemented with numerical simulations in which the flow around the particles is fully resolved. These simulations show that the regular chain spacing arises from the balance between long-range attractive and short-range repulsive hydrodynamic interactions, caused by the vortices in the steady streaming flow. We further show that these vortices induce an additional attractive interaction at very short range when $A_r/D\gtrsim 0.7$, which stabilizes the multiple-particle-wide bands. Finally, we give a comprehensive overview of the parameter space where we illustrate the different regions using our experimental data.

Human beings are easily impacted by social influences, due to their social nature. As an essential manifestation of social influences, conformity is associated with the frequency witnessed in others' behavior, involving normative conformity and informational conformity according to the reaction of individuals. The former comes from the fear of a normative environment, while the latter means most behaviors are followed due to information asymmetry. Normative conformity significantly enhances network reciprocity, producing optimal cooperation at a moderate proportion, which induces within-cluster behavioral homogeneity and between-cluster behavioral diversity. On the contrary, informational conformity has an inhibitory effect on the evolution of cooperation for a low proportion of the conformity population, which contributes to the formation of defectors' clusters. The symmetry and duality of the two types of conformity on cooperation evolution provide an interesting and unexplored approach for future research, revealing the mechanism of conformity in evolutionary games.

We propose a hybrid model governed by the Blume-Emery-Griffiths (BEG) Hamiltonian with a mean-field-like interaction, where the spins are randomly quenched such that some of them are “pure” Ising and the others admit the BEG set of states. It is found, by varying the concentration of the Ising spins, that the model displays different phase portraits in concentration-temperature parameter space, within the canonical and the microcanonical ensembles. Phenomenological indications that these portraits are rich and rather unusual are provided.

Fractional diffusion and Fokker-Planck equations are widely used tools to describe anomalous diffusion in a large variety of complex systems. The equivalent formulations in terms of Caputo or Riemann-Liouville fractional derivatives can be derived as continuum limits of continuous-time random walks and are associated with the Mittag-Leffler relaxation of Fourier modes, interpolating between a short-time stretched exponential and a long-time inverse power-law scaling. More recently, a number of other integrodifferential operators have been proposed, including the Caputo-Fabrizio and Atangana-Baleanu forms. Moreover, the conformable derivative has been introduced. We study here the dynamics of the associated generalized Fokker-Planck equations from the perspective of the moments, the time-averaged mean-squared displacements, and the autocovariance functions. We also study generalized Langevin equations based on these generalized operators. The differences between the Fokker-Planck and Langevin equations with different integrodifferential operators are discussed and compared with the dynamic behavior of established models of scaled Brownian motion and fractional Brownian motion. We demonstrate that the integrodifferential operators with exponential and Mittag-Leffler kernels are not suitable to be introduced to Fokker-Planck and Langevin equations for the physically relevant diffusion scenarios discussed in our paper. The conformable and Caputo Langevin equations are unveiled to share similar properties with scaled and fractional Brownian motion, respectively.

Proton-decoupled deuterium NMR spectra were obtained for an asymmetric liquid crystal dimer 1-(4-cyanobiphenyl-4′-yloxy)-6-(4-cyanobiphenyl-4′-yl)hexane (CB6OCB) containing a single -${\mathrm{CD}}_2$- group. The sample has two nematic liquid crystal phases: a twist-bend nematic, $N_{TB}$, at the lowest temperature followed by a uniaxial nematic, $N_U$, on increasing the temperature. Proton decoupling reduces the linewidths of the peaks in the deuterium spectrum from kHz to $\sim 100\mathrm{Hz}$, enabling quadrupolar splittings, $\mathrm{\Delta }\nu$, to be obtained with enhanced precision as well as the dipolar coupling between deuterium nuclei within the ${\mathrm{CD}}_2$ group, hence enhancing the information content.

The Mpemba effect refers to the counterintuitive effect where a system which is initially further from the final steady state equilibrates faster than an identical system that is initially closer. The closeness to the final state is defined in terms of a distance measure. For driven granular systems, the Mpemba effect has been illustrated in terms of an ad hoc measure of mean kinetic energy as the distance function. In this paper, by studying four different distance measures based on the mean kinetic energies as well as velocity distribution, we show that the Mpemba effect depends on the definition of the measures.

The numerical determination of solitary states is an important topic for such research areas as Bose-Einstein condensates, nonlinear optics, plasma physics, and so on. In this paper, we propose a data-driven approach for identifying solitons based on dynamical solutions of real-time differential equations. Our approach combines a machine-learning architecture called the complex-valued neural operator (CNO) with an energy-restricted gradient optimization. The CNO serves as a generalization of the traditional neural operator to the complex domain, and constructs a smooth mapping between the initial and final states; the energy-restricted optimization facilitates the search for solitons by constraining the energy space. We concretely demonstrate this approach on the quasi-one-dimensional Bose-Einstein condensate with homogeneous and inhomogeneous nonlinearities. Our work offers an idea for data-driven effective modeling and studies of solitary waves in nonlinear physical systems.

The thermodynamic uncertainty relation (TUR) provides a universal entropic bound for the precision of the fluctuation of the charge transfer, for example, for a class of continuous-time stochastic processes. However, its extension to general nonequilibrium dynamics is still an unsolved problem. We derive TUR for an arbitrary finite time from exchange fluctuation theorem under a geometric necessary and sufficient condition. We also generally show a necessary and sufficient condition of multidimensional TUR in a unified manner. As a nontrivial practical consequence, we obtain universal scaling relations among the mean and variance of the charge transfer in short time regime. In this manner, we can deepen our understanding of a link between two important rigorous relations, i.e., the fluctuation theorem and the thermodynamic uncertainty relation.

This paper provides an analysis of the time evolution of a many-particle system starting out of equilibrium with its control parameter fixed at a value corresponding to a many-body energy-level crossing (degeneracy). We prove theorems concerning ergodicity, equilibration, and thermalization. For certain conditions, the occupancy of symmetrically equivalent basis states has different time-averaged probabilities. This nonergodicity remains in equilibrium. If the symmetrically equivalent states have opposite parity in relation to some physical property, then a left and right particle number imbalance averaged in time is nonzero. This imbalance does not occur for all initial basis states. In addition, the Hilbert space of the system is not fragmented; however, there is a subspace spanned by favored basis states, where the system is most likely to be found. Therefore, our results reveal what appears to be a unique mechanism for a weak eigenstates-thermalization-hypothesis breakdown, where the degenerate eigenstates can work as nonthermal eigenstates. To illustrate these findings, we consider the Hubbard Hamiltonian. In this case, ergodicity breaking produces a left and right magnetization imbalance, where the time-averaged probability of finding a spin-$\sigma$ electron on one side of the crystal lattice is greater than on the other side. This imbalance is not associated with electrical charge; thus the conductance is preserved. The potential use in technology is discussed.

We study the stationary states of an overdamped active Brownian particle (ABP) in a harmonic trap in two dimensions via mathematical calculations and numerical simulations. In addition to translational diffusion, the ABP self-propels with a certain velocity, whose magnitude is constant, but its direction is subject to Brownian rotation. In the limit where translational diffusion is negligible, the stationary distribution of the particle's position shows a transition between two different shapes, one with maximum and the other with minimum density at the center, as the trap stiffness is increased. We show that this nonintuitive behavior is captured by the relevant Fokker-Planck equation, which, under minimal assumptions, predicts a continuous phase transition-like change between the two different shapes. As the translational diffusion coefficient is increased, both these distributions converge into the equilibrium, Boltzmann form. Our simulations support the analytical predictions and also show that the probability distribution of the orientation angle of the self-propulsion velocity undergoes a transition from unimodal to bimodal forms in this limit. We also extend our simulations to a three-dimensional trap and find similar behavior.

We present a detailed study by Monte Carlo simulations and finite-size scaling analysis of the phase diagram and ordered bulk phases for the three-dimensional Blume-Capel antiferromagnet in the space of temperature and magnetic and crystal fields (or two chemical potentials in an equivalent lattice-gas model with two particle species and vacancies). The phase diagram consists of surfaces of second- and first-order transitions that enclose a “volume” of ordered phases in the phase space. At relatively high temperatures, these surfaces join smoothly along a line of tricritical points, and at zero magnetic field we obtain good agreement with known values for tricritical exponent ratios [Y. Deng and H. W. J. Blöte, Phys. Rev. E 70, 046111 (2004)]. In limited field regions at lower temperatures (symmetric under reversal of the magnetic field), the tricritical line for this three-dimensional model bifurcates into lines of critical endpoints and critical points, connected by a surface of weak first-order transitions inside the region of ordered phases. This phenomenon is not seen in the two-dimensional version of the same model. We confirm the location of the bifurcation as previously reported [Y.-L. Wang and J. D. Kimel, J. Appl. Phys. 69, 6176 (1991)], and we identify the phases separated by this first-order surface as antiferromagnetically (three-dimensional checker-board) ordered with different vacancy densities. We visualize the phases by real-space snapshots and by structure factors in the three-dimensional space of wave vectors.

Swarmalators are phase oscillators that cluster in space, like fireflies flashing in a swarm to attract mates. Interactions between particles, which tend to synchronize their phases and align their motion, decrease with the distance and phase difference between them, coupling the spatial and phase dynamics. In this work, we explore the effects of inducing phase frustration on a system of swarmalators that move on a one-dimensional ring. Our model is inspired by the well-known Kuramoto-Sakaguchi equations. We find, numerically and analytically, the ordered and disordered states that emerge in the system. The active states, not present in the model without frustration, resemble states found previously in numerical studies for the two-dimensional swarmalators system. One of these states, in particular, shows similarities to turbulence generated in a flattened media. We show that all ordered states can be generated for any values of the coupling constants by tuning the phase frustration parameters only. Moreover, many of these combinations display multistability.

We study the $(1+1)$ focusing nonlinear Schrödinger equation for an initial condition with compactly supported parabolic profile and phase depending quadratically on the spatial coordinate. In the absence of dispersion, using the natural class of self-similar solutions, we provide a criterion for blowup in finite time, generalizing a result by Talanov et al. In the presence of dispersion, we numerically show that the same criterion determines, even beyond the semiclassical regime, whether the solution relaxes or develops a high-order rogue wave, whose onset time is predicted by the corresponding dispersionless catastrophe time. The sign of the chirp appears to determine the prevailing scenario among two competing mechanisms for rogue wave formation. For negative values, the numerical simulations are suggestive of the dispersive regularization of a gradient catastrophe described by Bertola and Tovbis for a different class of smooth, bell-shaped initial data. As the chirp becomes positive, the rogue wave seems to result from the interaction of counterpropagating dispersive dam break flows, as in the box problem recently studied by El, Khamis, and Tovbis. As the chirp and amplitude of the initial profile are relatively easy to manipulate in optical devices and water tank wave generators, we expect our observation to be relevant for experiments in nonlinear optics and fluid dynamics.

Compartmental models are the most widely used framework for modeling infectious diseases. These models have been continuously refined to incorporate all the realistic mechanisms that can shape the course of an epidemic outbreak. Building on a compartmental model that accounts for early detection and isolation of infectious individuals through testing, in this article we focus on the viability of detection processes under limited availability of testing resources, and we study how the latter impacts on the detection rate. Our results show that, in addition to the well-known epidemic transition at ${\mathcal{R}}_0=1$, a second transition occurs at ${\mathcal{R}}_0^{★}>1$ pinpointing the collapse of the detection system and, as a consequence, the switch from a regime of mitigation to a regime in which the pathogen spreads freely. We characterize the epidemic phase diagram of the model as a function of the relevant control parameters: the basic reproduction number, the maximum detection capacity of the system, and the fraction of individuals in shelter. Our analysis thus provides a valuable tool for estimating the detection resources and the level of confinement needed to face epidemic outbreaks.

Active filamentlike systems propelling along their backbone exist across scales ranging from motor-driven biofilaments to worms and robotic chains. In macroscopic active filaments such as a chain of robots, in contrast to their microscopic counterparts, inertial effects on their motion cannot be ignored. Nonetheless, the consequences of the interplay between inertia and flexibility on the shape and dynamics of active filaments remain unexplored. Here we examine inertial effects on a flexible tangentially driven active polymer model pertinent to the above examples and we determine the conditions under which inertia becomes important. Performing Langevin dynamics simulations of active polymers with underdamped and overdamped dynamics for a wide range of contour lengths and activities, we uncover striking inertial effects on conformation and dynamics for high levels of activities. Inertial collisions increase the persistence length of active polymers and remarkably alter their scaling behavior. In stark contrast to passive polymers, inertia leaves its fingerprint at long times by an enhanced diffusion of the center of mass. We rationalize inertia-induced enhanced dynamics by analytical calculations of center-of-mass velocity correlations, applicable to any active polymer model, which reveal significant contributions from active force fluctuations convoluted by inertial relaxation.

We analyze the mixing properties of a floating stirrer driven electromagnetically in a thin layer of electrolyte, consisting of two free-floating magnets with opposite polarities connected by a rigid coupling. The magnetic rotor is set in circular motion using Lorentz forces created due to the interaction of the magnetic field of the rotor with dc currents actuated in logic sequence. We identify a coherent structure similar to a tripolar vortex whose central vortex rotates in the same direction of the rotor promoting chaotic mixing of the fluid in the laminar regime ($\mathrm{Re}=45$). Dyed water visualization and particle image velocimetry were performed to characterize experimentally the mixing and flow dynamics at the surface of the electrolyte layer. A quasitwo-dimensional numerical simulation based on the immersed boundary method, which incorporates the fluid-solid interaction and reproduces the experimental observations satisfactorily, was carried out. Optimal mixing conditions are determined through the exponential growth of the material interfaces, which are established mainly by varying the distance separating the magnets of the rotor.

Nonthermal velocity distributions with much greater tails than a Maxwellian have been observed for radical atoms in plasmas for a long time. Historically, such velocity distributions have been modeled by a two-temperature Maxwell distribution. In this paper, I propose a model based on collisional energy cascade, which has been studied in the field of granular materials. In the collisional energy cascade, a particle ensemble undergoes energy input at the high-energy region, entropy production by elastic collisions among particles, and energy dissipation. For radical atoms, energy input may be caused by the Franck-Condon energy of molecular dissociation or charge-exchange collision with hot ions, and the input energy is eventually dissipated by collisions with the walls. I show that the steady-state velocity distribution in the collisional energy cascade is approximated by the generalized Mittag-Leffler distribution, which is a one-parameter extension of the Maxwell distribution. This parameter indicates the degree of the nonthermality and is related to the relative importance of energy dissipation over entropy production. This model is compared with a direct molecular dynamics simulation for simplified gaseous systems with energy input, as well as some experimentally observed velocity distributions of light radicals in plasmas.

Recent validation experiments on laser irradiation of polymer foils with and without implanted golden nanoparticles are discussed. First we analyze characteristics of craters, formed in the target after its interaction with the laser beam. Preliminary experimental results show significant production of deuterons when both the energy of laser pulse and concentration of nanoparticles are high enough. We consider the deuteron production via the nuclear transmutation reactions $p+C\to d+X$ where protons are accelerated by the Coulomb field generated in the target plasma. We argue that maximal proton energy can be above threshold values for these reactions and the deuteron yield may noticeably increase due to presence of nanoparticles.