Source author record

M. Jarrell

M. Jarrell 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

39works
8topics
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

39 published item(s)

preprint2020arXiv

Ab initio typical medium theory of substitutional disorder

By merging single-site typical medium theory with density functional theory we introduce a self-consistent framework for electronic structure calculations of materials with substitutional disorder which takes into account Anderson localization. The scheme and details of the implementation are presented and applied to the hypothetical alloy Li$_{c}$Be$_{1-c}$, and the results are compared with those obtained with the coherent potential approximation. Furthermore we demonstrate that Anderson localization suppresses ferromagnetic order for a very low concentration of (i) carbon impurities substituting oxygen in MgO$_{1-c}$C$_{c}$, and (ii) manganese impurities substituting magnesium in Mg$_{1-c}$Mn$_c$O for the low-spin magnetic configuration.

preprint2016arXiv

Generalized Multiband Typical Medium Dynamical Cluster Approximation: Application to (Ga,Mn)N

We generalize the multiband typical medium dynamical cluster approximation and the formalism introduced by Blackman, Esterling and Berk so that it can deal with localization in multiband disordered systems with both diagonal and off-diagonal disorder with complicated potentials. We also introduce a new ansatz for the angle resolved typical density of states that greatly improves the numerical stability of the method while preserving the independence of scattering events at different frequencies. Starting from the first-principles effective Hamiltonian, we apply this method to the diluted magnetic semiconductor Ga$_{1-x}$Mn$_x$N, and find the impurity band is completely localized for Mn concentrations $x<0.03$ while for $0.03 <x<0.10$ the impurity band has delocalized states but the chemical potential resides at or above the mobility edge. So, the system is always insulating within the experimental compositional limit ($x\approx 0.10$) due to Anderson localization. However, for $0.03 <x<0.10$ hole doping could make the system metallic allowing double exchange mediated, or enhanced, ferromagnetism. The developed method is expected to have a large impact on first-principles studies of Anderson localization.

preprint2016arXiv

The Bose-Hubbard model on a triangular lattice with diamond ring-exchange

Ring-exchange interactions have been proposed as a possible mechanism for a Bose-liquid phase at zero temperature, a phase that is compressible with no superfluidity. Using the Stochastic Green Function algorithm (SGF), we study the effect of these interactions for bosons on a two-dimensional triangular lattice. We show that the supersolid phase, that is known to exist in the ground state for a wide range of densities, is rapidly destroyed as the ring-exchange interactions are turned on. We establish the ground-state phase diagram of the system, which is characterized by the absence of the expected Bose-liquid phase.

preprint2015arXiv

Competing magnetic states, disorder, and the magnetic character of Fe3Ga4

The physical properties of metamagnetic Fe$_3$Ga$_4$ single crystals are investigated to explore the sensitivity of the magnetic states to temperature, magnetic field, and sample history. The data reveal a moderate anisotropy in the magnetization and the metamagnetic critical field along with features in the specific heat at the magnetic transitions $T_1=68$ K and $T_2=360$ K. Both $T_1$ and $T_2$ are found to be sensitive to the annealing conditions of the crystals suggesting that disorder affects the competition between the ferromagnetic (FM) and antiferromagnetic (AFM) states. Resistivity measurements reveal metallic transport with a sharp anomaly associated with the transition at $T_2$. The Hall effect is dominated by the anomalous contribution which rivals that of magnetic semiconductors in magnitude ($-5 μΩ$ cm at 2 T and 350 K) and undergoes a change of sign upon cooling into the low temperature FM state. The temperature and field dependence of the Hall effect indicate that the magnetism is likely to be highly itinerant in character and that a significant change in the electronic structure accompanies the magnetic transitions. We observe a contribution from the topological Hall effect in the AFM phase suggesting a non-coplanar contribution to the magnetism. Electronic structure calculations predict an AFM ground state with a wavevector parallel to the crystallographic $c$-axis preferred over the experimentally measured FM state by $\approx$ 50 meV per unit cell. However, supercell calculations with a small density of Fe-antisite defects introduced tend to stabilize the FM over the AFM state indicating that antisite defects may be the cause of the sensitivity to sample synthesis conditions.

preprint2015arXiv

