Source author record

N. M. R. Peres

N. M. R. Peres appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

77works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

77 published item(s)

preprint2022arXiv

Contribution to Excitonic Linewidth from Free Carrier-Exciton Scattering in Layered Materials: The example of h-BN

Scattering of excitons by free carriers is a phenomena which is especially important when considering moderately to heavily doped semiconductors in low temperature experiments, where the interaction of excitons with acoustic and optical phonons is reduced. In this paper, we consider the scattering of excitons by free carriers in monolayer hexagonal Boron Nitride encapsulated by a dielectric medium. We describe the excitonic states by variational wave functions, modeling the electrostatic interaction via the Rytova--Keldysh potential. Making the distinction between elastic and inelastic scattering, the relevance of each transition between excitonic states is also considered. Finally, we discuss the contribution of free carrier scattering to the excitonic linewidth, analyzing both its temperature and carrier density dependence.

preprint2022arXiv

Curvature-induced repulsive effect on the lateral Casimir-Polder--van der Waals force

We consider a perfectly conducting infinite cylinder with radius $R$, and investigate the Casimir-Polder (CP) and van der Waals (vdW) interactions with a neutral polarizable particle constrained to move in a plane distant $x_0>R$ from the axis of the cylinder. We show that when the relative curvature $x_0/R \lesssim 6.44$, this particle, under the action of the lateral CP force (which is the projection of the CP force onto the mentioned plane), is attracted to the point on the plane which is closest to the cylinder surface. On the other hand, when $x_0/R \gtrsim 6.44$, we also show that, for certain particle orientations and anisotropy, the lateral CP force can move the particle away from the cylinder. This repulsive behavior of such a component of the CP force reveals a nontrivial dependence of the CP interaction with the surface geometry, specifically of the relative curvature. In the vdW regime, we show that a similar nontrivial repulsive behavior occurs, but for the relative curvature $x_0/R \gtrsim 2.18$, which means that this effect requires a smaller cylinder curvature in the vdW regime than in the CP one. In addition, we also show that there are classical counterparts of these effects, involving a neutral particle with a permanent electric dipole moment. The prediction of such geometric effects on this force may be relevant for a better controlling of the interaction between a particle and a curved surface in classical and quantum physics.

preprint2022arXiv

Excitonic instability in transition metal dichalcogenides

When transition-metal dichalcogenide monolayers lack inversion symmetry, their low-energy single particle spectrum can described by tilted massive Dirac Hamiltonians. The so-called Janus materials fall into that category. Inversion symmetry can also be broken by the application of out-of-plane electric fields, or by the mere presence of a substrate. Here we explore the properties of excitons in TMDC monolayers lacking inversion symmetry. We find that exciton binding energies can be larger than the electronic band gap, making such materials promising candidates to host the elusive exciton insulator phase. We also investigate the excitonic contribution to their optical conductivity and discuss the associated optical selection rules.

preprint2022arXiv

Laser induced enhanced coupling between photons and squeezed magnons in antiferromagnets

In this paper we consider a honeycomb antiferromagnet subject to an external laser field. Obtaining a time-independent effective Hamiltonian, we find that the external laser renormalizes the exchange interaction between the in-plane components of the spin-operators, and induces a synthetic Dzyaloshinskii-Moria interaction (DMI) between second neighbors. The former allows the control of the magnon dispersion's bandwidth and the latter breaks time-reversal symmetry inducing non-reciprocity in momentum space. The eigen-excitations of the system correspond to squeezed magnons whose squeezing parameters depend on the properties of the laser. When studying how these spin excitations couple with cavity photons, we obtain a coupling strength which can be enhanced by an order of magnitude via careful tuning of the laser's intensity, when compared to the case where the laser is absent. The transmission plots through the cavity are presented, allowing the mapping of the magnons' dispersion relation.

preprint2022arXiv

Theoretical methods for excitonic physics in two-dimensional materials

In this tutorial we introduce the reader to several theoretical methods of determining the exciton wave functions and the corresponding eigenenergies. The methods covered are either analytical, semi-analytical, or numeric. We make explicit all the details associated with the different methods, thus allowing newcomers to do research on their own, without experiencing a steep learning curve. The tutorial starts with a variational method and ends with a simple semi-analytical approach to solve the Bethe-Salpeter equation in two-dimensional (2D) gapped materials. For the first methods addressed in this tutorial, we focus on a single layer of hexagonal Boron Nitride (hBN) and of transition metal dichalcogenide (TMD), as these are exemplary materials in the field of 2D excitons. For explaining the Bethe- Salpeter method we choose the biased bilayer graphene, which presents a tunnable band gap. The system has the right amount of complexity (without being excessive). This allows the presentation of the solution of the Bethe-Salpeter equation in a context that can be easily generalized to more complex systems or to apply it to simpler models.

preprint2022arXiv

Tunable Excitons in Rhombohedral Trilayer Graphene

Trilayer graphene is receiving an increasing level of attention due to its stacking--dependent magnetoelectric and optoelectric properties, and its more robust ferromagnetism relative to monolayer and bilayer variants. Additionally, rhombohedral stacked trilayer graphene presents the possibility of easily opening a gap via either an external electric field perpendicular to the layers, or via the application of external strain. In this paper, we consider an external electric field to open a bandgap in rhombohedral trilayer graphene and study the excitonic optical response of the system. This is done via the combination of a tight binding model with the Bethe--Salpeter equation, solved semi--analytically and requiring only a simple numerical quadrature. We then discuss the valley--dependent optical selection rules, followed by the computation of the excitonic linear optical conductivity for the case of a rhombohedral graphene trilayer encapsulated in hexagonal boron nitride. The tunability of the excitonic resonances via an external field is also discussed, together with the increasing localization of the excitonic states as the field increases.

preprint2022arXiv

Two limit cases of twisted hBN bilayers and their excitonic response

In this paper we discuss the optical response due to the excitonic effect of two types of hBN bilayers: AB and AA'. Understanding the properties of these bilayers is of great utility to the study of twisted bilayers at arbitrary angles, since these two configurations correspond to the limit cases of 0 and 60 degree rotation. To obtain the excitonic response we present a method to solve a four-band Bethe-Salpeter equation, by casting it into a 1D problem, thus greatly reducing the numerical burden of the calculation when compared with strictly 2D methods. We find results in good agreement with ab initio calculations already published in the literature for the AA' bilayer, and predict the excitonic conductivity of the AB bilayer, which remains largely unstudied. The main difference in the conductivity of these two types of bilayers is the appearance of a small, yet well resolved, resonance between two larger ones in the AB configuration. This resonance is due to a mainly interlayer exciton, and is absent in the AA' bilayer. Also, the conductivity of the AB bilayer is due to both intralayer and interlayer excitons and is dominated by p-states, while intralayer states are the relevant ones for the AA' configuration, like in a monolayer. The effect of introducing a bias in the AA' bilayer is also discussed.

preprint2021arXiv

Absorption and optical selection rules of tunable excitons in biased bilayer graphene

