An introduction to numerical methods in low-dimensional quantum systems

This is an introductory course to the Lanczos Method and Density Matrix Renormalization Group Algorithms (DMRG), two among the leading numerical techniques applied in studies of low-dimensional quantum models. The idea of studying the models on clusters of a finite size in order to extract their physical properties is briefly discussed. The important role played by the model symmetries is also examined. Special emphasis is given to the DMRG.

Competition between local potentials and attractive particle-particle interactions in superlattices

Naturally occuring or man-made systems displaying periodic spatial modulations of their properties on a nanoscale constitute superlattices. Such modulated structures are important both as prototypes of simple nanotechnological devices and as particular examples of emerging spatial inhomogeneity in interacting many-electron systems. Here we investigate the effect different types of modulation of the system parameters have on the ground-state energy and the charge-density distribution of the system. The superlattices are described by the inhomogeneous attractive Hubbard model, and the calculations are performed by density-functional and density-matrix renormalization group techniques. We find that modulations in local electric potentials are much more effective in shaping the system's properties than modulations in the attractive on-site interaction. This is the same conclusion we previously [M.F. Silva, N.A. Lima, A.L. Malvezzi, K. Capelle, Phys. Rev. B 71 (2005) 125130.] obtained for repulsive interactions, suggesting that it is not an artifact of a specific state, but a general property of modulated structures. (c) 2007 Elsevier Ltd. All rights reserved.

Absence of true long-range orbital order in a two-leg Kondo ladder

We investigate, through the density-matrix renormalization group and the Lanczos technique, the possibility of a two-leg Kondo ladder presenting an incommensurate orbital order. Our results indicate staggered short-range orbital order at half-filling. Away from half-filling our data are consistent with incommensurate quasi-long-range orbital order. We also observed that an interaction between the localized spins enhances the rung-rung current correlations.

Effects of nanoscale spatial inhomogeneity in strongly correlated systems

We calculate ground-state energies and density distributions of Hubbard superlattices characterized by periodic modulations of the on-site interaction and the on-site potential. Both density-matrix renormalization group and density-functional methods are employed and compared. We find that small variations in the on-site potential v(i) can simulate, cancel, or even overcompensate effects due to much larger variations in the on-site interaction U-i. Our findings highlight the importance of nanoscale spatial inhomogeneity in strongly correlated systems, and call for a reexamination of model calculations assuming spatial homogeneity.

Spin, charge, and orbital correlations in the one-dimensional t(2g)-orbital Hubbard model

We present the zero-temperature phase diagram of the one-dimensional t(2g)-orbital Hubbard model, obtained using the density-matrix renormalization group and Lanczos techniques. Emphasis is given to the case of the electron density n=5 corresponding to five electrons per site, while several other cases for electron densities between n=3 and 6 are also studied. At n=5, our results indicate a first-order transition between a paramagnetic (PM) insulator phase, with power-law slowly decaying correlations, and a fully polarized ferromagnetic (FM) state by tuning the Hund's coupling. The results also suggest a transition from the n=5 PM insulator phase to a metallic regime by changing the electron density, either via hole or electron doping. The behavior of the spin, charge, and orbital correlation functions in the FM and PM states are also described in the text and discussed. The robustness of these two states against varying parameters suggests that they may be of relevance in quasi-one-dimensional Co-oxide materials, or even in higher dimensional cobaltite systems as well.

Modulation of charge-density waves by superlattice structures

