Source author record

Carsten Honerkamp

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

37 published item(s)

preprint2022arXiv

Self-consistent formulations for stochastic nonlinear neuronal dynamics

Neural dynamics is often investigated with tools from bifurcation theory. However, many neuron models are stochastic, mimicking fluctuations in the input from unknown parts of the brain or the spiking nature of signals. Noise changes the dynamics with respect to the deterministic model; in particular bifurcation theory cannot be applied. We formulate stochastic neuronal dynamics in the Martin-Siggia-Rose de Dominicis-Janssen (MSRDJ) formalism and present the fluctuation expansion of the effective action and the functional renormalization group (fRG) as two systematic ways to incorporate corrections to the mean dynamics and time-dependent statistics due to fluctuations in the presence of nonlinear neuronal gain. To formulate self-consistency equations, we derive a fundamental link between the effective action in the Onsager-Machlup(OM) formalism, which allows the study of phase transitions, and the MSRDJ effective action, which is computationally advantageous. These results in particular allow the derivation of an OM effective action for systems with non-Gaussian noise. This approach naturally leads to effective deterministic equations for the first moment of the stochastic system; they explain how nonlinearities and noise cooperate to produce memory effects. Moreover, the MSRDJ formulation yields an effective linear system that has identical power spectra and linear response. Starting from the better known loopwise approximation, we then discuss the use of the fRG as a method to obtain self-consistency beyond the mean. We present a new efficient truncation scheme for the hierarchy of flow equations for the vertex functions by adapting the Blaizot, Méndez and Wschebor approximation from the derivative expansion to the vertex expansion. The methods are presented by means of the simplest possible example of a stochastic differential equation that has generic features of neuronal dynamics.

preprint2022arXiv

Understanding the three-dimensional quantum Hall effect in generic multi-Weyl semimetals