Finite Cluster Typical Medium Theory for Disordered Electronic Systems

We use the recently developed typical medium dynamical cluster (TMDCA) approach~[Ekuma \etal,~\textit{Phys. Rev. B \textbf{89}, 081107 (2014)}] to perform a detailed study of the Anderson localization transition in three dimensions for the Box, Gaussian, Lorentzian, and Binary disorder distributions, and benchmark them with exact numerical results. Utilizing the nonlocal hybridization function and the momentum resolved typical spectra to characterize the localization transition in three dimensions, we demonstrate the importance of both spatial correlations and a typical environment for the proper characterization of the localization transition in all the disorder distributions studied. As a function of increasing cluster size, the TMDCA systematically recovers the re-entrance behavior of the mobility edge for disorder distributions with finite variance, obtaining the correct critical disorder strengths, and shows that the order parameter critical exponent for the Anderson localization transition is universal. The TMDCA is computationally efficient, requiring only a small cluster to obtain qualitative and quantitative data in good agreement with numerical exact results at a fraction of the computational cost. Our results demonstrate that the TMDCA provides a consistent and systematic description of the Anderson localization transition.

preprint2015arXiv

Metal-Insulator-Transition in a Weakly interacting Disordered Electron System

The interplay of interactions and disorder is studied using the Anderson-Hubbard model within the typical medium dynamical cluster approximation. Treating the interacting, non-local cluster self-energy ($Σ_c[{\cal \tilde{G}}](i,j\neq i)$) up to second order in the perturbation expansion of interactions, $U^2$, with a systematic incorporation of non-local spatial correlations and diagonal disorder, we explore the initial effects of electron interactions ($U$) in three dimensions. We find that the critical disorder strength ($W_c^U$), required to localize all states, increases with increasing $U$; implying that the metallic phase is stabilized by interactions. Using our results, we predict a soft pseudogap at the intermediate $W$ close to $W_c^U$ and demonstrate that the mobility edge ($ω_ε$) is preserved as long as the chemical potential, $μ$, is at or beyond the mobility edge energy.

preprint2014arXiv

A Typical Medium Dynamical Cluster Approximation for the Study of Anderson Localization in Three Dimensions

We develop a systematic typical medium dynamical cluster approximation that provides a proper description of the Anderson localization transition in three dimensions (3D). Our method successfully captures the localization phenomenon both in the low and large disorder regimes, and allows us to study the localization in different momenta cells, which renders the discovery that the Anderson localization transition occurs in a cell-selective fashion. As a function of cluster size, our method systematically recovers the re-entrance behavior of the mobility edge and obtains the correct critical disorder strength for Anderson localization in 3D.

preprint2014arXiv

Dual-fermion approach to interacting disordered fermion systems

We generalize the recently introduced dual fermion (DF) formalism for disordered fermion systems by including the effect of interactions. For an interacting disordered system the contributions to the full vertex function have to be separated into elastic and inelastic scattering processes, and addressed differently when constructing the DF diagrams. By applying our approach to the Anderson-Falicov-Kimball model and systematically restoring the nonlocal correlations in the DF lattice calculation, we show a significant improvement over the Dynamical Mean-Field Theory and the Coherent Potential Approximation for both one-particle and two-particle quantities.

preprint2014arXiv

Electronic Structure and Spectra of CuO

We report the electronic structure of monoclinic CuO as obtained from first principles calculations utilizing density functional theory plus effective Coulomb interaction (DFT + U) method. In contrast to standard DFT calculations taking into account electronic correlations in DFT + U gave antiferromagnetic insulator with energy gap and magnetic moment values in good agreement with experimental data. The electronic states around the Fermi level are formed by partially filled Cu 3$d_{x^2-y^2}$ orbitals with significant admixture of O 2$p$ states. Theoretical spectra are calculated using DFT + U electronic structure method and their comparison with experimental photoemission and optical spectra show very good agreement.

preprint2014arXiv

Spectral changes in layered $f$-electron systems induced by Kondo hole substitution in the boundary-layer

