Source author record

Mona Berciu

Mona Berciu 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

49works
9topics
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

49 published item(s)

preprint2022arXiv

Dressing due to correlations strongly reduces the effect of electron-phonon coupling

We investigate the difference between the coupling of a bare carrier to phonons versus the coupling of a correlations-dressed quasiparticle to phonons, and show that latter may be weak even if the former is strong. Specifically, we analyze the effect of the hole-phonon coupling on the dispersion of the quasiparticle that forms when a single hole is doped into a cuprate layer. To model this, we start from the three-band Emery model supplemented by the Peierls modulation of the $p$-$d$ and $p$-$p$ hoppings due to the motion of O ions. We then project onto the strongly correlated $U_{dd}\rightarrow \infty$ limit where charge fluctuations are frozen on the Cu site. The resulting effective Hamiltonian describes the motion of a doped hole on the O sublattice, and its interactions with Cu spins and O phonons. We show that even though the hole-phonon coupling is moderate to strong, it leads to only a very minor increase of the quasiparticle's effective mass as compared to its mass in the absence of coupling to phonons, consistent with a weak coupling to phonons of the correlations-dressed quasiparticle. We explain the reasons for this suppression, revealing why it is expected to happen in any systems with strong correlations.

preprint2022arXiv

Role of long-range coupling on the properties of single polarons in models with dual electron-phonon couplings

We use the Variational Exact Diagonalization to investigate the single polaron properties for four different dual models, combining a short-range off-diagonal (Peierls) plus a longer-range diagonal (Holstein or breathing-mode) coupling. This allows us to investigate the sensitivity of various polaron properties both to the range of the diagonal coupling and to the specific diagonal coupling chosen. We find strong sensitivity to the range for all duals models as the adiabatic limit is approached, however considerable sensitivity is observed for some quantities even for large phonon frequencies. Also, the strong dependence of the results on the specific form of the diagonal coupling is observed everywhere in the parameter space. Taken together, these results suggest that careful consideration must be given to the specific coupling and its proper range when quantitative comparisons with experiments are sought.

preprint2022arXiv

Theory of dispersive optical phonons in resonant inelastic x-ray scattering experiments

The community currently lacks a complete understanding of how resonant inelastic x-ray scattering (RIXS) experiments probe the electron-phonon ($e$-ph) interaction in solids. For example, most theoretical models of this process have focused on dispersionless Einstein phonons. Using a recently developed momentum average (MA) variational approximation for computing RIXS spectra of band insulators, we examine the influence of both electron and phonon dispersion in the intermediate state of the scattering process. We find that the inclusion of either, and their mutual interplay, introduces significant momentum variations in the RIXS intensity, even for momentum-independent electron-phonon coupling. The phonon dispersion also induces nontrivial changes in the excitation line shapes, which can have a quantitative impact on the data analysis. These results highlight the considerable challenges of interpreting RIXS data in actual materials.

preprint2020arXiv

Peierls versus Holstein models for describing electron-phonon coupling in perovskites

We use the Momentum Average approximation together with perturbative approaches, in the appropriate limits, to study the single polaron physics on a perovskite lattice inspired by BaBiO3. We investigate electron-phonon coupling of the Peierls type whereby the motion of ions modulates the values of the hopping integrals between sites, and show that it cannot be mapped onto the simpler one-band Holstein model in the whole parameter space. This is because the dispersion of the Peierls polaron has sharp transitions where the ground-state momentum jumps between high-symmetry points in the Brillouin zone, whereas the Holstein polaron always has the same ground-state momentum. These results imply that careful consideration is required to choose the appropriate model for carrier-lattice coupling in such complex lattices.

preprint2020arXiv

The critical nature of the Ni spin state in doped NdNiO$_2$

Superconductivity with $T_c \approx 15K$ was recently found in doped NdNiO$_2$. The Ni$^{1+}$O$_2$ layers are expected to be Mott insulators so hole doping should produce Ni$^{2+}$ with $S=1$, incompatible with robust superconductivity. We show that the NiO$_2$ layers fall inside a ``critical'' region where the large $pd$ hybridization favors a singlet $^1\!A_1$ hole-doped state like in CuO$_2$. However, we find that the superexchange is about one order smaller than in cuprates, thus a magnon ``glue'' is very unlikely and another mechanism needs to be found.

preprint2019arXiv

Effect of the electron-lattice coupling on the charge and magnetic order in rare-earth nickelates

