Controlling a time-dependent force applied to single molecules or colloidal particles is crucial for many types of experiments. Since in optical tweezers the primary controlled variable is the position of the trap, imposing a target force requires an active feedback process. We analyze this feedback process for the paradigmatic case of a nonequilibrium steady state generated by a dichotomous force protocol, first theoretically for a colloidal particle in a harmonic trap and then with both simulations and experiments for a long DNA hairpin. For the first setup, we find there is an optimal feedback gain separating monotonic from oscillatory response, whereas a too strong feedback leads to an instability. For the DNA molecule, reaching the target force requires substantial feedback gain since weak feedback cannot overcome the tendency to relax towards the equilibrium force.

It is argued that two-dimensional $\text{U}\left(N\right)$ spin models for any $N$ undergo a Berezinskii-Kosterlitz-Thouless (BKT)-like phase transition, similarly to the famous $XY$ model. This conclusion follows from the Berezinskii-like calculation of the two-point correlation function in $\text{U}\left(N\right)$ models, approximate renormalization group analysis, and numerical investigations of the $\text{U}\left(2\right)$ model. It is shown, via Monte Carlo simulations, that the universality class of the $\text{U}\left(2\right)$ model coincides with that of the $XY$ model. Moreover, preliminary numerical results point out that two-dimensional $\text{SU}\left(N\right)$ spin models with the fundamental and adjoint terms and $N>4$ exhibit two phase transitions of BKT type, similarly to $Z\left(N\right)$ vector models.

We report the emergence of a collective dynamical state, namely, the phase-flip chimera, from an ensemble of identical nonlinear oscillators that are coupled indirectly via the dynamical variables from a common environment, which in turn are nonlocally coupled. The phase-flip chimera is characterized by the coexistence of two adjacent out-of-phase synchronized coherent domains interspersed by an incoherent domain, in which the nearby oscillators are in out-of-phase synchronized states. Attractors of the coherent domains are either from the same or from different basins of attractions, depending on whether they are periodic or chaotic. The conventional chimera precedes the phase-flip chimera in general. Further, the phase-flip chimera emerges after the completely synchronized evolution of the ensemble, in contrast to conventional chimeras, which emerge as an intermediate state between completely incoherent and coherent states. We have also characterized the observed dynamical transitions using the strength of incoherence, probability distribution of the correlation coefficient, and framework of the master stability function.

We study the growth kinetics of glassy correlations in a structural glass by monitoring the evolution, within mode-coupling theory, of a suitably defined three-point function ${\chi }_C(t,t_w)$ with time $t$ and waiting time $t_w$. From the complete wave-vector-dependent equations of motion for domain growth, we pass to a schematic limit to obtain a numerically tractable form. We find that the peak value ${\chi }_C^P$ of ${\chi }_C(t,t_w)$, which can be viewed as a correlation volume, grows as $t_w^{0.5}$, and the relaxation time as $t_w^{0.8}$, following a quench to a point deep in the glassy state. These results constitute a theoretical explanation of the simulation findings of Parisi [J. Phys. Chem. B 103, 4128 (1999)] and Kob and Barrat [Phys. Rev. Lett. 78, 4581 (1997)], and they are also in qualitative agreement with Parsaeian and Castillo [Phys. Rev. E 78, 060105(R) (2008)]. On the other hand, if the quench is to a point on the liquid side, the correlation volume grows to saturation. We present a similar calculation for the growth kinetics in a $p$-spin spin glass mean-field model where we find a slower growth, ${\chi }_C^P\sim t_w^{0.13}$. Further, we show that a shear rate $\dot{\gamma }$ cuts off the growth of glassy correlations when $t_w\sim 1/\dot{\gamma }$ for quench in the glassy regime and $t_w=min(t_r,1/\dot{\gamma })$ in the liquid, where $t_r$ is the relaxation time of the unsheared liquid. The relaxation time of the steady-state fluid in this case is $\propto {\dot{\gamma }}^{-0.8}$.

