Peter Kopietz

Accumulation of magnetoelastic bosons in yttrium iron garnet: kinetic theory and wave vector resolved Brillouin light scattering

We derive and solve quantum kinetic equations describing the accumulation of magnetoelastic bosons in an overpopulated magnon gas realized in a thin film of the magnetic insulator yttrium iron garnet. We show that in the presence of a magnon condensate, there is a non-equilibrium steady state in which incoherent magnetoelastic bosons accumulate in a narrow region in momentum space for energies slightly below the bottom of the magnon spectrum. The results of our calculations agree quite well with Brillouin light scattering measurements of the stationary non-equilibrium state of magnons and magnetoelastic bosons in yttrium iron garnet.

Critical spin dynamics of Heisenberg ferromagnets revisited

We calculate the dynamic structure factor $S (\boldsymbol{k},ω)$ in the paramagnetic regime of quantum Heisenberg ferromagnets for temperatures $T$ close to the critical temperature $T_c$ using our recently developed functional renormalization group approach to quantum spin systems. In $d=3$ dimensions we find that for small momenta $\boldsymbol{k}$ and frequencies $ω$ the dynamic structure factor assumes the scaling form $S(\boldsymbol{k},ω) = (τT G (\boldsymbol{k})/π)Φ(kξ, ωτ)$ , where $G (\boldsymbol{k})$ is the static spin-spin correlation function, $ξ$ is the correlation length, and the characteristic time-scale $τ$ is proportional to $ξ^{5/2}$ . We explicitly calculate the dynamic scaling function $Φ(x,y)$ and find satisfactory agreement with neutron scattering experiments probing the critical spin dynamics in EuO and EuS. Precisely at the critical point where $ξ= \infty$ our result for the dynamic structure factor can be written as $S (\boldsymbol{k},ω) = (πω_k)^{-1} T_c G (\boldsymbol{k}) Ψ_c (ω/ω_k)$ , where $ω_k \propto k^{5/2}$ . We find that $Ψ_c(ν)$ vanishes as $ν^{-13/5}$ for large $ν$ , and as $ν^{3/5}$ for small $ν$ . While the large-frequency behavior of $Ψ_c (ν)$ is consistent with calculations based on mode-coupling theory and with perturbative renormalization group calculations to second order in $ε= 6-d$ , our result for small frequencies disagrees with previous calculations. We argue that up until now neither experiments nor numerical simulations are sufficiently accurate to determine the low-frequency behavior of $Ψ_c (ν)$ . We also calculate the low-temperature behavior of $S ( \boldsymbol{k},ω)$ in one- and two dimensional ferromagnets and find that it satisfies dynamic scaling with exponent $z=2$ and exhibits a pseudogap for small frequencies.

Phonon renormalization and Pomeranchuk instability in the Holstein model

The Holstein model with dispersionless Einstein phonons is one of the simplest models describing electron-phonon interactions in condensed matter. A naive extrapolation of perturbation theory in powers of the relevant dimensionless electron-phonon coupling $λ_0$ suggests that at zero temperature the model exhibits a Pomeranchuk instability characterized by a divergent uniform compressibility at a critical value of $λ_0$ of order unity. In this work, we re-examine this problem using modern functional renormalization group (RG) methods. For dimensions $d > 3$ we find that the RG flow of the Holstein model indeed exhibits a tricritical fixed point associated with a Pomeranchuk instability. This non-Gaussian fixed point is ultraviolet stable and is closely related to the well-known ultraviolet stable fixed point of $ϕ^3$ -theory above six dimensions. To realize the Pomeranchuk critical point in the Holstein model at fixed density both the electron-phonon coupling $λ_0$ and the adiabatic ratio $ω_0 / ε_F$ have to be fine-tuned to assume critical values of order unity, where $ω_0$ is the phonon frequency and $ε_F$ is the Fermi energy. On the other hand, for dimensions $d \leq 3$ we find that the RG flow of the Holstein model does not have any critical fixed points. This rules out a quantum critical point associated with a Pomeranchuk instability in $d \leq 3$ .

Spin functional renormalization group for dimerized quantum spin systems

We investigate dimerized quantum spin systems using the spin functional renormalization group approach proposed by Krieg and Kopietz [Phys. Rev. B 99, 060403(R) (2019)] which directly focuses on the physical spin correlation functions and avoids the representation of the spins in terms of fermionic or bosonic auxiliary operators. Starting from decoupled dimers as initial condition for the renormalization group flow equations, we obtain the spectrum of the triplet excitations as well as the magnetization in the quantum paramagnetic, ferromagnetic, and thermally disordered phases at all temperatures. Moreover, we compute the full phase diagram of a weakly coupled dimerized spin system in three dimensions, including the correct mean field critical exponents at the two quantum critical points.