We investigate the impact of electron-lattice coupling on the stability of various magnetic orders in rare-earth nickelates. We use the Hartree-Fock approximation, at zero temperature, to study an effective, two-band model with correlations characterized by a Hubbard $U$ and a Hund's $J$. This is coupled to breathing-mode distortions of the octahedral oxygen cages, described semi-classically, with a Holstein term. We analyze the effect of the various parameters on the resulting phase diagram, in particular on the charge disproportionation and on the magnetic order. We confirm that the coupling to the lattice cooperates with Hund's coupling and thus encourages charge disproportionation. We also find that it favors the fully disproportionated, 4-site periodic magnetic order of type $\Uparrow 0 \Downarrow 0$. Other convergent magnetic phases, such as the collinear $\uparrow\uparrow\downarrow\downarrow$ and non-collinear $\uparrow\rightarrow\downarrow\leftarrow$ states, do not couple to the lattice because of their lack of charge disproportionation. Novel phases, e.g. with charge disproportionation but no magnetic order, are also found to be stabilized in specific conditions.

preprint2019arXiv

Relevance of Cu-3d multiplet structure in models of high Tc cuprates

We revisit the problem of the spectra of two holes in a CuO$_{2}$ layer, modeled as a Cu-d$^{8}$ impurity with full multiplet structure coupled to a full O-2p band as an approximation to the local electronic structure of a hole doped cuprate. Unlike previous studies that treated the O band as a featureless bath, we describe it with a realistic tight binding model. While our results are in qualitative agreement with previous work, we find considerable quantitative changes when using the proper O-2p band structure. We also find (i) that only the ligand O-2p orbitals play an essential role, within this impurity model; (ii) that the three-orbital Emery model provides an accurate description for the subspace with $^{1}\!A_1$ symmetry, which includes the ground-state in the relevant region of the phase diagram; (iii) that this ground-state has only $\sim 50\%$ overlap with a Zhang-Rice singlet; (iv) that there are other low-energy states, in subspaces with different symmetries, that are absent from the three-orbital Emery model and its one-band descendants. These states play an important role in describing the elementary excitations of doped cuprates.

preprint2016arXiv

Discontinuous polaron transition in a two-band model

We present exact diagonalization and momentum average approximation (MA) results for the single polaron properties of a one-dimensional two-band model with phonon-modulated hopping. At strong electron-phonon coupling, we find a novel type of sharp transition, where the polaron ground state momentum jumps discontinuously from $k=π$ to $k=0$. The nature and origin of this transition is investigated and compared to that of the Su-Schrieffer-Heeger (SSH) model, where a sharp but smooth transition was previously reported. We argue that such discontinuous transitions are a consequence of the multi-band nature of the model, and are unlikely to be observed in one-band models. We also show that MA describes qualitatively and even quantitatively accurately this polaron and its transition. Given its computationally efficient generalization to higher dimensions, MA thus promises to allow for accurate studies of electron-phonon coupling in multi-band models in higher dimensions.

preprint2016arXiv

Green's function variational approach to orbital polarons in KCuF3

We develop an $e_g$ orbital, $t$-$J$-like model of a single charge doped into a two-dimensional plane with ferromagnetic spin order and alternating orbital order, and present its solution by Green's functions in the variational approximation framework. The model is designed to represent the orbital physics within ferromagnetic $(a,b)$ planes of KCuF$_3$ and K$_2$CuF$_4$. The variational approximation (VA) relies on the systematic generation of equations of motion for the Green's function, taking into account the real-space constraints coming from the exclusion of doubly occupied sites. This method is compared to the firmly established self-consistent Born approximation, and to the variational cluster approximation (VCA) which relies on the itinerant regime of the model. We find that the present variational approximation captures the essential aspects of the spectral weight distribution of the coherent quasiparticle state and gives a result similar to the VCA, while also reproducing well the momentum dependence of the spectral moments. In contrast, the spectral function obtained within the self-consistent Born approximation is more incoherent and its quasiparticle is heavier, at strong effective couplings, than observed with VCA and VA.

preprint2016arXiv

Local environment effects on a charge carrier injected into an Ising chain

We present numerically exact results for the spectral function of a single charge carrier that is injected into an Ising chain at finite temperature $T$. Both ferromagnetic and antiferromagnetic coupling between the Ising spins are considered. The interaction between the carrier and the Ising spins is assumed to be on-site and of Ising type, as well. We find that the carrier's spectral function exhibits a distinctive fine structure of resonances that are due to the temporary entrapping of the carrier inside small magnetic domains. The connection to models of disordered binary alloys where similar effects can occur is explained. We use these results to construct an accurate (quasi)analytic approximation for low and medium-$T$ spectral functions.

preprint2016arXiv

Low-temperature evolution of the spectral weight of a spin-up carrier moving in a ferromagnetic background

We derive the lowest-temperature correction to the self-energy of a spin-up particle injected in a ferromagnetic background. The background is modeled with both Heisenberg and Ising Hamiltonians so that differences due to gapless vs. gapped magnons can be understood. Beside the expected thermal broadening of the quasiparticle peak as it becomes a resonance inside a continuum, we also find that spectral weight is transferred to regions lying outside this continuum. We explain the origin of this spectral weight transfer and its low-temperature evolution.

