Source author record

Enrico Arrigoni

Enrico Arrigoni 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

30works
6topics
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

30 published item(s)

preprint2020arXiv

Pseudogap Anderson impurity model out of equilibrium: A master equation tensor network approach

We study equilibrium and nonequilibrium properties of the single-impurity Anderson model with a power-law pseudogap in the density of states. In equilibrium, the model is known to display a quantum phase transition from a generalized Kondo to a local moment phase. In the present work, we focus on the extension of these phases beyond equilibrium, i.e. under the influence of a bias voltage. Within the auxiliary master equation approach combined with a scheme based on matrix product states (MPS) we are able to directly address the current-carrying steady state. Starting with the equilibrium situation, we first corroborate our results by comparing with a direct numerical evaluation of ground state spectral properties of the system by MPS. Here, a scheme to locate the phase boundary by extrapolating the power-law exponent of the self energy produces a very good agreement with previous results obtained by the numerical renormalization group. Our nonequilibrium study as a function of the applied bias voltage is then carried out for two points on either side of the phase boundary. In the Kondo regime the resonance in the spectral function is splitted as a function of the increasing bias voltage. The local moment regime, instead, displays a dip in the spectrum near the position of the chemical potentials. Similar features are observed in the corresponding self energies. The Kondo split peaks approximately obey a power-law behavior as a function of frequency, whose exponents depend only slightly on voltage. Finally, the differential conductance in the Kondo regime shows a peculiar maximum at finite voltages, whose height, however, is below the accuracy level.

preprint2016arXiv

Lindblad-Driven Discretized Leads for Non-Equilibrium Steady-State Transport in Quantum Impurity Models: Recovering the Continuum Limit

The description of interacting quantum impurity models in steady-state nonequilibrium is an open challenge for computational many-particle methods: the numerical requirement of using a finite number of lead levels and the physical requirement of describing a truly open quantum system are seemingly incompatible. One possibility to bridge this gap is the use of Lindblad-driven discretized leads (LDDL): one couples auxiliary continuous reservoirs to the discretized lead levels and represents these additional reservoirs by Lindblad terms in the Liouville equation. For quadratic models governed by Lindbladian dynamics, we present an elementary approach for obtaining correlation functions analytically. In a second part, we use this approach to explicitly discuss the conditions under which the continuum limit of the LDDL approach recovers the correct representation of thermal reservoirs. As an analytically solvable example, the nonequilibrium resonant level model is studied in greater detail. Lastly, we present ideas towards a numerical evaluation of the suggested Lindblad equation for interacting impurities based on matrix product states. In particular, we present a reformulation of the Lindblad equation, which has the useful property that the leads can be mapped onto a chain where both the Hamiltonian dynamics and the Lindblad driving are local at the same time. Moreover, we discuss the possibility to combine the Lindblad approach with a logarithmic discretization needed for the exploration of exponentially small energy scales.

preprint2016arXiv

Non-equilibrium variational cluster perturbation theory: quench dynamics of the quantum Ising model

We introduce a variational implementation of cluster perturbation theory (CPT) to address the dynamics of spin systems driven out of equilibrium. We benchmark the method with the quantum Ising model subject to a sudden quench of the transverse magnetic field across the transition or within a phase. We treat both the one-dimensional case, for which an exact solution is available, as well the two-dimensional one, for which has to resort to numerical results. Comparison with exact results shows that the approach provides a quite accurate description of the real-time dynamics up to a characteristic time scale $τ$ that increses with the size of the cluster used for CPT. In addition, and not surprisingly $τ$ is small for quenches across the equilibrium phase transition, but can be quite larger for quenches within the ordered or disordered phases.

preprint2016arXiv

Resonance Effects in Correlated Multilayer Heterostructures

We study the occurrence of negative differential conductance induced by resonance effects in a model for a multilayer heterostructure. In particular, we consider a system consisting of several correlated and non-correlated monoatomic layers, sandwiched between two metallic leads. The geometry confines electrons in wells within the heterostructures, which are connected to each other and to the leads by tunneling processes. The non-equilibrium situation is produced by applying a bias-voltage to the leads. Our results show that for specific values of the parameters resonance tunneling takes place. We investigate in detail its influence on the current-voltage characteristics. Our results are obtained via non-equilibrium real-space dynamical mean-field theory. As an impurity solver we use the so-called auxiliary master equation approach, which addresses the impurity problem within an auxiliary system consisting of a correlated impurity, a small number of uncorrelated bath sites, and two Markovian environments described by a generalized master equation.