The linear stability of plane-parallel flow of an incompressible viscous fluid over a saturated porous layer is studied to model the instability of water flow in a river over aquatic plants. The saturated porous layer is bounded from below by a rigid plate and the pure fluid layer has a free, undeformable upper boundary. A small inclination of the layers is imposed to simulate the riverbed slope. The layers are inclined at a small angle to the horizon. The problem is studied within two models: the Brinkman model with the boundary conditions by Ochoa-Tapia and Whitaker at the interface, and the Darcy-Forchheimer model with the conditions by Beavers and Joseph. The neutral curves and critical Reynolds numbers are calculated for various porous layer permeabilities and relative thicknesses of the porous layer. The results obtained within the two models are compared and analyzed.

The merger of two like-signed vortices is a well-studied problem, but in a turbulent flow, we may often have more than two like-signed vortices interacting. We study the merger of three or more identical corotating vortices initially arranged on the vertices of a regular polygon. At low to moderate Reynolds numbers, we find an additional stage in the merger process, absent in the merger of two vortices, where an annular vortical structure is formed and is long lived. Vortex merger is slowed down significantly due to this. Such annular vortices are known at far higher Reynolds numbers in studies of tropical cyclones, which have been noticed to a break down into individual vortices. In the preannular stage, vortical structures in a viscous flow are found here to tilt and realign in a manner similar to the inviscid case, but the pronounced filaments visible in the latter are practically absent in the former. Five or fewer vortices initially elongate radially, and then reorient their long axis closer to the azimuthal direction so as to form an annulus. With six or more vortices, the initial alignment is already azimuthal. Interestingly at higher Reynolds numbers, the merger of an odd number of vortices is found to proceed very differently from that of an even number. The former process is rapid and chaotic whereas the latter proceeds more slowly via pairing events. The annular vortex takes the form of a generalized Lamb-Oseen vortex (GLO), and diffuses inward until it forms a standard Lamb-Oseen vortex. For lower Reynolds number, the numerical (fully nonlinear) evolution of the GLO vortex follows exactly the analytical evolution until merger. At higher Reynolds numbers, the annulus goes through instabilities whose nonlinear stages show a pronounced difference between even and odd mode disturbances. Here again, the odd mode causes an early collapse of the annulus via decaying turbulence into a single central vortex, whereas the even mode disturbance causes a more orderly progression into a single vortex. Results from linear stability analysis agree with the nonlinear simulations, and predict the frequencies of the most unstable modes better than they predict the growth rates. It is hoped that the present findings, that multiple vortex merger is qualitatively different from the merger of two vortices, will motivate studies on how multiple vortex interactions affect the inverse cascade in two-dimensional turbulence.

In this paper we exactly solve, within the grand canonical ensemble, a minimal spin model with the hybrid phase transition. We call the model diffusion based because its Hamiltonian can be recovered from a simple dynamic procedure, which can be seen as an equilibrium statistical mechanics representation of a biased random walk. We outline the derivation of the phase diagram of the model, in which the triple point has the hallmarks of the hybrid transition: discontinuity in the average magnetization and algebraically diverging susceptibilities. At this point, two second-order transition curves meet in equilibrium with the first-order curve, resulting in a prototypical mixed-order behavior.

Lévy walks define a fundamental concept in random walk theory that allows one to model diffusive spreading faster than Brownian motion. They have many applications across different disciplines. However, so far the derivation of a diffusion equation for an $n$-dimensional correlated Lévy walk remained elusive. Starting from a fractional Klein-Kramers equation here we use a moment method combined with a Cattaneo approximation to derive a fractional diffusion equation for superdiffusive short-range auto-correlated Lévy walks in the large time limit, and we solve it. Our derivation discloses different dynamical mechanisms leading to correlated Lévy walk diffusion in terms of quantities that can be measured experimentally.

