We investigate the system of two coupled one-dimensional Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) equations which possess the global conservation law. Such system of equations has been recently derived for the quasiparticle densities in the two-band fermionic model with the particle-number conserving dissipative protocol. As standard FKPP equation the studied system of equations has one unstable and one stable homogeneous solution with travelling wave switching between them. We demonstrate that the conservation law enforces the synchronization of travelling waves for both densities and determine their minimal possible velocity. Surprisingly, we find the existence of jumps of the minimal velocity as function of control parameters. We obtain that the minimal velocity of the coupled FKPP equations may significantly exceed the minimal velocity for a single FKPP equation in a wide range of control parameters.

We present a theoretical study of terahertz radiation-induced transitions between attractive and repulsive Fermi polaron states in monolayers of transition metal dichalcogenides. Going beyond the simple few-particle trion picture, we develop a many-body description that explicitly accounts for correlations with the Fermi sea of resident charge carriers. We calculate the rate of the direct optical conversion process, showing that it features a characteristic frequency dependence near the threshold due to final-state electron-exciton scattering related to the trion correlation with the Fermi sea hole. Furthermore, we demonstrate that intense terahertz pulses can significantly heat the electron gas via Drude absorption enabling an additional, indirect conversion mechanism through collisions between hot electrons and polarons, which exhibits a strong exponential dependence on temperature. Our results reveal the important role of many-body correlations and thermal effects in the terahertz-driven dynamics of excitonic complexes in two-dimensional semiconductors.

We study the statistical properties of the complex generalization of Wigner time delay $\tau_\text{W}$ for sub-unitary wave chaotic scattering systems. We first demonstrate theoretically that the mean value of the $\text{Re}[\tau_\text{W}]$ distribution function for a system with uniform absorption strength $\eta$ is equal to the fraction of scattering matrix poles with imaginary parts exceeding $\eta$. The theory is tested experimentally with an ensemble of microwave graphs with either one or two scattering channels, and showing broken time-reversal invariance and variable uniform attenuation. The experimental results are in excellent agreement with the developed theory. The tails of the distributions of both real and imaginary time delay are measured and are also found to agree with theory. The results are applicable to any practical realization of a wave chaotic scattering system in the short-wavelength limit.

We develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for mass flow and pressure on a metric graph, where turbulent losses of momentum are modeled by phenomenological Darcy-Weisbach friction. Mass flow balance is maintained through the boundary conditions at the network nodes, where natural gas is injected or withdrawn from the system. Gas flow through the network is controlled by compressors boosting pressure at the inlet of the adjoint pipe. Our operator splitting numerical scheme is unconditionally stable and it is second order accurate in space and time. The scheme is explicit, and it is formulated to work with general networks with loops. We test the scheme over range of regimes and network configurations, also comparing its performance with performance of two other state of the art implicit schemes.

We consider disordered tight-binding models which Green's functions obey the self-consistent cavity equations . Based on these equations and the replica representation, we derive an analytical expression for the fractal dimension D_{1} that distinguishes between the extended ergodic, D_{1}=1, and extended non-ergodic (multifractal), 0

We study the response of an isolated quantum system governed by the Hamiltonian drawn from the Gaussian Rosenzweig-Porter random matrix ensemble to a perturbation controlled by a small parameter. We focus on the density of states, local density of states and the eigenfunction amplitude overlap correlation functions which are calculated exactly using the mapping to the supersymmetric nonlinear sigma model. We show that the susceptibility of eigenfunction fidelity to the parameter of perturbation can be expressed in terms of these correlation functions and is strongly peaked at the localization transition: It is independent of the effective disorder strength in the ergodic phase, grows exponentially with increasing disorder in the fractal phase and decreases exponentially in the localized phase. As a function of the matrix size, the fidelity susceptibility remains constant in the ergodic phase and increases in the fractal and in the localized phases at modestly strong disorder. We show that there is a critical disorder strength inside the insulating phase such that for disorder stronger than the critical the fidelity susceptibility decreases with increasing the system size. The overall behavior is very similar to the one observed numerically in a recent work by Sels and Polkovnikov [Phys. Rev. E 104, 054105 (2021)] for the normalized fidelity susceptibility in a disordered XXZ spin chain.

We develop a semi-quantitative theory of electron pairing and resulting superconductivity in bulk "poor conductors" in which Fermi energy $E_F$ is located in the region of localized states not so far from the Anderson mobility edge $E_c$. We review the existing theories and experimental data and argue that a large class of disordered films is described by this model. Our theoretical analysis is based on the analytical treatment of pairing correlations, described in the basis of the exact single-particle eigenstates of the 3D Anderson model, which we combine with numerical data on eigenfunction correlations. Fractal nature of critical wavefunction's correlations is shown to be crucial for the physics of these systems. We identify three distinct phases: 'critical' superconductive state formed at $E_F=E_c$, superconducting state with a strong pseudogap, realized due to pairing of weakly localized electrons and insulating state realized at $E_F$ still deeper inside localized band. The 'critical' superconducting phase is characterized by the enhancement of the transition temperature with respect to BCS result, by the inhomogeneous spatial distribution of superconductive order parameter and local density of states. The major new feature of the pseudo-gaped state is the presence of two independent energy scales: superconducting gap $\Delta$, that is due to many-body correlations and a new "pseudogap" energy scale $\Delta_P$ which characterizes typical binding energy of localized electron pairs and leads to the insulating behavior of the resistivity as a function of temperature above superconductive $T_c$. Two gap nature of the "pseudo-gaped superconductor" is shown to lead to a number of unusual physical properties.