Biased bilayer graphene, with its easily tunable band gap, presents itself as the ideal system to explore the excitonic effect in graphene based systems. In this paper we study the excitonic optical response of such a system by combining a tight binding model with the solution of the Bethe-Salpeter equation, the latter being solved in a semi-analytical manner, requiring a single numerical quadrature, thus allowing for a transparent calculation. With our approach we start by analytically obtaining the optical selection rules, followed by the computation of the absorption spectrum for the case of a biased bilayer encapsulated in hexagonal boron nitride, a system which has been the subject of a recent experimental study. An excellent agreement is seen when we compare our theoretical prediction with the experimental data.

preprint2021arXiv

Enhancing the hybridization of plasmons in graphene with 2D superconductor collective modes

We explore ways in which the close proximity between graphene sheets and monolayers of 2D superconductors can lead to hybridization between their collective excitations. We consider heterostructures formed by combinations of graphene sheets and 2D superconductor monolayers. The broad range of energies in which the graphene plasmon can exist, together with its tunability, makes such heterostrucutres promising platforms for probing the many-body physics of superconductors. We show that the hybridization between the graphene plasmon and the Bardasis-Schrieffer mode of a 2D superconductor results in clear signatures on the near-field reflection coefficient of the heterostructure, which in principle can be observed in scanning near-field microscopy experiments.

preprint2021arXiv

Harnessing Ultraconfined Graphene Plasmons to Probe the Electrodynamics of Superconductors

We show that the Higgs mode of a superconductor, which is usually challenging to observe by far-field optics, can be made clearly visible using near-field optics by harnessing ultraconfined graphene plasmons. As near-field sources we investigate two examples: graphene plasmons and quantum emitters. In both cases the coupling to the Higgs mode is clearly visible. In the case of the graphene plasmons, the coupling is signaled by a clear anticrossing stemming from the interaction of graphene plasmons with the Higgs mode of the superconductor. In the case of the quantum emitters, the Higgs mode is observable through the Purcell effect. When combining the superconductor, graphene, and the quantum emitters, a number of experimental knobs become available for unveiling and studying the electrodynamics of superconductors.

preprint2020arXiv

Analytical Quantitative Semi-Classical Approach to the LoSurdo-Stark Effect and Ionization in 2D excitons

Using a semi-classical approach, we derive a fully analytical expression for the ionization rate of excitons in two-dimensional materials due to an external static electric field, which eliminates the need for complicated numerical calculations. Our formula shows quantitative agreement with more sophisticated numerical methods based on the exterior complex scaling approach, which solves a non-hermitian eigenvalue problem yielding complex energy eigenvalues, where the imaginary part describes the ionization rate. Results for excitons in hexagonal boron nitride and the $A-$exciton in transition metal dichalcogenides are given as a simple examples. The extension of the theory to include spin-orbit-split excitons in transition metal dichalcogenides is trivial.

preprint2020arXiv

Exciton-polariton mediated interaction between two nitrogen-vacancy color centers in diamond using two-dimensional transition metal dichalcogenides

In this paper, starting from a quantum master equation, we discuss the interaction between two negatively charged Nitrogen-vacancy color centers in diamond via exciton-polaritons propagating in a two-dimensional transition metal dichalcogenide layer in close proximity to a diamond crystal. We focus on the optical 1.945 eV transition and model the Nitrogen-vacancy color centers as two-level (artificial) atoms. We find that the interaction parameters and the energy levels renormalization constants are extremely sensitive to the distance of the Nitrogen-vacancy centers to the transition metal dichalcogenide layer. Analytical expressions are obtained for the spectrum of the exciton-polaritons and for the damping constants entering the Lindblad equation. The conditions for occurrence of exciton mediated superradiance are discussed.

preprint2020arXiv

Excitonic Magneto-Optical Kerr Effect in 2D Transition Metal Dichalcogenides Induced by Spin Proximity

In this paper we develop the excitonic theory of Kerr rotation angle in a two-dimensional (2D) transition metal dichalcogenide at zero magnetic field. The finite Kerr angle is induced by the interplay between spin-orbit splitting and proximity exchange coupling due to the presence of a ferromagnet. We compare the excitonic effect with the single particle theory approach. We show that the excitonic properties of the 2D material lead to a dramatic change in the frequency dependence of the optical response function. We also find that the excitonic corrections enhance the optical response by a factor of two in the case of MoS2 in proximity to a Cobalt thin film.

preprint2020arXiv

Ionisation rate and Stark shift of a one-dimensional model of the Hydrogen molecular ion

In this paper we study the ionization rate and the Stark shift of a one-dimensional model of the H$_{2}^{+}$ ion. The finding of these two quantities is reduced to the solutions of a complex eigenvalue problem. We solve this problem both numerically and analytically. In the latter case we consider the regime of small external electrostatic fields and small internuclear distances. We find an excellent agreement between the ionization rate computed with the two approaches, even when the approximate result is pushed beyond its expected validity. The ionization rate is very sensitive to small changes of the external electrostatic field, spanning many orders of magnitude for small changes of the intensity of the external field. The dependence of the ionization on the internuclear distance is also studied, as this has a direct connection with experimental methods in molecular physics. It is shown that for large distances the ionization rate saturates, which is a direct consequence of the behavior of the energy eigenvalue with the internuclear distance. The Stark shift is computed and from it we extract the static polarizability of H$_{2}^{+}$ and compare our results with those found by other authors using more sophisticated methods.

preprint2020arXiv

Magneto-optical Kerr effect in spin split two-dimensional massive Dirac materials

Two-dimensional (2D) massive Dirac electrons possess a finite Berry curvature, with Chern number $\pm 1/2$, that entails both a quantized dc Hall response and a subgap full-quarter Kerr rotation. The observation of these effects in 2D massive Dirac materials such as gapped graphene, hexagonal boron nitride or transition metal dichalcogenides (TMDs) is obscured by the fact that Dirac cones come in pairs with opposite sign Berry curvatures, leading to a vanishing Chern number. Here, we show that the presence of spin-orbit interactions, combined with an exchange spin splitting induced either by diluted magnetic impurities or by proximity to a ferromagnetic insulator, gives origin to a net magneto-optical Kerr effect in such systems. We focus on the case of TMD monolayers and study the dependence of Kerr rotation on frequency and exchange spin splitting. The role of the substrate is included in the theory and found to critically affect the results. Our calculations indicate that state-of-the-art magneto-optical Kerr spectroscopy can detect a single magnetic impurity in diluted magnetic TMDs.

preprint2020arXiv

Non-reciprocal magnons in a two dimensional crystal with off-plane magnetization

Non reciprocal spin waves have a chiral asymmetry so that their energy is different for two opposite wave vectors. They are found in atomically thin ferromagnetic overlayers with in plane magnetization and are linked to the anti-symmetric Dzyaloshinskii-Moriya surface exchange. We use an itinerant fermion theory based on first principles calculations to predict that non-reciprocal magnons can occur in Fe$_3$GeTe$_2$, the first stand alone metallic two dimensional crystal with off-plane magnetization. We find that both the energy and lifetime of magnons are non-reciprocal and we predict that acoustic magnons can have lifetimes up to hundreds of picoseconds, orders of magnitude larger than in other conducting magnets.

preprint2020arXiv

Topological magnons in CrI$_3$ monolayers: an itinerant fermion description