preprint2016arXiv

Non-Zhang-Rice singlet character of the first ionization state of T-CuO

We argue that tetragonal CuO (T-CuO) has the potential to finally settle long-standing modelling issues for cuprate physics. We compare the one-hole quasiparticle (qp) dispersion of T-CuO to that of cuprates, in the framework of the strongly-correlated ($U_{dd}\rightarrow \infty$) limit of the three-band Emery model. Unlike in CuO$_2$, magnetic frustration in T-CuO breaks the $C_4$ rotational symmetry and leads to strong deviations from the Zhang-Rice singlet picture in parts of the reciprocal space. Our results are consistent with angle-resolved photoemission spectroscopy data but in sharp contradiction to those of a one-band model previously suggested for them. These differences identify T-CuO as an ideal material to test a variety of scenarios proposed for explaining cuprate phenomenology.

preprint2016arXiv

On the equivalence of models with similar low-energy quasiparticles

We use a Metropolis algorithm to calculate the finite temperature spectral weight of three related models that have identical quasiparticles at $T=0$, if the exchange favors the appearance of a ferromagnetic background. The low-energy behavior of two of the models remains equivalent at finite temperature, however that of the third does not because its low-energy behavior is controlled by rare events due to thermal fluctuations, which transfer spectral weight well below the $T=0$ quasiparticle peaks and generate a pseudogap-like phenomenology. Our results demonstrate that having $T=0$ spectra with similar quasiparticles is not a sufficient condition to ensure that two models are equivalent, {\em i.e.} that their low-energy properties are similar. We also argue that the pseudogap-like phenomenology is quite generic for models of $t$-$J$ type, appearing in any dimension and for carriers injected into both ferromagnetic and antiferromagnetic backgrounds.

preprint2016arXiv

The Green function variational approximation: Significance of physical constraints

We present a calculation of the spectral properties of a single charge doped at a Cu($3d$) site of the Cu-F plane in KCuF$_{3}$. The problem is treated by generating the equations of motion for the Green's function by means of subsequent Dyson expansions and solving the resulting set of equations. This method, dubbed the variational approximation, is both very dependable and flexible, since it is a systematic expansion with precise control over elementary physical processes. It allows for deep insight into the underlying physics of polaron formation as well as for inclusion of many physical constraints, such as excluding crossing diagrams and double occupation constraint, which are not included in the Self-Consistent Born Approximation. Here we examine the role and importance of such constraints by analyzing various spectral functions obtained in second order VA.

preprint2015arXiv

Differences in the quasiparticle dynamics for one-band and three-band cuprate models

We study the quasiparticles of the one-band $t$-$J$ and $t$-$t'$-$t"$-$J$ models using a variational approximation that includes spin fluctuations in the vicinity of the hole. We explain why the spin fluctuations and the longer range hopping have complementary contributions to the quasiparticle dynamics, and thus why both are essential to obtain a dispersion in agreement with that measured experimentally. This is very different from the three-band Emery model in the strongly-correlated limit, where the same variational approximation shows that spin fluctuations have a minor effect on the quasiparticle dynamics. This difference proves that these one-band and three-band models describe qualitatively different quasiparticles, and therefore they they cannot both be suitable to describe the physics of underdoped cuprates.

preprint2014arXiv

Large electropositive cations as surfactants for the growth of polar epitaxial films

Using density functional theory (DFT) we demonstrate that the adsorption of large cations such as potassium or cesium facilitates the epitaxial growth of polar LaAlO$_3$ (LAO) on SrTiO$_3$ (STO). The low ionization potential of K favors efficient electron transfer to the STO conduction band and results in a 2D electron gas which exactly compensates for the diverging potential with increasing layer thickness. For large cations like K or Cs, DFT total energy considerations show that they remain adsorbed on the LAO surface and do not enter substitutionally into LAO. These results suggest a novel scheme for growing clean LAO/STO interface systems, and polar systems in general, by performing the growth process in the presence of large, low ionization potential alkali metal ions.

preprint2014arXiv

Phonon-assisted carrier motion on the Wannier-Stark ladder

It is well known that at zero temperature and in the absence of electron-phonon coupling, the presence of an electric field leads to localization of carriers residing in a single band of finite bandwidth. We implement the Self-Consistent Born Approximation (SCBA) to study the effect of weak electron-phonon coupling on the motion of a carrier in a biased system. At moderate and strong electron-phonon coupling we supplement the SCBA, describing the string of phonons left behind by the carrier, with the Momentum Average (MA) approximation to describe the phonon cloud that accompanies the resulting polaron. We find that coupling to the lattice delocalizes the carrier, as expected, although long-lived resonances resulting from the Wannier-Stark states of the polaron may appear in the spectrum in certain regions of the parameter space. The approach we propose here can also be used to implement and check the validity of simple variational approximations.