preprint2015arXiv

Auxiliary master equation approach within matrix product states: Spectral properties of the nonequilibrium Anderson impurity model

Within the recently introduced auxiliary master equation approach it is possible to address steady state properties of strongly correlated impurity models, small molecules or clusters efficiently and with high accuracy. It is particularly suited for dynamical mean field theory in the nonequilibrium as well as in the equilibrium case. The method is based on the solution of an auxiliary open quantum system, which can be made quickly equivalent to the original impurity problem. In its first implementation a Krylov space method was employed. Here, we aim at extending the capabilities of the approach by adopting matrix product states for the solution of the corresponding auxiliary quantum master equation. This allows for a drastic increase in accuracy and permits us to access the Kondo regime for large values of the interaction. In particular, we investigate the nonequilibrium steady state of a single impurity Anderson model and focus on the spectral properties for temperatures T below the Kondo temperature T_K and for small bias voltages phi. For the two cases considered, with T=T_K/4 and T=T_K/10 we find a clear splitting of the Kondo resonance into a two-peak structure for phi close above T_K. In the equilibrium case (phi=0) and for T=T_K/4, the obtained spectral function essentially coincides with the one from numerical renormalization group.

preprint2015arXiv

Current characteristics of a one-dimensional Hubbard chain: The role of correlation and dissipation

We study the electronic transport in an infinite one-dimensional Hubbard chain, driven by a homogeneous electric field. The physical chain is coupled to fermion bath chains, in order to account for dissipation and to prevent the occurrence of Bloch Oscillations. The steady state current is computed in the frame of Keldysh Green's functions in Cluster Perturbation Theory. The current characteristics are dominated by resonant-tunneling-like structures, which can be traced back to Wannier-Stark resonances due to anti-ferromagnetic correlations. The same current characteristic occurs in a non-interacting Wannier-Stark model with alternating on-site energies. Non-local effects of the self energy can be accounted for the observed physical behaviour.

preprint2015arXiv

Master equation based steady-state cluster perturbation theory

A simple and efficient approximation scheme to study electronic transport characteristics of strongly correlated nano devices, molecular junctions or heterostructures out of equilibrium is provided by steady-state cluster perturbation theory. In this work, we improve the starting point of this perturbative, nonequilibrium Green's function based method. Specifically, we employ an improved unperturbed (so-called reference) state $\hatρ^S$, constructed as the steady-state of a quantum master equation within the Born-Markov approximation. This resulting hybrid method inherits beneficial aspects of both, the quantum master equation as well as the nonequilibrium Green's function technique. We benchmark the new scheme on two experimentally relevant systems in the single-electron transistor regime: An electron-electron interaction based quantum diode and a triple quantum dot ring junction, which both feature negative differential conductance. The results of the new method improve significantly with respect to the plain quantum maste equation treatment at modest additional computational cost.

preprint2015arXiv

Nonequilibrium, spatio-temporal formation of the Kondo screening-cloud on a lattice

We study the nonequilibrium formation of a spin screening cloud that accompanies the quenching of a local magnetic moment immersed in a Fermi sea at zero temperature. Based on high precision density matrix renormalization group results for the interacting single impurity Anderson model we discuss the real time evolution after a quantum quench in the impurity-reservoir hybridization using time evolving block decimation. We report emergent length and time scales in the spatiotemporal structure of non-local correlation functions in the spin- and the charge density channel. At equilibrium, our data for the correlation functions and the extracted length-scales show good agreement with existing results, as do local time dependent observables at the impurity. In the time-dependent data, we identify a major signal which defines a "light cone" moving at the Fermi velocity and a ferromagnetic component in its wake. Inside the light cone we find that the structure of the nonequilibrium correlation functions emerges on two time scales. Initially, the qualitative structure of the correlation functions develops rapidly at the lattice Fermi velocity. Subsequently the spin correlations converge to the equilibrium results on a much larger time scale. This process sets a dynamic energy scale, which we identify to be proportional to the Kondo temperature. Outside the light cone we observe two different power law decays of the correlation functions in space, with time and interaction strength independent exponents.

