Search

Valence bond (VB) states as the formation mechanism of Cooper pairs, eventually leading to high-temperature superconductivity, remain a controversial topic. Although various VB-like states find variational relevance in the description of specific spin models and quantum spin liquids, in the realm of many-body fermionic Hamiltonians, the evidence for such states as ground states wave functions remains elusive, challenging the valence-bond pairing mechanism. Here, we present evidence of a VB ground state with pairing tendencies, particularly at finite doping. We achieved this for the generic two-orbital Hubbard model in low dimension, where the VB states can be associated with the presence of the topological order manifested by edge states. Utilizing density-matrix renormalization group calculations, the study reveals key properties akin to those observed in superconductors' phase diagrams, such as pairing restricted to the regime of small but nonzero doping, presence of coherent pairs, and density oscillations in the charge sector.

We report neutron scattering, pressure-dependent AC calorimetry, and AC magnetic susceptibility measurements of triangular lattice NaYbSe$_2$. We observe a continuum of scattering, which is reproduced by matrix product simulations, and no phase transition is detected in any bulk measurements. Comparison to heat capacity simulations suggest the material is within the Heisenberg spin liquid phase. AC Susceptibility shows a significant 23~mK downturn, indicating a gap in the magnetic spectrum. The combination of a gap with no detectable magnetic order, comparison to theoretical models, and comparison to other $A$YbSe$_2$ compounds all strongly indicate NaYbSe$_2$ is within the quantum spin liquid phase. The gap also allows us to rule out a gapless Dirac spin liquid, with a gapped $\mathbb{Z}_2$ liquid the most natural explanation., Comment: 5 pages, 4 figures; 7 pages and 13 figures supplemental materials

One of the most famous quantum systems with topological properties, the spin $\mathcal{S}=1$ antiferromagnetic Heisenberg chain, is well-known to display exotic $\mathcal{S}=1/2$ edge states. However, this spin model has not been analyzed from the more general perspective of strongly correlated systems varying the electron-electron interaction strength. Here, we report the investigation of the emergence of the Haldane edge in a system of interacting electrons -- the two-orbital Hubbard model -- with increasing repulsion strength $U$ and Hund interaction $J_\mathrm{H}$. We show that interactions not only form the magnetic moments but also form a topologically nontrivial fermionic many-body ground-state with zero-energy edge states. Specifically, upon increasing the strength of the Hubbard repulsion and Hund exchange, we identify a sharp transition point separating topologically trivial and nontrivial ground-states. Surprisingly, such a behaviour appears already at rather small values of the interaction, in a regime where the magnetic moments are barely developed.

Spectroscopy experiments are routinely used to characterize the behavior of strongly correlated systems. An in-depth understanding of the different spectral features is thus essential. Here, we show that the spectrum of the multiorbital Hubbard model exhibits unique Hund \ms{bands} that occur at energies given only by the Hund coupling $J_\mathrm{H}$, as distinct from the Hubbard satellites following the interaction $U$. We focus on experimentally relevant single-particle and optical spectra that we calculate for a model related to iron chalcogenide ladders. The calculations are performed via the density-matrix renormalization group and Lanczos methods. The generality of the implications is verified by considering a generic multiorbital model within dynamical mean-field theory.

While geometrically frustrated quantum magnets are known for a variety of exotic spin states that are of great interests of understanding emergent phenomena as well as enabling revolutionary quantum technologies, most of them are necessarily good insulators which are difficult to be integrated with modern electrical circuit that relies on moving charge carriers. The grand challenge of converting fluctuations and excitations of frustrated moments into electronic responses is finding ways to introduce charge carriers that interact with the localized spins without destroying the spin states. Here, we show that, by designing a Bi2Ir2O7/Dy2Ti2O7 heterostructure, the breaking of the spin ice rule in insulating Dy2Ti2O7 can lead to a charge response in the Bi2Ir2O7 conducting layer that can be detected as anomalous magnetoresistance. These results demonstrate a novel and feasible interfacial approach for electronically probing exotic spin states in insulating magnets, laying out a blueprint for the metallization of frustrated quantum magnets.

Topological phases of matter are among the most intriguing research directions in Condensed Matter Physics. It is known that superconductivity induced on a topological insulator's surface can lead to exotic Majorana modes, the main ingredient of many proposed quantum computation schemes. In this context, the iron-based high critical temperature superconductors are a promising platform to host such an exotic phenomenon in real condensed-matter compounds. The Coulomb interaction is commonly believed to be vital for the magnetic and superconducting properties of these systems. This work bridges these two perspectives and shows that the Coulomb interaction can also drive a canonical superconductor with orbital degrees of freedom into the topological state. Namely, we show that above a critical value of the Hubbard interaction the system simultaneously develops spiral spin order, a highly unusual triplet amplitude in superconductivity, and, remarkably, Majorana fermions at the edges of the system.