preprint2014arXiv

Single polaron properties for double-well electron-phonon coupling

We show that in crystals where light ions are symmetrically intercalated between heavy ions, the electron-phonon coupling for carriers located at the light sites cannot be described by a Holstein model. We introduce the double-well electron-phonon coupling model to describe the most interesting parameter regime in such systems, and study it in the single carrier limit using the momentum average approximation. For sufficiently strong coupling, a small polaron with a robust phonon cloud appears at low energies. While some of its properties are similar to those of a Holstein polaron, we highlight some crucial differences. These prove that the physics of the double-well electron-phonon coupling model cannot be reproduced with a linear Holstein model.

preprint2014arXiv

Strongly bound yet light bipolarons for double-well electron-phonon coupling

We use the Momentum Average approximation (MA) to study the ground-state properties of strongly bound bipolarons in the double-well electron-phonon (el-ph) coupling model, which describes certain intercalated lattices where the linear term in the el-ph coupling vanishes due tosymmetry. We show that this model predicts the existence of strongly bound yet lightweight bipolarons in some regions of the parameter space. This provides a novel mechanism for the appearance of such bipolarons, in addition to long-range el-ph coupling and special lattice geometries.

preprint2014arXiv

The dynamics of a doped hole in cuprates is not controlled by spin fluctuations

Twenty seven years after the discovery of high-temperature superconductivity \cite{BedMu}, consensus on its theoretical explanation is still absent. To a good extent, this is due to the difficulty of studying strongly correlated systems near half-filling, needed to understand the behaviour of one or few holes doped into a CuO$_2$ layer. To simplify this task it is customary to replace three-band models \cite{Emery} describing the doping holes as entering the O $2p$ orbitals of these charge-transfer insulators \cite{ZSA} with much simpler one-band Hubbard or $tJ$ models \cite{rev1,rev2}. Here we challenge this approach, showing that not only is the dynamics of a doped hole easier to understand in models that explicitly include the O orbitals, but also that our solution contradicts the long-held belief that the quantum spin fluctuations of the antiferromagnetic (AFM) background play a key role in determining this dynamics. Indeed, we show that the correct, experimentally observed dispersion is generically obtained for a hole moving on the O sublattice, and coupled to a Néel lattice of spins without spin fluctuations. This marks a significant conceptual change in our understanding of the relevant phenomenology and opens the way to studying few-holes dynamics without finite-size effect issues \cite{BayoB}, to understand the actual strength of the "magnetic glue".

preprint2013arXiv

A variational approach for calculating Auger electron spectra: going beyond the impurity approximation

We propose a novel variational method to calculate the two-hole propagators relevant for Auger spectroscopy in transition metal oxides. This method can be thought of as an intermediary step between the full solution (which is difficult to generalize to systems with partially filled bands) and the impurity approximation. Like the former, our solution has full translational invariance, and like the latter, it can be generalized to certain types of systems with partially filled bands. Here we compare both our variational approximation and the impurity approximation against the exact solution for a simple one-dimensional model with filled bands. We show that when the energies of the eigenstates residing primarily on the transition metal ions do not overlap with those of the eigenstates residing primarily on Oxygen ions, both approximations are valid but the variational approach is superior.

preprint2013arXiv

Binding carriers to a non-magnetic impurity in a two-dimensional square Ising antiferromagnet

A hole in a two-dimensional Ising antiferromagnet was believed to be infinitely heavy due to the string of wrongly oriented spins it creates as it moves, which should trap it near its original location. Trugman showed that, in fact, the hole acquires a finite effective mass due to contributions from so-called {\it Trugman loop} processes, where the hole goes nearly twice around closed loops, first creating and then removing wrongly-oriented spins, and ending up at a different lattice site. This generates effective second- and third-nearest-neighbor hopping terms which keep the quasiparticle on the sublattice it was created on. Here, we investigate the trapping of the quasiparticle near a single attractive non-magnetic impurity placed at one lattice site. We consider the two cases with the quasiparticle and impurity being on the same versus on different sublattices. The main result is that even though the quasiparticle can not see the bare disorder in the latter case, the coupling to magnons generates an effective renormalized disorder on its own sublattice which is strong enough to lead to bound states, which however have a very different spectrum than when the quasiparticle and impurity are on the same sublattice.

preprint2013arXiv

Charge disproportionation without charge transfer in the rare-earth nickelates as a possible mechanism for the metal-insulator transition

