We present exact multiparameter families of soliton solutions for two- and three-component Manakov equations in the defocusing regime. Existence diagrams for such solutions in the space of parameters are presented. Fundamental soliton solutions exist only in finite areas on the plane of parameters. Within these areas, the solutions demonstrate rich spatiotemporal dynamics. The complexity increases in the case of three-component solutions. The fundamental solutions are dark solitons with complex oscillating patterns in the individual wave components. At the boundaries of existence, the solutions are transformed into plain (nonoscillating) vector dark solitons. The superposition of two dark solitons in the solution adds more frequencies in the patterns of oscillating dynamics. These solutions admit degeneracy when the eigenvalues of fundamental solitons in the superposition coincide.

Previous studies of nonlinear oscillator networks have shown that amplitude death (AD) occurs after tuning oscillator parameters and coupling properties. Here, we identify regimes where the opposite occurs and show that a local defect (or impurity) in network connectivity leads to AD suppression in situations where identically coupled oscillators cannot. The critical impurity strength value leading to oscillation restoration is an explicit function of network size and system parameters. In contrast to homogeneous coupling, network size plays a crucial role in reducing this critical value. This behavior can be traced back to the steady-state destabilization through a Hopf's bifurcation, which occurs for impurity strengths below this threshold. This effect is illustrated across different mean-field coupled networks and is supported by simulations and theoretical analysis. Since local inhomogeneities are ubiquitous and often unavoidable, such imperfections can be an unexpected source of oscillation control.

The roughness of a fracture surface records a crack's complex path through a material and can affect the resultant frictional or fluid transport properties of the broken medium. For brittle fractures, some of the most prominent surface features are long, step-like discontinuities called step lines. In heterogeneous materials, the mean crack surface roughness created by these step lines is well captured by a simple, one-dimensional ballistic annihilation model, which assumes the creation of these steps is a random processes with a single probability that depends on the heterogeneity of the material, and that their destruction occurs via pairwise interactions. Here, through an exhaustive study of experimentally generated crack surfaces in brittle hydrogels, we examine step interactions and show that interaction outcomes depend on the geometry of the incoming steps. The rules that govern step interactions can be categorized into three unique classes and are fully described, providing a complete framework for predicting fracture roughness.

Holes in mechanical metamaterials, quasilocalized plastic events in amorphous solids, and bound dislocations in a hexatic matter are different mechanisms of generic stress relaxation in solids. Regardless of the specific mechanism, these and other local stress relaxation modes are quadrupolar in nature, forming the foundation for stress screening in solids, similar to polarization fields in electrostatic media. We propose a geometric theory for stress screening in generalized solids based on this observation. The theory includes a hierarchy of screening modes, each characterized by internal length scales, and is partially analogous to theories of electrostatic screening such as dielectrics and Debye-Hückel theory. Additionally, our formalism suggests that the hexatic phase, traditionally defined by structural properties, can also be defined by mechanical properties and may exist in amorphous materials.

A wave near an isolated turning point is typically assumed to have an Airy function profile with respect to the separation distance. This description is incomplete, however, and is insufficient to describe the behavior of more realistic wave fields that are not simple plane waves. Asymptotic matching to a prescribed incoming wave field generically introduces a phase front curvature term that changes the characteristic wave behavior from the Airy function to that of the hyperbolic umbilic function. This function, which is one of the seven classic “elementary” functions from catastrophe theory along with the Airy function, can be understood intuitively as the solution for a linearly focused Gaussian beam propagating in a linearly varying density profile, as we show. The morphology of the caustic lines that govern the intensity maxima of the diffraction pattern as one alters the density length scale of the plasma, the focal length of the incident beam, and also the injection angle of the incident beam are presented in detail. This morphology includes a Goos-Hänchen shift and focal shift at oblique incidence that do not appear in a reduced ray-based description of the caustic. The enhancement of the intensity swelling factor for a focused wave compared to the typical Airy solution is highlighted, and the impact of a finite lens aperture is discussed. Collisional damping and finite beam waist are included in the model and appear as complex components to the arguments of the hyperbolic umbilic function. The observations presented here on the behavior of waves near turning points should aid the development of improved reduced wave models to be used, for example, in designing modern nuclear fusion experiments.

