BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. A. Jafari

S. A. Jafari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.mes-hallcond-mat.mtrl-scicond-mat.supr-congr-qchep-thcond-mat.dis-nncond-mat.othercond-mat.quant-gascond-mat.stat-mechmath-phmath.MPphysics.opticsquant-ph

Catalog footprint

What is connected

33works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

33 published item(s)

preprint2022arXiv

Undamped inter-valley paramagnons in doped graphene

We predict the existence of an undamped collective spin excitation in doped graphene in the paramagnetic regime, referred to as paramagnons. Since the electrons and the holes involved in this collective mode reside in different valleys of the band structure, the momentum of these inter-valley paramagnons is given by the separation of the valleys in momentum space. The energy of the inter-valley paramagnons lies in the void region below the continuum of inter-band single-particle electron-hole excitations that appears when graphene is doped. The paramagnons are undamped due to the lack of electron-hole excitations in this void region. Their energy strongly depends on doping concentration, which can help to identify them in future experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Generalization of Lieb-Wu wave function inspired by one-dimensional ionic Hubbard model

With the ionic Hubbard model (IHM) in mind, we construct a non-trivial generalization of the Bethe ansatz (BA) wave function which naturally generalizes the Lieb-Wu wave function with an ionic parameter $Δ$, and reduces to Lieb-Wu solution in the limit $Δ\to 0$. The resulting two-particle scattering matrix satisfies the Yang-Baxter equation. To the extent that the unit cells with more than two electrons (Choy-Haldane issue) are avoided on average, our wave function represents an effective soluiton for the one-dimensional IHM. The Choy-Haldane issue restricts the validity of our solution to low-filling and large $U\gtrsim 4$. This regime is attainable in cold atom realizations of the IHM. For this regime, we numerically solve the generalized Bethe equations and compute the ground state energy in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Synthetic non-Abelian gauge fields and gravitomagnetic effects in tilted Dirac cone systems

In planar tilted Dirac cone systems, the tilt parameter can be made space-dependent by either a perpendicular displacement field, or by chemical substitution in certain systems. We show that the symmetric partial derivative of the tilt parameter generates non-Abelian synthetic gauge fields in these systems. The small velocity limit of these gauge forces corresponds to Rashba and Dresselhaus spin-orbit couplings. At the classical level, the same symmetric spatial derivatives of tilt contribute to conservative, Lorentz-type and friction-like forces. The velocity dependent forces are odd with respect to tilt and therefore have opposite signs in the two valleys when the system is inversion symmetric. Furthermore, toggling the chemical potential between the valence and conduction bands reverses the sign of the all these classical forces, which indicates these forces couple to the electric charge of the carriers. As such, these forces are natural extensions of the electric and magnetic forces in the particular geometry of the tilted Dirac cone systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Undamped transverse electric mode in undoped two-dimensional tilted Dirac cone materials

Transverse electric (TE) modes can not propagate through the conducting solids. This is because the continuum of particle-hole excitations of conductors contaminates with the TE mode and dampes it out. But in solids hosting tilted Dirac cone (TDC) that admit a description in terms of a modified Minkowski spacetime, the new spacetime structure remedies this issue and therefore a tilted Dirac cone material (TDM) supports the propagation of an undamped TE mode which is sustained by density fluctuations. The resulting TE mode propagates at fermionic velocities which strongly confines the mode to the surface of the two-dimensional (2D) TDM.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Electrically charged Andreev modes in two-dimensional tilted Dirac cone systems

In a graphene-based Josephson junction, the Andreev reflection can become specular which gives rise to propagating Andreev modes. These propagating Andreev modes are essentially charge neutral and therefore they transfer energy but not electric charge. One main result of this work is that when the Dirac theory of graphene is deformed into a tilted Dirac cone, the breaking of charge conjugation symmetry of the Dirac equation renders the resulting Andreev modes electrically charged. We calculate an otherwise zero charge conductance arising solely from the tilt parameters $\vecζ=(ζ_x,ζ_y)$. The distinguishing feature of such a form of charge transport from the charge transport caused by normal electrons is their dependence on the phase difference $ϕ$ of the two superconductors which can be experimentally extracted by employing a flux bias. Another result concerns the enhancement of Josephson current in a regime where instead of propagating Andreev modes, localized Andreev levels are formed. In this regime, we find enhancement by orders of magnitude of the Josephson current when the tilt parameter is brought closer and closer to $ζ=1$ limit. We elucidate that, the enhancement is due to a combination of two effects: (i) enhancement of number of transmission channels by flattening of the band upon tilting to $ζ\approx 1$, and (ii) a non-trivial dependence on the angle $θ$ of the the tilt vector $\vecζ$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Optical conductivity of triple point fermions

