Showing 1–46 of 46 results for author: Kozik, E

Strange Metal to Insulator Transitions in the Lowest Landau Level

Authors: Ben Currie, Evgeny Kozik

Abstract: We study the microscopic model of electrons in the partially-filled lowest Landau level interacting via the Coulomb potential by the diagrammatic theory within the GW approximation. In a wide range of filling fractions and temperatures, we find a homogeneous non-Fermi liquid (nFL) state similar to that found in the Sachdev-Ye-Kitaev (SYK) model, with logarithmic corrections to the anomalous dimens… ▽ More We study the microscopic model of electrons in the partially-filled lowest Landau level interacting via the Coulomb potential by the diagrammatic theory within the GW approximation. In a wide range of filling fractions and temperatures, we find a homogeneous non-Fermi liquid (nFL) state similar to that found in the Sachdev-Ye-Kitaev (SYK) model, with logarithmic corrections to the anomalous dimension. In addition, the phase diagram is qualitatively similiar to that of SYK: a first-order transition terminating at a critical end-point separates the nFL phase from a band insulator that corresponds to the fully-filled Landau level. This critical point, as well as that of the SYK model -- whose critical exponents we determine more precisely -- are shown to both belong to the Van der Waals universality class. The possibility of a charge density wave (CDW) instability is also investigated, and we find the homogeneous nFL state to extend down to the ground state for fillings $0.2 \lesssim νにゅー\lesssim 0.8$, while a CDW appears outside this range of fillings at sufficiently low temperatures. Our results suggest that the SYK-like nFL state should be a generic feature of the partially-filled lowest Landau level at intermediate temperatures. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Comments: 11+6 pages, 6 figures

Combinatorial summation of Feynman diagrams: Equation of state of the 2D SU(N) Hubbard model

Authors: Evgeny Kozik

Abstract: Feynman's diagrammatic series is a common language for a formally exact theoretical description of systems of infinitely-many interacting quantum particles, as well as a foundation for precision computational techniques. Here we introduce a universal framework for efficient summation of connected or skeleton Feynman diagrams for generic quantum many-body systems. It is based on an explicit combina… ▽ More Feynman's diagrammatic series is a common language for a formally exact theoretical description of systems of infinitely-many interacting quantum particles, as well as a foundation for precision computational techniques. Here we introduce a universal framework for efficient summation of connected or skeleton Feynman diagrams for generic quantum many-body systems. It is based on an explicit combinatorial construction of the sum of the integrands by dynamic programming, at a computational cost that can be made only exponential in the diagram order on a classical computer and potentially polynomial on a quantum computer. We illustrate the technique by an unbiased diagrammatic Monte Carlo calculation of the equation of state of the $2D$ $SU(N)$ Hubbard model in an experimentally relevant regime, which has remained challenging for state-of-the-art numerical methods. △ Less

Submitted 9 April, 2024; v1 submitted 24 September, 2023; originally announced September 2023.

Comments: 10 pages (a section on nonsymmetrised summation is added), 5 figures

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

Authors: Aaram J. Kim, Evgeny Kozik

Abstract: 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 e… ▽ More 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. △ Less

Submitted 30 December, 2022; originally announced December 2022.

Comments: 6 pages, 5 figures

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

Authors: Connor Lenihan, Aaram J. Kim, Fedor Šimkovic IV., Evgeny Kozik

Abstract: 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… ▽ More 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. △ Less

Submitted 30 December, 2021; originally announced December 2021.

Physical and unphysical regimes of self-consistent many-body perturbation theory

Authors: K. Van Houcke, E. Kozik, R. Rossi, Y. Deng, F. Werner

Abstract: In the standard framework of self-consistent many-body perturbation theory, the skeleton series for the self-energy is truncated at a finite order $\mathcal{N}$ and plugged into the Dyson equation, which is then solved for the propagator $G_{\mathcal{N}}\,$. We consider two examples of fermionic models, the Hubbard atom at half filling and its zero space-time dimensional simplified version. First,… ▽ More In the standard framework of self-consistent many-body perturbation theory, the skeleton series for the self-energy is truncated at a finite order $\mathcal{N}$ and plugged into the Dyson equation, which is then solved for the propagator $G_{\mathcal{N}}\,$. We consider two examples of fermionic models, the Hubbard atom at half filling and its zero space-time dimensional simplified version. First, we show that $G_{\mathcal{N}}\,$ converges when $\mathcal{N}\to\infty$ to a limit $G_\infty\,$, which coincides with the exact physical propagator $G_{\rm exact}\ $ at small enough coupling, while $G_\infty \neq G_{\rm exact}\ $ at strong coupling. This follows from the findings of [Kozik, Ferrero and Georges, PRL 114, 156402 (2015)] and an additional subtle mathematical mechanism elucidated here. Second, we demonstrate that it is possible to discriminate between the $G_\infty=G_{\rm exact}\ $ and $G_\infty\neq G_{\rm exact}\ $ regimes thanks to a criterion which does not require the knowledge of $G_{\rm exact}\ $, as proposed in [Rossi et al., PRB 93, 161102(R) (2016)]. △ Less

Submitted 27 May, 2024; v1 submitted 8 February, 2021; originally announced February 2021.

Journal ref: SciPost Phys. 16, 133 (2024)

Strange Metal Solution in the Diagrammatic Theory for the $2d$ Hubbard Model

Authors: Aaram J. Kim, Philipp Werner, Evgeny Kozik

Abstract: We show that the numerically exact bold-line diagrammatic theory for the $2d$ Hubbard model exhibits a non-Fermi-liquid (NFL) strange metal state, which is connected to the SYK NFL in the strong-interaction limit. The solution for the doped system features the expected phenomenology with the NFL near half-filling at strong couplings and in a wide temperature range enclosed by the atomic state at h… ▽ More We show that the numerically exact bold-line diagrammatic theory for the $2d$ Hubbard model exhibits a non-Fermi-liquid (NFL) strange metal state, which is connected to the SYK NFL in the strong-interaction limit. The solution for the doped system features the expected phenomenology with the NFL near half-filling at strong couplings and in a wide temperature range enclosed by the atomic state at high temperatures and a Fermi liquid at low temperatures. We demonstrate, however, that this behavior in the weakly doped regime is due to the unphysical branch of the Luttinger-Ward functional. On the other hand, our analysis shows that the NFL physics is realized at larger doping. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: 7 pages, 5 figures

