Source author record

Evgeny Kozik

Evgeny Kozik appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Misleading convergence of the skeleton diagrammatic technique: when the correct solution can be found

Convergence of the self-consistent skeleton diagrammatic technique (SDT) -- in which the full Green's function is determined through summation of Feynman diagrams in terms of itself -- to the wrong answer has been associated with the existence of non-perturbative branches of the Luttinger-Ward functional. Although it has been possible to detect misleading convergence without the knowledge of the exact result, the SDT has remained inapplicable in the regimes where this happens. We show that misleading convergence does not always preclude recovering the exact solution. In addition to the established mechanism, convergence of the SDT to the wrong answer can stem from divergence of the inherent diagrammatic series, which allows us to recover the exact solution by a modified SDT protocol based on controlled analytic continuation. We illustrate this approach by its application to the analytically solvable (0+0)d Hubbard model, the Hubbard atom, and the 2d Hubbard model in a challenging strong-coupling regime, for which the SDT is solved with controlled accuracy by the diagrammatic Monte Carlo (DiagMC) technique.

preprint2021arXiv

Evaluating second-order phase transitions with Diagrammatic Monte Carlo: Néel Transition in the doped three-dimensional Hubbard model

Diagrammatic Monte Carlo -- the technique for numerically exact summation of all Feynman diagrams to high orders -- offers a unique unbiased probe of continuous phase transitions. Being formulated directly in the thermodynamic limit, the diagrammatic series is bound to diverge and is not resummable at the transition due to the non-analyticity of physical observables. This enables the detection of the transition with controlled error bars from an analysis of the series coefficients alone, avoiding the challenge of evaluating physical observables near the transition. We demonstrate this technique by the example of the Néel transition in the $3d$ Hubbard model. At half-filling and higher temperatures, the method matches the accuracy of state-of-the-art finite-size techniques, but surpasses it at low temperatures and allows us to map the phase diagram in the doped regime, where finite-size techniques struggle from the fermion sign problem. At low temperatures and sufficient doping, the transition to an incommensurate spin density wave state is observed.

preprint2020arXiv

Extended crossover from Fermi liquid to quasi-antiferromagnet in the half-filled 2D Hubbard model

The ground state of the Hubbard model with nearest-neighbor hopping on the square lattice at half filling is known to be that of an antiferromagnetic (AFM) band insulator for any on-site repulsion. At finite temperature, the absence of long-range order makes the question of how the interaction-driven insulator is realized nontrivial. We address this problem with controlled accuracy in the thermodynamic limit using self-energy diagrammatic determinant Monte Carlo and dynamical cluster approximation methods and show that development of long-range AFM correlations drives an extended crossover from Fermi liquid to insulating behavior in the parameter regime that precludes a metal-to-insulator transition. The intermediate crossover state is best described as a non-Fermi liquid with a partially gapped Fermi surface.

preprint2020arXiv

Spin and charge correlations across the metal-to-insulator crossover in the half-filled $2d$ Hubbard model

The $2d$ Hubbard model with nearest-neighbour hopping on the square lattice and an average of one electron per site is known to undergo an extended crossover from metallic to insulating behavior driven by proliferating antiferromagnetic correlations. We study signatures of this crossover in spin and charge correlation functions and present results obtained with controlled accuracy using diagrammatic Monte Carlo in the range of parameters amenable to experimental verification with ultracold atoms in optical lattices. The qualitative changes in charge and spin correlations associated with the crossover are observed at well-separated temperature scales, which encase the intermediary regime of non-Fermi-liquid character, where local magnetic moments are formed and non-local fluctuations in both channels are essential.

preprint2016arXiv

Fulde-Ferrell-Larkin-Ovchinnikov pairing as leading instability on the square lattice

We study attractively interacting spin-1/2 fermions on the square lattice subject to a spin population imbalance. Using unbiased diagrammatic Monte Carlo simulations we find an extended region in the parameter space where the Fermi liquid is unstable towards formation of Cooper pairs with non-zero center-of-mass momentum, known as the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state. In contrast to earlier mean-field and quasi-classical studies we provide quantitative and well-controlled predictions on the existence and location of the relevant Fermi-liquid instabilities. The highest temperature where the FFLO instability can be observed is about half of the superfluid transition temperature in the unpolarized system.

preprint2015arXiv

Emergent BCS regime of the two-dimensional fermionic Hubbard model: ground-state phase diagram

