Source author record

J. P. F. LeBlanc

J. P. F. LeBlanc 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

24works
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

24 published item(s)

preprint2026arXiv

NonHermitian Topological Phases in a Hermitian Modified Bosonic Kitaev Chain

We present a modification to the bosonic Kitaev chain that, despite being Hermitian, supports both nonHermitian skin effect and nontrivial topological edge modes in its excitation Hamiltonian. We establish an exact mapping between the excitation Hamiltonian of our system and a nonHermitian Su-Schrieffer-Heeger (SSH) model, which allows for a completely analytical characterization of its topology. In particular, topological phase transition points separating a topologically trivial and nontrivial regime were identified analytically by the appropriate winding number invariant and the presence of zero energy modes. Similarly to the regular bosonic Kitaev chain, the nonHermitian skin effect and some (but not all) topological edge modes are quickly destroyed at nonzero bosonic onsite potential (harmonic oscillator frequency). Remarkably, however, disorder partially recovers some of these features. This work thus demonstrates the potential of a modified bosonic Kitaev chain as a platform to generate rich nonHermitian topological phenomena from a completely Hermitian system's perspective. Lastly, we suggest a possible experimental realization of the model, which could allow for total control over the parameter space.

preprint2022arXiv

One- and two-particle properties of the weakly interacting two-dimensional Hubbard model in proximity to the van Hove singularity

We study the weak-coupling limit of the $t-t^\prime-U$ Hubbard model on a two-dimensional square lattice using a direct perturbative approach. Aided by symbolic computational tools, we compute the longitudinal density-density correlation functions in the $χ_{\uparrow \uparrow}$ and $χ_{\uparrow \downarrow}$ basis from which we can obtain the dynamical spin and charge susceptibilities at arbitrary doping and temperature. We find that for non-zero $t^\prime$, the zero frequency commensurate $\mathbf{q} = (π, π)$ spin and charge excitations are each strongest at different densities and we observe a clear behavioral change that appears tied to the van Hove singularity of the non-interacting dispersion upon which the perturbative expansion is built. We find a strongly reduced compressibility in the vicinity of the van Hove singularity as well as a behavioral change in the double occupancy. For finite $t^\prime$, the observed van Hove singularity occurs away from half-filling leading us to conclude that that this reduction in compressibility is distinct from Mott insulating physics that one expects in the strong-coupling regime. We compute the full dynamical spin and charge excitations and observe distinct structure for electron and hole doped scenarios in agreement with experiments on cuprate materials. Finally, we observe a peculiar splitting in spin and charge excitations in the vicinity of the van Hove singularity, the origin of which is traced to a splitting near the bottom of the band.

preprint2020arXiv

Algorithmic approach to diagrammatic expansions for real-frequency evaluation of susceptibility functions

We systematically generate the perturbative expansion for the two-particle spin susceptibility in the Feynman diagrammatic formalism and apply this expansion to a model system - the single-band Hubbard model on a square lattice. We make use of algorithmic Matsubara integration (AMI) [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 99, 035120 (2019)] to analytically evaluate Matsubara frequency summations, allowing us to symbolically impose analytic continuation to the real frequency axis. We minimize our computational expense by applying graph invariant transformations [Amir Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 101, 125109 (2020)]. We highlight extensions of the random-phase approximation and T-matrix methods that, due to AMI, become tractable. We present results for weak interaction strength where the direct perturbative expansion is convergent, and verify our results on the Matsubara axis by comparison to other numerical methods. By examining the spin susceptibility as a function of real-frequency via an order-by-order expansion we can identify precisely what role higher order corrections play on spin susceptibility and demonstrate the utility and limitations of our approach.

preprint2020arXiv

Extended crossover from Fermi liquid to quasi-antiferromagnet in the half-filled 2D Hubbard model

The ground state of the Hubbard model with nearest-neighbor hopping on the square lattice at half filling is known to be that of an antiferromagnetic (AFM) band insulator for any on-site repulsion. At finite temperature, the absence of long-range order makes the question of how the interaction-driven insulator is realized nontrivial. We address this problem with controlled accuracy in the thermodynamic limit using self-energy diagrammatic determinant Monte Carlo and dynamical cluster approximation methods and show that development of long-range AFM correlations drives an extended crossover from Fermi liquid to insulating behavior in the parameter regime that precludes a metal-to-insulator transition. The intermediate crossover state is best described as a non-Fermi liquid with a partially gapped Fermi surface.

preprint2019arXiv