Magnons dominate the magnetic response of the recently discovered insulating ferromagnetic two dimensional crystals such as CrI$_3$. Because of the arrangement of the Cr spins in a honeycomb lattice, magnons in CrI$_3$ bear a strong resemblance with electronic quasiparticles in graphene. Neutron scattering experiments carried out in bulk CrI$_3$ show the existence of a gap at the Dirac points, that has been conjectured to have a topological nature. Here we propose a theory for magnons in ferromagnetic CrI$_3$ monolayers based on an itinerant fermion picture, with a Hamiltonian derived from first principles. We obtain the magnon dispersion for 2D CrI$_3$ with a gap at the Dirac points with the same Berry curvature in both valleys. For CrI$_3$ ribbons, we find chiral in-gap edge states. Analysis of the magnon wave functions in momentum space further confirms their topological nature. Importantly, our approach does not require to define a spin Hamiltonian, and can be applied to both insulating and conducting 2D materials with any type of magnetic order.

preprint2020arXiv

Topological Photonic Tamm-States and the Su-Schrieffer-Heeger Model

In this paper we study the formation of topological Tamm states at the interface between a semi-infinite one-dimensional photonic-crystal and a metal. We show that when the system is topologically non-trivial there is a single Tamm state in each of the band-gaps, whereas if it is topologically trivial the band-gaps host no Tamm states. We connect the disappearance of the Tamm states with a topological transition from a topologically non-trivial system to a topologically trivial one. This topological transition is driven by the modification of the dielectric functions in the unit cell. Our interpretation is further supported by an exact mapping between the solutions of Maxwell's equations and the existence of a tight-binding representation of those solutions. We show that the tight-binding representation of the 1D photonic crystal, based on Maxwell's equations, corresponds to a Su-Schrieffer-Heeger-type model (SSH-model) for each set of pairs of bands. Expanding this representation near the band edge we show that the system can be described by a Dirac-like Hamiltonian. It allows one to characterize the topology associated with the solution of Maxwell's equations via the winding number. In addition, for the infinite system, we provide an analytical expression for the photonic bands from which the band-gaps can be computed.

preprint2019arXiv

Excitons in Phosphorene: A Semi-Analytical Perturbative Approach

In this paper we develop a semi-analytical perturbation-theory approach to the calculation of the energy levels (binding energies) and wave functions of excitons in phosphorene. Our method gives both the exciton wave function in real and reciprocal spaces with the same ease. This latter aspect is important for the calculation of the nonlinear optical properties of phosphorene. We find that our results are in agreement with calculations based both on the Bethe-Salpeter equation and on Monte Carlo simulations, which are computationally much more demanding. Our approach thus introduces a simple, viable, and accurate method to address the problem of excitons in anisotropic two-dimensional materials.

preprint2019arXiv

Quantization of graphene plasmons

In this article we perform the quantization of graphene plasmons using both a macroscopic approach based on the classical average electromagnetic energy and a quantum hydrodynamic model, in which graphene charge carriers are modeled as a charged fluid. Both models allow to take into account the dispersion of graphenes optical response, with the hydrodynamic model also allowing for the inclusion of non-local effects. Using both methods, the electromagnetic field mode-functions, and the respective frequencies, are determined for two different graphene structures. we show how to quantize graphene plasmons, considering that graphene is a dispersive medium, and taking into account both local and nonlocal descriptions. It is found that the dispersion of graphene's optical response leads to a non-trivial normalization condition for the mode-functions. The obtained mode-functions are then used to calculate the decay of an emitter, represented by a dipole, via the excitation of graphene surface plasmon-polaritons. The obtained results are compared with the total spontaneous decay rate of the emitter and a near perfect match is found in the relevant spectral range. It is found that non-local effects in graphene's conductivity, become relevant for the emission rate for small Fermi energies and small distances between the dipole and the graphene sheet.

preprint2019arXiv

The magnetic Purcell effect: the case of an emitter near an antiferromagnet

In this paper we discuss the magnetic Purcell effect of a magnetic dipole near a semi-infinite antiferromagnet. Contrary to the electric Purcell effect, the magnetic one is not so well studied in the literature. We derive the dispersion relation of the surface wave existing at an antiferromagnetic-dielectric interface from the calculation of the reflection coefficient of the structure. After characterizing the surface wave we quantize the electromagnetic vector potential of the surface wave. This allow us to discuss the magnetic Purcell effect via Fermi golden rule.

preprint2019arXiv

Twisted bilayer graphene: low-energy physics, electronic and optical properties

Van der Waals (vdW) heterostructures ---formed by stacking or growing two-dimensional (2D) crystals on top of each other--- have emerged as a new promising route to tailor and engineer the properties of 2D materials. Twisted bilayer graphene (tBLG), a simple vdW structure where the interference between two misaligned graphene lattices leads to the formation of a moiré pattern, is a test bed to study the effects of the interaction and misalignment between layers, key players for determining the electronic properties of these stackings. In this chapter, we present in a pedagogical way the general theory used to describe lattice mismatched and misaligned vdW structures. We apply it to the study of tBLG in the limit of small rotations and see how the coupling between the two layers leads both to an angle dependent renormalization of graphene's Fermi velocity and appearance of low-energy van Hove singularities. The optical response of this system is then addressed by computing the optical conductivity and the dispersion relation of tBLG surface plasmon-polaritons.

preprint2016arXiv

Graphene field-effect transistor array with integrated electrolytic gates scaled to 200 mm

Ten years have passed since the beginning of graphene research. In this period we have witnessed breakthroughs both in fundamental and applied research. However, the development of graphene devices for mass production has not yet reached the same level of progress. The architecture of graphene field-effect transistors (FET) has not significantly changed, and the integration of devices at the wafer scale has generally not been sought. Currently, whenever an electrolyte-gated FET (EGFET) is used, an external, cumbersome, out-of-plane gate electrode is required. Here, an alternative architecture for graphene EGFET is presented. In this architecture, source, drain, and gate are in the same plane, eliminating the need for an external gate electrode and the use of an additional reservoir to confine the electrolyte inside the transistor active zone. This planar structure with an integrated gate allows for wafer-scale fabrication of high-performance graphene EGFETs, with carrier mobility up to 1800 cm2 V-1 s-1. As a proof-of principle, a chemical sensor was achieved. It is shown that the sensor can discriminate between saline solutions of different concentrations. The proposed architecture will facilitate the mass production of graphene sensors, materializing the potential of previous achievements in fundamental and applied graphene research.

preprint2016arXiv

Graphene Plasmons in Triangular Wedges and Grooves

The ability to effectively guide electromagnetic radiation below the diffraction limit is of the utmost importance in the prospect of all-optical plasmonic circuitry. Here, we propose an alternative solution to conventional metal-based plasmonics by exploiting the deep subwavelength confinement and tunability of graphene plasmons guided along the apex of a graphene-covered dielectric wedge or groove. In particular, we present a quasi-analytic model to describe the plasmonic eigenmodes in such a system, including the complete determination of their spectrum and corresponding induced potential and electric field distributions. We have found that the dispersion of wedge/groove graphene plasmons follows the same functional dependence as their flat-graphene plasmons counterparts, but now scaled by a (purely) geometric factor in which all the information about the system's geometry is contained. We believe our results pave the way for the development of novel custom-tailored photonic devices for subwavelength waveguiding and localization of light based on recently discovered 2D materials.

preprint2016arXiv

