We study the effects of the Lifshitz transition from closed to open Fermi surface in dirty topological insulators with the nematic superconductivity near the critical temperature. We solve linearized Gor'kov equations and find that the nematic superconductor with an open Fermi surface has a lower critical temperature and more susceptible to the disorder than the superconductor with the closed Fermi surface. We propose that correspondence between the critical temperature and stability against the disorder is the general feature of the superconductivity. We investigate the effects of the Lifshitz transition on the competition between superconducting phases in a topological insulator. Open Fermi surface is beneficial for the nematic order parameter $\Delta_4$ in competition with orbital-triplet $\Delta_2$ and disfavors nematic state over the s-wave order parameter. We study Meissner currents in both clean and dirty limits. We found that transition from closed to open Fermi surface increases anisotropy of Meissner currents. Finite disorder suppresses superconducting density stronger than critical temperature. We compare our results with the existing experimental data.

Harrow-Hassidim-Lloyd algorithm (HHL) allows for the exponentially faster solution of a system of linear equations. However, this algorithm requires the postselection of an ancilla qubit to obtain the solution. This postselection makes the algorithm result probabilistic. Here we show conditions when the HHL algorithm can work without postselection of ancilla qubit. We derive expectation values for an observable $M$ on the HHL outcome state when ancilla qubit is measured in $\ket{0}$ and $\ket{1}$ and show condition for postselection-free HHL running. We provide an explicit example of a practically-interesting input matrix and an observable, which satisfy postselection-free HHL condition. Our work can improve the performance of the HHL-based algorithms.

It has been shown that the Kohn--Luttinger superconductivity mechanism interplaying with other types of ordering can be implemented in systems with a hexagonal lattice. A number of unusual properties of such systems in the normal phase have also been considered. Our previous results on Kohn--Luttinger superconductivity with $p$-, $d$-, and $f$-wave pairing in monolayer and AB bilayer graphene, obtained disregarding the effect of substrate potential and impurities, have been presented in the first part. Then, the interplay of the superconducting Kohn--Luttinger state with the spin density wave state in actual AB, AA, and twisted bilayer graphene has been discussed in detail. In the last parts, a number of anomalous properties in the normal phase and the appearance of nematic superconductivity alongside with the spin density wave in the twisted bilayer graphene have been presented.

We study the effects of the Lifshitz transition from closed to open Fermi surface in dirty topological insulators with the nematic superconductivity near the critical temperature. We solve linearized Gor'kov equations and find that the nematic superconductor with an open Fermi surface has a lower critical temperature and more susceptible to the disorder than the superconductor with the closed Fermi surface. We propose that correspondence between the critical temperature and stability against the disorder is the general feature of the superconductivity. We investigate the effects of the Lifshitz transition on the competition between superconducting phases in a topological insulator. Open Fermi surface is beneficial for the nematic order parameter $\Delta_4$ in competition with orbital-triplet $\Delta_2$ and disfavors nematic state over the s-wave order parameter. We study Meissner currents in both clean and dirty limits. We found that transition from closed to open Fermi surface increases anisotropy of Meissner currents. Finite disorder suppresses superconducting density stronger than critical temperature. We compare our results with the existing experimental data.

A family of conservative schemes for the axisymmetric contact smoothed particle hydrodynamics (CSPH) method, which ensure the accuracy and stability in modeling of complex multi-material flows of compressible media, is introduced. Among these schemes, the most convenient ones are considered. Simulations with the proposed schemes may be also improved by embedding of MUSCL reconstruction into a numerical scheme, or by correcting the kernel gradient as was proposed earlier for the Cartesian case. Verification of the proposed method is performed on several test problems: Sod's cylindrical test, Taylor bar test, and Sedov's point explosion. The conservative properties of the scheme are demonstrated. Finally, a set of simulations on air shock wave weakening by a breakaway sand barrier is performed and compared to experimental results.

The ecologically best way to produce nanoparticles (NP) is based on laser ablation in liquid (LAL). In the considered here case the LAL means that a gold target is irradiated through transparent water. During and after irradiation the heated material from surface of a target forms a plume which expands into liquid. In this paper we study a reach set of physical processes mixed with complicated hydrodynamic phenomena which all accompany LAL. These theoretical and simulation investigations are very important for practical applications. Laser pulses with different durations $\tau_L$ covering 5-th orders of magnitudes range from 0.1 ps to 0.5 ns and large absorbed fluences $F_{abs}$ near optical breakdown of liquid are compared. It is shown that the trajectory of the contact boundary with liquid at the middle and late stages after passing of the instant of maximum intensity of the longest pulse are rather similar for very different pulse durations (of course at comparable energies $F_{abs});$ we consider the pulses with a Gaussian temporal shape $I\propto \exp(-t^2/\tau_L^2).$ We follow how hot (few eV range) dense gold plasma expands, cools down, intersects a saturation curve, and condenses into NPs. These NPs appear first inside the water-gold diffusively mixed intermediate layer where gold vapor has the lowest temperature. Later in time pressure around the gold-water contact drops down below critical pressure for water. Thus NPs find themselves in gaseous water bubble where density of water gradually decreases to $10^{-4}-10^{-5}$ g/cm$\!^3$ at the instant of maximum expansion of a bubble.