As a low-energy effective theory on non-symmorphic lattices, we consider a generic triple point fermion Hamiltonian which is parameterized by an angular parameter $λ$. We find strong $λ$ dependence in both Drude and interband optical absorption of these systems. The deviation of the $T^2$ coefficient of the Drude weight from Dirac/Weyl fermions can be used as a quick way to optically distinguish the triple point degeneracies from the Dirac/Weyl degeneracies. At the particular $λ=π/6$ point, we find that the "helicity" reversal optical transition matrix element is identically zero. But deviating from this point, the helicity reversal emerges as an absorption channel.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Perpendicular Andreev reflection: Solid state signature of black hole horizon

Beenakker noticed that the peculiar band structure of Dirac fermions in 2D solids allows for the specular Andreev reflection in these systems which has no analogue in other 2D electron systems. An interesting deformation of the Dirac equation in the solid state is to tilt it which has now materials realization. In this work we report another peculiar feature of the Andreev reflection in tilted 2D Dirac cone systems. The tilt deformation of the Dirac equation is characterized by two parameters ${\boldsymbolζ}=ζ(\cosθ,\sinθ)$. We show that when the tilt parameter is tuned to its "horizon value" $ζ=1$, irrespective of the incidence angle of electrons, the Andreev reflected hole is always reflected {\em perpendicular} to the interface. Furthermore at the horizon value, if the tilt direction is perpendicular to the interface, the conductance will be energy independent for the entire sub gap energies. For generic tilt $\boldsymbolζ$, the tilt dependence of the conductance line shape can be used to extract information about the direction and magnitude of the tilt.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exactly solvable spin chain models corresponding to BDI class of topological superconductors

We present an exactly solvable extension of the quantum XY chain with longer range multi-spin interactions. Topological phase transitions of the model are classified in terms of the number of Majorana zero modes which are in turn related to an integer winding number. We further find a general relation between the winding number and the number of Majorana fermions at the ends of an open chain. The present class of exactly solvable models belong to the BDI class in the Altland-Zirnbauer classification (A. Altland, M. R. Zirnbauer, Phys. Rev. B (1997) 55 1142) of topological superconductors. We show that time reversal (TR) symmetry of the spin variables translates into a peculiar particle-hole transformation in the language of Jordan-Wigner (JW) fermions that is accompanied by a π shift in the wave vector ($π$PH). Presence of $π$PH symmetry restricts the winding number of TR symmetric extensions of XY to odd integers. The πPH operator may serve in further detailed classification of topological superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantitative measure for the spin-charge separation in two dimensional Hubbard model

We introduce a qauantitative measure of spin-charge separation, $ζ(t)$ which is based on the difference between the fluctuations with respect to background of the spin and charge profiles at any time t and is suitable for studying the non-equilibrium dynamics of excitations in strongly correlated systems. This quantity is not only a direct measure of the spin-charge separation in strongly correlated systems, but its long time behaviour can further serve as a possible order parameter for the interaction induced (Mott) insulating state. Within the nu- merically exact diagonzalization we calculate this quantity for the two dimensional Hubbard model away from Half filling. Our quantitative measure in chain, ladder and two-dimensional geometries gives the same order of magnitude for the quantity of spin-charge separation. Furthermore from the temporal behaviour of $ζ(t)$ a threshold time can be identified that provides clues onto the breakdown of underlying Mott insulating phase.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Localized magnetic states in three dimensional Dirac solids

Formation of localized magnetic states in a metallic host is a classic problem ofcondensed matter physics formalized by P. W. Anderson within the so called single impurity Anderson model (SIAM). The general picture in a host of a simple one-band metal is that a large Hubbard $U$ in the impurity orbital is pre-requisite for the formation of localized magnetic states. In recent years three dimensional (3D) Dirac solids have emerged the hallmark of which is strong spin-orbit interaction. In this work we show that such a strong spin-orbit interaction allows to form localized magnetic states even with small values of Hubbard $U$. This opens up the fascinating possibility of forming magnetic states with $s$ or $p$ orbital impurities -- different from traditional paradigms of $d$ or $f$ orbital based magnetic moments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strong coupling approach to Mott transition of massless and massive Dirac fermions on honeycomb lattice