Modeling the Excitation of Graphene Plasmons in Periodic Grids of Graphene Ribbons: An Analytical Approach

We study electromagnetic scattering and subsequent plasmonic excitations in periodic grids of graphene ribbons. To address this problem, we develop an analytical method to describe the plasmon-assisted absorption of electromagnetic radiation by a periodic structure of graphene ribbons forming a diffraction grating for THz and mid-IR light. The major advantage of this method lies in its ability to accurately describe the excitation of graphene surface plasmons (GSPs) in one-dimensional (1D) graphene gratings without the use of both time-consuming, and computationally-demanding full-wave numerical simulations. We thus provide analytical expressions for the reflectance, transmittance and plasmon-enhanced absorbance spectra, which can be readily evaluated in any personal laptop with little-to-none programming. We also introduce a semi-analytical method to benchmark our previous results and further compare the theoretical data with spectra taken from experiments, to which we observe a very good agreement. These theoretical tools may therefore be applied to design new experiments and cutting-edge nanophotonic devices based on graphene plasmonics.

preprint2016arXiv

Multiple negative differential conductance regions and inelastic phonon assisted tunneling in graphene-hBN-graphene structures

In this paper we study in detail the effect of the rotational alignment between a hexagonal boron nitride (hBN) slab and the graphene layers in the vertical current of a a graphene-hBN-graphene device. We show how for small rotational angles, the transference of momentum by the hBN crystal lattice leads to multiple peaks in the I-V curve of the device, giving origin to multiple regions displaying negative differential conductance. We also study the effect of scattering by phonons in the vertical current an see how the opening up of inelastic tunneling events allowed by spontaneous emission of optical phonons leads to sharp peaks in the second derivative of the current.

preprint2016arXiv

Numerical calculation of the Casimir-Polder interaction between a graphene sheet with vacancies and an atom