preprint2021arXiv

Zero-magnon sound in quantum Heisenberg ferromagnets

Using a functional renormalization group approach we show that at low temperatures and in a certain range of magnetic fields the longitudinal dynamic structure factor of quantum Heisenberg ferromagnets in dimensions $D\leq 2$ exhibits a well-defined quasi-particle peak with linear dispersion that we identify as zero-magnon sound. In $D>2$ the larger phase space available for the decay into transverse spin waves leads only to a broad hump centered at zero frequency whose width scales linearly in momentum.

preprint2020arXiv

Collisionless kinetic theory for parametrically pumped magnons

We discuss collisionless kinetic equations describing the non-equilibrium dynamics of magnons in a ferromagnet exposed to an oscillating microwave field. Previously, this problem has been treated within the so-called "S-theory" where the collision integral in the kinetic equation for the magnon distribution is either neglected or taken into account phenomenologically via an effective relaxation time. However, the possibility of magnon condensation has not been included in S-theory. Moreover, the momentum integrations appearing in the magnon self-energies are usually decoupled by retaining only the term where the loop momentum is equal to the external momentum. In this work we critically examine the accuracy of these approximations and develop the proper extensions of S-theory. We show that these extensions can significantly modify the time evolution of the magnon distribution.

preprint2016arXiv

Multilogarithmic velocity renormalization in graphene

We reexamine the effect of long-range Coulomb interactions on the quasiparticle velocity in graphene. Using a nonperturbative functional renormalization group approach with partial bosonization in the forward scattering channel and momentum transfer cutoff scheme, we calculate the quasiparticle velocity, $v (k)$ , and the quasiparticle residue, $Z$ , with frequency-dependent polarization. One of our most striking results is that $v ( k ) \propto \ln [ C_k (α) / k ]$ where the momentum- and interaction-dependent cutoff scale $C_k (α)$ vanishes logarithmically for $k \rightarrow 0$ . Here $k$ is measured with respect to one of the charge neutrality (Dirac) points and $α=2.2$ is the strength of dimensionless bare interaction. Moreover, we also demonstrate that the so-obtained multilogarithmic singularity is reconcilable with the perturbative expansion of $v (k)$ in powers of the bare interaction.

preprint2016arXiv

Physical dipoles and second order perturbation theory for dipolar fermions in two dimensions

In two dimensions the Fourier transform of the interaction between two point dipoles has a term which grows linearly in the modulus $| \mathbf{\textit{q}} |$ of the momentum . As a consequence, in second order perturbation theory the self-energy of two-dimensional dipolar fermions is ultraviolet divergent. We show that for electric dipoles this divergence can be avoided if one takes into account that physical dipoles consist of two opposite charges which are separated by a finite distance. Using this regularization, we calculate the self-energy, the renormalized chemical potential, and the renormalized Fermi surface of dipolar fermions in two dimensions in second order perturbation theory. We find that in the Fermi liquid phase the second order corrections weaken first order effects.

Hard-core boson approach to the spin-1/2 triangular-lattice antiferromagnet Cs $_2$ CuCl $_4$ at finite temperatures in magnetic fields higher than the saturation field

We study the high magnetic field regime of the antiferromagnetic insulator Cs $_2$ CuCl $_4$ by expressing the spin-1/2 operators in the relevant Heisenberg model in terms of hard-core bosons and implementing the hard-core constraint via an infinite on-site interaction. We focus on the case where the external magnetic field exceeds the saturation field $B_{c}\approx8.5\;\mathrm{T}$ and is oriented along the crystallographic $a$ axis perpendicular to the lattice plane. Because in this case the excited states are separated by an energy gap from the ground state, we may use the self-consistent ladder approximation to take the strong correlations due to the hard-core constraint into account. In Cs $_2$ CuCl $_4$ there are additional interactions besides the hard-core interaction which we treat in self-consistent Hartree-Fock approximation. We calculate the spectral function of the hard-core bosons from which we obtain the in-plane components of the dynamic structure factor, the magnetic susceptibility, and the specific heat. Our results for the specific heat are in good agreement with the available experimental data. We conclude that the self-consistent ladder approximation in combination with a self-consistent Hartree-Fock decoupling of the non-hard-core interactions gives an accurate description of the physical properties of gapped hard-core bosons in two dimensions at finite temperatures.

