We study the Brownian motion of a particle in a bounded circular two-dimensional domain in search for a stationary target on the boundary of the domain. The process switches between two modes: one where it performs a two-dimensional diffusion inside the circle and one where it diffuses along the one-dimensional boundary. During the process, the Brownian particle resets to its initial position with a constant rate $r$. The Fokker-Planck formalism allows us to calculate the mean time to absorption (MTA) as well as the optimal resetting rate for which the MTA is minimized. From the derived analytical results the parameter regions where resetting reduces the search time can be specified. We also provide a numerical method for the verification of our results.

The assumption that the local order parameter is related to an underlying spatially smooth auxiliary field, $u(\vec{r},t)$, is a common feature in theoretical approaches to non-conserved order parameter phase separation dynamics. In particular, the ansatz that $u(\vec{r},t)$ is a Gaussian random field leads to predictions for the decay of the autocorrelation function which are consistent with observations, but distinct from predictions using alternative theoretical approaches. In this paper, the auxiliary field is obtained directly from simulations of the time-dependent Ginzburg-Landau equation in two and three dimensions. The results show that $u(\vec{r},t)$ is equivalent to the distance to the nearest interface. In two dimensions, the probability distribution, $P\left(u\right)$, is well approximated as Gaussian except for small values of $u/L\left(t\right)$, where $L\left(t\right)$ is the characteristic length-scale of the patterns. The behavior of $P\left(u\right)$ in three dimensions is more complicated; the non-Gaussian region for small $u/L\left(t\right)$ is much larger than that in two dimensions but the tails of $P\left(u\right)$ begin to approach a Gaussian form at intermediate times. However, at later times, the tails of the probability distribution appear to decay faster than a Gaussian distribution.

We show that non-Markovian effects of the reservoirs can be used as a resource to extract work from an Otto cycle. The state transformation under non-Markovian dynamics is achieved via a two-step process, namely an isothermal process using a Markovian reservoir followed by an adiabatic process. From second law of thermodynamics, we show that the maximum amount of extractable work from the state prepared under the non-Markovian dynamics quantifies a lower bound of non-Markovianity. We illustrate our ideas with an explicit example of non-Markovian evolution.

We study simple models of intermittency, involving switching between two states, within the dynamical large-deviation formalism. Singularities appear in the formalism when switching is cooperative or when its basic time scale diverges. In the first case the unbiased trajectory distribution undergoes a symmetry breaking, leading to a change in shape of the large-deviation rate function for a particular dynamical observable. In the second case the symmetry of the unbiased trajectory distribution remains unbroken. Comparison of these models suggests that singularities of the dynamical large-deviation formalism can signal the dynamical equivalent of an equilibrium phase transition but do not necessarily do so.

A random search process in a networked environment is governed by the time it takes to visit every node, termed the cover time. Often, a networked process does not proceed in isolation but competes with many instances of itself within the same environment. A key unanswered question is how to optimize this process: How many concurrent searchers can a topology support before the benefits of parallelism are outweighed by competition for space? Here, we introduce the searcher-averaged parallel cover time (APCT) to quantify these economies of scale. We show that the APCT of the networked symmetric exclusion process is optimized at a searcher density that is well predicted by the spectral gap. Furthermore, we find that nonequilibrium processes, realized through the addition of bias, can support significantly increased density optima. Our results suggest alternative hybrid strategies of serial and parallel search for efficient information gathering in social interaction and biological transport networks.

Concentrated colloidal dispersions have been regarded as models for molecular glasses. One of the many ways to compare the behavior in these two different systems is by comparing the structural recovery or the physical aging behavior. However, recent investigations from our group to examine structural recovery in thermosensitive colloidal dispersions have shown contrasting results between the colloidal and the molecular glasses. The differences in the behaviors of the two systems have led us to pose this question: Is structural recovery behavior in colloidal glasses truly distinct from that of molecular glasses or is the conventional experimental condition (isobaric temperature-jumps) in determining the structural recovery in molecular glasses different from the experimental condition in the colloidal experiments (concentration- or volume fraction-jumps); i.e., are colloidal glasses inherently different from molecular glasses or not? To address the question, we resort to model calculations of structural recovery in a molecular glass under constant volume (isochoric) conditions following temperature only- and simultaneous volume- and temperature-jumps, which are closer to the volume fraction-jump conditions used in the thermosensitive-colloidal experiments. The current model predictions are then compared with the signatures of structural recovery under the conventional isobaric state in a molecular glass and with structural recovery behavior in colloidal glasses following volume fraction-jumps. We show that the results obtained from the experiments conducted by our group were contrasting to classical molecular glass behavior because the basis of our comparisons were incorrect (the histories were not analogous). The present calculations (with analogous histories) are qualitatively closer to the colloidal behavior. The signatures of “intrinsic isotherms” and “asymmetry of approach” in the current isochoric model predictions are quite different from those in the classical isobaric conditions while the “memory” signatures remain essentially the same. While there are qualitative similarities between the current isochoric model predictions and results from colloidal glasses, it appears from the calculations that the origins of these are different. The isochoric histories in the molecular glasses have compensating effects of pressure and departure from equilibrium which determines the structure dependence on mobility of the molecules. On the other hand, in the colloids it simply appears that the volume fraction-jump conditions simply do not exhibit such structure mobility dependence. The determining interplay of thermodynamic phase variables in colloidal and molecular systems might be very different or at least their correlations are yet to be ascertained. This topic requires further investigation to bring the similarities and differences between molecular and colloidal glass formers into fuller clarity.