The quantum Hall effect in three-dimensional Weyl semimetal (WSM) receives significant attention for the emergence of the Fermi loop where the underlying two-dimensional Hall conductivity, namely, sheet Hall conductivity, shows quantized plateaus. Considering the tilted lattice models for multi Weyl semimetals (mWSMs), we systematically study the Landau levels (LLs) and magneto-Hall conductivity in the presence of parallel and perpendicular (with respect to the Weyl node's separation) magnetic field, i.e., $\mathbf{ B}\parallel z$ and $\mathbf{B}\parallel x$, to explore the impact of tilting and non-linearity in the dispersion. We make use of two (single) node low-energy models to qualitatively explain the emergence of mid-gap chiral (linear crossing of chiral) LLs on the lattice for $\mathbf{ B}\parallel z$ ($\mathbf{ B}\parallel x$). Remarkably, we find that the sheet Hall conductivity becomes quantized for $\mathbf{ B}\parallel z$ even when two Weyl nodes project onto a single Fermi point in two opposite surfaces, forming a Fermi loop with $k_z$ as the good quantum number. On the other hand, the Fermi loop, connecting two distinct Fermi points in two opposite surfaces, with $k_x$ being the good quantum number, causes the quantization in sheet Hall conductivity for $\mathbf{ B}\parallel x$. The quantization is almost lost (perfectly remained) in the type-II phase for $\mathbf{ B}\parallel x$ ($\mathbf{ B}\parallel z$). Interestingly, the jump profiles between the adjacent quantized plateaus change with the topological charge for both of the above cases. The momentum-integrated three-dimensional Hall conductivity is not quantized; however, it bears the signature of chiral LLs as resulting in the linear dependence on $μ$ for small $μ$. The linear zone (its slope) reduces (increases) as the tilt (topological charge) of the underlying WSM increases.

preprint2021arXiv

Spin-fluctuation-induced pairing in twisted bilayer graphene

We investigate the interplay of magnetic fluctuations and Cooper pairing in twisted bilayer graphene from a purely microscopic model within a large-scale tight-binding approach resolving the Ångström scale. For local onsite repulsive interactions and using the random-phase approximation for spin fluctuations, we derive a microscopic effective pairing interaction that we use for self-consistent solutions of the Bogoliubov-de-Gennes equations of superconductivity. We study the predominant pairing types as function of interaction strength, temperature and band filling. For large regions of this parameter space, we find chiral $d$-wave pairing regimes, spontaneously breaking time-reversal symmetry, separated by magnetic instabilities at integer band fillings. Interestingly, the $d$-wave pairing is strongly concentrated in the AA regions of the moiré unit cell and exhibits phase windings of integer multiples of $2π$ around these superconducting islands, i.e. pinned vortices. The spontaneous circulating current creates a distinctive magnetic field pattern. This signature of the chiral pairing should be measurable by state-of-the-art experimental techniques.

preprint2020arXiv

A Long-Range Ising Model of a Barabási-Albert Network

Networks that have power-law connectivity, commonly referred to as the scale-free networks, are an important class of complex networks. A heterogeneous mean-field approximation has been previously proposed for the Ising model of the Barabási-Albert model of scale-free networks with classical spins on the nodes wherein it was shown that the critical temperature for such a system scales logarithmically with network size. For finite sizes, there is no criticality for such a system and hence no true phase transition in terms of singular behavior. Further, in the thermodynamic limit, the mean-field prediction of an infinite critical temperature for the system may exclude any true phase transition even then. Nevertheless, with an eye on potential applications of the model on biological systems that are generally finite, one may still try to find approximations that describe the relevant observables quantitatively. Here we present an alternative, approximate formulation for the description of the Ising model of a Barabási-Albert Network. Using the classical definition of magnetization, we show that Ising models on a network can be well-approximated by a long-range interacting homogeneous Ising model wherein each node of the network couples to all other spins with a strength determined by the mean degree of the Barabási-Albert Network. In such an effective long-range Ising model of a Barabási-Albert Network, the critical temperature is directly proportional to the number of preferentially attached links added to grow the network. The proposed model describes the magnetization of the majority of the sites with average or smaller than average degree better compared to the heterogeneous mean-field approximation. The long-range Ising model is the only homogeneous description of Barabási-Albert networks that we know of.

preprint2020arXiv

Functional renormalization group for a large moiré unit cell

Layers of two-dimensional materials arranged at a twist angle with respect to each other lead to enlarged unit cells with potentially strongly altered band structures, offering a new arena for novel and engineered many-body ground states. For the exploration of these, renormalization group methods are an appropriate, flexible tool that take into account the mutual influence of competing tendencies. Here we show that, within reasonable, non-trivial approximations, the functional renormalization group known from simpler two-dimensional systems can be employed for the large-unit cell moiré superlattices with more than 10.000 bands, remedying the need to employ ad hoc restrictions to effective low-energy theories of a few bands and/or effective continuum theories. This provides a description on the atomic scale, allowing one to absorb available ab-initio information on the model parameters and therefore lending the analysis a more concrete quantitative character. For the case of twisted bilayer graphene models, we explore the leading ordering tendencies depending on the band filling and the range of interactions. The results indicate a delicate balance between distinct magnetically ordered ground states, as well as the occurrence of a charge modulation within the moiré unit cell for sufficiently non-local repulsive interaction.

preprint2020arXiv

Pseudogap opening in the two-dimensional Hubbard model: A functional renormalization group analysis

Using the recently introduced multiloop extension of the functional renormalization group, we compute the frequency- and momentum-dependent self-energy of the two-dimensional Hubbard model at half filling and weak coupling. We show that, in the truncated-unity approach for the vertex, it is essential to adopt the Schwinger-Dyson form of the self-energy flow equation in order to capture the pseudogap opening. We provide an analytic understanding of the key role played by the flow scheme in correctly accounting for the impact of the antiferromagnetic fluctuations. For the resulting pseudogap, we present a detailed numerical analysis of its evolution with temperature, interaction strength, and loop order.

preprint2020arXiv

Quantifying the fragility of unprotected quadratic band crossing points

We examine a basic lattice model of interacting fermions that exhibits quadratic band crossing points (QBCPs) in the non-interacting limit. In particular, we consider spinless fermions on the honeycomb lattice with nearest neighbor hopping $t$ and third-nearest neighbor hopping $t''$, which exhibits fine-tuned QBCPs at the corners of the Brillouin zone for ${t'' = t/2}$. In this situation, the density of states remains finite at the Fermi level of the half-filled band and repulsive nearest-neighbor interactions $V$ lead to a charge-density-wave (CDW) instability at infinitesimally small $V$ in the random-phase approximation or mean-field theory. We examine the fragility of the QBCPs against dispersion renormalizations in the ${t\mbox{-}t''\mbox{-}V}$ model using perturbation theory, and find that the $t''$-value needed for the QBCPs increases with $V$ due to the hopping renormalization. However, the instability toward CDW formation always requires a nonzero threshold interaction strength, i.e., one cannot fine-tune $t''$ to recover the QBCPs in the interacting system. These perturbative arguments are supported by quantum Monte Carlo simulations for which we carefully compare the corresponding threshold scales at and beyond the QBCP fine-tuning point. From this analysis, we thus gain a quantitative microscopic understanding of the fragility of the QBCPs in this basic interacting fermion system.

preprint2020arXiv

Quantitative functional renormalization-group description of the two-dimensional Hubbard model

Using a leading algorithmic implementation of the functional renormalization group (fRG) for interacting fermions on two-dimensional lattices, we provide a detailed analysis of its quantitative reliability for the Hubbard model. In particular, we show that the recently introduced multiloop extension of the fRG flow equations for the self-energy and two-particle vertex allows for a precise match with the parquet approximation also for two-dimensional lattice problems. The refinement with respect to previous fRG-based computation schemes relies on an accurate treatment of the frequency and momentum dependences of the two-particle vertex, which combines a proper inclusion of the high-frequency asymptotics with the so-called truncated unity fRG for the momentum dependence. The adoption of the latter scheme requires, as an essential step, a consistent modification of the flow equation of the self-energy. We quantitatively compare our fRG results for the self-energy and momentum-dependent susceptibilities and the corresponding solution of the parquet approximation to determinant quantum Monte Carlo data, demonstrating that the fRG is remarkably accurate up to moderate interaction strengths. The presented methodological improvements illustrate how fRG flows can be brought to a quantitative level for two-dimensional problems, providing a solid basis for the application to more general systems.

preprint2020arXiv

Truncated Unity Parquet Solver

We present an implementation of a truncated unity parquet solver (TUPS) which solves the parquet equations using a truncated form-factor basis for the fermionic momenta. This way fluctuations from different scattering channels are treated on an equal footing. The essentially linear scaling of computational costs in the number of untruncated bosonic momenta allows us to treat system sizes of up to 76x76 discrete lattice momenta, unprecedented by previous unbiased methods that include the frequency dependence of the vertex. With TUPS, we provide the first numerical evidence that the parquet approximation might indeed respect the Mermin-Wagner theorem and further systematically analyze the convergence with respect to the number of form factors. Using a single form factor seems to qualitatively describe the physics of the half-filled Hubbard model correctly, including the pseudogap behaviour. Quantitatively, using a single or a few form factors only is not sufficient at lower temperatures or stronger coupling.

preprint2015arXiv

Correlated spinless fermions on the honeycomb lattice revisited

We investigate the quantum many-body instabilities of the extended Hubbard model for spinless fermions on the honeycomb lattice with repulsive nearest-neighbor and 2nd nearest-neighbor density-density interactions. Recent exact diagonalization and infinite density matrix renormalization group results suggest that a putative topological Mott insulator phase driven by the 2nd nearest-neighbor repulsion is suppressed, while other numerically exact approaches support the topological Mott insulator scenario. In the present work, we employ the functional renormalization group (fRG) for correlated fermionic systems. Our fRG results hint at a strong suppression of the scattering processes stabilizing the topological Mott insulator. From analyzing the effects of fermionic fluctuations, we obtain a phase diagram which is the result of the competition of various charge ordering instabilities.

preprint2015arXiv

Low-energy effective interactions beyond the constrained random-phase approximation by the functional renormalization group

In the derivation of low-energy effective models for solids targeting the bands near the Fermi level, the constrained random phase approximation (cRPA) has become an appreciated tool to compute the effective interactions. The Wick-ordered constrained functional renormalization group (cfRG) generalizes the cRPA approach by including all interaction channels in an unbiased way. Here we present applications of the cfRG to two simple multi-band systems and compare the resulting effective interactions to the cRPA. First we consider a multiband model for monolayer graphene, where we integrate out the $σ$-bands to get an effective theory for $π$-bands. It turns out that terms beyond cRPA are strongly suppressed by the different $xy$-plane reflection symmetry of the bands. In our model the cfRG-corrections to cRPA become visible when one disturbs this symmetry difference slightly, however without qualitative changes. This study shows that the embedding or layering of two-dimensional electronic systems can alter the effective interaction parameters beyond what is expected from screening considerations. The second example is a one-dimensional model for a diatomic system reminiscent of a CuO chain, where we consider an effective theory for Cu 3d-like orbitals. Here the fRG data shows relevant and qualitative corrections compared to the cRPA results. We argue that the new interaction terms affect the magnetic properties of the low-energy model.

preprint2015arXiv

Triplet pairing driven by Hund's coupling in doped monolayer MoS$_2$

We investigate superconducting pairing driven by electron-electron interactions in a theoretical model for monolayer MoS$_2$ with the temperature-flow functional renormalization group(fRG). At low doping, the dominant instability is toward odd-parity pairing with $f$-wave Mo-nearest-neighbor structure. We compute the fRG phase diagram versus electron doping below the van Hove filling of the conduction band. In the superconducting regime, the critical temperature grows with doping, comparable to the experiments. Near van Hove filling the system favors a ferromagnetic state. We demonstrate that the triplet pairing is driven by ferromagnetic fluctuations and that the multiorbital nature of the conduction band as well as the Hund's coupling appear crucial in making the physics of MoS$_2$ different from e.g. doped graphene.

preprint2014arXiv

Electronic Instabilities of the AA-Honeycomb Bilayer

We use a functional renormalization group approach to study the instabilities due to electron-electron interactions in a bilayer honeycomb lattice model with AA stacking, as it might be relevant for layered graphene with this structure. Starting with a tight- binding description for the four $π$-bands, we integrate out the modes of the dispersion by successively lowering an infrared cutoff and determine the leading tendencies in the effective interactions. The antiferromagnetic spin-density wave is an expected instability for dominant local repulsion among the electrons, but for nonlocal interaction terms also other instabilities occur. We discuss the phase diagrams depending on the model parameters. We compare our results to single-layer graphene and the more common AB-stacked bilayer, both qualitatively and quantitatively.

preprint2014arXiv

Functional renormalization group for commensurate antiferromagnets: Beyond the mean-field picture

We present a functional renormalization group (fRG) formalism for interacting fermions on lattices that captures the flow into states with commensurate spin-density wave order. During the flow, the growth of the order parameter is fed back into the flow of the interactions and all modes can be integrated out. This extends previous fRG flows in the symmetric phase that run into a divergence at a nonzero RG scale, i.e., that have to be stopped at the ordering scale. We use the corresponding Ward identity to check the accuracy of the results. We apply our new method to a model with two Fermi pockets that have perfect particle-hole nesting. The results obtained from the fRG are compared with those in random phase approximation.

preprint2014arXiv

Instabilities on graphene's honeycomb lattice with electron-phonon interactions

We study the impact of electron-phonon interactions on the many-body instabilities of electrons on the honeycomb lattice and their interplay with repulsive local and non-local Coulomb interactions at charge neutrality. To that end, we consider in-plane optical phonon modes with wavevectors close to the $Γ$ point as well as to the $K, -K$ points and calculate the effective phonon-mediated electron-electron interaction by integrating out the phonon modes. Ordering tendencies are studied by means of a momentum-resolved functional renormalization group approach allowing for an unbiased investigation of the appearing instabilities. In the case of an exclusive and supercritical phonon-mediated interaction, we find a Kekulé and a nematic bond ordering tendency being favored over the $s$-wave superconducting state. The competition between the different phonon-induced orderings clearly shows a repulsive interaction between phonons at small and large wavevector transfers. We further discuss the influence of phonon-mediated interactions on electronically-driven instabilities induced by onsite, nearest neighbor and next-to-nearest neighbor density-density interactions. We find an extension of the parameter regime of the spin density wave order going along with an increase of the critical scales where ordering occurs, and a suppression of competing orders.

preprint2014arXiv

Unconventional pairing and electronic dimerization instabilities in the doped Kitaev-Heisenberg model

We study the quantum many-body instabilities of the $t -J_{\mathrm{K}} - J_{\mathrm{H}}$ Kitaev-Heisenberg Hamiltonian on the honeycomb lattice as a minimal model for a doped spin-orbit Mott insulator. This spin-$1/2$ model is believed to describe the magnetic properties of the layered transition-metal oxide Na$_2$IrO$_3$. We determine the ground-state of the system with finite charge-carrier density from the functional renormalization group (fRG) for correlated fermionic systems. To this end, we derive fRG flow-equations adapted to the lack of full spin-rotational invariance in the fermionic interactions, here represented by the highly frustrated and anisotropic Kitaev exchange term. Additionally employing a set of Ward identities for the Kitaev-Heisenberg model, the numerical solution of the flow equations suggests a rich phase diagram emerging upon doping charge carriers into the ground-state manifold ($\mathbb{Z}_2$ quantum spin liquids and magnetically ordered phases). We corroborate superconducting triplet $p$-wave instabilities driven by ferromagnetic exchange and various singlet pairing phases. For filling $δ> 1/4$, the $p$-wave pairing gives rise to a topological state with protected Majorana edge-modes. For antiferromagnetic Kitaev and ferromagnetic Heisenberg exchange we obtain bond-order instabilities at van Hove filling supported by nesting and density-of-states enhancement, yielding dimerization patterns of the electronic degrees of freedom on the honeycomb lattice. Further, our flow equations are applicable to a wider class of model Hamiltonians.

preprint2013arXiv

An alternative functional renormalization group approach to the single impurity Anderson model

We present an alternative functional renormalization group (fRG) approach to the single-impurity Anderson model at finite temperatures. Starting with the exact self-energy and interaction vertex of a small system ('core') containing a correlated site, we switch on the hybridization with a non-interacting bath in the fRG-flow and calculate spectra of the correlated site. Different truncations of the RG-flow-equations and choices of the core are compared and discussed. Furthermore we calculate the linear conductance and the magnetic susceptibility as functions of temperature and interaction strength. The signatures of Kondo physics arising in the flow are compared with numerical renormalization group results.

preprint2013arXiv

Interplay between Point-Group Symmetries and the Choice of the Bloch Basis in Multiband Models

We analyze the point-group symmetries of generic multiband tight-binding models with respect to the transformation properties of the effective interactions. While the vertex functions in the orbital language may transform non-trivially under point-group operations, their point-group behavior in the band language can be simplified by choosing a suitable Bloch basis.We first give two analytically accessible examples. Then we show that, for a large class of models, a natural Bloch basis exists, in which the vertex functions in the band language transform trivially under all point-group operations. As a consequence, the point-group symmetries can be used to reduce the computational effort in perturbative many-particle approaches such as the functional renormalization group.

preprint2013arXiv

Multiorbital effects in the functional renormalization group: A weak-coupling study of the Emery model

We perform an instability analysis of the Emery three-band model at hole doping and weak coupling within a channel-decomposed functional renormalization group flow proposed in Phys. Rev. B 79, 195125 (2009). In our approach, momentum dependences are taken into account with improved precision compared to previous studies of related models. Around a generic parameter set, we find a strong competition of antiferromagnetic and d-wave Cooper instabilities with a smooth behavior under a variation of doping and additional hopping parameters. For increasingly incommensurate ordering tendencies in the magnetic channel, the d-wave pairing gap is deformed at its maxima. Comparing our results for the Emery model to those obtained for the two-dimensional one-band Hubbard model with effective parameters, we find that, despite considerable qualitative agreement, multi-orbital effects have a significant impact on a quantitative level.

preprint2013arXiv

Two-particle-correlations in a functional renormalization group scheme using a dynamical mean-field theory approach

We apply a recently introduced hybridization-flow functional renormalization group scheme for Anderson-like impurity models as an impurity solver in a dynamical mean-field theory (DMFT) approach to lattice Hubbard models. We present how this scheme is capable of reproducing metallic and insulating solutions of the lattice model. Our setup also offers a numerically rather inexpensive method to calculate two-particle correlation functions. For the paramagnetic Hubbard-model on the Bethe lattice in infinite dimensions we calculate the local two-particle-vertex for the metallic and the insulating phase. Then we go to a two-site cluster-DMFT-scheme for the two-dimensional Hubbard-model that includes short-range antiferromagnetic fluctuations and obtain the local and non-local two-particle-vertex-functions. We discuss the rich frequency structures of these vertices and compare with the vertex in the single-site solution.

preprint2012arXiv

Antiferromagnetism in the Hubbard Model on the Bernal-stacked Honeycomb Bilayer

Using a combination of quantum Monte Carlo simulations, functional renormalization group calculations and mean-field theory, we study the Hubbard model on the Bernal-stacked honeycomb bilayer at half-filling as a model system for bilayer graphene. The free bands consisting of two Fermi points with quadratic dispersions lead to a finite density of states at the Fermi level, which triggers an antiferromagnetic instability that spontaneously breaks sublattice and spin rotational symmetry once local Coulomb repulsions are introduced. Our results reveal an inhomogeneous participation of the spin moments in the ordered ground state, with enhanced moments at the three-fold coordinated sites. Furthermore, we find the antiferromagnetic ground state to be robust with respect to enhanced interlayer couplings and extended Coulomb interactions.

preprint2012arXiv

Critical scales in anisotropic spin systems from functional renormalization

We apply a recently developed functional renormalization group (fRG) scheme for quantum spin systems to the spin-1/2 antiferromagnetic XXZ model on a two-dimensional square lattice. Based on an auxiliary fermion representation we derive flow equations which allow a resummation of the perturbation series in the spin-spin interactions. Spin susceptibilities are calculated for different values of the anisotropy parameter. The phase transition between planar and axial ordering at the isotropic point is reproduced correctly. The results for the critical scales from the fRG as quantitative measures for the ordering temperatures are in good agreement with the exact solution in the Ising limit. On the easy-plane side, the deviations from critical temperatures obtained with quantum Monte Carlo are larger but still acceptable. However, at the isotropic point the Mermin-Wagner theorem is violated such that a precise description of the behavior in the vicinity of the phase transition is not possible. We discuss possible reasons for these discrepanies.

preprint2012arXiv

Effective three-particle interactions in low-energy models for multiband systems

We discuss different approximations for effective low-energy interactions in multi-band models for weakly correlated electrons. In the study of Fermi surface instabilities of the conduction band(s), the standard approximation consists only keeping those terms in the bare interactions that couple only to the conduction band(s), while corrections due to virtual excitations into bands away from the Fermi surface are typically neglected. Here, using a functional renormalization group approach, we present an improved truncation for the treatment of the effective interactions in the conduction band that keeps track of the generated three-particle interactions (six-point term) and hence allows one to include important aspects of these virtual interband excitations. Within a simplified two-patch treatment of the conduction band, we demonstrate that these corrections can have a rather strong effect in parts of the phase diagram by changing the critical scales for various orderings and the phase boundaries.

preprint2012arXiv

Functional renormalization group approach to correlated fermion systems

Numerous correlated electron systems exhibit a strongly scale-dependent behavior. Upon lowering the energy scale, collective phenomena, bound states, and new effective degrees of freedom emerge. Typical examples include (i) competing magnetic, charge, and pairing instabilities in two-dimensional electron systems, (ii) the interplay of electronic excitations and order parameter fluctuations near thermal and quantum phase transitions in metals, (iii) correlation effects such as Luttinger liquid behavior and the Kondo effect showing up in linear and non-equilibrium transport through quantum wires and quantum dots. The functional renormalization group is a flexible and unbiased tool for dealing with such scale-dependent behavior. Its starting point is an exact functional flow equation, which yields the gradual evolution from a microscopic model action to the final effective action as a function of a continuously decreasing energy scale. Expanding in powers of the fields one obtains an exact hierarchy of flow equations for vertex functions. Truncations of this hierarchy have led to powerful new approximation schemes. This review is a comprehensive introduction to the functional renormalization group method for interacting Fermi systems. We present a self-contained derivation of the exact flow equations and describe frequently used truncation schemes. Reviewing selected applications we then show how approximations based on the functional renormalization group can be fruitfully used to improve our understanding of correlated fermion systems.

preprint2012arXiv

Instabilities of quadratic band crossing points

Using a functional renormalization group approach, we study interaction-driven instabilities in quadratic band crossing point two-orbital models in two dimensions, extending a previous study of Sun et al. [1]. The wavevector-dependence of the Bloch eigenvectors of the free Hamiltonian causes interesting instabilities toward spin nematic, quantum anomalous Hall and quantum spin Hall states. In contrast with other known examples of interaction-driven topological insulators, in the system studied here, the QSH state occurs at arbitrarily small interaction strength and for rather simple intra- and inter-orbital repulsions.

preprint2012arXiv

Interacting electrons on trilayer honeycomb lattices

Few-layer graphene systems come in various stacking orders. Considering tight-binding models for electrons on stacked honeycomb layers, this gives rise to a variety of low-energy band structures near the charge neutrality point. Depending on the stacking order these band structures enhance or reduce the role of electron-electron interactions. Here, we investigate the instabilities of interacting electrons on honeycomb multilayers with a focus on trilayers with ABA and ABC stackings theoretically by means of the functional renormalization group. We find different types of competing instabilities and identify the leading ordering tendencies in the different regions of the phase diagram for a range of local and non-local short-ranged interactions. The dominant instabilities turn out to be toward an antiferromagnetic spin-density wave (SDW), a charge density wave and toward quantum spin Hall (QSH) order. Ab-initio values for the interaction parameters put the systems at the border between SDW and QSH regimes. Furthermore, we discuss the energy scales for the interaction-induced gaps of this model study and put them into context with the scales for single-layer and Bernal-stacked bilayer honeycomb lattices. This yields a comprehensive picture of the possible interaction-induced ground states of few-layer graphene.

preprint2012arXiv

Multiband Effects on Superconducting Instabilities Driven by Electron-Electron Interactions

We explore multiband effects on d-wave superconducting instabilities driven by electron-electron interactions. Our models on the two-dimensional square lattice consist of a main band with an extended Fermi surface and predominant weight from $d_{x^2-y^2}$ orbitals, whose orbital character is influenced by the admixture of other energetically neighbored orbitals. Using a functional renormalization group description of the superconducting instabilities of the system and different levels of approximations, we study how the energy scale for pairing and hence the critical temperature is affected by the band structure. We find that a reduction of orbital admixture as a function of the orbital energies can cause a $T_c$ enhancement although the Fermi surface becomes more curved and hence less favorable for antiferromagnetic spin fluctuations. While our study does not allow a quantitative understanding of the $T_c$ differences in realistic high-$T_c$ cuprate systems, it may reveal an underlying mechanism contributing to the actual material trends.

preprint2012arXiv

Renormalization group flow for fermions into antiferromagnetically ordered phases: Method and mean-field models

We present a functional renormalization group flow for many-fermion lattice models into phases with broken spin-rotational symmetry. The flow is expressed purely in terms of fermionic vertex functions. The symmetry breaking is seeded by a small initial anomalous self-energy which grows at the transition scale and which prevents a runaway-flow at nonzero scales. Focusing on the case of commensurate antiferromagnetism we discuss how the interaction vertex can be parametrized efficiently. For reduced models with long-range bare interactions we show the results of standard mean-field theory are reproduced by the fRG and how anisotropies in the spin sector change the flows. We then describe a more efficient decomposition of the interaction vertex that should allow for the treatment of more general models.

preprint2012arXiv

Self-energy feedback and frequency-dependent interactions in the functional renormalization group flow for the two-dimensional Hubbard model

We study the impact of including the self-energy feedback and frequency-dependent interactions on functional renormalization group grows for the two-dimensional Hubbard model on the square lattice at weak to moderate coupling strength. Previous studies using the functional renormalization group had ignored these two ingredients to large extent, and the question is how much the flows to strong coupling analyzed by this method depend on these approximations. Here we include the imaginary part of the self-energy on the imaginary axis and the frequency-dependence of the running interactions on a frequency mesh of 10 frequencies on the Matsubara axis. We find that i) the critical scales for the flows to strong coupling are shifted downwards by a factor that is usually of order one but can get larger in specific parameter regions, and ii) that the leading channel in this flow does not depend strongly on whether self-energies and frequency-dependence is included or not. We also discuss the main features of the self-energies developing during the flows.