We study a model for the metal-insulator (MI) transition in the rare-earth nickelates RNiO$_3$, based upon a negative charge transfer energy and coupling to a rock-salt like lattice distortion of the NiO$_6$ octahedra. Using exact diagonalization and the Hartree-Fock approximation we demonstrate that electrons couple strongly to these distortions. For small distortions the system is metallic, with ground state of predominantly $d^8\ligand$ character, where $\ligand$ denotes a ligand hole. For sufficiently large distortions ($δd_{\rm Ni-O} \sim 0.05 - 0.10Å$), however, a gap opens at the Fermi energy as the system enters a periodically distorted state alternating along the three crystallographic axes, with $(d^8\ligand^2)_{S=0}(d^8)_{S=1}$ character, where $S$ is the total spin. Thus the MI transition may be viewed as being driven by an internal volume "collapse" where the NiO$_6$ octahedra with two ligand holes shrink around their central Ni, while the remaining octahedra expand accordingly, resulting in the ($1/2,1/2,1/2$) superstructure observed in x-ray diffraction in the insulating phase. This insulating state is an example of a new type of charge ordering achieved without any actual movement of the charge.

preprint2013arXiv

Hidden Quasiparticles and Incoherent Photoemission Spectra in Na2IrO3

We study two Heisenberg-Kitaev t-J-like models on a honeycomb lattice, focusing on the zigzag magnetic phase of Na$_2$IrO$_3$, and investigate hole motion by exact diagonalization and variational methods. The spectral functions are quite distinct from those of cuprates and are dominated by large incoherent spectral weight at high energy, almost independent of the microscopic parameters --- a universal and generic feature for zigzag magnetic correlations. We explain why quasiparticles at low energy are strongly suppressed in the photoemission spectra and determine an analog of a pseudogap. We point out that the qualitative features of the predominantly incoherent spectra obtained within the two different models for the zigzag phase are similar, and they have remarkable similarity to recently reported angular resolved photoemission spectra for Na$_2$IrO$_3$.

preprint2013arXiv

Investigating polaron transitions with polar molecules

We determine the phase diagram of a polaron model with mixed breathing-mode and Su-Schrieffer-Heeger couplings and show that it has two sharp transitions, in contrast to pure models which exhibit one (for Su-Schrieffer-Heeger coupling) or no (for breathing-mode coupling) transition. Our results indicate that the physics of realistic mixed polaron models is much richer than that of simplified models. We then show that ultracold molecules trapped in optical lattices can be used to study precisely this mixed Hamiltonian, and that the relative contributions of the two couplings can be tuned with external electric fields. The parameters of current experimental set-ups place them in the region where one of the transitions occurs. We propose a scheme to measure the polaron dispersion using stimulated Raman spectroscopy.

preprint2013arXiv

The effect of dispersive optical phonons on the behaviour of a Holstein polaron

We use the approximation-free Bold Diagrammatic Monte Carlo technique to study the effects of a finite dispersion of the optical phonon mode on the properties of the Holstein polaron, especially its effective mass. For weak electron-phonon coupling the effect is very small, but it becomes significant for moderate and large electron-phonon coupling. The effective mass is found to increase (decrease) if the phonon dispersion has a negative (positive) curvature at the centre of the Brillouin zone.

preprint2012arXiv

A perturbational study of the lifetime of a Holstein polaron in the presence of weak disorder

Using the momentum average (MA) approximation, we find an analytical expression for the disorder-averaged Green's function of a Holstein polaron in a three-dimensional simple cubic lattice with random on-site energies. The on-site disorder is assumed to be weak compared to the kinetic energy of the polaron, and is treated perturbationally. Within this scheme, the states at the bottom of the polaron band are found to have an infinite lifetime, signaling a failure of perturbation theory at these energies. The higher-energy polaron states have a finite lifetime. We study this lifetime and the disorder-induced energy shift of these eigenstates for various strengths of disorder and electron-phonon coupling. We compare our findings to the predictions of Fermi's golden rule and the average T-matrix method, and find a significant quantitative discrepancy at strong electron-phonon coupling, where the polaron lifetime is much shorter than Fermi's golden rule prediction. We attribute this to the renormalization of the on-site potential by the electron-phonon coupling.

preprint2012arXiv

Magnon-mediated interactions between fermions depend strongly on the lattice structure

We propose two new methods to calculate exactly the spectrum of two spin-${1\over 2}$ charge carriers moving in a ferromagnetic background, at zero temperature. We find that if the spins are located on a different sublattice than that on which the fermions move, magnon-mediated effective interactions are very strong and can bind the fermions into low-energy bipolarons with triplet character. This never happens in models where spins and charge carriers share the same lattice, whether they are in the same band or in different bands. This proves that effective one-lattice models do not describe correctly the low-energy part of the two-carrier spectrum of a two-sublattice model, even though they may describe the low-energy single-carrier spectrum appropriately.

preprint2012arXiv

The role of the lattice structure in determining the magnon-mediated interactions between charge carriers doped into a magnetically ordered background