Synchronization processes are ubiquitous despite the many connectivity patterns that complex systems can show. Usually, the emergence of synchrony is a macroscopic observable; however, the microscopic details of the system, as, e.g., the underlying network of interactions, is many times partially or totally unknown. We already know that different interaction structures can give rise to a common functionality, understood as a common macroscopic observable. Building upon this fact, here we propose network transformations that keep the collective behavior of a large system of Kuramoto oscillators invariant. We derive a method based on information theory principles, that allows us to adjust the weights of the structural interactions to map random homogeneous in-degree networks into random heterogeneous networks and vice versa, keeping synchronization values invariant. The results of the proposed transformations reveal an interesting principle; heterogeneous networks can be mapped to homogeneous ones with local information, but the reverse process needs to exploit higher-order information. The formalism provides analytical insight to tackle real complex scenarios when dealing with uncertainty in the measurements of the underlying connectivity structure.

We develop a theoretical model to describe two-beam energy exchange in a hybrid photorefractive cholesteric cell with a short-pitch helix oriented parallel to the cell substrates (so-called uniformly lying helix configuration). Weak and strong light beams incident on the hybrid cell interfere and induce a periodic space-charge field in the photorefractive substrate of the cell, which penetrates into the cholesteric liquid crystal (LC). Due to the flexoelectro-optic effect an interaction of the photorefractive field with the LC flexopolarization causes the spatially periodic modulation of the helix axis in the plane parallel to the cell substrates. Coupling of a weak signal beam with a strong pump beam at the LC permittivity grating, induced by the periodically tilted helix axis, leads to the energy gain of the weak signal beam. Dependence of the signal beam gain coefficient on the parameters of the short-pitch cholesteric LC is studied.

A plasma photonic crystal consists of a plasma density grating which is created in underdense plasma by counterpropagating laser beams. When a high-power laser pulse impinges the crystal, it might be reflected or transmitted. So far only one type of pulse polarization, namely the so-called $s$ wave (or TE mode) was investigated (when the electric field vector is perpendicular to the plane of incidence). Here, when investigating also so-called $p$ waves (or TM modes, where the magnetic field vector is perpendicular to the plane of incidence), it is detected that the transmission and reflection properties of the plasma photonic crystal depend on polarization. A simple analytic model of the crystal allows one to make precise predictions. The first conclusion is that in some operational regime the crystal can act as a plasma polarizer for high-intensity laser pulses. Also, differences in phase velocities for grazing incidence between $s$ and $p$ polarization are found. Thus, secondly, the crystal can be utilized as a waveplate, e.g., transforming linearly polarized laser light into circular polarization. All these processes extend to laser intensities beyond the damage intensities of so far used solid state devices.

The paper presents a derivation of analytical components of $S$ matrices for arbitrary planar diffractive structures and metasurfaces in the Fourier domain. The attained general formulas for $S$-matrix components can be applied within both formulations in the Cartesian and curvilinear metric. A numerical method based on these results can benefit from all previous improvements of the Fourier domain methods. In addition, we provide expressions for $S$-matrix calculation in the case of periodically corrugated layers of two-dimensional materials, which are valid for arbitrary corrugation depth-to-period ratios. As an example, the derived equations are used to simulate resonant grating excitation of graphene plasmons and the impact of a silica interlayer on corresponding reflection curves.

This paper considers two nonequilibrium model systems described by linear Fokker-Planck equations for the time-dependent velocity distribution functions that yield steady state Kappa distributions for specific system parameters. The first system describes the time evolution of a charged test particle in a constant temperature heat bath of a second charged particle. The time dependence of the distribution function of the test particle is given by a Fokker-Planck equation with drift and diffusion coefficients for Coulomb collisions as well as a diffusion coefficient for wave-particle interactions. A second system involves the Fokker-Planck equation for electrons dilutely dispersed in a constant temperature heat bath of atoms or ions and subject to an external time-independent uniform electric field. The momentum transfer cross section for collisions between the two components is assumed to be a power law in reduced speed. The time-dependent Fokker-Planck equations for both model systems are solved with a numerical finite difference method and the approach to equilibrium is rationalized with the Kullback-Leibler relative entropy. For particular choices of the system parameters for both models, the steady distribution is found to be a Kappa distribution. Kappa distributions were introduced as an empirical fitting function that well describe the nonequilibrium features of the distribution functions of electrons and ions in space science as measured by satellite instruments. The calculation of the Kappa distribution from the Fokker-Planck equations provides a direct physically based dynamical approach in contrast to the nonextensive entropy formalism by Tsallis [J. Stat. Phys. 53, 479 (1988)].