We investigate the effect of disorder on the dynamical spectrum of layered $f$-electron systems. With random dilution of $f$-sites in a single Kondo insulating layer, we explore the range and extent to which Kondo hole incoherence can penetrate into adjacent layers. We consider three cases of neighboring layers: band insulator, Kondo insulator and simple metal. The disorder-induced spectral weight transfer, used here for quantification of the proximity effect, decays algebraically with distance from the boundary layer. Further, we show that the spectral weight transfer is highly dependent on the frequency range considered as well as the presence of interactions in the clean adjacent layers. The changes in the low frequency spectrum are very similar when the adjacent layers are either metallic or Kondo insulating, and hence are independent of interactions. In stark contrast, a distinct picture emerges for the spectral weight transfers across large energy scales. The spectral weight transfer over all energy scales is much higher when the adjacent layers are non-interacting as compared to when they are strongly interacting Kondo insulators. Thus, over all scales, interactions screen the disorder effects significantly. We discuss the possibility of a crossover from non-Fermi liquid to Fermi liquid behavior upon increasing the ratio of clean to disordered layers in particle-hole asymmetric systems.

preprint2014arXiv

Study of off-diagonal disorder using the typical medium dynamical cluster approximation

We generalize the typical medium dynamical cluster approximation (TMDCA) and the local Blackman, Esterling, and Berk (BEB) method for systems with off-diagonal disorder. Using our extended formalism we perform a systematic study of the effects of non-local disorder-induced correlations and of off-diagonal disorder on the density of states and the mobility edge of the Anderson localized states. We apply our method to the three-dimensional Anderson model with configuration dependent hopping and find fast convergence with modest cluster sizes. Our results are in good agreement with the data obtained using exact diagonalization, and the transfer matrix and kernel polynomial methods.

preprint2013arXiv

Mean-field embedding of the dual fermion approach for correlated electron systems

To reduce the rapidly growing computational cost of the dual fermion lattice calculation with increasing system size, we introduce two embedding schemes. One is the real fermion embedding, and the other is the dual fermion embedding. Our numerical tests show that the real fermion and dual fermion embedding approaches converge to essentially the same result. The application on the Anderson disorder and Hubbard models shows that these embedding algorithms converge more quickly with system size as compared to the conventional dual fermion method, for the calculation of both single-particle and two-particle quantities.

preprint2013arXiv

Phase Stability in the Two dimensional Anisotropic Boson Hubbard Hamiltonian