We use two recently proposed methods to calculate exactly the spectrum of two spin-${1\over 2}$ charge carriers moving in a ferromagnetic background, at zero temperature, for three types of models. By comparing the low-energy states in both the one-carrier and the two-carrier sectors, we analyze whether complex models with multiple sublattices can be accurately described by simpler Hamiltonians, such as one-band models. We find that while this is possible in the one-particle sector, the magnon-mediated interactions which are key to properly describe the two-carrier states of the complex model are not reproduced by the simpler models. We argue that this is true not just for ferromagnetic, but also for antiferromagnetic backgrounds. Our results question the ability of simple one-band models to accurately describe the low-energy physics of cuprate layers.

preprint2012arXiv

Trapping of three-dimensional Holstein polarons by various impurities

We study the bound states of a three-dimensional Holstein polaron near various kinds of single impurities, using the momentum average approximation. We show that the electron-phonon coupling is responsible for a strong renormalization of the impurity potential, resulting in an effective potential with significant retardation effects, which describes essential physics ignored by "instantaneous" approximations. The accuracy of our approximation is gauged by comparison with results from Diagrammatic Monte Carlo for the case of an impurity that modifies the on-site energy of the electron. We also discuss impurities that modify the local strength of the electron-phonon coupling, as well as isotope substitutions that change both the electron-phonon coupling and the phonon frequency, and contrast and highlight the difference between these cases.

preprint2011arXiv

A computational approach to a doped antiferromagnet: correlations between two spin-polarons in the lightly doped CuO$_2$ plane

We extend the methods recently introduced in Phys. Rev. Lett. 106 036401 (2011) to investigate correlations between two spin-polarons in a quasi-two-dimensional CuO2 layer. The low-energy wavefunctions for two doped holes introduced in a half-filled CuO2 plane with 32 copper and 64 oxygen sites are calculated explicitly using an efficient yet accurate truncation scheme to model the antiferromagnet background. The energetics and wavefucntions show that the charges form three-spin polarons and the spin is carried by a disturbance around the three-spin polaron core. The low-energy band results from the competition between the kinetic energy and a local attractive potential which favors d(x^2-y^2) states. Lastly, we point out features that are expected to be robust for larger systems.

preprint2011arXiv

Aharonov-Bohm interference for a hole in a two-dimensional Ising antiferromagnet in a transverse magnetic field

We show that a proper consideration of the contribution of Trugman loops leads to a fairly low effective mass for a hole moving in a square lattice Ising antiferromagnet, if the bare hopping and the exchange energy scales are comparable. This contradicts the general view that because of the absence of spin fluctuations, this effective mass must be extremely large. Moreover, in the presence of a transverse magnetic field, we show that the effective hopping integrals acquire an unusual dependence on the magnetic field, through Aharonov-Bohm interference, in addition to significant retardation effects. The effect of the Aharonov-Bohm interference on the cyclotron frequency (for small magnetic fields) and the Hofstadter butterfly (for large magnetic fields) is analyzed.

preprint2011arXiv

Few-particle Green's functions for strongly correlated systems on infinite lattices

We show how few-particle Green's functions can be calculated efficiently for models with nearest-neighbor hopping, for infinite lattices in any dimension. As an example, for one dimensional spinless fermions with both nearest-neighbor and second nearest-neighbor interactions, we investigate the ground states for up to 5 fermions. This allows us not only to find the stability region of various bound complexes, but also to infer the phase diagram at small but finite concentrations.

preprint2011arXiv

High-spin polaron in lightly doped CuO$_2$ planes

We device and investigate numerically a minimal yet detailed spin polaron model that describes lightly doped CuO$_2$ layers. The low-energy physics of a hole is studied by total-spin-resolved exact diagonalization on clusters of up to 32 CuO$_2$ unit cells, revealing features missed by previous studies. In particular, spin-polaron states with total spin 3/2 are the lowest eigenstates in several regions of the Brillouin zone. In these regions, and also at other points the quasiparticle weight is identically zero, indicating orthogonal states to those represented in the one electron Green's function. This highlights the importance of proper treatment of spin fluctuations in the many-body background.

preprint2011arXiv

Optical conductivity of the Holstein polaron

The Momentum Average approximation is used to derive a new kind of non-perturbational analytical expression for the optical conductivity (OC) of a Holstein polaron at zero temperature. This provides insight into the shape of the OC, by linking it to the structure of the polaron's phonon cloud. Our method works in any dimension, properly enforces selection rules, can be systematically improved, and also generalizes to momentum-dependent couplings. Its accuracy is demonstrated by a comparison with the first detailed set of three-dimensional numerical OC results, obtained using the approximation-free diagrammatic Monte Carlo method.

preprint2011arXiv

Reply to "Comment on 'High-Spin Polaron in Lightly Doped CuO$_2$ Planes'"