Absorption refrigerators transfer thermal energy from a cold bath to a hot bath without input power by utilizing heat from an additional “work” reservoir. Particularly interesting is a three-level design for a quantum absorption refrigerator, which can be optimized to reach the maximal (Carnot) cooling efficiency. Previous studies of three-level chillers focused on the behavior of the averaged cooling current. Here, we go beyond that and study the full counting statistics of heat exchange in a three-level chiller model. We explain how to obtain the complete cumulant generating function of the refrigerator in a steady state, then derive a partial cumulant generating function, which yields closed-form expressions for both the averaged cooling current and its noise. Our analytical results and simulations are beneficial for the design of nanoscale engines and cooling systems far from equilibrium, with their performance optimized according to different criteria, efficiency, power, fluctuations, and dissipation.

We elucidate the impact of diffusive motion on the nature of anomalous dispersion in layered and fibrous disordered media. We consider two types of disorder characterized by quenched random velocities and quenched random retardation properties. Purely advective particle motion is ballistic in both disorder models. This changes dramatically in the presence of transverse diffusion, which leads to dimension-dependent disorder sampling. For $d\le 3$ dimensions, heavy-tailed velocity distributions render large-scale particle motion a correlated Lévy flight, while transport in the quenched random retardation model behaves as a biased continuous time random walk with correlated time increments.

In its classical version, the theory of large deviations makes quantitative statements about the probability of outliers when estimating time averages, if time series data are identically independently distributed. We study large-deviation probabilities (LDPs) for time averages in short- and long-range correlated Gaussian processes and show that long-range correlations lead to subexponential decay of LDPs. A particular deterministic intermittent map can, depending on a control parameter, also generate long-range correlated time series. We illustrate numerically, in agreement with the mathematical literature, that this type of intermittency leads to a power law decay of LDPs. The power law decay holds irrespective of whether the correlation time is finite or infinite, and hence irrespective of whether the central limit theorem applies or not.

We present a general dynamic finite-size scaling theory for the quantum dynamics after an abrupt quench, at both continuous and first-order quantum transitions. For continuous transitions, the scaling laws are naturally ruled by the critical exponents and the renormalization-group dimension of the perturbation at the transition. In the case of first-order transitions, it is possible to recover a universal scaling behavior, which is controlled by the size behavior of the energy gap between the lowest-energy levels. We discuss these findings in the framework of the paradigmatic quantum Ising ring, and support the dynamic scaling laws by numerical evidence.

Aspects of isomorph theory, Rosenfeld-Tarazona temperature scaling, and thermodynamic geometry are comparatively discussed on the basis of the Lennard-Jones potential. The first two approaches approximate the high-density fluid state well when the repulsive interparticle interactions become dominant, which is typically the case close to the freezing line. However, previous studies of Rosenfeld-Tarazona scaling for the isochoric heat capacity and its relation to isomorph theory reveal deviations for the temperature dependence. It turns out that a definition of a state region in which repulsive interactions dominate is required for achieving consistent results. The Riemannian thermodynamic scalar curvature $R$ allows for such a classification, indicating predominantly repulsive interactions by $R>0$. An analysis of the isomorphic character of the freezing line and the validity of Rosenfeld-Tarazona temperature scaling show that these approaches are consistent only in a small state region.

In this paper, we investigate the coupled Sasa-Satsuma equations, which describe the simultaneous propagation of two ultrashort pulses in the birefringent or two-mode fiber with the third-order dispersion, self-steepening, and stimulated Raman scattering effects. Darboux-dressing transformation is applied to obtain the dark-bright soliton and semirational rogue-wave solutions. Dark-bright one solitons with the single-hump, double-hump, and even breather-like structures are presented. Interactions between the double-peak breather and different kinds of dark-bright solitons are studied. We show that the double-peak (or single-peak) rogue wave can coexist and interact with different kinds of dark-bright solitons. Coexistence of the solitons with different velocities and rogue waves is also found. Numerical stabilities of the dark-bright solitons and semirational rogue waves are exhibited. It is expected that those localized wave phenomena can be experimentally observed and have potential applications.

We present a set of phase-space portraits illustrating the extraordinary oscillatory possibilities of the dynamical systems through the so-called generalized Landau scenario. In its simplest form the scenario develops in $N$ dimensions around a saddle-node pair of fixed points experiencing successive Hopf bifurcations up to exhausting their stable manifolds and generating $N-1$ different limit cycles. The oscillation modes associated with these cycles extend over a wide phase-space region by mixing ones within the others and by affecting both the transient trajectories and the periodic orbits themselves. A mathematical theory covering the mode-mixing mechanisms is lacking, and our aim is to provide an overview of their main qualitative features in order to stimulate research on it.