-
Possibilities for enhanced electron-phonon interactions and high-$T_c$ superconductivity in engineered bimetallic nano-structured superlattices
Authors:
Shinjan Mandal,
Shrihari Soundararajan,
Manish Jain,
H. R. Krishnamurthy
Abstract:
We explore theoretically the properties of engineered bimetallic nano-structured superlattices where an array of nano-clusters of a simple (single band) metal are embedded periodically inside another simple metal with a different work function. The exploration is done using a simplified tight-binding model with Coulomb interactions included, as well as density functional theory. Taking arrays of "…
▽ More
We explore theoretically the properties of engineered bimetallic nano-structured superlattices where an array of nano-clusters of a simple (single band) metal are embedded periodically inside another simple metal with a different work function. The exploration is done using a simplified tight-binding model with Coulomb interactions included, as well as density functional theory. Taking arrays of "Ag" clusters of fixed sizes and configurations (when unrelaxed) embedded periodically in an "Au" matrix as an example, we show that a significant enhancement of electron-phonon interactions ensues, implying possibilities for high-$T_c$ superconductivity. The enhancement stems from a strong coupling, via Coulomb interactions, between the dipolar charge distribution that forms at the Au-Ag interfaces and the breathing and other modes of vibration of the light Ag atoms caged inside the heavier Au matrix. The interface dipoles form because of the interplay between the mismatch of the local potential seen by the conduction electrons localised in Wannier orbitals at the Ag and Au sites (the Ag sites being slightly repulsive relative to the Au sites) and the (long-range) Coulomb repulsion between electrons occupying these Wannier orbitals. We also discuss the DC transport in such systems.
△ Less
Submitted 28 August, 2024;
originally announced August 2024.
-
Engineering ultra-strong electron-phonon coupling and nonclassical electron transport in crystalline gold with nanoscale interfaces
Authors:
Shreya Kumbhakar,
Tuhin Kumar Maji,
Binita Tongbram,
Shinjan Mandal,
Shri Hari Soundararaj,
Banashree Debnath,
T. Phanindra Sai,
Manish Jain,
H. R. Krishnamurthy,
Anshu Pandey,
Arindam Ghosh
Abstract:
Electrical resistivity in good metals, particularly noble metals such as gold (Au), silver (Ag), or copper, increases linearly with temperature ($T$) for $T > Θ_{\mathrm{D}}$, where $Θ_{\mathrm{D}}$ is the Debye temperature. This is because the coupling ($λ$) between the electrons and the lattice vibrations, or phonons, in these metals is rather weak with $λ\sim 0.1-0.2$, and a perturbative analys…
▽ More
Electrical resistivity in good metals, particularly noble metals such as gold (Au), silver (Ag), or copper, increases linearly with temperature ($T$) for $T > Θ_{\mathrm{D}}$, where $Θ_{\mathrm{D}}$ is the Debye temperature. This is because the coupling ($λ$) between the electrons and the lattice vibrations, or phonons, in these metals is rather weak with $λ\sim 0.1-0.2$, and a perturbative analysis suffices to explain the $T$-linear electron-phonon scattering rate. In this work, we outline a new nanostructuring strategy of crystalline Au where this foundational concept of metallic transport breaks down. We show that by embedding a distributed network of ultra-small Ag nanoparticles (AgNPs) of radius $\sim1-2$ nm inside a crystalline Au shell, an unprecedented enhancement in the electron-phonon interaction, with $λ$ as high as $\approx 20$, can be achieved. This is over hundred times that of bare Au or Ag, and ten times larger than any known metal. With increasing AgNP density, the electrical resistivity deviates from $T$-linearity, and approaches a saturation to the Mott-Ioffe-Regel scale $ρ_{\mathrm{MIR}}\sim h a /e^2$ for both disorder ($T\to 0$) and phonon ($T \gg Θ_{\mathrm{D}}$)-dependent components of resistivity (here, $a=0.3$~nm, is the lattice constant of Au). This giant electron-phonon interaction, which we suggest arises from the coulomb interaction-induced coupling of conduction electrons to the localized phonon modes at the buried Au-Ag hetero-interfaces, allows experimental access to a regime of nonclassical metallic transport that has never been probed before.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
Phonon Linewidths in Twisted Bilayer Graphene near Magic Angle
Authors:
Shinjan Mandal,
Indrajit Maity,
H. R. Krishnamurthy,
Manish Jain
Abstract:
We present a computational study of the phonon linewidths in twisted bilayer graphene arising from electron-phonon interactions and anharmonic effects. The electronic structure is calculated using distance-dependent transfer integrals based on the atomistic Slater-Koster tight-binding formalism, including electron-electron interactions treated at the Hartree level, and the phonons are calculated u…
▽ More
We present a computational study of the phonon linewidths in twisted bilayer graphene arising from electron-phonon interactions and anharmonic effects. The electronic structure is calculated using distance-dependent transfer integrals based on the atomistic Slater-Koster tight-binding formalism, including electron-electron interactions treated at the Hartree level, and the phonons are calculated using classical force fields. These ingredients are used to calculate the phonon linewidths arising from electron-phonon interactions. Furthermore, anharmonic effects on the linewidths are computed using the mode-projected velocity autocorrelation function obtained from classical molecular dynamics. We predict a moiré potential induced splitting of this mode, which arises due to contributions from high symmetry stacking regions. Our findings show that both electron-phonon and anharmonic effects have a significant impact on the linewidth of the Raman active G mode near the magic angle.
△ Less
Submitted 2 July, 2024; v1 submitted 15 February, 2024;
originally announced March 2024.
-
Electric field tunable superconductivity with competing orders in near magic-angle twisted bilayer graphene
Authors:
Ranit Dutta,
Ayan Ghosh,
Shinjan Mandal,
K. Watanabe,
T. Taniguchi,
H. R. Krishnamurthy,
Sumilan Banerjee,
Manish Jain,
Anindya Das
Abstract:
Superconductivity (SC) in twisted bilayer graphene (tBLG) has been explored by varying carrier concentrations, twist angles, and screening strength, with the aim of uncovering its origin and possible connections to strong electronic correlations in narrow bands and various resulting broken symmetries. However, the link between the tBLG band structure and the onset of SC and other orders largely re…
▽ More
Superconductivity (SC) in twisted bilayer graphene (tBLG) has been explored by varying carrier concentrations, twist angles, and screening strength, with the aim of uncovering its origin and possible connections to strong electronic correlations in narrow bands and various resulting broken symmetries. However, the link between the tBLG band structure and the onset of SC and other orders largely remains unclear. In this study, we address this crucial gap by examining in-situ band structure tuning of a near magic-angle ($θ\approx0.95^0$) tBLG device with displacement field ($D$) and reveal remarkable competition between SC and other broken symmetries. At zero $D$, the device exhibits superconducting signatures without the resistance peak at half-filling, a characteristic signature with a strong electronic correlation. As $D$ increases, the SC is suppressed, accompanied by the appearance of a resistance peak at half-filling. Hall density measurements reveal that at zero $D$, SC arises around the van Hove singularity (vHs) from an isospin or spin-valley unpolarized band. At higher $D$, the suppression of SC coincides with broken isospin symmetry near half-filling with lifted degeneracy ($g_d \sim 2$). Additionally, when SC is suppressed at higher $D$, density waves around the superconducting dome are seen. These findings, with our theoretical calculations, highlight the competition between SC and other orders, and the stabilization of the latter by the electric field enhanced nesting in the band structure of tBLG.
△ Less
Submitted 18 February, 2024;
originally announced February 2024.
-
Calculation of Gilbert damping and magnetic moment of inertia using torque-torque correlation model within ab initio Wannier framework
Authors:
Robin Bajaj,
Seung-Cheol Lee,
H. R. Krishnamurthy,
Satadeep Bhattacharjee,
Manish Jain
Abstract:
Magnetization dynamics in magnetic materials are well described by the modified semiclassical Landau-Lifshitz-Gilbert (LLG) equation, which includes the magnetic damping $α$ and the magnetic moment of inertia $\mathrm{I}$ tensors as key parameters. Both parameters are material-specific and physically represent the time scales of damping of precession and nutation in magnetization dynamics. $α$ and…
▽ More
Magnetization dynamics in magnetic materials are well described by the modified semiclassical Landau-Lifshitz-Gilbert (LLG) equation, which includes the magnetic damping $α$ and the magnetic moment of inertia $\mathrm{I}$ tensors as key parameters. Both parameters are material-specific and physically represent the time scales of damping of precession and nutation in magnetization dynamics. $α$ and $\mathrm{I}$ can be calculated quantum mechanically within the framework of the torque-torque correlation model. The quantities required for the calculation are torque matrix elements, the real and imaginary parts of the Green's function and its derivatives. Here, we calculate these parameters for the elemental magnets such as Fe, Co and Ni in an ab initio framework using density functional theory and Wannier functions. We also propose a method to calculate the torque matrix elements within the Wannier framework. We demonstrate the effectiveness of the method by comparing it with the experiments and the previous ab initio and empirical studies and show its potential to improve our understanding of spin dynamics and to facilitate the design of spintronic devices.
△ Less
Submitted 1 January, 2024;
originally announced January 2024.
-
Controlling Umklapp scattering in bilayer graphene moir'e superlattice
Authors:
Mohit Kumar Jat,
Shubhankar Mishra,
Harsimran Kaur Mann,
Robin Bajaj,
Kenji Watanabe,
Takashi Taniguchi,
H. R. Krishnamurthy,
Manish Jain,
Aveek Bid
Abstract:
In this Letter, we present experimental findings on electron-electron scattering in a two-dimensional moir'e heterostructure with tunable Fermi wave vector, reciprocal lattice vector, and band gap. We achieve this in high-mobility aligned heterostructures of bilayer graphene (BLG) and hBN. Around half-filling, the primary contribution to the resistance of BLG/hBN aligned superlattices arises from…
▽ More
In this Letter, we present experimental findings on electron-electron scattering in a two-dimensional moir'e heterostructure with tunable Fermi wave vector, reciprocal lattice vector, and band gap. We achieve this in high-mobility aligned heterostructures of bilayer graphene (BLG) and hBN. Around half-filling, the primary contribution to the resistance of BLG/hBN aligned superlattices arises from electron-electron Umklapp (Uee) scattering, making the resistance of graphene/hBN moir'e devices significantly larger than that of non-aligned devices (where Uee is forbidden). We quantify the strength of the Uee scattering and find that it follows a universal scaling with Fermi energy and has a non-monotonic dependence on the charge carrier density. The Uee scattering is strongly electric field tunable and affected by layer-polarization of BLG. It has a strong particle-hole asymmetry - the resistance when the chemical potential is in the conduction band is significantly lesser than when it is in the valence band, making the electron-doped regime more practical for potential applications.
△ Less
Submitted 15 February, 2024; v1 submitted 13 October, 2023;
originally announced October 2023.
-
Higher-order Bragg gaps in the electronic band structure of bilayer graphene renormalized by recursive supermoiré potential
Authors:
Mohit Kumar Jat,
Priya Tiwari,
Robin Bajaj,
Ishita Shitut,
Shinjan Mandal,
Kenji Watanabe,
Takashi Taniguchi,
H. R. Krishnamurthy,
Manish Jain,
Aveek Bid
Abstract:
This letter presents our findings on the recursive band gap engineering of chiral fermions in bilayer graphene doubly aligned with hBN. By utilizing two interfering moiré potentials, we generate a supermoiré pattern which renormalizes the electronic bands of the pristine bilayer graphene, resulting in higher-order fractal gaps even at very low energies. These Bragg gaps can be mapped using a uniqu…
▽ More
This letter presents our findings on the recursive band gap engineering of chiral fermions in bilayer graphene doubly aligned with hBN. By utilizing two interfering moiré potentials, we generate a supermoiré pattern which renormalizes the electronic bands of the pristine bilayer graphene, resulting in higher-order fractal gaps even at very low energies. These Bragg gaps can be mapped using a unique linear combination of periodic areas within the system. To validate our findings, we used electronic transport measurements to identify the position of these gaps as functions of the carrier density and establish their agreement with the predicted carrier densities and corresponding quantum numbers obtained using the continuum model. Our work provides direct experimental evidence of the quantization of the area of quasi-Brillouin zones in supermoiré systems. It fills essential gaps in understanding the band structure engineering of Dirac fermions by a recursive doubly periodic superlattice potential.
△ Less
Submitted 4 April, 2023;
originally announced April 2023.
-
Exciton fine structure in twisted transition metal dichalcogenide heterostructures
Authors:
Sudipta Kundu,
Tomer Amit,
H. R. Krishnamurthy,
Manish Jain,
Sivan Refaely-Abramson
Abstract:
Moiré superlattices of transition metal dichalcogenide (TMD) heterostructures give rise to rich excitonic phenomena associated with the interlayer twist angle and induced changes in the involved quantum states. Theoretical calculations of excitons in such systems are typically based on model moiré potentials to mitigate the computational cost. However, an ab initio understanding of the electron-ho…
▽ More
Moiré superlattices of transition metal dichalcogenide (TMD) heterostructures give rise to rich excitonic phenomena associated with the interlayer twist angle and induced changes in the involved quantum states. Theoretical calculations of excitons in such systems are typically based on model moiré potentials to mitigate the computational cost. However, an ab initio understanding of the electron-hole coupling dominating the excitations is crucial to realize the twist-induced modifications of the optical selection rules. In this work we use many-body perturbation theory to compute and analyze the relation between twist angle and exciton properties in twisted TMD heterostructures. We present a general approach for unfolding excitonic states from the moiré Brillouin zone onto the Brillouin zones of the separate layers. Applying this method to a twisted MoS$_2$/MoSe$_2$ bilayer, we find that the optical excitation spectrum is dominated by mixed transitions between electrons and holes with different momenta in the separate monolayers, leading to unexpected and angle-dependent hybridization between interlayer and intralayer excitons. Our findings offer a design pathway for tuning exciton layer-localization in TMD heterostructures as a function of twist angle.
△ Less
Submitted 12 September, 2022;
originally announced September 2022.
-
Atomic reconstruction and flat bands in strain engineered transition metal dichalcogenide bilayer moiré systems
Authors:
Sudipta Kundu,
Indrajit Maity,
Robin Bajaj,
H. R. Krishnamurthy,
Manish Jain
Abstract:
Strain-induced lattice mismatch leads to moiré patterns in homobilayer transition metal dichalcogenides (TMDs). We investigate the structural and electronic properties of such strained moiré patterns in TMD homobilayers. The moiré patterns in strained TMDs consist of several stacking domains which are separated by tensile solitons. Relaxation of these systems distributes the strain unevenly in the…
▽ More
Strain-induced lattice mismatch leads to moiré patterns in homobilayer transition metal dichalcogenides (TMDs). We investigate the structural and electronic properties of such strained moiré patterns in TMD homobilayers. The moiré patterns in strained TMDs consist of several stacking domains which are separated by tensile solitons. Relaxation of these systems distributes the strain unevenly in the moiré superlattice, with the maximum strain energy concentrating at the highest energy stackings. The order parameter distribution shows the formation of aster topological defects at the same sites. In contrast, twisted TMDs host shear solitons at the domain walls, and the order parameter distribution in these systems shows the formation of vortex defects. The strained moiré systems also show the emergence of several well-separated flat bands at both the valence and conduction band edges, and we observe a significant reduction in the band gap. The flat bands in these strained moiré superlattices provide platforms for studying the Hubbard model on a triangular lattice as well as the ionic Hubbard model on a honeycomb lattice. Furthermore, we study the localization of the wave functions corresponding to these flat bands. The wave functions localize at different stackings compared to twisted TMDs, and our results are in excellent agreement with recent spectroscopic experiments [1].
△ Less
Submitted 2 August, 2021;
originally announced August 2021.
-
Moiré induced topology and flat bands in twisted bilayer WSe$_2$: A first-principles study
Authors:
Sudipta Kundu,
Mit H. Naik,
H. R. Krishnamurthy,
Manish Jain
Abstract:
We study the influence of strong spin-orbit interaction on the formation of flat bands in relaxed twisted bilayer WSe$_2$. Flat bands, well separated in energy, emerge at the band edges for twist angles ($θ$) near 0$^{\circ}$ and 60$^{\circ}$. For $θ$ near 0$^{\circ}$, the interlayer hybridization together with a moiré potential determines the electronic structure. The bands near the valence band…
▽ More
We study the influence of strong spin-orbit interaction on the formation of flat bands in relaxed twisted bilayer WSe$_2$. Flat bands, well separated in energy, emerge at the band edges for twist angles ($θ$) near 0$^{\circ}$ and 60$^{\circ}$. For $θ$ near 0$^{\circ}$, the interlayer hybridization together with a moiré potential determines the electronic structure. The bands near the valence band edge have nontrivial topology, with Chern numbers equal to +1 or $-$1. We propose that the nontrivial topology of the first band can be probed experimentally for twist angles less than a critical angle of 3.5$^{\circ}$. For $θ$ near 60$^{\circ}$, the flattening of the bands arising from the K point of the unit cell Brillouin zone is a result of atomic rearrangements in the individual layers. Our findings on the flat bands and the localization of their wavefunctions for both ranges of $θ$ match well with recent experimental observations.
△ Less
Submitted 22 February, 2022; v1 submitted 12 March, 2021;
originally announced March 2021.
-
Unconventional superconductivity in a strongly correlated band-insulator without doping
Authors:
Anwesha Chattopadhyay,
H. R. Krishnamurthy,
Arti Garg
Abstract:
We present a novel route for attaining unconventional superconductivity (SC) in a strongly correlated system without doping. In a simple model of a correlated band insulator (BI) at half-filling we demonstrate, based on a generalization of the projected wavefunctions method, that SC emerges when e-e interactions and the bare band-gap are both much larger than the kinetic energy, provided the syste…
▽ More
We present a novel route for attaining unconventional superconductivity (SC) in a strongly correlated system without doping. In a simple model of a correlated band insulator (BI) at half-filling we demonstrate, based on a generalization of the projected wavefunctions method, that SC emerges when e-e interactions and the bare band-gap are both much larger than the kinetic energy, provided the system has sufficient frustration against the magnetic order. As the interactions are tuned, SC appears sandwiched between the correlated BI followed by a paramagnetic metal on one side, and a ferrimagnetic metal, antiferromagnetic (AF) half-metal, and AF Mott insulator phases on the other side.
△ Less
Submitted 9 July, 2020; v1 submitted 7 July, 2020;
originally announced July 2020.
-
Excess entropy and breakdown of semiclassical description of thermoelectricity in twisted bilayer graphene close to half filling
Authors:
Bhaskar Ghawri,
Phanibhusan S. Mahapatra,
Shinjan Mandal,
Aditya Jayaraman,
Manjari Garg,
Kenji Watanabe,
Takashi Taniguchi,
H. R. Krishnamurthy,
Manish Jain,
Sumilan Banerjee,
U. Chandni,
Arindam Ghosh
Abstract:
In moiré systems with twisted bilayer graphene (tBLG), the amplification of Coulomb correlation effects at low twist angles ($θ$) is a result of nearly flat low-energy electronic bands and divergent density of states (DOS) at van Hove singularities (vHS). This not only causes superconductivity, Mott insulating states, and quantum anomalous Hall effect close to the critical (or magic) angle…
▽ More
In moiré systems with twisted bilayer graphene (tBLG), the amplification of Coulomb correlation effects at low twist angles ($θ$) is a result of nearly flat low-energy electronic bands and divergent density of states (DOS) at van Hove singularities (vHS). This not only causes superconductivity, Mott insulating states, and quantum anomalous Hall effect close to the critical (or magic) angle $θ= θ_{c} \approx 1.1^\circ$, but also unconventional metallic states that are claimed to exhibit non-Fermi liquid (NFL) excitations. However, unlike superconductivity and the correlation-induced gap in the DOS, unambiguous signatures of NFL effects in the metallic state remain experimentally elusive. Here we report simultaneous measurement of electrical resistivity ($ρ$) and thermoelectric power ($S$) in tBLG at $θ\approx 1.6^\circ$. We observe an emergent violation of the semiclassical Mott relation in the form of excess $S$ close to half-filling. The excess $S$ ($\approx 2$ $μ$V/K at low temperature $T \sim 10$ K) persists up to $\approx 40$ K, and is accompanied by metallic $T$-linear $ρ$ with transport scattering rate ($τ^{-1}$) of near-Planckian magnitude $τ^{-1} \sim k_{B}T/\hbar$. The combination of non-trivial electrical transport and violation of Mott relation provides compelling evidence of NFL physics intrinsic to tBLG, at small twist angle and half-filling.
△ Less
Submitted 26 April, 2020;
originally announced April 2020.
-
Reconstruction of moiré lattices in twisted transition metal dichalcogenide bilayers
Authors:
Indrajit Maity,
Prabal K. Maiti,
H. R. Krishnamurthy,
Manish Jain
Abstract:
An important step in understanding the exotic electronic, vibrational, and optical properties of the moiré lattices is the inclusion of the effects of structural relaxation of the un-relaxed moiré lattices. Here, we propose novel structures for twisted bilayer of transition metal dichalcogenides (TMDs). For $θ\gtrsim 58.4^{\circ}$, we show a dramatic reconstruction of the moiré lattices, leading t…
▽ More
An important step in understanding the exotic electronic, vibrational, and optical properties of the moiré lattices is the inclusion of the effects of structural relaxation of the un-relaxed moiré lattices. Here, we propose novel structures for twisted bilayer of transition metal dichalcogenides (TMDs). For $θ\gtrsim 58.4^{\circ}$, we show a dramatic reconstruction of the moiré lattices, leading to a trimerization of the unfavorable stackings. We show that the development of curved domain walls due to the three-fold symmetry of the stacking energy landscape is responsible for such lattice reconstruction. Furthermore, we show that the lattice reconstruction notably changes the electronic band-structure. This includes the occurrence of flat bands near the edges of the conduction as well as valence bands, with the valence band maximum, in particular, corresponding to localized states enclosed by the trimer. We also find possibilities for other complicated, entropy stabilized, lattice reconstructed structures.
△ Less
Submitted 18 February, 2021; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Correlation driven metallic and half-metallic phases in a band insulator
Authors:
Soumen Bag,
Arti Garg,
H. R. Krishnamurthy
Abstract:
We demonstrate, using dynamical mean-field theory with the hybridization expansion continuous time quantum montecarlo impurity solver, a rich phase diagram with {\em correlation driven metallic and half-metallic phases} in a simple model of a correlated band insulator, namely, the half-filled ionic Hubbard model (IHM) with first {\em and} second neighbor hopping ($t$ and $t'$), an on-site repulsio…
▽ More
We demonstrate, using dynamical mean-field theory with the hybridization expansion continuous time quantum montecarlo impurity solver, a rich phase diagram with {\em correlation driven metallic and half-metallic phases} in a simple model of a correlated band insulator, namely, the half-filled ionic Hubbard model (IHM) with first {\em and} second neighbor hopping ($t$ and $t'$), an on-site repulsion $U$, and a staggered potential $Δ$. Without $t'$ the IHM has a direct transition from a paramagnetic band insulator (BI) to an antiferromagnetic Mott insulator (AFI) phase as $U$ increases. For weak to intermediate correlations, $t'$ frustrates the AF order, leading to a paramagnetic metal (PM) phase, a ferrimagnetic metal (FM) phase and an anti-ferromagnetic half-metal (AFHM) phase in which electrons with one spin orientation, say up-spin, have gapless excitations while the down-spin electrons are gapped. For $t'$ less than a threshold $ t_1$, there is a direct, first-order, BI to AFI transition as $U$ increases, as for $t'=0$; for $t_4< t' < Δ/2$, the BI to AFI transition occurs via an intervening PM phase. For $t' > Δ/2$, there is no BI phase, and the system has a PM to AFI transition as $U$ increases. In an intermediate-range $t_2 < t' < t_3$, as $U$ increases the system undergoes four transitions, in the sequence BI $\rightarrow$ PM $\rightarrow$ FM $\rightarrow$ AFHM $\rightarrow$ AFI; the FM phase is absent in the ranges of $t'$ on either side, implying three transitions. The BI-PM, FM-AFHM and AFHM-AFI transitions, and a part of the PM-FM transition are continuous, while the rest of the transitions are first order in nature. The PM, FM and the AFHM phases have, respectively, spin symmetric, partially polarized and fully polarized electron [hole] pockets around the ($\pmπ/2$, $\pmπ/2$) [($\pm π, 0$), ($0. \pm π$)] points in the Brillouin zone.
△ Less
Submitted 9 September, 2019;
originally announced September 2019.
-
Phonons in Twisted Transition Metal Dichalcogenide Bilayers ("Twistnonics"): Ultra-soft Phasons, and a transition from Superlubric to Pinned Phase
Authors:
Indrajit Maity,
Mit H. Naik,
Prabal K Maiti,
H. R. Krishnamurthy,
Manish Jain
Abstract:
The tunability of the interlayer coupling by twisting one layer with respect to another layer of two-dimensional materials provides a unique way to manipulate the phonons and related properties. We refer to this engineering of phononic properties as "Twistnonics". We study the effects of twisting on low-frequency shear (SM) and layer breathing (LBM) modes in transition metal dichalcogenide (TMD) b…
▽ More
The tunability of the interlayer coupling by twisting one layer with respect to another layer of two-dimensional materials provides a unique way to manipulate the phonons and related properties. We refer to this engineering of phononic properties as "Twistnonics". We study the effects of twisting on low-frequency shear (SM) and layer breathing (LBM) modes in transition metal dichalcogenide (TMD) bilayer using atomistic classical simulations. We show that these low-frequency modes are extremely sensitive to twist and can be used to infer the twist angle. We find unique "ultra-soft" phason modes (frequency $\lesssim 1\ \mathrm{cm^{-1}}$, comparable to acoustic modes) for any non-zero twist, corresponding to an \textit{effective} translation of the moir{é} lattice by relative displacement of the constituent layers in a non-trivial way. Unlike the acoustic modes, the velocity of the phason modes is quite sensitive to twist angle. As twist angle decreases, ($θ\lesssim 3^{\circ},\ \gtrsim 57^{\circ}$) the ultra-soft modes represent the acoustic modes of the "emergent" soft moir{é} scale lattice. Also, new high-frequency SMs appear, identical to those in stable bilayer TMD ($θ= 0\degree/60\degree$), due to the overwhelming growth of stable stacking regions in relaxed twisted structures. Furthermore, we find remarkably different structural relaxation as $θ\to 0^{\circ}$, $\to 60^{\circ}$ due to sub-lattice symmetry breaking. Our study reveals the possibility of an intriguing $θ$ dependent superlubric to pinning behavior and of the existence of ultra-soft modes in \textit{all} two-dimensional (2D) materials.
△ Less
Submitted 20 December, 2019; v1 submitted 27 May, 2019;
originally announced May 2019.
-
Is there a superconducting phase in the half-filled ionic Hubbard model ?
Authors:
Anwesha Chattopadhyay,
Soumen Bag,
H. R. Krishnamurthy,
Arti Garg
Abstract:
We investigate the ionic Hubbard model (IHM) at half-filling in the limit of strong correlations and large ionic potential. The low energy effective Hamiltonian in this limit, obtained by a similarity transformation, is a modified $t-J$ model with effective second neighbour hopping terms. We explore the possibilities of d-wave pairing and extended s-wave pairing superconducting (SC) phases on a tw…
▽ More
We investigate the ionic Hubbard model (IHM) at half-filling in the limit of strong correlations and large ionic potential. The low energy effective Hamiltonian in this limit, obtained by a similarity transformation, is a modified $t-J$ model with effective second neighbour hopping terms. We explore the possibilities of d-wave pairing and extended s-wave pairing superconducting (SC) phases on a two dimensional square lattice at zero temperature within a Gutzwiller projected renormalized mean field theory. In the sector of solutions that forbid spin ordering, the system shows a finite non-zero d-wave as well as extended s-wave pairing amplitude for $Δ\sim U \gg t$. The width of the superconducting phase in $U-Δ$ regime shrinks with increase in $U$ and $Δ$, though the extended s-wave pairing phase is higher in energy than the d-wave pairing superconducting phase. But in a spin resolved renormalized mean field calculation, which allows for an antiferromagnetic (AF) order along with the d-wave or extended s-wave pairing, the SC phase is no longer viable and the system shows a direct transition from an AF ordered phase to a paramagnetic band insulator. Except for a thin sliver of a half-metallic AF phase close to the AF transition point, most of the AF ordered phase is a Mott insulator. We benchmarked the AF Mott insulator to band insulator transition within the Gutzwiller projected renormalized mean field theory against the dynamical mean field theory (DMFT) solved using continuous time quantum Monte-Carlo (CTQMC). Our work suggests that the ground state phase diagram of the IHM at half-filling in the limit of extreme correlations does not have any SC phase. The SC phase seen in the paramagnetic sector is a metastable phase, being higher in energy than the AF Mott insulator phase.
△ Less
Submitted 15 February, 2019;
originally announced February 2019.
-
Reversal of particle-hole scattering-rate asymmetry in Anderson impurity model
Authors:
R. Zitko,
H. R. Krishnamurthy,
B. Sriram Shastry
Abstract:
We study the particle-hole asymmetry of the scattering rate in strongly correlated electron systems by examining the cubic $ω^3$ and $ωT^2$ terms in the imaginary part of the self-energy of the Anderson impurity model. We show that the sign is opposite in the weak-coupling and strong-coupling limits, explaining the differences found in theoretical approaches taking the respective limits as the sta…
▽ More
We study the particle-hole asymmetry of the scattering rate in strongly correlated electron systems by examining the cubic $ω^3$ and $ωT^2$ terms in the imaginary part of the self-energy of the Anderson impurity model. We show that the sign is opposite in the weak-coupling and strong-coupling limits, explaining the differences found in theoretical approaches taking the respective limits as the starting points. The sign change in fact precisely delineates the cross-over between the weak and strong correlation regimes of the model. For weak interaction $U$ the sign reversal occurs for small values of the doping $δ=1-n$, while for interaction of order $U \approx 2 Γ$, $Γ$ being the hybridization strength, the cross-over curve rapidly shifts to the large-doping range. This curve based on the impurity dynamics is genuinely different from other cross-over curves defined through impurity thermodynamic and static properties.
△ Less
Submitted 10 October, 2018; v1 submitted 30 July, 2018;
originally announced July 2018.
-
Surprises in the t-J model: Implications for cuprates
Authors:
Aabhaas Vineet Mallik,
Gaurav Kumar Gupta,
Vijay B. Shenoy,
H. R. Krishnamurthy
Abstract:
The t-J model is a paradigmatic model for the study of strongly correlated electron systems. In particular, it has been argued that it is an appropriate model to describe the cuprate high-Tc superconductors. It turns out that a comprehensive understanding of the gamut of physics encoded by the t-J model is still an open problem. In recent years some remarkable experiments on the cuprates, for exam…
▽ More
The t-J model is a paradigmatic model for the study of strongly correlated electron systems. In particular, it has been argued that it is an appropriate model to describe the cuprate high-Tc superconductors. It turns out that a comprehensive understanding of the gamut of physics encoded by the t-J model is still an open problem. In recent years some remarkable experiments on the cuprates, for example, discovery of nodeless superconductivity in underdoped samples (PNAS 109, 18332 (2012)), discovery of s-wave like gap in the pseudogap phase (Phys. Rev. Lett. 111, 107001 (2013)), and observation of polar Kerr effect (PKE) (Phys. Rev. Lett. 112, 047003 (2014)), have thrown up new challenges for this model. Here, we present results demonstrating that, within the slave-particle formulation of the t-J model, the d-wave superconductor is unstable at low doping to its own anti-symmetric phase mode fluctuations when the effect of fluctuations is treated self-consistently. We then show that this instability gives way to a time reversal symmetry broken d + is-SC in the underdoped region which has superfluid stiffness consistent with Uemura relation, even with a large pair amplitude. We show that our results are consistent with existing experiments on cuprates and suggest that Josephson (SQUID interferometry) experiments can clearly distinguish the d+is-SC from a host of other possibilities alluded to be contributing to the physics of underdoped cuprates. We also comment on other theoretical studies vis-a-vis ours.
△ Less
Submitted 7 May, 2018;
originally announced May 2018.
-
Seebeck Coefficient of a Single van der Waals Junction in Twisted Bilayer Graphene
Authors:
Phanibhusan S. Mahapatra,
Kingshuk Sarkar,
H. R. Krishnamurthy,
Subroto Mukerjee,
Arindam Ghosh
Abstract:
When two planar atomic membranes are placed within the van der Waals distance, the charge and heat transport across the interface are coupled by the rules of momentum conservation and structural commensurability, leading to outstanding thermoelectric properties. Here we show that an effective "interlayer phonon drag" determines the Seebeck coefficient (S) across the van der Waals gap formed in twi…
▽ More
When two planar atomic membranes are placed within the van der Waals distance, the charge and heat transport across the interface are coupled by the rules of momentum conservation and structural commensurability, leading to outstanding thermoelectric properties. Here we show that an effective "interlayer phonon drag" determines the Seebeck coefficient (S) across the van der Waals gap formed in twisted bilayer graphene (tBLG). The cross-plane thermovoltage, which is non-monotonic in both temperature and density, is generated through scattering of electrons by the out-of-plane layer breathing (ZO'/ZA2) phonon modes and differs dramatically from the expected Landauer-Buttiker formalism in conventional tunnel junctions. The tunability of the cross-plane Seebeck effect in van der Waals junctions may be valuable in creating a new genre of versatile thermoelectric systems with layered solids.
△ Less
Submitted 4 January, 2018;
originally announced January 2018.
-
Ken Wilson - A Tribute: Some recollections and a few thoughts on education
Authors:
H. R. Krishnamurthy
Abstract:
I had the marvelous good fortune to be Ken Wilson's graduate student at the Physics Department, Cornell University, from 1972 to 1976. In this article, I present some recollections of how this came about, my interactions with Ken, and Cornell during this period; and acknowledge my debt to Ken, and to John Wilkins and Michael Fisher, who I was privileged to have as my main mentors at Cornell. I end…
▽ More
I had the marvelous good fortune to be Ken Wilson's graduate student at the Physics Department, Cornell University, from 1972 to 1976. In this article, I present some recollections of how this came about, my interactions with Ken, and Cornell during this period; and acknowledge my debt to Ken, and to John Wilkins and Michael Fisher, who I was privileged to have as my main mentors at Cornell. I end with some thoughts on the challenges of reforming education, a subject that was one of Ken's major preoccupations in the second half of his professional life.
△ Less
Submitted 31 December, 2016;
originally announced January 2017.
-
Constant Matrix Element Approximation to Time-Resolved Angle-Resolved Photoemission Spectroscopy
Authors:
J. K. Freericks,
H. R. Krishnamurthy
Abstract:
We describe different approximations associated with employing a constant matrix element for the coupling of light to multiband electrons in the context of time-resolved angle-resolved photoemission spectroscopy (TR-ARPES). In particular, we demonstrate that the constant matrix approximation only holds for one of the coordinate representations, and changing to other bases, requires including nonco…
▽ More
We describe different approximations associated with employing a constant matrix element for the coupling of light to multiband electrons in the context of time-resolved angle-resolved photoemission spectroscopy (TR-ARPES). In particular, we demonstrate that the constant matrix approximation only holds for one of the coordinate representations, and changing to other bases, requires including nonconstant corrections to the matrix element. We also discuss some simplifying approximations, where a constant matrix element is employed in multiple bases, and the consequences of this further approximation (especially with respect to the TR-ARPES signal no longer being nonnegative). We also discuss issues related to gauge invariance of the final spectra.
△ Less
Submitted 8 October, 2016;
originally announced October 2016.
-
Generalized gradient expansion for inhomogeneous dynamical mean-field theory: Application to ultracold atoms in a harmonic trap
Authors:
J. K. Freericks,
Shuyang Han,
Karlis Mikelsons,
H. R. Krishnamurthy
Abstract:
We develop a generalized gradient expansion of the inhomogeneous dynamical mean-field theory method for determining properties of ultracold atoms in a trap. This approach goes beyond the well-known local density approximation and at higher temperatures, in the normal phase, it shows why the local density approximation works so well, since the local density and generalized gradient approximations a…
▽ More
We develop a generalized gradient expansion of the inhomogeneous dynamical mean-field theory method for determining properties of ultracold atoms in a trap. This approach goes beyond the well-known local density approximation and at higher temperatures, in the normal phase, it shows why the local density approximation works so well, since the local density and generalized gradient approximations are essentially indistinguishable from each other (and from the exact solution within full inhomogeneous dynamical mean-field theory). But because the generalized gradient expansion only involves nearest-neighbor corrections, it does not work as well at low temperatures, when the systems enter into ordered phases. This is primarily due to the problem that ordered phases often satisfy some global constraints which determine the spatial ordering pattern, and the local density and generalized gradient approximations are not able to impose those kinds of constraints; they also overestimate the tendency to order. The theory is applied to phase separation of different mass fermionic mixtures represented by the Falicov-Kimball model and to determining the entropy per particle of a fermionic system represented by the Hubbard model. The generalized gradient approximation is a useful diagnostic for the accuracy of the local density approximation---when both methods agree, they are likely accurate, when they disagree, neither is likely to be correct.
△ Less
Submitted 1 June, 2016;
originally announced June 2016.
-
The role of average time dependence on the relaxation of excited electron populations in nonequilibrium many-body physics
Authors:
A. F. Kemper,
H. R. Krishnamurthy,
J. K. Freericks
Abstract:
We examine the exact equation of motion for the relaxation of populations of strongly correlated electrons after a nonequilibrium excitation by a pulsed field, and prove that the populations do not change when the Green's functions have no average time dependence. We show how the average time dependence enters into the equation of motion to lowest order and describe what governs the relaxation pro…
▽ More
We examine the exact equation of motion for the relaxation of populations of strongly correlated electrons after a nonequilibrium excitation by a pulsed field, and prove that the populations do not change when the Green's functions have no average time dependence. We show how the average time dependence enters into the equation of motion to lowest order and describe what governs the relaxation process of the electron populations in the long-time limit. While this result may appear, on the surface, to be required by any steady-state solution, the proof is nontrivial, and provides new critical insight into how nonequilibrium populations relax, which goes beyond the assumption that they thermalize via a simple relaxation rate determined by the imaginary part of the self-energy, or that they can be described by a quasi-equilibrium condition with a Fermi-Dirac distribution and a time-dependent temperature. We also discuss the implications of this result to approximate theories, which may not satisfy the exact relation in the equation of motion.
△ Less
Submitted 25 October, 2016; v1 submitted 30 April, 2016;
originally announced May 2016.
-
Crucial role of Internal Collective Modes in Underdoped Cuprates
Authors:
Aabhaas V. Mallik,
Umesh K. Yadav,
Amal Medhi,
H. R. Krishnamurthy,
Vijay B. Shenoy
Abstract:
The enigmatic cuprate superconductors have attracted resurgent interest with several recent reports and discussions of competing orders in the underdoped side. Motivated by this, here we address the natural question of fragility of the d-wave superconducting state in underdoped cuprates. Using a combination of theoretical approaches we study t-J like models, and discover an - as yet unexplored - i…
▽ More
The enigmatic cuprate superconductors have attracted resurgent interest with several recent reports and discussions of competing orders in the underdoped side. Motivated by this, here we address the natural question of fragility of the d-wave superconducting state in underdoped cuprates. Using a combination of theoretical approaches we study t-J like models, and discover an - as yet unexplored - instability that is brought about by an "internal" (anti-symmetric mode) fluctuation of the d-wave state. This new theoretical result is in good agreement with recent STM and ARPES studies of cuprates. We also suggest experimental directions to uncover this physics.
△ Less
Submitted 7 February, 2017; v1 submitted 31 March, 2016;
originally announced March 2016.
-
Renormalization Group Study of a Fragile Fermi liquid in $1+ε$ dimensions
Authors:
Peizhi Mai,
H. R. Krishnamurthy,
B. Sriram Shastry
Abstract:
We present a calculation of the low energy Greens function in $1+ε$ dimensions using the method of extended poor man's scaling, developed here. We compute the wave function renormalization $Z(ω)$ and also the decay rate near the Fermi energy. Despite the lack of $ω^2$ damping characteristic of 3-dimensional Fermi liquids, we show that quasiparticles do exist in $1+ε$ dimensions, in the sense that…
▽ More
We present a calculation of the low energy Greens function in $1+ε$ dimensions using the method of extended poor man's scaling, developed here. We compute the wave function renormalization $Z(ω)$ and also the decay rate near the Fermi energy. Despite the lack of $ω^2$ damping characteristic of 3-dimensional Fermi liquids, we show that quasiparticles do exist in $1+ε$ dimensions, in the sense that the quasiparticle weight $Z$ is finite and that the damping rate is smaller than the energy. We explicitly compute the crossover from this behavior to a 1-dimensional type Tomonaga-Luttinger liquid behavior at higher energies.
△ Less
Submitted 14 February, 2016; v1 submitted 11 February, 2016;
originally announced February 2016.
-
Phase Diagram of the Half-Filled Ionic Hubbard Model
Authors:
Soumen Bag,
Arti Garg,
H. R. Krishnamurthy
Abstract:
We study the phase diagram of the ionic Hubbard model (IHM) at half-filling using dynamical mean field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered potential $Δ$ and the on-site Hubbard U. In both the methods we find that for a finite…
▽ More
We study the phase diagram of the ionic Hubbard model (IHM) at half-filling using dynamical mean field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered potential $Δ$ and the on-site Hubbard U. In both the methods we find that for a finite $Δ$ and at zero temperature, anti-ferromagnetic (AFM) order sets in beyond a threshold $U=U_{AF}$ via a first order phase transition below which the system is a paramagnetic band insulator. Both the methods show a clear evidence for a transition to a half-metal phase just after the AFM order is turned on, followed by the formation of an AFM insulator on further increasing U. We show that the results obtained within both the methods have good qualitative and quantitative consistency in the intermediate to strong coupling regime. On increasing the temperature, the AFM order is lost via a first order phase transition at a transition temperature $T_{AF}(U, Δ)$ within both the methods, for weak to intermediate values of U/t. But in the strongly correlated regime, where the effective low energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) transition from the AFM phase to the paramagnetic phase, but the CTQMC does. As a result, at any finite temperature T, DMFT+CTQMC shows a second phase transition (not seen within DMFT+IPT) on increasing U beyond $U_{AF}$. At $U_N > U_{AF}$, when the Neel temperature $T_N$ for the effective Heisenberg model becomes lower than T, the AFM order is lost via a second order transition. In the 3-dimensonal parameter space of $(U/t,T/t,Δ/t)$, there is a line of tricritical points that separates the surfaces of first and second order phase transitions.
△ Less
Submitted 18 May, 2015; v1 submitted 30 March, 2015;
originally announced March 2015.
-
Feshbach modulation spectroscopy
Authors:
A. Dirks,
K. Mikelsons,
H. R. Krishnamurthy,
J. K. Freericks
Abstract:
In the vicinity of a Feshbach resonance, a system of ultracold atoms on an optical lattice undergoes rich physical transformations which involve molecule formation and hopping of molecules on the lattice and thus goes beyond a single-band Hubbard model description. We propose to probe the behavior of this system with a harmonic modulation of the magnetic field, and thus of the scattering length, a…
▽ More
In the vicinity of a Feshbach resonance, a system of ultracold atoms on an optical lattice undergoes rich physical transformations which involve molecule formation and hopping of molecules on the lattice and thus goes beyond a single-band Hubbard model description. We propose to probe the behavior of this system with a harmonic modulation of the magnetic field, and thus of the scattering length, across the Feshbach resonance, as an alternative to lattice-depth modulation spectroscopy. In the regime in which the single-band Hubbard model is still valid, we provide simulation data for this type of spectroscopy. The method may uncover a route towards the efficient creation of ultracold molecules and provides an alternate means for lattice modulation spectroscopy.
△ Less
Submitted 1 August, 2014;
originally announced August 2014.
-
Gauge invariance in the theoretical description of time-resolved angle-resolved pump/probe photoemission spectroscopy
Authors:
J. K. Freericks,
H. R. Krishnamurthy,
M. A. Sentef,
T. P. Devereaux
Abstract:
Nonequilibrium calculations in the presence of an electric field are usually performed in a gauge, and need to be transformed to reveal the gauge-invariant observables. In this work, we discuss the issue of gauge invariance in the context of time-resolved angle-resolved pump/probe photoemission. If the probe is applied while the pump is still on, one must ensure that the calculations of the observ…
▽ More
Nonequilibrium calculations in the presence of an electric field are usually performed in a gauge, and need to be transformed to reveal the gauge-invariant observables. In this work, we discuss the issue of gauge invariance in the context of time-resolved angle-resolved pump/probe photoemission. If the probe is applied while the pump is still on, one must ensure that the calculations of the observed photocurrent are gauge invariant. We also discuss the requirement of the photoemission signal to be positive and the relationship of this constraint to gauge invariance. We end by discussing some technical details related to the perturbative derivation of the photoemission spectra, which involve processes where the pump pulse photoexcites electrons due to nonequilibrium effects.
△ Less
Submitted 28 March, 2014;
originally announced March 2014.
-
Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory
Authors:
A. Dirks,
K. Mikelsons,
H. R. Krishnamurthy,
J. K. Freericks
Abstract:
We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and…
▽ More
We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse $10^d\times 10^d$ matrix, with $d > 6$. As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in non-equilibrium setups without mass transport.
△ Less
Submitted 23 September, 2013;
originally announced September 2013.
-
Nonequilibrium "melting" of a charge density wave insulator via an ultrafast laser pulse
Authors:
W. Shen,
Yizhi Ge,
A. Y. Liu,
H. R. Krishnamurthy,
T. P. Devereaux,
J. K. Freericks
Abstract:
We employ an exact solution of the simplest model for pump-probe time-resolved photoemission spectroscopy in charge-density-wave systems to show how, in nonequilibrium the gap in the density of states disappears while the charge density remains modulated, and then the gap reforms after the pulse has passed. This nonequilibrium scenario qualitatively describes the common short-time experimental fea…
▽ More
We employ an exact solution of the simplest model for pump-probe time-resolved photoemission spectroscopy in charge-density-wave systems to show how, in nonequilibrium the gap in the density of states disappears while the charge density remains modulated, and then the gap reforms after the pulse has passed. This nonequilibrium scenario qualitatively describes the common short-time experimental features in TaS2 and TbTe3 indicating a quasiuniversality for nonequilibrium "melting" with qualitative features that can be easily understood within a simple picture.
△ Less
Submitted 28 August, 2013;
originally announced August 2013.
-
Strong-coupling expansion for ultracold bosons in an optical lattice at finite temperatures in the presence of superfluidity
Authors:
Manjari Gupta,
H. R. Krishnamurthy,
J. K. Freericks
Abstract:
We develop a strong-coupling ($t \ll U$) expansion technique for calculating the density profile for bosonic atoms trapped in an optical lattice with an overall harmonic trap at finite temperature and finite on site interaction in the presence of superfluid regions. Our results match well with quantum Monte Carlo simulations at finite temperature. We also show that the superfluid order parameter n…
▽ More
We develop a strong-coupling ($t \ll U$) expansion technique for calculating the density profile for bosonic atoms trapped in an optical lattice with an overall harmonic trap at finite temperature and finite on site interaction in the presence of superfluid regions. Our results match well with quantum Monte Carlo simulations at finite temperature. We also show that the superfluid order parameter never vanishes in the trap due to proximity effect. Our calculations for the scaled density in the vacuum to superfluid transition agree well with the experimental data for appropriate temperatures. We present calculations for the entropy per particle as a function of temperature which can be used to calibrate the temperature in experiments. We also discuss issues connected with the demonstration of universal quantum critical scaling in the experiments.
△ Less
Submitted 21 August, 2013; v1 submitted 19 August, 2013;
originally announced August 2013.
-
Doping a correlated band insulator: A new route to half metallic behaviour
Authors:
Arti Garg,
H. R. Krishnamurthy,
Mohit Randeria
Abstract:
We demonstrate in a simple model the surprising result that turning on an on-site Coulomb interaction U in a doped band insulator leads to the formation of a half-metallic state. In the undoped system, we show that increasing U leads to a first order transition between a paramagnetic, band insulator and an antiferomagnetic Mott insulator at a finite value U_{AF}. Upon doping, the system exhibits h…
▽ More
We demonstrate in a simple model the surprising result that turning on an on-site Coulomb interaction U in a doped band insulator leads to the formation of a half-metallic state. In the undoped system, we show that increasing U leads to a first order transition between a paramagnetic, band insulator and an antiferomagnetic Mott insulator at a finite value U_{AF}. Upon doping, the system exhibits half metallic ferrimagnetism over a wide range of doping and interaction strengths on either side of U_{AF}. Our results, based on dynamical mean field theory, suggest a novel route to half-metallic behavior and provide motivation for experiments on new materials for spintronics.
△ Less
Submitted 10 July, 2013;
originally announced July 2013.
-
Theoretical Description of Coherent Doublon Creation via Lattice Modulation Spectroscopy
Authors:
A. Dirks,
K. Mikelsons,
H. R. Krishnamurthy,
J. K. Freericks
Abstract:
Using a recently developed strong-coupling method, we present a comprehensive theory for doublon production processes in modulation spectroscopy of a three-dimensional system of ultracold fermionic atoms in an optical lattice with a trap. The theoretical predictions compare well to the experimental time traces of doublon production. For experimentally feasible conditions, we provide a quantitative…
▽ More
Using a recently developed strong-coupling method, we present a comprehensive theory for doublon production processes in modulation spectroscopy of a three-dimensional system of ultracold fermionic atoms in an optical lattice with a trap. The theoretical predictions compare well to the experimental time traces of doublon production. For experimentally feasible conditions, we provide a quantitative prediction for the presence of a nonlinear "two-photon" excitation at strong modulation amplitudes.
△ Less
Submitted 28 January, 2014; v1 submitted 29 April, 2013;
originally announced April 2013.
-
Quasi-universal transient behavior of a nonequilibrium Mott insulator driven by an electric field
Authors:
K. Mikelsons,
J. K. Freericks,
H. R. Krishnamurthy
Abstract:
We use a self-consistent strong-coupling expansion for the self-energy (perturbation theory in the hopping) to describe the nonequilibrium dynamics of strongly correlated lattice fermions. We study the three-dimensional homogeneous Fermi-Hubbard model driven by an external electric field showing that the damping of the ensuing Bloch oscillations depends on the direction of the field, and that for…
▽ More
We use a self-consistent strong-coupling expansion for the self-energy (perturbation theory in the hopping) to describe the nonequilibrium dynamics of strongly correlated lattice fermions. We study the three-dimensional homogeneous Fermi-Hubbard model driven by an external electric field showing that the damping of the ensuing Bloch oscillations depends on the direction of the field, and that for a broad range of field strengths, a long-lived transient prethermalized state emerges. This long-lived transient regime implies that thermal equilibrium may be out of reach of the time scales accessible in present cold atom experiments, but shows that an interesting new quasi-universal transient state exists in nonequilibrium governed by a thermalized kinetic energy but not a thermalized potential energy. In addition, when the field strength is equal in magnitude to the interaction between atoms, the system undergoes a rapid thermalization, characterized by a different quasi-universal behavior of the current and spectral function for different values of the hopping.
△ Less
Submitted 21 October, 2012;
originally announced October 2012.
-
Spectral moment sum rules for the retarded Green's function and self-energy of the inhomogeneous Bose-Hubbard model in equilibrium and nonequilibrium
Authors:
J. K. Freericks,
V. Turkowski,
H. R. Krishnamurthy,
M. Knap
Abstract:
We derive expressions for the zeroth and the first three spectral moment sum rules for the retarded Green's function and for the zeroth and the first spectral moment sum rules for the retarded self-energy of the inhomogeneous Bose-Hubbard model in nonequilibrium, when the local on-site repulsion and the chemical potential are time-dependent, and in the presence of an external time-dependent electr…
▽ More
We derive expressions for the zeroth and the first three spectral moment sum rules for the retarded Green's function and for the zeroth and the first spectral moment sum rules for the retarded self-energy of the inhomogeneous Bose-Hubbard model in nonequilibrium, when the local on-site repulsion and the chemical potential are time-dependent, and in the presence of an external time-dependent electromagnetic field. We also evaluate these moments for equilibrium, where all time dependence and external fields vanish. Unlike similar sum rules for the Fermi-Hubbard model, in the Bose-Hubbard model case, the sum rules often depend on expectation values that cannot be determined simply from parameters in the Hamiltonian like the interaction strength and chemical potential, but require knowledge of equal time many-body expectation values from some other source. We show how one can approximately evaluate these expectation values for the Mott-insulating phase in a systematic strong-coupling expansion in powers of the hopping divided by the interaction. We compare the exact moments to moments of spectral functions calculated from a variety of different approximations and use them to benchmark their accuracy.
△ Less
Submitted 19 August, 2012;
originally announced August 2012.
-
Synchronous and Asynchronous Mott Transitions in Topological Insulator Ribbons
Authors:
Amal Medhi,
Vijay B. Shenoy,
H. R. Krishnamurthy
Abstract:
We address how the nature of linearly dispersing edge states of two dimensional (2D) topological insulators evolves with increasing electron-electron correlation engendered by a Hubbard like on-site repulsion $U$ in finite ribbons of two models of topological band insulators. Using an inhomogeneous cluster slave rotor mean-field method developed here, we show that electronic correlations drive the…
▽ More
We address how the nature of linearly dispersing edge states of two dimensional (2D) topological insulators evolves with increasing electron-electron correlation engendered by a Hubbard like on-site repulsion $U$ in finite ribbons of two models of topological band insulators. Using an inhomogeneous cluster slave rotor mean-field method developed here, we show that electronic correlations drive the topologically nontrivial phase into a Mott insulating phase via two different routes. In a synchronous transition, the entire ribbon attains a Mott insulating state at one critical $U$ that depends weakly on the width of the ribbon. In the second, asynchronous route, Mott localization first occurs on the edge layers at a smaller critical value of electronic interaction which then propagates into the bulk as $U$ is further increased until all layers of the ribbon become Mott localized. We show that the kind of Mott transition that takes place is determined by certain properties of the linearly dispersing edge states which characterize the topological resilience to Mott localization.
△ Less
Submitted 19 December, 2011;
originally announced December 2011.
-
Two-channel Kondo physics in odd impurity chains
Authors:
Andrew K. Mitchell,
David E. Logan,
H. R. Krishnamurthy
Abstract:
We study odd-membered chains of spin-(1/2) impurities, with each end connected to its own metallic lead. For antiferromagnetic exchange coupling, universal two-channel Kondo (2CK) physics is shown to arise at low energies. Two overscreening mechanisms are found to occur depending on coupling strength, with distinct signatures in physical properties. For strong inter-impurity coupling, a residual c…
▽ More
We study odd-membered chains of spin-(1/2) impurities, with each end connected to its own metallic lead. For antiferromagnetic exchange coupling, universal two-channel Kondo (2CK) physics is shown to arise at low energies. Two overscreening mechanisms are found to occur depending on coupling strength, with distinct signatures in physical properties. For strong inter-impurity coupling, a residual chain spin-(1/2) moment experiences a renormalized effective coupling to the leads; while in the weak-coupling regime, Kondo coupling is mediated via incipient single-channel Kondo singlet formation. We also investigate models where the leads are tunnel-coupled to the impurity chain, permitting variable dot filling under applied gate voltages. Effective low-energy models for each regime of filling are derived, and for even-fillings where the chain ground state is a spin singlet, an orbital 2CK effect is found to be operative. Provided mirror symmetry is preserved, 2CK physics is shown to be wholly robust to variable dot filling; in particular the single-particle spectrum at the Fermi level, and hence the low-temperature zero-bias conductance, is always pinned to half-unitarity. We derive a Friedel-Luttinger sum rule and from it show that, in contrast to a Fermi liquid, the Luttinger integral is non-zero and determined solely by the `excess' dot charge as controlled by gate voltage. The relevance of the work to real quantum dot devices, where inter-lead charge-transfer processes fatal to 2CK physics are present, is also discussed. Physical arguments and numerical renormalization group techniques are used to obtain a detailed understanding of these problems.
△ Less
Submitted 28 June, 2011; v1 submitted 25 March, 2011;
originally announced March 2011.
-
Novel effects of localization due to `intrinsic disorder' in the `two-fluid' model for manganites
Authors:
Prabuddha Sanyal,
V. B. Shenoy,
H. R. Krishnamurthy,
T. V. Ramakrishnan
Abstract:
We discuss the effects of a novel polaronic disorder in the recently proposed two-fluid model for manganites. Using effective field theory as well as direct numerical simulations, we show that this disorder can have dramatic effects in terms of the transition from ferromagnetic insulator to ferromagnetic metal upon hole-doping, including an Anderson localized regime where variable range hopping ma…
▽ More
We discuss the effects of a novel polaronic disorder in the recently proposed two-fluid model for manganites. Using effective field theory as well as direct numerical simulations, we show that this disorder can have dramatic effects in terms of the transition from ferromagnetic insulator to ferromagnetic metal upon hole-doping, including an Anderson localized regime where variable range hopping may be observed.
△ Less
Submitted 7 January, 2011;
originally announced January 2011.
-
Long Range Coulomb Interactions and Nanoscale Electronic Inhomogeneities in Correlated Oxides
Authors:
Vijay B. Shenoy,
Tribikram Gupta,
H. R. Krishnamurthy,
T. V. Ramakrishnan
Abstract:
Electronic, magnetic or structural inhomogeneities ranging in size from nanoscopic to mesoscopic scales seem endemic, and are possibly generic, to colossal magnetoresistance manganites and other transition metal oxides. We show here that an extension, to include long range Coulomb interactions, of a quantum two-fluid $\ell-b$ model proposed recently for manganites [Phys. Rev. Lett., {\bf 92}, 15…
▽ More
Electronic, magnetic or structural inhomogeneities ranging in size from nanoscopic to mesoscopic scales seem endemic, and are possibly generic, to colossal magnetoresistance manganites and other transition metal oxides. We show here that an extension, to include long range Coulomb interactions, of a quantum two-fluid $\ell-b$ model proposed recently for manganites [Phys. Rev. Lett., {\bf 92}, 157203 (2004)] leads to an excellent description of such inhomogeneities. In the $\ell-b$ model two very different kinds of electronic states, one localized and polaronic ($\ell$), and the other extended or broad band ($b$) co-exist. For model parameters appropriate to manganites, and even within a simple dynamical mean-filed theory (DMFT) framework, it describes many of the unusual phenomena seen in manganites, including colossal magnetoresistance (CMR), qualitatively and quantitatively. However, in the absence of long ranged Coulomb interaction, a system described by such a model would actually phase separate, into macroscopic regions of $l$ and $b$ electrons respectively. As we show in this paper, in the presence of Coulomb interactions, the {\em macroscopic} phase separation gets suppressed, and instead nanometer scale regions of polarons interspersed with band electron puddles appear, constituting a new kind of quantum Coulomb glass. Our work points to an interplay of strong correlations, long range Coulomb interaction and dopant ion disorder as the origin of nanoscale inhomogeneities, rather than disorder frustrated phase competition as is generally believed. Based on this, we argue that the observed micrometer(meso)-scale inhomogeneities owe their existence to extrinsic causes, eg. strain due to cracks and defects. We suggest possible experiments to validate our speculation.
△ Less
Submitted 12 April, 2009;
originally announced April 2009.
-
Strong-coupling expansion for the momentum distribution of the Bose Hubbard model with benchmarking against exact numerical results
Authors:
J. K. Freericks,
H. R. Krishnamurthy,
Yasuyuki Kato,
Naoki Kawashima,
Nandini Trivedi
Abstract:
A strong-coupling expansion for the Green's functions, self-energies and correlation functions of the Bose Hubbard model is developed. We illustrate the general formalism, which includes all possible inhomogeneous effects in the formalism, such as disorder, or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the…
▽ More
A strong-coupling expansion for the Green's functions, self-energies and correlation functions of the Bose Hubbard model is developed. We illustrate the general formalism, which includes all possible inhomogeneous effects in the formalism, such as disorder, or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the bosons in the Mott phase for an infinite homogeneous periodic system at zero temperature through third-order in the hopping. By using scaling theory for the critical behavior at zero momentum and at the critical value of the hopping for the Mott insulator to superfluid transition along with a generalization of the RPA-like form for the momentum distribution, we are able to extrapolate the series to infinite order and produce very accurate quantitative results for the momentum distribution in a simple functional form for one, two, and three dimensions; the accuracy is better in higher dimensions and is on the order of a few percent relative error everywhere except close to the critical value of the hopping divided by the on-site repulsion. In addition, we find simple phenomenological expressions for the Mott phase lobes in two and three dimensions which are much more accurate than the truncated strong-coupling expansions and any other analytic approximation we are aware of. The strong-coupling expansions and scaling theory results are benchmarked against numerically exact QMC simulations in two and three dimensions and against DMRG calculations in one dimension. These analytic expressions will be useful for quick comparison of experimental results to theory and in many cases can bypass the need for expensive numerical simulations.
△ Less
Submitted 19 February, 2009;
originally announced February 2009.
-
Synthesis, Structure and Properties of Boron and Nitrogen Doped Graphene
Authors:
L. S. Panchakarla,
K. S. Subrahmanyam,
S. K. Saha,
A. Govindaraj,
H. R. Krishnamurthy,
U. V. Waghmare,
C. N. R. Rao
Abstract:
Two-dimensional graphene exhibits many fascinating properties such as ballistic electronic conduction and quantum Hall effect at room temperature.1-4 Graphene doped electrochemically or through charge-transfer with electron-donor and -acceptor molecules,5-7 shows marked changes in electronic structure, with characteristic signatures in the Raman spectra.5-10 Substitutional doping, universally us…
▽ More
Two-dimensional graphene exhibits many fascinating properties such as ballistic electronic conduction and quantum Hall effect at room temperature.1-4 Graphene doped electrochemically or through charge-transfer with electron-donor and -acceptor molecules,5-7 shows marked changes in electronic structure, with characteristic signatures in the Raman spectra.5-10 Substitutional doping, universally used in tuning properties of semiconductors, could be a powerful tool to control the electronic properties of graphene. Here, we present the structure and properties of boron and nitrogen doped graphenes, obtained by more than one method involving arc discharge in appropriate gaseous atmosphere, by using modified graphite electrode or by the transformation of nano-diamond. Using a combination of experiment and firstprinciples theory, we demonstrate systematic changes in the carrier-concentration and electronic structure of graphenes with B/N-doping, accompanied by stiffening of the Gband and intensification of the defect related D-band in the Raman spectra. Such n/p - type graphenes obtained without external fields or chemical agents should find device applications.
△ Less
Submitted 18 February, 2009;
originally announced February 2009.
-
Perfect dc conductance of a finite width Mott insulator sandwiched between metallic leads at zero temperature: a quantum emergent phenomenon in strongly correlated multilayers
Authors:
H. Zenia,
J. K. Freericks,
H. R. Krishnamurthy,
Th. Pruschke
Abstract:
Using inhomogeneous dynamical mean-field theory, we argue that the normal-metal proximity effect forces any finite number of "barrier" planes that are described by the (paramagnetic) Hubbard model and sandwiched between semi-infinite metallic leads to always be a Fermi liquid at T=0. This then implies that the inhomogeneous system restores lattice periodicity at zero frequency, has a well-define…
▽ More
Using inhomogeneous dynamical mean-field theory, we argue that the normal-metal proximity effect forces any finite number of "barrier" planes that are described by the (paramagnetic) Hubbard model and sandwiched between semi-infinite metallic leads to always be a Fermi liquid at T=0. This then implies that the inhomogeneous system restores lattice periodicity at zero frequency, has a well-defined Fermi surface, and should display perfect (ballistic) conductivity or "transparency". These results are, however, fragile with respect to finite frequency, V, T, disorder, or magnetism, all of which restore the expected quantum tunneling regime through a finite-width Mott insulator. Our formal results are complemented by numerical renormalization group studies on small thickness barriers that illustrate under what circumstances this behavior might be seen in real experimental systems.
△ Less
Submitted 29 September, 2008;
originally announced September 2008.
-
Theoretical description of time-resolved pump/probe photoemission in TaS\_2: a single-band DFT+DMFT(NRG) study within the quasiequilibrium approximation
Authors:
J. K. Freericks,
H. R. Krishnamurthy,
Yizhi Ge,
A. Y. Liu,
Th. Pruschke
Abstract:
In this work, we theoretically examine recent pump/probe photoemission experiments on the strongly correlated charge-density-wave insulator TaS\_2. We describe the general nonequilibrium many-body formulation of time-resolved photoemission in the sudden approximation, and then solve the problem using dynamical mean-field theory with the numerical renormalization group and a bare density of state…
▽ More
In this work, we theoretically examine recent pump/probe photoemission experiments on the strongly correlated charge-density-wave insulator TaS\_2. We describe the general nonequilibrium many-body formulation of time-resolved photoemission in the sudden approximation, and then solve the problem using dynamical mean-field theory with the numerical renormalization group and a bare density of states calculated from density functional theory including the charge-density-wave distortion of the ion cores and spin-orbit coupling
We find a number of interesting results: (i) the bare band structure actually has more dispersion in the perpendicular direction than in the two-dimensional planes; (ii) the DMFT approach can produce upper and lower Hubbard bands that resemble those in the experiment, but the upper bands will overlap in energy with other higher energy bands; (iii) the effect of the finite width of the probe pulse is minimal on the shape of the photoemission spectra; and (iv) the quasiequilibrium approximation does not fully describe the behavior in this system.
△ Less
Submitted 13 September, 2008;
originally announced September 2008.
-
Theoretical description of time-resolved photoemission spectroscopy: application to pump-probe experiments
Authors:
J. K. Freericks,
H. R. Krishnamurthy,
Th. Pruschke
Abstract:
The theory for time-resolved photoemission spectroscopy as applied to pump-probe experiments is developed and solved for the generic case of a strongly correlated material. The formal development incorporates all of the nonequilibrium effects of the pump pulse and the finite time width of the probe pulse. While our formal development is completely general, in our numerical illustration for the H…
▽ More
The theory for time-resolved photoemission spectroscopy as applied to pump-probe experiments is developed and solved for the generic case of a strongly correlated material. The formal development incorporates all of the nonequilibrium effects of the pump pulse and the finite time width of the probe pulse. While our formal development is completely general, in our numerical illustration for the Hubbard model, we assume the pump pulse drives the electrons into a nonequilibrium configuration, which rapidly thermalizes to create a hot (quasi-equilibrium) electronic system, and we then study the effects of windowing that arise from the finite width of the probe pulse. We find sharp features in the spectra are broadened, particularly the quasiparticle peak of strongly correlated metals at low temperature.
△ Less
Submitted 29 June, 2008;
originally announced June 2008.
-
Canonical Trajectories and Critical Coupling of the Bose-Hubbard Hamiltonian in a Harmonic Trap
Authors:
G. G. Batrouni,
H. R. Krishnamurthy,
K. W. Mahmud,
V. G. Rousseau,
R. T. Scalettar
Abstract:
Quantum Monte Carlo (QMC) simulations and the Local Density Approximation (LDA) are used to map the constant particle number (canonical) trajectories of the Bose Hubbard Hamiltonian confined in a harmonic trap onto the $(μ/U,t/U)$ phase diagram of the uniform system. Generically, these curves do not intercept the tips of the Mott insulator (MI) lobes of the uniform system. This observation neces…
▽ More
Quantum Monte Carlo (QMC) simulations and the Local Density Approximation (LDA) are used to map the constant particle number (canonical) trajectories of the Bose Hubbard Hamiltonian confined in a harmonic trap onto the $(μ/U,t/U)$ phase diagram of the uniform system. Generically, these curves do not intercept the tips of the Mott insulator (MI) lobes of the uniform system. This observation necessitates a clarification of the appropriate comparison between critical couplings obtained in experiments on trapped systems with those obtained in QMC simulations. The density profiles and visibility are also obtained along these trajectories. Density profiles from QMC in the confined case are compared with LDA results.
△ Less
Submitted 30 March, 2008; v1 submitted 19 March, 2008;
originally announced March 2008.
-
Spectral properties in the charge density wave phase of the half-filled Falicov-Kimball Model
Authors:
S. R. Hassan,
H. R. Krishnamurthy
Abstract:
We study the spectral properties of charge density wave (CDW) phase of the half-filled spinless Falicov-Kimball model within the framework of the Dynamical Mean Field Theory. We present detailed results for the spectral function in the CDW phase as function of temperature and $U$. We show how the proximity of the non-fermi liquid phase affects the CDW phase, and show that there is a region in th…
▽ More
We study the spectral properties of charge density wave (CDW) phase of the half-filled spinless Falicov-Kimball model within the framework of the Dynamical Mean Field Theory. We present detailed results for the spectral function in the CDW phase as function of temperature and $U$. We show how the proximity of the non-fermi liquid phase affects the CDW phase, and show that there is a region in the phase diagram where we get a CDW phase without a gap in the spectral function. This is a radical deviation from the mean-field prediction where the gap is proportional to the order parameter.
△ Less
Submitted 26 September, 2007;
originally announced September 2007.
-
Unusual doping and temperature dependence of photoemission spectra from manganites
Authors:
Prabuddha Sanyal,
Subhra Sen Gupta,
Nandan Pakhira,
H. R. Krishnamurthy,
D. D. Sarma,
T. V. Ramakrishnan
Abstract:
A recent, major, puzzle in the core-level photoemission spectra of doped manganites is the observation of a 1-2 eV wide shoulder with intensity varying with temperature T as the square of the magnetization over a T scale of order 200K, an order of magnitude less than electronic energies. This is addressed and resolved here, by extending a recently proposed two electron fluid l-b model for these…
▽ More
A recent, major, puzzle in the core-level photoemission spectra of doped manganites is the observation of a 1-2 eV wide shoulder with intensity varying with temperature T as the square of the magnetization over a T scale of order 200K, an order of magnitude less than electronic energies. This is addressed and resolved here, by extending a recently proposed two electron fluid l-b model for these systems to include core-hole effects. The shoulder arises from a rapid redistribution of e_g electron density, as a function of T, between the highly localized (l) and band-like (b) states. Furthermore, our theory leads to a correspondence between spectral changes due to increasing doping and decreasing T, as experimentally observed.
△ Less
Submitted 30 April, 2007;
originally announced April 2007.
-
Frequency shift of optical phonons in doped graphene layers
Authors:
Srijan Kumar Saha,
U. V. Waghmare,
H. R. Krishnamurthy,
A. K. Sood
Abstract:
We use first-principles density-functional calculations to determine the frequency shift of the A$'_1$-${\bf K}$ phonon (Raman D band) in monolayer graphene, as a function of the charge doping. A detailed DFT study on the electron-phonon coupling and the phonon line width of E$_{2g}$-${\bm Γ}$ and A$'_1$-${\bf K}$ phonons are also performed for graphene multi-layers. Furthermore, we explain the…
▽ More
We use first-principles density-functional calculations to determine the frequency shift of the A$'_1$-${\bf K}$ phonon (Raman D band) in monolayer graphene, as a function of the charge doping. A detailed DFT study on the electron-phonon coupling and the phonon line width of E$_{2g}$-${\bm Γ}$ and A$'_1$-${\bf K}$ phonons are also performed for graphene multi-layers. Furthermore, we explain the experimentally observed '1/(Number of Layers)' behaviour of the Raman G band position after including the dynamic response treated within time dependent perturbation theory.
△ Less
Submitted 15 March, 2007; v1 submitted 27 February, 2007;
originally announced February 2007.
-
Dynamics of capacitively coupled double quantum dots
Authors:
Martin R. Galpin,
David E. Logan,
H. R. Krishnamurthy
Abstract:
We consider a double dot system of equivalent, capacitively coupled semiconducting quantum dots, each coupled to its own lead, in a regime where there are two electrons on the double dot. Employing the numerical renormalization group, we focus here on single-particle dynamics and the zero-bias conductance, considering in particular the rich range of behaviour arising as the interdot coupling is…
▽ More
We consider a double dot system of equivalent, capacitively coupled semiconducting quantum dots, each coupled to its own lead, in a regime where there are two electrons on the double dot. Employing the numerical renormalization group, we focus here on single-particle dynamics and the zero-bias conductance, considering in particular the rich range of behaviour arising as the interdot coupling is progressively increased through the strong coupling (SC) phase, from the spin-Kondo regime, across the SU(4) point to the charge-Kondo regime; and then towards and through the quantum phase transition to a charge-ordered (CO) phase. We first consider the two-self-energy description required to describe the broken symmetry CO phase, and implications thereof for the non-Fermi liquid nature of this phase. Numerical results for single-particle dynamics on all frequency scales are then considered, with particular emphasis on universality and scaling of low-energy dynamics throughout the SC phase. The role of symmetry breaking perturbations is also briefly discussed.
△ Less
Submitted 8 August, 2006;
originally announced August 2006.
-
Renormalization group study of capacitively coupled double quantum dots
Authors:
Martin R. Galpin,
David E. Logan,
H. R. Krishnamurthy
Abstract:
The numerical renormalization group is employed to study a double quantum (DQD) dot system consisting of two equivalent single-level dots, each coupled to its own lead and with a mutual capacitive coupling embodied in an interdot interaction U', in addition to the intradot Coulomb interaction U. We focus on the regime with two electrons on the DQD, and the evolution of the system on increasing U…
▽ More
The numerical renormalization group is employed to study a double quantum (DQD) dot system consisting of two equivalent single-level dots, each coupled to its own lead and with a mutual capacitive coupling embodied in an interdot interaction U', in addition to the intradot Coulomb interaction U. We focus on the regime with two electrons on the DQD, and the evolution of the system on increasing U'/U. The spin-Kondo effect arising for U'=0 (SU(2) x SU(2)) is found to persist robustly with increasing U'/U, before a rapid but continuous crossover to (a) the SU(4) point U'=U where charge and spin degrees of freedom are entangled and the Kondo scale strongly enhanced; and then (b) a charge-Kondo state, in which a charge-pseudospin is quenched on coupling to the leads/conduction channels. A quantum phase transition of Kosterlitz-Thouless type then occurs from this Fermi liquid, strong coupling (SC) phase, to a broken symmetry, non-Fermi liquid charge ordered (CO) phase at a critical U'_c. Our emphasis in this paper is on the structure, stability and flows between the underlying RG fixed points, on the overall phase diagram in the (U,U')-plane and evolution of the characteristic low-energy Kondo scale inherent to the SC phase; and on static physical properties such as spin- and charge-susceptibilities (staggered and uniform), including universality and scaling behaviour in the strongly correlated regime. Some exact results for associated Wilson ratios are also obtained.
△ Less
Submitted 7 August, 2006;
originally announced August 2006.