In arXiv:1108.5413v1, Lee and Lee use several comparisons to argue that the physics of the three-band model found (Phys. Rev. Lett. 106 036401 (2011)) can be explained in the one-band model's framework (Phys. Rev. Lett. 91 057001 (2003)). While superficial similarities exist between the two sets of results, for reasons discussed in this reply, we disagree that they describe the same physics.

preprint2010arXiv

Computational scheme for the spin-1/2 square lattice Heisenberg antiferromagnet based on an octapartite description of the square lattice

We introduce a novel parametrization scheme for the Hilbert space of a spin-1/2 Heisenberg antiferromagnet (AFM) based on an octapartite description of the square lattice. Our formulation provides an efficient yet systematic way to model the singlet ground-state wave function within a truncated basis. We demonstrate its effectiveness by using exact diagonalization to study systems with up to 64 spins. At significantly reduced computational cost, we obtain ground-state energies within less than 1% of the lowest values published. The non-iterative nature of this formulation and the fact that spatial symmetries remain unexploited lead to opportunities for extensions such as the incorporation into established iterative methods and the application to interacting models, for example in the study of the propagation of charged fermions in a 2D antiferromagnet.

preprint2010arXiv

Heat transport in quantum spin chains: the relevance of integrability

We investigate heat transport in various quantum spin chains, using the projector operator technique. We find that anomalous heat transport is linked not to the integrability of the Hamiltonian, but to whether it can be mapped to a model of non-interacting fermions. Our results also suggest how seemingly anomalous transport may occur at low temperatures in a much wider class of models.

preprint2010arXiv

Holstein magneto-polarons: from Landau levels to Hofstadter butterflies

We study the Holstein polaron in transverse magnetic field using non-perturbational methods. At strong fields and large coupling, we show that the polaron has a Hofstadter spectrum, however very distorted and of lower symmetry than that of a (heavier) bare particle. For weak magnetic fields, we identify non-perturbational behaviour of the Landau levels not previously known.

preprint2010arXiv

Kubo formula for finite size systems

We demonstrate that the proper calculation of the linear response for finite-size systems can only be performed if the coupling to the leads/baths is explicitly taken into consideration. We exemplify this by obtaining a Kubo-type formula for heat transport in a finite-size system coupled to two thermal baths, kept at different temperatures. We show that the proper calculation results in a well-behaved response, without the singular contributions from degenerate states encountered when Kubo formulae for infinite-size systems are inappropriately used for finite-size systems.

preprint2010arXiv

Momentum average approximation for models with boson-modulated hopping: Role of closed loops in the dynamical generation of a finite quasiparticle mass

We generalize the momentum average approximation to study the properties of single polarons in models with boson affected hopping, where the fermion-boson scattering depends explicitly on both the fermion's and the boson's momentum. As a specific example, we investigate the Edwards fermion-boson model in both one and two dimensions. In one dimension, this allows us to compare our results with exact diagonalization results, to validate the accuracy of our approximation. The generalization to two-dimensional lattices allows us to calculate the polaron's quasiparticle weight and dispersion throughout the Brillouin zone and to demonstrate the importance of Trugman loops in generating a finite effective mass even when the free fermion has an infinite mass.

preprint2010arXiv

The end points in the dispersion of Holstein polarons

We investigate the existence of end points in the dispersion of Holstein polarons in various dimensions, using the Momentum Average approximation which has proved to be very accurate for this model. An end point separates momenta for which the lowest-energy state is a discrete level, i.e., an infinitely-lived polaron, from those where the lowest-energy feature is a continuum in which the "polaron'" is signalled by a resonance with a finite lifetime. While such end points are known to not appear in 1D, we show here that they are generic in 3D if the particle-boson coupling is not too strong. The 2D case is "critical": a pure 2D Holstein model has no end points, like the 1D case. However, any amount of interlayer hopping leads to 3D-like behavior. As a result, such end points are expected to appear in the spectra of layered, quasi-2D systems described by Holstein models. Generalizations to other models are also briefly discussed.

preprint2009arXiv

Ferromagnetic spin-polaron on complex lattices

We present a simpler derivation of the exact solution of a spin-polaron in a ferromagnet and generalize it to complex lattices and/or longer range exchange interactions. As a specific example, we analyze a two-dimensional MnO$_2$-like lattice (as in the ferromagnetic layers in LaMnO$_3$) and discuss the properties of the resulting spin-polaron in various regimes. At strong couplings the solution is reminiscent of the Zhang-Rice singlet, however the electronic wavefunction involved in the singlet is dependent on the momentum of the singlet, and multiple bands may appear.

preprint2009arXiv

Holstein polaron in the presence of disorder