We present a comprehensive study of the spin excitations - as measured by the dynamical spin structure factor $S(q,\omega)$ - of the so-called block-magnetic state of low-dimensional orbital-selective Mott insulators. We realize this state via both a multi-orbital Hubbard model and a generalized Kondo-Heisenberg Hamiltonian. Due to various competing energy scales present in the models, the system develops periodic ferromagnetic islands of various shapes and sizes, which are antiferromagnetically coupled. The 2$\times$2 particular case was already found experimentally in the ladder material BaFe$_2$Se$_3$ that becomes superconducting under pressure. Here we discuss the electronic density as well as Hubbard and Hund coupling dependence of $S(q,\omega)$ using density matrix renormalization group method. Several interesting features were identified: (1) An acoustic (dispersive spin-wave) mode develops. (2) The spin-wave bandwidth establishes a new energy scale that is strongly dependent on the size of the magnetic island and becomes abnormally small for large clusters. (3) Optical (dispersionless spin excitation) modes are present for all block states studied here. In addition, a variety of phenomenological spin Hamiltonians have been investigated but none matches entirely our results that were obtained primarily at intermediate Hubbard $U$ strengths. Our comprehensive analysis provides theoretical guidance and motivation to crystal growers to search for appropriate candidate materials to realize the block states, and to neutron scattering experimentalists to confirm the exotic dynamical magnetic properties unveiled here, with a rich mixture of acoustic and optical features.

Competing interactions in Quantum Materials induce novel states of matter such as frustrated magnets, an extensive field of research both from the theoretical and experimental perspectives. Here, we show that competing energy scales present in the low-dimensional orbital-selective Mott phase (OSMP) induce an exotic magnetic order, never reported before. Earlier neutron scattering experiments on iron-based 123 ladder materials, where OSMP is relevant, already confirmed our previous theoretical prediction of block-magnetism (magnetic order of the form $\uparrow\uparrow\downarrow\downarrow$). Now we argue that another novel phase can be stabilized in multi-orbital Hubbard models, the {\it block-spiral state}. In this state, the magnetic islands form a spiral propagating through the chain but with the blocks maintaining their identity, namely rigidly rotating. This new spiral state is stabilized without any apparent frustration, the common avenue to generate spiral arrangements in multiferroics. By examining the behaviour of the electronic degrees of freedom, parity breaking quasiparticles are revealed. Finally, a simple phenomenological model that accurately captures the macroscopic spin spiral arrangement is also introduced, and fingerprints for the neutron scattering experimental detection of our new state are provided.

Inelastic neutron scattering recently confirmed the theoretical prediction of a $\uparrow\uparrow\downarrow\downarrow$-magnetic state along the legs of quasi-one-dimensional (quasi-1D) iron-based ladders in the orbital-selective Mott phase (OSMP). We show here that electron-doping of the OSMP induces a whole class of novel block-states with a variety of periodicities beyond the previously reported $\pi/2$ pattern. We discuss the magnetic phase diagram of the OSMP regime that could be tested by neutrons once appropriate quasi-1D quantum materials with the appropriate dopings are identified.

Iron-based superconductors display a variety of magnetic phases originating in the competition between electronic, orbital, and spin degrees of freedom. Previous theoretical investigations of the multi-orbital Hubbard model in one dimension revealed the existence of an orbital-selective Mott phase (OSMP) with block spin order. Recent inelastic neutron scattering (INS) experiments on the BaFe$_2$Se$_3$ ladder compound confirmed the relevance of the block-OSMP. Moreover, the powder INS spectrum reveled an unexpected structure, containing both low-energy acoustic and high-energy optical modes. Here we present the theoretical prediction for the dynamical spin structure factor within a block-OSMP regime using the density-matrix renormalization group method. In agreement with experiments we find two dominant features: low-energy dispersive and high-energy dispersionless modes. We argue that the former represents the spin-wave-like dynamics of the block ferromagnetic islands, while the latter is attributed to a novel type of local on-site spin excitations controlled by the Hund coupling.

We present a method for computing resonant inelastic x-ray scattering (RIXS) spectra in one-dimensional systems using the density matrix renormalization group (DMRG) method. By using DMRG to address the problem, we shift the computational bottleneck from the memory requirements associated with exact diagonalization (ED) calculations to the computational time associated with the DMRG algorithm. This approach is then used to obtain RIXS spectra on cluster sizes well beyond state-of-the-art ED techniques. Using this new procedure, we compute the low-energy magnetic excitations observed in Cu $L$-edge RIXS for the challenging corner shared CuO$_4$ chains, both for large multi-orbital clusters and downfolded $t$-$J$ chains. We are able to directly compare results obtained from both models defined in clusters with identical momentum resolution. In the strong coupling limit, we find that the downfolded $t$-$J$ model captures the features of the magnetic excitations probed by RIXS after a uniform scaling of the spectra is taken into account., Comment: Main text: 6 pages, 4 figures; Supplementary Material: 4 pages, 3 figures