Homotopic Action: A Pathway to Convergent Diagrammatic Theories

Authors: Aaram J. Kim, Nikolay V. Prokof'ev, Boris V. Svistunov, Evgeny Kozik

Abstract: 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 homot… ▽ More 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. △ Less

Submitted 23 June, 2021; v1 submitted 11 October, 2020; originally announced October 2020.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. Lett. 126, 257001 (2021)

Dual Boson Diagrammatic Monte Carlo Approach Applied to the Extended Hubbard Model

Authors: M. Vandelli, V. Harkov, E. A. Stepanov, J. Gukelberger, E. Kozik, A. Rubio, A. I. Lichtenstein

Abstract: In this work we introduce the Dual Boson Diagrammatic Monte Carlo technique for strongly interacting electronic systems. This method combines the strength of dynamical mean-filed theory for non-perturbative description of local correlations with the systematic account of non-local corrections in the Dual Boson theory by the diagrammatic Monte Carlo approach. It allows us to get a numerically exact… ▽ More In this work we introduce the Dual Boson Diagrammatic Monte Carlo technique for strongly interacting electronic systems. This method combines the strength of dynamical mean-filed theory for non-perturbative description of local correlations with the systematic account of non-local corrections in the Dual Boson theory by the diagrammatic Monte Carlo approach. It allows us to get a numerically exact solution of the dual boson theory at the two-particle local vertex level for the extended Hubbard model. We show that it can be efficiently applied to description of single particle observables in a wide range of interaction strengths. We compare our exact results for the self-energy with the ladder Dual Boson approach and determine a physical regime, where description of collective electronic effects requires more accurate consideration beyond the ladder approximation. Additionally, we find that the order-by-order analysis of the perturbative diagrammatic series for the single-particle Green's function allows to estimate the transition point to the charge density wave phase. △ Less

Submitted 19 October, 2020; v1 submitted 9 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. B 102, 195109 (2020)

Tracking the Footprints of Spin Fluctuations: A MultiMethod, MultiMessenger Study of the Two-Dimensional Hubbard Model

Authors: Thomas Schäfer, Nils Wentzell, Fedor Šimkovic IV, Yuan-Yao He, Cornelia Hille, Marcel Klett, Christian J. Eckhardt, Behnam Arzhang, Viktor Harkov, François-Marie Le Régent, Alfred Kirsch, Yan Wang, Aaram J. Kim, Evgeny Kozik, Evgeny A. Stepanov, Anna Kauch, Sabine Andergassen, Philipp Hansmann, Daniel Rohe, Yuri M. Vilk, James P. F. LeBlanc, Shiwei Zhang, A. -M. S. Tremblay, Michel Ferrero, Olivier Parcollet , et al. (1 additional authors not shown)

Abstract: The Hubbard model represents the fundamental model for interacting quantum systems and electronic correlations. Using the two-dimensional half-filled Hubbard model at weak coupling as a testing ground, we perform a comparative study of a comprehensive set of state of the art quantum many-body methods. Upon cooling into its insulating antiferromagnetic ground-state, the model hosts a rich sequence… ▽ More The Hubbard model represents the fundamental model for interacting quantum systems and electronic correlations. Using the two-dimensional half-filled Hubbard model at weak coupling as a testing ground, we perform a comparative study of a comprehensive set of state of the art quantum many-body methods. Upon cooling into its insulating antiferromagnetic ground-state, the model hosts a rich sequence of distinct physical regimes with crossovers between a high-temperature incoherent regime, an intermediate temperature metallic regime and a low-temperature insulating regime with a pseudogap created by antiferromagnetic fluctuations. We assess the ability of each method to properly address these physical regimes and crossovers through the computation of several observables probing both quasiparticle properties and magnetic correlations, with two numerically exact methods (diagrammatic and determinantal quantum Monte Carlo) serving as a benchmark. By combining computational results and analytical insights, we elucidate the nature and role of spin fluctuations in each of these regimes. Based on this analysis, we explain how quasiparticles can coexist with increasingly long-range antiferromagnetic correlations, and why dynamical mean-field theory is found to provide a remarkably accurate approximation of local quantities in the metallic regime. We also critically discuss whether imaginary time methods are able to capture the non-Fermi liquid singularities of this fully nested system. △ Less

Submitted 29 March, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 55 pages, 39 figures, 288 references

Journal ref: Phys. Rev. X 11, 011058 (2021)

GW space-time method: Energy band-gap of solid hydrogen

Authors: Sam Azadi, Arkadiy Davydov, Evgeny Kozik

Abstract: We implement the GW space-time method at finite temperatures, in which the Green's function G and the screened Coulomb interaction W are represented in the real space on a suitable mesh and in imaginary time in terms of Chebyshev polynomials, paying particular attention to controlling systematic errors of the representation. Having validated the technique by the canonical application to silicon an… ▽ More We implement the GW space-time method at finite temperatures, in which the Green's function G and the screened Coulomb interaction W are represented in the real space on a suitable mesh and in imaginary time in terms of Chebyshev polynomials, paying particular attention to controlling systematic errors of the representation. Having validated the technique by the canonical application to silicon and germanium, we apply it to calculation of band gaps in hexagonal solid hydrogen with the bare Green's function obtained from density functional approximation and the interaction screened within the random phase approximation (RPA). The band gap results, obtained from the asymptotic decay of the full Green's function without resorting to analytic continuation, suggest that the solid hydrogen above 270 GPa can not adopt the hexagonal-closed-pack (hcp) structure. The demonstrated ability of the method to store the full G and W functions in memory with sufficient accuracy is crucial for its subsequent extensions to include higher orders of the diagrammatic series by means of diagrammatic Monte Carlo algorithms. △ Less

Submitted 6 December, 2021; v1 submitted 11 June, 2020; originally announced June 2020.

Entropy in the non-Fermi-liquid regime of the doped $2d$ Hubbard model

