Superfluid Edge Dislocation: Transverse Quantum Fluid.

Recently, it was argued [Kuklov et al., Phys. Rev. Lett. 128, 255301 (2022)PRLTAO0031-900710.1103/PhysRevLett.128.255301] that unusual features associated with the superflow-through-solid effect observed in solid ^{4}He can be explained by unique properties of dilute distribution of superfluid edge dislocations. We demonstrate that stability of supercurrents controlled by quantum phase slips (instantons), and other exotic infrared properties of the superfluid dislocations readily follow from a one-dimensional quantum liquid distinguished by an effectively infinite compressibility (in the absence of Peierls potential) associated with the edge dislocation's ability to climb. This establishes a new class of quasi-one-dimensional superfluid states that remain stable and long-range ordered despite their dimensionality. Our theory is consistent with the existing experimental data, and we propose an experiment to test the mass-current-pressure characteristic prediction.

Superconducting Transition Temperature of the Bose One-Component Plasma.

We present results of numerically exact simulations of the Bose one-component plasma, i.e., a Bose gas with pairwise Coulomb interactions among particles and a uniform neutralizing background. We compute the superconducting transition temperature for a wide range of densities, in two and three dimensions, for both continuous and lattice versions of the model. The Coulomb potential causes the weakly interacting limit to be approached at high density, but gives rise to no qualitatively different behavior, vis-à-vis the superfluid transition, with respect to short-ranged interactions. Our results are of direct relevance to quantitative studies of bipolaron mechanisms of (high-temperature) superconductivity.

Dynamic Response of an Electron Gas: Towards the Exact Exchange-Correlation Kernel.

Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. We introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using a series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.

Supertransport by Superclimbing Dislocations in

The unique superflow-through-solid effect observed in solid ^{4}He and attributed to the quasi-one-dimensional superfluidity along the dislocation cores exhibits two extraordinary features: (i) an exponentially strong suppression of the flow by a moderate increase in pressure and (ii) an unusual temperature dependence of the flow rate with no analogy to any known system and in contradiction with the standard Luttinger liquid paradigm. Based on ab initio and model simulations, we argue that the two features are closely related: Thermal fluctuations of the shape of a superclimbing edge dislocation induce large, correlated, and asymmetric stress fields acting on the superfluid core. The critical flux is most sensitive to strong rare fluctuations and hereby acquires a sharp temperature dependence observed in experiments.

Anti-drude metal of bosons.

In the absence of frustration, interacting bosons in their ground state in one or two dimensions exist either in the superfluid or insulating phases. Superfluidity corresponds to frictionless flow of the matter field, and in optical conductivity is revealed through a distinct δ-functional peak at zero frequency with the amplitude known as the Drude weight. This characteristic low-frequency feature is instead absent in insulating phases, defined by zero static optical conductivity. Here we demonstrate that bosonic particles in disordered one dimensional chains can also exist in a conducting, non-superfluid, phase when their hopping is of the dipolar type, often viewed as short-ranged in one dimension. This phase is characterized by finite static optical conductivity, followed by a broad anti-Drude peak at finite frequencies. Off-diagonal correlations are also unconventional: they feature an integrable algebraic decay for arbitrarily large values of disorder. These results do not fit the description of any known quantum phase, and strongly suggest the existence of an unusual conducting state of bosonic matter in the ground state.

Homotopic Action: A Pathway to Convergent Diagrammatic Theories.

The major obstacle preventing Feynman diagrammatic expansions from accurately solving many-fermion systems in strongly correlated regimes is the series slow convergence or divergence problem. Several techniques have been proposed to address this issue: series resummation by conformal mapping, changing the nature of the starting point of the expansion by shifted action tools, and applying the homotopy analysis method to the Dyson-Schwinger equation. They emerge as dissimilar mathematical procedures aimed at different aspects of the problem. The proposed homotopic action offers a universal and systematic framework for unifying the existing-and generating new-methods and ideas to formulate a physical system in terms of a convergent diagrammatic series. It eliminates the need for resummation, allows one to introduce effective interactions, enables a controlled ultraviolet regularization of continuous-space theories, and reduces the intrinsic polynomial complexity of the diagrammatic Monte Carlo method. We illustrate this approach by an application to the Hubbard model.

Fermi blockade of the strong electron-phonon interaction in modelled optimally doped high temperature superconductors.

We study how manifestations of strong electron-phonon interaction depend on the carrier concentration by solving the two-dimensional Holstein model for the spin-polarized fermions using an approximation free bold-line diagrammatic Monte Carlo method. We show that the strong electron-phonon interaction, obviously present at very small Fermion concentration, is masked by the Fermi blockade effects and Migdal's theorem to the extent that it manifests itself as moderate one at large carriers densities. Suppression of strong electron-phonon interaction fingerprints is in agreement with experimental observations in doped high temperature superconductors.