Magneto-caloric effects, quantum critical points, and the Berezinsky-Kosterlitz-Thouless transition in 2D coupled spin dimer systems

Spin dimer systems are a promising playground for the detailed study of quantum phase transitions. Using the magnetic field as the tuning parameter it is in principle possible to observe a crossover from the characteristic scaling near critical points to the behavior of a finite temperature phase transition. In this work we study two-dimensional coupled spin dimer systems by comparing numerical quantum Monte Carlo simulations with analytical calculations of the susceptibility, the magneto-caloric effect, and the helicity modulus. The magneto-caloric behavior of the magnetization with temperature can be used to determine the critical fields with high accuracy, but the critical scaling does not show the expected logarithmic corrections. The zeros of the cooling rate are an excellent indicator of the competition between quantum criticality and vortex physics, but they are not directly associated with the quantum phase transition or the finite temperature Berezinsky-Kosterlitz-Thouless transition. The results give a unified picture of the full quantum and finite temperature phase diagram.

Non-analytic magnetic field dependence of quasi-particle properties of two-dimensional metals

We show that in a weak external magnetic field H the quasi-particle residue and the renormalized electron Lande factor of two-dimensional Fermi liquids exhibit a non-analytic magnetic field dependence proportional to |H| which is due to electron-electron interactions. We explicitly calculate the corresponding prefactors to second order in the interaction and show that they are determined by low-energy scattering processes involving only momenta close to the Fermi surface. Experimentally, these non-analytic terms can be detected from measurements of the magnetic field dependence of the density of states and the magnetoconductivity.

Non-perturbative renormalization group calculation of the quasi-particle velocity and the dielectric function of graphene

Using a non-perturbative functional renormalization group approach we calculate the renormalized quasi-particle velocity $v (k)$ and the static dielectric function $ε( k )$ of suspended graphene as functions of an external momentum $k$ . Our numerical result for $v (k )$ can be fitted by $v ( k ) / v_F = A + B \ln ( Λ_0 / k )$ , where $v_F$ is the bare Fermi velocity, $Λ_0$ is an ultraviolet cutoff, and $A = 1.37$ , $B =0.51$ for the physically relevant value ( $e^2/v_F =2.2$ ) of the coupling constant. In contrast to calculations based on the static random-phase approximation, we find that $ε(k )$ approaches unity for $k \rightarrow 0$ . Our result for $v (k )$ agrees very well with a recent measurement by Elias et al. [Nat. Phys. 7, 701 (2011)].

Rayleigh-Jeans condensation of pumped magnons in thin film ferromagnets

We show that the formation of a magnon condensate in thin ferromagnetic films can be explained within the framework of a classical stochastic non-Markovian Landau-Lifshitz-Gilbert equation where the properties of the random magnetic field and the dissipation are determined by the underlying phonon dynamics. We have numerically solved this equation for a tangentially magnetized yttrium-iron garnet film in the presence of a parallel parametric pumping field. We obtain a complete description of all stages of the nonequilibrium time evolution of the magnon gas which is in excellent agreement with experiments. Our calculation proves that the experimentally observed condensation of magnons in yttrium-iron garnet at room temperature is a purely classical phenomenon which should be called Rayleigh-Jeans rather than Bose-Einstein condensation.

Second-order interaction corrections to the Fermi surface and the quasiparticle properties of dipolar fermions in three dimensions

We calculate the renormalized Fermi surface and the quasiparticle properties in the Fermi liquid phase of three-dimensional dipolar fermions to second order in the dipole-dipole interaction. Using parameters relevant to an ultracold gas of erbium atoms, we find that the second-order corrections typically renormalize the Hartree-Fock results by less than one percent. On the other hand, if we use the second-order correction to the compressibility to estimate the regime of stability of the system, the point of instability is already reached for a significantly smaller interaction strength than in the Hartree-Fock approximation.

Summing parquet diagrams using the functional renormalization group: X-ray problem revisited

We present a simple method for summing so-called parquet diagrams of fermionic many-body systems with competing instabilities using the functional renormalization group. Our method is based on partial bosonization of the interaction utilizing multi-channel Hubbard-Stratonovich transformations. A simple truncation of the resulting flow equations, retaining only the frequency-independent parts of the two-point and three-point vertices amounts to solving coupled Bethe-Salpeter equations for the effective interaction to leading logarithmic order. We apply our method by revisiting the X-ray problem and deriving the singular frequency dependence of the X-ray response function and the particle-particle susceptibility. Our method is quite general and should be useful in many-body problems involving strong fluctuations in several scattering channels.