preprint2011arXiv

Effective interactions in multi-band systems from constrained summations

In correlated electron materials, the application of many-body techniques for the study of interaction effects or unconventional superconductivity often requires the formulation of an effective low-energy model that contains only the relevant bands near the Fermi level. However the bands away from the Fermi level are known to renormalize the low-energy interactions substantially. Here we compare different schemes to derive low-energy effective theories for interacting electrons in solids. The frequently used constrained random phase approximation (cRPA) is identified as a particular resummation of higher-order interaction terms that includes important virtual corrections. We then propose an adapted functional renormalization group (fRG) scheme that includes the cPRA, but also allows one to go beyond the cRPA approximation. We study a simple two-band model in order to explore the differences between the approximations.

preprint2011arXiv

Instabilities of interacting electrons on the honeycomb bilayer

We investigate the instabilities of interacting electrons on the honeycomb bilayer by means of the functional renormalization group for a range of interactions up to the third-nearest neighbor. Besides a novel instability toward a gapless charge-density wave we find that using interaction parameters as determined by ab-initio calculations for graphene and graphite puts the system close to the boundary between antiferromagnetic and quantum spin Hall instabilities. Importantly, the energy scales for these instabilities are large such that imperfections and deviations from the basic model are expected to play a major role in real bilayer graphene, where interaction effects seem to be seen only at smaller scales. We therefore analyze how reducing the critical scale and small doping of the layers affect the instabilities.

