It is a well established notion that animals can detect the Earth's magnetic field, while the biophysical origin of such magnetoreception is still elusive. Recently, a magnetic receptor Drosophila CG8198 (MagR) with a rodlike protein complex is reported [S. Qin et al., Nat. Mater. 15, 217 (2016)] to act like a compass needle to guide the magnetic orientation of animals. This view, however, is challenged [M. Meister, Elife 5, e17210 (2016)] by arguing that thermal fluctuations beat the Zeeman coupling of the proteins's magnetic moment with the rather weak geomagnetic field ($\sim 25–65$ $\mu \mathrm{T}$). In this work, we show that the spin-mechanical interaction at the atomic scale gives rise to a high blocking temperature which allows a good alignment of the protein's magnetic moment with the Earth's magnetic field at room temperature. Our results provide a promising route to resolve the debate on the thermal behaviors of MagR, and may stimulate a broad interest in spin-mechanical couplings down to atomistic levels.

Tilted smectic liquid crystal phases such as the smectic-C phase seen in calamitic liquid crystals are usually treated using the assumption of biaxial orthorhombic symmetry. However, the smectic-C phase has monoclinic symmetry, thereby allowing disassociation of the principal optic and dielectric axes based on symmetry and invariance principles. This is demonstrated here by comparing optical and dielectric measurements for two materials with highly first-order direct transitions from nematic to smectic-C phases. The results show a high difference between the orientations of the principal axes sets, which is interpreted as the existence of two distinct cone angles for optical and dielectric frequencies. Both materials exhibit an increasing degree of monoclinic behavior with decreasing temperature. Due to fast switching speeds, ferroelectric smectic-C* materials are important for fast modulators and LCoS devices, where the dielectric biaxiality influences device operation.

A biaxial surface potential ${\mathrm{\Phi }}_s$ of smectic-$C^{*}$ surface-stabilized ferroelectric liquid crystals (SSFLCs) is introduced in this paper to explain the experimentally observed electric-field dependence of polarization ${\tilde{P}}_{\mathrm{cell}}\left(E\right)$, in particular the shape of the static hysteresis loops. Our potential consists of three independent parts. The first nonpolar part ${\mathrm{\Phi }}_n$ describes the deviation of the prime director $n$ (which is the most probable orientation of the long molecular axes) from the easy alignment axis $R$, which is located in the boundary surface plane. It is introduced in the same manner as the uniaxial Rapini potential. The second part ${\mathrm{\Phi }}_p$ of the potential is a polar term associated with the presence of the polar axis in a FLC. The third part ${\mathrm{\Phi }}_m$ relates to the inherent FLC biaxiality, which has not been taken into consideration previously. The ${\mathrm{\Phi }}_m$ part takes into account the deviations of the secondary director $m$ (which is the most probable orientation of the short molecular axes) from the normal to the boundary surface. The overall surface potential ${\mathrm{\Phi }}_s$, which is a sum of ${\mathrm{\Phi }}_n,{\mathrm{\Phi }}_p$, and ${\mathrm{\Phi }}_m$, allows one to model the conditions when either one, two, or three minima of the SSFLC cell free energy are realized depending on the biaxiality extent. A monodomain or polydomain structure, as well as the bistability or monostability of SSFLC cells, depends on the number of free-energy minima, as confirmed experimentally. In this paper, we analyze the biaxiality impact on the FLC alignment. We also answer the question of whether the bistable or monostable structure can be formed in an SSFLC cell. Our approach is essentially based on a consideration of the biaxial surface potential, while the uniaxial surface potential cannot adequately describe the experimental observations in the FLC.

We provide a theoretical framework to analyze the properties of frontal collisions of two growing interfaces considering different short-range interactions between them. Due to their roughness, the collision events spread in time and form rough domain boundaries, which defines collision interfaces in time and space. We show that statistical properties of such interfaces depend on the kinetics of the growing interfaces before collision, but are independent of the details of their interaction and of their fluctuations during the collision. Those properties exhibit dynamic scaling with exponents related to the growth kinetics, but their distributions may be nonuniversal. Our results are supported by simulations of lattice models with irreversible dynamics and local interactions. Relations to first passage processes are discussed and a possible application to grain-boundary formation in two-dimensional materials is suggested.