Supervised quantum learning is an emergent multidisciplinary domain bridging between variational quantum algorithms and classical machine learning. Here, we study experimentally a hybrid classifier model accelerated by a quantum simulator - a linear array of four superconducting transmon artificial atoms - trained to solve multilabel classification and image recognition problems. We train a quantum circuit on simple binary and multi-label tasks, achieving classification accuracy around 95%, and a hybrid model with data re-uploading with accuracy around 90% when recognizing handwritten decimal digits. Finally, we analyze the inference time in experimental conditions and compare the performance of the studied quantum model with known classical solutions.

We develop the microscopic theory for the attenuation of out-of-plane phonons in stressed flexible two-dimensional crystalline materials. We demonstrate that the presence of nonzero tension strongly reduces the relative magnitude of the attenuation and, consequently, results in parametrical narrowing of the phononspectral line. We predict the specific power-law dependence of the spectral-line width on temperature and tension. We speculate that suppression of the phonon attenuation by nonzero tension might be responsible for high quality factors of mechanical nanoresonators based on flexural two-dimensional materials.

It has recently been suggested that at the post-inflationary stage of the mixed Higgs-$R^2$ model of inflation efficient particle production can arise from the tachyonic instability of the Higgs field. It might complete the preheating of the Universe if appropriate conditions are satisfied, especially in the Higgs-like regime. In this paper, we study this behavior in more depth, including the conditions for occurrence, analytical estimates for the maximal efficiency, and the necessary degree of fine-tuning among the model parameters to complete preheating by this effect. We find that the parameter sets that cause the most efficient tachyonic instabilities obey simple laws in both the Higgs-like regime and the $R^2$-like regime, respectively. We then estimate the efficiency of this instability. In particular, even in the deep $R^2$-like regime with a small non-minimal coupling, this effect is strong enough to complete preheating although a severe fine-tuning is required among the model parameters. We also estimate how much fine-tuning is needed to complete preheating by this effect. It is shown that the fine-tuning of parameters for the sufficient particle production is at least < \mathcal{O}(0.1) in the deep Higgs-like regime with a large scalaron mass, while it is more severe $\sim {\cal O}(10^{-4})-{\cal O}(10^{-5})$ in the $R^2$-like regime with a small non-minimal coupling.

We consider the static and dynamic phases in a Rosenzweig-Porter (RP) random matrix ensemble with the tailed distribution of off-diagonal matrix elements of the form of the large-deviation ansatz. We present a general theory of survival probability in such a random-matrix model and show that the {\it averaged} survival probability may decay with time as the simple exponent, as the stretch-exponent and as a power-law or slower. Correspondingly, we identify the exponential, the stretch-exponential and the frozen-dynamics phases. As an example, we consider the mapping of the Anderson model on Random Regular Graph (RRG) onto the "multifractal" RP model and find exact values of the stretch-exponent $\kappa$ depending on box-distributed disorder in the thermodynamic limit. As another example we consider the logarithmically-normal RP (LN-RP) random matrix ensemble and find analytically its phase diagram and the exponent $\kappa$. In addition, our theory allows to compute the shift of apparent phase transition lines at a finite system size and show that in the case of RP associated with RRG and LN-RP with the same symmetry of distribution function of hopping, a finite-size multifractal "phase" emerges near the tricritical point which is also the point of localization transition.

We analyze a two-field inflationary model consisting of the Ricci scalar squared ($R^2$) term and the standard Higgs field non-minimally coupled to gravity in addition to the Einstein $R$ term. Detailed analysis of the power spectrum of this model with mass hierarchy is presented, and we find that one can describe this model as an effective single-field model in the slow-roll regime with a modified sound speed. The scalar spectral index predicted by this model coincides with those given by the $R^2$ inflation and the Higgs inflation implying that there is a close relation between this model and the $R^2$ inflation already in the original (Jordan) frame. For a typical value of the self-coupling of the standard Higgs field at the high energy scale of inflation, the role of the Higgs field in parameter space involved is to modify the scalaron mass, so that the original mass parameter in the $R^2$ inflation can deviate from its standard value when non-minimal coupling between the Ricci scalar and the Higgs field is large enough.