Phase transitions in the Hubbard model and ionic Hubbard model at half-filling on the honeycomb lattice are investigated in the strong coupling perturbation theory which corresponds to an expansion in powers of the hopping $t$ around the atomic limit. Within this formulation we find analytic expressions for the single-particle spectrum, whereby the calculation of the insulating gap is reduced to a simple root finding problem. This enables high precision determination of the insulating gap that does not require any extrapolation procedure. The critical value of Mott transition on the honeycomb lattice is obtained to be $U_c\approx 2.38 t$. Studying the ionic Hubbard model at the lowest order, we find two insulating states, one with Mott character at large $U$ and another with single-particle gap character at large ionic potential, $Δ$. The present approach gives a critical gapless state at $U=2Δ$ at lowest order. By systematically improving on the perturbation expansion, the density of states around this critical gapless phase reduces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gapless chiral excitons in thin films of topological insulators

In a nanoscopic thin film of a strong topological insulator (TI) the Coulomb interaction in the channel that exchanges the two electrons with the same chirality in two different planes of the slab takes advantage of the minus sign resulting from such "exchange" and gives rise to a bound state between the positive energy states in one surface and the negative energy states in the opposite surfaces. Therefore particle and hole pairs in the {\em undoped} Dirac cone of the TI thin film form an inter-surface spin-singlet state that lies below the continuum of free particle-hole pairs. This mode is similar to the excitons of semiconductors, albeit formed between the electron and hole pairs from two different two-dimensional surfaces. For low-momenta the dispersion relation characterizing this collective mode is linear. Experimental comparison of two slabs with different thicknesses can capture the exponential dependence of the present effect on the slab thickness.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum transport through 3D Dirac materials

Bismuth and its alloys provide a paradigm to realize three dimensional materials whose low-energy effective theory is given by Dirac equation in 3+1 dimensions. We study the quantum transport properties of three dimensional Dirac materials within the framework of Landauer-Büttiker formalism. Charge carriers in normal metal satisfying the Schrödinger equation, can be split into four-component with appropriate matching conditions at the boundary with the three dimensional Dirac material (3DDM). We calculate the conductance and the Fano factor of an interface separating 3DDM from a normal metal, as well as the conductance through a slab of 3DDM. Under certain circumstances the 3DDM appears transparent to electrons hitting the 3DDM. We find that electrons hitting the metal-3DDM interface from metallic side can enter 3DDM in a reversed spin state as soon as their angle of incidence deviates from the the direction perpendicular to interface. However the presence of a second interface completely cancels this effect.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stable local moments of vacancies and hollow-site impurities in graphene

Taking into account the possibility of a p-wave hybridization function $V(\vec k)$ of ad-atom with Dirac electrons in graphene -- which arises for vacancies and hollow-site impurities -- we study the nature of magnetic moment formation within the single impurity Anderson model (SIAM). Compared to the s-wave hybridization function, we find that the local moments formed within the Hartree mean field are robust against the change in the parameters of the model. Further we investigate the stability of the local moments with respect to quantum fluctuations by going beyond the Hartree approximation. We find that for parameter regimes where local moments formed by top-site ad-atoms are completely washed out by quantum fluctuations, those formed by vacancies (or hollow-site impurities) survive the quantum fluctuations captured by post-Hartree approximation. Hence vacancies and hollow-site ad-atoms are suitable candidates to produce stable local moments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Valence Bond Phases in $S=1/2$ Kane-Mele-Heisenberg Model

The phase diagram of Kane-Mele-Heisenberg (KMH) model in classical limit~\cite{zare}, contains disordered regions in the coupling space, as the result of to competition among different terms in the Hamiltonian, leading to frustration in finding a unique ground state. In this work we explore the nature of these phase in the quantum limit, for a $S=1/2$. Employing exact diagonalization (ED) in $S_z$ and nearest neighbor valence bond (NNVB) bases, bond and plaquette valence bond mean field theories, We show that the disordered regions are divided into ordered quantum states in the form of plaquette valence bond crystal(PVBC) and staggered dimerized (SD) phases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kondo resonance from vacancies in graphene

Using the slave-rotor mean-field theory, we study the formation of Kondo resonance in graphene induced from vacancies. At the Dirac neutrality point, we find acritical hybridization strength beyond which the Kondo resonance takes place, despite vanishing density of states. We also find that the line-shapes of the Kondo peak are entirely different from Lorentzian shown in normal metallic hosts. The dependence of Kondo temperature on the effective hybridization parameter of the spinon field turns out to be linear in contrast to the quadratic dependence in normal metals. Upon doping we find strong electron-hole asymmetry for the dependence of Kondo temperature on the chemical potential.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