We study the dynamical evolution toward steady state of the stochastic nonequilibrium model known as the totally asymmetric simple exclusion process, in both uniform and nonuniform (staggered) one-dimensional systems with open boundaries. Domain-wall theory and numerical simulations are used and, where pertinent, their results are compared to existing mean-field predictions and exact solutions where available. For uniform chains we find that the inclusion of fluctuations inherent to the domain-wall formulation plays a crucial role in providing good agreement with simulations, which is severely lacking in the corresponding mean-field predictions. For alternating-bond chains the domain-wall predictions for the features of the phase diagram in the parameter space of injection and ejection rates turn out to be realized only in an incipient and quantitatively approximate way. Nevertheless, significant quantitative agreement can be found between several additional domain-wall theory predictions and numerics.

Rings of oscillators with delayed pulse coupling are studied analytically, numerically, and experimentally. The basic regimes observed in such rings are rotating waves with constant interspike intervals and phase lags between the neighbors. We show that these rotating waves may destabilize leading to the so-called jittering waves. For these regimes, the interspike intervals are no more equal but form a periodic sequence in time. Analytic criterion for the emergence of jittering waves is derived and confirmed by the numerical and experimental data. The obtained results contribute to the hypothesis that the multijitter instability is universal in systems with pulse coupling.

The dynamics of a pair of harmonic oscillators represented by three-dimensional fields coupled with a repulsive cubic nonlinearity is investigated through direct simulations of the respective field equations and with the help of the finite-mode Galerkin approximation (GA), which represents the two interacting fields by a superposition of $3+3$ harmonic-oscillator $p$-wave eigenfunctions with orbital and magnetic quantum numbers $l=1$ and $m=1$, 0, $-1$. The system can be implemented in binary Bose-Einstein condensates, demonstrating the potential of the atomic condensates to emulate various complex modes predicted by classical field theories. First, the GA very accurately predicts a broadly degenerate set of the system's ground states in the $p$-wave manifold, in the form of complexes built of a dipole coaxial with another dipole or vortex, as well as complexes built of mutually orthogonal dipoles. Next, pairs of noncoaxial vortices and/or dipoles, including pairs of mutually perpendicular vortices, develop remarkably stable dynamical regimes, which feature periodic exchange of the angular momentum and periodic switching between dipoles and vortices. For a moderately strong nonlinearity, simulations of the coupled-field equations agree very well with results produced by the GA, demonstrating that the dynamics is accurately spanned by the set of six modes limited to $l=1$.

We develop the theory of sparse multiplex networks with partially overlapping links based on their local treelikeness. This theory enables us to find the giant mutually connected component in a two-layer multiplex network with arbitrary correlations between connections of different types. We find that correlations between the overlapping and nonoverlapping links markedly change the phase diagram of the system, leading to multiple hybrid phase transitions. For assortative correlations we observe recurrent hybrid phase transitions.

Dielectric properties are evaluated by means of molecular dynamics simulations on two model systems made up of dipolar molecules. One of them mimics methanol, whereas the other differs from the former only in the ability to form hydrogen bonds. Static dielectric properties such as the permittivity and the Kirkwood factor are evaluated, and results are analyzed by considering the distribution of relative orientations between molecular dipoles. Dipole moment–time correlation functions are also evaluated. The relevance of contributions associated with autocorrelations of molecular dipoles and with cross-correlations between dipoles belonging to different molecules has been investigated. For methanol, the Debye approximation for the overall dipole moment correlation function is not valid at room temperature. The model applies when hydrogen bonds are suppressed, but it fails upon cooling the nonassociated liquid. Important differences between relaxation times associated with dipole auto- versus cross-correlations as well as their relative relevance are at the root of the Debye model breakdown.