Many experimentally accessible, finite-sized interacting quantum systems are most appropriately described by the canonical ensemble of statistical mechanics. Conventional numerical simulation methods either approximate them as being coupled to a particle bath or use projective algorithms which may suffer from nonoptimal scaling with system size or large algorithmic prefactors. In this paper, we introduce a highly stable, recursive auxiliary field quantum Monte Carlo approach that can directly simulate systems in the canonical ensemble. We apply the method to the fermion Hubbard model in one and two spatial dimensions in a regime known to exhibit a significant “sign” problem and find improved performance over existing approaches including rapid convergence to ground-state expectation values. The effects of excitations above the ground state are quantified using an estimator-agnostic approach including studying the temperature dependence of the purity and overlap fidelity of the canonical and grand canonical density matrices. As an important application, we show that thermometry approaches often exploited in ultracold atoms that employ an analysis of the velocity distribution in the grand canonical ensemble may be subject to errors leading to an underestimation of extracted temperatures with respect to the Fermi temperature.

First-passage properties of continuous stochastic processes confined in a one-dimensional interval are well described. However, for jump processes (discrete random walks), the characterization of the corresponding observables remains elusive, despite their relevance in various contexts. Here we derive exact asymptotic expressions for the leftward, rightward, and complete exit-time distributions from the interval $[0,x]$ for symmetric jump processes starting from $x_0=0$, in the large $x$ and large time limit. We show that both the leftward probability $F_{\underline{0},x}\left(n\right)$ to exit through 0 at step $n$ and rightward probability $F_{0,\underline{x}}\left(n\right)$ to exit through $x$ at step $n$ exhibit a universal behavior dictated by the large-distance decay of the jump distribution parametrized by the Levy exponent $\mu$. In particular, we exhaustively describe the $n\ll {(x/a_{\mu })}^{\mu }$ and $n\gg {(x/a_{\mu })}^{\mu }$ limits and obtain explicit results in both regimes. Our results finally provide exact asymptotics for exit-time distributions of jump processes in regimes where continuous limits do not apply.

We have studied the performance of a measurement-based quantum Otto engine (QOE) in a working system of two spins coupled by Heisenberg anisotropic interaction. A nonselective quantum measurement fuels the engine. We have calculated thermodynamic quantities of the cycle in terms of the transition probabilities between the instantaneous energy eigenstates, and also between the instantaneous energy eigenstates and the basis states of the measurement, when the unitary stages of the cycle operate for a finite time $\tau$. The efficiency attains a large value in the limit of $\tau \to 0$ and then gradually reaches the adiabatic value in a long-time limit $\tau \to \infty$. For finite values of $\tau$ and for anisotropic interaction, an oscillatory behavior of the efficiency of the engine is observed. This oscillation can be interpreted in terms of interference between the relevant transition amplitudes in the unitary stages of the engine cycle. Therefore, for a suitable choice of timing of the unitary processes in the short time regime, the engine can have a higher work output and less heat absorption, such that it works more efficiently than a quasistatic engine. In the case of an always-on heat bath, in a very short time, the bath has a negligible effect on its performance.

A prominent type of collective dynamics in networks of coupled oscillators is the coexistence of coherently and incoherently oscillating domains known as chimera states. Chimera states exhibit various macroscopic dynamics with different motions of the Kuramoto order parameter. Stationary, periodic and quasiperiodic chimeras are known to occur in two-population networks of identical phase oscillators. In a three-population network of identical Kuramoto-Sakaguchi phase oscillators, stationary and periodic symmetric chimeras were previously studied on a reduced manifold in which two populations behaved identically [Phys. Rev. E 82, 016216 (2010)]. In this paper, we study the full phase space dynamics of such three-population networks. We demonstrate the existence of macroscopic chaotic chimera attractors that exhibit aperiodic antiphase dynamics of the order parameters. We observe these chaotic chimera states in both finite-sized systems and the thermodynamic limit outside the Ott-Antonsen manifold. The chaotic chimera states coexist with a stable chimera solution on the Ott-Antonsen manifold that displays periodic antiphase oscillation of the two incoherent populations and with a symmetric stationary chimera solution, resulting in tristability of chimera states. Of these three coexisting chimera states, only the symmetric stationary chimera solution exists in the symmetry-reduced manifold.

It is well known that tree-based theories can describe the properties of undirected clustered networks with extremely accurate results [S. Melnik et al., Phys. Rev. E 83, 036112 (2011)]. It is reasonable to suggest that a motif-based theory would be superior to a tree one, since additional neighbor correlations are encapsulated in the motif structure. In this paper, we examine bond percolation on random and real world networks using belief propagation in conjunction with edge-disjoint motif covers. We derive exact message passing expressions for cliques and chordless cycles of finite size. Our theoretical model gives good agreement with Monte Carlo simulation and offers a simple, yet substantial improvement on traditional message passing, showing that this approach is suitable to study the properties of random and empirical networks.