preprint2015arXiv

Quasiparticle excitations in steady state transport across a correlated layer

In this work we investigate the spectral and transport properties of a single correlated layer attached to two metallic leads, with particular focus on the low-energy physics. A steady state current is driven across the layer by applying a bias voltage between the leads. Extending previous work we introduce a nonzero temperature in the leads, which enables us to study the influence of quasiparticle excitations on the transport characteristics in detail. Even though the system is clearly three dimensional we obtain current-voltage curves that closely resemble those of single quantum dots. Furthermore, a splitting of the quasiparticle excitation with bias voltage is observed in the spectral function.

preprint2015arXiv

Transport Through a Correlated Interface: Auxiliary Master Equation Approach

We present improvements of a recently introduced numerical method [Arrigoni etal, Phys. Rev. Lett. 110, 086403 (2013)] to compute steady state properties of strongly correlated electronic systems out of equilibrium. The method can be considered as a non-equilibrium generalization of exact diagonalization based dynamical mean-field theory (DMFT). The key modification for the non-equilibrium situation consists in addressing the DMFT impurity problem within an auxiliary system consisting of the correlated impurity, $N_b$ uncorrelated bath sites and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of $N_b$ than previously in DMFT. This increases the accuracy of the results and is crucial for a correct description of the physical behavior of the system in the relevant parameter range including a semi-quantitative description of the Kondo regime. To illustrate the approach we consider a monoatomic layer of correlated orbitals, described by the single-band Hubbard model, attached to two metallic leads. The non-equilibrium situation is driven by a bias-voltage applied to the leads. For this system, we investigate the spectral function and the steady state current-voltage characteristics in the weakly as well as in the strongly interacting limit. In particular we investigate the non-equilibrium behavior of quasi-particle excitations within the Mott gap of the correlated layer. We find for low bias voltage Kondo like behavior in the vicinity of the insulating phase. In particular we observe a splitting of the Kondo resonance as a function of the bias voltage.

preprint2014arXiv

Auxiliary master equation approach to non-equilibrium correlated impurities

We present a numerical method for the study of correlated quantum impurity problems out of equilibrium, which is particularly suited to address steady state properties within Dynamical Mean Field Theory. The approach, recently introduced in [Arrigoni et al., Phys. Rev. Lett. 110, 086403 (2013)], is based upon a mapping of the original impurity problem onto an auxiliary open quantum system, consisting of the interacting impurity coupled to bath sites as well as to a Markovian environment. The dynamics of the auxiliary system is governed by a Lindblad master equation whose parameters are used to optimize the mapping. The accuracy of the results can be readily estimated and systematically improved by increasing the number of auxiliary bath sites, or by introducing a linear correction. Here, we focus on a detailed discussion of the proposed approach including technical remarks. To solve for the Green's functions of the auxiliary impurity problem, a non-hermitian Lanczos diagonalization is applied. As a benchmark, results for the steady state current-voltage characteristics of the single impurity Anderson model are presented. Furthermore, the bias dependence of the single particle spectral function and the splitting of the Kondo resonance are discussed. In its present form the method is fast, efficient and features a controlled accuracy.

preprint2014arXiv

Effects of electronic correlations and magnetic field on a molecular ring out of equilibrium

We study effects of electron-electron interactions on the steady-state characteristics of a hexagonal molecular ring in a magnetic field, as a model for a benzene molecular junction. The system is driven out of equilibrium by applying a bias voltage across two metallic leads. We employ a model Hamiltonian approach to evaluate the effects of on-site as well as nearest-neighbour density-density type interactions in a physically relevant parameter regime. Results for the steady-state current, charge density and magnetization in three different junction setups (para, meta and ortho) are presented. Our findings indicate that interactions beyond the mean-field level renormalize voltage thresholds as well as current plateaus. Electron-electron interactions lead to substantial charge redistribution as compared to the mean-field results. We identify a strong response of the circular current on the electronic structure of the metallic leads. Our results are obtained by steady-state Cluster Perturbation Theory, a systematically improvable approximation to study interacting molecular junctions out of equilibrium, even in magnetic fields. Within this framework general expressions for the current, charge density and magnetization in the steady-state are derived. The method is flexible and fast and can straight-forwardly be applied to effective models as obtained from ab-initio calculations.