preprint2011arXiv

Mechanism for a Pairing State with Time-Reversal Symmetry Breaking in Iron-Based Superconductors

The multipocket Fermi surfaces of iron-based superconductors promote pairing states with both s_{+-}-wave and d_{x^2-y^2}-wave symmetry. We argue that the competition between these two order parameters could lead to a time-reversal-symmetry breaking state with s+id-pairing symmetry in the iron-based superconductors, and propose serveral scenarios in which this phase may be found. To understand the emergence of such a pairing state on a more rigorous footing, we start from a microscopic 5-orbital description representative for the pnictides. Using a combined approach of functional renormalization group and mean-field analysis, we identify the microscopic parameters of the s+id-pairing state. There, we find the most promising region for s+id-pairing in the electron doped regime with an enhanced pnictogen height.

preprint2010arXiv

Effective low-energy Hamiltonians for interacting nanostructures

We present a functional renormalization group (fRG) treatment of trigonal graphene nanodiscs and composites thereof, modeled by finite-size Hubbard-like Hamiltonians with honeycomb lattice structure. At half filling, the noninteracting spectrum of these structures contains a certain number of half-filled states at the Fermi level. For the case of trigonal nanodiscs, including interactions between these degenerate states was argued to lead to a large ground state spin with potential spintronics applications. Here we perform a systematic fRG flow where the excited single-particle states are integrated out with a decreasing energy cutoff, yielding a renormalized low-energy Hamiltonian for the zero-energy states that includes effects of the excited levels. The numerical implementation corroborates the results obtained with a simpler Hartree-Fock treatment of the interaction effects within the zero-energy states only. In particular, for trigonal nanodiscs the degeneracy of the one-particle-states with zero-energy turns out to be very robust against influences of the higher levels. As an explanation, we give a general argument that within this fRG scheme the zero-energy degeneracy remains unsplit under quite general conditions and for any size of the trigonal nanodisc. We furthermore discuss the differences in the effective Hamiltonian and their ground states of single nanodiscs and composite bow-tie-shaped systems.