We study the magnetic and charge dynamical response of a Hubbard model in a two-leg ladder geometry using the density matrix renormalization group (DMRG) method and the random phase approximation within the fluctuation-exchange approximation (RPA+FLEX). Our calculations reveal that RPA+FLEX can capture the main features of the magnetic response from weak up to intermediate Hubbard repulsion for doped ladders, when compared with the numerically exact DMRG results. However, while at weak Hubbard repulsion both the spin and charge spectra can be understood in terms of weakly-interacting electron-hole excitations across the Fermi surface, at intermediate coupling DMRG shows gapped spin excitations at large momentum transfer that remain gapless within the RPA+FLEX approximation. For the charge response, RPA+FLEX can only reproduce the main features of the DMRG spectra at weak coupling and high doping levels, while it shows an incoherent character away from this limit. Overall, our analysis shows that RPA+FLEX works surprisingly well for spin excitations at weak and intermediate Hubbard $U$ values even in the difficult low-dimensional geometry such as a two-leg ladder. Finally, we discuss the implications of our results for neutron scattering and resonant inelastic x-ray scattering experiments on two-leg ladder cuprate compounds., Comment: 9 pages, 8 figures, submitted to Physical Review B

Magnetic interactions are widely believed to play a crucial role in the microscopic mechanism leading to high critical temperature superconductivity. It is therefore important to study the signatures of pairing in the magnetic excitation spectrum of simple models known to show unconventional superconducting tendencies. Using the Density Matrix Renormalization Group technique, we calculate the dynamical spin structure factor $S({\bf k},\omega)$ of a generalized $t-U-J$ Hubbard model away from half-filling in a two-leg ladder geometry. The addition of $J$ enhances pairing tendencies. We analyze quantitatively the signatures of pairing in the magnetic excitation spectra. We found that the superconducting pair-correlation strength, that can be estimated independently from ground state properties, is closely correlated with the integrated low-energy magnetic spectral weight in the vicinity of $(\pi,\pi)$. In this wave-vector region, robust spin incommensurate features develop with increasing doping. The branch of the spectrum with rung direction wave-vector $k_{rung}=0$ does not change substantially with doping where pairing dominates, and thus plays a minor role. We discuss the implications of our results for neutron scattering experiments, where the spin excitation dynamics of hole-doped quasi-one dimensional magnetic materials can be measured, and also address implications for recent resonant inelastic X-ray scattering experiments., Comment: 9 pages, 9 figures, submitted to Physical Review B

We present a time-dependent density-matrix renormalization group investigation of the quantum distillation process within the Fermi--Hubbard model on a quasi-1D ladder geometry. The term distillation refers to the dynamical, spatial separation of singlons and doublons in the sudden expansion of interacting particles in an optical lattice, i.e., the release of a cloud of atoms from a trapping potential. Remarkably, quantum distillation can lead to a contraction of the doublon cloud, resulting in an increased density of the doublons in the core region compared to the initial state. As a main result, we show that this phenomenon is not limited to chains that were previously studied. Interestingly, there are additional dynamical processes on the two-leg ladder such as density oscillations and selftrapping of defects that lead to a less efficient distillation process. An investigation of the time evolution starting from product states provides an explanation for this behaviour. Initial product states are also considered, since in optical lattice experiments such states are often used as the initial setup. We propose configurations that lead to a fast and efficient quantum distillation.

The recent discovery of superconductivity under high pressure in the ladder compound BaFe$_2$S$_3$ has opened a new field of research in iron-based superconductors with focus on quasi one-dimensional geometries. In this publication, using the Density Matrix Renormalization Group technique, we study a two-orbital Hubbard model defined in one dimensional chains. Our main result is the presence of hole binding tendencies at intermediate Hubbard $U$ repulsion and robust Hund coupling $J_H/U=0.25$. Binding does not occur neither in weak coupling nor at very strong coupling. The pair-pair correlations that are dominant near half-filling, or of similar strength as the charge and spin correlation channels, involve hole-pair operators that are spin singlets, use nearest-neighbor sites, and employ different orbitals for each hole. The Hund coupling strength, presence of robust magnetic moments, and antiferromagnetic correlations among them are important for the binding tendencies found here.

Iron telluride FeTe is known to display bicollinear magnetic order at low temperatures together with a monoclinic lattice distortion. Because the bicollinear order can involve two different wavevectors $(\pi/2,\pi/2)$ and $(\pi/2,-\pi/2)$, symmetry considerations allow for the possible stabilization of a nematic state with short-range bicollinear order coupled to monoclinic lattice distortions at a $T_S$ higher than the temperature $T_N$ where long-range bicollinear order fully develops. As a concrete example, the three-orbitals spin-fermion model for iron telluride is studied with an additional coupling $\tilde\lambda_{12}$ between the monoclinic lattice strain and an orbital-nematic order parameter with $B_{2g}$ symmetry. Monte Carlo simulations show that with increasing $\tilde\lambda_{12}$ the first-order transition characteristic of FeTe splits and bicollinear nematicity is stabilized in a (narrow) temperature range. In this new regime the lattice is monoclinically distorted and short-range spin and orbital order breaks rotational invariance. A discussion of possible realizations of this exotic state is provided.