preprint2014arXiv

Magneto-elastic modes and lifetime of magnons in thin yttrium-iron garnet films

We calculate the effects of the spin-lattice coupling on the magnon spectrum of thin ferromagnetic films consisting of the magnetic insulator yttrium-iron garnet. The magnon-phonon hybridisation generates a characteristic minimum in the spin dynamic structure factor which quantitatively agrees with recent Brillouin light scattering experiments. We also show that at room temperature the phonon contribution to the magnon damping exhibits a rather complicated momentum dependence: In the exchange regime the magnon damping is dominated by Cherenkov type scattering processes, while in the long-wavelength dipolar regime these processes are subdominant and the magnon damping is two orders of magnitude smaller. We supplement our calculations by actual measurements of the magnon relaxation in the dipolar regime. Our theory provides a simple explanation of a recent experiment probing the different temperatures of the magnon and phonon gases in yttrium-iron garnet.

preprint2014arXiv

Singular spin-wave theory and scattering continua in the cone state of Cs_2CuCl_4

For temperatures below 0.6 K the geometrically frustrated layered quantum antiferromagnet Cs $_2$ CuCl $_4$ in a magnetic field perpendicular to the layers orders magnetically in a so-called cone state, where the magnetic moments have a finite component in the field direction while their projection onto the layers forms a spiral. Modeling this system by a two-dimensional spatially anisotropic quantum Heisenberg antiferromagnet with Dzyaloshinskii-Moriya interaction, we find that even for vanishing temperature the usual spin-wave expansion is plagued by infrared divergencies which are due to the coupling between longitudinal and transverse spin fluctuations in the cone state. Similar divergencies appear also in the ground state of the interacting Bose gas in two and three dimensions. Using known results for the correlation functions of the interacting Bose gas, we present a non-perturbative expression for the dynamic structure factor in the cone state of Cs $_2$ CuCl $_4$ . We show that in this state the spectral line shape of spin fluctuations exhibits singular scattering continua which can be understood in terms of the well-known anomalous longitudinal fluctuations in the ground state of the two-dimensional Bose gas.

preprint2013arXiv

Majorana spin liquid and dimensional reduction in Cs2CuCl4

The low-temperature behavior of the magnetic insulator Cs2CuCl4 can be modeled by an anisotropic triangular lattice spin-1/2 Heisenberg antiferromagnet with two different exchange couplings J and J' = J/3. We show that in a wide range of magnetic fields the experimentally observed field dependence of the crossover temperature T_c for spin-liquid behavior can be explained within a mean-field theory based on the representation of spin operators in terms of Majorana fermions. We also show that for small magnetic fields the specific heat and the spin susceptibility both exhibit a maximum as a function of temperature at T_c = J/2. In the spin-liquid regime, the Majorana fermions can only propagate along the direction of the strongest bond, in agreement with the dimensional reduction scenario advanced by Balents [Nature (London) 464, 199 (2010)].

preprint2013arXiv

Solution of the Anderson impurity model via the functional renormalization group

We show that the functional renormalization group is a numerically cheap method to obtain the low-energy behavior of the Anderson impurity model describing a localized interacting electron coupled to a bath of conduction electrons. Our approach uses an external magnetic field as flow parameter, partial bosonization of the transverse spin fluctuations, and frequency-independent interaction vertices which are fixed by Ward identities. We calculate the quasi-particle residue and the spin susceptibility in the particle-hole symmetric case and obtain excellent agreement with the Bethe ansatz results for arbitrary strengths of the interaction.

preprint2013arXiv

Spontaneous ferromagnetism in the spinor Bose gas with Rashba spin-orbit coupling

We show that in the two-component Bose gas with Rashba spin-orbit coupling an arbitrarily small attractive interaction between bosons with opposite spin induces spontaneous ferromagnetism below a finite critical temperature $T_c$ . In the ferromagnetic phase the single-particle spectrum exhibits a unique minimum in momentum space in the direction of the magnetization. For sufficiently small temperatures below $T_c$ the bosons eventually condense into the unique state at the bottom of the spectrum, forming a ferromagnetic Bose-Einstein condensate.

Damping of phase fluctuations in superfluid Bose gases