We investigate the effect of a time-reversal-breaking impurity term (of strength ${\lambda }_d$) on both the equilibrium and nonequilibrium critical properties of entanglement entropy (EE) in a three-spin-interacting transverse Ising model, which can be mapped to a $p$-wave superconducting chain with next-nearest-neighbor hopping and interaction. Importantly, we find that the logarithmic scaling of the EE with block size remains unaffected by the application of the impurity term, although, the coefficient (i.e., central charge) varies logarithmically with the impurity strength for a lower range of ${\lambda }_d$ and eventually saturates with an exponential damping factor $[\sim exp(-{\lambda }_d)]$ for the phase boundaries shared with the phase containing two Majorana edge modes. On the other hand, it receives a linear correction in term of ${\lambda }_d$ for an another phase boundary. Finally, we focus to study the effect of the impurity in the time evolution of the EE for the critical quenching case where the impurity term is applied only to the final Hamiltonian. Interestingly, it has been shown that for all the phase boundaries, contrary to the equilibrium case, the saturation value of the EE increases logarithmically with the strength of impurity in a certain regime of ${\lambda }_d$ and finally, for higher values of ${\lambda }_d$, it increases very slowly dictated by an exponential damping factor. The impurity-induced behavior of EE might bear some deep underlying connection to thermalization.

Work fluctuations and work probability distributions are fundamentally different in systems with short-ranged versus long-ranged correlations. Specifically, in systems with long-ranged correlations the work distribution is extraordinarily broad compared to systems with short-ranged correlations. This difference profoundly affects the possible applicability of fluctuation theorems like the Jarzynski fluctuation theorem. The Heisenberg ferromagnet, well below its Curie temperature, is a system with long-ranged correlations in very low magnetic fields due to the presence of Goldstone modes. As the magnetic field is increased the correlations gradually become short ranged. Hence, such a ferromagnet is an ideal system for elucidating the changes of the work probability distribution as one goes from a domain with long-ranged correlations to a domain with short-ranged correlations by tuning the magnetic field. A quantitative analysis of this crossover behavior of the work probability distribution and the associated fluctuations is presented.

The relevance of the concept of Fisher information is increasing in both statistical physics and quantum computing. From a statistical mechanical standpoint, the application of Fisher information in the kinetic theory of gases is characterized by its decrease along the solutions of the Boltzmann equation for Maxwellian molecules in the two-dimensional case. From a quantum mechanical standpoint, the output state in Grover's quantum search algorithm follows a geodesic path obtained from the Fubini-Study metric on the manifold of Hilbert-space rays. Additionally, Grover's algorithm is specified by constant Fisher information. In this paper, we present an information geometric characterization of the oscillatory or monotonic behavior of statistically parametrized squared probability amplitudes originating from special functional forms of the Fisher information function: constant, exponential decay, and power-law decay. Furthermore, for each case, we compute both the computational speed and the availability loss of the corresponding physical processes by exploiting a convenient Riemannian geometrization of useful thermodynamical concepts. Finally, we briefly comment on the possibility of using the proposed methods of information geometry to help identify a suitable trade-off between speed and thermodynamic efficiency in quantum search algorithms.

