-
Dynamical mean-field theory of the Anderson-Hubbard model with local and non-local disorder in tensor formulation
Authors:
A. Weh,
Y. Zhang,
A. Östlin,
H. Terletska,
D. Bauernfeind,
K. -M. Tam,
H. G. Evertz,
K. Byczuk,
D. Vollhardt,
L. Chioncel
Abstract:
To explore correlated electrons in the presence of local and non-local disorder, the Blackman-Esterling-Berk method for averaging over off-diagonal disorder is implemented into dynamical mean-field theory using tensor notation. The impurity model combining disorder and correlations is solved using the recently developed fork tensor-product state solver, which allows one to calculate the single par…
▽ More
To explore correlated electrons in the presence of local and non-local disorder, the Blackman-Esterling-Berk method for averaging over off-diagonal disorder is implemented into dynamical mean-field theory using tensor notation. The impurity model combining disorder and correlations is solved using the recently developed fork tensor-product state solver, which allows one to calculate the single particle spectral functions on the real-frequency axis. In the absence of off-diagonal hopping, we establish exact bounds of the spectral function of the non-interacting Bethe lattice with coordination number $Z$. In the presence of interaction, the Mott insulating paramagnetic phase of the one-band Hubbard model is computed at zero temperature in alloys with site- and off-diagonal disorder. When the Hubbard $U$ parameter is increased, transitions from an alloy band-insulator through a correlated metal into a Mott insulating phase are found to take place.
△ Less
Submitted 14 May, 2021;
originally announced May 2021.
-
Impact ionization and multiple photon absorptions in the two-dimensional photoexcited Hubbard model
Authors:
F. Maislinger,
H. G. Evertz
Abstract:
We study the non-equilibrium response of a 4x3 Hubbard model at U=8 under the influence of a short electric field pulse, with the main focus on multiple photon excitations and on the change of double occupancy after the pulse. The behavior mainly depends on the driving frequency of the electric field. The largest change of double occupancy occurs during the pulse. For frequencies below the Mott ga…
▽ More
We study the non-equilibrium response of a 4x3 Hubbard model at U=8 under the influence of a short electric field pulse, with the main focus on multiple photon excitations and on the change of double occupancy after the pulse. The behavior mainly depends on the driving frequency of the electric field. The largest change of double occupancy occurs during the pulse. For frequencies below the Mott gap, we observe multiphoton excitations at large field intensities. For frequencies beyond the gap energy, there is a region where Auger recombination reduces the double occupancy after the pulse. Impact ionization (Multi Exciton Generation), namely a growing double occupancy after the pulse, occurs for frequencies larger than twice the Mott gap. From the Loschmidt amplitude we compute the eigenstate spectrum of the quantum state after the pulse, observing multiple distinct photon excitation peaks, in line with expectations from a quasiparticle picture. We introduce a technique with which we analyze the time evolution of double occupancy in each peak individually. The long-term behavior of the double occupancy almost only depends on the absorbed energy, and we explore the connection of this property to the Eigenstate Thermalization Hypothesis.
△ Less
Submitted 3 August, 2020; v1 submitted 31 July, 2020;
originally announced July 2020.
-
Spectral properties of heterostructures containing half-metallic ferromagnets in the presence of local many-body correlations
Authors:
A. Weh,
J. Otsuki,
H. Schnait,
H. G. Evertz,
U. Eckern,
A. I. Lichtenstein,
L. Chioncel
Abstract:
In this work, we investigate models for bulk, bi- and multilayers containing half-metallic ferromagnets (HMFs), at zero and at finite temperature, in order to elucidate the effects of strong electronic correlations on the spectral properties (density of states). Our focus is on the evolution of the finite-temperature many-body induced tails in the half-metallic gap. To this end, the dynamical mean…
▽ More
In this work, we investigate models for bulk, bi- and multilayers containing half-metallic ferromagnets (HMFs), at zero and at finite temperature, in order to elucidate the effects of strong electronic correlations on the spectral properties (density of states). Our focus is on the evolution of the finite-temperature many-body induced tails in the half-metallic gap. To this end, the dynamical mean-field theory (DMFT) is employed. For the bulk, a Bethe lattice model is solved using a matrix product states based impurity solver at zero temperature and a continuous-time quantum Monte Carlo (CT-QMC) solver at finite temperature. We demonstrate numerically, in agreement with the analytical result, that the tails vanish at the Fermi level at zero temperature. In order to study multilayers, taken to be square lattices within the layers, we use the real-space DMFT extension with the CT-QMC impurity solver. For bilayers formed by the HMF with a band or correlated insulator, we find that charge fluctuations between the layers enhance the finite temperature tails. In addition, in the presence of inter-layer hopping, a coherent quasiparticle peak forms in the otherwise correlated insulator. In the multilayer heterostructure setup, we find that by suitably choosing the model parameters, the tails at the HMF/Mott insulator interface can be reduced significantly, and that a high spin polarization is conceivable, even in the presence of long-ranged electrostatic interactions.
△ Less
Submitted 30 October, 2020; v1 submitted 6 July, 2020;
originally announced July 2020.
-
Comparison of MPS based real time evolution algorithms for Anderson Impurity Models
Authors:
Daniel Bauernfeind,
Markus Aichhorn,
Hans Gerd Evertz
Abstract:
We perform a detailed comparison of two Matrix Product States (MPS) based time evolution algorithms for Anderson Impurity Models. To describe the bath, we use both the star-geometry as well as the commonly employed Wilson chain geometry. For each bath geometry, we use either the Time Dependent Variational Principle (TDVP) or the Time Evolving Block Decimation (TEBD) to perform the time evolution.…
▽ More
We perform a detailed comparison of two Matrix Product States (MPS) based time evolution algorithms for Anderson Impurity Models. To describe the bath, we use both the star-geometry as well as the commonly employed Wilson chain geometry. For each bath geometry, we use either the Time Dependent Variational Principle (TDVP) or the Time Evolving Block Decimation (TEBD) to perform the time evolution. To apply TEBD for the star-geometry, we use a specially adapted algorithm that can deal with the long-range coupling terms. Analyzing the major sources of errors, one expects them to be proportional to the system size for all algorithms. Surprisingly, we find errors independent of system size except for TEBD in chain geometry. Additionally, we show that the right combination of bath representation and time evolution algorithm is important. While TDVP in chain geometry is a very precise approach, TEBD in star geometry is much faster, such that for a given accuracy it is superior to TDVP in chain geometry. This makes the adapted version of TEBD in star geometry the most efficient method to solve impurity problems.
△ Less
Submitted 21 June, 2019;
originally announced June 2019.
-
Dynamical Mean-Field Theory on the Real-Frequency Axis: p-d Hybridizations and Atomic Physics in SrMnO$_3$
Authors:
Daniel Bauernfeind,
Robert Triebl,
Manuel Zingl,
Markus Aichhorn,
Hans Gerd Evertz
Abstract:
We investigate the electronic structure of SrMnO$_3$ with Density Functional Theory (DFT) plus Dynamical Mean-Field Theory (DMFT). Within this scheme the selection of the correlated subspace and the construction of the corresponding Wannier functions is a crucial step. Due to the crystal field splitting of the Mn-$3d$ orbitals and their separation from the O-2$p$ bands, SrMnO$_3$ is a material whe…
▽ More
We investigate the electronic structure of SrMnO$_3$ with Density Functional Theory (DFT) plus Dynamical Mean-Field Theory (DMFT). Within this scheme the selection of the correlated subspace and the construction of the corresponding Wannier functions is a crucial step. Due to the crystal field splitting of the Mn-$3d$ orbitals and their separation from the O-2$p$ bands, SrMnO$_3$ is a material where on first sight a 3-band $d$-only model should be sufficient. However, in the present work we demonstrate that the resulting spectrum is considerably influenced by the number of correlated orbitals and the number of bands included in the Wannier function construction. For example, in a $d$-$dp$ model we observe a splitting of the \tg lower Hubbard band into a more complex spectral structure, not observable in $d$-only models. To illustrate these high-frequency differences we employ the recently developed Fork Tensor Product State (FTPS) impurity solver, as it provides the necessary spectral resolution on the real-frequency axis. We find that the spectral structure of a 5-band $d$-$dp$ model is in good agreement with PES and XAS experiments. Our results demonstrate that the FTPS solver is capable of performing full 5-band DMFT calculations directly on the real-frequency axis.
△ Less
Submitted 19 June, 2018; v1 submitted 21 December, 2017;
originally announced December 2017.
-
Interaction effects in a microscopic quantum wire model with strong spin-orbit interaction
Authors:
Georg W. Winkler,
Martin Ganahl,
Dirk Schuricht,
Hans Gerd Evertz,
Sabine Andergassen
Abstract:
We investigate the effect of strong interactions on the spectral properties of quantum wires with strong Rashba spin-orbit interaction in a magnetic field, using a combination of Matrix Product State and bosonization techniques. Quantum wires with strong Rashba spin-orbit interaction and magnetic field exhibit a partial gap in one-half of the conducting modes. Such systems have attracted wide-spre…
▽ More
We investigate the effect of strong interactions on the spectral properties of quantum wires with strong Rashba spin-orbit interaction in a magnetic field, using a combination of Matrix Product State and bosonization techniques. Quantum wires with strong Rashba spin-orbit interaction and magnetic field exhibit a partial gap in one-half of the conducting modes. Such systems have attracted wide-spread experimental and theoretical attention due to their unusual physical properties, among which are spin-dependent transport, or a topological superconducting phase when under the proximity effect of an s-wave superconductor. As a microscopic model for the quantum wire we study an extended Hubbard model with spin-orbit interaction and Zeeman field. We obtain spin resolved spectral densities from the real-time evolution of excitations, and calculate the phase diagram. We find that interactions increase the pseudo gap at $k = 0$ and thus also enhance the Majorana-supporting phase and stabilize the helical spin order. Furthermore, we calculate the optical conductivity and compare it with the low energy spiral Luttinger Liquid result, obtained from field theoretical calculations. With interactions, the optical conductivity is dominated by an excotic excitation of a bound soliton-antisoliton pair known as a breather state. We visualize the oscillating motion of the breather state, which could provide the route to their experimental detection in e.g. cold atom experiments.
△ Less
Submitted 9 April, 2017; v1 submitted 13 January, 2017;
originally announced January 2017.
-
Fork Tensor Product States - Efficient Three Orbital Real Time DMFT Solver
Authors:
Daniel Bauernfeind,
Manuel Zingl,
Robert Triebl,
Markus Aichhorn,
Hans Gerd Evertz
Abstract:
We present a tensor network especially suited for multi-orbital Anderson impurity models and as an impurity solver for multi-orbital dynamical mean-field theory (DMFT). The solver works directly on the real-frequency axis and yields very high spectral resolution at all frequencies. We use a large number $\left(\mathcal{O}(100)\right)$ of bath sites, and therefore achieve an accurate representation…
▽ More
We present a tensor network especially suited for multi-orbital Anderson impurity models and as an impurity solver for multi-orbital dynamical mean-field theory (DMFT). The solver works directly on the real-frequency axis and yields very high spectral resolution at all frequencies. We use a large number $\left(\mathcal{O}(100)\right)$ of bath sites, and therefore achieve an accurate representation of the bath. The solver can treat full rotationally invariant interactions with reasonable numerical effort. We show the efficiency and accuracy of the method by a benchmark for the testbed material SrVO$_3$. There we observe multiplet structures in the high-energy spectrum which are almost impossible to resolve by other multi-orbital methods. The resulting structure of the Hubbard bands can be described as a broadened atomic spectrum with rescaled interaction parameters. Additional features emerge when $U$ is increased. The impurity solver offers a new route to the calculation of precise real-frequency spectral functions of correlated materials.
△ Less
Submitted 14 June, 2017; v1 submitted 16 December, 2016;
originally announced December 2016.
-
Universal front propagation in the quantum Ising chain with domain-wall initial states
Authors:
Viktor Eisler,
Florian Maislinger,
Hans Gerd Evertz
Abstract:
We study the melting of domain walls in the ferromagnetic phase of the transverse Ising chain, created by flipping the order-parameter spins along one-half of the chain. If the initial state is excited by a local operator in terms of Jordan-Wigner fermions, the resulting longitudinal magnetization profiles have a universal character. Namely, after proper rescalings, the profiles in the noncritical…
▽ More
We study the melting of domain walls in the ferromagnetic phase of the transverse Ising chain, created by flipping the order-parameter spins along one-half of the chain. If the initial state is excited by a local operator in terms of Jordan-Wigner fermions, the resulting longitudinal magnetization profiles have a universal character. Namely, after proper rescalings, the profiles in the noncritical Ising chain become identical to those obtained for a critical free-fermion chain starting from a step-like initial state. The relation holds exactly in the entire ferromagnetic phase of the Ising chain and can even be extended to the zero-field XY model by a duality argument. In contrast, for domain-wall excitations that are highly non-local in the fermionic variables, the universality of the magnetization profiles is lost. Nevertheless, for both cases we observe that the entanglement entropy asymptotically saturates at the ground-state value, suggesting a simple form of the steady state.
△ Less
Submitted 8 December, 2016; v1 submitted 5 October, 2016;
originally announced October 2016.
-
Quasi-soliton scattering in quantum spin chains
Authors:
Rogier Vlijm,
Martin Ganahl,
Davide Fioretto,
Michael Brockmann,
Masudul Haque,
Hans Gerd Evertz,
Jean-Sébastien Caux
Abstract:
The quantum scattering of magnon bound states in the anisotropic Heisenberg spin chain is shown to display features similar to the scattering of solitons in classical exactly solvable models. Localized colliding Gaussian wave packets of bound magnons are constructed from string solutions of the Bethe equations and subsequently evolved in time, relying on an algebraic Bethe ansatz based framework f…
▽ More
The quantum scattering of magnon bound states in the anisotropic Heisenberg spin chain is shown to display features similar to the scattering of solitons in classical exactly solvable models. Localized colliding Gaussian wave packets of bound magnons are constructed from string solutions of the Bethe equations and subsequently evolved in time, relying on an algebraic Bethe ansatz based framework for the computation of local expectation values in real space-time. The local magnetization profile shows the trajectories of colliding wave packets of bound magnons, which obtain a spatial displacement upon scattering. Analytic predictions on the displacements for various values of anisotropy and string lengths are derived from scattering theory and Bethe ansatz phase shifts, matching time evolution fits on the displacements. The time evolved block decimation (TEBD) algorithm allows for the study of scattering displacements from spin-block states, showing similar scattering displacement features.
△ Less
Submitted 17 August, 2015; v1 submitted 30 July, 2015;
originally announced July 2015.
-
Auxiliary master equation approach within matrix product states: Spectral properties of the nonequilibrium Anderson impurity model
Authors:
Antonius Dorda,
Martin Ganahl,
Hans Gerd Evertz,
Wolfgang von der Linden,
Enrico Arrigoni
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 December, 2015; v1 submitted 10 July, 2015;
originally announced July 2015.
-
Nonequilibrium, spatio-temporal formation of the Kondo screening-cloud on a lattice
Authors:
Martin Nuss,
Martin Ganahl,
Enrico Arrigoni,
Wolfgang von der Linden,
Hans Gerd Evertz
Abstract:
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 evol…
▽ More
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.
△ Less
Submitted 3 March, 2015; v1 submitted 2 September, 2014;
originally announced September 2014.
-
Efficient DMFT impurity solver using real-time dynamics with Matrix Product States
Authors:
Martin Ganahl,
Markus Aichhorn,
Patrik Thunström,
Karsten Held,
Hans Gerd Evertz,
Frank Verstraete
Abstract:
We propose to calculate spectral functions of quantum impurity models using the Time Evolving Block Decimation (TEBD) for Matrix Product States. The resolution of the spectral function is improved by a so-called linear prediction approach. We apply the method as an impurity solver within the Dynamical Mean Field Theory (DMFT) for the single- and two-band Hubbard model on the Bethe lattice. For the…
▽ More
We propose to calculate spectral functions of quantum impurity models using the Time Evolving Block Decimation (TEBD) for Matrix Product States. The resolution of the spectral function is improved by a so-called linear prediction approach. We apply the method as an impurity solver within the Dynamical Mean Field Theory (DMFT) for the single- and two-band Hubbard model on the Bethe lattice. For the single-band model we observe sharp features at the inner edges of the Hubbard bands. A finite size scaling shows that they remain present in the thermodynamic limit. We analyze the real time-dependence of the double occupation after adding a single electron and observe oscillations at the same energy as the sharp feature in the Hubbard band, indicating a long-lived coherent superposition of states that correspond to the Kondo peak and the side peaks. For a two-band Hubbard model we observe an even richer structure in the Hubbard bands, which cannot be related to a multiplet structure of the impurity, in addition to sharp excitations at the band edges of a type similar to the single-band case.
△ Less
Submitted 6 September, 2015; v1 submitted 26 May, 2014;
originally announced May 2014.
-
Chebyshev expansion for Impurity Models using Matrix Product States
Authors:
Martin Ganahl,
Patrik Thunström,
Frank Verstraete,
Karsten Held,
Hans Gerd Evertz
Abstract:
We improve a recently developed expansion technique for calculating real frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost and, in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than…
▽ More
We improve a recently developed expansion technique for calculating real frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost and, in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory (DMFT). Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group (DDMRG), at lower cost. We also propose an alternative expansion, of 1-exp(-tau H) instead of the usual H, which opens the possibility of using established methods for the time evolution of matrix product states to calculate spectral functions directly.
△ Less
Submitted 5 March, 2014;
originally announced March 2014.
-
Quantum Monte Carlo Simulations for Stacked Spin Ladder Systems Containing Low Concentrations of Non-Magnetic Impurities; Application to the Low Temperature Broadening of NMR-Spectra in SrCu2O3
Authors:
Robert Achleitner,
Hans Gerd Evertz,
Masatoshi Imada,
Ralf Gamillscheg,
Peter Mohn
Abstract:
We present a Quantum Monte Carlo study for Heisenberg spin-1/2 two leg ladder systems doped with non magnetic impurities. The simulations are applied to the doped spin ladder compound SrCu2O3 doped with Zn, where a large broadening of the 65Cu NMR lines has been observed in experiment at low temperatures but far above the Neel temperature. We find that interladder-couplings with a sizeable couplin…
▽ More
We present a Quantum Monte Carlo study for Heisenberg spin-1/2 two leg ladder systems doped with non magnetic impurities. The simulations are applied to the doped spin ladder compound SrCu2O3 doped with Zn, where a large broadening of the 65Cu NMR lines has been observed in experiment at low temperatures but far above the Neel temperature. We find that interladder-couplings with a sizeable coupling in stacking direction are required to describe the line broadening, which can not be explained by considering a single ladder only. Around a single impurity, spin correlations cause an exponentially decaying antiferromagnetic local magnetization in a magnetic field. We develop an effective model for the local magnetization of systems with many randomly distributed impurities, with few parameters which can be extracted out of QMC calculations with a single impurity. The broadening arises from a drag effect, where the magnetization around an impurity works as an effective field for spins on the neighboring ladders, causing a non-exponentially decaying magnetization cloud around the impurity. Our results show that even for impurity concentrations as small as x = 0.001 and x = 0.0025, the broadening effect is large, in good quantitative agreement with experiment. We also develop a simple model for the effective interaction of two impurity spins.
△ Less
Submitted 18 December, 2013; v1 submitted 7 October, 2013;
originally announced October 2013.
-
Quantum Bowling: Particle-hole transmutation in one-dimensional strongly interacting lattice models
Authors:
M. Ganahl,
M. Haque,
H. G. Evertz
Abstract:
We study the scattering of a soliton-like propagating particle with a wall of bound particles, in several strongly interacting one-dimensional lattice models with discrete degrees of freedom. We consider spin-polarized fermions (anisotropic Heisenberg spin chain), the fermionic Hubbard model, and the Bose Hubbard model, using precise numerical time dependent Density Matrix Renormalization Group te…
▽ More
We study the scattering of a soliton-like propagating particle with a wall of bound particles, in several strongly interacting one-dimensional lattice models with discrete degrees of freedom. We consider spin-polarized fermions (anisotropic Heisenberg spin chain), the fermionic Hubbard model, and the Bose Hubbard model, using precise numerical time dependent Density Matrix Renormalization Group techniques. We show that in all integrable models studied, there is no reflection. Instead, an incoming particle experiences particle-hole transmutation upon entry and exit of the wall, and travels inside the wall as a hole, analoguous to Klein tunneling, even though the dispersion is highly nonlinear and there is no external potential. {\em Two} particles are added to the wall on the incoming side and removed on the opposite side. For spin-polarized fermions a single transmitted particle thus shifts the wall by two lattice sites, in complete contrast to classical physics. For both Hubbard models, the wall shifts by one doubly occupied single site. In the nonintegrable models studied, the same process occurs in linear superposition with backscattering events. We demonstrate a corresponding fermionic quantum Newton's cradle and a metamaterial with "tachyonic" modes travelling faster than in an empty system. We present a possible atomic scale signal counter for spintronics. Our scenario should be realizable in future cold atom experiments.
△ Less
Submitted 11 February, 2013;
originally announced February 2013.
-
Steady-state and quench dependent relaxation of a quantum dot coupled to one-dimensional leads
Authors:
Martin Nuss,
Martin Ganahl,
Hans Gerd Evertz,
Enrico Arrigoni,
Wolfgang von der Linden
Abstract:
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…
▽ More
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.
△ Less
Submitted 30 July, 2013; v1 submitted 14 January, 2013;
originally announced January 2013.
-
Steady-state spectra, current and stability diagram of a quantum dot: a non-equilibrium Variational Cluster Approach
Authors:
Martin Nuss,
Christoph Heil,
Martin Ganahl,
Michael Knap,
Hans Gerd Evertz,
Enrico Arrigoni,
Wolfgang von der Linden
Abstract:
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-equ…
▽ More
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)].
△ Less
Submitted 16 October, 2012; v1 submitted 24 July, 2012;
originally announced July 2012.
-
Time Evolution within a Comoving Window: Scaling of signal fronts and magnetization plateaus after a local quench in quantum spin chains
Authors:
V. Zauner,
M. Ganahl,
H. G. Evertz,
T. Nishino
Abstract:
We present a modification of Matrix Product State time evolution to simulate the propagation of signal fronts on infinite one-dimensional systems. We restrict the calculation to a window moving along with a signal, which by the Lieb-Robinson bound is contained within a light cone. Signal fronts can be studied unperturbed and with high precision for much longer times than on finite systems. Entangl…
▽ More
We present a modification of Matrix Product State time evolution to simulate the propagation of signal fronts on infinite one-dimensional systems. We restrict the calculation to a window moving along with a signal, which by the Lieb-Robinson bound is contained within a light cone. Signal fronts can be studied unperturbed and with high precision for much longer times than on finite systems. Entanglement inside the window is naturally small, greatly lowering computational effort. We investigate the time evolution of the transverse field Ising (TFI) model and of the S=1/2 XXZ antiferromagnet in their symmetry broken phases after several different local quantum quenches.
In both models, we observe distinct magnetization plateaus at the signal front for very large times, resembling those previously observed for the particle density of tight binding (TB) fermions. We show that the normalized difference to the magnetization of the ground state exhibits similar scaling behaviour as the density of TB fermions. In the XXZ model there is an additional internal structure of the signal front due to pairing, and wider plateaus with tight binding scaling exponents for the normalized excess magnetization. We also observe parameter dependent interaction effects between individual plateaus, resulting in a slight spatial compression of the plateau widths.
In the TFI model, we additionally find that for an initial Jordan-Wigner domain wall state, the complete time evolution of the normalized excess longitudinal magnetization agrees exactly with the particle density of TB fermions.
△ Less
Submitted 3 July, 2015; v1 submitted 3 July, 2012;
originally announced July 2012.
-
Worm algorithms for the 3-state Potts model with magnetic field and chemical potential
Authors:
Ydalia Delgado,
Hans Gerd Evertz,
Christof Gattringer
Abstract:
We discuss worm algorithms for the 3-state Potts model with external field and chemical potential. The complex phase problem of this system can be overcome by using a flux representation where the new degrees of freedom are dimer and monomer variables. Working with this representation we discuss two different generalizations of the conventional Prokof'ev-Svistunov algorithm suitable for Monte Carl…
▽ More
We discuss worm algorithms for the 3-state Potts model with external field and chemical potential. The complex phase problem of this system can be overcome by using a flux representation where the new degrees of freedom are dimer and monomer variables. Working with this representation we discuss two different generalizations of the conventional Prokof'ev-Svistunov algorithm suitable for Monte Carlo simulations of the model at arbitrary chemical potential and evaluate their performance.
△ Less
Submitted 20 February, 2012;
originally announced February 2012.
-
Observing complex bound states in the spin-1/2 Heisenberg XXZ chain using local quantum quenches
Authors:
Martin Ganahl,
Elias Rabel,
Fabian H. L. Essler,
Hans Gerd Evertz
Abstract:
We consider the non-equilibrium evolution in the spin-1/2 XXZ Heisenberg chain for fixed magnetization after a local quantum quench. This model is equivalent to interacting spinless fermions. Initially an infinite magnetic field is applied to n consecutive sites and the ground state is calculated. At time t=0 the field is switched off and the time evolution of observables such as the z-component o…
▽ More
We consider the non-equilibrium evolution in the spin-1/2 XXZ Heisenberg chain for fixed magnetization after a local quantum quench. This model is equivalent to interacting spinless fermions. Initially an infinite magnetic field is applied to n consecutive sites and the ground state is calculated. At time t=0 the field is switched off and the time evolution of observables such as the z-component of spin is computed using the Time Evolving Block Decimation (TEBD) algorithm. We find that the observables exhibit strong signatures of linearly propagating spinon and bound state excitations. These persist even when integrability-breaking perturbations are included. Since bound states ("strings") are notoriously difficult to observe using conventional probes such as inelastic neutron scattering, we conclude that local quantum quenches are an ideal setting for studying their properties. We comment on implications of our results for cold atom experiments.
△ Less
Submitted 19 December, 2011; v1 submitted 14 December, 2011;
originally announced December 2011.
-
The QCD Phase Diagram with Effective Theories
Authors:
Ydalia Delgado,
Hans Gerd Evertz,
Christof Gattringer,
Daniel Göschl
Abstract:
We study two effective theories for QCD at non-zero temperature and finite chemical potential, using local Polyakov loops as the degrees of freedom. The sign problem is solved by exactly mapping the partition function to a sum over flux and monomer variables with only real and positive weights, making the two models accessible to Monte Carlo simulation techniques. We use generalized worm algorithm…
▽ More
We study two effective theories for QCD at non-zero temperature and finite chemical potential, using local Polyakov loops as the degrees of freedom. The sign problem is solved by exactly mapping the partition function to a sum over flux and monomer variables with only real and positive weights, making the two models accessible to Monte Carlo simulation techniques. We use generalized worm algorithms and a local Metropolis update to perform the simulations and determine the phase diagram as a function of the temperature and the chemical potential.
△ Less
Submitted 3 November, 2011;
originally announced November 2011.
-
Solving the Complex Phase Problem in a QCD Related Model
Authors:
Ydalia Delgado,
Hans Gerd Evertz,
Christof Gattringer
Abstract:
We discuss an effective theory for QCD at finite chemical potential and non-zero temperature, where QCD is reduced to its center degrees of freedom. The effective action can be mapped to a flux representation, where the complex phase problem is solved and the theory accessible to Monte Carlo techniques. In this work, we use a generalized Prokof'ev-Svistunov worm algorithm to perform the simulation…
▽ More
We discuss an effective theory for QCD at finite chemical potential and non-zero temperature, where QCD is reduced to its center degrees of freedom. The effective action can be mapped to a flux representation, where the complex phase problem is solved and the theory accessible to Monte Carlo techniques. In this work, we use a generalized Prokof'ev-Svistunov worm algorithm to perform the simulations and determine the phase diagram as a function of temperature, quark mass and chemical potential. It turns out that the transition is qualitatively as expected for QCD.
△ Less
Submitted 31 October, 2011;
originally announced October 2011.
-
Superconductivity and charge order of confined Fermi systems
Authors:
E. Assmann,
S. Chiesa,
G. G. Batrouni,
H. G. Evertz,
R. T. Scalettar
Abstract:
The low-temperature properties of the two-dimensional attractive Hubbard model are strongly influenced by the fermion density. Away from half-filling, there is a finite-temperature transition to a phase with s-wave pairing order. However, the critical temperature is suppressed to zero at half-filling, where long-range charge-density-wave order also appears, degenerate with superconductivity. This…
▽ More
The low-temperature properties of the two-dimensional attractive Hubbard model are strongly influenced by the fermion density. Away from half-filling, there is a finite-temperature transition to a phase with s-wave pairing order. However, the critical temperature is suppressed to zero at half-filling, where long-range charge-density-wave order also appears, degenerate with superconductivity. This paper presents Determinant Quantum Monte Carlo simulations of the attractive Hubbard model in the presence of a confining potential V which makes the fermion density ρ inhomogeneous across the lattice. Pair correlations are shown to be large at low temperatures in regions of the trapped system with incommensurate filling, and to exhibit a minimum as the local density ρ(i) passes through one fermion per site. In this ring of ρ=1, charge order is enhanced. A comparison is made between treating V within the local-density approximation (LDA) and in an ab initio manner. It is argued that certain sharp features of the LDA result at integer filling do not survive the proximity of doped sites. The critical temperature of confined systems of fixed characteristic density is estimated.
△ Less
Submitted 2 September, 2011; v1 submitted 31 August, 2011;
originally announced August 2011.
-
The QCD phase diagram according to the center group
Authors:
Ydalia Delgado Mercado,
Hans Gerd Evertz,
Christof Gattringer
Abstract:
We study an effective theory for QCD at finite temperature and density which contains the leading center symmetric and center symmetry breaking terms. The effective theory is studied in a flux representation where the complex phase problem is absent and the model becomes accessible to Monte Carlo techniques also at finite chemical potential. We simulate the system using a generalized Prokof'ev-Svi…
▽ More
We study an effective theory for QCD at finite temperature and density which contains the leading center symmetric and center symmetry breaking terms. The effective theory is studied in a flux representation where the complex phase problem is absent and the model becomes accessible to Monte Carlo techniques also at finite chemical potential. We simulate the system using a generalized Prokof'ev-Svistunov worm algorithm and compare the results to a low temperature expansion. The phase diagram is determined as a function of temperature, chemical potential and quark mass. Shape and quark mass dependence of the phase boundaries are as expected for QCD. The transition into the deconfined phase is smooth throughout, without any discontinuities or critical points.
△ Less
Submitted 12 May, 2011; v1 submitted 15 February, 2011;
originally announced February 2011.
-
The ALPS project release 2.0: Open source software for strongly correlated systems
Authors:
B. Bauer,
L. D. Carr,
H. G. Evertz,
A. Feiguin,
J. Freire,
S. Fuchs,
L. Gamper,
J. Gukelberger,
E. Gull,
S. Guertler,
A. Hehn,
R. Igarashi,
S. V. Isakov,
D. Koop,
P. N. Ma,
P. Mates,
H. Matsuo,
O. Parcollet,
G. Pawlowski,
J. D. Picon,
L. Pollet,
E. Santos,
V. W. Scarola,
U. Schollwöck,
C. Silva
, et al. (7 additional authors not shown)
Abstract:
We present release 2.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. The code development is centered on common XML and HDF5 data formats, libraries…
▽ More
We present release 2.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. The code development is centered on common XML and HDF5 data formats, libraries to simplify and speed up code development, common evaluation and plotting tools, and simulation programs. The programs enable non-experts to start carrying out serial or parallel numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), the density matrix renormalization group (DMRG) both in a static version and a dynamic time-evolving block decimation (TEBD) code, and quantum Monte Carlo solvers for dynamical mean field theory (DMFT). The ALPS libraries provide a powerful framework for programers to develop their own applications, which, for instance, greatly simplify the steps of porting a serial code onto a parallel, distributed memory machine. Major changes in release 2.0 include the use of HDF5 for binary data, evaluation tools in Python, support for the Windows operating system, the use of CMake as build system and binary installation packages for Mac OS X and Windows, and integration with the VisTrails workflow provenance tool. The software is available from our web server at http://alps.comp-phys.org/.
△ Less
Submitted 23 May, 2011; v1 submitted 13 January, 2011;
originally announced January 2011.
-
Effect of the exchange hole on the Gutzwiller approximation in one dimension
Authors:
Balazs Hetenyi,
Hans Gerd Evertz,
Wolfgang von der Linden
Abstract:
The Gutzwiller approximate solution to the Gutzwiller wavefunction yields exact results for the Gutzwiller wavefunction in the infinite dimensional limit. Implicit in the Gutzwiller approximation is an approximate local form of the fermion exchange hole. This approximate form is the same for all dimensions but is incorrect except in infinite dimensions. We implement the correct form for the exch…
▽ More
The Gutzwiller approximate solution to the Gutzwiller wavefunction yields exact results for the Gutzwiller wavefunction in the infinite dimensional limit. Implicit in the Gutzwiller approximation is an approximate local form of the fermion exchange hole. This approximate form is the same for all dimensions but is incorrect except in infinite dimensions. We implement the correct form for the exchange hole into the Gutzwiller approximation. We perform calculations on the one-dimensional Hubbard model at half-filling. They indicate that the implementation of the exchange hole already brings the Gutzwiller approximation into very close quantitative agreement with the results of the full Gutzwiller wavefunction. Metallicity as well as anti-ferromagnetism are recovered.
△ Less
Submitted 27 July, 2009; v1 submitted 16 June, 2009;
originally announced June 2009.
-
Excitation spectra of strongly correlated lattice bosons and polaritons
Authors:
Peter Pippan,
Hans Gerd Evertz,
Martin Hohenadler
Abstract:
Spectral properties of the Bose-Hubbard model and a recently proposed coupled-cavity model are studied by means of quantum Monte Carlo simulations in one dimension. Both models exhibit a quantum phase transition from a Mott insulator to a superfluid phase. The dynamic structure factor $S(k,ω)$ and the single-particle spectrum $A(k,ω)$ are calculated, focusing on the parameter region around the p…
▽ More
Spectral properties of the Bose-Hubbard model and a recently proposed coupled-cavity model are studied by means of quantum Monte Carlo simulations in one dimension. Both models exhibit a quantum phase transition from a Mott insulator to a superfluid phase. The dynamic structure factor $S(k,ω)$ and the single-particle spectrum $A(k,ω)$ are calculated, focusing on the parameter region around the phase transition from the Mott insulator with density one to the superfluid phase, where correlations are important. The strongly interacting nature of the superfluid phase manifests itself in terms of additional gapped modes in the spectra. Comparison is made to recent analytical work on the Bose-Hubbard model. Despite some subtle differences due to the polaritonic particles in the cavity model, the gross features are found to be very similar to the Bose-Hubbard case. For the polariton model, emergent particle-hole symmetry near the Mott lobe tip is demonstrated, and temperature and detuning effects are analyzed. A scaling analysis for the generic transition suggests mean field exponents, in accordance with field theory results.
△ Less
Submitted 9 September, 2009; v1 submitted 8 April, 2009;
originally announced April 2009.
-
Loop updates for variational and projector quantum Monte Carlo simulations in the valence-bond basis
Authors:
A. W. Sandvik,
H. G. Evertz
Abstract:
We show how efficient loop updates, originally developed for Monte Carlo simulations of quantum spin systems at finite temperature, can be combined with a ground-state projector scheme and variational calculations in the valence bond basis. The methods are formulated in a combined space of spin z-components and valence bonds. Compared to schemes formulated purely in the valence bond basis, the com…
▽ More
We show how efficient loop updates, originally developed for Monte Carlo simulations of quantum spin systems at finite temperature, can be combined with a ground-state projector scheme and variational calculations in the valence bond basis. The methods are formulated in a combined space of spin z-components and valence bonds. Compared to schemes formulated purely in the valence bond basis, the computational effort is reduced from up to O(N^2) to O(N) for variational calculations, where N is the system size, and from O(m^2) to O(m) for projector simulations, where m>> N is the projection power. These improvements enable access to ground states of significantly larger lattices than previously. We demonstrate the efficiency of the approach by calculating the sublattice magnetization M_s of the two-dimensional Heisenberg model to high precision, using systems with up to 256*256 spins. Extrapolating the results to the thermodynamic limit gives M_s=0.30743(1). We also discuss optimized variational amplitude-product states, which were used as trial states in the projector simulations, and compare results of projecting different types of trial states.
△ Less
Submitted 24 May, 2010; v1 submitted 4 July, 2008;
originally announced July 2008.
-
Parametrized canonical transformation for the Hubbard-model at arbitrary interaction strength
Authors:
Balazs Hetenyi,
Hans Gerd Evertz
Abstract:
The $t-J$ and Heisenberg models are truncated expansions of a canonically transformed Hubbard model coinciding with it at $U\to \infty$. We show that a modified canonical transformation applied to the Hubbard model leads to alternative models of the same form, but whose convergence properties with respect to the expansion are more favourable, resulting in a good description of the half-filled gr…
▽ More
The $t-J$ and Heisenberg models are truncated expansions of a canonically transformed Hubbard model coinciding with it at $U\to \infty$. We show that a modified canonical transformation applied to the Hubbard model leads to alternative models of the same form, but whose convergence properties with respect to the expansion are more favourable, resulting in a good description of the half-filled ground state even at $0<U\leq1$. We investigate the transformed Hamiltonian and observables for metallic and insulating variational wave-functions.
△ Less
Submitted 27 July, 2008; v1 submitted 24 June, 2008;
originally announced June 2008.
-
Efficient Matrix Product State Method for periodic boundary conditions
Authors:
Peter Pippan,
Steven R. White,
Hans Gerd Evertz
Abstract:
We introduce an efficient method to calculate the ground state of one-dimensional lattice models with periodic boundary conditions. The method works in the representation of Matrix Product States (MPS), related to the Density Matrix Renormalization Group (DMRG) method. It improves on a previous approach by Verstraete et al. We introduce a factorization procedure for long products of MPS matrices…
▽ More
We introduce an efficient method to calculate the ground state of one-dimensional lattice models with periodic boundary conditions. The method works in the representation of Matrix Product States (MPS), related to the Density Matrix Renormalization Group (DMRG) method. It improves on a previous approach by Verstraete et al. We introduce a factorization procedure for long products of MPS matrices, which reduces the computational effort from m^5 to m^3, where m is the matrix dimension, and m ~ 100 - 1000 in typical cases. We test the method on the S=1/2 and S=1 Heisenberg chains. It is also applicable to non-translationally invariant cases. The new method makes ground state calculations with periodic boundary conditions about as efficient as traditional DMRG calculations for systems with open boundaries.
△ Less
Submitted 16 February, 2010; v1 submitted 13 January, 2008;
originally announced January 2008.
-
Perfect Tempering
Authors:
M. Daghofer,
M. Konegger,
H. G. Evertz,
W. von der Linden
Abstract:
Multimodal structures in the sampling density (e.g. two competing phases) can be a serious problem for traditional Markov Chain Monte Carlo (MCMC), because correct sampling of the different structures can only be guaranteed for infinite sampling time. Samples may not decouple from the initial configuration for a long time and autocorrelation times may be hard to determine.
We analyze a suitabl…
▽ More
Multimodal structures in the sampling density (e.g. two competing phases) can be a serious problem for traditional Markov Chain Monte Carlo (MCMC), because correct sampling of the different structures can only be guaranteed for infinite sampling time. Samples may not decouple from the initial configuration for a long time and autocorrelation times may be hard to determine.
We analyze a suitable modification (C. J. Geyer and E. A. Thompson, J. Amer. Statist. Assoc., 90, 909, 1995) of the simulated tempering idea (E. Marinari and G. Parisi, Europhys. Lett. 19, 451, 1992), which has orders of magnitude smaller autocorrelation times for multimodal sampling densities and which samples all peaks of multimodal structures according to their weight. The method generates exact, i.e. uncorrelated, samples and thus gives access to reliable error estimates. Exact tempering is applicable to arbitrary (continuous or discreet) sampling densities and moreover presents a possibility to calculate integrals over the density (e.g. the partition function for the Boltzmann distribution), which are not accessible by usual MCMC.
△ Less
Submitted 22 May, 2006; v1 submitted 19 December, 2005;
originally announced December 2005.
-
Aspects of the FM Kondo Model: From Unbiased MC Simulations to Back-of-an-Envelope Explanations
Authors:
M. Daghofer,
W. Koller,
A. Pruell,
H. G. Evertz,
W. von der Linden
Abstract:
Effective models are derived from the ferromagnetic Kondo lattice model with classical corespins, which greatly reduce the numerical effort. Results for these models are presented. They indicate that double exchange gives the correct order of magnitude and the correct doping dependence of the Curie temperature. Furthermore, we find that the jump in the particle density previously interpreted as…
▽ More
Effective models are derived from the ferromagnetic Kondo lattice model with classical corespins, which greatly reduce the numerical effort. Results for these models are presented. They indicate that double exchange gives the correct order of magnitude and the correct doping dependence of the Curie temperature. Furthermore, we find that the jump in the particle density previously interpreted as phase separation is rather explained by ferromagnetic polarons.
△ Less
Submitted 19 December, 2005;
originally announced December 2005.
-
Creation and destruction of a spin gap in weakly coupled quarter-filled ladders
Authors:
B. Edegger,
H. G. Evertz,
R. M. Noack
Abstract:
We investigate weakly coupled quarter-filled ladders with model parameters relevant for NaV_2O_5 using density-matrix renormalization group calculations on an extended Hubbard model coupled to the lattice. NaV2O5 exhibits super-antiferroelectric charge order with a zigzag pattern on each ladder. We show that this order, with a periodicity of four ladders, causes a spin dimerization along the lad…
▽ More
We investigate weakly coupled quarter-filled ladders with model parameters relevant for NaV_2O_5 using density-matrix renormalization group calculations on an extended Hubbard model coupled to the lattice. NaV2O5 exhibits super-antiferroelectric charge order with a zigzag pattern on each ladder. We show that this order, with a periodicity of four ladders, causes a spin dimerization along the ladder and a corresponding spin gap of the same magnitude as that observed experimentally. The spin gap is destroyed again at large charge order due to a restructuring of the spins. An analysis of an effective spin model predicts a re-creation of the gap by inter-ladder singlets when the charge order increases further.
△ Less
Submitted 13 March, 2006; v1 submitted 12 October, 2005;
originally announced October 2005.
-
Charge order induced by electron-lattice interaction in NaV2O5
Authors:
B. Edegger,
H. G. Evertz,
R. M. Noack
Abstract:
We present Density Matrix Renormalization Group calculations of the ground-state properties of quarter-filled ladders including static electron-lattice coupling. Isolated ladders and two coupled ladders are considered, with model parameters obtained from band-structure calculations for $α^\prime$-NaV$_2$O$_5$. The relevant Holstein coupling to the lattice causes static out-of-plane lattice disto…
▽ More
We present Density Matrix Renormalization Group calculations of the ground-state properties of quarter-filled ladders including static electron-lattice coupling. Isolated ladders and two coupled ladders are considered, with model parameters obtained from band-structure calculations for $α^\prime$-NaV$_2$O$_5$. The relevant Holstein coupling to the lattice causes static out-of-plane lattice distortions, which appear concurrently with a charge-ordered state and which exhibit the same zigzag pattern observed in experiments. The inclusion of electron-lattice coupling drastically reduces the critical nearest-neighbor Coulomb repulsion $V_c$ needed to obtain the charge-ordered state. No spin gap is present in the ordered phase. The charge ordering is driven by the Coulomb repulsion and the electron-lattice interaction. With electron-lattice interaction, coupling two ladders has virtually no effect on $V_c$ or on the characteristics of the charge-ordered phase. At $V=0.46\eV$, a value consistent with previous estimates, the lattice distortion, charge gap, charge order parameter, and the effective spin coupling are in good agreement with experimental data for NaV$_2$O_5$.
△ Less
Submitted 30 May, 2005;
originally announced May 2005.
-
Ferromagnetic polarons in the one-dimensional ferromagnetic Kondo model with quantum mechanical S=3/2 core spins
Authors:
D. R. Neuber,
M. Daghofer,
R. M. Noack,
H. G. Evertz,
W. von der Linden
Abstract:
We present an extensive numerical study of the ferromagnetic Kondo lattice model with quantum mechanical S=3/2 core spins. We treat one orbital per site in one dimension using the density matrix renormalization group and include on-site Coulomb repulsion between the electrons. We examine parameters relevant to manganites, treating the range of low to intermediate doping, 0<x<0.5. In particular,…
▽ More
We present an extensive numerical study of the ferromagnetic Kondo lattice model with quantum mechanical S=3/2 core spins. We treat one orbital per site in one dimension using the density matrix renormalization group and include on-site Coulomb repulsion between the electrons. We examine parameters relevant to manganites, treating the range of low to intermediate doping, 0<x<0.5. In particular, we investigate whether quantum fluctuations favor phase separation over the ferromagnetic polarons observed in a model with classical core spins. We obtain very good agreement of the quantum model with previous results for the classical model, finding separated polarons which are repulsive at short distance for finite t_2g superexchange J'. Taking on-site Coulomb repulsion into account, we observe phase separation for small but finite superexchange J', while for larger J' polarons are favored in accordance with simple energy considerations previously applied to classical spins. We discuss the interpretation of compressibilities and present a phase diagram with respect to doping and the t_2g superexchange parameter J' with and without Coulomb repulsion.
△ Less
Submitted 2 December, 2005; v1 submitted 11 January, 2005;
originally announced January 2005.
-
Ferromagnetic Polarons, Phase Separation, Stripes and Polaron Lattice: The two- and three-dimensional ferromagnetic Kondo Model
Authors:
M. Daghofer,
W. Koller,
H. G. Evertz,
W. von der Linden
Abstract:
We investigate the two- and three-dimensional ferromagnetic Kondo lattice model by unbiased Monte Carlo simulations. A phase diagram for the two-dimensional model is presented, in which the stability of magnetic order and ferromagnetic polarons is examined with respect to the antiferromagnetic superexchange J' and temperature. The Monte Carlo simulations reveal that J' > 0.02 strengthens individ…
▽ More
We investigate the two- and three-dimensional ferromagnetic Kondo lattice model by unbiased Monte Carlo simulations. A phase diagram for the two-dimensional model is presented, in which the stability of magnetic order and ferromagnetic polarons is examined with respect to the antiferromagnetic superexchange J' and temperature. The Monte Carlo simulations reveal that J' > 0.02 strengthens individual polarons while small J' < 0.02 favors larger clusters and phase separation except for small doping. Lowering the temperature stabilizes ferromagnetic polarons for realistic J' > 0.01, while phase separation is only favored for very small J' < 0.01. Our Monte Carlo simulations show that low temperatures can lead to diagonal or vertical stripes depending on J'.
Simulations for three-dimensional systems yield ferromagnetic polarons, which form a `polaron lattice' at higher doping levels 0.2 < x < 0.23, when independent polarons do no longer fit into the system. No tendency to phase separation is observed in three dimensions.
△ Less
Submitted 12 October, 2004;
originally announced October 2004.
-
Single-particle spectral function of quarter-filled ladder systems
Authors:
M. Aichhorn,
E. Ya. Sherman,
H. G. Evertz
Abstract:
We study the single-particle properties of quarter-filled ladder systems such as sodium vanadate by means of a recently developed generalization of the variational cluster perturbation theory to extended Hubbard models. We find a homogeneous antiferromagnetic insulating phase for nearest-neighbor repulsions smaller than a critical value, without any metallic phase for small repulsions. Different…
▽ More
We study the single-particle properties of quarter-filled ladder systems such as sodium vanadate by means of a recently developed generalization of the variational cluster perturbation theory to extended Hubbard models. We find a homogeneous antiferromagnetic insulating phase for nearest-neighbor repulsions smaller than a critical value, without any metallic phase for small repulsions. Different from C-DMFT and LDA considerations, the inclusion of diagonal hopping within a ladder has little effect on the bonding bands, while flattening and shifting the antibonding bands. In the low-temperature charge-ordered phase, the spectrum depends on whether the ordering is driven by the Coulomb repulsion or by the coupling to a static lattice distortion. The small change of the experimentally observed gap upon charge ordering implies that the lattice coupling plays an important role in this ordering. Inter-ladder coupling is straightforward to include within our method. We show that it has only a minor effect on the spectral function. The numerically calculated spectra show good agreement with experimental angle-resolved photo-emission data.
△ Less
Submitted 20 October, 2005; v1 submitted 7 September, 2004;
originally announced September 2004.
-
Finite-temperature investigation of quarter filled ladder systems
Authors:
C. Gabriel,
E. Sherman,
T. C. Lang,
M. Aichhorn,
H. G. Evertz
Abstract:
We investigate charge ordering in a quarter-filled ladder at finite temperature by determinantal Quantum Monte Carlo. The sign problem is moderate in a wide range of model parameters relevant for NaV2O5. The charge order parameter exhibits a crossover as a function of inverse temperature on finite systems. Above a critical nearest neighbor Coulomb repulsion, the correlation length grows exponent…
▽ More
We investigate charge ordering in a quarter-filled ladder at finite temperature by determinantal Quantum Monte Carlo. The sign problem is moderate in a wide range of model parameters relevant for NaV2O5. The charge order parameter exhibits a crossover as a function of inverse temperature on finite systems. Above a critical nearest neighbor Coulomb repulsion, the correlation length grows exponentially with inverse temperature, indicative of the ordered phase at T=0. We find a clear single-particle gap manifesting itself in a flat n(μ) dependence at large nearest neighbor Coulomb repulsion.
△ Less
Submitted 29 June, 2004;
originally announced June 2004.
-
Single particle spectrum of the flux phase in the FM Kondo Model
Authors:
M. Daghofer,
W. Koller,
W. von der Linden,
H. G. Evertz
Abstract:
We investigate the 2D ferromagnetic Kondo lattice model for manganites with classical corespins at Hund's rule coupling J_H=6, with antiferromagnetic superexchange 0.03 < J' < 0.05. We employ canonical and grand canonical unbiased Monte Carlo simulations and find paramagnetism, weak ferromagnetism and the Flux phase, depending on doping and on J'. The observed single particle spectrum in the flu…
▽ More
We investigate the 2D ferromagnetic Kondo lattice model for manganites with classical corespins at Hund's rule coupling J_H=6, with antiferromagnetic superexchange 0.03 < J' < 0.05. We employ canonical and grand canonical unbiased Monte Carlo simulations and find paramagnetism, weak ferromagnetism and the Flux phase, depending on doping and on J'. The observed single particle spectrum in the flux phase differs from the idealized infinite lattice case, but agrees well with an idealized finite lattice case with thermal fluctuations.
△ Less
Submitted 25 June, 2004;
originally announced June 2004.
-
Optical properties, electron-phonon coupling, and Raman scattering of vanadium ladder compounds
Authors:
J. Spitaler,
E. Ya. Sherman,
H. G. Evertz,
C. Ambrosch-Draxl
Abstract:
The electronic structure of two V-based ladder compounds, the quarter-filled NaV$_2$O$_5$ in the symmetric phase and the iso-structural half-filled CaV$_2$O$_5$ is investigated by ab initio calculations. Based on the bandstructure we determine the dielectric tensor $ε(ω)$ of these systems in a wide energy range. The frequencies and eigenvectors of the fully symmetric A$_{g}$ phonon modes and the…
▽ More
The electronic structure of two V-based ladder compounds, the quarter-filled NaV$_2$O$_5$ in the symmetric phase and the iso-structural half-filled CaV$_2$O$_5$ is investigated by ab initio calculations. Based on the bandstructure we determine the dielectric tensor $ε(ω)$ of these systems in a wide energy range. The frequencies and eigenvectors of the fully symmetric A$_{g}$ phonon modes and the corresponding electron-phonon and spin-phonon coupling parameters are also calculated from first-principles. We determine the Raman scattering intensities of the A$_g$ phonon modes as a function of polarization and frequency of the exciting light.
All results, i.e. shape and magnitude of the dielectric function, phonon frequencies and Raman intensities show very good agreement with available experimental data.
△ Less
Submitted 20 September, 2004; v1 submitted 3 April, 2004;
originally announced April 2004.
-
Charge ordering in extended Hubbard models: Variational cluster approach
Authors:
M. Aichhorn,
H. G. Evertz,
W. von der Linden,
M. Potthoff
Abstract:
We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalisation of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For…
▽ More
We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalisation of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For one dimension, we show that quantum Monte Carlo and density-matrix renormalization-group results can be reproduced with very good accuracy. Moreover we apply the method to the two-dimensional extended Hubbard model on a square lattice. In contrast to the one-dimensional case, a first order phase transition between spin density wave phase and charge density wave phase is found as function of the nearest-neighbor interaction at onsite interactions U>=3t. The single-particle spectral function is calculated for both the one-dimensional and the two-dimensional system.
△ Less
Submitted 24 February, 2004;
originally announced February 2004.
-
Holstein polaron in two and three dimensions by quantum Monte Carlo
Authors:
Martin Hohenadler,
Hans Gerd Evertz,
Wolfgang von der Linden
Abstract:
A recently developed quantum Monte Carlo approach to the Holstein model with one electron [PRB 69, 024301 (2004)] is extended to two and three dimensional lattices. A moderate sign problem occurs, which is found to diminish with increasing system size in all dimensions, and not to affect simulations significantly. We present an extensive study of the influence of temperature, system size, dimens…
▽ More
A recently developed quantum Monte Carlo approach to the Holstein model with one electron [PRB 69, 024301 (2004)] is extended to two and three dimensional lattices. A moderate sign problem occurs, which is found to diminish with increasing system size in all dimensions, and not to affect simulations significantly. We present an extensive study of the influence of temperature, system size, dimensionality and model parameters on the small-polaron cross over. Results are extrapolated to remove the error due to the Trotter discretization, which significantly improves the accuracy. Comparison with existing work and other quantum Monte Carlo methods is made. The method can be extended to the many-electron case.
△ Less
Submitted 30 January, 2004;
originally announced January 2004.
-
Charge ordering in quarter-filled ladder systems coupled to the lattice
Authors:
M. Aichhorn,
M. Hohenadler,
E. Ya. Sherman,
J. Spitaler,
C. Ambrosch-Draxl,
H. G. Evertz
Abstract:
We investigate charge ordering in the presence of electron-phonon coupling for quarter-filled ladder systems by using Exact Diagonalization. As an example we consider NaV2O5 using model parameters obtained from first-principles band-structure calculations. The relevant Holstein coupling to the lattice considerably reduces the critical value of the nearest-neighbor Coulomb repulsion at which form…
▽ More
We investigate charge ordering in the presence of electron-phonon coupling for quarter-filled ladder systems by using Exact Diagonalization. As an example we consider NaV2O5 using model parameters obtained from first-principles band-structure calculations. The relevant Holstein coupling to the lattice considerably reduces the critical value of the nearest-neighbor Coulomb repulsion at which formation of the zig-zag charge-ordered state occurs, which is then accompanied by a static lattice distortion. Energy and length of a kink-like excitation on the background of the distorted lattice are calculated. Spin and charge spectra on ladders with and without static distortion are obtained, and the charge gap and the effective spin-spin exchange parameter J are extracted. J agrees well with experimental results. Analysis of the dynamical Holstein model, restricted to a small number of phonons, shows that low frequency lattice vibrations increase the charge order, accompanied by dynamically produced zig-zag lattice distortions.
△ Less
Submitted 26 May, 2004; v1 submitted 19 December, 2003;
originally announced December 2003.
-
Polaronic Aspects of the two-dimensional Ferromagnetic Kondo Model
Authors:
M. Daghofer,
W. Koller,
H. G. Evertz,
W. von der Linden
Abstract:
The 2D ferromagnetic Kondo model with classical corespins is studied for a hole doping up to x = 0.125 via unbiased Monte Carlo (MC) simulations. A canonical algorithm for finite temperatures is developed.
We show that with realistic parameters for the manganites and at low temperatures, the double-exchange mechanism does not lead to phase separation on a two-dimensional lattice but rather sta…
▽ More
The 2D ferromagnetic Kondo model with classical corespins is studied for a hole doping up to x = 0.125 via unbiased Monte Carlo (MC) simulations. A canonical algorithm for finite temperatures is developed.
We show that with realistic parameters for the manganites and at low temperatures, the double-exchange mechanism does not lead to phase separation on a two-dimensional lattice but rather stabilizes individual ferromagnetic polarons. A detailed analysis of unbiased MC results reveals that the polarons can be treated as independent particles.
It is found that a simple polaron model perfectly describes the physics of the FM Kondo model. The ferromagnetic polaron picture provides an obvious explanation for the pseudogap in the one-particle spectral function $A_k(ω)$ observed in the FM Kondo model.
△ Less
Submitted 5 August, 2004; v1 submitted 16 December, 2003;
originally announced December 2003.
-
Quantum Monte Carlo and variational approaches to the Holstein model
Authors:
Martin Hohenadler,
Hans Gerd Evertz,
Wolfgang von der Linden
Abstract:
Based on the canonical Lang-Firsov transformation of the Hamiltonian we develop a very efficient quantum Monte Carlo algorithm for the Holstein model with one electron. Separation of the fermionic degrees of freedom by a reweighting of the probability distribution leads to a dramatic reduction in computational effort. A principal component representation of the phonon degrees of freedom allows t…
▽ More
Based on the canonical Lang-Firsov transformation of the Hamiltonian we develop a very efficient quantum Monte Carlo algorithm for the Holstein model with one electron. Separation of the fermionic degrees of freedom by a reweighting of the probability distribution leads to a dramatic reduction in computational effort. A principal component representation of the phonon degrees of freedom allows to sample completely uncorrelated phonon configurations. The combination of these elements enables us to perform efficient simulations for a wide range of temperature, phonon frequency and electron-phonon coupling on clusters large enough to avoid finite-size effects. The algorithm is tested in one dimension and the data are compared with exact-diagonalization results and with existing work. Moreover, the ideas presented here can also be applied to the many-electron case. In the one-electron case considered here, the physics of the Holstein model can be described by a simple variational approach.
△ Less
Submitted 24 January, 2004; v1 submitted 16 May, 2003;
originally announced May 2003.
-
Magnetic Polarons in the 1D FM Kondo Model
Authors:
Winfried Koller,
Alexander Prüll,
Hans Gerd Evertz,
Wolfgang von der Linden
Abstract:
The ferromagnetic Kondo model with classical corespins is studied via unbiased Monte-Carlo simulations. We show that with realistic parameters for the manganites and at low temperatures, the double-exchange mechanism does not lead to phase separation in one-dimensional chains but rather stabilizes individual ferromagnetic polarons. Within the ferromagnetic polaron picture, the pseudogap in the o…
▽ More
The ferromagnetic Kondo model with classical corespins is studied via unbiased Monte-Carlo simulations. We show that with realistic parameters for the manganites and at low temperatures, the double-exchange mechanism does not lead to phase separation in one-dimensional chains but rather stabilizes individual ferromagnetic polarons. Within the ferromagnetic polaron picture, the pseudogap in the one-particle spectral function A_k(ω) can easily be explained. Ferromagnetic polarons also clear up a seeming failure of the double-exchange mechanism in explaining the comparable bandwidths in the ferromagnetic and paramagnetic phase. For our analysis, we extend a simplified model, the finite temperature uniform hopping approach (UHA), to include polarons. It can easily be evaluated numerically and provides a simple quantitative understanding of the physical features of the ferromagnetic Kondo model.
△ Less
Submitted 23 January, 2003; v1 submitted 20 January, 2003;
originally announced January 2003.
-
Uniform hopping approach to the FM Kondo Model at finite temperature
Authors:
Winfried Koller,
Alexander Pruell,
Hans Gerd Evertz,
Wolfgang von der Linden
Abstract:
We study the ferromagnetic Kondo model with classical corespins via unbiased Monte-Carlo simulations and derive a simplified model for the treatment of the corespins at any temperature. Our simplified model captures the main aspects of the Kondo model and can easily be evaluated both numerically and analytically. It provides a better qualitative understanding of the physical features of the Kond…
▽ More
We study the ferromagnetic Kondo model with classical corespins via unbiased Monte-Carlo simulations and derive a simplified model for the treatment of the corespins at any temperature. Our simplified model captures the main aspects of the Kondo model and can easily be evaluated both numerically and analytically. It provides a better qualitative understanding of the physical features of the Kondo model and rationalizes the Monte-Carlo results, including the spectral density A_k(omega) of a 1D chain with nearest neighbor Coulomb repulsion. By calculating the specific heat and the susceptibility of systems up to size 16^3, we determine the Curie temperature of the 3D one-orbital double-exchange model, which agrees with experimental values.
△ Less
Submitted 29 November, 2002; v1 submitted 28 November, 2002;
originally announced November 2002.
-
Low Temperature Lanczos Method
Authors:
Markus Aichhorn,
Maria Daghofer,
Hans Gerd Evertz,
Wolfgang von der Linden
Abstract:
We present a modified finite temperature Lanczos method for the evaluation of dynamical and static quantities of strongly correlated electron systems that complements the finite temperature method (FTLM) introduced by Jaklic and Prelovsek for low temperatures. Together they allow accurate calculations at any temperature with moderate effort. As an example we calculate the static spin correlation…
▽ More
We present a modified finite temperature Lanczos method for the evaluation of dynamical and static quantities of strongly correlated electron systems that complements the finite temperature method (FTLM) introduced by Jaklic and Prelovsek for low temperatures. Together they allow accurate calculations at any temperature with moderate effort. As an example we calculate the static spin correlation function and the regular part of the optical conductivity of the one dimensional Hubbard model at half-filling and show in detail the connection between the ground state and finite temperature method. By using Cluster Perturbation Theory (CPT), the finite temperature spectral function is extended to the infinite system, clearly exhibiting the effects of spin-charge separation.
△ Less
Submitted 8 November, 2002;
originally announced November 2002.
-
Effective spinless fermions in the strong coupling Kondo model
Authors:
Winfried Koller,
Alexander Pruell,
Hans Gerd Evertz,
Wolfgang von der Linden
Abstract:
Starting from the two-orbital Kondo-lattice model with classical t_2g spins, an effective spinless fermion model is derived for strong Hund coupling J_H with a projection technique. The model is studied by Monte Carlo simulations and analytically using a uniform hopping approximation. The results for the spinless fermion model are in remarkable agreement with those of the original Kondo-lattice…
▽ More
Starting from the two-orbital Kondo-lattice model with classical t_2g spins, an effective spinless fermion model is derived for strong Hund coupling J_H with a projection technique. The model is studied by Monte Carlo simulations and analytically using a uniform hopping approximation. The results for the spinless fermion model are in remarkable agreement with those of the original Kondo-lattice model, independent of the carrier concentration, and even for moderate Hund coupling J_H. Phase separation, the phase diagram in uniform hopping approximation, as well as spectral properties including the formation of a pseudo-gap are discussed for both the Kondo-lattice and the effective spinless fermion model in one and three dimensions.
△ Less
Submitted 21 November, 2002; v1 submitted 5 June, 2002;
originally announced June 2002.
-
Diamagnetic Properties of Doped Antiferromagnets
Authors:
D. Veberic,
P. Prelovsek,
H. G. Evertz
Abstract:
Finite-temperature diamagnetic properties of doped antiferromagnets as modeled by the two-dimensional t-J model were investigated by numerical studies of small model systems. Two numerical methods were used: the worldline quantum Monte Carlo method with a loop cluster algorithm (QMC) and the finite-temperature Lanczos method (FTLM), yielding consistent results. The diamagnetic susceptibility int…
▽ More
Finite-temperature diamagnetic properties of doped antiferromagnets as modeled by the two-dimensional t-J model were investigated by numerical studies of small model systems. Two numerical methods were used: the worldline quantum Monte Carlo method with a loop cluster algorithm (QMC) and the finite-temperature Lanczos method (FTLM), yielding consistent results. The diamagnetic susceptibility introduced by coupling of the magnetic field to the orbital current reveals an anomalous temperature dependence, changing character from diamagnetic to paramagnetic at intermediate temperatures.
△ Less
Submitted 25 March, 2002;
originally announced March 2002.