Using Popov's hydrodynamic approach we derive an effective Euclidean action for the long-wavelength phase fluctuations of superfluid Bose gases in D dimensions. We then use this action to calculate the damping of phase fluctuations at zero temperature as a function of D. For D >1 and wavevectors | k | << 2 mc (where m is the mass of the bosons and c is the sound velocity) we find that the damping in units of the phonon energy E_k = c | k | is to leading order gamma_k / E_k = A_D (k_0^D / 2 pi rho) (| k | / k_0)^{2 D -2}, where rho is the boson density and k_0 =2 mc is the inverse healing length. For D -> 1 the numerical coefficient A_D vanishes and the damping is proportional to an additional power of |k | /k_0; a self-consistent calculation yields in this case gamma_k / E_k = 1.32 (k_0 / 2 pi rho)^{1/2} |k | / k_0. In one dimension, we also calculate the entire spectral function of phase fluctuations.

Functional renormalization group approach to the Ising-nematic quantum critical point of two-dimensional metals

Using functional renormalization group methods, we study an effective low-energy model describing the Ising-nematic quantum critical point in two-dimensional metals. We treat both gapless fermionic and bosonic degrees of freedom on equal footing and explicitly calculate the momentum and frequency dependent effective interaction between the fermions mediated by the bosonic fluctuations. Following earlier work by S.-S. Lee for a one-patch model, Metlitski and Sachdev [Phys. Rev. B {\bf{82}}, 075127] recently found within a field-theoretical approach that certain three-loop diagrams strongly modify the one-loop results, and that the conventional 1/N expansion breaks down in this problem. We show that the singular three-loop diagrams considered by Metlitski and Sachdev are included in a rather simple truncation of the functional renormalization group flow equations for this model involving only irreducible vertices with two and three external legs. Our approximate solution of these flow equations explicitly yields the vertex corrections of this problem and allows us to calculate the anomalous dimension $η_ψ$ of the fermion field.

Thermalization of magnons in yttrium-iron garnet: nonequilibrium functional renormalization group approach

Using a nonequilibrium functional renormalization group (FRG) approach we calculate the time evolution of the momentum distribution of a magnon gas in contact with a thermal phonon bath. As a cutoff for the FRG procedure we use a hybridization parameter Λ giving rise to an artificial damping of the phonons. Within our truncation of the FRG flow equations the time evolution of the magnon distribution is obtained from a rate equation involving cutoff-dependent nonequilibrium self-energies, which in turn satisfy FRG flow equations depending on cutoff-dependent transition rates. Our approach goes beyond the Born collision approximation and takes the feedback of the magnons on the phonons into account. We use our method to calculate the thermalization of a quasi two-dimensional magnon gas in the magnetic insulator yttrium-iron garnet after a highly excited initial state has been generated by an external microwave field. We obtain good agreement with recent experiments.

Time-dependent spin-wave theory

We generalize the spin-wave expansion in powers of the inverse spin to time-dependent quantum spin models describing rotating magnets or magnets in time-dependent external fields. We show that in these cases, the spin operators should be projected onto properly defined rotating reference frames before the spin components are bosonized using the Holstein-Primakoff transformation. As a first application of our approach, we calculate the reorganization of the magnetic state due to Bose-Einstein condensation of magnons in the magnetic insulator yttrium-iron garnet; we predict a characteristic dip in the magnetization which should be measurable in experiments.

Bose-Einstein condensation at finite momentum and magnon condensation in thin film ferromagnets

We use the Gross-Pitaevskii equation to determine the spatial structure of the condensate density of interacting bosons whose energy dispersion epsilon_k has two degenerate minima at finite wave-vectors q. We show that in general the Fourier transform of the condensate density has finite amplitudes for all integer multiples of q. If the interaction is such that many Fourier components contribute, the Bose condensate is localized at the sites of a one-dimensional lattice with spacing 2pi/q; in this case Bose-Einstein condensation resembles the transition from a liquid to a crystalline solid. We use our results to investigate the spatial structure of the Bose condensate formed by magnons in thin films of ferromagnets with dipole-dipole interactions.

Elastic constants and ultrasonic attenuation in the cone state of the frustrated antiferromagnet Cs_2CuCl_4

