-
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.
-
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 10 September, 2024; v1 submitted 24 September, 2023;
originally announced September 2023.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Strong and Electromagnetic Forces in Heavy Ion Collisions
Authors:
Mariola Klusek-Gawenda,
Ewa Kozik,
Andrzej Rybicki,
Iwona Sputowska,
Antoni Szczurek
Abstract:
The interplay between the strong and electromagnetic force in high energy nucleus-nucleus collisions was studied experimentally and theoretically in our earlier works. This effect appeared to result in very large distortions in spectra of charged pions produced in the collision. It was also found to bring new, independent information on the space-time evolution of the non-perturbative process of p…
▽ More
The interplay between the strong and electromagnetic force in high energy nucleus-nucleus collisions was studied experimentally and theoretically in our earlier works. This effect appeared to result in very large distortions in spectra of charged pions produced in the collision. It was also found to bring new, independent information on the space-time evolution of the non-perturbative process of particle production.
In this paper, we present our new results on the influence of the spectator-induced electromagnetic force on spectra of charged particles produced in two different Pb-induced reactions. For the first time, we also address the topic of p+A collisions in view of obtaining information about their centrality and nuclear break-up, both subjects being of importance in the context of the new p+A data collected at the LHC.
△ Less
Submitted 26 March, 2013;
originally announced March 2013.
-
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.
-
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.
-
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.
-
Sequential and simultaneous emission of fragments from p+Al collisions
Authors:
M. Fidelus,
D. Filges,
F. Goldenbaum,
H. Hodde,
A. Jany,
L. Jarczyk,
B. Kamys,
M. Kistryn,
St. Kistryn,
St. Kliczewski,
E. Kozik,
P. Kulessa,
H. Machner,
A. Magiera,
B. Piskor-Ignatowicz,
K. Pysz,
Z. Rudy,
Sushil K. Sharma,
R. Siudak,
M. Wojciechowski
Abstract:
The energy and angular dependencies of double differential cross sections d^2sigma/dOmegadE were measured for p,d,t,^{3,4,6}He,^{6,7,8}Li,{7,9,10}Be, and ^{10,11,12}B produced in collisions of 1.2, 1.9, and 2.5 GeV protons with an Al target. The spectra and angular distributions of Li, Be, and B isotopes indicate a presence of two contributions: an isotropic, low energy one which is attributed to…
▽ More
The energy and angular dependencies of double differential cross sections d^2sigma/dOmegadE were measured for p,d,t,^{3,4,6}He,^{6,7,8}Li,{7,9,10}Be, and ^{10,11,12}B produced in collisions of 1.2, 1.9, and 2.5 GeV protons with an Al target. The spectra and angular distributions of Li, Be, and B isotopes indicate a presence of two contributions: an isotropic, low energy one which is attributed to the evaporation of particles from excited remnants of the intranuclear cascade, and an anisotropic part which is interpreted to be due to multifragmentation of the remnants. It was found that such a hypothesis leads to very good description of spectra and angular distributions of all intermediate mass fragments (^6He-^{12}B) using the critical value of the excitation energy per nucleon as a single parameter, varying slowly with the beam energy.
△ Less
Submitted 14 September, 2011;
originally announced September 2011.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Evidence of Ni break-up from total production cross sections in p+Ni collisions
Authors:
A. Budzanowski,
M. Fidelus,
D. Filges,
F. Goldenbaum,
H. Hodde,
L. Jarczyk,
B. Kamys,
M. Kistryn,
St. Kistryn,
St. Kliczewski,
A. Kowalczyk,
E. Kozik,
P. Kulessa,
H. Machner,
A. Magiera,
B. Piskor-Ignatowicz,
K. Pysz,
Z. Rudy,
R. Siudak,
M. Wojciechowski
Abstract:
The total production cross sections of light charged particles (LCPs), intermediate mass fragments (IMFs) and heavy reaction products of p+Ni collisions available in the literature have been compared with predictions of a two-step model in the proton beam energy range from reaction threshold up to approximately 3 GeV. Model cross sections were calculated assuming, that the reaction proceeds via…
▽ More
The total production cross sections of light charged particles (LCPs), intermediate mass fragments (IMFs) and heavy reaction products of p+Ni collisions available in the literature have been compared with predictions of a two-step model in the proton beam energy range from reaction threshold up to approximately 3 GeV. Model cross sections were calculated assuming, that the reaction proceeds via an intranuclear cascade of nucleon-nucleon collisions followed by evaporation of particles from an equilibrated, heavy target residuum. The shape of the excitation functions was well described by model calculations for all reaction products. The magnitude of the cross sections was reasonably well reproduced for heavy reaction products, i.e. for nuclei heavier than Al, but the cross sections for lighter products were systematically underestimated. This fact was used as an argument in favor of a significant break-up contribution to the reaction mechanism. The present conclusions are supported by recently published results of investigations of differential cross sections in p+Ni collisions, which showed that hypothesis of the break-up of target nucleus is indispensable for a good reproduction of d2 sigma/dOmega dE for LCPs and IMFs.
△ Less
Submitted 9 October, 2009;
originally announced October 2009.
-
Variation of nonequilibrium processes in p+Ni system with beam energy
Authors:
A. Budzanowski,
M. Fidelus,
D. Filges,
F. Goldenbaum,
H. Hodde,
L. Jarczyk,
B. Kamys,
M. Kistryn,
St. Kistryn,
St. Kliczewski,
A. Kowalczyk,
E. Kozik,
P. Kulessa,
H. Machner,
A. Magiera,
B. Piskor-Ignatowicz,
K. Pysz,
Z. Rudy,
R. Siudak,
M. Wojciechowski
Abstract:
The energy and angular dependence of double differential cross sections dsigma/dOmega dE were measured for p, d, t, 3,4He, 6,7Li, 7,9Be, and 10,11B produced in collisions of 0.175 GeV protons with Ni target. The analysis of measured dfferential cross sections allowed to extract total production cross sections for ejectiles listed above. The shape of the spectra and angular distributions indicate…
▽ More
The energy and angular dependence of double differential cross sections dsigma/dOmega dE were measured for p, d, t, 3,4He, 6,7Li, 7,9Be, and 10,11B produced in collisions of 0.175 GeV protons with Ni target. The analysis of measured dfferential cross sections allowed to extract total production cross sections for ejectiles listed above. The shape of the spectra and angular distributions indicate the presence of other nonequilibrium processes besides the emission of nucleons from the intranuclear cascade, and besides the evaporation of various particles from remnants of intranuclear cascade. These nonequilibrium processes consist of coalescence of nucleons into light charged particles during the intranuclear cascade, of the fireball emission which contributes to the cross sections of protons and deuterons, and of the break-up of the target nucleus which leads to the emission of intermediate mass fragments. All such processes were found earlier at beam energies 1.2, 1.9, and 2.5 GeV for Ni as well as for Au targets, however, significant differences in properties of these processes at high and low beam energy are observed in the present study.
△ Less
Submitted 8 September, 2009;
originally announced September 2009.
-
Comparison of nonequilibrium processes in p+Ni and p+Au collisions at GeV energies
Authors:
A. Budzanowski,
M. Fidelus,
D. Filges,
F. Goldenbaum,
H. Hodde,
L. Jarczyk,
B. Kamys,
M. Kistryn,
St. Kistryn,
St. Kliczewski,
A. Kowalczyk,
E. Kozik,
P. Kulessa,
H. Machner,
A. Magiera,
B. Piskor-Ignatowicz,
K. Pysz,
Z. Rudy,
R. Siudak,
M. Wojciechowski
Abstract:
The energy and angular dependence of double differential cross sections d2sigma/dOmega dE were measured for p, d, t, 3,4,6He, 6,7,8Li, 7,9,10Be, 10,11B, and C produced in collisions of 1.2, 1.9, and 2.5 GeV protons with a Ni target. The shape of the spectra and angular distributions does almost not change whereas the absolute value of the cross sections increases by a factor about 1.7 for all ej…
▽ More
The energy and angular dependence of double differential cross sections d2sigma/dOmega dE were measured for p, d, t, 3,4,6He, 6,7,8Li, 7,9,10Be, 10,11B, and C produced in collisions of 1.2, 1.9, and 2.5 GeV protons with a Ni target. The shape of the spectra and angular distributions does almost not change whereas the absolute value of the cross sections increases by a factor about 1.7 for all ejectiles in this beam energy range. It was found that energy and angular dependencies of the cross sections cannot be reproduced by the microscopic model of intranuclear cascade with coalescence of nucleons and the statistical model for evaporation of particles from excited, equilibrated residual nuclei. The inclusion of nonequilibrium processes, described by a phenomenological model of the emission from fast and hot moving sources, resulting from break-up of the target nucleus by impinging proton, leads to very good reproduction of data. Cross sections of these processes are quite large, exhausting approximately half of the total production cross sections. Due to good reproduction of energy and angular dependencies of d2sigma/dOmega dE by model calculation it was possible to determine total production cross sections for all studied ejectiles. Results obtained in this work point to the analogous reaction mechanism for proton induced reactions on Ni target as that observed previously for Au target in the same beam energy range.
△ Less
Submitted 31 August, 2009;
originally announced August 2009.
-
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.
-
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.
-
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.
-
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.
-
Competition of coalescence and "fireball" processes in nonequilibrium emission of light charged particles from p+Au collisions
Authors:
A. Budzanowski,
M. Fidelus,
D. Filges,
F. Goldenbaum,
H. Hodde,
L. Jarczyk,
B. Kamys,
M. Kistryn,
St. Kistryn,
St. Kliczewski,
A. Kowalczyk,
E. Kozik,
P. Kulessa,
H. Machner,
A. Magiera,
B. Piskor-Ignatowicz,
K. Pysz,
Z. Rudy,
R. Siudak,
M. Wojciechowski
Abstract:
The energy and angular dependence of double differential cross sections was measured for p,d,t,He,Li,Be, and B isotopes produced in collisions of 1.2 and 1.9 GeV protons with Au target. The shape of the spectra and angular distributions almost does not change in the beam energy range from 1.2 to 2.5 GeV, however, the absolute value of the cross sections increases for all ejectiles. A phenomenolo…
▽ More
The energy and angular dependence of double differential cross sections was measured for p,d,t,He,Li,Be, and B isotopes produced in collisions of 1.2 and 1.9 GeV protons with Au target. The shape of the spectra and angular distributions almost does not change in the beam energy range from 1.2 to 2.5 GeV, however, the absolute value of the cross sections increases for all ejectiles. A phenomenological model of two emitting, moving sources reproduces very well spectra and angular distributions of intermediate mass fragments. Double differential cross sections for light charged particles (LCP) were analyzed in the frame of the microscopic model of intranuclear cascade (INC) with coalescence of nucleons and statistical model for evaporation of particles from excited residual nuclei. Energy and angular dependencies of data agree satisfactorily neither with predictions of microscopic intranuclear cascade calculations for protons, nor with coalescence calculations for other LCP. Phenomenological inclusion of another reaction mechanism - emission of LCP from a "fireball", i.e., fast and hot moving source - combined with the microscopic model calculations of INC, coalescence and evaporation of particles leads to very good description of the data. It was found that nonequilibrium processes are very important for production of LCP. They exhaust 40-80% of the total cross sections - depending on the emitted particles. Coalescence and "fireball" emission give comparable contributions to the cross sections with exception of 3He data where coalescence clearly dominates. The ratio of sum of all nonequilibrium processes to those proceeding through stage of statistical equilibrium does almost not change in the beam energy range from 1.2 GeV to 2.5 GeV for all light charged particles.
△ Less
Submitted 29 January, 2008;
originally announced January 2008.
-
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.
-
Non-equilibrium emission of complex fragments from p+Au collisions at 2.5 GeV proton beam energy
Authors:
A. Budzanowski,
D. Filges,
F. Goldenbaum,
A. Heczko,
H. Hodde,
L. Jarczyk,
B. Kamys,
M. Kistryn,
St. Kistryn,
St. Kliczewski,
A. Kowalczyk,
M. Smoluchowski,
E. Kozik,
P. Kulessa,
H. Machner,
A. Magiera,
W. Migdal,
N. Paul,
B. Piskor-Ignatowicz,
M. Puchala,
K. Pysz,
Z. Rudy,
R. Siudak,
M. Wojciechowski
Abstract:
Energy and angular dependence of double differential cross sections d$^2σ$/d$Ω$dE was measured for reactions induced by 2.5 GeV protons on Au target with isotopic identification of light products (H, He, Li, Be, and B) and with elemental identification of heavier intermediate mass fragments (C, N, O, F, Ne, Na, Mg, and Al). It was found that two different reaction mechanisms give comparable cont…
▽ More
Energy and angular dependence of double differential cross sections d$^2σ$/d$Ω$dE was measured for reactions induced by 2.5 GeV protons on Au target with isotopic identification of light products (H, He, Li, Be, and B) and with elemental identification of heavier intermediate mass fragments (C, N, O, F, Ne, Na, Mg, and Al). It was found that two different reaction mechanisms give comparable contributions to the cross sections. The intranuclear cascade of nucleon-nucleon collisions followed by evaporation from an equilibrated residuum describes low energy part of the energy distributions whereas another reaction mechanism is responsible for high energy part of the spectra of composite particles. Phenomenological model description of the differential cross sections by isotropic emission from two moving sources led to a very good description of all measured data. Values of the extracted parameters of the emitting sources are compatible with the hypothesis claiming that the high energy particles emerge from pre-equilibrium processes consisting in a breakup of the target into three groups of nucleons; small, fast and hot fireball of $\sim$ 8 nucleons, and two larger, excited prefragments, which emits the light charged particles and intermediate mass fragments. The smaller of them contains $\sim$ 20 nucleons and moves with velocity larger than the CM velocity of the proton projectile and the target. The heavier prefragment behaves similarly as the heavy residuum of the intranuclear cascade of nucleon-nucleon collisions. %The mass and charge dependence of the total production cross %sections was extracted from the above analysis for all observed %reaction products. This dependence follows the power low behavior %(A$^{-τ}$ or Z$^{-τ}$).
△ Less
Submitted 18 September, 2007;
originally announced September 2007.
-
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.
-
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.
-
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.
-
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.