Non-local, inhomogeneous and retarded response observed in experiments is reproduced by introducing the Inhomogeneous Momentum Average (IMA) method to study single polaron problems with disorder in the on-site potential and/or spatial variations of the electron-phonon couplings and/or phonon frequencies. We show that the electron-phonon coupling gives rise to an additional inhomogeneous, strongly retarded potential, which makes instant approximations questionable. The accuracy of IMA is demonstrated by comparison with results from the approximation free Diagrammatic Monte Carlo (DMC) method. Its simplicity allows for easy study of many problems that were previously unaccessible. As an example, we show how inhomogeneities in the electron-phonon coupling lead to nonlocal, retarded response in scanning tunneling microscopy (STM) images.

preprint2009arXiv

Polaron formation in the presence of Rashba spin-orbit coupling: implications for spintronics

We study the effects of the Rashba spin-orbit coupling on the polaron formation, using a suitable generalization of the Momentum Average approximation. Contrary to previous investigations of this problem, we find that the spin-orbit interaction decreases the effective electron-phonon coupling. It is thus possible to lower the effective mass of the polaron by increasing the spin-orbit coupling. We also show that when the spin-orbit coupling is large as compared to the phonon energy, the polaron retains only one of the spin polarized bands in its coherent spectrum. This has major implications for the propagation of spin-polarized currents in such materials, and thus for spintronic applications.

preprint2009arXiv

The effect of impurities on spin polarized Zeeman bound states in dilute magnetic semiconductor-superconductor hybrids

We investigate the effect of single and multiple impurities on the Zeeman-localized, spin polarized bound states in dilute magnetic semiconductor hybrid system. Such bound states appear whenever a dilute magnetic semiconductor showing giant Zeeman effect is exposed to an external magnetic field showing nanoscale inhomogeneity. We consider the specific example of a superconductor-dilute magnetic semiconductor hybrid, calculate the energy spectrum and the wave functions of the bound states in the presence of a single impurity, and monitor the evolution of the bound state as a function of the impurity strength and impurity location with respect to the center of the Zeeman trapping potential. Our results have important experimental implications as they predict robust spin textures even for than than ideal samples. We find that for all realistic impurity strengths the Zeeman bound state survives the presence of the impurity. We also investigate the effect of a large number of impurities and perform ensemble averages with respect to the impurity locations. We find that the spin polarized Zeeman bound states are very robust, and they remain bound to the external field inhomogeneity throughout the experimentally relevant region of impurity concentration and scattering strength.

preprint2008arXiv

Electronic polarons and bipolarons in Fe-based superconductors: a pairing mechanism

Superconductivity is a fascinating example of how "more is different". It is due to electrons binding into bosonic Cooper pairs, which exhibit coherent behavior across a macroscopic sample. Finding the mechanism responsible for this binding is one of the more difficult tasks of condensed matter physics. For conventional superconductors the solution was given by the BCS theory as being due to exchange of phonons. For the cuprate high-Tc superconductors a widely-accepted explanation is still missing despite intense effort. The recently discovered Fe-based high-Tc superconductors pose now a new challenge. We present here a quantum mechanical theory for pnictides describing the influence of the large electronic polarizability of the heavy anions. We demonstrate that its inclusion results in electronic polarons as the low-energy quasi-particles and also unveils a pairing mechanism for these electronic polarons.

preprint2007arXiv

Single polaron properties of the breathing-mode Hamiltonian

We investigate numerically various properties of the one-dimensional (1D) breathing-mode polaron. We use an extension of a variational scheme to compute the energies and wave-functions of the two lowest-energy eigenstates for any momentum, as well as a scheme to compute directly the polaron Greens function. We contrast these results with results for the 1D Holstein polaron. In particular, we find that the crossover from a large to a small polaron is significantly sharper. Unlike for the Holstein model, at moderate and large couplings the breathing-mode polaron dispersion has non-monotonic dependence on the polaron momentum k. Neither of these aspects is revealed by a previous study based on the self-consistent Born approximation.

preprint2007arXiv

Spectral weight transfer in a disorder-broadened Landau level

In the absence of disorder, the degeneracy of a Landau level (LL) is $N=BA/ϕ_0$, where $B$ is the magnetic field, $A$ is the area of the sample and $ϕ_0=h/e$ is the magnetic flux quantum. With disorder, localized states appear at the top and bottom of the broadened LL, while states in the center of the LL (the critical region) remain delocalized. This well-known phenomenology is sufficient to explain most aspects of the Integer Quantum Hall Effect (IQHE) [1]. One unnoticed issue is where the new states appear as the magnetic field is increased. Here we demonstrate that they appear predominantly inside the critical region. This leads to a certain ``spectral ordering'' of the localized states that explains the stripes observed in measurements of the local inverse compressibility [2-3], of two-terminal conductance [4], and of Hall and longitudinal resistances [5] without invoking interactions as done in previous work [6-8].