We study the effects of charge degrees of freedom on the spin excitation dynamics in quasi-one dimensional magnetic materials. Using the density matrix renormalization group method, we calculate the dynamical spin structure factor of the Hubbard model at half electronic filling on a chain and on a ladder geometry, and compare the results with those obtained using the Heisenberg model, where charge degrees of freedom are considered frozen. For both chains and two-leg ladders, we find that the Hubbard model spectrum qualitatively resembles the Heisenberg spectrum -- with low-energy peaks resembling spinonic excitations -- already at intermediate on-site repulsion as small as $U/t \sim 2-3$, although ratios of peak intensities at different momenta continue evolving with increasing $U/t$ converging only slowly to the Heisenberg limit. We discuss the implications of these results for neutron scattering experiments and we propose criteria to establish the values of $U/t$ of quasi-one dimensional systems described by one-orbital Hubbard models from experimental information., Comment: 10 pages, 9 figures

We study non-local correlations in a three-orbital Hubbard model defined on an extended one-dimensional chain using determinant quantum Monte Carlo and density matrix renormalization group methods. We focus on a parameter with robust Hund's coupling, which produces an orbital selective Mott phase (OSMP) at intermediate values of the Hubbard U, as well as an orbitally ordered ferromagnetic insulating state at stronger coupling. An examination of the orbital and spin-correlation functions indicates that the orbital ordering occurs before the onset of magnetic correlations in this parameter regime as a function of temperature. In the OSMP, we find that the self-energy for the itinerant electrons is momentum dependent, indicating a degree of non-local correlations while the localized electrons have largely momentum independent self-energies. These non-local correlations also produce relative shifts of the hole-like and electron-like bands within our model. The overall momentum dependence of these quantities is strongly suppressed in the orbitally-ordered insulating phase., Comment: we have 12 pages and 9 figures

Polycrystalline samples of the quarter-doped manganites $R_{0.75}$Ca$_{0.25}$MnO$_3$ ($R$ = Y, Tb, Dy, Ho, and Er) were studied by X-ray diffraction and AC/DC susceptibility measurements. All five samples are orthorhombic and exhibit similar magnetic properties: enhanced ferromagnetism below $T_1$ ($\sim80$ K) and a spin glass (SG) state below $T_{SG}$ ($\sim30$ K). With increasing $R^{3+}$ ionic size, both $T_1$ and $T_{SG}$ generally increase. The single crystal neutron diffraction results on Tb$_{0.75}$Ca$_{0.25}$MnO$_3$ revealed that the SG state is mainly composed of a short-range ordered version of a novel canted (i.e. noncollinear) antiferromagnetic spin state. Furthermore, calculations based on the double exchange model for quarter-doped manganites reveal that this new magnetic phase provides a transition state between the ferromagnetic state and the theoretically predicted spin-orthogonal stripe phase., Comment: 10 pages, 10 figures. Submitted to Physical Review B

Using the time-dependent density-matrix renormalization group (tDMRG), we study the time evolution of electron wave packets in one-dimensional (1D) metal-superconductor heterostructures. The results show Andreev reflection at the interface, as expected. By combining these results with the well-known single-spin-species electron-hole transformation in the Hubbard model, we predict an analogous spin Andreev reflection in metal-Mott insulator heterostructures. This effect is numerically confirmed using 1D tDMRG, but it is expected to be present also in higher dimensions, as well as in more general Hamiltonians. We present an intuitive picture of the spin reflection, analogous to that of Andreev reflection at metal-superconductors interfaces. This allows us to discuss a novel antiferromagnetic proximity effect. Possible experimental realizations are discussed., Comment: K.A.A. and J.R. contributed equally to this work. 5+2 pages, 4+1 figures

Using time-dependent density-matrix renormalization group, we study the time evolution of electronic wave packets in the one-dimensional extended Hubbard model with on-site and nearest neighbor repulsion, U and V, respectively. As expected, the wave packets separate into spin-only and charge-only excitations (spin-charge separation). Charge and spin velocities exhibit non-monotonic dependence on V. For small and intermediate values of V, both velocities increase with V. However, the charge velocity exhibits a stronger dependence than that of the spin, leading to a more pronounced spin-charge separation. Charge fractionalization, on the other hand, is weakly affected by V. The results are explained in terms of Luttinger liquid theory in the weak-coupling limit, and an effective model in the strong-coupling regime., Comment: 6 pages, 4 figures

Recent studies have reported the existence of an epitaxially-stabilized tetragonal-like ('T-like') monoclinic phase in BiFeO3 thin-films with high levels of compressive strain. While their structural and ferroelectric properties are different than those of rhombohedral-like ('R-like') films with lower levels of strain, little information exists on magnetic properties. Here, we report a detailed neutron scattering study of a nearly phase-pure film of T-like BiFeO3. By tracking the temperature dependence and relative intensity of several superstructure peaks in the reciprocal lattice cell, we confirm antiferromagnetism with largely G-type character and TN = 324 K, significantly below a structural phase transition at 375 K, contrary to previous reports. Evidence for a second transition, possibly a minority magnetic phase with C-type character is also reported with TN = 260 K. The co-existence of the two magnetic phases in T-like BiFeO3 and the difference in ordering temperatures between R-like and T-like systems is explained through simple Fe-O-Fe bond distance considerations., Comment: 9 pages, 7 figures. Slight modifications to introductory text and presentation of data. Supplementary Material section added