preprint2013arXiv

Nonequilibrium self-energy functional theory

The self-energy functional theory (SFT) is generalized to describe the real-time dynamics of correlated lattice-fermion models far from thermal equilibrium. This is achieved by means of a reformulation of the original equilibrium theory in terms of double-time Green's functions on the Keldysh-Matsubara contour. We construct a functional which is stationary at the physical (nonequilibrium) self-energy and which yields the grand potential of the initial thermal state at the physical point. Non-perturbative approximations can be defined by specifying a reference system that serves to generate trial self-energies. These self-energies are varied by varying the reference system's one-particle parameters on the Keldysh-Matsubara contour. In case of thermal equilibrium, the new approach reduces to the conventional SFT. However, "unphysical" variations, i.e., variations that are different on the upper and the lower branch of the Keldysh contour, must be considered to fix the time-dependence of the optimal physical parameters via the variational principle. Functional derivatives in the nonequilibrium SFT Euler equation are carried out analytically to derive conditional equations for the variational parameters that are accessible to a numerical evaluation via a time-propagation scheme. Approximations constructed by means of the nonequilibrium SFT are shown to be inherently causal, internally consistent and to respect macroscopic conservation laws resulting from gauge symmetries of the Hamiltonian. This comprises the nonequilibrium dynamical mean-field theory but also dynamical-impurity and variational-cluster approximations that are specified by reference systems with a finite number of degrees of freedom. In this way, non-perturbative and consistent approximations can be set up, the numerical evaluation of which is accessible to an exact-diagonalization approach.

preprint2013arXiv

Steady-state and quench dependent relaxation of a quantum dot coupled to one-dimensional leads

We study the time evolution and steady state of the charge current in a single-impurity Anderson model, using matrix product states techniques. A nonequilibrium situation is imposed by applying a bias voltage across one-dimensional tight-binding leads. Focusing on particle-hole symmetry, we extract current-voltage characteristics from universal low-bias up to high-bias regimes, where band effects start to play a dominant role. We discuss three quenches, which after strongly quench-dependent transients yield the same steady-state current. Among these quenches we identify those favorable for extracting steady-state observables. The period of short-time oscillations is shown to compare well to real-time renormalization group results for a simpler model of spinless fermions. We find indications that many-body effects play an important role at high-bias voltage and finite bandwidth of the metallic leads. The growth of entanglement entropy after a certain time scale (proportional to the inverse of Delta) is the major limiting factor for calculating the time evolution. We show that the magnitude of the steady-state current positively correlates with entanglement entropy. The role of high-energy states for the steady-state current is explored by considering a damping term in the time evolution.

preprint2013arXiv

Vibration-mediated correlation effects in the transport properties of a benzene molecule

We theoretically analyze correlation effects on the transport properties of a benzene molecule that are mediated by interactions between the motion of the nuclei and the transmitted charge. We focus on the lowest-lying molecular vibrational mode which allows us to derive an analytic expression for the current. The results provide transparent interpretations of various features of the highly nonlinear current-voltage characteristics, which is experimentally accessible through resonant inelastic electron-tunneling spectroscopy.

preprint2012arXiv

Characterization of Mott-insulating and superfluid phases in the one-dimensional Bose--Hubbard model