External kink modes, believed to be the drive of the $\beta$-limiting resistive wall mode, are strongly stabilized by the presence of a separatrix. We thus propose a novel mechanism explaining the appearance of long-wavelength global instabilities in free boundary high-$\beta$ diverted tokamaks, retrieving the experimental observables within a physical framework dramatically simpler than most of the models employed for the description of such phenomena. It is shown that the magnetohydrodynamic stability is worsened by the synergy of $\beta$ and plasma resistivity, with wall effects significantly screened in an ideal, i.e., with vanishing resistivity, plasma with separatrix. Stability can be improved by toroidal flows, depending on the proximity to the resistive marginal boundary. The analysis is performed in tokamak toroidal geometry, and includes averaged curvature and essential separatrix effects.

Complex system stability can be studied via linear stability analysis using random matrix theory (RMT) or via feasibility (requiring positive equilibrium abundances). Both approaches highlight the importance of interaction structure. Here we show, analytically and numerically, how RMT and feasibility approaches can be complementary. In generalized Lotka-Volterra (GLV) models with random interaction matrices, feasibility increases when predator-prey interactions increase; increasing competition/mutualism has the opposite effect. These changes have crucial impact on the stability of the GLV model.

Although the cooperative dynamics emerging from a network of interacting players has been exhaustively investigated, it is not yet fully understood when and how network reciprocity drives cooperation transitions. In this work, we investigate the critical behavior of evolutionary social dilemmas on structured populations by using the framework of master equations and Monte Carlo simulations. The developed theory describes the existence of absorbing, quasiabsorbing, and mixed strategy states and the transition nature, continuous or discontinuous, between the states as the parameters of the system change. In particular, when the decision-making process is deterministic, in the limit of zero effective temperature of the Fermi function, we find that the copying probabilities are discontinuous functions of the system's parameters and of the network degrees sequence. This may induce abrupt changes in the final state for any system size, in excellent agreement with the Monte Carlo simulation results. Our analysis also reveals the existence of continuous and discontinuous phase transitions for large systems as the temperature increases, which is explained in the mean-field approximation. Interestingly, for some game parameters, we find optimal “social temperatures” maximizing or minimizing the cooperation frequency or density.

Employing the quantum magnetohydrodynamic (QMHD) model, the basic properties of magnetosonic waves were investigated in a magnetorotating quantum plasma. The contemplated system considered a combined effects of quantum tunneling and degeneracy forces, dissipation influence, spin magnetization, in addition to the Coriolis force. Fast and slow magnetosonic modes were obtained and examined in the linear regime. Their frequencies are significantly modified due to the rotating parameters (frequency and angle) in addition to quantum correction effects. The nonlinear Korteweg–de Vries–Burger equation was derived using the reductive perturbation approach in a small amplitude limit. The aspects of magnetosonic shock profiles were explored analytically by applying the Bernoulli equation approach and numerically using the Runge-Kutta method. The regarded plasma parameters due to the investigated effects were found to play major roles in specifying the nature of monotonic and oscillatory shock waves' structures and their features. Our results may be applicable in magnetorotating quantum plasma in astrophysical environments such as neutron stars and white dwarfs.

Size-invariant shape transformation is a technique of changing the shape of a domain while preserving its sizes under the Lebesgue measure. In quantum-confined systems, this transformation leads to so-called quantum shape effects in the physical properties of confined particles associated with the Dirichlet spectrum of the confining medium. Here we show that the geometric couplings between levels generated by the size-invariant shape transformations cause nonuniform scaling in the eigenspectra. In particular, the nonuniform level scaling, in the direction of increasing quantum shape effect, is characterized by two distinct spectral features: lowering of the first eigenvalue (ground-state reduction) and changing of the spectral gaps (energy level splitting or degeneracy formation depending on the symmetries). We explain the ground-state reduction by the increase in local breadth (i.e., parts of the domain becoming less confined) that is associated with the sphericity of these local portions of the domain. We accurately quantify the sphericity using two different measures: the radius of the inscribed $n$-sphere and the Hausdorff distance. Due to Rayleigh-Faber-Krahn inequality, the greater the sphericity, the lower the first eigenvalue. Then level splitting or degeneracy, depending on the symmetries of the initial configuration, becomes a direct consequence of size invariance dictating the eigenvalues to have the same asymptotic behavior due to Weyl law. Such level splittings may be interpreted as geometric analogs of Stark and Zeeman effects. Furthermore, we find that the ground-state reduction causes a quantum thermal avalanche which is the underlying reason for the peculiar effect of spontaneous transitions to lower entropy states in systems exhibiting the quantum shape effect. Unusual spectral characteristics of size-preserving transformations can assist in designing confinement geometries that could lead to classically inconceivable quantum thermal machines.