A detailed description of the time-step-targetting time evolution method within the DMRG algorithm is presented. The focus of this publication is on the implementation of the algorithm, and on its generic application. The case of one-site excitations within a Hubbard model is analyzed as a test for the algorithm, using open chains and two-leg ladder geometries. The accuracy of the procedure in the case of the recently discussed holon-doublon photo excitations of Mott insulators is also analyzed. Performance and parallelization issues are discussed. In addition, the full open-source code is provided as supplementary material., Comment: RevTeX4, 9 pages, 5 figures

Electric-field controlled exchange bias in a heterostructure composed of the ferromagnetic manganite La0.7Sr0.3MO3 and the ferroelectric antiferromagnetic BiFeO3 has recently been demonstrated experimentally. By means of a microscopic model Hamiltonian we provide a possible explanation of the origin of this magnetoelectric coupling. We find, in agreement with experimental results, a net ferromagnetic moment at the BiFeO3 interface. The induced ferromagnetic moment is the result of the competition between the e_g-electrons double exchange and the t_2g-spins antiferromagnetic superexchange that dominate in bulk BiFeO3. The balance of these simultaneous ferromagnetic and antiferromagnetic tendencies is strongly affected by the interfacial electronic charge density which, in turn, can be controlled by the BiFeO3 ferroelectric polarization., Comment: 4 pages, including 4 figures

We calculate the non-equilibrium electronic transport properties of a one-dimensional interacting chain at half filling, coupled to non-interacting leads. The interacting chain is initially in a Mott insulator state that is driven out of equilibrium by applying a strong bias voltage between the leads. For bias voltages above a certain threshold we observe the breakdown of the Mott insulator state and the establishment of a steady-state electronic current through the system. Based on extensive time-dependent density matrix renormalization group simulations, we show that this steady-state current always has the same functional dependence on voltage, independent of the microscopic details of the model and relate the value of the threshold to the Lieb-Wu gap. We frame our results in terms of the Landau-Zener dielectric breakdown picture. Finally, we also discuss the real-time evolution of the current, and characterize the current-carrying state resulting from the breakdown of the Mott insulator by computing the double occupancy, the spin structure factor, and the entanglement entropy., Comment: 12 pages RevTex4, 12 eps figures, as published, minor revisions

Using the Density Matrix Renormalization Group and two-leg ladders, we investigate an electronic two-orbital Hubbard model including plaquette diagonal hopping amplitudes. Our goal is to search for regimes where charges added to the undoped state form pairs, presumably a precursor of a superconducting state.For the electronic density $\rho=2$, i.e. the undoped limit, our investigations show a robust $(\pi,0)$ antiferromagnetic ground state, as in previous investigations. Doping away from $\rho=2$ and for large values of the Hund coupling $J$, a ferromagnetic region is found to be stable. Moreover, when the interorbital on-site Hubbard repulsion is smaller than the Hund coupling, i.e. for $U'

We use neutron scattering to investigate the doping evolution of the magnetic correlations in the single-layer manganite $\rm Pr_{\it 1-x}Ca_{\it 1+x}MnO_4$, away from the $x=0.5$ composition where the CE-type commensurate antiferromagnetic (AF) structure is stable. We find that short-range incommensurate spin correlations develop as the system is electron doped ($x<0.5$), which coexist with the CE-type AF order. This suggests that electron doping in this system induces an inhomogeneous electronic self-organization, where commensurate AF patches with $x=0.5$ are separated by electron-rich domain walls with short range magnetic correlations. This behavior is strikingly different than for the perovksite $\rm Pr_{\it 1-x}Ca_{\it x}MnO_3$, where the long-range CE-type commensurate AF structure is stable over a wide range of electron or hole doping around $x=0.5$., Comment: 4 figures. To be published in Phys. Rev. Lett

The spatial length of the Kondo screening is still a controversial issue related to Kondo physics. While renormalization group and Bethe Anzats solutions have provided detailed information about the thermodynamics of magnetic impurities, they are insufficient to study the effect on the surrounding electrons, i.e., the spatial range of the correlations created by the Kondo effect between the localized magnetic moment and the conduction electrons. The objective of this work is to present a quantitative way of measuring the extension of these correlations by studying their effect directly on the local density of states (LDOS) at arbitrary distances from the impurity. The numerical techniques used, the Embedded Cluster Approximation, the Finite U Slave Bosons, and Numerical Renormalization Group, calculate the Green functions in real space. With this information, one can calculate how the local density of states away from the impurity is modified by its presence, below and above the Kondo temperature, and then estimate the range of the disturbances in the non-interacting Fermi sea due to the Kondo effect, and how it changes with the Kondo temperature $T_{\rm K}$. The results obtained agree with results obtained through spin-spin correlations, showing that the LDOS captures the phenomenology of the Kondo cloud as well. To the best of our knowledge, it is the first time that the LDOS is used to estimate the extension of the Kondo cloud.

One of the main open problems in the field of transport in strongly interacting nanostructures is the understanding of currents beyond the linear response regime. In this work, we consider the single-impurity Anderson model and use the adaptive time-dependent density matrix renormalization group (tDMRG) method to compute real-time currents out of equilibrium. We first focus on the particle-hole symmetric point where Kondo correlations are the strongest and then extend the study of the nonequilibrium transport to the mixed-valence regime. As a main result, we present accurate data for the current-voltage characteristics of this model., Comment: 4 pages Revtex4, 4 eps figures, additional references, revised version, to appear in Phys. Rev. B