Authors: Connor Lenihan, Aaram J. Kim, Fedor Šimkovic IV., Evgeny Kozik

Abstract: We study thermodynamic properties of the doped Hubbard model on the square lattice in the regime of strong charge and spin fluctuations at low temperatures near the metal-to-insulator crossover and obtain results with controlled accuracy using the diagrammatic Monte Carlo method directly in the thermodynamic limit. The behavior of the entropy reveals a non-Fermi-liquid state at sufficiently high i… ▽ More We study thermodynamic properties of the doped Hubbard model on the square lattice in the regime of strong charge and spin fluctuations at low temperatures near the metal-to-insulator crossover and obtain results with controlled accuracy using the diagrammatic Monte Carlo method directly in the thermodynamic limit. The behavior of the entropy reveals a non-Fermi-liquid state at sufficiently high interactions near half-filling: a maximum in the entropy at non-zero doping develops as the coupling strength is increased, along with an inflection point, evidencing a metal to non-Fermi-liquid crossover. The specific heat exhibits additional distinctive features of a non-Fermi-liquid state. Measurements of the entropy can therefore be used as a probe of the state of the system in quantum simulation experiments with ultracold atoms in optical lattices. △ Less

Submitted 23 December, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 126, 105701 (2021)

Superfluid ground state phase diagram of the two-dimensional Hubbard model in the emergent Bardeen-Cooper-Schrieffer regime

Authors: Fedor Šimkovic IV, Youjin Deng, Evgeny Kozik

Abstract: In nonperturbative regimes, the superfluid instability in the two-dimensional Hubbard model can be described by an emergent BCS theory with small effective pairing constants. We compute the effective couplings using a controlled bold-line diagrammatic Monte Carlo approach, which stochastically sums all skeleton Feynman diagrams dressed in the one- and two-particle channels to high expansion orders… ▽ More In nonperturbative regimes, the superfluid instability in the two-dimensional Hubbard model can be described by an emergent BCS theory with small effective pairing constants. We compute the effective couplings using a controlled bold-line diagrammatic Monte Carlo approach, which stochastically sums all skeleton Feynman diagrams dressed in the one- and two-particle channels to high expansion orders, and map out the resulting superfluid ground-state phase diagram in a range of next-nearest-neighbor hopping $0 \leq t^{\prime} \leq 0.3t$, interaction strength $0 \leq U \leq 3t$, and lattice filling $0 \leq n \leq 2$. The phase diagram is dramatically transformed in the hole-doped region and becomes particularly rich at larger doping and $t'$. At $t'=0.3$, the weak-coupling picture with the dominant triplet pairing sharply peaked at $n \approx 0.73$ due to the Van Hove singularity is replaced by a plateau of the singlet $d_{x^2-y^2}$ paring, while for $U \gtrsim 3t$ the effective couplings are consistent with the $d_{x^2-y^2}$ high-temperature superconductivity in the hole-doped region near cuprates' optimal doping. △ Less

Submitted 20 July, 2021; v1 submitted 30 December, 2019; originally announced December 2019.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 104, 020507 (2021)

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

Authors: Aaram J. Kim, Fedor Simkovic IV, Evgeny Kozik

Abstract: 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 diagrammat… ▽ More 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. △ Less

Submitted 18 March, 2020; v1 submitted 30 May, 2019; originally announced May 2019.

Comments: 6 pages, 5 figures, published version

Journal ref: Phys. Rev. Lett. 124, 117602 (2020)

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

Authors: Fedor Šimkovic IV, J. P. F. LeBlanc, Aaram J. Kim, Youjin Deng, N. V. Prokof'ev, B. V. Svistunov, Evgeny Kozik

Abstract: 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 thermody… ▽ More 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. △ Less

Submitted 1 February, 2020; v1 submitted 30 December, 2018; originally announced December 2018.

Comments: 6 pages, 4 figures, with supplemental material: 2 pages, 3 figures

Journal ref: Phys. Rev. Lett. 124, 017003 (2020)

Determinant Monte Carlo for irreducible Feynman diagrams in the strongly correlated regime

Authors: Fedor Simkovic IV., Evgeny Kozik

Abstract: We develop a numerically exact method for the summation of irreducible Feynman diagrams for fermionic self-energy in the thermodynamic limit. The technique, based on the Diagrammatic Determinant Monte Carlo and its recent extension to connected diagrams, allows us to reach high ($\sim 10$) orders of the weak-coupling expansion for the self-energy of the two-dimensional Hubbard model. Access to hig… ▽ More We develop a numerically exact method for the summation of irreducible Feynman diagrams for fermionic self-energy in the thermodynamic limit. The technique, based on the Diagrammatic Determinant Monte Carlo and its recent extension to connected diagrams, allows us to reach high ($\sim 10$) orders of the weak-coupling expansion for the self-energy of the two-dimensional Hubbard model. Access to high orders reveals a non-trivial analytic structure of the self-energy and enables its controlled reconstruction with arbitrary momentum resolution in the nonperturbative regime of essentially strong correlations, which has recently been reached with ultracold atoms in optical lattices. △ Less

Submitted 4 September, 2019; v1 submitted 28 December, 2017; originally announced December 2017.

Comments: 5 pages, 4 figures, published version

Journal ref: Phys. Rev. B 100, 121102 (2019)

Magnetic Correlations in the Two-dimensional Repulsive Fermi Hubbard Model

Authors: Fedor Šimkovic IV., Youjin Deng, N. V. Prokof'ev, B. V. Svistunov, I. Tupitsyn, Evgeny Kozik