We consider a quasi-PT-symmetric system of two resonators, one of which interacts with a finite-size environment. The interaction with the environment leads to energy losses in the resonators, and the finite size of the environment leads to a non-Markovian dynamics of the relaxation process. We demonstrate that non-Markovian processes in the quasi-PT-symmetric system can make the states of the system infinitely living, loss-protected states, even in the absence of gain. There is a critical value of the interaction between the resonator and the environment below which any state of the system is loss-protected. When the interaction magnitude is greater than the critical value, depending on the coupling strength between the resonators, either one or both states are unprotected. We show that the boundaries of regions with different numbers of protected states are determined by the relaxation rates in the quasi-PT-symmetric system, calculated in the Markovian approximation. By changing the coupling strength between the resonators and the interaction magnitude between the resonator and the environment, the system switches between modes with two, one, or no loss-protected states. This makes it possible to realize stable PT-symmetric devices based on purely dissipative systems. The obtained results are applicable to quantum systems with single excitations, allowing the concept of PT symmetry to be extended to such systems.

The quantum regression theorem is a powerful tool for calculating the muli-time correlators of operators of open quantum systems which dynamics can be described in Markovian approximation. It enables to obtain the closed system of equation for the multi-time correlators. However, the scope of the quantum regression theorem is limited by a particular time order of the operators in multi-time correlators and does not include out-of-time-ordered correlators. In this work, we obtain an adjoint master equation for multi-time correlators that is applicable to out-of-time-ordered correlators. We show that this equation can be derived for various approaches to description of the dynamics of open quantum systems, such as the global or local approach. We show that the adjoint master equation for multi-time correlators is self-consistent. Namely, the final equation does not depend on how the operators are grouped inside the correlator, and it coincides with the quantum regression theorem for the particular time ordering of the operators.

We develop a microscopic theory for the dynamics of quantum fluids of light, deriving an effective kinetic equation in momentum space that takes the form of the convection-diffusion equation. In the particular case of two-dimensional systems with parabolic dispersion, it reduces to the Bateman--Burgers equation. The hydrodynamic analogy unifies nonlinear wave phenomena, such as shock wave formation and turbulence, with non-equilibrium Bose--Einstein condensation of photons and polaritons in optical cavities. We introduce the Reynolds number $(\textit{Re})$ and demonstrate that the condensation threshold corresponds exactly to a critical Reynolds number of unity $(\textit{Re}=1)$, beyond which $(\textit{Re} > 1)$ a shock-like front emerges in the momentum space, characterized by the Bose--Einstein distribution for the particle density in states with high momentum.

A nonlocal quantum-field model is constructed for the system of hydrodynamic equations for incompressible viscous fluid (the stochastic Navier--Stokes (NS) equation and the continuity equation). This model is studied by the following two mutually parallel methods: the Wilson--Polchinski functional renormalization group method (FRG), which is based on the exact functional equation for the generating functional of amputated connected Green's functions (ACGF), and the Heemskerk--Polchinski holographic renormalization group method (HRG), which is based on the functional Hamilton--Jacobi (HJ) equation for the holographic boundary action. Both functional equations are equivalent to infinite hierarchies of integro-differential equations (coupled in the FRG case) for the corresponding families of Green's functions (GF). The RG-flow equations can be derived explicitly for two-particle functions. Because the HRG-flow equation is closed (contains only a two-particle GF), the explicit analytic solutions are obtained for the two-particle GF (in terms of the modified Bessel functions $I$ and $K$) in the framework of the minimal holographic model and its simple generalization, and these solutions have a remarkable property of minimal dependence on the details of the random force correlator (the function of the energy pumping into the system). The restrictions due to the time-gauged Galilean symmetry present in this theory, the problem of choosing the pumping function, and some generalizations of the standard RG-flow procedures are discussed in detail. Finally, the question of whether the HRG-solutions can be used to solve the FRG-flow equation for the two-particle GF (in particular, the relationship between the regulators in the two methods) is studied.