Correlations between particles can lead to subtle and sometimes counterintuitive phenomena. We analyze one such case, occurring during the sudden expansion of fermions in a lattice when the initial state has a strong admixture of double occupancies. We promote the notion of quantum distillation: during the expansion, and in the presence of strongly repulsive interactions, doublons group together, forming a nearly ideal band insulator, which is metastable with a low entropy. We propose that this effect could be used for cooling purposes in experiments with two-component Fermi gases., Comment: Final version as published, minor revisions, more discussion on the cooling proposal, includes auxiliary material (EPAPS), 7 pages Revtex 4, 12 eps figures

A two orbital model for the new Fe-based superconductors is studied using the Lanczos method as well as pairing mean-field approximations. Our main goals are (i) to provide a comprehensive analysis of this model using numerical techniques with focus on half-filling and on the state with two more electrons than half-filling and (ii) to investigate the nodal structure of the mean-field superconducting state and compare the results with angle-resolved photoemission data. In particular, we provide evidence that at half-filling spin 'stripes', as observed experimentally, dominate over competing states. Depending on parameters, the state with two more electrons added to half filling is either triplet or singlet. Since experiments suggest spin singlet pairs, our focus is on this state. Under rotation, it transforms as the B_2g representation of the D_4h group. We also show that the s+/- pairing operator transforms as A_1g and becomes dominant only in an unphysical regime of the model where the undoped state is an insulator. For robust values of the effective electronic attraction producing the Cooper pairs, assumption compatible with recent angle-resolved photoemission (ARPES) results that suggesting small Cooper-pair size, the nodes of the two-orbital model are found to be located only at the electron pockets. Since recent ARPES efforts have searched for nodes at the hole pockets or only in a few directions at the electron pockets, our results for the nodal distribution may help to guide future experiments. More in general, the investigations reported here aim to establish several of the properties of the two orbital model. Only a detailed comparison with experiments will clarify how far this simple model present a valid description of the Fe pnictides.

We apply the adaptive time-dependent Density Matrix Renormalization Group method (tDMRG) to the study of transport properties of quantum-dot systems connected to metallic leads. Finite-size effects make the usual tDMRG description of the Kondo regime a numerically demanding task. We show that such effects can be attenuated by describing the leads by "Wilson chains", in which the hopping matrix elements decay exponentially away from the impurity ($t_n \propto \Lambda^{-n/2}$). For a given system size and in the linear response regime, results for $\Lambda > 1$ show several improvements over the undamped, $\Lambda=1$ case: perfect conductance is obtained deeper in the strongly interacting regime and current plateaus remain well defined for longer time scales. Similar improvements were obtained in the finite-bias regime up to bias voltages of the order of the Kondo temperature. These results show that, with the proposed modification, the tDMRG characterization of Kondo correlations in the transport properties can be substantially improved, while it turns out to be sufficient to work with much smaller system sizes. We discuss the numerical cost of this approach with respect to the necessary system sizes and the entanglement growth during the time-evolution., Comment: 10 pages, 9 figures. Updated: new section on finite bias added

This work proposes a new approach to study transport properties of highly correlated local structures. The method, dubbed the Logarithmic Discretization Embedded Cluster Approximation (LDECA), consists of diagonalizing a finite cluster containing the many-body terms of the Hamiltonian and embedding it into the rest of the system, combined with Wilson's idea of a logarithmic discretization of the representation of the Hamiltonian. The physics associated with both one embedded dot and a double-dot side-coupled to leads is discussed in detail. In the former case, the results perfectly agree with Bethe ansatz data, while in the latter, the physics obtained is framed in the conceptual background of a two-stage Kondo problem. A many-body formalism provides a solid theoretical foundation to the method. We argue that LDECA is well suited to study complicated problems such as transport through molecules or quantum dot structures with complex ground states., Comment: 17 pages, 13 figures

A two-orbital model for Fe-pnictide superconductors is investigated using computational techniques on two-dimensional square clusters. The hopping amplitudes are derived from orbital overlap integrals, or by band structure fits, and the spin frustrating effect of the plaquette-diagonal Fe-Fe hopping is remarked. A spin 'striped' state is stable in a broad range of couplings in the undoped regime, in agreement with neutron scattering. Adding two electrons to the undoped ground state of a small cluster, the dominant pairing operators are found. Depending on parameters, two pairing operators were identified: they involve inter-xz-yz orbital combinations forming spin singlets or triplets, transforming according to the B_2g and A_2g representations of the D_4h group, respectively., Comment: 4 pages, 4 figures, terminology for pairing symmetries improved, accepted for Phys. Rev. Lett

We study the real-time dynamics of a pair hole/doubly-occupied-site, namely a holon and a doublon, in a 1D Hubbard insulator with on-site and nearest-neighbor Coulomb repulsion. Our analysis shows that the pair is long-lived and the expected decay mechanism to underlying spin excitations is actually inefficient. For a nonzero inter-site Coulomb repulsion, we observe that part of the wave-function remains in a bound state. Our study also provides insight on the holon-doublon propagation in real space. Due to the one-dimensional nature of the problem, these particles move in opposite directions even in the absence of an applied electric field. The potential relevance of our results to solar cell applications is discussed., Comment: To Appear in Phys. Rev. Lett