We use strong-coupling perturbation theory, the variational cluster approach (VCA), and the dynamical density-matrix renormalization group (DDMRG) method to investigate static and dynamical properties of the one-dimensional Bose--Hubbard model in both the Mott-insulating and superfluid phases. From the von Neumann entanglement entropy we determine the central charge and the transition points for the first two Mott lobes. Our DMRG results for the ground-state energy, momentum distribution function, boson correlation function decay, Mott gap, and single particle-spectral function are reproduced very well by the strong-coupling expansion to fifth order, and by VCA with clusters up to 12 sites as long as the ratio between the hopping amplitude and on-site repulsion, t/U, is smaller than 0.15 and 0.25, respectively. In addition, in the superfluid phase VCA captures well the ground-state energy and the sound velocity of the linear phonon modes. This comparison provides an authoritative estimate for the range of applicability of these methods. In strong-coupling theory for the Mott phase, the dynamical structure factor is obtained from the solution of an effective single-particle problem with an attractive potential. The resulting resonances show up as double-peak structure close to the Brillouin zone boundary. These high-energy features also appear in the superfluid phase which is characterized by a pronounced phonon mode at small momenta and energies, as predicted by Bogoliubov and field theory. In one dimension, there are no traces of an amplitude mode in the dynamical single-particle and two-particle correlation functions.

preprint2012arXiv

Nonequilibrium Dynamical Mean Field Theory: an auxiliary Quantum Master Equation approach

We introduce a versatile method to compute electronic steady state properties of strongly correlated extended quantum systems out of equilibrium. The approach is based on dynamical mean-field theory (DMFT), in which the original system is mapped onto an auxiliary non-equilibrium impurity problem imbedded in a Markovian environment. The steady state Green's function of the auxiliary system is solved by full diagonalization of the corresponding Lindblad equation. The approach can be regarded as the nontrivial extension of the exact-diagonalization based DMFT to the non-equilibrium case. As a first application, we consider an interacting Hubbard layer attached to two metallic leads and present results for the steady-state current and the non-equilibrium density of states.

preprint2012arXiv

Steady-state spectra, current and stability diagram of a quantum dot: a non-equilibrium Variational Cluster Approach

We calculate steady-state properties of a strongly correlated quantum dot under voltage bias by means of non-equilibrium Cluster Perturbation Theory and the non-equilibrium Variational Cluster Approach, respectively. Results for the steady-state current are benchmarked against data from accurate Matrix Product State based time evolution. We show that for low to medium interaction strength, non-equilibrium Cluster Perturbation Theory already yields good results, while for higher interaction strength the self-consistent feedback of the non-equilibrium Variational Cluster Approach significantly enhances the accuracy. We report the current-voltage characteristics for different interaction strengths. Furthermore we investigate the non-equilibrium local density of states of the quantum dot and illustrate that within the variational approach a linear splitting and broadening of the Kondo resonance is predicted which depends on interaction strength. Calculations with applied gate voltage, away from particle hole symmetry, reveal that the maximum current is reached at the crossover from the Kondo regime to the doubly-occupied or empty quantum dot. Obtained stability diagrams compare very well to recent experimental data [Phys. Rev. B, 84, 245316 (2011)].

preprint2012arXiv

Variational cluster approach to the single impurity Anderson model

We study the single impurity Anderson model by means of cluster perturbation theory and the variational cluster approach (VCA). An expression for the VCA grand potential for a system in a non interacting bath is presented. Results for the single-particle dynamics in different parameter regimes are shown to be in good agreement with established renormalization group results. We aim at a broad and comprehensive overview of the capabilities and shortcomings of the methods. We address the question to what extent the elusive low energy properties of the model are reproducible within the framework of VCA. These are furthermore benchmarked against continuous time quantum Monte Carlo calculations. We also discuss results obtained by an alternative, i.e. self consistent formulation of VCA, which was introduced recently in the context of non equilibrium systems.

preprint2011arXiv

Emission characteristics of laser-driven dissipative coupled-cavity systems

We consider a laser-driven and dissipative system of two coupled cavities with Jaynes-Cummings nonlinearity. In particular, we investigate both incoherent and coherent laser driving, corresponding to different experimental situations. We employ Arnoldi time evolution as a numerical tool to solve exactly the many-body master equation describing the non-equilibrium quantum system. We evaluate the fluorescence spectrum and the spectrum of the second-order correlation function of the emitted light field. Finally, we relate the measured spectra of the dissipative quantum system to excitations of the corresponding non-dissipative quantum system. Our results demonstrate how to interpret spectra obtained from dissipative quantum systems and specify what information is contained therein.

preprint2011arXiv

Extended self-energy functional approach for strongly-correlated lattice bosons in the superfluid phase