We study the zero-temperature many-body properties of twisted bilayer graphene with a twist angle equal to the so-called `first magic angle'. The system low-energy single-electron spectrum consists of four (eight, if spin label is accounted) weakly-dispersing partially degenerate bands, each band accommodating one electron per Moir{\'{e}} cell per spin projection. This weak dispersion makes electrons particularly susceptible to the effects of interactions. Introducing several excitonic order parameters with spin-density-wave-like structure, we demonstrate that (i)~the band degeneracy is partially lifted by the interaction, and (ii)~the details of the low-energy spectrum becomes doping-dependent. For example, at or near the undoped state, interactions separate the eight bands into two quartets (one quartet is almost filled, the other is almost empty), while for two electrons per Moir\'{e} cell, the quartets are pulled apart, and doublets emerge. When the doping is equal to one or three electrons per cell, the doublets split into singlets. Hole doping produces similar effects. As a result, electronic properties (e.g., the density of states at the Fermi energy) demonstrate oscillating dependence on the doping concentration. This allows us to reproduce qualitatively the behavior of the conductance observed recently in experiments [Cao et al., Nature {\bf 556}, 80 (2018)]. Near half-filling, the electronic spectrum loses hexagonal symmetry indicating the appearance of a many-body nematic state.

We numerically study the Kitaev honeycomb model with the additional XX Ising interaction between the nearest and the next nearest neighbors (Kitaev-Ising-$J_1$-$J_2$ model), by using the density matrix renormalization group (DMRG) method. Such additional interaction correspond to the nearest and diagonal interactions on the square lattice. Phase diagram of the bare Kitaev model consist of low entangled commensurate magnetic phases and entagled Kitaev spin liquid. Anisotropic Ising interaction allows the entangled incommensurate magnetic phases in the phase diagram, which previously was predicted only for more complex type of interactions. We study the scaling law of the entanglement entropy and the bond dimension of the matrix product state with the size of the system. In addition, we propose an optimization algorithm to prevent DMRG from getting stuck in the low-entangled phases.

Ultrafast emission rates obtained from quantum emitters coupled to plasmonic nanoantennas have recently opened fundamentally new possibilities in quantum information and sensing applications. Plasmonic nanoantennas greatly improve the brightness of quantum emitters by dramatically shortening their fluorescence lifetimes. Gap plasmonic nanocavities that support strongly confined modes are of particular interest for such applications. We demonstrate single-photon emission from nitrogen-vacancy (NV) centers in nanodiamonds coupled to nanosized gap plasmonic cavities with internal mode volumes about 10 000 times smaller than the cubic vacuum wavelength. The resulting structures features sub-nanosecond NV excited-state lifetimes and detected photon rates up to 50 million counts per second. Analysis of the fluorescence saturation allows the extraction of the multi-order excitation rate enhancement provided by the nanoantenna. Efficiency analysis shows that the NV center is producing up to 0.25 billion photons per second in the far-field.

In frames of the nonlocal and nonpolynomial quantum theory of the one component scalar field in $D$-dimensional spacetime, stated by Gariy Vladimirovich Efimov, the expansion of the $\mathcal{S}$-matrix is revisited for different interaction Lagrangians and for some kinds of Gaussian propagators modified by different ultraviolet form factors $F$ which depend on some length parameter $l$. The expansion of the $\mathcal{S}$-matrix is of the form of a grand canonical partition function of some $D+N$-dimensional ($N\geq 1$) classical gas with interaction. The toy model of the realistic quantum field theory (QFT) is considered where the $\mathcal{S}$-matrix is calculated in closed form. Then, the functional Schwinger-Dyson and Schrödinger equations for the $\mathcal{S}$-matrix in Efimov representation are derived. These equations play a central role in the present paper. The functional Schwinger-Dyson and Schrödinger equations in Efimov representation do not involve explicit functional derivatives but involve a shift of the field which is the $\mathcal{S}$-matrix argument. The asymptotic solutions of the Schwinger-Dyson equation are obtained in different limits. Also, the solution is found in one heuristic case allowing us to study qualitatively the behavior of the $\mathcal{S}$-matrix for an arbitrary finite value of its argument. Self-consistency equations, which arise during the process of derivation, are of a great interest. Finally, in the light of the discussion of QFT functional equations, ultraviolet form factors and extra dimensions, the connection with functional (in terms of the Wilson-Polchinski and Wetterich-Morris functional equations) and holographic renormalization groups (in terms of the functional Hamilton-Jacobi equation) is made. In addition the Hamilton-Jacobi equation is formulated in an unconventional way.

Superconducting fluxonium qubits provide a promising alternative to transmons on the path toward large-scale superconductor-based quantum computing due to their better coherence and larger anharmonicity. A major challenge for multi-qubit fluxonium devices is the experimental demonstration of a scalable crosstalk-free multi-qubit architecture with high fidelity single-qubit and two-qubit gates, single-shot readout and state initialization. Here, we present a two-qubit fluxonium-based quantum processor with a tunable coupler element following our theoretical proposal [DOI: https://doi.org/10.1063/5.0064800]. We experimentally demonstrate fSim-type and controlled-Z gates with $99.55\%$ and $99.23\%$ fidelities, respectively. The residual ZZ interaction is suppressed down to the few kHz level. Using a galvanically coupled flux control line, we implement high fidelity single-qubit gates and ground state initialization with a single arbitrary waveform generator channel per qubit.

We discover an unexpected behavior in a hybrid system composed of cavity strongly coupled to molecules and subjected to high intensity coherent pumping. We show that if the frequency of the pumping wave is close to polariton transitions in the hybrid system, non monotone dependence of fluorescence and scattering amplitudes on pump intensity with a narrow resonance-like response occurs. We demonstrate that this phenomenon occurs due to hybridization of lower and upper polaritons with substantially different excitation numbers caused by the pumping affecting polariton states of the system. This occurs when the Rabi interaction with coherent field is comparable to the field-matter coupling constant which in turn needs to be sufficient for the manifestation of strong coupling. This non-monotonic dependence of the fluorescence and scattering amplitudes on pump rate intensity pave the way for creation of nonlinear optical devices.

Polariton thermalization is a key process in achieving light-matter Bose--Einstein condensation, spanning from solid-state semiconductor microcavities at cryogenic temperatures to surface plasmon nanocavities with molecules at room temperature. Originated from the matter component of polariton states, the microscopic mechanisms of thermalization are closely tied to specific material properties. In this work, we investigate polariton thermalization in strongly-coupled molecular systems. We develop a microscopic theory addressing polariton thermalization through electron-phonon interactions (known as exciton-vibration coupling) with low-energy molecular vibrations. This theory presents a simple analytical method to calculate the temperature-dependent polariton thermalization rate, utilizing experimentally accessible spectral properties of bare molecules, such as the Stokes shift and temperature-dependent linewidth of photoluminescence, in conjunction with well-known parameters of optical cavities. Our findings demonstrate qualitative agreement with recent experimental reports of nonequilibrium polariton condensation in both ground and excited states, and explain the thermalization bottleneck effect observed at low temperatures. This study showcases the significance of vibrational degrees of freedom in polariton condensation and offers practical guidance for future experiments, including the selection of suitable material systems and cavity designs.

We study the effects of magnetization on the properties of the doped topological insulator with nematic superconductivity. We found that the direction of the in-plane magnetization fixes the direction of the nematicity in the system. The chiral state is more favorable than the nematic state for large values of out-of-plane magnetization. Overall, the critical temperature of the nematic state is resilient against magnetization. We explore the spectrum of the system with the pinned direction of the nematic order parameter $\Delta_{y}$ in details. Without magnetization, there is a full gap in the spectrum. At strong enough out-of-plane $m_z$ or orthogonal in-plane $m_x$ magnetization, the spectrum is closed at the nodal points that are split by the magnetization. Flat Majorana surface states connect such split bulk nodal points. Parallel magnetization $m_y$ lifts nodal points and opens a full gap in the spectrum. We discuss relevant experiments and propose experimental verifications of our theory.

We report an experimental and numerical investigation of the fragmentation mechanisms of micrometer-sized metal droplet irradiated by ultrashort laser pulses. The results of the experiment show that the fast one-side heating of such a droplet may lead to either symmetric or asymmetric expansion followed by different fragmentation scenarios. To unveil the underlying processes leading to fragmentation we perform simulation of liquid-tin droplet expansion produced by the initial conditions similar to those in experiment using the smoothed particle hydrodynamics (SPH) method. Simulation demonstrates that a thin heated surface layer generates a ultrashort shock wave propagating from the frontal side to rear side of the droplet. Convergence of such shock wave followed by a rarefaction tale to the droplet center results in the cavitation of material inside the central region by the strong tensile stress. Reflection of the shock wave from the rear side of droplet produces another region of highly stretched material where the spallation may occur producing a thin spallation layer moving with a velocity higher than expansion of the central shell after cavitation. It is shown both experimentally and numerically that the threshold laser intensity necessary for the spallation is higher than the threshold required to induce cavitation in the central region of droplet. Thus, the regime of asymmetrical expansion is realized if the laser intensity exceeds the spallation threshold. The transverse and longitudinal expansion velocities obtained in SPH simulations of different regimes of expansion are agreed well with our experimental data.