Abstract: The repulsive Fermi Hubbard model on the square lattice has a rich phase diagram near half-filling (corresponding to the particle density per lattice site $n=1$): for $n=1$ the ground state is an antiferromagnetic insulator, at $0.6 < n \lesssim 0.8$, it is a $d_{x^2-y^2}$-wave superfluid (at least for moderately strong interactions $U \lesssim 4t$ in terms of the hopping $t$), and the region… ▽ More The repulsive Fermi Hubbard model on the square lattice has a rich phase diagram near half-filling (corresponding to the particle density per lattice site $n=1$): for $n=1$ the ground state is an antiferromagnetic insulator, at $0.6 < n \lesssim 0.8$, it is a $d_{x^2-y^2}$-wave superfluid (at least for moderately strong interactions $U \lesssim 4t$ in terms of the hopping $t$), and the region $1-n \ll 1$ is most likely subject to phase separation. Much of this physics is preempted at finite temperatures and to an extent driven by strong magnetic fluctuations, their quantitative characteristics and how they change with the doping level being much less understood. Experiments on ultra-cold atoms have recently gained access to this interesting fluctuation regime, which is now under extensive investigation. In this work we employ a self-consistent skeleton diagrammatic approach to quantify the characteristic temperature scale $T_{M}(n)$ for the onset of magnetic fluctuations with a large correlation length and identify their nature. Our results suggest that the strongest fluctuations---and hence highest $T_{M}$ and easiest experimental access to this regime---are observed at $U/t \approx 4-6$. △ Less

Submitted 22 June, 2017; originally announced June 2017.

Comments: 5 pages, 6 figures

Journal ref: Phys. Rev. B 96, 081117 (2017)

Diagrammatic Monte Carlo approach for diagrammatic extensions of dynamical mean-field theory -- convergence analysis of the dual fermion technique

Authors: Jan Gukelberger, Evgeny Kozik, Hartmut Hafermann

Abstract: The dual-fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in th… ▽ More The dual-fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in the two-particle scattering vertex to a leading-order or ladder-type approximation. In this work we compute the dual-fermion expansion for the two-dimensional Hubbard model including all diagram topologies with two-particle interactions to high orders by means of a stochastic diagrammatic Monte Carlo algorithm. We benchmark the obtained self-energy against numerically exact Diagrammatic Determinant Monte Carlo simulations to systematically assess convergence of the dual-fermion series and the validity of these approximations. We observe that, from high temperatures down to the vicinity of the DMFT Néel transition, the dual-fermion series converges very quickly to the exact solution in the whole range of Hubbard interactions considered ($4 \leq U/t \leq 12$), implying that contributions from higher-order vertices are small. As the temperature is lowered further, we observe slower series convergence, convergence to incorrect solutions, and ultimately divergence. This happens in a regime where magnetic correlations become significant. We find however that the self-consistent particle-hole ladder approximation yields reasonable and often even highly accurate results in this regime. △ Less

Submitted 8 May, 2017; v1 submitted 22 November, 2016; originally announced November 2016.

Comments: 25 pages, 20 figures; v2: improved presentation and figures, new overview section

Journal ref: Phys. Rev. B 96, 035152 (2017)

Controlling Feynman diagrammatic expansions: physical nature of the pseudo gap in the two-dimensional Hubbard model

Authors: Wei Wu, Michel Ferrero, Antoine Georges, Evgeny Kozik