Polaron Mobility in the "Beyond Quasiparticles" Regime.

In a number of physical situations, from polarons to Dirac liquids and to non-Fermi liquids, one encounters the "beyond quasiparticles" regime, in which the inelastic scattering rate exceeds the thermal energy of quasiparticles. Transport in this regime cannot be described by the kinetic equation. We employ the diagrammatic Monte Carlo method to study the mobility of a Fröhlich polaron in this regime and discover a number of nonperturbative effects: a strong violation of the Mott-Ioffe-Regel criterion at intermediate and strong couplings, a mobility minimum at Tâˆ¼Ω in the strong-coupling limit (Ω is the optical mode frequency), a substantial delay in the onset of an exponential dependence of the mobility for T<Ω at intermediate coupling, and complete smearing of the Drude peak at strong coupling. These effects should be taken into account when interpreting mobility data in materials with strong electron-phonon coupling.

Supersolid Stripe Crystal from Finite-Range Interactions on a Lattice.

Strong, long-range interactions present a unique challenge for the theoretical investigation of quantum many-body lattice models, due to the generation of large numbers of competing states at low energy. Here, we investigate a class of extended bosonic Hubbard models with off-site terms interpolating between short and infinite range, thus allowing for an exact numerical solution for all interaction strengths. We predict a novel type of stripe crystal at strong coupling. Most interestingly, for intermediate interaction strengths we demonstrate that the stripes can turn superfluid, thus leading to a self-assembled array of quasi-one-dimensional superfluids. These bosonic superstripes turn into an isotropic supersolid with decreasing the interaction strength. The mechanism for stripe formation is based on cluster self-assembling in the corresponding classical ground state, reminiscent of classical soft-matter models of polymers, different from recently proposed mechanisms for cold gases of alkali or dipolar magnetic atoms.

Restoring a smooth function from its noisy integrals.

Numerical (and experimental) data analysis often requires the restoration of a smooth function from a set of sampled integrals over finite bins. We present the bin hierarchy method that efficiently computes the maximally smooth function from the sampled integrals using essentially all the information contained in the data. We perform extensive tests with different classes of functions and levels of data quality, including Monte Carlo data suffering from a severe sign problem and physical data for the Green's function of the Fröhlich polaron.

Stability of Dirac Liquids with Strong Coulomb Interaction.

We develop and apply the diagrammatic Monte Carlo technique to address the problem of the stability of the Dirac liquid state (in a graphene-type system) against the strong long-range part of the Coulomb interaction. So far, all attempts to deal with this problem in the field-theoretical framework were limited either to perturbative or random phase approximation and functional renormalization group treatments, with diametrically opposite conclusions. Our calculations aim at the approximation-free solution with controlled accuracy by computing vertex corrections from higher-order skeleton diagrams and establishing the renormalization group flow of the effective Coulomb coupling constant. We unambiguously show that with increasing the system size L (up to ln(L)∼40), the coupling constant always flows towards zero; i.e., the two-dimensional Dirac liquid is an asymptotically free T=0 state with divergent Fermi velocity.

Quantum Walk in Degenerate Spin Environments.

We study the propagation of a hole in degenerate (paramagnetic) spin environments. This canonical problem has important connections to a number of physical systems, and is perfectly suited for experimental realization with ultracold atoms in an optical lattice. At the short-to-intermediate time scale that we can access using a stochastic-series-type numeric scheme, the propagation turns out to be distinctly nondiffusive with the probability distribution featuring minima in both space and time due to quantum interference, yet the motion is not ballistic, except at the beginning. We discuss possible scenarios for long-term evolution that could be explored with an unprecedented degree of detail in experiments with single-atom resolved imaging.

Spin-Ice State of the Quantum Heisenberg Antiferromagnet on the Pyrochlore Lattice.

We study the low-temperature physics of the SU(2)-symmetric spin-1/2 Heisenberg antiferromagnet on a pyrochlore lattice and find "fingerprint" evidence for the thermal spin-ice state in this frustrated quantum magnet. Our conclusions are based on the results of bold diagrammatic Monte Carlo simulations, with good convergence of the skeleton series down to the temperature T/J=1/6. The identification of the spin-ice state is done through a remarkably accurate microscopic correspondence for the static structure factor between the quantum Heisenberg, classical Heisenberg, and Ising models at all accessible temperatures, and the characteristic bowtie pattern with pinch points observed at T/J=1/6. The dynamic structure factor at real frequencies (obtained by the analytic continuation of numerical data) is consistent with diffusive spinon dynamics at the pinch points.

p-Wave superfluidity by spin-nematic Fermi surface deformation.