In an external magnetic field perpendicular to the plane of the layers, the quasi two-dimensional frustrated antiferromagnet Cs_2CuCl_4 exhibits a magnetically ordered "cone state" at low temperatures. In this state the component of the magnetic moments in field direction is finite, while their projection onto the plane of the layers forms a spiral. We present both theoretical and experimental results for the magnetic field dependence of the elastic constants and the ultrasonic attenuation rate in the cone state. Our theoretical analysis is based on the usual spin-wave expansion around the classical ground state of a Heisenberg model on an anisotropic triangular lattice with Dzyaloshinskii-Moriya interactions. Magnon-phonon interactions are modeled by expanding the exchange interactions up to second order in powers of the phonon coordinates. As long as the external magnetic field is not too close to the critical field where the cone state becomes unstable, we obtain reasonable agreement between theory and experiment, suggesting that at least in this regime magnons behave as well-defined quasiparticles. We also show that the assumption of well-defined magnons implies that at long wavelengths the ultrasonic attenuation rate in the cone state of Cs_2CuCl_4 is proportional to the fourth power of the phonon momentum.

Non-equilibrium time evolution of bosons from the functional renormalization group

We develop a functional renormalization group approach to obtain the time evolution of the momentum distribution function of interacting bosons out of equilibrium. Using an external out-scattering rate as flow parameter, we derive formally exact renormalization group flow equations for the non-equilibrium self-energies in the Keldysh basis. A simple perturbative truncation of these flow equations leads to an approximate solution of the quantum Boltzmann equation which does not suffer from secular terms and gives accurate results even for long times. We demonstrate this explicitly within a simple exactly solvable toy model describing a quartic oscillator with off-diagonal pairing terms.

Quantum criticality of dipolar spin chains

We show that a chain of Heisenberg spins interacting with long-range dipolar forces in a magnetic field h perpendicular to the chain exhibits a quantum critical point belonging to the two-dimensional Ising universality class. Within linear spin-wave theory the magnon dispersion for small momenta k is [Delta^2 + v_k^2 k^2]^{1/2}, where Delta^2 \propto |h - h_c| and v_k^2 \propto |ln k|. For fields close to h_c linear spin-wave theory breaks down and we investigate the system using density-matrix and functional renormalization group methods. The Ginzburg regime where non-Gaussian fluctuations are important is found to be rather narrow on the ordered side of the transition, and very broad on the disordered side.

Functional renormalization group approach to interacting bosons at zero temperature

We investigate the single-particle spectral density of interacting bosons within the non-perturbative functional renormalization group technique. The flow equations for a Bose gas are derived in a scheme which treats the two-particle density-density correlations exactly but neglects irreducible correlations among three and more particles. These flow equations are solved within a truncation which allows to extract the complete frequency and momentum structure of the normal and anomalous self-energies. Both the asymptotic small momentum regime, where perturbation regime fails, as well as the perturbative regime at larger momenta are well described within a single unified approach. The self-energies do not exhibit any infrared divergences, satisfy the U(1) symmetry constraints, and are in accordance with the Nepomnyashchy relation which states that the anomalous self-energy vanishes at zero momentum and zero frequency. From the self-energies we extract the single-particle spectral density of the two-dimensional Bose gas. The dispersion is found to be of the Bogoliubov form and shows the crossover from linear Goldstone modes to the quadratic behavior of quasi-free bosons. The damping of the quasiparticles is found to be in accordance with the standard Beliaev damping. We furthermore recover the exact asymptotic limit of the propagators derived by Gavoret and Nozieres and discuss the nature of the non-analyticities of the self-energies in the very small momentum regime.

Parametric pumping and kinetics of magnons in dipolar ferromagnets

The time evolution of magnons subject to a time-dependent microwave field is usually described within the so-called "S-theory", where kinetic equations for the distribution function are obtained within the time-dependent Hartree-Fock approximation. To explain the recent observation of "Bose-Einstein condensation of magnons" in an external microwave field [Demokritov et al., Nature 443, 430 (2006)], we extend the "S-theory" to include the Gross-Pitaevskii equation for the time-dependent expectation values of the magnon creation and annihilation operators. We explicitly solve the resulting coupled equations within a simple approximation where only a single condensed mode is retained. We also re-examine the usual derivation of an effective boson model from a realistic spin model for yttrium-iron garnet films and argue that in the parallel pumping geometry (where both the static and the time-dependent magnetic field are parallel to the macroscopic magnetization) the time-dependent Zeemann energy cannot give rise to magnon condensation.

Spectral function of the Anderson impurity model at finite temperatures

Using the functional renormalization group (FRG) and the numerical renormalization group (NRG), we calculate the spectral function of the Anderson impurity model at zero and finite temperatures. In our FRG scheme spin fluctuations are treated non-perturbatively via a suitable Hubbard-Stratonovich field, but vertex corrections are neglected. A comparison with our highly accurate NRG results shows that this FRG scheme gives a quantitatively good description of the spectral line-shape at zero and finite temperatures both in the weak and strong coupling regimes, although at zero temperature the FRG is not able to reproduce the known exponential narrowing of the Kondo resonance at strong coupling.