The two dimensional square lattice hard-core boson Hubbard model with near neighbor interactions has a `checkerboard' charge density wave insulating phase at half-filling and sufficiently large intersite repulsion. When doped, rather than forming a supersolid phase in which long range charge density wave correlations coexist with a condensation of superfluid defects, the system instead phase separates. However, it is known that there are other lattice geometries and interaction patterns for which such coexistence takes place. In this paper we explore the possibility that anisotropic hopping or anisotropic near neighbor repulsion might similarly stabilize the square lattice supersolid. By considering the charge density wave structure factor and superfluid density for different ratios of interaction strength and hybridization in the $\hat x$ and $\hat y$ directions, we conclude that phase separation still occurs.

preprint2013arXiv

Re-examining the electronic structure of germanium: A first-principle study

We report results from an efficient, robust, ab-initio method for self-consistent calculations of electronic and structural properties of Ge. Our non-relativistic calculations employed a generalized gradient approximation (GGA) potential and the linear combination of atomic orbitals (LCAO) formalism. The distinctive feature of our computations stem from the use of Bagayoko-Zhao-Williams-Ekuma-Franklin (BZW-EF) method. Our results are in agreement with experimental ones where the latter are available. In particular, our theoretical, indirect band gap of 0.65 eV, at the experimental lattice constant of 5.66 Å, is in excellent agreement with experiment. Our predicted, equilibrium lattice constant is 5.63 Å, with a corresponding indirect band gap of 0.65 eV and a bulk modulus of 80 GPa. We also calculated the effective masses in various directions with respect to the $Γ$ point.

preprint2012arXiv

Dual Fermion Method for Disordered Electronic Systems

While the coherent potential approximation (CPA) is the prevalent method for the study of disordered electronic systems, it fails to capture non-local correlations and Anderson localization. To incorporate such effects, we extend the dual fermion approach to disordered non-interacting systems using the replica method. Results for single- and two- particle quantities show good agreement with cluster extensions of the CPA; moreover, weak localization is captured. As a natural extension of the CPA, our method presents an alternative to the existing cluster theories. It can be used in various applications, including the study of disordered interacting systems, or for the description of non-local effects in electronic structure calculations.

preprint2012arXiv

Electronic, structural, and elastic properties of metal nitrides XN (X = Sc, Y): A first principle study

We utilized a simple, robust, first principle method, based on basis set optimization with the BZW-EF method, to study the electronic and related properties of transition metal mono-nitrides: ScN and YN. We solved the KS system of equations self-consistently within the linear combination of atomic orbitals (LCAO) formalism. It is shown that the band gap and low energy conduction bands, as well as elastic and structural properties, can be calculated with a reasonable accuracy when the LCAO formalism is used to obtain an optimal basis. Our calculated, indirect electronic band gap (E$^\mathrm{Γ-X}_g$) is 0.79 (LDA) and 0.88 eV (GGA) for ScN. In the case of YN, we predict an indirect band gap (E$^\mathrm{Γ-X}_g$) of 1.09 (LDA) and 1.15 eV (GGA). We also calculated the equilibrium lattice constants, the bulk moduli (B$_{o}$), effective masses, and elastic constants for both systems. Our calculated values are in excellent agreement with experimental ones where the latter are available.

preprint2012arXiv

First principle electronic, structural, elastic, and optical properties of strontium titanate

We report self-consistent ab-initio electronic, structural, elastic, and optical properties of cubic SrTiO$_{3}$ perovskite. Our non-relativistic calculations employed a generalized gradient approximation (GGA) potential and the linear combination of atomic orbitals (LCAO) formalism. The distinctive feature of our computations stem from solving self-consistently the system of equations describing the GGA, using the Bagayoko-Zhao-Williams (BZW) method. Our results are in agreement with experimental ones where the later are available. In particular, our theoretical, indirect band gap of 3.24 eV, at the experimental lattice constant of 3.91 Å, is in excellent agreement with experiment. Our predicted, equilibrium lattice constant is 3.92 Å, with a corresponding indirect band gap of 3.21 eV and bulk modulus of 183 GPa.

preprint2012arXiv

First Principle Local Density Approximation Description of the Electronic Properties of Ferroelectric Sodium Nitrite

The electronic structure of the ferroelectric crystal, NaNO$_2$, is studied by means of first-principles, local density calculations. Our ab-initio, non-relativistic calculations employed a local density functional approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). Following the Bagayoko, Zhao, Williams, method, as enhanced by Ekuma, and Franklin (BZW-EF), we solved self-consistently both the Kohn-Sham equation and the equation giving the ground state charge density in terms of the wave functions of the occupied states. We found an indirect band gap of 2.83 eV, from W to R. Our calculated direct gaps are 2.90, 2.98, 3.02, 3.22, and 3.51 eV at R, W, X, Γ, and T, respectively. The band structure and density of states show high localization, typical of a molecular solid. The partial density of states shows that the valence bands are formed only by complex anionic states. These results are in excellent agreement with experiment. So are the calculated densities of states. Our calculated electron effective masses of 1.18, 0.63, and 0.73 mo in the Γ-X, Γ-R, and Γ-W directions, respectively, show the highly anisotropic nature of this material.

preprint2012arXiv

Physical Properties of $Ba_2 Mn_2 Sb_2 O$ Single Crystals

We report both experimental and theoretical investigations of the physical properties of Ba$_\mathrm{2}$Mn$_\mathrm{2}$Sb$_\mathrm{2}$O single crystals. This material exhibits a hexagonal structure with lattice constants: a = 4.7029(15) Å and c = 19.9401(27) Å, as obtained from powder X-ray diffraction measurements, and in agreement with structural optimization through density functional theory (DFT) calculations. The magnetic susceptibility and specific heat show anomalies at T$_\mathrm{N}$ = 60 K, consistent with antiferromagnetic ordering. However, the magnitude of T$_\mathrm{N}$ is significantly smaller than the Curie-Weiss temperature ($\mid$$\mathrm{Θ_{CW}}$$\mid$ $\approx$ 560 K), suggesting a magnetic system of reduced dimensionality. The temperature dependence of both the in-plane and out-of-plane resistivity changes from an activated at $T$ $>$ T$_\mathrm{x}$ $\sim$ 200 K to a logarithmic at $T$ $<$ T$_\mathrm{x}$. Correspondingly, the magnetic susceptibility displays a bump at T$_\mathrm{x}$. DFT calculations at the DFT + U level support the experimental observation of an antiferromagnetic ground state.

preprint2011arXiv

Dual Fermion Dynamical Cluster Approach for Strongly Correlated Systems

We have designed a new multi-scale approach for Strongly Correlated Systems by combining the Dynamical Cluster Approximation (DCA) and the recently introduced dual-fermion formalism. This approach employs an exact mapping from a real lattice to a DCA cluster of linear size Lc embedded in a dual fermion lattice. Short-length-scale physics is addressed by the DCA cluster calculation, while longer-length-scale physics is addressed diagrammatically using dual fermions. The bare and dressed dual Fermionic Green functions scale as O(1/Lc) so perturbation theory on the dual lattice converges very quickly. E.g., the dual Fermion self-energy calculated with simple second order perturbation theory is of order O(1/Lc^3), with third order and three body corrections down by an additional factor of O(1/Lc^2).

preprint2011arXiv

Role of the van Hove Singularity in the Quantum Criticality of the Hubbard Model

A quantum critical point (QCP), separating the non-Fermi liquid region from the Fermi liquid, exists in the phase diagram of the 2D Hubbard model [Vidhyadhiraja et. al, Phys. Rev. Lett. 102, 206407 (2009)]. Due to the vanishing of the critical temperature associated with a phase separation transition, the QCP is characterized by a vanishing quasiparticle weight. Near the QCP, the pairing is enhanced since the real part of the bare d-wave p-p susceptibility exhibits algebraic divergence with decreasing temperature, replacing the logarithmic divergence found in a Fermi liquid [Yang et. al, Phys. Rev. Lett. 106, 047004 (2011)]. In this paper we explore the single-particle and transport properties near the QCP. We focus mainly on a van Hove singularity (vHS) coming from the relatively flat dispersion that crosses the Fermi level near the quantum critical filling. The flat part of the dispersion orthogonal to the antinodal direction remains pinned near the Fermi level for a range of doping that increases when we include a negative next-near-neighbor hopping t' in the model. For comparison, we calculate the bare d-wave pairing susceptibility for non-interacting models with the usual two-dimensional tight binding dispersion and a hypothetical quartic dispersion. We find that neither model yields a vHS that completely describes the critical algebraic behavior of the bare d-wave pairing susceptibility. The resistivity, thermal conductivity, thermopower, and the Wiedemann-Franz Law are examined in the Fermi liquid, marginal Fermi liquid, and pseudo-gap doping regions. A negative next-near-neighbor hopping t' increases the doping region with marginal Fermi liquid character. Both T and negative t' are relevant variables for the QCP, and both the transport and the motion of the vHS with filling suggest that they are qualitatively similar in their effect.

preprint2011arXiv

Solving the Parquet Equations for the Hubbard Model beyond Weak Coupling

We find that imposing the crossing symmetry in the iteration process considerably extends the range of convergence for solutions of the parquet equations for the Hubbard model. When the crossing symmetry is not imposed, the convergence of both simple iteration and more complicated continuous loading (homotopy) methods are limited to high temperatures and weak interactions. We modify the algorithm to impose the crossing symmetry without increasing the computational complexity. We also imposed time reversal and a subset of the point group symmetries, but they did not further improve the convergence. We elaborate the details of the latency hiding scheme which can significantly improve the performance in the computational implementation. With these modifications, stable solutions for the parquet equations can be obtained by iteration more quickly even for values of the interaction that are a significant fraction of the bandwidth and for temperatures that are much smaller than the bandwidth. This may represent a crucial step towards the solution of two-particle field theories for correlated electron models.

preprint2011arXiv

Using off-diagonal confinement as a cooling method

In a recent letter [Phys. Rev. Lett. 104, 167201 (2010)] we proposed a new confining method for ultracold atoms on optical lattices, based on off-diagonal confinement (ODC). This method was shown to have distinct advantages over the conventional diagonal confinement (DC) that makes use of a trapping potential, including the existence of pure Mott phases and highly populated condensates. In this paper we show that the ODC method can also lead to temperatures that are smaller than with the conventional DC method, depending on the control parameters. We determine these parameters using exact diagonalizations for the hard-core case, then we extend our results to the soft-core case by performing quantum Monte Carlo (QMC) simulations for both DC and ODC systems at fixed temperatures, and analysing the corresponding entropies. We also propose a method for measuring the entropy in QMC simulations.

preprint2010arXiv

A dynamical mean-field approximation study of a tight-binding model for Ga$_{1-x}$Mn$_{x}$As

The magneto-optical properties of the ferromagnetic semiconductor Ga$_{1-x}$Mn$_{x}$As are studied within the dynamical mean-field approximation (DMFA). A material-specific multiband $sp^{3}$ tight-binding Hamiltonian is employed for the dispersion of the GaAs host. The calculated density of states shows an impurity band and a distorted valence band for large and moderate values of magnetic coupling, respectively. Upon using the more realistic band structure, the ferromagnetic transition temperature is significantly closer to the experimental results than the previous predictions of $k\cdot p$ models. The optical conductivity shows a Drude-like peak at low frequencies which is suppressed by increasing of the magnetic coupling.

preprint2010arXiv

Cluster solver for dynamical mean-field theory with linear scaling in inverse temperature

Dynamical mean field theory and its cluster extensions provide a very useful approach for examining phase transitions in model Hamiltonians, and, in combination with electronic structure theory, constitute powerful methods to treat strongly correlated materials. The key advantage to the technique is that, unlike competing real space methods, the sign problem is well controlled in the Hirsch-Fye (HF) quantum Monte Carlo used as an exact cluster solver. However, an important computational bottleneck remains; the HF method scales as the cube of the inverse temperature, $β$. This often makes simulations at low temperatures extremely challenging. We present here a new method based on determinant quantum Monte Carlo which scales linearly in $β$, and demonstrate that the sign problem is identical to HF.

preprint2010arXiv

Phonons and the coherence scale of models of heavy fermions

We consider models of heavy fermions in the strong coupling or local moment limit and include phonon degrees of freedom on the conduction electrons. Due to the large mass or low coherence temperature of the heavy fermion state, it is shown that such a regime is dominated by vertex corrections which leads to the complete failure of the Migdal theorem. Even at weak electron-phonon couplings, binding of the conduction electrons competes with the Kondo effect and substantially reduces the coherence temperature, ultimately leading to the Kondo breakdown. Those results are obtained using a combination of the slave boson method and Migdal-Eliashberg approximation as well as the dynamical mean-field theory approximation.

preprint2010arXiv

Pure Mott phases in confined ultracold atomic systems

We propose a novel scheme for confining atoms to optical lattices by engineering a spatially-inhomogeneous hopping matrix element in the Hubbard-model (HM) description, a situation we term off-diagonal confinement (ODC). We show, via an exact numerical solution of the boson HM with ODC, that this scheme possesses distinct advantages over the conventional method of confining atoms using an additional trapping potential, including the presence of incompressible Mott phases at commensurate filling and a phase diagram that is similar to the uniform HM. The experimental implementation of ODC will thus allow a more faithful realization of correlated phases of interest in cold atom experiments.

preprint2010arXiv

Quantum Criticality and Incipient Phase Separation in the Thermodynamic Properties of the Hubbard Model

Transport measurements on the cuprates suggest the presence of a quantum critical point hiding underneath the superconducting dome near optimal hole doping. We provide numerical evidence in support of this scenario via a dynamical cluster quantum Monte Carlo study of the extended two-dimensional Hubbard model. Single particle quantities, such as the spectral function, the quasiparticle weight and the entropy, display a crossover between two distinct ground states: a Fermi liquid at low filling and a non-Fermi liquid with a pseudogap at high filling. Both states are found to cross over to a marginal Fermi-liquid state at higher temperatures. For finite next-nearest-neighbor hopping t' we find a classical critical point at temperature T_c. This classical critical point is found to be associated with a phase separation transition between a compressible Mott gas and an incompressible Mott liquid corresponding to the Fermi liquid and the pseudogap state, respectively. Since the critical temperature T_c extrapolates to zero as t' vanishes, we conclude that a quantum critical point connects the Fermi-liquid to the pseudogap region, and that the marginal-Fermi-liquid behavior in its vicinity is the analogous of the supercritical region in the liquid-gas transition.

preprint2010arXiv

Quantum Criticality Due to Incipient Phase Separation in the Two-dimensional Hubbard Model

We investigate the two-dimensional Hubbard model with next-nearest-neighbor hopping, t', using the dynamical cluster approximation. We confirm the existence of a first-order phase-separation transition terminating at a second-order critical point at filling n_c(t') and temperature T_ps(t'). We find that as t' approaches zero, T_ps(t') vanishes and n_c(t') approaches the filling associated with the quantum critical point separating the Fermi liquid from the pseudogap phase. We propose that the quantum critical point under the superconducting dome is the zero-temperature limit of the line of second-order critical points.

preprint2010arXiv

Quantum Monte Carlo study of the two-dimensional fermion Hubbard Model

We report large scale determinant Quantum Monte Carlo calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian at half-filling. The sharp Fermi surface of the non-interacting limit is significantly broadened by the electronic correlations, but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak coupling to the strong coupling Heisenberg limit. Our lattices provide improved resolution of the Green's function in momentum space, allowing a more quantitative comparison with time-of-flight optical lattice experiments.

preprint2010arXiv

The effect of spin-orbit interaction and attractive Coulomb potential on the magnetic properties of Ga$_{1-x}$Mn$_{x}$As

We employ the dynamical mean-field approximation to study the magnetic properties of a model relevant for the dilute magnetic semiconductors. Our model includes the spin-orbit coupling on the hole bands, the exchange interaction, and the attractive Coulomb potential between the negatively charged magnetic ions and the itinerant holes. The inclusion of the Coulomb potential significantly renormalizes the exchange coupling and enhances the ferromagnetic transition temperature for a wide range of couplings. We also explore the effect of the spin-orbit interaction by using two different values of the ratio of the effective masses of the heavy and light holes. We show that in the regime of small $J_{c}$-$V$ the spin-orbit interaction enhances $T_{c}$, while for large enough values of $J_{c}$-$V$ magnetic frustration reduces $T_c$ to values comparable to the previously calculated strong coupling limit.

preprint2010arXiv

Validity of the spin-susceptibility "glue" approximation for pairing in the two-dimensional Hubbard model

We examine the validity of the weak coupling spin-susceptibility "glue" approximation (SSGA) in a two-dimensional Hubbard model for cuprates. For comparison, we employ the well-established dynamical cluster approximation (DCA) with a quantum Monte Carlo algorithm as a cluster solver. We compare the leading eigenvalues and corresponding eigenfunctions of the DCA and SSGA pairing matrices. For realistic model parameters, we find that the SSGA fails to capture the leading pairing symmetries seen in the DCA. Furthermore, when the SSGA is improved through the addition of a term with d-wave symmetry, the strength of this additional term is found to be larger than that of the "glue" approximation.

preprint2009arXiv

Analytic Continuation of Quantum Monte Carlo Data by Stochastic Analytical Inference

We present an algorithm for the analytic continuation of imaginary-time quantum Monte Carlo data which is strictly based on principles of Bayesian statistical inference. Within this framework we are able to obtain an explicit expression for the calculation of a weighted average over possible energy spectra, which can be evaluated by standard Monte Carlo simulations, yielding as by-product also the distribution function as function of the regularization parameter. Our algorithm thus avoids the usual ad-hoc assumptions introduced in similar algortihms to fix the regularization parameter. We apply the algorithm to imaginary-time quantum Monte Carlo data and compare the resulting energy spectra with those from a standard maximum entropy calculation.

preprint2009arXiv

Parquet approximation for the 4x4 Hubbard cluster

We present a numerical solution of the parquet approximation (PA), a conserving diagrammatic approach which is self-consistent at both the single-particle and the two-particle levels. The fully irreducible vertex is approximated by the bare interaction thus producing the simplest approximation that one can perform with the set of equations involved in the formalism. The method is applied to the Hubbard model on a half-filled 4x4 cluster. Results are compared to those obtained from Determinant Quantum Monte Carlo (DQMC), FLuctuation EXchange (FLEX), and self-consistent second-order approximation methods. This comparison shows a satisfactory agreement with DQMC and a significant improvement over the FLEX or the self-consistent second-order approximation.

preprint2009arXiv

The screening of 4f moments and delocalization in the compressed light rare earths

Spin and charge susceptibilities and the 4f^n, 4f^{n-1}, and 4f^{n+1} configuration weights are calculated for compressed Ce (n=1), Pr (n=2), and Nd (n=3) metals using dynamical mean field theory combined with the local-density approximation. At ambient and larger volumes these trivalent rare earths are pinned at sharp 4f^n configurations, their 4f moments assume atomic-limiting values, are unscreened, and the 4f charge fluctuations are small indicating little f state density near the Fermi level. Under compresssion there is dramatic screening of the moments and an associated increase in both the 4f charge fluctuations and static charge susceptibility. These changes are coincident with growing weights of the 4f^{n-1} configurations, which it is argued are better measures of delocalization than the 4f^{n+1} weights which are compromised by an increase in the number of 4f electrons caused by rising 6s, 6p bands. This process is continuous and prolonged as a function of volume, with strikingly similarity among the three rare earths, aside from the effects moderating and shifting to smaller volumes for the heavier members. The observed alpha-gamma collapse in Ce occurs over the large-volume half of this evolution, the Pr analog at smaller volumes, and Nd has no collapse.

preprint2006arXiv

Multi-scale Extensions to Quantum Cluster Methods for Strongly Correlated Electron Systems

A numerically implementable Multi-scale Many-Body approach to strongly correlated electron systems is introduced. An extension to quantum cluster methods, it approximates correlations on any given length-scale commensurate with the strength of the correlations on the respective scale. Short length-scales are treated explicitly, long ones are addressed at a mean-field level and intermediate length-regime correlations are assumed to be weak and are approximated diagrammatically. To illustrate and test this method, we apply it to the one dimensional Hubbard model. The resulting multi-scale self-energy provides a very good quantitative agreement with substantially more numerically expensive, explicit Quantum Monte-Carlo calculations.

preprint1995arXiv

Anomalous Normal-State Properties of High-T$_c$ Superconductors -- Intrinsic Properties of Strongly Correlated Electron Systems?

A systematic study of optical and transport properties of the Hubbard model, based on Metzner and Vollhardt's dynamical mean-field approximation, is reviewed. This model shows interesting anomalous properties that are, in our opinion, ubiquitous to single-band strongly correlated systems (for all spatial dimensions greater than one), and also compare qualitatively with many anomalous transport features of the high-T$_c$ cuprates. This anomalous behavior of the normal-state properties is traced to a ``collective single-band Kondo effect,'' in which a quasiparticle resonance forms at the Fermi level as the temperature is lowered, ultimately yielding a strongly renormalized Fermi liquid at zero temperature.

preprint1995arXiv

Ground-State Dynamical Correlation Functions: An Approach from Density Matrix Renormalization Group Method

A numerical approach to ground-state dynamical correlation functions from Density Matrix Renormalization Group (DMRG) is developed. Using sum rules, moments of a dynamic correlation function can be calculated with DMRG, and with the moments the dynamic correlation function can be obtained by the maximum entropy method. We apply this method to one-dimensional spinless fermion system, which can be converted to the spin 1/2 Heisenberg model in a special case. The dynamical density-density correlation function is obtained.

preprint1995arXiv

Spectral Properties and Bandstructure of Correlated Electron Systems

We present $\vec{k}$-dependent one-particle spectra and corresponding effective bandstructures for the $2d$ Hubbard model calculated within the dynamical molecular field theory (DMFT). This method has proven to yield highly nontrivial results for a variety of quantities but the question remains open to what extent it is applicable to relevant physical situations. To address this problem we compare our results for spectral functions to those obtained by QMC simulations. The good agreement supports our notion that the DMFT is indeed a sensible ansatz for correlated models even in to $d=2$.