We study attractively interacting fermions on a square lattice with dispersion relations exhibiting strong spin-dependent anisotropy. The resulting Fermi surface mismatch suppresses the s-wave BCS-type instability, clearing the way for unconventional types of order. Unbiased sampling of the Feynman diagrammatic series using diagrammatic Monte Carlo methods reveals a rich phase diagram in the regime of intermediate coupling strength. Instead of a proposed Cooper-pair Bose metal phase [A. E. Feiguin and M. P. A. Fisher, Phys. Rev. Lett. 103, 025303 (2009)], we find an incommensurate density wave at strong anisotropy and two different p-wave superfluid states with unconventional symmetry at intermediate anisotropy.

Diagrammatic Monte Carlo method for many-polaron problems.

We introduce the first bold diagrammatic Monte Carlo approach to deal with polaron problems at a finite electron density nonperturbatively, i.e., by including vertex corrections to high orders. Using the Holstein model on a square lattice as a prototypical example, we demonstrate that our method is capable of providing accurate results in the thermodynamic limit in all regimes from a renormalized Fermi liquid to a single polaron, across the nonadiabatic region where Fermi and Debye energies are of the same order of magnitude. By accounting for vertex corrections, the accuracy of the theoretical description is increased by orders of magnitude relative to the lowest-order self-consistent Born approximation employed in most studies. We also find that for the electron-phonon coupling typical for real materials, the quasiparticle effective mass increases and the quasiparticle residue decreases with increasing the electron density at constant electron-phonon coupling strength.

Critical exponents of the superfluid-Bose-glass transition in three dimensions.

Recent experimental and numerical studies of the critical-temperature exponent ϕ for the superfluid-Bose-glass universality in three-dimensional systems report strong violations of the key quantum critical relation, ϕ=νz, where z and ν are the dynamic and correlation-length exponents, respectively; these studies question the conventional scaling laws for this quantum critical point. Using Monte Carlo simulations of the disordered Bose-Hubbard model, we demonstrate that previous work on the superfluid-to-normal-fluid transition-temperature dependence on the chemical potential (or the magnetic field, in spin systems), T_{c}∝(µ-µ_{c})^{ϕ}, was misinterpreting transient behavior on approach to the fluctuation region with the genuine critical law. When the model parameters are modified to have a broad quantum critical region, simulations of both quantum and classical models reveal that the ϕ=νz law [with ϕ=2.7(2), z=3, and ν=0.88(5)] holds true, resolving the ϕ-exponent "crisis."

Universal conductivity in a two-dimensional superfluid-to-insulator quantum critical system.

We compute the universal conductivity of the (2+1)-dimensional XY universality class, which is realized for a superfluid-to-Mott insulator quantum phase transition at constant density. Based on large-scale Monte Carlo simulations of the classical (2+1)-dimensional J-current model and the two-dimensional Bose-Hubbard model, we can precisely determine the conductivity on the quantum critical plateau, σ(∞) = 0.359(4)σQ with σQ the conductivity quantum. The universal conductivity curve is the standard example with the lowest number of components where the bottoms-up AdS/CFT correspondence from string theory can be tested and made to use [R. C. Myers, S. Sachdev, and A. Singh, Phys. Rev. D 83, 066017 (2011)]. For the first time, the shape of the σ(iω(n)) - σ(∞) function in the Matsubara representation is accurate enough for a conclusive comparison and establishes the particlelike nature of charge transport. We find that the holographic gauge-gravity duality theory for transport properties can be made compatible with the data if temperature of the horizon of the black brane is different from the temperature of the conformal field theory. The requirements for measuring the universal conductivity in a cold gas experiment are also determined by our calculation.

Universal properties of the Higgs resonance in (2+1)-dimensional U(1) critical systems.

We present spectral functions for the magnitude squared of the order parameter in the scaling limit of the two-dimensional superfluid to Mott insulator quantum phase transition at constant density, which has emergent particle-hole symmetry and Lorentz invariance. The universal functions for the superfluid, Mott insulator, and normal liquid phases reveal a low-frequency resonance which is relatively sharp and is followed by a damped oscillation (in the first two phases only) before saturating to the quantum critical plateau. The counterintuitive resonance feature in the insulating and normal phases calls for deeper understanding of collective modes in the strongly coupled (2+1)-dimensional relativistic field theory. Our results are derived from analytically continued correlation functions obtained from path-integral Monte Carlo simulations of the Bose-Hubbard model.

Regularization of diagrammatic series with zero convergence radius.

The divergence of perturbative expansions which occurs for the vast majority of macroscopic systems and follows from Dyson's collapse argument prevents the direct use of Feynman's diagrammatic technique for controllable studies of strongly interacting systems. We show how the problem of divergence can be solved by replacing the original model with a convergent sequence of successive approximations which have a convergent perturbative series while maintaining the diagrammatic structure. As an instructive model, we consider the zero-dimensional |ψ|4 theory.