Abstract: We introduce a method for summing Feynman's perturbation series based on diagrammatic Monte Carlo that significantly improves its convergence properties. This allows us to investigate in a controllable manner the pseudogap regime of the Hubbard model and to study the nodal/antinodal dichotomy at low doping and intermediate coupling. Marked differences from the weak coupling scenario are manifest,… ▽ More We introduce a method for summing Feynman's perturbation series based on diagrammatic Monte Carlo that significantly improves its convergence properties. This allows us to investigate in a controllable manner the pseudogap regime of the Hubbard model and to study the nodal/antinodal dichotomy at low doping and intermediate coupling. Marked differences from the weak coupling scenario are manifest, such as a higher degree of incoherence at the antinodes than at the `hot spots'. Our results show that the pseudogap and reduction of quasiparticle coherence at the antinode is due to antiferromagnetic spin correlations centered around the commensurate $(πぱい,πぱい)$ wavevector. In contrast, the dominant source of scattering at the node is associated with incommensurate momentum transfer. Umklapp scattering is found to play a key role in the nodal/antinodal dichotomy. △ Less

Submitted 30 August, 2016; originally announced August 2016.

Comments: 5+ pages, 4 figures

Journal ref: Phys. Rev. B 96, 041105 (2017)

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

Authors: Fedor Simkovic, Xuan-Wen Liu, Youjin Deng, Evgeny Kozik

Abstract: 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 riche… ▽ More 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. △ Less

Submitted 14 December, 2015; originally announced December 2015.

Comments: 11 Pages, 9 Figures

Journal ref: Phys. Rev. B 94, 085106 (2016)

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

Authors: Jan Gukelberger, Sebastian Lienert, Evgeny Kozik, Lode Pollet, Matthias Troyer

Abstract: 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 e… ▽ More 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. △ Less

Submitted 19 June, 2016; v1 submitted 16 September, 2015; originally announced September 2015.

Comments: 7 pages, 4 figures; v2: improved references and discussion, added calculations with larger cutoff order that corroborate our earlier results

Journal ref: Phys. Rev. B 94, 075157 (2016)

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Authors: J. P. F. LeBlanc, Andrey E. Antipov, Federico Becca, Ireneusz W. Bulik, Garnet Kin-Lic Chan, Chia-Min Chung, Youjin Deng, Michel Ferrero, Thomas M. Henderson, Carlos A. Jiménez-Hoyos, E. Kozik, Xuan-Wen Liu, Andrew J. Millis, N. V. Prokof'ev, Mingpu Qin, Gustavo E. Scuseria, Hao Shi, B. V. Svistunov, Luca F. Tocchio, I. S. Tupitsyn, Steven R. White, Shiwei Zhang, Bo-Xiao Zheng, Zhenyue Zhu, Emanuel Gull

Abstract: Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare… ▽ More Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods. △ Less

Submitted 15 December, 2015; v1 submitted 9 May, 2015; originally announced May 2015.

Journal ref: Phys. Rev. X 5, 041041 (2015)

p-wave Superfluidity by Spin-Nematic Fermi Surface Deformation

Authors: Jan Gukelberger, Evgeny Kozik, Lode Pollet, Nikolay Prokof'ev, Manfred Sigrist, Boris Svistunov, Matthias Troyer

Abstract: 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 regim… ▽ More 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. △ Less

Submitted 21 January, 2015; v1 submitted 13 August, 2014; originally announced August 2014.

Comments: published version, supplemental material included as appendix, 4+4 pages, 5+7 figures

Journal ref: Phys. Rev. Lett. 113, 195301 (2014)

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

Authors: Youjin Deng, Evgeny Kozik, Nikolay V. Prokof'ev, Boris V. Svistunov

Abstract: 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-… ▽ More 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$. △ Less

Submitted 5 July, 2015; v1 submitted 9 August, 2014; originally announced August 2014.

Comments: 5 pages, 5 figures

Journal ref: EPL 110 (2015) 57001

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

Authors: Evgeny Kozik, Michel Ferrero, Antoine Georges

Abstract: 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… ▽ More 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. △ Less

Submitted 17 April, 2015; v1 submitted 21 July, 2014; originally announced July 2014.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Lett. 114, 156402 (2015)

Contact and Momentum Distribution of the Unitary Fermi Gas

Authors: R. Rossi, T. Ohgoe, E. Kozik, N. Prokof'ev, B. Svistunov, K. Van Houcke, F. Werner

Abstract: A key quantity in strongly-interacting resonant Fermi gases is the contact $\mathcal{C}$, which characterizes numerous properties such as the momentum distribution at large momenta or the pair correlation function at short distances. The temperature dependence of $\mathcal{C}$ was measured at unitarity, where existing theoretical predictions differ substantially even at the qualitative level. We r… ▽ More A key quantity in strongly-interacting resonant Fermi gases is the contact $\mathcal{C}$, which characterizes numerous properties such as the momentum distribution at large momenta or the pair correlation function at short distances. The temperature dependence of $\mathcal{C}$ was measured at unitarity, where existing theoretical predictions differ substantially even at the qualitative level. We report accurate data for the contact and the momentum distribution of the unitary gas in the normal phase, obtained by Bold Diagrammatic Monte Carlo and Borel resummation. Our results agree with experimental data within error bars and provide crucial benchmarks for the development of advanced theoretical treatments and precision measurements. △ Less

Submitted 10 September, 2018; v1 submitted 25 March, 2013; originally announced March 2013.

Comments: Updated data using conformal-Borel resummation. Added observation of non-Fermi liquid behavior (Fig. 6). Supplemental Material in ancillary files

Journal ref: Phys. Rev. Lett. 121, 130406 (2018)

Néel temperature and thermodynamics of the half-filled 3D Hubbard model by Diagrammatic Determinant Monte Carlo

Authors: E. Kozik, E. Burovski, V. W. Scarola, M. Troyer

Abstract: We study thermodynamics of the 3D Hubbard model at half filling on approach to the Néel transition by means of large-scale unbiased Diagrammatic Determinant Monte Carlo simulations. We obtain the transition temperature in the strongly correlated regime, as well as temperature dependence of energy, entropy, double occupancy, and the nearest-neighbor spin correlation function. Our results improve th… ▽ More We study thermodynamics of the 3D Hubbard model at half filling on approach to the Néel transition by means of large-scale unbiased Diagrammatic Determinant Monte Carlo simulations. We obtain the transition temperature in the strongly correlated regime, as well as temperature dependence of energy, entropy, double occupancy, and the nearest-neighbor spin correlation function. Our results improve the accuracy of previous unbiased studies and present accurate benchmarks in the ongoing effort to realize the antiferromagnetic state of matter with ultracold atoms in optical lattices. △ Less

Submitted 24 May, 2013; v1 submitted 12 December, 2012; originally announced December 2012.

Comments: 17 pages, 12 figures, 5 tables

Journal ref: Phys. Rev. B 87, 205102 (2013)

Feynman diagrams versus Fermi-gas Feynman emulator

Authors: K. Van Houcke, F. Werner, E. Kozik, N. Prokofev, B. Svistunov, M. J. H. Ku, A. T. Sommer, L. W. Cheuk, A. Schirotzek, M. W. Zwierlein

Abstract: Precise understanding of strongly interacting fermions, from electrons in modern materials to nuclear matter, presents a major goal in modern physics. However, the theoretical description of interacting Fermi systems is usually plagued by the intricate quantum statistics at play. Here we present a cross-validation between a new theoretical approach, Bold Diagrammatic Monte Carlo (BDMC), and precis… ▽ More Precise understanding of strongly interacting fermions, from electrons in modern materials to nuclear matter, presents a major goal in modern physics. However, the theoretical description of interacting Fermi systems is usually plagued by the intricate quantum statistics at play. Here we present a cross-validation between a new theoretical approach, Bold Diagrammatic Monte Carlo (BDMC), and precision experiments on ultra-cold atoms. Specifically, we compute and measure with unprecedented accuracy the normal-state equation of state of the unitary gas, a prototypical example of a strongly correlated fermionic system. Excellent agreement demonstrates that a series of Feynman diagrams can be controllably resummed in a non-perturbative regime using BDMC. This opens the door to the solution of some of the most challenging problems across many areas of physics. △ Less

Submitted 19 March, 2012; v1 submitted 17 October, 2011; originally announced October 2011.

Journal ref: Nature Phys. 8, 366 (2012)

Diagrammatic Quantum Monte Carlo solution of the two-dimensional Cooperon-Fermion model

Authors: Kai-Yu Yang, E. Kozik, Xin Wang, M. Troyer

Abstract: We investigate the two-dimensional cooperon-fermion model in the correlated regime with a new continuous-time diagrammatic determinant quantum Monte Carlo (DDQMC) algorithm. We estimate the transition temperature $T_{c}$, examine the effectively reduced band gap and cooperon mass, and find that delocalization of the cooperons enhances the diamagnetism. When applied to diamagnetism of the pseudogap… ▽ More We investigate the two-dimensional cooperon-fermion model in the correlated regime with a new continuous-time diagrammatic determinant quantum Monte Carlo (DDQMC) algorithm. We estimate the transition temperature $T_{c}$, examine the effectively reduced band gap and cooperon mass, and find that delocalization of the cooperons enhances the diamagnetism. When applied to diamagnetism of the pseudogap phase in high-$T_{c}$ cuprates, we obtain results in a qualitative agreement with recent torque magnetization measurements. △ Less

Submitted 15 March, 2011; originally announced March 2011.

Comments: 8 pages, 11 figures

Journal ref: Phys. Rev. B 83, 214516 (2011)

Comment on "Dispersive bottleneck delaying thermalization of turbulent Bose-Einstein Condensates" by Krstulovic and Brachet [arXiv:1007.4441]

Authors: Evgeny Kozik

Abstract: We reveal the connection of the recent numerical observations of Krstulovic and Brachet [arXiv:1007.4441] with the general theory of relaxation kinetics of the strongly non-equilibrium Bose-Einstein condensates. We reveal the connection of the recent numerical observations of Krstulovic and Brachet [arXiv:1007.4441] with the general theory of relaxation kinetics of the strongly non-equilibrium Bose-Einstein condensates. △ Less

Submitted 1 September, 2011; v1 submitted 14 February, 2011; originally announced February 2011.

Comments: comment on arXiv:1007.4441, published version, minor stylistic changes

Journal ref: Phys. Rev. Lett. 107, 099601 (2011)

Thermodynamics of the 3D Hubbard model on approach to the Neel transition

Authors: S. Fuchs, E. Gull, L. Pollet, E. Burovski, E. Kozik, T. Pruschke, M. Troyer

Abstract: We study the thermodynamic properties of the 3D Hubbard model for temperatures down to the Neel temperature using cluster dynamical mean-field theory. In particular we calculate the energy, entropy, density, double occupancy and nearest-neighbor spin correlations as a function of chemical potential, temperature and repulsion strength. To make contact with cold-gas experiments, we also compute prop… ▽ More We study the thermodynamic properties of the 3D Hubbard model for temperatures down to the Neel temperature using cluster dynamical mean-field theory. In particular we calculate the energy, entropy, density, double occupancy and nearest-neighbor spin correlations as a function of chemical potential, temperature and repulsion strength. To make contact with cold-gas experiments, we also compute properties of the system subject to an external trap in the local density approximation. We find that an entropy per particle $S/N \approx 0.65(6)$ at $U/t=8$ is sufficient to achieve a Neel state in the center of the trap, substantially higher than the entropy required in a homogeneous system. Precursors to antiferromagnetism can clearly be observed in nearest-neighbor spin correlators. △ Less

Submitted 14 September, 2010; originally announced September 2010.

Comments: 4 pages, 6 figures

Journal ref: Phys. Rev. Lett. 106, 030401 (2011)

The 17/5 spectrum of the Kelvin-wave cascade

Authors: Evgeny Kozik, Boris Svistunov

Abstract: Direct numeric simulation of the Biot-Savart equation readily resolves the 17/5 spectrum of the Kelvin-wave cascade from the 11/3 spectrum of the non-local (in the wavenumber space) cascade scenario by L'vov and Nazarenko. This result is a clear-cut visualisation of the unphysical nature of the 11/3 solution, which was established earlier on the grounds of symmetry. Direct numeric simulation of the Biot-Savart equation readily resolves the 17/5 spectrum of the Kelvin-wave cascade from the 11/3 spectrum of the non-local (in the wavenumber space) cascade scenario by L'vov and Nazarenko. This result is a clear-cut visualisation of the unphysical nature of the 11/3 solution, which was established earlier on the grounds of symmetry. △ Less

Submitted 28 July, 2010; originally announced July 2010.

Comments: 2 pages, 1 figure

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

Authors: Evgeny Kozik, Boris Svistunov

Abstract: 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 b… ▽ More 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. △ Less

Submitted 12 October, 2010; v1 submitted 9 June, 2010; originally announced June 2010.

Comments: comment on arXiv:1005.4575, version accepted in JLTP with minimal changes: abstract added

Journal ref: J. Low. Temp. Phys. 161, 603 (2010)

Geometric Symmetries in Superfluid Vortex Dynamics

Authors: Evgeny Kozik, Boris Svistunov

Abstract: 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… ▽ More 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. △ Less

Submitted 22 August, 2010; v1 submitted 2 June, 2010; originally announced June 2010.

Comments: 4 REVTeX pages, minor stylistic changes, references to recent related preprints added

Journal ref: Phys. Rev. B 82, 140510(R) (2010)

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

Authors: Robert Jördens, Leticia Tarruell, Daniel Greif, Thomas Uehlinger, Niels Strohmaier, Henning Moritz, Tilman Esslinger, Lorenzo De Leo, Corinna Kollath, Antoine Georges, Vito Scarola, Lode Pollet, Evgeni Burovski, Evgeny Kozik, Matthias Troyer

Abstract: 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 t… ▽ More 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. △ Less

Submitted 7 May, 2010; v1 submitted 18 December, 2009; originally announced December 2009.

Comments: 5 pages, 3 figures, 1 table

Journal ref: Phys. Rev. Lett. 104, 180401 (2010)

Diagrammatic Monte Carlo for Correlated Fermions

Authors: E. Kozik, K. Van Houcke, E. Gull, L. Pollet, N. Prokof'ev, B. Svistunov, M. Troyer

Abstract: We show that Monte Carlo sampling of the Feynman diagrammatic series (DiagMC) can be used for tackling hard fermionic quantum many-body problems in the thermodynamic limit by presenting accurate results for the repulsive Hubbard model in the correlated Fermi liquid regime. Sampling Feynman's diagrammatic series for the single-particle self-energy we can study moderate values of the on-site repul… ▽ More We show that Monte Carlo sampling of the Feynman diagrammatic series (DiagMC) can be used for tackling hard fermionic quantum many-body problems in the thermodynamic limit by presenting accurate results for the repulsive Hubbard model in the correlated Fermi liquid regime. Sampling Feynman's diagrammatic series for the single-particle self-energy we can study moderate values of the on-site repulsion ($U/t \sim 4$) and temperatures down to $T/t=1/40$. We compare our results with high temperature series expansion and with single-site and cluster dynamical mean-field theory. △ Less

Submitted 22 January, 2010; v1 submitted 5 July, 2009; originally announced July 2009.

Comments: 4 pages, 5 figures, stylistic changes

Journal ref: EPL 90, 10004 (2010)

Theory of decay of superfluid turbulence in the low-temperature limit

Authors: E. V. Kozik, B. V. Svistunov

Abstract: We review the theory of relaxational kinetics of superfluid turbulence--a tangle of quantized vortex lines--in the limit of very low temperatures when the motion of vortices is conservative. While certain important aspects of the decay depend on whether the tangle is non-structured, like the one in the Kibble-Zurek picture, or essentially polarized, like the one that emulates the Richardson-Kolm… ▽ More We review the theory of relaxational kinetics of superfluid turbulence--a tangle of quantized vortex lines--in the limit of very low temperatures when the motion of vortices is conservative. While certain important aspects of the decay depend on whether the tangle is non-structured, like the one in the Kibble-Zurek picture, or essentially polarized, like the one that emulates the Richardson-Kolmogorov regime of classical turbulence, there are common fundamental features. In both cases, there exists an asymptotic range in the wavenumber space where the energy flux is supported by the cascade of Kelvin waves (kelvons)--precessing distortions propagating along the vortex filaments. At large enough wavenumbers, the Kelvin-wave cascade is supported by three-kelvon elastic scattering. At T=0 the dissipative cutoff of the Kelvin-wave cascade is due to the emission of phonons, in which an elementary process converts two kelvons with almost opposite momenta into one bulk phonon. Along with the standard set of conservation laws, a crucial role in the theory is played by the fact of integrability of the local induction approximation (LIA) controlled by the parameter Λらむだ= \ln (λらむだ/a_0), with λらむだthe characteristic kelvon wavelength and a_0 the vortex core radius. While excluding a straightforward onset of the pure three-kelvon cascade, the integrability of LIA does not plug the cascade because of the availability of the kinetic channels associated with vortex line reconnections. We argue that the crossover from Richardson-Kolmogorov to the Kelvin-wave cascade is due to eventual dominance of local induction of a single line over the collective induction of polarized eddies, which causes the breakdown of classical-fluid regime and gives rise to a reconnection-driven inertial range. △ Less

Submitted 8 April, 2009; originally announced April 2009.

Comments: Review article, 52 pages, 8 figures

Journal ref: J. Low Temp. Phys. 156, 215 (2009)

Critical Temperature Curve in the BEC-BCS Crossover

Authors: Evgeni Burovski, Evgeny Kozik, Nikolay Prokof'ev, Boris Svistunov, Matthias Troyer

Abstract: The strongly-correlated regime of the BCS-BEC crossover can be realized by diluting a system of two-component fermions with a short-range attractive interaction. We investigate this system via a novel continuous-space-time diagrammatic determinant Monte Carlo method and determine the universal curve $T_c/εいぷしろん_F$ for the transition temperature between the normal and the superfluid states as a functi… ▽ More The strongly-correlated regime of the BCS-BEC crossover can be realized by diluting a system of two-component fermions with a short-range attractive interaction. We investigate this system via a novel continuous-space-time diagrammatic determinant Monte Carlo method and determine the universal curve $T_c/εいぷしろん_F$ for the transition temperature between the normal and the superfluid states as a function of the scattering length with the maximum on the BEC side. At unitarity, we confirm that $T_c/εいぷしろん_F = 0.152(7)$. △ Less

Submitted 4 September, 2008; v1 submitted 20 May, 2008; originally announced May 2008.

Comments: revtex, 4 pages, 3 figures, published version

Journal ref: Phys. Rev. Lett. 101, 090402 (2008)

Diagrammatic Monte Carlo

Authors: Kris Van Houcke, Evgeny Kozik, Nikolay Prokof'ev, Boris Svistunov

Abstract: Diagrammatic Monte Carlo (DiagMC) is a numeric technique that allows one to calculate quantities specified in terms of diagrammatic expansions, the latter being a standard tool of many-body quantum statistics. The sign problem that is typically fatal to Monte Carlo approaches, appears to be manageable with DiagMC. Starting with a general introduction to the principles of DiagMC, we present a det… ▽ More Diagrammatic Monte Carlo (DiagMC) is a numeric technique that allows one to calculate quantities specified in terms of diagrammatic expansions, the latter being a standard tool of many-body quantum statistics. The sign problem that is typically fatal to Monte Carlo approaches, appears to be manageable with DiagMC. Starting with a general introduction to the principles of DiagMC, we present a detailed description of the DiagMC scheme for interacting fermions (Hubbard model), as well as the first illustrative results for the equations of state. △ Less

Submitted 20 February, 2008; originally announced February 2008.

Comments: 15 pages, 7 figures. To appear in: Computer Simulation Studies in Condensed Matter Physics XXI, Eds. D.P. Landau, S.P. Lewis, and H.B. Schuttler (Springer Verlag, Heidelberg, Berlin 2008)

Journal ref: Phys. Procedia 6, 95 (2010)

Scanning Superfluid-Turbulence Cascade by Its Low-Temperature Cutoff

Authors: Evgeny Kozik, Boris Svistunov

Abstract: On the basis of recently proposed scenario of the transformation of the Kolmogorov cascade into the Kelvin-wave cascade, we develop a theory of low-temperature cutoff. The theory predicts a specific behavior of the quantized vortex line density, $L$, controlled by the frictional coefficient, $αあるふぁ(T) \ll 1$, responsible for the cutoff. The curve $\ln L(\ln αあるふぁ)$ is found to directly reflect the struc… ▽ More On the basis of recently proposed scenario of the transformation of the Kolmogorov cascade into the Kelvin-wave cascade, we develop a theory of low-temperature cutoff. The theory predicts a specific behavior of the quantized vortex line density, $L$, controlled by the frictional coefficient, $αあるふぁ(T) \ll 1$, responsible for the cutoff. The curve $\ln L(\ln αあるふぁ)$ is found to directly reflect the structure of the cascade, revealing four qualitatively distinct wavenumber regions. Excellent agreement with recent experiment by Walmsley {\it et al.} [arXiv:0710.1033]--in which $L(T)$ has been measured down to $T \sim 0.08 $K--implies that the scenario of low-temperature superfluid turbulence is now experimentally validated, and allows to quantify the Kelvin-wave cascade spectrum. △ Less

Submitted 2 March, 2008; v1 submitted 24 October, 2007; originally announced October 2007.

Comments: 4 pages, 2 figures, v2: extended introduction, the controversy with the scenario by L'vov et al. [13] is discussed in conclusion

Journal ref: Phys. Rev. Lett. 100, 195302 (2008)

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

Authors: Evgeny Kozik, Boris Svistunov

Abstract: 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 K… ▽ More 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. △ Less

Submitted 24 October, 2007; v1 submitted 1 March, 2007; originally announced March 2007.

Comments: 4 RevTex pages, 1 figure. Quantitative analysis of the regime 2 has been revised

Journal ref: Phys. Rev. B 77, 060502 (2008)

On-site number statistics of ultracold lattice bosons

Authors: Barbara Capogrosso-Sansone, Evgeny Kozik, Nikolay Prokof'ev, Boris Svistunov

Abstract: We study on-site occupation number fluctuations in a system of interacting bosons in an optical lattice. The ground-state distribution is obtained analytically in the limiting cases of strong and weak interaction, and by means of exact Monte Carlo simulations in the strongly correlated regime. As the interaction is increased, the distribution evolves from Poissonian in the non-interacting gas to… ▽ More We study on-site occupation number fluctuations in a system of interacting bosons in an optical lattice. The ground-state distribution is obtained analytically in the limiting cases of strong and weak interaction, and by means of exact Monte Carlo simulations in the strongly correlated regime. As the interaction is increased, the distribution evolves from Poissonian in the non-interacting gas to a sharply peaked distribution in the Mott-insulator (MI) regime. In the special case of large occupation numbers, we demonstrate analytically and check numerically that there exists a wide interval of interaction strength, in which the on-site number fluctuations remain Gaussian and are gradually squeezed until they are of order unity near the superfluid (SF)-MI transition. Recently, the on-site number statistics were studied experimentally in a wide range of lattice potential depths [Phys. Rev. Lett. \textbf{96}, 090401 (2006)]. In our simulations, we are able to directly reproduce experimental conditions using temperature as the only free parameter. Pronounced temperature dependence suggests that measurements of on-site atom number fluctuations can be employed as a reliable method of thermometry in both SF and MI regimes. △ Less

Submitted 23 September, 2006; originally announced September 2006.

Comments: 9 pages, 4 figures

Journal ref: Phys. Rev. A 75, 013619 (2007)

Vortex-Phonon Interaction in the Kosterlitz-Thouless Theory

Authors: Evgeny Kozik, Nikolay Prokof'ev, Boris Svistunov

Abstract: 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-… ▽ More 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. △ Less

Submitted 27 October, 2010; v1 submitted 4 November, 2005; originally announced November 2005.

Comments: 4 pages, RevTex

Journal ref: Phys. Rev. B 73, 092501 (2006)

Vortex--Phonon Interaction

Authors: Evgeny Kozik, Boris Svistunov

Abstract: 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… ▽ More 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. △ Less

Submitted 10 November, 2005; v1 submitted 1 May, 2005; originally announced May 2005.

Comments: 4 pages, RevTex, a remark regarding elastic and inelastic phonon scattering is added, to be published in PRB

Journal ref: Phys. Rev. B 72, 172505 (2005)

Superfluid Interfaces in Quantum Solids

Authors: Evgeni Burovski, Evgeni Kozik, Anatoly Kuklov, Nikolay Prokof'ev, Boris Svistunov

Abstract: One scenario for the non-classical moment of inertia of solid He-4 discovered by Kim and Chan [Nature 427, 225 (2004)] is the superfluidity of micro-crystallite interfaces. On the basis of the most simple model of a quantum crystal--the checkerboard lattice solid--we show that the superfluidity of interfaces between solid domains can exist in a wide range of parameters. At strong enough inter-pa… ▽ More One scenario for the non-classical moment of inertia of solid He-4 discovered by Kim and Chan [Nature 427, 225 (2004)] is the superfluidity of micro-crystallite interfaces. On the basis of the most simple model of a quantum crystal--the checkerboard lattice solid--we show that the superfluidity of interfaces between solid domains can exist in a wide range of parameters. At strong enough inter-particle interaction, a superfluid interface becomes an insulator via a quantum phase transition. Under the conditions of particle-hole symmetry, the transition is of the standard U(1) universality class in 3D, while in 2D the onset of superfluidity is accompanied by the interface roughening, driven by fractionally charged topological excitations. △ Less

Submitted 22 December, 2004; originally announced December 2004.

Comments: 4 revtex4 pages

Journal ref: Phys. Rev. Lett., vol. 94, p. 165301 (2005)

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

Authors: Evgeny Kozik, Boris Svistunov

Abstract: 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… ▽ More 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). △ Less

Submitted 8 October, 2004; v1 submitted 11 August, 2004; originally announced August 2004.

Comments: 4 pages, RevTex

Journal ref: Phys. Rev. Lett. 94, 025301 (2005)

Kelvin Wave Cascade and Decay of Superfluid Turbulence

Authors: Evgeny Kozik, Boris Svistunov

Abstract: Kelvin waves (kelvons)--the distortion waves on vortex lines--play a key part in the relaxation of superfluid turbulence at low temperatures. We present a weak-turbulence theory of kelvons. We show that non-trivial kinetics arises only beyond the local-induction approximation and is governed by three-kelvon collisions; corresponding kinetic equation is derived. On the basis of the kinetic equati… ▽ More Kelvin waves (kelvons)--the distortion waves on vortex lines--play a key part in the relaxation of superfluid turbulence at low temperatures. We present a weak-turbulence theory of kelvons. We show that non-trivial kinetics arises only beyond the local-induction approximation and is governed by three-kelvon collisions; corresponding kinetic equation is derived. On the basis of the kinetic equation, we prove the existence of Kolmogorov cascade and find its spectrum. The qualitative analysis is corroborated by numeric study of the kinetic equation. The application of the results to the theory of superfluid turbulence is discussed. △ Less

Submitted 10 August, 2003; originally announced August 2003.

Comments: 4 pages, RevTeX