Among the various numerical techniques to study the physics of strongly correlated quantum many-body systems, the self-energy functional approach (SFA) has become increasingly important. In its previous form, however, SFA is not applicable to Bose-Einstein condensation or superfluidity. In this paper we show how to overcome this shortcoming. To this end we identify an appropriate quantity, which we term $D$, that represents the correlation correction of the condensate order parameter, as it does the self-energy for the Green's function. An appropriate functional is derived, which is stationary at the exact physical realizations of $D$ and of the self-energy. Its derivation is based on a functional-integral representation of the grand potential followed by an appropriate sequence of Legendre transformations. The approach is not perturbative and therefore applicable to a wide range of models with local interactions. We show that the variational cluster approach based on the extended self-energy functional is equivalent to the "pseudoparticle" approach introduced in Phys. Rev. B, 83, 134507 (2011). We present results for the superfluid density in the two-dimensional Bose-Hubbard model, which show a remarkable agreement with those of Quantum-Monte-Carlo calculations.

preprint2011arXiv

Nonequilibrium steady state for strongly-correlated many-body systems: variational cluster approach

A numerical approach is presented that allows to compute nonequilibrium steady state properties of strongly correlated quantum many-body systems. The method is imbedded in the Keldysh Green's function formalism and is based upon the idea of the variational cluster approach as far as the treatment of strong correlations is concerned. It appears that the variational aspect is crucial as it allows for a suitable optimization of a "reference" system to the nonequilibrium target state. The approach is neither perturbative in the many-body interaction nor in the field, that drives the system out of equilibrium, and it allows to study strong perturbations and nonlinear responses of systems in which also the correlated region is spatially extended. We apply the presented approach to non-linear transport across a strongly correlated quantum wire described by the fermionic Hubbard model. We illustrate how the method bridges to cluster dynamical mean-field theory upon coupling two baths containing and increasing number of uncorrelated sites.

preprint2011arXiv

Polaritonic properties of the Jaynes-Cummings lattice model in two dimensions

Light-matter systems allow to realize a strongly correlated phase where photons are present. In these systems strong correlations are achieved by optical nonlinearities, which appear due to the coupling of photons to atomic-like structures. This leads to intriguing effects, such as the quantum phase transition from the Mott to the superfluid phase. Here, we address the two-dimensional Jaynes-Cummings lattice model. We evaluate the boundary of the quantum phase transition and study polaritonic properties. In order to be able to characterize polaritons, we investigate the spectral properties of both photons as well as two-level excitations. Based on this information we introduce polariton quasiparticles as appropriate wavevector, band index, and filling dependent superpositions of photons and two-level excitations. Finally, we analyze the contributions of the individual constituents to the polariton quasiparticles.

preprint2011arXiv

Time-reversal symmetry breaking phase in the Hubbard model: a VCA study

We study the stability of the time-reversal symmetry breaking staggered flux phase of a single band Hubbard model, within the Variational Cluster Approach (VCA). For intermediate and small values of the interaction $U$, we find metastable solutions for the staggered flux phase, with a maximum current per bond at $U\approx 3.2$. However, allowing for antiferromagnetic and superconducting long-range order it turns out that in the region at and close to half filling the antiferromagnetic phase is the most favorable energetically. The effect of nearest-neighbour interaction is also considered. Our results show that a negative nearest-neighbour interaction and finite doping favors the stability of the staggered-flux phase. We also present preliminary results for the three-band Hubbard model obtained with a restricted set of variational parameters. For this case, no spontaneous time-reversal symmetry breaking phase is found in our calculations.

preprint2011arXiv

Variational cluster approach for strongly correlated lattice bosons in the superfluid phase

We extend the variational cluster approach to deal with strongly correlated lattice bosons in the superfluid phase. To this end, we reformulate the approach within a pseudoparticle formalism, whereby cluster excitations are described by particlelike excitations. The approximation amounts to solving a multicomponent noninteracting bosonic system by means of a multimode Bogoliubov approximation. A source-and-drain term is introduced in order to break U(1) symmetry at the cluster level. We provide an expression for the grand potential, the single-particle normal and anomalous Green's functions, the condensate density, and other static quantities. As a first nontrivial application of the method we choose the two-dimensional Bose-Hubbard model and evaluate results in both the Mott and the superfluid phases. Our results show an excellent agreement with quantum Monte Carlo calculations.