preprint2010arXiv

Exact diagonalization study of the trionic crossover and the trion liquid in the attractive three-component Hubbard model

We investigate the trion formation and the effective trionic properties in the attractive Hubbard model with three fermionic colors using exact diagonalization. The crossover to the trionic regime with colorless compound fermions upon increasing strength of the onsite attraction parameter U features smoothly evolving ground state properties and exhibits clear similarities to the BCS/BEC-crossover for two colors. In the excitation spectrum, there is a clear gap opening between a band of well-defined trions and excitations of broken-up trions at U_c ~ 1.8t. This picture remains the same away from the SU(3)-symmetric point. The spatial pairing correlations for colored Cooper pairs are compatible with a power-law at small attractions and change to an exponential decay above the trionic crossover. Furthermore, we show that the effective trionic liquid for U > U_c can be well modeled with spinless 'heavy' fermions interacting with a strong nearest neighbor repulsion.

preprint2009arXiv

Pairing in the iron arsenides: a functional RG treatment

We study the phase diagram of a microscopic model for the superconducting iron arsenides by means of a functional renormalization group. Our treatment establishes a connection between a strongly simplified two-patch model by Chubukov et al. and a five-band- analysis by Wang et al.. For a wide parameter range, the dominant pairing instability occurs in the extended s-wave channel. The results clearly show the relevance of pair scattering between electron and hole pockets. We also give arguments that the phase transition between the antiferromagnetic phase for the undoped system and the superconducting phase may be first order.

preprint2006arXiv

Phonons and d-wave pairing in the two-dimensional Hubbard model

We analyze the influence of phonons on the d-wave pairing instability in the Hubbard model on the two-dimensional square lattice at weak to moderate interaction U, using a functional renormalization group scheme with frequency-dependent interaction vertices. As measured by the pairing scale, the B1g buckling mode enhances the pairing, while other phonon modes decrease the pairing. When various phonon modes are included together, the net effect on the scale is small. However, in situations where d-wave superconductivity and other tendencies, e.g. antiferromagnetism, are closely competing, the combined effect of different phonons may be able to tip the balance towards pairing.

preprint1997arXiv

Andreev reflection in unitary and non-unitary triplet states

The quasiparticle reflection and transmission properties at normal conductor-superconductor interfaces are discussed for unitary and non-unitary spin triplet pairing states recently proposed for Sr_2RuO_4. We find resonance peaks in the Andreev reflection amplitude, which are related to surface bound states in the superconductor. They lead to conductance peak features below the quasiparticle gap in the superconductor. The position of the peaks depends on the angle of the incident quasiparticles in a different way for unitary and non-unitary states. Based on this observation we propose a possible experiment which allows to distinguish between the two kinds of superconducting states.