A significant part of the phase diagram of the two-dimensional fermionic Hubbard model for moderate interactions and filling factors ($U < 4, \, n<0.7$) is governed by effective Fermi liquid physics with weak BCS-type instabilities. We access this regime in a controlled way by a combination of the bold-line diagrammatic Monte Carlo method with an additional ladder-diagram summation trick and semi-analytic treatment of the weak instability in the Cooper channel. We obtain the corresponding ground-state phase diagram in the $(n,U)$ plane describing the competition between the $p-$ and $d-$wave superfluid states. We also claim the values of the dimensionless BCS coupling constants controlling the superfluid $T_c$ at the phase boundaries, which prove to be very small up to $U=4, n = 0.6$.

preprint2015arXiv

Ground state phase diagram of the repulsive fermionic $t-t^{\prime}$ Hubbard model on the square lattice from weak-coupling

We obtain a complete and exact in the weak-coupling limit ($U \rightarrow 0$) ground state phase diagram of the repulsive fermionic Hubbard model on the square lattice for filling factors $0 < n < 2$ and next-nearest-neighbour hopping amplitudes $0 \le t^{\prime} \le 0.5$. Phases are distinguished by the symmetry and the number of nodes of the superfluid order parameter. The phase diagram is richer than may be expected and typically features states with a high --- higher than that of the fundamental mode of the corresponding irreducible representation --- number of nodes. The effective coupling strength in the Cooper channel $λ$, which determines the critical temperature $T_c$ of the superfluid transition, is calculated in the whole parameter space and regions with high values of $λ$ are identified. It is shown that besides the expected increase of $λ$ near the Van Hove singularity line, joining the ferromagnetic and antiferromagnetic points, another region with high values of $λ$ can be found at quarter filling and $t^{\prime}=0.5$ due to the presence of a line of nesting at $t^{\prime} \ge 0.5$. The results can serve as benchmarks for controlled non-perturbative methods and guide the ongoing search for high-$T_c$ superconductivity in the Hubbard model.

preprint2015arXiv

Non-existence of the Luttinger-Ward functional and misleading convergence of skeleton diagrammatic series for Hubbard-like models

The Luttinger-Ward functional $Φ[\mathbf{G}]$, which expresses the thermodynamic grand potential in terms of the interacting single-particle Green's function $\mathbf{G}$, is found to be ill-defined for fermionic models with the Hubbard on-site interaction. In particular, we show that the self-energy $\mathbfΣ[\mathbf{G}] \propto δΦ[\mathbf{G}]/δ\mathbf{G}$ is not a single-valued functional of $\mathbf{G}$: in addition to the physical solution for $\mathbfΣ[\mathbf{G}]$, there exists at least one qualitatively distinct unphysical branch. This result is demonstrated for several models: the Hubbard atom, the Anderson impurity model, and the full two-dimensional Hubbard model. Despite this pathology, the skeleton Feynman diagrammatic series for $\mathbfΣ$ in terms of $\mathbf{G}$ is found to converge at least for moderately low temperatures. However, at strong interactions, its convergence is to the unphysical branch. This reveals a new scenario of breaking down of diagrammatic expansions. In contrast, the bare series in terms of the non-interacting Green's function $\mathbf{G}_0$ converges to the correct physical branch of $\mathbfΣ$ in all cases currently accessible by diagrammatic Monte Carlo. Besides their conceptual importance, these observations have important implications for techniques based on the explicit summation of diagrammatic series.

preprint2015arXiv

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.

preprint2010arXiv

Comment on "Symmetries and Interaction Coefficients of Kelvin waves" [arXiv:1005.4575] by Lebedev and L'vov

We comment on the claim by Lebedev and L'vov [arXiv:1005.4575] that the symmetry with respect to a tilt of a quantized vortex line does not yet prohibit coupling between Kelvin waves and the large-scale slope of the line. Ironically, the counterexample of an effective scattering vertex in the local induction approximation (LIA) attempted by Lebedev and L'vov invalidates their logic all by itself being a notoriously known example of how symmetries impose stringent constraints on kelvon kinetics---not only the coupling in question but the kinetics in general are absent within LIA. We further explain that the mistake arises from confusing symmetry properties of a specific mathematical representation in terms of the canonical vortex position field w(z) = x(z) + iy(z), which explicitly breaks the tilt symmetry due to an arbitrary choice of the z-axis, with those of the real physical system recovered in final expressions.

preprint2010arXiv