Fluctuation diagnostics of the finite temperature quasi-antiferromagnetic regime of the 2D Hubbard model

We study the finite temperature Fermi-liquid to non-Fermi-liquid crossover in the 2D Hubbard model for a range of dopings using the self-consistent ladder dual fermion method. We consider relatively high temperatures where we identify a suppression of the density of states near the Fermi level caused by a quasi-antiferromagnetic behaviour that is itself characterized by a long, but finite, correlation length scale. We perform fluctuation diagnostics to decompose the single-particle self energy into scattering $q$-vector and bosonic frequency contributions. Within this framework we find that the key contributions to the single-particle self energy that give non-Fermi-liquid character, even at weak coupling, are caused by relatively sharp $q=(π,π)$ spin fluctuations, while the decomposition in the bosonic frequency channel shows a complicated dependence on the relative strengths of zero, positive and negative frequency contributions. Finally, variation in density suggests that the tendency towards non-Fermi-liquid behavior is not substantially different for electron or hole doped systems.

preprint2019arXiv

Optimal grouping of arbitrary diagrammatic expansions via analytic pole structure

We present a general method to optimize the evaluation of Feynman diagrammatic expansions, which requires the automated symbolic assignment of momentum/energy conserving variables to each diagram. With this symbolic representation, we utilize the pole structure of each diagram to automatically sort the Feynman diagrams into groups that are likely to contain nearly equal or nearly cancelling diagrams, and we show that for some systems this cancellation is exact. This allows for a potentially massive cancellation during the numerical integration of internal momenta variables, leading to an optimal suppression of the `sign problem' and hence reducing the computational cost. Although we define these groups using a frequency space representation, the equality or cancellation of diagrams within the group remains valid in other representations such as imaginary time used in standard diagrammatic Monte Carlo. As an application of the approach we apply this method, combined with algorithmic Matsubara integration (AMI) [Phys. Rev. B 99, 035120 (2019)] and Monte Carlo methods, to the Hubbard model self-energy expansion on a 2D square lattice up to sixth order which we evaluate and compare with existing benchmarks.

preprint2016arXiv

Parquet decomposition calculations of the electronic self-energy

The parquet decomposition of the self-energy into classes of diagrams, those associated with specific scattering processes, can be exploited for different scopes. In this work, the parquet decomposition is used to unravel the underlying physics of non-perturbative numerical calculations. We show the specific example of dynamical mean field theory (DMFT) and its cluster extensions (DCA) applied to the Hubbard model at half-filling and with hole doping: These techniques allow for a simultaneous determination of two-particle vertex functions and self-energies, and hence, for an essentially "exact" parquet decomposition at the single-site or at the cluster level. Our calculations show that the self-energies in the underdoped regime are dominated by spin scattering processes, consistent with the conclusions obtained by means of the fluctuation diagnostics approach [Phys. Rev. Lett. 114, 236402 (2015)]. However, differently from the latter approach, the parquet procedure displays important changes with increasing interaction: Even for relatively moderate couplings, well before the Mott transition, singularities appear in different terms, with the notable exception of the predominant spin-channel. We explain precisely how these singularities, which partly limit the utility of the parquet decomposition, and - more generally - of parquet-based algorithms, are never found in the fluctuation diagnostics procedure. Finally, by a more refined analysis, we link the occurrence of the parquet singularities in our calculations to a progressive suppression of charge fluctuations and the formation of an RVB state, which are typical hallmarks of a pseudogap state in DCA.

preprint2015arXiv

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

preprint2015arXiv

Superconducting Fluctuations in the Normal State of the Two-Dimensional Hubbard Model

We compute the two-particle quantities relevant for superconducting correlations in the two-dimensional Hubbard model within the dynamical cluster approximation. In the normal state we identify the parameter regime in density, interaction, and second-nearest-neighbor hopping strength that maximizes the $d_{x^2-y^2}$ superconducting transition temperature. We find in all cases that the optimal transition temperature occurs at intermediate coupling strength, and is suppressed at strong and weak interaction strengths. Similarly, superconducting fluctuations are strongest at intermediate doping and suppressed towards large doping and half-filling. We find a change in sign of the vertex contributions to $d_{xy}$ superconductivity from repulsive near half filling to attractive at large doping. $p$-wave superconductivity is not found at the parameters we study, and $s$-wave contributions are always repulsive. For negative second-nearest-neighbor hopping the optimal transition temperature shifts towards the electron-doped side in opposition to the van Hove singularity which moves towards hole doping. We surmise that an increase of the local interaction of the electron-doped compounds would increase $T_c$.

preprint2014arXiv

Dielectric screening of surface states in a topological insulator

Hexagonal warping provides an anisotropy to the dispersion curves of the helical Dirac fermions that exist at the surface of a topological insulator. A sub-dominant quadratic in momentum term leads to an asymmetry between conduction and valence band. A gap can also be opened through magnetic doping. We show how these various modifications to the Dirac spectrum change the polarization function of the surface states and employ our results to discuss their effect on the plasmons. In the long wavelength limit, the plasmon dispersion retains its square root dependence on its momentum, $\boldsymbol{q}$, but its slope is modified and it can acquire a weak dependence on the direction of $\boldsymbol{q}$. Further, we find the existence of several plasmon branches, one which is damped for all values of $\boldsymbol{q}$, and extract the plasmon scattering rate for a representative case.

preprint2014arXiv

Fluctuation diagnostics of the electron self-energy: Origin of the pseudogap physics

We demonstrate how to identify which physical processes dominate the low-energy spectral functions of correlated electron systems. We obtain an unambiguous classification through an analysis of the equation of motion for the electron self-energy in its charge, spin and particle-particle representations. Our procedure is then employed to clarify the controversial physics responsible for the appearance of the pseudogap in correlated systems. We illustrate our method by examining the attractive and repulsive Hubbard model in two-dimensions. In the latter, spin fluctuations are identified as the origin of the pseudogap, and we also explain why $d-$wave pairing fluctuations play a marginal role in suppressing the low-energy spectral weight, independent of their actual strength.

preprint2014arXiv

Signatures of a momentum independent pseudogap in the electronic density of states and Raman spectroscopy of the underdoped cuprates

We propose a hybridization phenomenology to describe the pseudogap state of the underdoped cuprates. We show how a momentum independent pseudogap opens asymmetrically from the Fermi-surface but symmetric to the zeroes of the hybridized bonding dispersion, which results in false d-wave characteristics of the pseudogap at the Fermi level. By comparing against a d-wave form factor we illustrate the difficulty in identifying a momentum independent order in momentum averaged quantities such as the electronic Raman response. We identify a suppression in the single-particle density of states which produces a hump feature which should be observable experimentally in tunnelling $dI/dV$ spectra and distinguishes the s-wave and d-wave ordering scenarios.

preprint2013arXiv

Effect of Dynamical Screening on Single Particle Spectral Features of Uniaxially Strained Graphene: Tuning the Plasmaron Ring

Electronic screening renormalizes the linear bands of graphene and in the vicinity of the Dirac point, creates a diamond shaped structure in the quasiparticle spectral density. This is a result of electron-plasmon scattering processes which produce a finite momentum feature referred to as the `plasmaron ring'. In this work we explore the effects of uniaxial strain on these spectral features with the aim of understanding how strain modifies correlations. We derive and calculate the spectral density to the G$_0$W-RPA level which allows us to identify the dispersive behaviour of the diamond geometry, and thus electron-plasmon scattering, for variation in electron-electron coupling strength and magnitude of applied strain. We find that the application of strain changes the geometry (in momentum) of the electron-plasmon scattering and that renormalizations beyond simple geometrical scalings further enhance this effect. These results suggest that the properties of the plasmaron ring can be tuned through the application of uniaxial strain, effectively producing a larger fine structure constant without the need to change the sample substrate.