Ward identities for the Anderson impurity model: derivation via functional methods and the exact renormalization group

Using functional methods and the exact renormalization group we derive Ward identities for the Anderson impurity model. In particular, we present a non-perturbative proof of the Yamada-Yosida identities relating certain coefficients in the low-energy expansion of the self-energy to thermodynamic particle number and spin susceptibilities of the impurity. Our proof underlines the relation of the Yamada-Yosida identities to the U(1) x U(1) symmetry associated with particle number and spin conservation in a magnetic field.

preprint2009arXiv

Condensate density of interacting bosons: a functional renormalization group approach

We calculate the temperature dependent condensate density $ρ^0 (T)$ of interacting bosons in three dimensions using the functional renormalization group (FRG). From the numerical solution of suitably truncated FRG flow equations for the irreducible vertices we obtain $ρ^0 (T)$ for arbitrary temperatures. We carefully extrapolate our numerical results to the critical point and determine the order parameter exponent $β\approx 0.32$ , in reasonable agreement with the expected value $0.345$ associated with the XY-universality class. We also calculate the condensate density in two dimensions at zero temperature using a truncation of the FRG flow equations based on the derivative expansion including cubic and quartic terms in the expansion of the effective potential in powers of the density. As compared with the widely used quadratic approximation for the effective potential, the coupling constants associated with the cubic and quartic terms increase the result for the condensate density by a few percent. However, the cubic and quartic coupling constants flow to rather large values, which sheds some doubt on FRG calculations based on a low order polynomial approximation for the effective potential.

preprint2009arXiv

Spectral function and quasi-particle damping of interacting bosons in two dimensions

We employ the functional renormalization group to study dynamical properties of the two-dimensional Bose gas. Our approach is free of infrared divergences, which plague the usual diagrammatic approaches, and is consistent with the exact Nepomnyashchy identity, which states that the anomalous self-energy vanishes at zero frequency and momentum. We recover the correct infrared behavior of the propagators and present explicit results for the spectral line-shape, from which we extract the quasi-particle dispersion and damping.

preprint2008arXiv

Functional renormalization group in the broken symmetry phase: momentum dependence and two-parameter scaling of the self-energy

We include spontaneous symmetry breaking into the functional renormalization group (RG) equations for the irreducible vertices of Ginzburg-Landau theories by augmenting these equations by a flow equation for the order parameter, which is determined from the requirement that at each RG step the vertex with one external leg vanishes identically. Using this strategy, we propose a simple truncation of the coupled RG flow equations for the vertices in the broken symmetry phase of the Ising universality class in D dimensions. Our truncation yields the full momentum dependence of the self-energy Sigma (k) and interpolates between lowest order perturbation theory at large momenta k and the critical scaling regime for small k. Close to the critical point, our method yields the self-energy in the scaling form Sigma (k) = k_c^2 sigma^{-} (k | xi, k / k_c), where xi is the order parameter correlation length, k_c is the Ginzburg scale, and sigma^{-} (x, y) is a dimensionless two-parameter scaling function for the broken symmetry phase which we explicitly calculate within our truncation.

preprint2007arXiv

Damping of zero sound in Luttinger liquids

We calculate the damping gamma_q of collective density oscillations (zero sound) in a one-dimensional Fermi gas with dimensionless forward scattering interaction F and quadratic energy dispersion k^2 / 2 m at zero temperature. For wave-vectors | q| /k_F small compared with F we find to leading order gamma_q = v_F^{-1} m^{-2} Y (F) | q |^3, where v_F is the Fermi velocity, k_F is the Fermi wave-vector, and Y (F) is proportional to F^3 for small F. We also show that zero-sound damping leads to a finite maximum proportional to |k - k_F |^{-2 + 2 eta} of the charge peak in the single-particle spectral function, where eta is the anomalous dimension. Our prediction agrees with photoemission data for the blue bronze K_{0.3}MoO_3.

preprint1998arXiv

Vertex corrections in gauge theories for two-dimensional condensed matter systems