We study transitions in log-correlated random energy models (logREMs) that are related to the violation of a Seiberg bound in Liouville field theory (LFT): the binding transition and the termination point transition (a.k.a., pre-freezing). By means of LFT-logREM mapping, replica symmetry breaking and traveling-wave equation techniques, we unify both transitions in a two-parameter diagram, which describes the free-energy large deviations of logREMs with a deterministic background log potential, or equivalently, the joint moments of the free energy and Gibbs measure in logREMs without background potential. Under the LFT-logREM mapping, the transitions correspond to the competition of discrete and continuous terms in a four-point correlation function. Our results provide a statistical interpretation of a peculiar nonlocality of the operator product expansion in LFT. The results are rederived by a traveling-wave equation calculation, which shows that the features of LFT responsible for the transitions are reproduced in a simple model of diffusion with absorption. We examine also the problem by a replica symmetry breaking analysis. It complements the previous methods and reveals a rich large deviation structure of the free energy of logREMs with a deterministic background log potential. Many results are verified in the integrable circular logREM, by a replica-Coulomb gas integral approach. The related problem of common length (overlap) distribution is also considered. We provide a traveling-wave equation derivation of the LFT predictions announced in a precedent work.

Strong-coupling statistical thermodynamics is formulated as the Hamiltonian dynamics of an observed system interacting with another unobserved system (a bath). It is shown that the entropy production functional of stochastic thermodynamics, defined as the log ratio of forward and backward system path probabilities, is in a one-to-one relation with the log ratios of the joint initial conditions of the system and the bath. A version of strong-coupling statistical thermodynamics where the system-bath interaction vanishes at the beginning and at the end of a process is, as is also weak-coupling stochastic thermodynamics, related to the bath initially in equilibrium by itself. The heat is then the change of bath energy over the process, and it is discussed when this heat is a functional of the system history alone. The version of strong-coupling statistical thermodynamics introduced by Seifert and Jarzynski is related to the bath initially in conditional equilibrium with respect to the system. This leads to heat as another functional of the system history which needs to be determined by thermodynamic integration. The log ratio of forward and backward system path probabilities in a stochastic process is finally related to log ratios of the initial conditions of a combined system and bath. It is shown that the entropy production formulas of stochastic processes under a general class of time reversals are given by the differences of bath energies in a larger underlying Hamiltonian system. The paper highlights the centrality of time reversal in stochastic thermodynamics, also in the case of strong coupling.

The generalized nonparaxial nonlinear Schrödinger (NLS) equation in optical fibers filled with chiral materials is reduced to the higher-order integrable Hirota equation. Based on the modified Darboux transformation method, the nonparaxial chiral optical rogue waves are constructed from the scalar model with modulated coefficients. We show that the parameters of nonparaxiality, third-order dispersion, and differential gain or loss term are the main keys to control the amplitude, linear, and nonlinear effects in the model. Moreover, the influence of nonparaxiality, optical activity, and walk-off effect are also evidenced under the defocusing and focusing regimes of the vector nonparaxial NLS equations with constant and modulated coefficients. Through an algorithm scheme of wider applicability on nonparaxial beam propagation methods, the most influential effect and the simultaneous controllability of combined effects are underlined, showing their properties and their potential applications in optical fibers and in a variety of complex dynamical systems.

We derive the transmission probabilities and delay times, and identify quasibound state structures in an open quantum system consisting of three Gaussian potential energy peaks, a system whose classical scattering dynamics we show to be chaotic. Such open quantum systems can serve as models for nanoscale quantum devices and their wave dynamics are similar to electromagnetic wave dynamics in optical microcavities. We use a quantum web to determine energy regimes for which the system exhibits the quantum manifestations of chaos, and we show that the classical scattering dynamics contains a significant amount of chaos. We also derive an exact expression for the non-Hermitian Hamiltonian whose eigenvalues give quasibound state energies and lifetimes of the system.

The spread of an infectious disease is known to change people's behavior, which in turn affects the spread of disease. Adaptive network models that account for both epidemic and behavioral change have found oscillations, but in an extremely narrow region of the parameter space, which contrasts with intuition and available data. In this paper we propose a simple susceptible-infected-susceptible epidemic model on an adaptive network with time-delayed rewiring, and show that oscillatory solutions are now present in a wide region of the parameter space. Altering the transmission or rewiring rates reveals the presence of an endemic bubble—an enclosed region of the parameter space where oscillations are observed.