RKKY interaction in heavily vacant graphene

Dirac electrons in clean graphene can mediate the interactions between two localized magnetic moments. The functional form of the RKKY interaction in pristine graphene is specified by two main features: (i) an atomic scale oscillatory part determined by a wave vector $\vec Q$ connecting the two valleys. Furthermore with doping another longer range oscillation appears which arise from the existence of an extended Fermi surface characterized by a single momentum scale $k_F$. (ii) $R^α$ decay in large distances where the exponent $α=-3$ is a distinct feature of undoped Dirac sea (with a linear dispersion relation) in two dimensions. In this work, we investigate the effect of a few percent vacancies on the above properties. Depending on the doping level, if the chemical potential lies on the linear part of the density of states, the exponent $α$ remains close to -3. Otherwise $α$ reduces towards more negative values which means that the combined effect of vacancies and the randomness in their positions makes it harder for the carriers of the medium to mediate the magnetic interaction. Addition of a few percent of vacancies diminishes the atomic scale oscillations of the RKKY interaction signaling the destruction of two-valley structure of the parent graphene material. Surprisingly by allowing the chemical potential to vary, we find that the longer-range oscillations expected to arise from the existence of a $k_F$ scale in the vacant graphene are absent. This may indicate possible non-Fermi liquid behavior by "alloying" graphene with vacancies. The complete absence of oscillations in heavily vacant graphene can be considered an advantage for applications as a uniform sign of the exchange interaction is desirable for magnetic ordering.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collective excitations and the nature of Mott transition in undoped gapped graphene

Particle-hole continuum (PHC) for massive Dirac fermions in presence of short range interactions, provides an unprecedented opportunity for formation of two collective split-off states, one in the singlet and the other in the triplet (spin-1) channel in undoped system. Both poles are close in energy and are separated from thec continuum of free particle-hole excitations by an energy scale of the order of gap parameter $Δ$. They both disperse linearly with two different velocities reminiscent of spin-charge separation in Luttinger liquids. When the strength of Hubbard interactions is stronger than a critical value, the velocity of singlet excitation which we interpret as a charge boson composite becomes zero, and renders the system a Mott insulator. Beyond this critical point, the low-energy sector is left with a linearly dispersing triplet mode -- a characteristic of a Mott insulator. The velocity of triplet mode at the Mott criticality is twice the velocity of underlying Dirac fermions. The phase transition line in the space of $U$ and $Δ$ is in qualitative agreement with a more involved dynamical mean field theory (DMFT) calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear optical response in gapped graphene

We present a formulation for the nonlinear optical response in gapped graphene, where the low-energy single-particle spectrum is modeled by massive Dirac theory. As a representative example of the formulation presented here, we obtain closed form formula for the third harmonic generation (THG) in gapped graphene. It turns out that the covariant form of the low-energy theory gives rise to a peculiar logarithmic singularities in the nonlinear optical spectra. The universal functional dependence of the response function on dimension-less quantities indicates that the optical nonlinearity can be largely enhanced by tuning the gap to smaller values.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Collective Spin and Charge Excitations in Planar Aromatic Molecules

Employing high accuracy fixed node diffusion Monte Carlo (DMC) method we calculated the lowest triplet collective excitation (spin gap), as well as an upper bound for the singlet excitations (charge gap) in a series of charge neutral planar non-ladder aromatic compounds. Both excitation energies lie below the continuum of particle-hole excitation energies obtained from Hartree-Fock orbitals. Hence they can be interpreted as genuine bound states in the particle-hole channel. Assuming a resonating valence bond (RVB) ground state which has been recently suggested for $sp^2$ bonded systems [ M. Marchi, {\em et. al.}, Phys. Rev. Lett. {\bf 107}, 086807 (2011)], offers a unified description of both excited states as two-spinon and doublon-holon bound states. We corroborate our interpretation, by Exact diagonalization study of a minimal model on finite honeycomb clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Equations of motion method for triplet excitation operators in graphene

Particle-hole continuum in Dirac sea of graphene has a unique window underneath,which in principle leaves a room for bound state formation in the triplet particle hole channel [Phys. Rev. Lett. {\bf 89}, 016402 (2002)]. In this work, we construct appropriate triplet particle-hole operators, and using a repulsive Hubbard type effective interaction, we employ equations of motions to derive approximate eigen-value equation for such triplet operators. While the secular equation for the spin density fluctuations gives rise to an equation which is second order in the strength of the short range interaction, the explicit construction of the triplet operators obtained here shows that in terms of these operators, the second order can be factorized to two first order equations, one of which gives rise to a solution below the particle-hole continuum of Dirac electrons in undoped graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Josephson super-current in graphene-superconductor junction

