Hierarchical Parallelisation of Functional Renormalisation Group Calculations -- hp-fRG

The functional renormalisation group (fRG) has evolved into a versatile tool in condensed matter theory for studying important aspects of correlated electron systems. Practical applications of the method often involve a high numerical effort, motivating the question in how far High Performance Computing (HPC) can leverage the approach. In this work we report on a multi-level parallelisation of the underlying computational machinery and show that this can speed up the code by several orders of magnitude. This in turn can extend the applicability of the method to otherwise inaccessible cases. We exploit three levels of parallelisation: Distributed computing by means of Message Passing (MPI), shared-memory computing using OpenMP, and vectorisation by means of SIMD units (single-instruction-multiple-data). Results are provided for two distinct High Performance Computing (HPC) platforms, namely the IBM-based BlueGene/Q system JUQUEEN and an Intel Sandy-Bridge-based development cluster. We discuss how certain issues and obstacles were overcome in the course of adapting the code. Most importantly, we conclude that this vast improvement can actually be accomplished by introducing only moderate changes to the code, such that this strategy may serve as a guideline for other researcher to likewise improve the efficiency of their codes

Fermionic renormalization group flow into phases with broken discrete symmetry: charge-density wave mean-field model

We generalize the application of the functional renormalization group (fRG) method for the fermionic flow into the symmetry-broken phase to finite temperatures. We apply the scheme to the case of a broken discrete symmetry: the charge-density wave (CDW) mean-field model at half filling. We show how an arbitrarily small initial CDW order parameter starts to grow at the CDW instability and how it flows to the correct final value, suppressing the divergence of the effective interaction in the fRG flow. The effective interaction peaks at the instability and saturates at low energy scales or temperatures. The relation to the mean-field treatment, differences compared to the flow for a broken continuous symmetry, and the prospects of the new method are discussed.Comment: 10 pages, 7 figures. V2: minor corrections, journal references and DOI adde

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.Comment: 15 pages, 15 figures, version as publishe

Pseudogap at hot spots in the two-dimensional Hubbard model at weak coupling

We analyze the interaction-induced renormalization of single-particle excitations in the two-dimensional Hubbard model at weak coupling using the Wick-ordered version of the functional renormalization group. The self energy is computed for real frequencies by integrating a flow equation with renormalized two-particle interactions. In the vicinity of hot spots, that is points where the Fermi surface intersects the umklapp surface, self energy effects beyond the usual quasi-particle renormalizations and damping occur near instabilities of the normal, metallic phase. Strongly enhanced renormalized interactions between particles at different hot spots generate a pronounced low-energy peak in the imaginary part of the self energy, leading to a pseudogap-like double-peak structure in the spectral function for single-particle excitations.Comment: 14 pages, 7 figure

Interaction flow method for many-fermion systems

We propose an interaction flow scheme that sums up the perturbation expansion of many-particle systems by successively increasing the interaction strength. It combines the unbiasedness of renormalization group methods with the simplicity of straight-forward perturbation theory. Applying the scheme to fermions in one dimension and to the two-dimensional Hubbard model we find that at one-loop level and low temperatures there is ample agreement with previous one-loop renormalization group approaches. We furthermore present results for the momentum-dependence of spin, charge and pairing interactions in the two-dimensional Hubbard model.Comment: 14 pages, 14 figure

Quasi-particle functional Renormalisation Group calculations in the two-dimensional t-t'-Hubbard model

We extend and apply a recently introduced quasi-particle functional renormalisation group scheme to the two-dimensional Hubbard model with next-nearest-neighbour hopping and away from half filling. We confirm the generation of superconducting correlations in some regions of the phase diagram, but also find that the inclusion of self-energy feedback by means of a decreasing quasi-particle weight can suppress superconducting tendencies more than anti-ferromagnetic correlations by which they are generated. As a supplement, we provide sample results for the self-energy in second-order perturbation theory and address some conceptual matters

Pair fluctuation induced pseudogap in the normal phase of the two-dimensional attractive Hubbard model at weak coupling

One-particle spectral properties in the normal phase of the two-dimensional attractive Hubbard model are investigated in the weak coupling regime using the non-selfconsistent T-matrix approximation. The corresponding equations are evaluated numerically directly on the real frequency axis. For temperatures sufficiently close to the superconducting transition temperature a pseudogap in the one-particle spectral function is observed, which can be assigned to the increasing importance of pair fluctuations.Comment: 22 pages, 13 figure

Scaling study of the pion electroproduction cross sections and the pion form factor

The $^{1}$H($e,e^\prime \pi^+$)n cross section was measured for a range of four-momentum transfer up to $Q^2$=3.91 GeV$^2$ at values of the invariant mass, $W$, above the resonance region. The $Q^2$-dependence of the longitudinal component is consistent with the $Q^2$-scaling prediction for hard exclusive processes. This suggests that perturbative QCD concepts are applicable at rather low values of $Q^2$. Pion form factor results, while consistent with the $Q^2$-scaling prediction, are inconsistent in magnitude with perturbative QCD calculations. The extraction of Generalized Parton Distributions from hard exclusive processes assumes the dominance of the longitudinal term. However, transverse contributions to the cross section are still significant at $Q^2$=3.91 GeV$^2$.Comment: 6 pages, 3 figure

Measurement of Nuclear Transparency for the A(e,e' pi^+) Reaction

We have measured the nuclear transparency of the A(e,e' pi^+) process in ^{2}H,^{12}C, ^{27}Al, ^{63}Cu and ^{197}Au targets. These measurements were performed at the Jefferson Laboratory over a four momentum transfer squared range Q^2 = 1.1 - 4.7 (GeV/c)^2. The nuclear transparency was extracted as the super-ratio of $(\sigma_A/\sigma_H)$ from data to a model of pion-electroproduction from nuclei without pi-N final state interactions. The Q^2 and atomic number dependence of the nuclear transparency both show deviations from traditional nuclear physics expectations, and are consistent with calculations that include the quantum chromodynamical phenomenon of color transparency.Comment: 5 pages, 3 figs Changes to figure 2 and 3 (error band updated and theory curves updated