We discuss the interplay between electronic correlations and an underlying superlattice structure in determining the period of charge density waves (CDW's), by considering a one-dimensional Hubbard model with a repeated (nonrandom) pattern of repulsive (U > 0) and free (U=0) sites. Density matrix renormalization group diagonalization of finite systems (up to 120 sites) is used to calculate the charge-density correlation function and structure factor in the ground state. The modulation period can still be predicted through effective Fermi wave vectors k(F)(*) and densities, and we have found that it is much more sensitive to electron (or hole) doping, both because of the narrow range of densities needed to go from q(*)=0 to pi, but also due to sharp 2k(F)(*)-4k(F)(*) transitions; these features render CDW's more versatile for actual applications in heterostructures than in homogeneous systems.

Multiperiodic magnetic structures in Hubbard superlattices

We consider fermions in one-dimensional superlattices (SL's), modeled by site-dependent Hubbard-U couplings arranged in a repeated pattern of repulsive (i.e., U>0) and free (U=0) sites. Density matrix renormalization group diagonalization of finite systems is used to calculate the local moment and the magnetic structure factor in the ground state. We have found four regimes for magnetic behavior: uniform local moments forming a spin-density wave (SDW), floppy local moments with short-ranged correlations, local moments on repulsive sites forming long-period SDW's superimposed with short-ranged correlations, and local moments on repulsive sites solely with long-period SDW's; the boundaries between these regimes depend on the range of electronic densities ρ and on the SL aspect ratio. Above a critical electronic density, ρ↑↓, the SDW period oscillates both with ρ and with the spacer thickness. The former oscillation allows one to reproduce all SDW wave vectors within a small range of electronic densities, unlike the homogeneous system. The latter oscillation is related to the exchange oscillation observed in magnetic multilayers. A crossover between regimes of thin to thick layers has also been observed.

Scaling limit of virtual states of triatomic systems

Natural scales determine the physics of quantum few-body systems with short-range interactions. Thus, the scaling limit is found when the ratio between the scattering length and the interaction range tends to infinity, while the ratio between the physical scales are kept fixed. From the formal point of view, the relation of the scaling limit and the renormalization aspects of a few-body model with a zero-range interaction, through the derivation of subtracted three-body T-matrix equations that are renormalization-group invariant.

Relating a gluon mass scale to an infrared fixed point in pure gauge QCD

The condition for the global minimum of the vacuum energy for a non-Abelian gauge theory with a dynamically generated gauge boson mass scale which implies the existence of a nontrivial IR fixed point of the theory was shown. Thus, this vacuum energy depends on the dynamical masses through the nonperturbative propagators of the theory. The results show that the freezing of the QCD coupling constant observed in the calculations can be a natural consequence of the onset of a gluon mass scale, giving strong support to their claim.

The temperature dependence of the QCD running coupling

We study the running of the QCD coupling with the momentum squared (Q 2) and the temperature scales in the high temperature limit (T > Tc), using a mass dependent renormalization scheme to build the Renormalization Group Equations. The approach used guaranty gauge invariance, through the use of the Hard Thermal Loop approximation, and independence of the vertex chosen to renormalize the coupling. In general, the dependence of the coupling with the temperature is not logarithmical, although in the region Q2 ∼ T2 the logarithm approximation is reasonable. Finally, as known from Debye screening, color charge is screened in the coupling. The number of flavors, however, is anti-screened.

Gluon mass generation in the PT-BFM scheme

In this article we study the general structure and special properties of the Schwinger-Dyson equation for the gluon propagator constructed with the pinch technique, together with the question of how to obtain infrared finite solutions, associated with the generation of an effective gluon mass. Exploiting the known all-order correspondence between the pinch technique and the background field method, we demonstrate that, contrary to the standard formulation, the non-perturbative gluon self-energy is transverse order-by-order in the dressed loop expansion, and separately for gluonic and ghost contributions. We next present a comprehensive review of several subtle issues relevant to the search of infrared finite solutions, paying particular attention to the role of the seagull graph in enforcing transversality, the necessity of introducing massless poles in the three-gluon vertex, and the incorporation of the correct renormalization group properties. In addition, we present a method for regulating the seagull-type contributions based on dimensional regularization; its applicability depends crucially on the asymptotic behavior of the solutions in the deep ultraviolet, and in particular on the anomalous dimension of the dynamically generated gluon mass. A linearized version of the truncated Schwinger-Dyson equation is derived, using a vertex that satisfies the required Ward identity and contains massless poles belonging to different Lorentz structures. The resulting integral equation is then solved numerically, the infrared and ultraviolet properties of the obtained solutions are examined in detail, and the allowed range for the effective gluon mass is determined. Various open questions and possible connections with different approaches in the literature are discussed. © SISSA 2006.

Analyzing dynamical gluon mass generation

We study the necessary conditions for obtaining infrared finite solutions from the Schwinger-Dyson equation governing the dynamics of the gluon propagator. The equation in question is set up in the Feynman gauge of the background field method, thus capturing a number of desirable features. Most notably, and in contradistinction to the standard formulation, the gluon self-energy is transverse order-by-order in the dressed loop expansion, and separately for gluonic and ghost contributions. Various subtle field-theoretic issues, such as renormalization group invariance and regularization of quadratic divergences, are briefly addressed. The infrared and ultraviolet properties of the obtained solutions are examined in detail, and the allowed range for the effective gluon mass is presented.

Thermoelectric effects in quantum dots

We report a numerical renormalization-group study of the thermoelectric effect in the single-electron transistor (SET) and side-coupled geometries. As expected, the computed thermal conductance and thermopower curves show signatures of the Kondo effect and of Fano interference. The thermopower curves are also affected by particle-hole asymmetry. © 2009 Elsevier B.V. All rights reserved.

Universal conductance for the Anderson model

We discuss the thermal dependence of the zero-bias electrical conductance for a quantum dot embedded in a quantum wire, or side-coupled to it. In the Kondo regime, the temperature-dependent conductances map linearly onto the conductance for the symmetric Anderson Hamiltonian. The mapping fits accurately numerical renormalization-group results for the conductance in each geometry. In the side-coupled geometry, the conductance is markedly affected by a gate potential applied to the wire; in the embedded geometry, it is not. © 2010 IOP Publishing Ltd.

As funções do grupo de renormalização no modelo com auto-interação quártica fermiônica

Modelos com interações quárticas fermiônicas tem sido estudadas para clarificar aspectos conceituais e possíveis aplicações em teoria quântica de campos. Neste trabalho apresentamos a estrutura do grupo de renormalização no modelo de Nambu-Jona-Lasinio até a ordem de 1-loop. O modelo é não renormalizável perturbativamente, no sentido usual de contagem de potência, mas é tratado como uma teoria efetiva, válida numa escala de energia onde p << ^, sendo p o momento externo do loop e ^ um parâmetro de escala de massa que caracteriza o acoplamento do vértice não renormalizável. Esclarecemos a estrutura tensorial dos vértices de interação e calculamos as funções do grupo de renormalização. A análise dos pontos fixos da teoria também é apresentada e discutida usando o formalismo de redução das constantes de acoplamento proposto por Zimmermann. Encontramos a baixas eneergias a origem como ponto fixo infravermelho estável e um ponto fixo não trivial ultravioleta estável, indicando a consistência perturbativa se o momento é pequeno.