In this work the Casimir{Polder interaction energy between a rubidium atom and a disordered graphene sheet is investigated beyond the Dirac cone approximation by means of accurate real-space calculations. As a model of defected graphene, we consider a tight-binding model of Π-electrons on a honeycomb lattice with a small concentration of point defects. The optical response of the graphene sheet is evaluated with full spectral resolution by means of exact Chebyshev polynomial expansions of the Kubo formula in large lattices with in excess of ten million atoms. At low temperatures, the optical response of defected graphene is found to display two qualitatively distinct behavior with a clear transition around non-zero Fermi energy, μ~μ*. In the vicinity of the Dirac point, the imaginary part of optical conductivity is negative for low frequencies while the real part is strongly suppressed. On the other hand, for high doping, it has the same features found in the Drude model within the Dirac cone approximation, namely, a Drude peak at small frequencies and a change of sign in the imaginary part above the interband threshold ω> 2μ. These characteristics translate into a non-monotonic behavior of the Casimir{Polder interaction energy with very small variation with doping in the vicinity of the neutrality point while having the same form of the interaction calculated with Drude's model at high electronic density.

preprint2016arXiv

Pumping electrons in graphene to the $\mathbf{M}$-point in the Brillouin zone: The emergence of anisotropic plasmons

We consider the existence of plasmons in a non-equilibrium situation where electrons from the valence band of graphene are pumped to states in the Brillouin zone around the $\mathbf{M}$-point by a high intensity UV electromagnetic field. The resulting out-of-equilibrium electron gas is later probed by a weak electromagnetic field of different frequency. We show that the optical properties of the system and the dispersion of the plasmons are strongly anisotropic, depending on the pumping radiation properties: its intensity, polarization, and frequency. This anisotropy has its roots in the saddle-like nature of the electronic dispersion relation around that particular point in the Brillouin zone. It is found that despite the strong anisotropy, the dispersion of the plasmons scales with the square root of the wave number but is characterized an effective Fermi energy, which depends on the properties of the pumping radiation. Our calculations go beyond the usual Dirac cone approximation taking the full band structure of graphene into account. This is a necessary condition for discussing plasmons at the $\mathbf{M}$-point in the Brillouin zone.

preprint2016arXiv

Scattering of surface plasmon-polaritons in a graphene multilayer photonic crystal with inhomogeneous doping

The propagation of a surface plasmon-polariton along a stack of doped graphene sheets is considered. This auxiliary problem is used to discuss: (i) the scattering of such a mode at an interface between the stack and the vacuum; (ii) the scattering at an interface where there is a sudden change of the electronic doping. The formalism is then extended to the {\it barrier problem}. In this system rich physics is found for the plasmonic mode, showing: total reflection, total transmission, Fabry-Pérot oscillations, and coupling to photonic modes.

preprint2015arXiv

Active magneto-optical control of spontaneous emission in graphene

We investigate the spontaneous emission rate of a two-level quantum emitter near a graphene-coated substrate under the influence of an external magnetic field or strain induced pseudo-magnetic field. We demonstrate that the application of the magnetic field can substantially increase or decrease the decay rate. We show that a suppression as large as 99$\%$ in the Purcell factor is achieved even for moderate magnetic fields. The emitter's lifetime is a discontinuous function of $|{\bf B}|$, which is a direct consequence of the occurrence of discrete Landau levels in graphene. We demonstrate that, in the near-field regime, the magnetic field enables an unprecedented control of the decay pathways into which the photon/polariton can be emitted. Our findings strongly suggest that a magnetic field could act as an efficient agent for on-demand, active control of light-matter interactions in graphene at the quantum level.

preprint2015arXiv

Anderson localization of light in disordered superlattices containing graphene layers

We theoretically investigate light propagation and Anderson localization in one-dimensional disordered superlattices composed of dielectric stacks with graphene sheets in between. Disorder is introduced either on graphene material parameters ({\it e.g.} Fermi energy) or on the widths of the dielectric stacks. We derive an analytic expression for the localization length $ξ$, and compare it to numerical simulations using transfer matrix technique; a very good agreement is found. We demonstrate that the presence of graphene may strongly attenuate the anomalously delocalised Breswter modes, and is at the origin of a periodic dependence of $ξ$ on frequency, in contrast to the usual asymptotic decay, $ξ\propto ω^{-2}$. By unveiling the effects of graphene on Anderson localization of light, we pave the way for new applications of graphene-based, disordered photonic devices in the THz spectral range.

preprint2015arXiv

Cloaking Resonant Scatterers and Tuning Electron Flow in Graphene

We consider resonant scatterers with large scattering cross-sections in graphene that are produced by a gated disk or a vacancy, and show that a gated ring can be engineered to produce an efficient electron cloak. We also demonstrate that this same scheme can be applied to tune the direction of electron flow. Our analysis is based on a partial-wave expansion of the electronic wave-functions in the continuum approximation, described by the Dirac equation. Using a symmetrized version of the massless Dirac equation, we derive a general condition for the cloaking of a scatterer by a potential with radial symmetry. We also perform tight-binding calculations to show that our findings are robust against the presence of disorder in the gate potential.

preprint2015arXiv

Electrically tunable resonant scattering in fluorinated bilayer graphene

Adatom-decorated graphene offers a promising new path towards spintronics in the ultrathin limit. We combine experiment and theory to investigate the electronic properties of dilutely fluorinated bilayer graphene, where the fluorine adatoms covalently bond to the top graphene layer. We show that fluorine adatoms give rise to resonant impurity states near the charge neutrality point of the bilayer, leading to strong scattering of charge carriers and hopping conduction inside a field-induced band gap. Remarkably, the application of an electric field across the layers is shown to tune the resonant scattering amplitude from fluorine adatoms by nearly twofold. The experimental observations are well explained by a theoretical analysis combining Boltzmann transport equations and fully quantum-mechanical methods. This paradigm can be generalized to many bilayer graphene-adatom materials, and we envision that the realization of electrically tunable resonance may be a key advantage in graphene-based spintronic devices.

preprint2015arXiv

Electronic transport across linear defects in graphene

We investigate the low-energy electronic transport across grain boundaries in graphene ribbons and infinite flakes. Using the recursive Green's function method, we calculate the electronic transmission across different types of grain boundaries in graphene ribbons. We show results for the charge density distribution and the current flow along the ribbon. We study linear defects at various angles with the ribbon direction, as well as overlaps of two monolayer ribbon domains forming a bilayer region. For a class of extended defect lines with periodicity 3, an analytic approach is developed to study transport in infinite flakes. This class of extended grain boundaries is particularly interesting, since the $K$ and $K'$ Dirac points are superposed.

preprint2015arXiv

Renormalization of nanoparticle polarizability in the vicinity of a graphene-covered interface

We study the electromagnetic properties of a metamaterial consisting of polarizable (nano)particles and a single graphene sheet placed at the interface between two dielectrics. We show that the particle's polarizability is renormalized because of the electromagnetic coupling to surface plasmons supported by graphene, which results in a dispersive behavior, different for the polarizability components corresponding to the induced dipole moment, parallel and perpendicular to the graphene sheet. In particular, this effect is predicted to take place for a metallic particle whose bare polarizability in the terahertz (THz) region is practically equal to the cube of its radius (times $4π\varepsilon_0$). This opens the possibility to excite surface plasmons in graphene and enhance its absorption in the THz range by simply using a monolayer of metallic particles randomly deposited on top of it, as we show by explicit calculations.

preprint2015arXiv

Strain induced edge magnetism at zigzag edge in graphene quantum dot

We study the temperature dependent magnetic susceptibility of a strained graphene quantum dot by using the determinant quantum Monte Carlo method. Within the Hubbard model on a honeycomb lattice, our unbiased numerical results show that a relative small interaction $U$ may lead to a edge ferromagnetic like behavior in the strained graphene quantum dot, and a possible room temperature transition is suggested. Around half filling, the ferromagnetic fluctuations at the zigzag edge is strengthened both markedly by the on-site Coulomb interaction and the strain, especially in low temperature region. The resultant strongly enhanced ferromagnetic like behavior may be important for the development of many applications.

preprint2014arXiv

Optical bistability of graphene in the terahertz range

We use an exact solution of the relaxation-time Boltzmann equation in a uniform AC electric field to describe the nonlinear optical response of graphene in the terahertz (THz). The cases of monolayer, bilayer and ABA-stacked trilayer graphene are considered, and the monolayer species is shown to be the most appropriate one to exploit the nonlinear free electron response. We find that a single layer of monolayer graphene shows an optical bistability in the THz range, within the electromagnetic power range attainable in practice. The current associated with the third harmonic generation is also computed.

preprint2014arXiv

Orbital symmetry fingerprints for magnetic adatoms in graphene

In this paper, we describe the formation of local resonances in graphene in the presence of magnetic adatoms containing localized orbitals of arbitrary symmetry, corresponding to any given angular momentum state. We show that quantum interference effects which are naturally inbuilt in the honeycomb lattice in combination with the specific orbital symmetry of the localized state lead to the formation of fingerprints in differential conductance curves. In the presence of Jahn-Teller distortion effects, which lift the orbital degeneracy of the adatoms, the orbital symmetries can lead to distinctive signatures in the local density of states. We show that those effects allow scanning tunneling probes to characterize adatoms and defects in graphene.

preprint2014arXiv

Terahertz response of patterned epitaxial graphene

We study the interaction between polarized terahertz (THz) radiation and micro-structured large-area graphene in transmission geometry. In order to efficiently couple the radiation into the two-dimensional material, a lateral periodic patterning of a closed graphene sheet by intercalation doping into stripes is chosen, yielding unequal transmittance of the radiation polarized parallel and perpendicular to the stripes. Indeed, a polarization contrast up to 20% is observed. The effect even increases up to 50% when removing graphene stripes in analogy to a wire grid polarizer. The polarization dependence is analyzed in a large frequency range from < 80 GHz to 3 THz, including the plasmon-polariton resonance. The results are in excellent agreement with theoretical calculations based on the electronic energy spectrum of graphene and the electrodynamics of the patterned structure.

preprint2013arXiv

A Primer on Surface Plasmon-Polaritons in Graphene

We discuss the properties of surface plasmons-polaritons in graphene and describe three possible ways of coupling electromagnetic radiation in the terahertz (THz) spectral range to this type of surface waves. (i) the attenuated total reflection (ATR) method using a prism in the Otto configuration, (ii) graphene micro-ribbon arrays or monolayers with modulated conductivity, (iii) a metal stripe on top of the graphene layer, and (iv) graphene-based gratings. The text provides a number of original results along with their detailed derivation and discussion.

preprint2013arXiv

Scattering by linear defects in graphene: a tight-binding approach

We develop an analytical scattering formalism for computing the transmittance through periodic defect lines within the tight-binding model of graphene. We first illustrate the method with a relatively simple case, the pentagon-only defect line. Afterwards, more complex defect lines are treated, namely the zz(558) and the zz(5757) ones. The formalism developed, only uses simple tight-binding concepts, reducing the problem to matrix manipulations which can be easily worked out by any computational algebraic calculator.

preprint2013arXiv

Strong light-matter interaction in systems described by a modified Dirac equation

The bulk states of some materials, such as topological insulators, are described by a modified Dirac equation. Such systems may have trivial and non-trivial phases. In this paper, we show that in the non-trivial phase a strong light-matter interaction exists in a two-dimensional system, which leads to an optical conductivity at least one order of magnitude larger than that of graphene.

preprint2012arXiv

Complete light absorption in graphene-metamaterial corrugated structures

We show that surface-plasmon polaritons excited in negative permittivity metamaterials having shallow periodic surface corrugation profiles can be explored to push the absorption of single and continuous sheets of graphene up to 100%. In the relaxation regime, the position of the plasmonic resonances of the hybrid system is determined by the plasma frequency of the metamaterial, allowing the frequency range for enhanced absorption to be set without the need of engineering graphene.

preprint2012arXiv

Confined magneto-optical waves in graphene

The electromagnetic mode spectrum of single-layer graphene subjected to a quantizing magnetic field is computed taking into account intraband and interband contributions to the magneto-optical conductivity. We find that a sequence of weakly decaying quasi-transverse-electric modes, separated by magnetoplasmon polariton modes, emerge due to the quantizing magnetic field. The characteristics of these modes are tuneable, by changing the magnetic field or the Fermi energy.

preprint2012arXiv

Continuum Model of the Twisted Bilayer

The continuum model of the twisted graphene bilayer (Phys. Rev. Lett. 99, 256802, 2007) is extended to include all types of commensurate structures. The essential ingredient of the model, the Fourier components of the spatially modulated hopping amplitudes, can be calculated analytically, for any type of commensurate structures in the low twist angle limit. We show that the Fourier components that could give rise to a gap in the SE-even structures discussed by Mele (Phys. Rev. B 81, 161405 2010) vanish linearly with angle, whereas the amplitudes that saturate to finite values, as $θ\to0$, ensure that all low angle structures share essentially the same physics. We extend our previous calculations beyond the validity of perturbation theory, to discuss the disappearance of Dirac cone structure at angles below 1^{\circ}.

preprint2012arXiv

Efficient graphene-based photodetector with two cavities

We present an efficient graphene-based photodetector with two Fabri-Pérot cavities. It is shown that the absorption can reach almost 100% around a given frequency, which is determined by the two-cavity lengths. It is also shown that hysteresis in the absorbance is possible, with the transmittance amplitude of the mirrors working as an external driving field. The role of non-linear contributions to the optical susceptibility of graphene is discussed.

preprint2012arXiv

Enhanced Optical Dichroism of Graphene Nanoribbons

The optical conductivity of graphene nanoribbons is analytical and exactly derived. It is shown that the absence of translation invariance along the transverse direction allows considerable intra-band absorption in a narrow frequency window that varies with the ribbon width, and lies in the THz range domain for ribbons 10-100nm wide. In this spectral region the absorption anisotropy can be as high as two orders of magnitude, which renders the medium strongly dichroic, and allows for a very high degree of polarization (up to ~85) with just a single layer of graphene. The effect is resilient to level broadening of the ribbon spectrum potentially induced by disorder. Using a cavity for impedance enhancement, or a stack of few layer nanoribbons, these values can reach almost 100%. This opens a potential prospect of employing graphene ribbon structures as efficient polarizers in the far IR and THz frequencies.

preprint2012arXiv

Exact solution for square-wave grating covered with graphene: Surface plasmon-polaritons in the THz range

We provide an analytical solution to the problem of scattering of electromagnetic radiation by a square-wave grating with a flat graphene sheet on top. We show that for deep groves there is a strong plasmonic response with light absorption in the graphene sheet reaching more than 45%, due to the excitation of surface plasmon-polaritons. The case of grating with a graphene sheet presenting an induced periodic modulation of the conductivity is also discussed.

preprint2012arXiv

Faraday effect in graphene enclosed in an optical cavity and the equation of motion method for the study of magneto-optical transport in solids

We show that by enclosing graphene in an optical cavity, giant Faraday rotations in the infrared regime are generated and measurable Faraday rotation angles in the visible range become possible. Explicit expressions for the Hall steps of the Faraday rotation angle are given for relevant regimes. In the context of this problem we develop an equation of motion (EOM) method for calculation of the magneto-optical properties of metals and semiconductors. It is shown that properly regularized EOM solutions are fully equivalent to the Kubo formula.

preprint2012arXiv

Light scattering by a medium with a spatially modulated optical conductivity: the case of graphene

We describe light scattering from a graphene sheet having a modulated optical conductivity. We show that such modulation enables the excitation of surface plasmon-polaritons by an electromagnetic wave impinging at normal incidence. The resulting surface plasmon-polaritons are responsible for a substantial increase of electromagnetic radiation absorption by the graphene sheet. The origin of the modulation can be due either to a periodic strain field or to adatoms (or absorbed molecules) with a modulated adsorption profile.

preprint2012arXiv

Scattering by linear defects in graphene: a continuum approach

We study the low-energy electronic transport across periodic extended defects in graphene. In the continuum low-energy limit, such defects act as infinitesimally thin stripes separating two regions where Dirac Hamiltonian governs the low-energy phenomena. The behavior of these systems is defined by the boundary condition imposed by the defect on the massless Dirac fermions. We demonstrate how this low-energy boundary condition can be computed from the tight-binding model of the defect line. For simplicity we consider defect lines oriented along the zigzag direction, which requires the consideration of only one copy of Dirac equation. Three defect lines of this kind are studied and shown to be mappable between them: the pentagon-only, the zz(558) and the zz(5757) defect lines. In addition, in this same limit, we calculate the conductance across such defect lines with size L, and find it to be proportional to k_FL at low temperatures.

preprint2011arXiv

Coulomb Drag and High Resistivity Behavior in Double Layer Graphene

We show that Coulomb drag in ultra-clean graphene double layers can be used for controlling the on/off ratio for current flow by tunning the external gate voltage. Hence, although graphene remains semi-metallic, the double layer graphene system can be tuned from conductive to a highly resistive state. We show that our results explain previous data of Coulomb drag in double layer graphene samples in disordered SiO2 substrates.

preprint2011arXiv

Field-effect tunneling transistor based on vertical graphene heterostructures

We report a bipolar field effect tunneling transistor that exploits to advantage the low density of states in graphene and its one atomic layer thickness. Our proof-of-concept devices are graphene heterostructures with atomically thin boron nitride acting as a tunnel barrier. They exhibit room temperature switching ratios ~50, a value that can be enhanced further by optimizing the device structure. These devices have potential for high frequency operation and large scale integration.

preprint2011arXiv

Mechanism for graphene-based optoelectronic switches by tuning surface plasmon-polaritons in monolayer graphene

It is shown that one can explore the optical conductivity of graphene, together with the ability of controlling its electronic density by an applied gate voltage, in order to achieve resonant coupling between an external electromagnetic radiation and surface plasmon-polaritons in the graphene layer. This opens the possibility of electrical control of the intensity of light reflected inside a prism placed on top of the graphene layer, by switching between the regimes of total reflection and total absorption. The predicted effect can be used to build graphene-based opto-electronic switches.

preprint2011arXiv

Observation of Intra- and Inter-band Transitions in the Optical Response of Graphene

The optical conductivity of freely suspended graphene was examined under non-equilibrium conditions using femtosecond pump-probe spectroscopy. We observed a conductivity transient that varied strongly with the electronic temperature, exhibiting a crossover from enhanced to decreased absorbance with increasing pump fluence. The response arises from a combination of bleaching of the inter-band transitions by Pauli blocking and induced absorption from the intra-band transitions of the carriers. The latter dominates at low electronic temperature, but, despite an increase in Drude scattering rate, is overwhelmed by the former at high electronic temperature. The time-evolution of the optical conductivity in all regimes can described in terms of a time-varying electronic temperature.

preprint2011arXiv

Solution of the quantum harmonic oscillator plus a delta-function potential at the origin: The oddness of its even-parity solutions

We derive the energy levels associated with the even-parity wave functions of the harmonic oscillator with an additional delta-function potential at the origin. Our results bring to the attention of students a non-trivial and analytical example of a modification of the usual harmonic oscillator potential, with emphasis on the modification of the boundary conditions at the origin. This problem calls the attention of the students to an inaccurate statement in quantum mechanics textbooks often found in the context of solution of the harmonic oscillator problem.

preprint2011arXiv

Stability of boron nitride bilayers: Ground state energies, interlayer distances, and tight-binding description

We have studied boron nitride monolayer and bilayer band structures. For bilayers, the ground state energies of the different five stackings are computed using DFT in order to determine the most stable configuration. Also, the interlayer distance for the five different types of stacking in which boron-nitride bilayers can be found is determined. Using a minimal tight binding model for the band structures of boron nitride bilayers, the hopping parameters and the onsite energies have been extracted by fitting a tight binding empirical model to the DFT results.

preprint2011arXiv

Transport properties of graphene with one-dimensional charge defects

We study the effect of extended charge defects in electronic transport properties of graphene. Extended defects are ubiquitous in chemically and epitaxially grown graphene samples due to internal strains associated with the lattice mismatch. We show that at low energies these defects interact quite strongly with the 2D Dirac fermions and have an important effect in the DC-conductivity of these materials.

preprint2011arXiv

Unified description of the dc conductivity of monolayer and bilayer graphene at finite densities based on resonant scatterers

We show that a coherent picture of the dc conductivity of monolayer and bilayer graphene at finite electronic densities emerges upon considering that strong short-range potentials are the main source of scattering in these two systems. The origin of the strong short-range potentials may lie in adsorbed hydrocarbons at the surface of graphene. The equivalence among results based on the partial-wave description of scattering, the Lippmann-Schwinger equation, and the T-matrix approach is established. Scattering due to resonant impurities close to the neutrality point is investigated via a numerical computation of the Kubo formula using a kernel polynomial method. We find that relevant adsorbate species originate impurity bands in monolayer and bilayer graphene close to the Dirac point. In the midgap region, a plateau of minimum conductivity of about $e^2/h$ (per layer) is induced by the resonant disorder. In bilayer graphene, a large adsorbate concentration can develop an energy gap between midgap and high-energy states. As a consequence, the conductivity plateau is supressed near the edges and a "conductivity gap" takes place. Finally, a scattering formalism for electrons in biased bilayer graphene, taking into account the degeneracy of the spectrum, is developed and the dc conductivity of that system is studied.

preprint2011arXiv

Zigzag graphene nanoribbon edge reconstruction with Stone-Wales defects

In this article, we study zigzag graphene nanoribbons with edges reconstructed with Stone-Wales defects, by means of an empirical (first-neighbor) tight-binding method, with parameters determined by ab-initio calculations of very narrow ribbons. We explore the characteristics of the electronic band structure with a focus on the nature of edge states. Edge reconstruction allows the appearance of a new type of edge states. They are dispersive, with non-zero amplitudes in both sub-lattices; furthermore, the amplitudes have two components that decrease with different decay lengths with the distance from the edge; at the Dirac points one of these lengths diverges, whereas the other remains finite, of the order of the lattice parameter. We trace this curious effect to the doubling of the unit cell along the edge, brought about by the edge reconstruction. In the presence of a magnetic field, the zero-energy Landau level is no longer degenerate with edge states as in the case of pristine zigzag ribbon.

preprint2010arXiv

Analytical study of non-linear transport across a semiconductor-metal junction

In this paper we study analytically a one-dimensional model for a semiconductor-metal junction. We study the formation of Tamm states and how they evolve when the semi-infinite semiconductor and metal are coupled together. The non-linear current, as a function of the bias voltage, is studied using the non-equilibrium Green's function method and the density matrix of the interface is given. The electronic occupation of the sites defining the interface has strong non-linearities as function of the bias voltage due to strong resonances present in the Green's functions of the junction sites. The surface Green's function is computed analytically by solving a quadratic matrix equation, which does not require adding a small imaginary constant to the energy. The wave function for the surface states is given.

preprint2010arXiv

Bilayer graphene: gap tunability and edge properties

Bilayer graphene -- two coupled single graphene layers stacked as in graphite -- provides the only known semiconductor with a gap that can be tuned externally through electric field effect. Here we use a tight binding approach to study how the gap changes with the applied electric field. Within a parallel plate capacitor model and taking into account screening of the external field, we describe real back gated and/or chemically doped bilayer devices. We show that a gap between zero and midinfrared energies can be induced and externally tuned in these devices, making bilayer graphene very appealing from the point of view of applications. However, applications to nanotechnology require careful treatment of the effect of sample boundaries. This being particularly true in graphene, where the presence of edge states at zero energy -- the Fermi level of the undoped system -- has been extensively reported. Here we show that also bilayer graphene supports surface states localized at zigzag edges. The presence of two layers, however, allows for a new type of edge state which shows an enhanced penetration into the bulk and gives rise to band crossing phenomenon inside the gap of the biased bilayer system.

preprint2010arXiv

Colloquium: The transport properties of graphene: An introduction

An introduction to the transport properties of graphene combining experimental results and theoretical analysis is presented. In the theoretical description simple intuitive models are used to illustrate important points on the transport properties of graphene. The concept of chirality, stemming from the massless Dirac nature of the low energy physics of the material, is shown to be instrumental in understanding its transport properties: the conductivity minimum, the electronic mobility, the effect of strain, the weak (anti-)localization, and the optical conductivity.

preprint2010arXiv

Dynamical polarizability of graphene beyond the Dirac cone approximation

We compute the dynamical polarizability of graphene beyond the usual Dirac cone approximation, integrating over the full Brillouin zone. We find deviations at $\hbarω=2t$ ($t$ the hopping parameter) which amount to a logarithmic singularity due to the van Hove singularity and derive an approximate analytical expression. Also at low energies, we find deviations from the results obtained from the Dirac cone approximation which manifest themselves in a peak spitting at arbitrary direction of the incoming wave vector $\q$. Consequences for the plasmon spectrum are discussed.

preprint2010arXiv

Electronic properties of a biased graphene bilayer

We study, within the tight-binding approximation, the electronic properties of a graphene bilayer in the presence of an external electric field applied perpendicular to the system -- \emph{biased bilayer}. The effect of the perpendicular electric field is included through a parallel plate capacitor model, with screening correction at the Hartree level. The full tight-binding description is compared with its 4-band and 2-band continuum approximations, and the 4-band model is shown to be always a suitable approximation for the conditions realized in experiments. The model is applied to real biased bilayer devices, either made out of SiC or exfoliated graphene, and good agreement with experimental results is found, indicating that the model is capturing the key ingredients, and that a finite gap is effectively being controlled externally. Analysis of experimental results regarding the electrical noise and cyclotron resonance further suggests that the model can be seen as a good starting point to understand the electronic properties of graphene bilayer. Also, we study the effect of electron-hole asymmetry terms, as the second-nearest-neighbor hopping energies $t'$ (in-plane) and $γ_{4}$ (inter-layer), and the on-site energy $Δ$.

preprint2010arXiv

Excitonic effects in the optical conductivity of gated graphene

We study the effect of electron-electron interactions in the optical conductivity of graphene under applied bias and derive a generalization of Elliot's formula, commonly used for semiconductors, for the optical intensity. We show that {\it excitonic resonances} are responsible for several features of the experimentally measured mid-infrared response of graphene such as the increase of the conductivity beyond the "universal" value above the Fermi blocked regime, the broadening of the absorption at the threshold, and the decrease of the optical conductivity at higher frequencies. Our results are also in agreement with {\it ab initio} calculations in the neutral regime.

preprint2010arXiv

Optical Properties of Strained Graphene

The optical conductivity of graphene strained uniaxially is studied within the Kubo-Greenwood formalism. Focusing on inter-band absorption, we analyze and quantify the breakdown of universal transparency in the visible region of the spectrum, and analytically characterize the transparency as a function of strain and polarization. Measuring transmittance as a function of incident polarization directly reflects the magnitude and direction of strain. Moreover, direction-dependent selection rules permit identification of the lattice orientation by monitoring the van-Hove transitions. These photoelastic effects in graphene can be explored towards atomically thin, broadband optical elements.

preprint2007arXiv

Fermi Liquid Theory of a Fermi Ring

We study the effect of electron-electron interactions in the electronic properties of a biased graphene bilayer. This system is a semiconductor with conduction and valence bands characterized by an unusual ``mexican-hat'' dispersion. We focus on the metallic regime where the chemical potential lies in the ``mexican-hat'' in the conduction band, leading to a topologically non-trivial Fermi surface in the shape of a ring. We show that due to the unusual topology of the Fermi surface electron-electron interactions are greatly enhanced. We discuss the possibility of an instability towards a ferromagnetic phase due to this enhancement. We compute the electronic polarization function in the random phase approximation and show that, while at low energies the system behaves as a Fermi liquid (albeit with peculiar Friedel oscillations), at high frequencies it shows a highly anomalous response when compare to ordinary metals.

preprint2007arXiv

Graphene bilayer with a twist: electronic structure

Electronic properties of bilayer and multilayer graphene have generally been interpreted in terms of AB or Bernal stacking. However, it is known that many types of stacking defects can occur in natural and synthetic graphite; rotation of the top layer is often seen in scanning tunneling microscopy (STM) studies of graphite. In this paper we consider a graphene bilayer with a relative small angle rotation between the layers and calculate the electronic structure near zero energy in a continuum approximation. Contrary to what happens in a AB stacked bilayer and in accord with observations in epitaxial graphene we find: (a) the low energy dispersion is linear, as in a single layer, but the Fermi velocity can be significantly smaller than the single layer value; (b) an external electric field, perpendicular to the layers, does not open an electronic gap

preprint2006arXiv

Electronic properties of graphene multilayers

We study the effects of disorder in the electronic properties of graphene multilayers, with special focus on the bilayer and the infinite stack. At low energies and long wavelengths, the electronic self-energies and density of states exhibit behavior with divergences near half-filling. As a consequence, the spectral functions and conductivities do not follow Landau's Fermi liquid theory. In particular, we show that the quasiparticle decay rate has a minimum as a function of energy, there is a universal minimum value for the in-plane conductivity of order e^2/h per plane and, unexpectedly, the c-axis conductivity is enhanced by disorder at low doping, leading to an enormous conductivity anisotropy at low temperatures.

preprint2004arXiv

Phase diagram and magnetic collective excitations of the Hubbard model in graphene sheets and layers

We discuss the magnetic phases of the Hubbard model for the honeycomb lattice both in two and three spatial dimensions. A ground state phase diagram is obtained depending on the interaction strength U and electronic density n. We find a first order phase transition between ferromagnetic regions where the spin is maximally polarized (Nagaoka ferromagnetism) and regions with smaller magnetization (weak ferromagnetism). When taking into account the possibility of spiral states, we find that the lowest critical U is obtained for an ordering momentum different from zero. The evolution of the ordering momentum with doping is discussed. The magnetic excitations (spin waves) in the antiferromagnetic insulating phase are calculated from the random-phase-approximation for the spin susceptibility. We also compute the spin fluctuation correction to the mean field magnetization by virtual emission/absorpion of spin waves. In the large $U$ limit, the renormalized magnetization agrees qualitatively with the Holstein-Primakoff theory of the Heisenberg antiferromagnet, although the latter approach produces a larger renormalization.

preprint2003arXiv

Spin waves in La_2CuO_4: band structure and correlation effects

We calculate the antiferromagnetic spin wave dispersion in the half-filled (electronic density n=1) Hubbard model for a two-dimensional square lattice, using the random phase approximation (RPA) in a broken symmetry (spin density wave) ground state. Our results for the spin wave dispersion, ω(\vec q), are compared with high-resolution inelastic neutron scattering performed on La_2CuO$_4. The effects of different band structures and different values of the on-site Coulomb interaction on the spin wave spectrum is studied. Particular attention is put on the high energy dispersion values ω(π/2,π/2) and ω(0,π).

preprint2002arXiv

Numerical and Monte Carlo Bethe ansatz method: 1D Heisenberg model

In this paper we present two new numerical methods for studying thermodynamic quantities of integrable models. As an example of the effectiveness of these two approaches, results from numerical solutions of all sets of Bethe ansatz equations, for small Heisenberg chains, and Monte Carlo simulations in quasi-momentum space, for a relatively larger chains, are presented. Our results agree with those obtained by thermodynamics Bethe ansatz (TBA) and Quantum Transfer Matrix (QTM).

preprint2000arXiv

Superconductivity in the SU(N) Anderson Lattice at U=\infty

We present a mean-field study of superconductivity in a generalized N-channel cubic Anderson lattice at U=\infty taking into account the effect of a nearest-neighbor attraction J. The condition U=\infty is implemented within the slave-boson formalism considering the slave bosons to be condensed. We consider the $f$-level occupancy ranging from the mixed valence regime to the Kondo limit and study the dependence of the critical temperature on the various model parameters for each of three possible Cooper pairing symmetries (extended s, d-wave and p-wave pairing) and find interesting crossovers. It is found that the d- and p- wave order parameters have, in general, very similar critical temperatures. The extended s-wave pairing seems to be relatively more stable for electronic densities per channel close to one and for large values of the superconducting interaction J.

preprint1999arXiv

Finite-Temperature Transport in Finite-Size Hubbard Rings in the Strong-Coupling Limit

We study the current, the curvature of levels, and the finite temperature charge stiffness, D(T,L), in the strongly correlated limit, U>>t, for Hubbard rings of L sites, with U the on-site Coulomb repulsion and t the hopping integral. Our study is done for finite-size systems and any band filling. Up to order t we derive our results following two independent approaches, namely, using the solution provided by the Bethe ansatz and the solution provided by an algebraic method, where the electronic operators are represented in a slave-fermion picture. We find that, in the U=\infty case, the finite-temperature charge stiffness is finite for electronic densities, n, smaller than one. These results are essencially those of spinless fermions in a lattice of size L, apart from small corrections coming from a statistical flux, due to the spin degrees of freedom. Up to order t, the Mott-Hubbard gap is Δ_{MH}=U-4t, and we find that D(T) is finite for n<1, but is zero at half-filling. This result comes from the effective flux felt by the holon excitations, which, due to the presence of doubly occupied sites, is renormalized to Φ^{eff}=ϕ(N_h-N_d)/(N_d+N_h), and which is zero at half-filling, with N_d and N_h being the number of doubly occupied and empty lattice sites, respectively. Further, for half-filling, the current transported by any eigenstate of the system is zero and, therefore, D(T) is also zero.

preprint1996arXiv

Pseudoparticle Description of the 1D Hubbard Model Electronic Transport Properties

We extend the pseudoparticle transport description of the Hubbard chain to all energy scales. In particular we compute the mean value of the electric current transported by any Bethe-ansatz state and the transport masses of the charge carriers. We present numerical results for the optical conductivity of the model at half-filling for values of U/t=3 and 4. We show that these are in good agreement with the pseudoparticle description of the finite-energy transitions involving new pseudoparticle energy bands.