The Kondo lattice model enlarged by an antiferromagnetic coupling $J_{\rm AF}$ between the localized spins is here investigated using computational techniques. Our results suggest the existence of a d-wave superconducting phase close to half-filling mediated by antiferromagnetic fluctuations. This establishes a closer connection between theory and heavy fermion experiments than currently provided by the standard Kondo lattice model with $J_{\rm AF}$=0., Comment: 4 pages, 4 figures, to appear in PRL

The one-particle spectral function of a state formed by superconducting (SC) clusters is studied via Monte Carlo techniques. The clusters have similar SC amplitudes but randomly distributed phases. This state is stabilized by the competition with anti-ferromagnetism, after quenched disorder is introduced. Fermi arcs between the critical temperature Tc and the cluster formation temperature scale T* are observed, similarly as in the pseudo-gap state of the cuprates. The arcs originate at metallic regions in between the neighboring clusters that present large SC phase differences., Comment: 4 figures, RevTeX4

A spin-tunnel-junction based on manganites, with La$_{1-x}$Sr$_x$MnO$_3$ (LSMO) as ferromagnetic metallic electrodes and the undoped parent compound LaMnO$_3$ (LMO) as insulating barrier, is here theoretically discussed using double exchange model Hamiltonians and numerical techniques. For an even number of LMO layers, the ground state is shown to have anti-parallel LSMO magnetic moments. This highly resistive, but fragile, state is easily destabilized by small magnetic fields, which orient the LSMO moments in the direction of the field. The magnetoresistance associated with this transition is very large, according to Monte Carlo and Density Matrix Renormalization Group studies. The influence of temperature, the case of an odd number of LMO layers, and the differences between LMO and SrTiO$_3$ as barriers are also addressed. General trends are discussed., Comment: 14 pages, 18 figures

We study the sudden expansion of strongly correlated fermions in a one-dimensional lattice, utilizing the time-dependent density-matrix renormalization group method. Our focus is on the behavior of experimental observables such as the density, the momentum distribution function, and the density and spin structure factors. As our main result, we show that correlations in the transient regime can be accurately described by equilibrium reference systems. In addition, we find that the expansion from a Mott insulator produces distinctive peaks in the momentum distribution function at |k| ~ pi/2, accompanied by the onset of power-law correlations., Comment: 8 pages Revtex4, 7 eps-figures, minor changes, additional references, as published

In this manuscript, recent theoretical investigations by the authors in the area of oxide multilayers are briefly reviewed. The calculations were carried out using model Hamiltonians and a variety of non-perturbative techniques. Moreover, new results are also included here. They correspond to the generation of a metallic state by mixing insulators in a multilayer geometry, using the Hubbard and Double Exchange models. For the latter, the resulting metallic state is also ferromagnetic. This illustrates how electron or hole doping via transfer of charge in multilayers can lead to the study of phase diagrams of transition metal oxides in the clean limit. Currently, these phase diagrams are much affected by the disordering standard chemical doping procedure, which introduces quenched disorder in the material., Comment: 14 pages, 9 figures. Invited article for a special issue of JPCM on Metal Oxide Thin Films; minor changes in the text

The numerical analysis of strongly interacting nanostructures requires powerful techniques. Recently developed methods, such as the time-dependent density matrix renormalization group (tDMRG) approach or the embedded-cluster approximation (ECA), rely on the numerical solution of clusters of finite size. For the interpretation of numerical results, it is therefore crucial to understand finite-size effects in detail. In this work, we present a careful finite-size analysis for the examples of one quantum dot, as well as three serially connected quantum dots. Depending on odd-even effects, physically quite different results may emerge from clusters that do not differ much in their size. We provide a solution to a recent controversy over results obtained with ECA for three quantum dots. In particular, using the optimum clusters discussed in this paper, the parameter range in which ECA can reliably be applied is increased, as we show for the case of three quantum dots. As a practical procedure, we propose that a comparison of results for static quantities against those of quasi-exact methods, such as the ground-state density matrix renormalization group (DMRG) method or exact diagonalization, serves to identify the optimum cluster type. In the examples studied here, we find that to observe signatures of the Kondo effect in finite systems, the best clusters involving dots and leads must have a total z-component of the spin equal to zero., Comment: 16 pages, 14 figures, revised version to appear in Eur. Phys. J. B, additional references