The dynamics of a pair of identical artificial microswimmers bound inside two harmonic traps, in a thin sheared fluid film, is numerically investigated. In a two-dimensional Oseen approximation, the hydrodynamic pair coupling is long-ranged and proportional to the particle radius to film thickness ratio. On increasing such ratio above a certain threshold, a transition occurs between a free regime, where each swimmer orbits in its own trap with random phase, and a strong synchronization regime, where the two swimmers strongly repel each other to an average distance larger than both the trap distance and their free orbit diameter. Moreover, the swimmers tend to synchronize their positions opposite the center of the system.

We investigate Marangoni instability in a thin liquid film resting on a substrate of low thermal conductivity and separated from the surrounding gas phase by a deformable free surface. Considering a nonmonotonic variation of surface tension with temperature, here we analytically derive the neutral stability curve for the monotonic and oscillatory modes of instability (for both the long-wave and short-wave perturbations) under the framework of linear stability analysis. For the long-wave instability, we derive a set of amplitude equations using the scaling $k\sim {\left(\mathrm{Bi}\right)}^{1/2}$, where $k$ is the wave number and Bi is the Biot number. Through this investigation, we demonstrate that for such a fluid layer upon heating from below, both monotonic and oscillatory instability can appear for a certain range of the dimensionless parameters, viz., Biot number $\left(\mathrm{Bi}\right)$, Galileo number $\left(\mathrm{Ga}\right)$, and inverse capillary number $\left(\mathrm{\Sigma }\right)$. Moreover, we unveil, through this study, the influential role of the above-mentioned parameters on the stability of the system and identify the critical values of these parameters above which instability initiates in the liquid layer.

The linear theory of Rayleigh-Taylor instability is developed for the case of a viscous fluid layer accelerated by a semi-infinite viscous fluid, considering that the top interface is a free surface. Effects of the surface tensions at both interfaces are taken into account. When viscous effects dominate on surface tensions, an interplay of two mechanisms determines opposite behaviors of the instability growth rate with the thickness of the heavy layer for an Atwood number $A_T=1$ and for sufficiently small values of $A_T$. In the former case, viscosity is a less effective stabilizing mechanism for the thinnest layers. However, the finite thickness of the heavy layer enhances its viscous effects that, in general, prevail on the viscous effects of the semi-infinite medium.

The flow of nematic liquid crystals (NLCs) in the presence of an electric field is typically characterized by the variation in its rheological properties due to transition in its molecular arrangements. Here, we bring out a nontrivial interplay of a consequent alteration in the resistive viscous effects and driving electrocapillary interactions, toward maneuvering the capillary filling dynamics over miniaturized scales. Considering a dynamic interplay of the relevant bulk and interfacial forces acting in tandem, our results converge nicely to previously reported experimental data. Finally, we attempt a scaling analysis to bring forth further insight to the reported observations. Our analysis paves the way for the development of microfluidic strategies with previously unexplored paradigms of interaction between electrical and fluidic phenomenon, providing with an augmented controllability on capillary filling as compared to tthose reported to be achievable by the existing strategies. This, in turn, holds utilitarian scopes in improved designs of functional capillarities in electro-optical systems, electrorheological utilities, electrokinetic flow control, as well as in interfacing and imaging systems for biomedical microdevices.

We introduce two-color counterrotating circularly polarized laser fields as a way to spectrally control high harmonic generation (HHG) from relativistic plasma mirrors. Through particle-in-cell simulations, we show that only a selected group of harmonic orders can appear owing to the symmetry of the laser fields and the related conservation laws. By adjusting the intensity ratio of the two driving field components, we demonstrate the overall HHG efficiency, the relative intensity of allowed neighboring harmonic orders, and that the polarization state of the harmonic source can be tuned. The HHG efficiency of this scheme can be as high as that driven by a linearly polarized laser field.