Within the tunneling Hamiltonian formulation for the eight-component spinors,the Josephson critical super-current has been calculated in a planar superconductor-normal graphene-superconductor junction. Coupling between superconductor regions and graphene is taken into account by a tunneling Hamiltonian which contains two types of tunneling, intra-valley and inter-valley tunneling. Within the present tunneling approach, we find that the contributions of two kinds of tunneling to the critical super-current, are completely separable. Therefore, it is possible to consider the effect of the inter-valley tunnelings in the critical super-current. The incorporation of these type of processes into the tunneling Hamiltonian, exposes a special feature of the graphene Josephson junctions. The effect of inter-valley tunneling appears in the length dependence plot of critical current in the form of oscillations. We also present the results for temperature dependence of critical super-current and compare with experimental results and other theoretical calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Neutral Triplet Collective Mode in Doped Graphene

Particle-hole continuum in Dirac sea of graphene has a unique window underneath, which provides a unique opportunity for emergence of a pole in the susceptibility of the {\em triplet} particle-hole channel in the entire Brillouin zone (BZ). Here we use random phase approximation (RPA) to study such collective mode at zero temperature, in a single layer of doped graphene. We find that due to the chiral nature of one-particle states, in undoped graphene, the wave function overlap factors do not lead to qualitative differences, while in doped graphene they will kill small momentum part of the branch of magnetic excitations by pushing it to touch the lower part of the continuum. The pole corresponding to magnetic excitations survives for for larger momenta in the BZ.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Plaquette valence-bond ordering in J_1-J_2 Heisenberg antiferromagnet on the honeycomb lattice

We study S=1/2 Heisenberg model on the honeycomb lattice with first and second neighbor antiferromagnetic exchange (J_{1}-J_{2} model), employing exact diagonalization in both S_z=0 basis and nearest neighbor singlet valence bond (NNVB) basis. We find that for 0.2<J_2/J_1<0.3, NNVB basis gives a proper description of the ground state in comparison with the exact results. By analyzing the dimer-dimer as well as plaquette-plaquette correlations and also defining appropriate structure factors, we investigate possible symmetry breaking states as the candidate for the ground state in the frustrated region. We provide body of evidences in favor of plaquette valence bond ordering for 0.2<J_2/J_1<0.3. By further increasing the ratio J_2/J_1, this state undergoes a transition to the staggered dimerized state.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Short range Coulomb correlations render massive Dirac fermions massless

Tight binding electrons on a honeycomb lattice are described by an effective Dirac theory at low energies. Lowering symmetry by an alternate ionic potential ($Δ$) generates a single-particle gap in the spectrum. We employ the dynamical mean field theory (DMFT) technique, to study the effect of on-site electron correlation ($U$) on massive Dirac fermions. For a fixed mass parameter $Δ$, we find that beyond a critical value $U_{c1}(Δ)$ massive Dirac fermions become massless. Further increasing $U$ beyond $U_{c2}(Δ)$, there will be another phase transition to the Mott insulating state. Therefore the competition between the single-particle gap parameter, $Δ$, and the Hubbard $U$ restores the semi-metallic nature of the parent Hamiltonian. The width of the intermediate semi-metallic regime shrinks by increasing the ionic potential. However, at small values of $Δ$, there is a wide interval of $U$ values for which the system remains semi-metal.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Anderson Transition in Disordered Bilayer Graphene

Employing the Kernel Polynomial method (KPM), we study the electronic properties of the graphene bilayers in the presence of diagonal disorder, within the tight-binding approximation. The KPM method enables us to calculate local density of states (LDOS) without need to exactly diagonalize the Hamiltonian. We use the geometrical averaging of the LDOS's at different lattice sites as a criterion to distinguish the localized states from extended ones. We find that bilayer graphene undergoes Anderson metal-insulator transition at a critical value of disorder strength.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Excitation Spectrum of One-dimensional Extended Ionic Hubbard Model