We calculate the self-energy of two-dimensional fermions that are coupled to transverse gauge fields, taking two-loop corrections into account. Given a bare gauge field propagator that diverges for small momentum transfers q as 1 / q^{eta}, 1 < eta < 2, the fermionic self-energy without vertex corrections vanishes for small frequencies omega as Sigma (omega) propto omega^{gamma with gamma = {frac{2}{1 + eta}} < 1. We show that inclusion of the leading radiative correction to the fermion - gauge field vertex leads to Sigma (omega) propto omega^{gamma} [ 1 - a_{eta} ln (omega_0 / omega) ], where a_η is a positive numerical constant and omega_0 is some finite energy scale. The negative logarithmic correction is consistent with the scenario that higher order vertex corrections push the exponent gamma to larger values.

preprint1996arXiv

Bosonization and the eikonal expansion: similarities and differences

We compare two non-perturbative techniques for calculating the single-particle Green's function of interacting Fermi systems with dominant forward scattering: our recently developed functional integral approach to bosonization in arbitrary dimensions, and the eikonal expansion. In both methods the Green's function is first calculated for a fixed configuration of a background field, and then averaged with respect to a suitably defined effective action. We show that, after linearization of the energy dispersion at the Fermi surface, both methods yield for Fermi liquids exactly the same non-perturbative expression for the quasi-particle residue. However, in the case of non-Fermi liquid behavior the low-energy behavior of the Green's function predicted by the eikonal method can be erroneous. In particular, for the Tomonaga-Luttinger model the eikonal method neither reproduces the correct scaling behavior of the spectral function, nor predicts the correct location of its singularities.

Anomalous scaling and spin-charge separation in coupled chains

We use a bosonization approach to show that the three dimensional Coulomb interaction in coupled metallic chains leads to a Luttinger liquid for vanishing inter-chain hopping $t_{\bot}$ , and to a Fermi liquid for any finite $t_{\bot}$ . However, for small $t_{\bot} \neq 0$ the Greens-function satisfies a homogeneity relation with a non-trivial exponent $γ_{cb}$ in a large intermediate regime. Our results offer a simple explanation for the large values of $γ_{cb}$ inferred from recent photoemission data from quasi one-dimensional conductors and might have some relevance for the understanding of the unusual properties of the high-temperature superconductors.

Bosonization of interacting fermions in arbitrary dimension beyond the Gaussian approximation

We use our recently developed functional bosonization approach to bosonize interacting fermions in arbitrary dimension $d$ beyond the Gaussian approximation. Even in $d=1$ the finite curvature of the energy dispersion at the Fermi surface gives rise to interactions between the bosons. In higher dimensions scattering processes describing momentum transfer between different patches on the Fermi surface (around-the-corner processes) are an additional source for corrections to the Gaussian approximation. We derive an explicit expression for the leading correction to the bosonized Hamiltonian and the irreducible self-energy of the bosonic propagator that takes the finite curvature as well as around-the-corner processes into account. In the special case that around-the-corner scattering is negligible, we show that the self-energy correction to the Gaussian propagator is negligible if the dimensionless quantities $( \frac{q_{c} }{ k_{F}} )^d F_{0} [ 1 + F_{0} ]^{-1} \fracμ{ν^α} | \frac{ \partial ν^α }{ \partial μ} |$ are small compared with unity for all patches $α$ . Here $q_{c}$ is the cutoff of the interaction in wave-vector space, $k_{F}$ is the Fermi wave-vector, $μ$ is the chemical potential, $F_{0}$ is the usual dimensionless Landau interaction-parameter, and $ν^α$ is the {\it{local}} density of states associated with patch $α$ . We also show that the well known cancellation between vertex- and self-energy corrections in one-dimensional systems, which is responsible for the fact that the random-phase approximation for the density-density correlation function is exact in $d=1$ , exists also in $d> 1$ , provided (1) the interaction cutoff $q_{c}$ is small compared with $k_{F}$ , and (2) the energy dispersion is locally linearized at the Fermi the Fermi surface. Finally, we suggest a new systematic method to calculate corrections to the RPA, which is based on the perturbative calculation of the irreducible bosonic self-energy arising from the non-Gaussian terms of the bosonized Hamiltonian.

Functional Bosonization of Interacting Fermions in Arbitrary Dimensions

We bosonize the long-wavelength excitations of interacting fermions in arbitrary dimension by directly applying a suitable Hubbard-Stratonowich transformation to the Grassmannian generating functional of the fermionic correlation functions. With this technique we derive a surprisingly simple expression for the single-particle Greens-function, which is valid for arbitrary interaction strength and can describe Fermi- as well as Luttinger liquids. Our approach sheds further light on the relation between bosonization and the random-phase approximation, and enables us to study screening in a non-perturbative way.