preprint2013arXiv

Equation of State of the Fermionic 2D Hubbard Model

We present results for the equation of state of the two-dimensional Hubbard model on an isotropic square lattice as obtained from a controlled and numerically exact large-cluster dynamical mean field simulation. Our results are obtained for large but finite systems and are extrapolated to infinite system size using a known finite size scaling relation. We present the energy, entropy, double occupancy and nearest-neighbour spin correlations extrapolated to the thermodynamic limit and discuss the implications of these calculations on pseudogap physics of the 2D-Hubbard model away from half filling. We find a strong behavioural shift in energy below a temperature $T^*$ which becomes more pronounced for larger clusters. Finally, we provide reference calculations and tables for the equation of state for values of doping away from half filling which are of interest to cold atom experiments.

preprint2012arXiv

Distinguishing Coulomb and electron-phonon interactions for massless Dirac fermions

While many physical properties of graphene can be understood qualitatively on the basis of bare Dirac bands, there is specific evidence that electron-electron (EE) and electron-phonon (EP) interactions can also play an important role. We discuss strategies for extracting separate images of the EE and EP interactions as they present themselves in the electron spectral density and related self-energies. While for momentum, $k$, equal to its Fermi value, $k_F$, a composite structure is obtained which can be difficult to separate into its two constituent parts, at smaller values of $k$ the spectral function shows distinct incoherent sidebands on the left and right of the main quasiparticle line. These image respectively the EE and EP interactions, each being most prominent in its own energy window. We employ a maximum entropy inversion technique on the self energy to reveal the electron-phonon spectral density separate from the excitation spectrum due to coulomb correlations. Our calculations show that this technique can provide important new insights into inelastic scattering processes in graphene.