The electron doping of undoped high-$T_c$ cuprates via the transfer of charge from manganites (or other oxides) using heterostructure geometries is here theoretically discussed. This possibility is mainly addressed via a detailed analysis of photoemission and diffusion voltage experiments, which locate the Fermi level of manganites above the bottom of the upper Hubbard band of some cuprate parent compounds. A diagram with the relative location of Fermi levels and gaps for several oxides is presented. The procedure discussed here is generic, allowing for the qualitative prediction of the charge flow direction at several oxide interfaces. The addition of electrons to antiferromagnetic Cu oxides may lead to a superconducting state at the interface with minimal quenched disorder. Model calculations using static and dynamical mean-field theory, supplemented by a Poisson equation formalism to address charge redistribution at the interface, support this view. The magnetic state of the manganites could be antiferromagnetic or ferromagnetic. The former is better to induce superconductivity than the latter, since the spin-polarized charge transfer will be detrimental to singlet superconductivity. It is concluded that in spite of the robust Hubbard gaps, the electron doping of undoped cuprates at interfaces appears possible, and its realization may open an exciting area of research in oxide heterostructures., Comment: 12 pages, 9 figures

The crossover between an impurity band (IB) and a valence band (VB) regime as a function of the magnetic impurity concentration in models for diluted magnetic semiconductors (DMS) is studied systematically by taking into consideration the Coulomb attraction between the carriers and the magnetic impurities. The density of states and the ferromagnetic transition temperature of a Spin-Fermion model applied to DMS are evaluated using Dynamical Mean-Field Theory (DMFT) and Monte Carlo (MC) calculations. It is shown that the addition of a square-well-like attractive potential can generate an IB at small enough Mn doping $x$ for values of the $p-d$ exchange $J$ that are not strong enough to generate one by themselves. We observe that the IB merges with the VB when $x >= x_c$ where $x_c$ is a function of $J$ and the Coulomb attraction strength $V$. Using MC calculations, we demonstrate that the range of the Coulomb attraction plays an important role. While the on-site attraction, that has been used in previous numerical simulations, effectively renormalizes $J$ for all values of $x$, an unphysical result, a nearest-neighbor range attraction renormalizes $J$ only at very low dopings, i.e., until the bound holes wave functions start to overlap. Thus, our results indicate that the Coulomb attraction can be neglected to study Mn doped GaSb, GaAs, and GaP in the relevant doping regimes, but it should be included in the case of Mn doped GaN that is expected to be in the IB regime., Comment: 8 pages, 4 Postscript figures, RevTex4

We present a numerical study of the magnetization process of frustrated quantum spin-S chains with S=1, 3/2, 2 as well as the classical limit. Using the exact diagonalization and density-matrix renormalization techniques, we provide evidence that a plateau at one third of the saturation magnetization exists in the magnetization curve of frustrated spin-S chains with S>1/2. Similar to the case of S=1/2, this plateau state breaks the translational symmetry of the Hamiltonian and realizes an up-up-down pattern in the spin component parallel to the external field. Our study further shows that this plateau exists both in the cases of an isotropic exchange and in the easy-axis regime for spin-S=1, 3/2, and 2, but is absent in classical frustrated spin chains with isotropic interactions. We discuss the magnetic phase diagram of frustrated spin-1 and spin-3/2 chains as well as other emergent features of the magnetization process such as kink singularities, jumps, and even-odd effects. A quantitative comparison of the one-third plateau in the easy-axis regime between spin-1 and spin-3/2 chains on the one hand and the classical frustrated chain on the other hand indicates that the critical frustration and the phase boundaries of this state rapidly approach the classical result as the spin S increases., Comment: 15 pages RevTex4, 13 figures

We investigate in detail, using both analytical and numerical tools, the decoherence of electron spins in quantum dots (QDs) coupled to a bath of nuclear spins in magnetic fields or with various initial bath polarizations, focusing on the longitudinal relaxation in low and moderate field/polarization regimes. An increase of the initial polarization of nuclear spin bath has the same effect on the decoherence process as an increase of the external magnetic field, namely, the decoherence dynamics changes from smooth decay to damped oscillations. This change can be observed experimentally for a single QD and for a double-QD setup. Our results indicate that substantial increase of the decoherence time requires very large bath polarizations, and the use of other methods (dynamical decoupling or control of the nuclear spins distribution) may be more practical for suppressing decoherence of QD-based qubits., Comment: Rev. Tex, 5 pages, 3 eps color figures, submitted to Phys. Rev. B

We show that the symmetry of the spin zigzag chain E phase of the orthorhombic perovskite manganites and nickelates allows for the existence of a finite ferroelectric polarization. The proposed microscopic mechanism is independent of spin-orbit coupling. We predict that the polarization induced by the E-type magnetic order can potentially be enhanced by up to two orders of magnitude with respect to that in the spiral magnetic phases of TbMnO3 and similar multiferroic compounds., Comment: 4 pages, 2 figures, somewhat changed emphases, accepted to PRL

We investigate the ground state phase diagram of the half-filled repulsive Hubbard model in two dimensions in the presence of a staggered potential $\Delta$, the so-called ionic Hubbard model, using cluster dynamical mean field theory. We find that for large Coulomb repulsion, $U\gg \Delta$, the system is a Mott insulator (MI). For weak to intermediate values of $\Delta$, on decreasing $U$, the Mott gap closes at a critical value $U_{c1}(\Delta)$ beyond which a correlated insulating phase with possible bond order (BO) is found. Further, this phase undergoes a first-order transition to a band insulator (BI) at $U_{c2}(\Delta)$ with a finite charge gap at the transition. For large $\Delta$, there is a direct first-order transition from a MI to a BI with a single metallic point at the phase boundary.