For the study of symmetry-breaking phenomena in neuronal networks, simplified versions of the FitzHugh-Nagumo model are widely used. In this paper, these phenomena are investigated in a network of FitzHugh-Nagumo oscillators taken in the form of the original model and it is found that it exhibits diverse partial synchronization patterns that are unobserved in the networks with simplified models. Apart from the classical chimera, we report a new type of chimera pattern whose incoherent clusters are characterized by spatial random swings among a few fixed periodic attractors. Another peculiar hybrid state is found that combines the features of this chimera state and a solitary state such that the main coherent cluster is interspersed with some nodes with identical solitary dynamics. In addition, oscillation death including chimera death emerges in this network. A reduced model of the network is derived to study oscillation death, which helps explaining the transition from spatial chaos to oscillation death via the chimera state with a solitary state. This study deepens our understanding of chimera patterns in neuronal networks.

We consider a class of Hubbard-Stratonovich transformations suitable for treating Hubbard interactions in the context of quantum Monte Carlo simulations. A tunable parameter $p$ allows us to continuously vary from a discrete Ising auxiliary field ($p=\infty$) to a compact auxiliary field that couples to electrons sinusoidally ($p=0$). In tests on the single-band square and triangular Hubbard models, we find that the severity of the sign problem decreases systematically with increasing $p$. Selecting $p$ finite, however, enables continuous sampling methods such as the Langevin or Hamiltonian Monte Carlo methods. We explore the tradeoffs between various simulation methods through numerical benchmarks.

The attenuation in the propagation of a plane wave in conducting media has been studied. We analyzed a wave motion suffering dissipation by the Joule effect in its propagation in a medium with global disorder. We solved the stochastic telegrapher's equation in the Fourier-Laplace representation allowing us to find the space penetration length of a plane wave in a complex conducting medium. Considering fluctuations in the loss of energy, we found a critical value $k_c$ for Fourier's modes, thus if $\left|k\right|<k_c$ the waves are localized. We showed that the penetration length is inversely proportional to $k_c$. Thus, the penetration length $L=k_c^{-1}$ becomes an important piece of information for describing wave propagation with Markovian and non-Markovian fluctuations in the rate of the absorption of energy ${\tau }^{-1}$. In addition, intermittent fluctuations in this rate have also been studied.

A simple two-dimensional statistical mechanical water model, called the rose model, was used in this work. We studied how a homogeneous constant electric field affects the properties of water. The rose model is a very simple model that helps explain the anomalous properties of water. Rose water molecules are represented as two-dimensional Lennard-Jones disks with potentials for orientation-dependent pairwise interactions mimicking formations of hydrogen bonds. The original model is modified by addition of charges for interaction with the electric field. We studied what kind of influence the electric field strength has on the model's properties. To determine the structure and thermodynamics of the rose model under the influence of the electric field we used Monte Carlo simulations. Under the influence of a weak electric field the anomalous properties and phase transitions of the water do not change. On the other hand, the strong fields shift the phase transition points as well as the position of the density maximum.

In a two-phase solid, we examine the growth of a preexisting macroscopic crack based on simulations of a random spring network model. We find that the enhancement in toughness, as well as strength, is strongly dependent on the ratio of elastic moduli as well as on the relative proportion of the phases. We find that the mechanism that leads to enhancement in toughness is not the same as that for enhancement in strength; however, the overall enhancement is similar in mode I and mixed-mode loading. Based on the crack paths, and the spread of the fracture process zone, we identify the type of fracture to transition from nucleation type, for close to single-phase compositions, whether hard or soft, to avalanche type for more mixed compositions. We also show that the associated avalanche distributions exhibit power-law statistics with different exponents for each phase. The significance of variations in the avalanche exponents with the relative proportion of phases and possible connections to the fracture types are discussed in detail.