Geometric Symmetries in Superfluid Vortex Dynamics

Dynamics of quantized vortex lines in a superfluid feature symmetries associated with the geometric character of the complex-valued field, $w(z)=x(z)+iy(z)$, describing the instant shape of the line. Along with a natural set of Noether's constants of motion, which---apart from their rather specific expressions in terms of $w(z)$---are nothing but components of the total linear and angular momenta of the fluid, the geometric symmetry brings about crucial consequences for kinetics of distortion waves on the vortex lines---the Kelvin waves. It is the geometric symmetry that renders Kelvin-wave cascade local in the wavenumber space. Similar considerations apply to other systems with purely geometric degrees of freedom.

preprint2010arXiv

Quantitative Determination of Temperature in the Approach to Magnetic Order of Ultracold Fermions in an Optical Lattice

We perform a quantitative simulation of the repulsive Fermi-Hubbard model using an ultracold gas trapped in an optical lattice. The entropy of the system is determined by comparing accurate measurements of the equilibrium double occupancy with theoretical calculations over a wide range of parameters. We demonstrate the applicability of both high-temperature series and dynamical mean-field theory to obtain quantitative agreement with the experimental data. The reliability of the entropy determination is confirmed by a comprehensive analysis of all systematic errors. In the center of the Mott insulating cloud we obtain an entropy per atom as low as 0.77kB which is about twice as large as the entropy at the Neel transition. The corresponding temperature depends on the atom number and for small fillings reaches values on the order of the tunneling energy.

preprint2010arXiv

Vortex-Phonon Interaction in the Kosterlitz-Thouless Theory

The "canonical" variables of the Kosterlitz-Thouless theory--fields $Φ_0({\bf r})$ and $ϕ({\bf r})$, generally believed to stand for vortices and phonons (or their XY equivalents, like spin waves, etc.) turn out to be neither vortices and phonons, nor, strictly speaking, {\it canonical} variables. The latter fact explains paradoxes of (i) absence of interaction between $Φ_0$ and $ϕ$, and (ii) non-physical contribution of small vortex pairs to long-range phase correlations. We resolve the paradoxes by explicitly relating $Φ_0$ and $ϕ$ to canonical vortex-pair and phonon variables.

preprint2007arXiv

Kolmogorov and Kelvin-Wave Cascades of Superfluid Turbulence at T=0: What is in Between?

As long as vorticity quantization remains irrelevant for the long-wave physics, superfluid turbulence supports a regime macroscopically identical to the Kolmogorov cascade of a normal liquid. At high enough wavenumbers, the energy flux in the wavelength space is carried by individual Kelvin-wave cascades on separate vortex lines. We analyze the transformation of the Kolmogorov cascade into the Kelvin-wave cascade, revealing a chain of three distinct intermediate cascades, supported by local-induction motion of the vortex lines, and distinguished by specific reconnection mechanisms. The most prominent qualitative feature predicted is unavoidable production of vortex rings of the size of the order of inter-vortex distance.

preprint2005arXiv

Vortex--Phonon Interaction

Kelvin waves (kelvons)--helical waves on quantized vortex lines--are the normal modes of vortices in a superfluid. At zero temperature, the only dissipative channel of vortex dynamics is phonon emission. Starting with the hydrodynamic action, we derive the Hamiltonian of vortex-phonon interaction, thereby reducing the problem of the interaction of Kelvin waves with sound to inelastic elementary excitation scattering. On the basis of this formalism, we calculate the rate of sound radiation by superfluid turbulence at zero temperature and estimate the value of short-wavelength cutoff of the turbulence spectrum.

preprint2004arXiv

Scale Separation Scheme for Simulating Superfluid Turbulence: Kelvin-Wave Cascade

A Kolmogorov-type cascade of Kelvin waves--the distortion waves on vortex lines--plays a key part in the relaxation of superfluid turbulence at low temperatures. We propose an efficient numeric scheme for simulating the Kelvin wave cascade on a single vortex line. The idea is likely to be generalizable for a full-scale simulation of different regimes of superfluid turbulence. With the new scheme, we are able to unambiguously resolve the cascade spectrum exponent, and thus to settle the controversy between recent simulations [1] and recently developed analytic theory [2]. [1] W.F. Vinen, M. Tsubota and A. Mitani, Phys. Rev. Lett. 91, 135301 (2003). [2] E.V. Kozik and B.V. Svistunov, Phys. Rev. Lett. 92, 035301 (2004).