Spatiotemporal thermal response and characteristics of net entropy production rate of a gold nanosphere (radius: 50–200 nm), subjected to a short-pulse, femtosecond laser is reported. In order to correctly illustrate the temperature history of laser-metal interaction(s) at picoseconds transient with a comprehensive single temperature definition in macroscale and to further understand how the thermophysical response of the single-phase lag (SPL) and dual-phase lag (DPL) frameworks (with various lag-ratios’) differs, governing energy equations derived from these benchmark non-Fourier frameworks are numerically solved and thermodynamic assessment under both the classical irreversible thermodynamics (CIT) as well as extended irreversible thermodynamics (EIT) frameworks is subsequently carried out. Under the frameworks of SPL and DPL with small lag ratio, thermophysical anomalies such as temperature overshooting characterized by adverse temperature gradient is observed to violate the local thermodynamic equilibrium (LTE) hypothesis. The EIT framework, however, justifies the compatibility of overshooting of temperature with the second law of thermodynamics under a nonequilibrium paradigm. The DPL framework with higher lag ratio was however observed to remain free from temperature overshooting and finds suitable consistency with LTE hypothesis. In order to solve the dimensional non-Fourier governing energy equation with volumetric laser-irradiation source term(s), the lattice Boltzmann method (LBM) is extended and a three-time level, fully implicit, second order accurate finite difference method (FDM) is illustrated. For all situations under observation, the LBM scheme is featured to be computationally superior to remaining FDM schemes. With detailed prediction of maximum temperature rise and the corresponding peaking time by all the numerical schemes, effects of the change of radius of the gold nanosphere, the magnitude of fluence of laser, and laser irradiation with multiple pulses on thermal energy transport and lagging behavior (if any) are further elucidated at different radial locations of the gold nanosphere. Last, efforts are further made to address the thermophysical characteristics when effective thermal conductivity (with temporal and size effects) is considered instead of the usual bulk thermal conductivity.

A first analysis of fluctuations of the light intensity of vertical-cavity surface-emitting lasers operating in a bistable regime reveals the presence of $1/f$ noise. In this regime the intensity fluctuates between two recently characterized states with residence times $\left\{{\tau }_1\right\}$ and $\left\{{\tau }_2\right\}$. We identify three distinct processes. One of them presents a coherence enhancement phenomenon, and in the other two the distribution of residence times in one of the states follows either a power law $P\left({\tau }_1\right)\sim {\tau }_1^{-2}$ or $P\left({\tau }_2\right)\sim {\tau }_2^{-2}$, and this is the cause of the $1/f$ shape in the spectral density of the intensity. The process at the coherence enhancement zone shows $1/f$ fluctuations in the light intensity and also in the time residence process. It is shown that the origin of these fluctuations is due to a power-law distribution in the time separation between pulses observed in the time residence series.

Here we study the relaxation of a chain consisting of three masses joined by nonlinear springs and periodic conditions when the stiffness is weakened. This system, when expressed in their normal coordinates, yields a softened Henon-Heiles system. By reducing the stiffness of one low-frequency vibrational mode, a faster relaxation is enabled. This is due to a reduction of the energy barrier heights along the softened normal mode as well as for a widening of the opening channels of the energy landscape in configurational space. The relaxation is for the most part exponential, and can be explained by a simple flux equation. Yet, for some initial conditions the relaxation follows as a power law, and in many cases there is a regime change from exponential to power-law decay. We pinpoint the initial conditions for the power-law decay, finding two regions of sticky states. For such states, quasiperiodic orbits are found since almost for all components of the initial momentum orientation, the system is trapped inside two pockets of configurational space. The softened Henon-Heiles model presented here is intended as the simplest model in order to understand the interplay of rigidity, nonlinear interactions and relaxation for nonequilibrium systems such as glass-forming melts or soft matter. Our softened system can be applied to model $\beta$ relaxation in glasses and suggest that local reorientational jumps can have an exponential and a nonexponential contribution for relaxation, the latter due to asymmetric molecules sticking in cages for certain orientations.