We study the statistical properties of the complex generalization of Wigner time delay $\tau_\text{W}$ for sub-unitary wave chaotic scattering systems. We first demonstrate theoretically that the mean value of the $\text{Re}[\tau_\text{W}]$ distribution function for a system with uniform absorption strength $\eta$ is equal to the fraction of scattering matrix poles with imaginary parts exceeding $\eta$. The theory is tested experimentally with an ensemble of microwave graphs with either one or two scattering channels, and showing broken time-reversal invariance and variable uniform attenuation. The experimental results are in excellent agreement with the developed theory. The tails of the distributions of both real and imaginary time delay are measured and are also found to agree with theory. The results are applicable to any practical realization of a wave chaotic scattering system in the short-wavelength limit.

We study quasiparticle energy relaxation at sub-kelvin temperatures by injecting hot electrons into an aluminium island and measuring the energy flux from electrons into phonons both in the superconducting and in the normal state. The data show strong reduction of the flux at low temperatures in the superconducting state, in qualitative agreement with the presented quasiclassical theory for clean superconductors. Quantitatively, the energy flux exceeds that from the theory both in the superconducting and in the normal state, possibly suggesting an enhanced or additional relaxation process.

We extend the calculation of the Unruh effect to the universality classes of quantum vacua obeying topologically protected invariance under anisotropic scaling ${\bf r} \rightarrow b {\bf r}$, $t \rightarrow b^z t$. Two situations are considered. The first one is related to the accelerated detector which detects the electron - hole pairs. The second one is related to the system in external electric field, when the electron - hole pairs are created due to the Schwinger process. As distinct from the Unruh effect in relativistic systems (where $z=1$) the calculated radiation is not thermal, but has properties of systems in the vicinity of quantum criticality. The vacuum obeying anisotropic scaling can be realized, in particular, in multilayer graphene with the rhombohedral stacking. Opportunities of the experimental realization of Unruh effect in this situation are discussed.

We perform full-scale numerical simulation of instability of weakly nonlinear waves on the surface of deep fluid. We show that the instability development leads to chaotization and formation of wave turbulence. We study instability both of propagating and standing waves. We studied separately pure capillary wave unstable due to three-wave interactions and pure gravity waves unstable due to four-wave interactions. The theoretical description of instabilities in all cases is included into the article. The numerical algorithm used in these and many other previous simulations performed by authors is described in details.

We study the electron-phonon relaxation in the model of a granular metal film, where the grains are formed by regularly arranged potential barriers of arbitrary transparency. The relaxation rate of Debye acoustic phonons is calculated, taking into account two mechanisms of electron-phonon scattering: the standard Fröhlich interaction of the lattice deformation with the electron density and the interaction mediated by the displacement of grain boundaries dragged by the lattice vibration. At the lowest temperatures, the electron-phonon cooling power follows the power-law temperature dependence typical for clean systems but with the prefactor growing as the transparency of the grain boundaries decreases.

It is known that scalar-tensor theory of gravity admits regular crossing of the phantom divide line $w_{\DE}=-1$ for dark energy, and existing viable models of present dark energy for its particular case -- $f(R)$ gravity -- possess one such crossing in the recent past, after the end of the matter dominated stage. It was recently noted that during the future evolution of these models the dark energy equation of state $w_{\DE}$ may oscillate with an arbitrary number of phantom divide crossings. In this paper we prove that the number of crossings can be infinite, present an analytical condition for the existence of this effect and investigate it numerically. With the increase of the present mass of the scalaron (a scalar particle appearing in $f(R)$ gravity) beyond the boundary of the appearance of such oscillations, their amplitude is shown to decrease very fast. As a result, the effect quickly becomes small and its beginning is shifted to the remote future.

We develop the pQCD description of the diffraction dissociation (DD) of longitudinal photons. We demonstrate that the longitudinal diffractive structure function does not factor into the flux of pomerons and the partonic structure function of the pomeron, thus defying the usually assumed Regge factorization. In contrast to DD of the transverse photons, DD of the longitudinal photons is strongly peaked at $\beta =1$. We comment on duality properties of DD in deep inelastic scattering.

Weyl semimetal is a solid material with isolated touching points between conduction and valence bands in its Brillouin zone -- Weyl points. Low energy excitations near these points exhibit a linear dispersion and act as relativistic massless particles. Weyl points are stable topological objects robust with respect to most perturbations. We study effects of weak disorder on the spectral and transport properties of Weyl semimetals in the limit of low energies. We use a model of Gaussian white-noise potential and apply dimensional regularization scheme near three dimensions to treat divergent terms in the perturbation theory. In the framework of self-consistent Born approximation, we find closed expressions for the average density of states and conductivity. Both quantities are analytic functions in the limit of zero energy. We also include interference terms beyond the self-consistent Born approximation up to the third order in the disorder strength. These interference corrections are stronger than the mean-field result and non-analytic as functions of energy. Our main result is the dependence of conductivity (in units $e^2/h$) on the electron concentration $\sigma = \sigma_0 - 0.891\, n^{1/3} + 0.115\, (n^{2/3}/\sigma_0) \ln|n|$.