preprint2012arXiv

Emergence of Plasmaronic Structure in the Near Field Optical Response of Graphene

The finite momentum optical response $σ({\boldsymbol{q}},ω)$ of graphene can be probed with the innovative technique of infrared nanoscopy where mid-infrared radiation is confined by an atomic force microscope cantilever tip. In contrast to conventional $q\sim 0$ optical absorption which primarily involves Dirac fermions with momentum near the Fermi momentum, $k\sim k_F$, for finite $q$, $σ({\boldsymbol{q}},ω)$ has the potential to provide information on many body renormalizations and collective phenomena which have been found at small $k< k_F$ near the Dirac point in electron-doped graphene. For electron-electron interactions, the low energy excitation spectrum characterizing the incoherent part of the quasiparticle spectral function of Dirac electrons with $k\sim k_F$ consists of a flat, small amplitude background which scales with chemical potential and Fermi momentum. However, probing of the states with $k$ near $k=0$ will reveal plasmarons, a collective state of a charge carrier and a plasmon. These collective modes in graphene have recently been seen in angle-resolved photoemission spectroscopy and here we describe how they manifest in near field optics.

preprint2012arXiv

Impact of Electron-Phonon Coupling on Near-Field Optical Spectra

The finite momentum transfer ($\boldsymbol{q}$) longitudinal optical response $σ^L(\boldsymbol{q},ω)$ of graphene has a peak at an energy $ω=\hbar v_F q$. This corresponds directly to a quasiparticle peak in the spectral density at momentum relative to the Fermi momentum $k_F -q$. Inclusion of coupling to a phonon mode at $ω_E$ results, for $ω<|ω_E|$, in a constant electron-phonon renormalization of the bare bands by a mass enhancement factor $(1+λ)$ and this is followed by a phonon kink at $ω_E$ where additional broadening begins. Here we study the corresponding changes in the optical quasiparticle peaks which we find to continue to directly track the renormalized quasiparticle energies until $q$ is large enough that the optical transitions begin to sample the phonon kink region of the dispersion curves where linearity in momentum is lost in the renormalized Dirac Fermion dispersion curves and the correspondence to a single quasiparticle energy is lost. Nevertheless there remains in $σ^L(\boldsymbol{q},ω)$ features analogous to the phonon kinks of the dispersion curves which are observable through variation of $q$ and $ω$.

preprint2012arXiv

Optical Self Energy in Graphene due to Correlations

In highly correlated systems one can define an optical self energy in analogy to its quasiparticle (QP) self energy counterpart. This quantity provides useful information on the nature of the excitations involved in inelastic scattering processes. Here we calculate the self energy of the intraband optical transitions in graphene originating in the electron-electron interaction (EEI) as well as electron-phonon interaction (EPI). Although optics involves an average over all momenta ($k$) of the charge carriers, the structure in the optical self energy is nevertheless found to mirror mainly that of the corresponding quasiparticles for $k$ equal to or near the Fermi momentum $k_F$. Consequently plasmaronic structures which are associated with momenta near the Dirac point at $k=0$ are not important in the intraband optical response. While the structure of the electron-phonon interaction (EPI) reflects the sharp peaks of the phonon density of states, the excitation spectrum associated with the electron-electron interaction is in comparison structureless and flat and extends over an energy range which scales linearly with the value of the chemical potential. Modulations seen on the edge of the interband optical conductivity as it rises towards its universal background value are traced to structure in the quasiparticle self energies around $k_F$ of the lower Dirac cone associated with the occupied states.

preprint2011arXiv

Effect of electron-phonon coupling on energy and density of states renormalizations of dynamically screened graphene