preprint2010arXiv

Benchmarking the variational cluster approach by means of the one-dimensional Bose-Hubbard model

Convergence properties of the variational cluster approach with respect to the variational parameter space, cluster size, and boundary conditions of the reference system are investigated and discussed for bosonic many-body systems. Specifically, the variational cluster approach is applied to the one-dimensional Bose-Hubbard model, which exhibits a quantum phase transition from Mott to superfluid phase. In order to benchmark the variational cluster approach, results for the phase boundary delimiting the first Mott lobe are compared with essentially exact density matrix renormalization group data. Furthermore, static quantities, such as the ground state energy and the one-particle density matrix are compared with high-order strong coupling perturbation theory results. For reference systems with open boundary conditions the variational parameter space is extended by an additional variational parameter which allows for a more uniform particle density on the reference system and thus drastically improves the results. It turns out that the variational cluster approach yields accurate results with relatively low computational effort for both spectral as well as static properties of the one-dimensional Bose-Hubbard model, even at the tip of the first Mott lobe where correlation effects are most pronounced.

preprint2010arXiv

Excitations in disordered bosonic optical lattices

Spectral excitations of ultracold gases of bosonic atoms trapped in one dimensional optical lattices with disorder are investigated by means of the variational cluster approach applied to the Bose-Hubbard model. In particular, qualitatively different disorder distributions typically employed in experiments are considered. The computed spectra exhibit a strong dependence on both the shape of the disorder distribution and the disorder strength. We compare alternative results for the Mott gap obtained from its formal definition and from the minimum peak distance, which is the quantity available from experiments.

preprint2010arXiv

Quantum phase transition and excitations of the Tavis-Cummings lattice model

The enormous progress in controlling quantum optical and atomic systems has prompted ideas for new experimental realizations of strongly correlated many-body systems operating with light. These systems consist of photons confined in optical cavities, which interact strongly with atoms or atomiclike structures. Due to the interaction between the two particle species optical nonlinearities appear, leading to a quantum phase transition from Mott to superfluid phase. Here, we address the Tavis-Cummings lattice model, which describes light-matter systems containing multiple atomiclike structures in each cavity. In particular, we investigate the phase boundary delimiting Mott from superfluid phase and the elementary excitations of the two-dimensional Tavis-Cummings lattice model in dependence of the number of atomiclike structures per cavity. In order to obtain the results we employ the variational cluster approach. We evaluate spectral functions and densities of states of both particle species, which allows us to characterize the fundamental excitations of light-matter systems. These excitations are termed polaritons and are superpositions of photons and atomic excitations. We introduce polariton quasiparticles as appropriate linear combinations of both particle species and analyze the weights of their constituents. Our results demonstrate the dependence of the quantum phase transition and the elementary excitations on the number of atomiclike structures per cavity and provide thus valuable insight into the physics of light-matter systems.

preprint2010arXiv

Spectral properties of coupled cavity arrays in one dimension

Spectral properties of coupled cavity arrays in one dimension are investigated by means of the variational cluster approach. Coupled cavity arrays consist of two distinct "particles," namely, photons and atomiclike excitations. Spectral functions are evaluated and discussed for both particle types. In addition, densities of states, momentum distributions and spatial correlation functions are presented. Based on this information, polariton "quasiparticles" are introduced as appropriate, wave vector and filling dependent linear combinations of photon and atomiclike particles. Spectral functions and densities of states are evaluated for the polariton quasiparticles, and the weights of their components are analyzed.

preprint2010arXiv

Spectral properties of strongly correlated bosons in two-dimensional optical lattices

Spectral properties of the two-dimensional Bose-Hubbard model, which emulates ultracold gases of atoms confined in optical lattices, are investigated by means of the variational cluster approach. The phase boundary of the quantum phase transition from Mott phase to superfluid phase is calculated and compared to recent work. Moreover the single-particle spectral functions in both the first and the second Mott lobe are presented and the corresponding densities of states and momentum distributions are evaluated. A qualitatively similar intensity distribution of the spectral weight can be observed for spectral functions in the first and the second Mott lobe.