We use Perturbative Continuous Unitary Transformations (PCUT) to study the one dimensional Extended Ionic Hubbard Model (EIHM) at half-filling in the band insulator region. The extended ionic Hubbard model, in addition to the usual ionic Hubbard model, includes an inter-site nearest-neighbor (n.n.) repulsion, $V$. We consider the ionic potential as unperturbed part of the Hamiltonian, while the hopping and interaction (quartic) terms are treated as perturbation. We calculate total energy and ionicity in the ground state. Above the ground state, (i) we calculate the single particle excitation spectrum by adding an electron or a hole to the system. (ii) the coherence-length and spectrum of electron-hole excitation are obtained. Our calculations reveal that for V=0, there are two triplet bound state modes and three singlet modes, two anti-bound states and one bound state, while for finite values of $V$ there are four excitonic bound states corresponding to two singlet and two triplet modes. The major role of on-site Coulomb repulsion $U$ is to split singlet and triplet collective excitation branches, while $V$ tends to pull the singlet branches below the continuum to make them bound states.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Substitutional doping of Cu in diamond: Mott physics with $p$ orbitals

Discovery of superconductivity in the impurity band formed by heavy doping of boron into diamond (C:B) as well as doping of boron into silicon (Si:B) has provided a rout for the possibility of new families of superconducting materials. Motivated by the special role played by copper atoms in high temperature superconducting materials where essentially Cu $d$ orbitals are responsible for a variety of correlation induced phases, in this paper we investigate the effect of substitutional doping of Cu into diamond. Our extensive first principle calculations averaged over various geometries based on density functional theory, indicates the formation of a mid-gap band, which mainly arises from the $t_{2g}$ and $4p$ states of Cu. For impurity concentrations of more than $\sim 1%, the effect of $2p$ bands of neighboring carbon atoms can be ignored. Based on our detailed analysis, we suggest a two band model for the mid-gap states consisting of a quarter-filled hole like $t_{2g}$ band, and a half-filled band of $4p$ states. Increasing the concentration of the Cu impurity beyond $\sim 5%, completely closes the spectral gap of the host diamond.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Classical Analogue of the Ionic Hubbard Model

In our earlier work [M. Hafez, {\em et al.}, Phys. Lett. A {\bf 373} (2009) 4479] we employed the flow equation method to obtain a classic effective model from a quantum mechanical parent Hamiltonian called, the ionic Hubbard model (IHM). The classical ionic Hubbard model (CIHM) obtained in this way contains solely Fermionic occupation numbers of two species corresponding to particles with $\up$ and $\down$ spin, respectively. In this paper, we employ the transfer matrix method to analytically solve the CIHM at finite temperature in one dimension. In the limit of zero temperature, we find two insulating phases at large and small Coulomb interaction strength, $U$, mediated with a gap-less metallic phase, resulting in two continuous metal-insulator transitions. Our results are further supported with Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fermi surface nesting and possibility of orbital ordering in FeO

We study FeO, a Mott insulator in GGA and GGA+U approximations. In the GGA we find a multi-band metallic state with remarkable inter-band nesting between two $t_{2g}$ bands of Fermi surface, which signals possible instability towards an orbital ordered insulating phase. Such broken symmetry state, although has lower energy than the underlying homogeneous metallic state, but the gap magnitude is less than the experimentally observed optical gap. Therefore we incorporate the calculated value of on-site Coulomb repulsion U on orbital ordered state. We find that symmetry breaking and Coulomb correlations cooperate together to stabilize the system and give an insulating orbital ordered state, with the gap magnitude very close to the experimental value. We propose this method as a possible indication of orbital ordering in LDA and GGA calculations. We check our method with known examples of LiVO$_2$ and LaMnO$_3$.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Flow Equations for the Ionic Hubbard Model

Taking the site-diagonal terms of the one-dimensional ionic Hubbard model (IHM) as $H_0$, we employ Continuous Unitary Transformations (CUT) to obtain a "classical" effective Hamiltonian in which hopping term has been integrated out. For this Hamiltonian spin gap and charge gap are calculated at half-filling and subject to periodic boundary conditions. Our calculations indicate two transition points. In fixed $Δ$, as $U$ increases from zero, there is a region in which both spin gap and charge gap are positive and identical; characteristic of band insulators. Upon further increasing $U$, first transition occurs at $U=U_{c_{1}}$, where spin and charge gaps both vanish and remain zero up to $U=U_{c_{2}}$. A gap-less state in charge and spin sectors characterizes a metal. For $U>U_{c_{2}}$ spin gap remains zero and charge gap becomes positive. This third region corresponds to a Mott insulator in which charge excitations are gaped, while spin excitations remain gap-less.

0 citations0 reviewsNo code linkedOpen access