Electronic screening strongly renormalizes the linear bands which occur near the Dirac crossing in graphene. The single bare Dirac crossing is split into two individual Dirac-like points, which are separated in energy but still at zero momentum relative to the K-point. A diamond-like structure occurs in between as a result of the formation of plasmarons. In this work we explore the combined effect of electron-electron and electron-phonon coupling on the renormalized energy dispersion, the spectral function and on the electronic density of states. We find that distinct signatures of the plasmaron structure are observable in the density of states with the split Dirac points presenting themselves as minima with quadratic dependence on energy about such points. By examining the slopes of both the density of states and the renormalized dispersion near the Fermi level, we illustrate how one can separate $k$-dependent and $ω$-dependent renormalizations and suggest how this might allow for the isolation of the renormalization due to the electron-phonon interaction from that of the electron-electron interaction.

preprint2011arXiv

Electronic structure in underdoped cuprates due to the emergence of a pseudogap

The phenomenological Green's function developed in the works of Yang, Rice and Zhang has been very successful in understanding many of the anomalous superconducting properties of the deeply underdoped cuprates. It is based on considerations of the resonating valence bond spin liquid approximation and is designed to describe the underdoped regime of the cuprates. Here we emphasize the region of doping, $x$, just below the quantum critical point at which the pseudogap develops. In addition to Luttinger hole pockets centered around the nodal direction, there are electron pockets near the antinodes which are connected to the hole pockets by gapped bridging contours. We determine the contours of nearest approach as would be measured in angular resolved photoemission experiments and emphasize signatures of the Fermi surface reconstruction from the large Fermi contour of Fermi liquid theory (which contains $1+x$ hole states) to the Luttinger pocket (which contains $x$ hole states). We find that the quasiparticle effective mass renormalization increases strongly towards the edge of the Luttinger pockets beyond which it diverges.

preprint2010arXiv

Effects of a particle-hole asymmetric pseudogap on Bogoliubov quasiparticles

We show that in the presence of a pseudogap, the spectral function in the superconducting state of the underdoped cuprates exhibits additional Bogoliubov quasiparticle peaks at both positive and negative energy which are revealed by the particle-hole asymmetry of the pseudogapped energy bands. This provides direct information on the unoccupied band via measurement of the occupied states. When sufficiently close, these Bogoliubov peaks will appear to merge with existing peaks leading to the anomalous observation, seen in experiment, that the carrier spectral density broadens with reduced temperature in the superconducting state. Using the resonating valence bond (RVB) spin liquid model in conjunction with recent angle-resolved photoemission spectroscopy (ARPES) data allows for an empirical determination of the temperature dependence of the pseudogap suggesting that it opens only very gradually below the pseudogap onset temperature $T^*$.

preprint2009arXiv

Effect of pseudogap formation on the penetration depth of underdoped high $T_c$ cuprates

The penetration depth is calculated over the entire doping range of the cuprate phase diagram with emphasis on the underdoped regime. Pseudogap formation on approaching the Mott transition, for doping below a quantum critical point, is described within a model based on the resonating valence bond spin liquid which provides an ansatz for the coherent piece of the Green's function. Fermi surface reconstruction, which is an essential element of the model, has a strong effect on the superfluid density at T=0 producing a sharp drop in magnitude, but does not change the slope of the linear low temperature variation. Comparison with recent data on Bi-based cuprates provides validation of the theory and shows that the effects of correlations, captured by Gutzwiller factors, are essential for a qualitative understanding of the data. We find that the Ferrell-Glover-Tinkham sum rule still holds and we compare our results with those for the Fermi arc and the nodal liquid models.

preprint2009arXiv

Signatures of Fermi surface reconstruction in Raman spectra of underdoped cuprates

We have calculated the Raman B$_{1g}$ and B$_{2g}$ spectra as a function of temperature, as well as doping, for the underdoped cuprates, using a model based on the resonating valence-bond spin-liquid. We discuss changes in intensity and peak position brought about by the presence of a pseudogap and the implied Fermi surface reconstruction, which are elements of this model. Signatures of Fermi surface reconstruction are evident as a sharp rise in the doping dependence of the antinodal to nodal peak ratio which occurs below the quantum critical point. The temperature dependence of the B$_{1g}$ polarization can be used to determine if the superconducting gap is limited to the Fermi pocket, as seen in angle resolved photoemission spectroscopy, or extends beyond. We find that the slope of the linear low energy B$_{2g}$ spectrum maintains its usual d-wave form, but with an effective gap which reflects the gap amplitude projected on the Fermi pocket. Our calculations capture the main qualitative features revealed in the extensive data set available on the HgBa$_2$CuO$_{4+δ}$ (Hg-1201) cuprate.