β relaxation, which is predicted by mode coupling theory and involves the localized motions of particles, initiates in a supercooled liquid and continues into glassy state. It correlates essentially with many fundamental properties of amorphous materials. Despite its importance, the underlying mechanisms leading to the β relaxation have remained elusive. As natural heterogeneity, the original distributed free volume has been supposed to be associated with the activation of β relaxation in amorphous solids. However, there has been no direct experimental proof for this hypothesis. Here we used a colloidal glass to directly observe the β relaxation and free-volume distribution. We found a spatial correlation between the β relaxation and free volume. The large free volume regions were observed to possess a low-energy cost of relaxation-induced strain, indicating that the large free volume region presenting a low-energy barrier for structural relaxation benefits the β relaxation.

The objective of this study is to investigate the influence of compressibility on Lagrangian statistics of vorticity and strain-rate interactions. The Lagrangian statistics are extracted from “almost” time-continuous data sets of direct numerical simulations of compressible decaying isotropic turbulence by employing a cubic spline-based Lagrangian particle tracker. We study the influence of compressibility on Lagrangian statistics of alignment in terms of compressibility parameters—turbulent Mach number, normalized dilatation-rate, and flow topology. In comparison to incompressible turbulence, we observe that the presence of compressibility in a flow field weakens the alignment tendency of vorticity toward the largest strain-rate eigenvector. Based on the Lagrangian statistics of alignment conditioned on dilatation and topology, we find that the weakened tendency of alignment observed in compressible turbulence is because of a special group of fluid particles that have an initially negligible dilatation-rate and are associated with stable-focus-stretching topology.

A methodology to simulate thermal fields in complex microflow geometries is proposed. For the flow fields, the regularized multiple-relaxation-time lattice Boltzmann method (LBM) is applied coupled with the diffusive-bounce-back boundary condition for wall boundaries. For the thermal fields, the regularized lattice Bhatnagar-Gross-Krook model is applied. For the thermal wall boundary condition, a newly developed boundary condition, which is a mixture of the diffuse scattering and constant temperature conditions, is applied. The proposed set of schemes is validated by reference data in the Fourier flows and square cylinder flows confined in a microchannel. The obtained results confirm that it is essential to apply the regularization to the thermal LBM for avoiding kinked temperature profiles in complex thermal flows. The proposed wall boundary condition is successful to obtain thermal jumps at the walls with good accuracy.

The problem of synchronization of finite-size chemical oscillators described by active inertial particles is addressed for situations in which they are immersed in a reacting nonstationary chaotic flow. Active substances in the fluid will be modeled by Lagrangian particles closely following the fluid streamlines. Their interaction with the active inertial particles as well as the properties of the fluid dynamics will result in modifying the synchronization state of the chemical oscillators. This behavior is studied in terms of the exchange rate between the Lagrangian and inertial particles, and the finite-time Lyapunov exponents characterizing the flow. The coherence of the population of oscillators is determined by means of the order parameter introduced by Kuramoto. The different dynamics observed for the inertial particles (chemical oscillators) and Lagrangian particles (describing chemicals in the flow) lead to nonlinear interactions and patches of synchronized regions within the domain.

An algebraic solution of the Dirac equation is reinvestigated. Slater-type spinor orbitals and their corresponding system of differential equations are defined in two- and four-component formalism. They describe the radial function in components of the wave function of the Dirac equation solution to high accuracy. They constitute the matrix elements arising in a generalized eigenvalue equation. These terms are evaluated through prolate spheroidal coordinates. The corresponding integrals are calculated by the numerical global-adaptive method taking into account the Gauss-Kronrod numerical integration extension. Sample calculations are performed using flexible basis sets generated with both signs of the relativistic angular momentum quantum number $\kappa$. Applications to one-electron atoms and diatomics are detailed. Variationally optimum values for orbital parameters are obtained at given nuclear separation. Methods discussed in this work are capable of yielding highly accurate relativistic two-center integrals for all ranges of orbital parameters. This work provides an efficient way to overcome the problems that arise in relativistic calculations.

In this Rapid Communication we express the partition function of the integrable elliptic solid-on-solid model with domain-wall boundary conditions as a single determinant. This representation appears naturally as the solution of a system of functional equations governing the model's partition function.