Source author record

Steven G. Louie

Steven G. Louie 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

91works
13topics
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

91 published item(s)

preprint2023arXiv

Evidence for excitonic insulator ground state in triangulene Kagome lattice

Electron-hole pair excitations in semiconductors have been predicted to be able to give rise to a highly correlated many-body ground state, the excitonic insulator (EI). Under appropriate conditions below a critical temperature (Tc), strongly bound electron-hole pairs spontaneously form and undergo a phase transition from a normal band insulator into an exciton condensate, transforming the parent material into a novel correlated insulator. Despite recent advances in spectroscopic tools, clear direct experimental evidence for the EI state has been elusive and is often obfuscated by accompanying electronic effects. Here we present the reticular bottom-up synthesis of a Kagome lattice of [4]triangulene, a two-dimensional (2D) covalent organic framework (COF) imbued with a deliberate excitonic instability excitons with binding energies larger than the bandgap arising from a pair of flat bands (FBs). Theoretical analyses based on first-principles calculations and scanning tunnelling spectroscopy (STS) reveal quasiparticle spectral signatures mixing valence (VB) and conduction (CB) characteristics of the FBs along with a non-trivial semiconducting gap that can only be explained by invoking many-body theory. Our findings spectroscopically corroborate the nature of a FB induced exciton insulator ground state and provide a robust yet highly tuneable platform for the exploration of correlated quasi boson physics in quantum materials.

preprint2022arXiv

Exciton-Driven Renormalization of Quasiparticle Band Structure in Monolayer MoS2

Optical excitation serves as a powerful approach to control the electronic structure of layered Van der Waals materials via many-body screening effects, induced by photoexcited free carriers, or via light-driven coherence, such as optical Stark and Bloch-Siegert effects. Although theoretical work has also pointed to an exotic mechanism of renormalizing band structure via excitonic correlations in bound electron-hole pairs (excitons), experimental observation of such exciton-driven band renormalization and the full extent of their implications is still lacking, largely due to the limitations of optical probes and the impact of screening effects. Here, by using extreme-ultraviolet time-resolved angle-resolved photoemission spectroscopy together with excitonic many-body theoretical calculations, we directly unmask the band renormalization effects driven by excitonic correlations in a monolayer semiconductor. We revealed a surprising bandgap opening, increased by 40 meV, and a simultaneous enhancement of band effective mass. Our findings unmask the novel exciton-driven mechanism towards the band engineering in photoexcited semiconducting materials, opening a new playground to manipulate the transient energy states in layered quantum materials via optical controls of excitonic many-body correlations.

preprint2022arXiv

Many-Body Effects in the X-ray Absorption Spectra of Liquid Water

X-ray absorption spectroscopy (XAS) is a powerful experimental technique to probe the local order in materials with core electron excitations. Experimental interpretation requires supporting theoretical calculations. For water, these calculations are very demanding and, to date, could only be done with major approximations that limited the accuracy of the calculated spectra. This prompted an intense debate on whether a substantial revision of the standard picture of tetrahedrally bonded water was necessary to improve the agreement of theory and experiment. Here, we report a new first-principles calculation of the XAS of water that avoids the approximations of prior work thanks to recent advances in electron excitation theory. The calculated XAS spectra, and their variation with changes of temperature and/or with isotope substitution, are in excellent quantitative agreement with experiments. The approach requires accurate quasi-particle wavefunctions beyond density functional theory approximations, accounts for the dynamics of quasi-particles and includes dynamic screening as well as renormalization effects due to the continuum of valence-level excitations. The three features observed in the experimental spectra are unambiguously attributed to excitonic effects. The pre-edge feature is associated to a bound intramolecular exciton, the main-edge feature is associated to an exciton localized within the coordination shell of the excited molecule, while the post-edge one is delocalized over more distant neighbors, as expected for a resonant state. The three features probe the local order at short, intermediate, and longer range relative to the excited molecule. The calculated spectra are fully consistent with a standard tetrahedral picture of water.

preprint2022arXiv

Nature of novel moiré exciton states in WSe$_2$/WS$_2$ heterobilayers

Moiré patterns of transition metal dichalcogenide (TMD) heterobilayers have proven to be an ideal platform to host unusual correlated electronic phases, emerging magnetism, and correlated exciton physics. While the existence of novel moiré excitonic states is established through optical measurements, the microscopic nature of these states is still poorly understood, often relying on empirically fit models. Here, combining large-scale first-principles GW-BSE calculations and micro-reflection spectroscopy, we identify the nature of the exciton resonances in WSe$_2$/WS$_2$ moiré superlattices, discovering a surprisingly rich set of moiré excitons that cannot be even qualitatively captured by prevailing continuum models. Our calculations reveal moiré excitons with distinct characters, including modulated Wannier excitons and previously unindentified intralayer charge-transfer excitons. Signatures of these distinct excitonic characters are confirmed experimentally via the unique carrier-density and magnetic-field dependences of different moiré exciton resonances. Our study highlights the highly non-trivial exciton states that can emerge in TMD moiré superlattices, and suggests novel ways of tuning many-body physics in moiré systems by engineering excited-states with specific spatial characters.

preprint2022arXiv

Optical and magneto-optical properties of ferromagnetic monolayer CrBr$_3$: A first-principles $GW$ and $GW$ plus Bethe-Salpeter equation study

The discovery of atomically thin two-dimensional (2D) magnetic semiconductors has triggered enormous research interest recently. In this work, we use first-principles many-body perturbation theory to study a prototypical 2D ferromagnetic semiconductor, monolayer chromium tribromide (CrBr$_3$). With broken time-reversal symmetry, spin-orbit coupling, and excitonic effects included through the full-spinor $GW$ and $GW$ plus Bethe-Salpeter equation ($GW$-BSE) methods, we compute the frequency-dependent layer polarizability tensor and dielectric function tensor that govern the optical and magneto-optical properties. In addition, we provide a detailed theoretical formalism for simulating magnetic circular dichroism, magneto-optical Kerr effect, and Faraday effect, demonstrating the approach with monolayer CrBr$_3$. Due to reduced dielectric screening in 2D and the localized nature of the Cr $d$ orbitals, we find strong self-energy effects on the quasiparticle band structure of monolayer CrBr$_3$ that give a 3.8 eV indirect band gap. Also, excitonic effects dominate the low-energy optical and magneto-optical responses in monolayer CrBr$_3$ where a large exciton binding energy of 2.3 eV is found for the lowest bright exciton state with excitation energy at 1.5 eV. We further find that the magneto-optical signals demonstrate strong dependence on the excitation frequency and substrate refractive index. Our theoretical framework for modeling optical and magneto-optical effects could serve as a powerful theoretical tool for future study of optoelectronic and spintronics devices consisting of van der Waals 2D magnets.

preprint2022arXiv

The electronic structure of $β$-HgS via $GW$ calculations

The electronic structure of the zincblende $β$-HgS is not well understood. Previous first-principles calculations using fully-relativistic density functional theory and many-body perturbation theory in the fully-relativistic $GW$ approach have predicted an inverted, topologically non-trivial ordering of these states, with the $s$-like $Γ_6$ state occupied. However, other calculations using the $GW$ approach in which spin-orbit coupling is added perturbatively ("$GW$+SOC") predict the $p$-$d$ hybridized $Γ_7$ and $Γ_8$ states to be occupied and the $Γ_6$ state to be unoccupied, suggesting that $β$-HgS is a topologically trivial small band gap semiconductor. In the present work, a plane-wave pseudopotential fully-relativistic $GW$ calculation finds a band ordering in agreement with the previous $GW$+SOC calculations. The calculated band gap is 0.10 eV and the electron effective mass is 0.07 $m_e$, in good agreement with experiment.

preprint2022arXiv

Theory of Exciton-Phonon Coupling

The effects of the electron-phonon interaction on optical excitations can be understood in terms of exciton-phonon coupling, and require a careful treatment in low-dimensional materials with strongly bound excitons or strong electron-hole interaction in general. Through phonon absorption and emission processes, the optically accessible excitons are scattered into otherwise optically dark finite-momentum exciton states. We derive a practical expression for the phonon-induced term of the exciton self-energy (denoted as the exciton-phonon self-energy) that gives the temperature dependence of the optical transition energies and their lifetime broadening resulting from the exciton's interaction with the phonons. We illustrate this theory on a two-dimensional model, and show that our expression for the exciton-phonon self-energy differs qualitatively from previous expressions found in the literature that neglect the exciton binding or electron-hole correlations.

preprint2022arXiv

Unconventional Excitonic States with Phonon Sidebands in Layered Silicon Diphosphide

Many-body interactions between quasiparticles (electrons, excitons, and phonons) have led to the emergence of new complex correlated states and are at the core of condensed matter physics and material science. In low-dimensional materials, unique electronic properties for these correlated states could significantly affect their optical properties. Herein, combining photoluminescence, optical reflection measurements and theoretical calculations, we demonstrate an unconventional excitonic state and its bound phonon sideband in layered silicon diphosphide (SiP$_2$), in which the bound electron-hole pair is composed of electrons confined within one-dimensional phosphorus$-$phosphorus chains and holes extended in two-dimensional SiP$_2$ layers. The excitonic state and the emergent phonon sideband show linear dichroism and large energy redshifts with increasing temperature. Within the $GW$ plus Bethe$-$Salpeter equation calculations and solving the generalized Holstein model non-perturbatively, we confirm that the observed sideband feature results from the correlated interaction between excitons and optical phonons. Such a layered material provides a new platform to study excitonic physics and many-particle effects.

preprint2022arXiv

Using dynamic mode decomposition to predict the dynamics of a two-time non-equilibrium Green's function

Computing the numerical solution of the Kadanoff-Baym equations, a set of nonlinear integral differential equations satisfied by two-time Green's functions derived from many-body perturbation theory for a quantum many-body system away from equilibrium, is a challenging task. Recently, we have successfully applied dynamic mode decomposition (DMD) to construct a data driven reduced order model that can be used to extrapolate the time-diagonal of a two-time Green's function from numerical solution of the KBE within a small time window. In this paper, we extend the previous work and use DMD to predict off-diagonal elements of the two-time Green's function. We partition the two-time Green's function into a number of one-time functions along the diagonal and subdiagonls of the two-time window as well as in horizontal and vertical directions. We use DMD to construct separate reduced order models to predict the dynamics of these one-time functions in a two-step procedure. We extrapolate along diagonal and several subdiagonals within a subdiagonal band of a two-time window in the first step. In the second step, we use DMD to extrapolate the Green's function outside of the sub-diagonal band. We demonstrate the efficiency and accuracy of this approach by applying it to a two-band Hubbard model problem.

preprint2021arXiv

Exciton Shift Currents: DC Conduction with Sub-bandgap Photo Excitations

Shift current is a DC current generated from nonlinear light-matter interaction in a non-centrosymmetric crystal and is considered a promising candidate for next generation photovoltaic devices. The mechanism for shift currents in real materials is, however, still not well understood, especially if electron-hole interactions are taken into account. Here, we employ a first-principles interacting Green's-function approach on the Keldysh contour to study photocurrents generated by nonlinear optical processes in real materials and discover a strong DC shift current at subbandgap excitation frequencies in monolayer GeS due to strongly bound excitons, as well as giant enhancement in the shift current coefficients at above bandgap photon frequencies. Our results suggest that atomically thin two-dimensional materials may be promising building blocks for next generation shift current devices with efficiencies beyond the Shockley-Queisser limit.

preprint2021arXiv

Solving the Bethe-Salpeter Equation on a Subspace: Approximations and Consequences for Low-dimensional Materials

It is well known that the ambient environment can dramatically renormalize the quasiparticle gap and exciton binding energies in low-dimensional materials, but the effect of the environment on the energy splitting of the spin-singlet and spin-triplet exciton states is less understood. A prominent effect is the renormalization of the exciton binding energy and optical strength (and hence the optical spectrum) through additional screening of the direct Coulomb term describing the attractive electron-hole interaction in the kernel of the Bethe-Salpeter equation (BSE). The repulsive exchange interaction responsible for the singlet-triplet slitting, on the other hand, is unscreened within formal many-body perturbation theory. However, Loren Benedict argued that in practical calculations restricted to a subspace of the full Hilbert space, the exchange interaction should be appropriately screened by states outside of the subspace, the so-called $S$ approximation \cite{Benedict2002}. Here, we systematically explore the accuracy of the $S$ approximation for different confined systems, including a molecule and heterostructures of semiconducting and metallic layered materials. We show that the $S$ approximation is actually exact in the limit of small exciton binding energies (i.e., small direct term) and can be used to significantly accelerate convergence of the exciton energies with respect to the number of empty states, provided that a particular effective screening consistent with the conventional Tamm-Dancoff approximation is employed. We further find that the singlet-triplet splitting in the energy of the excitons is largely unaffected by the external dielectric environment for most quasi-two-dimensional materials.

preprint2021arXiv

Spin Splitting of Dopant Edge States in Magnetic Zigzag Graphene Nanoribbons

Spin-ordered electronic states in hydrogen-terminated zigzag nanographene give rise to magnetic quantum phenomena that have sparked renewed interest in carbon-based spintronics. Zigzag graphene nanoribbons (ZGNRs), quasi one-dimensional semiconducting strips of graphene featuring two parallel zigzag edges along the main axis of the ribbon, are predicted to host intrinsic electronic edge states that are ferromagnetically ordered along the edges of the ribbon and antiferromagnetically coupled across its width. Despite recent advances in the bottom-up synthesis of atomically-precise ZGNRs, their unique electronic structure has thus far been obscured from direct observations by the innate chemical reactivity of spin-ordered edge states. Here we present a general technique for passivating the chemically highly reactive spin-polarized edge states by introducing a superlattice of substitutional nitrogen-dopants along the edges of a ZGNR. First-principles GW calculations and scanning tunneling spectroscopy reveal a giant spin splitting of the low-lying nitrogen lone-pair flat bands by a large exchange field (~850 Tesla) induced by the spin-polarized ferromagnetically ordered edges of ZGNRs. Our findings directly corroborate the nature of the predicted emergent magnetic order in ZGNRs and provide a robust platform for their exploration and functional integration into nanoscale sensing and logic devices.

preprint2020arXiv

Band gap renormalization, carrier mobilities, and the electron-phonon self-energy in crystalline naphthalene

Organic molecular crystals are expected to feature appreciable electron-phonon interactions that influence their electronic properties at zero and finite temperature. In this work, we report first-principles calculations and an analysis of the electron-phonon self-energy in naphthalene crystals. We compute the zero-point renormalization and temperature dependence of the fundamental band gap, and the resulting scattering lifetimes of electronic states near the valence- and conduction-band edges employing density functional theory. Further, our calculated phonon renormalization of the $GW$-corrected quasiparticle band structure predicts a fundamental band gap of 5 eV for naphthalene at room temperature, in good agreement with experiments. From our calculated phonon-induced electron lifetimes, we obtain the temperature-dependent mobilities of electrons and holes in good agreement with experimental measurements at room temperatures. Finally, we show that an approximate energy self-consistent computational scheme for the electron-phonon self-energy leads to the prediction of strong satellite bands in the electronic band structure. We find that a single calculation of the self-energy can reproduce the self-consistent results of the band gap renormalization and electrical mobilities for naphthalene, provided that the on-the-mass-shell approximation is used, i.e., if the self-energy is evaluated at the bare eigenvalues.

preprint2020arXiv

Comparison of GW band structure to semi-empirical approach for an FeSe monolayer

We present the G$_0$W$_0$ band structure, core levels, and deformation potential of monolayer FeSe in the paramagnetic phase based on a starting mean field of the Kohn Sham density functional theory (DFT) with the PBE functional. We find the GW correction increases the bandwidth of the states forming the $M$ pocket near the Fermi energy, while leaving the $Γ$ pocket roughly unchanged. We then compare the G$_0$W$_0$ quasiparticle band energies with the band structure from a simple empirical +A approach, which was recently proposed to capture the renormalization of the electron-phonon interaction going beyond DFT in FeSe, when used as a starting point in density functional perturbation theory (DFPT). We show that this empirical correction succeeds in approximating the GW non-local and dynamical self energy in monolayer FeSe and reproduces the GW band structure near the Fermi surface, the core energy levels, and the deformation potential (electron-phonon coupling).

preprint2020arXiv

Identifying substitutional oxygen as a prolific point defect in monolayer transition metal dichalcogenides with experiment and theory

Chalcogen vacancies are considered to be the most abundant point defects in two-dimensional (2D) transition-metal dichalcogenide (TMD) semiconductors, and predicted to result in deep in-gap states (IGS). As a result, important features in the optical response of 2D-TMDs have typically been attributed to chalcogen vacancies, with indirect support from Transmission Electron Microscopy (TEM) and Scanning Tunneling Microscopy (STM) images. However, TEM imaging measurements do not provide direct access to the electronic structure of individual defects; and while Scanning Tunneling Spectroscopy (STS) is a direct probe of local electronic structure, the interpretation of the chemical nature of atomically-resolved STM images of point defects in 2D-TMDs can be ambiguous. As a result, the assignment of point defects as vacancies or substitutional atoms of different kinds in 2D-TMDs, and their influence on their electronic properties, has been inconsistent and lacks consensus. Here, we combine low-temperature non-contact atomic force microscopy (nc-AFM), STS, and state-of-the-art ab initio density functional theory (DFT) and GW calculations to determine both the structure and electronic properties of the most abundant individual chalcogen-site defects common to 2D-TMDs. Surprisingly, we observe no IGS for any of the chalcogen defects probed. Our results and analysis strongly suggest that the common chalcogen defects in our 2D-TMDs, prepared and measured in standard environments, are substitutional oxygen rather than vacancies.

preprint2020arXiv

The bright side of defects in MoS$_2$ and WS$_2$ and a generalizable chemical treatment protocol for defect passivation

Structural defects are widely regarded as detrimental to the optoelectronic properties of monolayer transition metal dichalcogenides, leading to concerted efforts to eliminate defects via improved materials growth or post-growth passivation. Here, using steady-state and ultrafast optical spectroscopy, supported by ab initio calculations, we demonstrate that sulfur vacancy defects act as exciton traps. Current chemical treatments do not passivate these sites, leading to decreased mobility and trap-limited photoluminescence. We present a generalizable treatment protocol based on the use of passivating agents such as thiols or sulfides in combination with a Lewis acid to passivate sulfur vacancies in monolayer MoS$_2$ and WS$_2$, increasing photoluminescence up to 275 fold, while maintaining mobilities. Our findings suggest a route for simple and rational defect engineering strategies, where the passivating agent varies the electronic properties, thereby allowing the design of new heterostructures.

preprint2020arXiv

Topology Classification from Chiral Symmetry: Chiral Phase Index and Spin Correlations in Graphene Nanoribbons

Topology concepts have significantly deepened of our understanding in recent years of the electronic properties of one-dimensional (1D) nano structures such as the graphene nanoribbons. Controlling topological electronic properties of GNRs has been demonstrated in both theoretical studies and experimental realization. Most previous works rely on classification theory requiring both time reversal and spatial symmetry of a unit cell in the 1D bulk material that is commensurate to its boundary. To access boundary structures that lead to unit cell with no spatial symmetry and to generalize the theory, we propose here another classification scheme, using chiral symmetry, to arrive at a Z classification that is not only applicable to GNRs with arbitrary terminations, but also to any general 1D chiral structures. This theory, combining with Lieb's theorem, moreover enables access to the electron's spin degree of freedom, allowing for investigation of spin physics.

preprint2019arXiv

Reproducibility in $G_0W_0$ Calculations for Solids

Ab initio many-body perturbation theory within the $GW$ approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point'', are used to construct $G$ and $W$, and then perturbatively corrected by the resultant $GW$ self-energy. In practice, there are several ways to construct the $GW$ self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO$_2$, reported $GW$ fundamental gaps can vary by more than 1 eV. In this work, we address the convergence and key approximations in contemporary $G_0W_0$ calculations, including frequency-integration schemes and the treatment of the Coulomb divergence in the exact-exchange term. We study several systems,and compare three different $GW$ codes: BerkeleyGW, Abinit and Yambo. We demonstrate, for the first time, that the same quasiparticle energies for systems in the condensed phase can be obtained with different codes, and we provide a comprehensive assessment of implementations of the $GW$ approximation.

preprint2019arXiv

Visualizing Exotic Orbital Texture in the Single-Layer Mott Insulator 1T-TaSe2

Mott insulating behavior is induced by strong electron correlation and can lead to exotic states of matter such as unconventional superconductivity and quantum spin liquids. Recent advances in van der Waals material synthesis enable the exploration of novel Mott systems in the two-dimensional limit. Here we report characterization of the local electronic properties of single- and few-layer 1T-TaSe2 via spatial- and momentum-resolved spectroscopy involving scanning tunneling microscopy and angle-resolved photoemission. Our combined experimental and theoretical study indicates that electron correlation induces a robust Mott insulator state in single-layer 1T-TaSe2 that is accompanied by novel orbital texture. Inclusion of interlayer coupling weakens the insulating phase in 1T-TaSe2, as seen by strong reduction of its energy gap and quenching of its correlation-driven orbital texture in bilayer and trilayer 1T-TaSe2. Our results establish single-layer 1T-TaSe2 as a useful new platform for investigating strong correlation physics in two dimensions.

preprint2016arXiv

Accurate determination of the valence band edge in hard x-ray photoemission spectra using GW theory

We introduce a new method for determining accurate values of the valence-band maximum in x-ray photoemission spectra. Specifically, we align the sharpest peak in the valence-band region of the experimental spectrum with the corresponding feature of a theoretical valence-band density of states curve from ab initio GW theory calculations. This method is particularly useful for soft and hard x-ray photoemission studies of materials with a mixture of valence-band characters, where strong matrix element effects can render standard methods for extracting the valence-band maximum unreliable. We apply our method to hydrogen-terminated boron-doped diamond, which is a promising substrate material for novel solar cell devices. By carrying out photoemission experiments with variable light polarizations, we verify the accuracy of our analysis and the general validity of the method.

preprint2016arXiv

Antiferromagnetism enables electron-phonon coupling in iron-based superconductors

It is shown that a generic form of an antiferromagnetic wave function opens strong electron-phonon coupling channels in the iron-based superconductors. In the nonmagnetic state these channels exist locally on a single iron atom, but are canceled out between the two iron atoms in the primitive unit cell. Our findings are based on symmetry and the presence of an xz/yz Fermi surface near the M point and thus should be relevant for the known iron-based pnictide or chalcogenide superconductors.

preprint2016arXiv

Dispersion and lineshape of plasmon satellites in one, two and three dimensions

Using state-of-the-art many-body Green's function calculations based on the "GW plus cumulant" approach, we analyze the properties of plasmon satellites in the electron spectral function resulting from electron-plasmon interactions in one-, two- and three-dimensional systems. Specifically, we show how their dispersion relation, lineshape and linewidth are related to the properties of the constituent electrons and plasmons. To gain insight into the many-body processes giving rise to the formation of plasmon satellites, we connect the "GW plus cumulant" approach to a many-body wavefunction picture of electron-plasmon interactions and introduce the coupling-strength weighted electron-plasmon joint-density states as a powerful concept for understanding plasmon satellites.

preprint2016arXiv

Effective Mass in Bilayer Graphene at Low Carrier Densities: the Role of Potential Disorder and Electron-Electron Interaction

In a two-dimensional electron gas, the electron-electron interaction generally becomes stronger at lower carrier densities and renormalizes the Fermi liquid parameters such as the effective mass of carriers. We combine experiment and theory to study the effective masses of electrons and holes $m^*_e$ and $m^*_h$ in bilayer graphene in the low carrier density regime of order 1 * 10^11 cm^-2. Measurements use temperature-dependent low-field Shubnikov-de Haas (SdH) oscillations are observed in high-mobility hexagonal boron nitride (h-BN) supported samples. We find that while $m^*_e$ follows a tight-binding description in the whole density range, $m^*_h$ starts to drop rapidly below the tight-binding description at carrier density n = 6 * 10^11 cm^-2 and exhibits a strong suppression of 30% when n reaches 2 * 10^11 cm^-2. Contributions from electron-electron interaction alone, evaluated using several different approximations, cannot explain the experimental trend. Instead, the effect of potential fluctuation and the resulting electron-hole puddles play a crucial role. Calculations including both the electron-electron interaction and disorder effects explain the experimental data qualitatively and quantitatively. This study reveals an unusual disorder effect unique to two-dimensional semi-metallic systems.

preprint2016arXiv

Klein tunneling and supercollimation of pseudospin-1 electromagnetic waves

Pseudospin plays a central role in many novel physical properties of graphene and other artificial systems which have pseudospins of 1/2. Here we show that in certain photonic crystals (PCs) exhibiting conical dispersions at k = 0, the eigenmodes near the "Dirac-like point" can be described by an effective spin-orbit Hamiltonian with a pseudospin of 1, treating wave propagations in the upper cone, the lower cone and a flat band (corresponding to zero refractive index) within a unified framework. The 3-component spinor gives rise to boundary conditions distinct from those of pseudospin-1/2, leading to new wave transport behaviors as manifested in Klein tunneling and supercollimation. For example, collimation can be realized more easily with pseudospin-1 than pseudospin-1/2. The special wave scattering properties of pseudospin-1 photons, coupled with the discovery that the effective photonic "potential" can be varied by a simple change of length scale, may offer new ways to control photon transport.

preprint2016arXiv

Proposal for a bulk material based on a monolayer FeSe on SrTiO3 high-temperature superconductor

Motivated by the high superconducting transition temperature of monolayer FeSe on SrTiO3, we propose a potential three-dimensional high-temperature superconductor superlattice FeSe-SrTiO3 and study its structural stability and electronic structure using density functional theory. We find that the binding energy between the FeSe and SrTiO3 layers is about 0.7 eV per (Fe2Se2) unit and that it saturates already within a single TiO2 atomic layer of SrTiO3. In addition we analyzed the dynamical stability of the superlattice and compared it to the case of bulk SrTiO3.

preprint2016arXiv

Screening and Many-Body Effects in Two-Dimensional Crystals: Monolayer MoS$_{2}$

We present a systematic study of the variables affecting the electronic and optical properties of two-dimensional(2D) crystals within \textit{ab initio} GW and GW plus Bethe Salpeter Equation (GW-BSE) calculations. As a prototypical 2D transition metal dichalcogenide material, we focus our study on monolayer MoS${}_2$. We find that the reported variations in GW-BSE results in the literature for monolayer MoS${}_2$ and related systems arise from different treatments of the long-range Coulomb interaction in supercell calculations and convergence of k-grid sampling and cutoffs for various quantities such as the dielectric screening. In particular, the quasi-2D nature of the system gives rise to fast spatial variations in the screening environment, which are computationally challenging to resolve. We also show that common numerical treatments to remove the divergence in the Coulomb interaction can shift the exciton continuum leading to false convergence with respect to k-point sampling. Our findings apply to GW-BSE calculations on any low-dimensional semiconductors.

preprint2016arXiv

Temperature-induced topological phase transitions: promoted vs. suppressed non-trivial topology

We determine the topological phase diagram of BiTl(S$_{1-δ}$Se$_δ$)$_2$ as a function of doping and temperature from first-principles calculations. Due to electron\textendash phonon interaction, the bands are renormalized at finite temperature, allowing for a transition between the trivial ($Z_2=0$) and non-trivial ($Z_2=1$) topological phase. We find two distinct regions of the phase diagram with non-trivial topology. In BiTlS$_2$, the phonons promote the crystal to the topological phase at high temperature, while in BiTlSe$_2$, the topological phase exists only at low temperature. This behaviour is explained by the symmetry of the phonon coupling potential, whereby the even phonon modes (whose potential is even under inversion) promote the topological phase and the odd phonon modes promote the trivial phase.

preprint2016arXiv

Why do Collapsed Carbon Nanotubes Twist?

We study the collapsing and subsequent spontaneous twisting of a carbon nanotube by in-situ transmission electron microscopy. A custom-sized nanotube is first created in the microscope by selectively extracting shells from a parent multi-wall tube. The few-wall, large-diameter daughter nanotube is driven to collapse via mechanical stimulation, after which the ribbon-like collapsed tube spontaneously twists along its long axis. In-situ diffraction experiments fully characterize the uncollapsed and collapsed tubes. From the experimental observations and associated theoretical analysis, the origin of the twisting is determined to be compressive strain due to charge imbalance.

preprint2015arXiv

Automated construction of maximally localized Wannier functions: Optimized projection functions method

Maximally localized Wannier functions are widely used in electronic structure theory for analyses of bonding, electric polarization, orbital magnetization, and for interpolation. The state of the art method for their construction is based on the method of Marzari and Vanderbilt. One of the practical difficulties of this method is guessing functions (initial projections) that approximate the final Wannier functions. Here we present an approach based on optimized projection functions that can construct maximally localized Wannier functions without a guess. We describe and demonstrate this approach on several realistic examples.

preprint2015arXiv

First-principles theory of electron-spin fluctuation coupling and superconducting instabilities in iron selenide

We present first-principles calculations of the coupling of quasiparticles to spin fluctuations in iron selenide and discuss which types of superconducting instabilities this coupling gives rise to. We find that strong antiferromagnetic stripe-phase spin fluctuations lead to large coupling constants for superconducting gaps with $s_\pm$-symmetry, but these coupling constants are significantly reduced by other spin fluctuations with small wave vectors. An accurate description of this competition and an inclusion of band structure and Stoner parameter renormalization effects lead to a value of the coupling constant for an $s_\pm$-symmetric gap which can produce a superconducting transition temperature consistent with experimental measurements.

preprint2015arXiv

Large electron-phonon interactions from FeSe phonons in a monolayer

We show that electron-phonon coupling can induce strong electron pairing in an FeSe monolayer on a SrTiO3 substrate (experimental indications for superconducting Tc are between 65 and 109K). The role of the SrTiO3 substrate in increasing the coupling is two-fold. First, the interaction of the FeSe and TiO2 terminated face of SrTiO3 prevents the FeSe monolayer from undergoing a shear-type (orthorhombic, nematic) structural phase transition. Second, the substrate allows an anti-ferromagnetic ground state of FeSe which opens electron-phonon coupling channels within the monolayer that are prevented by symmetry in the non-magnetic phase. The spectral function for the electron-phonon coupling (alpha2F) in our calculations agrees well with inelastic tunneling data.

preprint2015arXiv

Nonanalyticity, Valley Quantum Phases, and Light-like Exciton Dispersion in Monolayer TMDs: Theory and First-Principles Calculations

Exciton dispersion as a function of center-of-mass momentum, \textbf{Q}, is essential to the understanding of exciton dynamics. We use the ab-initio GW-Bethe Salpeter equation method to calculate the dispersion of excitons in monolayer MoS$_2$ and find a nonanalytic light-like dispersion. This behavior arises from an unusual $|Q|$ term in both the intra- and inter-valley exchange of the electron-hole interaction, which concurrently gives rise to a valley quantum phase of winding number two. A simple effective Hamiltonian to $Q^2$ order with analytic solutions is derived to describe quantitatively these behaviors.

preprint2015arXiv

Observing Atomic Collapse Resonances in Artificial Nuclei on Graphene

Relativistic quantum mechanics predicts that when the charge of a superheavy atomic nucleus surpasses a certain threshold, the resulting strong Coulomb field causes an unusual atomic collapse state; this state exhibits an electron wave function component that falls toward the nucleus, as well as a positron component that escapes to infinity. In graphene, where charge carriers behave as massless relativistic particles, it has been predicted that highly charged impurities should exhibit resonances corresponding to these atomic collapse states. We have observed the formation of such resonances around artificial nuclei (clusters of charged calcium dimers) fabricated on gated graphene devices via atomic manipulation with a scanning tunneling microscope. The energy and spatial dependence of the atomic collapse state measured with scanning tunneling microscopy revealed unexpected behavior when occupied by electrons.

preprint2015arXiv

Probing the Role of Interlayer Coupling and Coulomb Interactions on Electronic Structure in Few-Layer MoSe2 Nanostructures

Despite the weak nature of interlayer forces in transition metal dichalcogenide (TMD) materials, their properties are highly dependent on the number of layers in the few-layer two-dimensional (2D) limit. Here, we present a combined scanning tunneling microscopy/spectroscopy and GW theoretical study of the electronic structure of high quality single- and few-layer MoSe2 grown on bilayer graphene. We find that the electronic (quasiparticle) bandgap, a fundamental parameter for transport and optical phenomena, decreases by nearly one electronvolt when going from one layer to three due to interlayer coupling and screening effects. Our results paint a clear picture of the evolution of the electronic wave function hybridization in the valleys of both the valence and conduction bands as the number of layers is changed. This demonstrates the importance of layer number and electron-electron interactions on van der Waals heterostructures, and helps to clarify how their electronic properties might be tuned in future 2D nanodevices.

preprint2015arXiv

Spectral Functions of the Uniform Electron Gas via Coupled-Cluster Theory and Comparison to the $GW$ and Related Approximations

We use, for the first time, ab initio coupled-cluster theory to compute the spectral function of the uniform electron gas at a Wigner-Seitz radius of $r_\mathrm{s}=4$. The coupled-cluster approximations we employ go significantly beyond the diagrammatic content of state-of-the-art $GW$ theory. We compare our calculations extensively to $GW$ and $GW$-plus-cumulant theory, illustrating the strengths and weaknesses of these methods in capturing the quasiparticle and satellite features of the electron gas. Our accurate calculations further allow us to address the long-standing debate over the occupied bandwidth of metallic sodium. Our findings indicate that the future application of coupled-cluster theory to condensed phase material spectra is highly promising.

preprint2015arXiv

Structure Preserving Parallel Algorithms for Solving the Bethe-Salpeter Eigenvalue Problem

The Bethe-Salpeter eigenvalue problem is a dense structured eigenvalue problem arising from discretized Bethe-Salpeter equation in the context of computing exciton energies and states. A computational challenge is that at least half of the eigenvalues and the associated eigenvectors are desired in practice. We establish the equivalence between Bethe-Salpeter eigenvalue problems and real Hamiltonian eigenvalue problems. Based on theoretical analysis, structure preserving algorithms for a class of Bethe-Salpeter eigenvalue problems are proposed. We also show that for this class of problems all eigenvalues obtained from the Tamm-Dancoff approximation are overestimated. In order to solve large scale problems of practical interest, we discuss parallel implementations of our algorithms targeting distributed memory systems. Several numerical examples are presented to demonstrate the efficiency and accuracy of our algorithms.

preprint2015arXiv

SU(4) symmetry breaking revealed by magneto-optical spectroscopy in epitaxial graphene

Refined infrared magnetotransmission experiments have been performed in magnetic fields B up to 35 T on a series of multilayer epitaxial graphene samples. Following the main optical transition involving the n=0 Landau level (LL), we observe a new absorption transition increasing in intensity with magnetic fields B>26 T. Our analysis shows that this is a signature of the breaking of the SU(4) symmetry of the n=0 LL. Using a quantitative model, we show that the only symmetry-breaking scheme consistent with our experiments is a charge density wave (CDW).

preprint2014arXiv

Ab initio quasiparticle bandstructure of ABA and ABC-stacked graphene trilayers

We obtain the quasiparticle band structure of ABA and ABC-stacked graphene trilayers through ab initio density functional theory (DFT) and many-body quasiparticle calculations within the GW approximation. To interpret our results, we fit the DFT and GW $π$ bands to a low energy tight-binding model, which is found to reproduce very well the observed features near the K point. The values of the extracted hopping parameters are reported and compared with available theoretical and experimental data. For both stackings, the self energy corrections lead to a renormalization of the Fermi velocity, an effect also observed in previous calculations on monolayer graphene. They also increase the separation between the higher energy bands, which is proportional to the nearest neighbor interlayer hopping parameter $γ_1$. Both features are brought to closer agreement with experiment through the self energy corrections. Finally, other effects, such as trigonal warping, electron-hole asymmetry and energy gaps are discussed in terms of the associated parameters.

preprint2014arXiv

Effect of spin fluctuations on quasiparticle excitations: first-principles theory and application to sodium and lithium

We present first-principles calculations for quasiparticle excitations in sodium and lithium including the effects of charge and spin fluctuations. We employ the Overhauser-Kukkonen form for the electron self energy arising from spin fluctuations and demonstrate that the coupling of electrons to spin fluctuations gives an important contribution to the quasiparticle lifetime, but does not significantly reduce the occupied bandwidth. Including correlation effects beyond the random-phase approximation in the screening from charge fluctuations yields good agreement with experiment.

preprint2014arXiv

Effects of self-consistency and plasmon-pole models on GW calculations for closed-shell molecules

We present theoretical calculations of quasiparticle energies in closed-shell molecules using the GW method. We compare three different approaches: a full-frequency $G_0W_0$ (FF-$G_0W_0$) method with density functional theory (DFT-PBE) used as a starting mean field; a full-frequency $GW_0$ (FF-$GW_0$) method where the interacting Green's function is approximated by replacing the DFT energies with self-consistent quasiparticle energies or Hartree-Fock energies; and a $G_0W_0$ method with a Hybertsen-Louie generalized plasmon-pole model (HL GPP-$G_0W_0$). While the latter two methods lead to good agreement with experimental ionization potentials and electron affinities for methane, ozone, and beryllium oxide molecules, FF-$G_0W_0$ results can differ by more than one electron volt from experiment. We trace this failure of the FF-$G_0W_0$ method to the occurrence of incorrect self-energy poles describing shake-up processes in the vicinity of the quasiparticle energies.

preprint2014arXiv

Electron supercollimation in graphene and Dirac fermion materials using one-dimensional disorder potentials

Electron supercollimation, in which a wavepacket is guided to move undistorted along a selected direction, is a highly desirable property that has yet been realized experimentally. Disorder in general is expected to inhibit supercollimation. Here, we report a counter-intuitive phenomenon of electron supercollimation by disorder in graphene and related Dirac fermion materials. We show that one can use one-dimensional disorder potentials to control electron wavepacket transport. This is distinct from known systems where an electron wavepacket would be further spread by disorder and hindered in the potential fluctuating direction. The predicted phenomenon has significant implications in the understanding and applications of electron transport in Dirac fermion materials.

preprint2014arXiv

First-principles DFT+GW study of oxygen vacancies in rutile TiO2

We perform first-principles calculations of the quasiparticle defect states, charge transition levels, and formation energies of oxygen vacancies in rutile titanium dioxide. The calculations are done within the recently developed combined DFT+GW formalism, including the necessary electrostatic corrections for the supercells with charged defects. We find the oxygen vacancy to be a negative U defect, where U is the defect electron addition energy. For the values of Fermi level below 2.8 eV (relative to the valence band maximum) we find the +2 charge state of the vacancy to be the most stable, while above 2.8 eV we find that the neutral charge state is the most stable.

preprint2014arXiv

Numerical integration for ab initio many-electron self energy calculations within the GW approximation

We present a numerical integration scheme for evaluating the convolution of a Green's function with a screened Coulomb potential on the real axis in the GW approximation of the self energy. Our scheme takes the zero broadening limit in Green's function first, replaces the numerator of the integrand with a piecewise polynomial approximation, and performs principal value integration on subintervals analytically. We give the error bound of our numerical integration scheme and show by numerical examples that it is more reliable and accurate than the standard quadrature rules such as the composite trapezoidal rule. We also discuss the benefit of using different self energy expressions to perform the numerical convolution at different frequencies.

preprint2014arXiv

Probing Excitonic Dark States in Single-layer Tungsten Disulfide

Transition metal dichalcogenide (TMDC) monolayer has recently emerged as an important two-dimensional semiconductor with promising potentials for electronic and optoelectronic devices. Unlike semi-metallic graphene, layered TMDC has a sizable band gap. More interestingly, when thinned down to a monolayer, TMDC transforms from an indirect bandgap to a direct bandgap semiconductor, exhibiting a number of intriguing optical phenomena such as valley selective circular dichroism, doping dependent charged excitons, and strong photocurrent responses. However, the fundamental mechanism underlying such a strong light-matter interaction is still under intensive investigation. The observed optical resonance was initially considered to be band-to-band transitions. In contrast, first-principle calculations predicted a much larger quasiparticle band gap size and an optical response that is dominated by excitonic effects. Here, we report experimental evidence of the exciton dominance mechanism by discovering a series of excitonic dark states in single-layer WS2 using two-photon excitation spectroscopy. In combination with GW-BSE theory, we find the excitons are Wannier excitons in nature but possess extraordinarily large binding energy (~0.7 eV), leading to a quasiparticle band gap of 2.7 eV. These strongly bound exciton states are observed stable even at room temperature. We reveal an exciton series in significant deviation from hydrogen models, with a novel inverse energy dependence on the orbital angular momentum. These excitonic energy levels are experimentally found robust against environmental perturbations. The discovery of excitonic dark states and exceptionally large binding energy not only sheds light on the importance of many-electron effects in this two-dimensional gapped system, but also holds exciting potentials for the device application of TMDC monolayers and their heterostructures.

preprint2014arXiv

Satellite structures in the spectral functions of semiconductor quantum well two-dimensional electron gases: a GW plus cumulant study

We present theoretical calculations for the spectral functions and single-particle densities of states of the two-dimensional electron gas in semiconductor quantum wells at different electron densities using the GW+cumulant method. We compare our results to GW only calculations and find significant differences in the description of the satellites between the two theories: while GW theory predicts the existence of a plasmaron excitation, no such excitation is found in GW+cumulant theory. We compare our results to experimental tunneling spectra from semiconductor quantum wells and find good agreement for the satellite properties.

preprint2014arXiv

Tunable Magnetism and Half-Metallicity in Hole-doped Monolayer GaSe

We find, through first-principles calculations, that hole doping induces a ferromagnetic phase transition in monolayer GaSe. Upon increasing hole density, the average spin magnetic moment per carrier increases and reaches a plateau near 1.0 $μ_{\rm{B}}$/carrier in a range of $3\times 10^{13}$/cm$^{2}$-$1\times 10^{14}$/cm$^{2}$ with the system in a half-metal state before the moment starts to descend abruptly. The predicted magnetism originates from an exchange splitting of electronic states at the top of the valence band where the density of states exhibits a sharp van Hove singularity in this quasi-two-dimensional system.

preprint2013arXiv

Atomically perfect torn graphene edges and their reversible reconstruction

The atomic structure of graphene edges is critical in determining the electrical, magnetic, and chemical properties of truncated graphene structures, notably nanoribbons. Unfortunately, graphene edges are typically far from ideal and suffer from atomic-scale defects, structural distortion, and unintended chemical functionalization, leading to unpredictable properties. Here we report that graphene edges fabricated by electron-beam-initiated mechanical rupture or tearing in high vacuum are clean and largely atomically perfect, oriented in either the armchair or zigzag direction. Via aberration-corrected transmission electron microscopy, we demonstrate reversible and extended pentagon-heptagon (5-7) reconstruction at zigzag edges, and explore experimentally and theoretically the dynamics of the transitions between configuration states. Good theoretical-experimental agreement is found for the flipping rates between 5-7 and 6-6 zigzag edge states. Our study demonstrates that simple ripping is remarkably effective in producing atomically clean, ideal terminations thus providing a valuable tool for realizing atomically-tailored graphene and facilitating meaningful experimental study.

preprint2013arXiv

Experimentally Engineering the Edge Termination of Graphene Nanoribbons

The edges of graphene nanoribbons (GNRs) have attracted much interest due to their potentially strong influence on GNR electronic and magnetic properties. Here we report the ability to engineer the microscopic edge termination of high quality GNRs via hydrogen plasma etching. Using a combination of high-resolution scanning tunneling microscopy and first-principles calculations, we have determined the exact atomic structure of plasma-etched GNR edges and established the chemical nature of terminating functional groups for zigzag, armchair and chiral edge orientations. We find that the edges of hydrogen-plasma-etched GNRs are generally flat, free of structural reconstructions and are terminated by hydrogen atoms with no rehybridization of the outermost carbon edge atoms. Both zigzag and chiral edges show the presence of edge states.

preprint2013arXiv

Imprint of transition metal d-orbitals on graphene Dirac cone

We investigate the influence of SiO2, Au, Ag, Cu, and Pt substrates on the Raman spectrum of graphene. Experiments reveal particularly strong modifications to the intensity, position, width, and shape of the Raman signal of graphene on platinum, compared to that of suspended graphene. The modifications strongly depend on the relative orientation of the graphene and platinum lattices. These observations are theoretically investigated and shown to originate from hybridization of electronic states in graphene and d-orbitals in platinum. It is expected that, quite generally, hybridization between graphene and any material with d-orbitals near the Fermi level will result in an imprint on the graphene Dirac cone which depends sensitively on the relative orientation of the respective lattices.

preprint2013arXiv

Intermolecular interactions and substrate effects for an adamantane monolayer on the Au(111) surface

We study theoretically and experimentally the infrared (IR) spectrum of an adamantane monolayer on a Au(111) surface. Using a new STM-based IR spectroscopy technique (IRSTM) we are able to measure both the nanoscale structure of an adamantane monolayer on Au(111) as well as its infrared spectrum, while DFT-based ab initio calculations allow us to interpret the microscopic vibrational dynamics revealed by our measurements. We find that the IR spectrum of an adamantane monolayer on Au(111) is substantially modified with respect to the gas-phase IR spectrum. The first modification is caused by the adamantane--adamantane interaction due to monolayer packing and it reduces the IR intensity of the 2912 cm$^{-1}$ peak (gas phase) by a factor of 3.5. The second modification originates from the adamantane--gold interaction and it increases the IR intensity of the 2938 cm$^{-1}$ peak (gas phase) by a factor of 2.6, and reduces its frequency by 276 cm$^{-1}$. We expect that the techniques described here can be used for an independent estimate of substrate effects and intermolecular interactions in other diamondoid molecules, and for other metallic substrates.

preprint2013arXiv

Optical Spectrum of MoS${}_2$: Many-body Effects and Diversity of Exciton States

We present first-principles calculations of the optical response of monolayer molybdenum disulfide employing the GW-Bethe Salpeter equation (GW-BSE) approach including self-energy, excitonic, and electron-phonon effects. We show that monolayer MoS2 possesses a large and diverse number of strongly bound excitonic states with novel k-space characteristics that were not previously seen experimentally or theoretically. The absorption spectrum is shown to be dominated by excitonic states with binding energy close to $1$~eV and by strong electron-phonon broadening in the visible to ultra-violet range. Our results explain recent experimental measurements and resolve inconsistencies between previous GW-BSE calculations.

preprint2013arXiv

Photoelectron spin-flipping and texture manipulation in a topological insulator

Recently discovered materials called three-dimensional topological insulators constitute examples of symmetry protected topological states in the absence of applied magnetic fields and cryogenic temperatures. A hallmark characteristic of these non-magnetic bulk insulators is the protected metallic electronic states confined to the material's surfaces. Electrons in these surface states are spin polarized with their spins governed by their direction of travel (linear momentum), resulting in a helical spin texture in momentum space. Spin- and angle-resolved photoemission spectroscopy (spin-ARPES) has been the only tool capable of directly observing this central feature with simultaneous energy, momentum, and spin sensitivity. By using an innovative photoelectron spectrometer with a high-flux laser-based light source, we discovered another surprising property of these surface electrons which behave like Dirac fermions. We found that the spin polarization of the resulting photoelectrons can be fully manipulated in all three dimensions through selection of the light polarization. These surprising effects are due to the spin-dependent interaction of the helical Dirac fermions with light, which originates from the strong spin-orbit coupling in the material. Our results illustrate unusual scenarios in which the spin polarization of photoelectrons is completely different from the spin state of electrons in the originating initial states. The results also provide the basis for a novel source of highly spin-polarized electrons with tunable polarization in three dimensions.

preprint2013arXiv

Physical origin of satellite in photoemission of doped graphene: An ab-initio GW plus cumulant study

We calculate the photoemission spectra of suspended and epitaxial doped graphene using an ab initio cumulant expansion of the Green's function based on the GW self-energy. Our results are compared to experiment and to standard GW calculations. For doped graphene on a silicon carbide substrate, we find, in contrast to earlier calculations, that the spectral function from GW only does not reproduce experimental satellite properties. However, ab initio GW plus cumulant theory combined with an accurate description of the substrate screening results in good agreement with experiment, but gives no plasmaron (i.e., no extra well-defined excitation satisfying Dyson's equation).

preprint2013arXiv

Surface Atom Motion to Move Iron Nanocrystals through Constrictions in Carbon Nanotubes under the Action of an Electric Current

Under the application of electrical currents, metal nanocrystals inside carbon nanotubes can be bodily transported. We examine experimentally and theoretically how an iron nanocrystal can pass through a constriction in the carbon nanotube with a smaller cross-sectional area than the nanocrystal itself. Remarkably, through in situ transmission electron imaging and diffraction, we find that, while passing through a constriction, the nanocrystal remains largely solid and crystalline and the carbon nanotube is unaffected. We account for this behavior by a pattern of iron atom motion and rearrangement on the surface of the nanocrystal. The nanocrystal motion can be described with a model whose parameters are nearly independent of the nanocrystal length, area, temperature, and electromigration force magnitude. We predict that metal nanocrystals can move through complex geometries and constrictions, with implications for both nanomechanics and tunable synthesis of metal nanoparticles.

preprint2013arXiv

Systematic Determination of Absolute Absorption Cross-section of Individual Carbon Nanotubes

Determination of optical absorption cross-section is always among the central importance of understanding a material. However its realization on individual nanostructures, such as carbon nanotubes, is experimentally challenging due to the small extinction signal using conventional transmission measurements. Here we develop a technique based on polarization manipulation to enhance the sensitivity of single-nanotube absorption spectroscopy by two-orders of magnitude. We systematically determine absorption cross-section over broad spectral range at single-tube level for more than 50 chirality-defined single-walled nanotubes. Our data reveals chirality-dependent one-dimensional photo-physics through the behaviours of exciton oscillator strength and lifetime. We also establish an empirical formula to predict absorption spectrum of any nanotube, which provides the foundation to determine quantum efficiencies in important photoluminescence and photovoltaic processes.

preprint2013arXiv

Theoretical study of solid iron nanocrystal movement inside a carbon nanotube

We use a first-principles based kinetic Monte Carlo simulation to study the movement of a solid iron nanocrystal inside a carbon nanotube driven by the electrical current. The origin of the iron nanocrystal movement is the electromigration force. Even though the iron nanocrystal appears to be moving as a whole, we find that the core atoms of the nanocrystal is completely stationary, and only the surface atoms are moving. Movement in the contact region with the carbon nanotube is driven by electromigration forces, and the movement on the remaining surfaces is driven by diffusion. Results of our calculations also provide a simple model which can predict the center of mass speed of the iron nanocrystal over a wide range of parameters. We find both qualitative and quantitative agreement of the iron nanocrystal center of mass speed with experimental data.

preprint2013arXiv

Theory of the Raman spectrum of rotated double-layer graphene

We study theoretically the Raman spectrum of the rotated double-layer graphene, consisting of two graphene layers rotated with respect to each other by an arbitrary angle θ. We find a relatively simple dependence of the Raman G peak intensity on the angle θ. On the other hand, the Raman 2D peak position, intensity, and width show a much more complicated dependence on the angle θ. We account for all of these effects, including dependence on the incoming photon energy, in good agreement with the experimental data. We find that it is sufficient to include the interaction between the graphene layers on the electronic degrees of freedom (resulting in the occurrence of van Hove singularities in the density of states). We assume that the phonon degrees of freedom are unaffected by the interaction between the layers. Furthermore, we decompose the Raman 2D peak into two components having very different linewidths; these widths are almost independent of the angle θ. The change in the intensity and the peak position of one of these two components gives insight into the dependence of the overall Raman 2D peak features as a function of the angle θ. Additionally, we show regions in the phonon spectrum giving rise to the Raman 2D peak signal. This work provides an insight into the interplay between the mechanical degree of freedom (angle θ) and the electronic degrees of freedom (singularities in the density of states) in rotated double-layer graphene. Additionally, this work provides a way to establish experimentally the value of the rotation angle θusing Raman spectroscopy measurement. This procedure becomes even more robust if one repeats the Raman spectroscopy measurement with a different incoming photon energy.

preprint2013arXiv

Tuning many-body interactions in graphene: The effects of doping on excitons and carrier lifetimes

The optical properties of graphene are strongly affected by electron-electron (e-e) and electron-hole (e-h) interactions. Here we tune these many-body interactions through varying the density of free charge carriers. Measurements from the infrared to the ultraviolet reveal significant changes in the optical conductivity of graphene for both electron and hole doping. The shift, broadening, and modification in shape of the saddle-point exciton resonance reflect strong screening of the many-body interactions by the carriers, as well as changes in quasi-particle lifetimes. Ab initio calculations by the GW Bethe-Salpeter equation (GW-BSE), which take into account modification of both the repulsive e-e and the attractive e-h interactions, provide excellent agreement with experiment. Understanding the optical properties and high-energy carrier dynamics of graphene over a wide range of doping is crucial for both fundamental graphene physics and for emerging applications of graphene in photonics.

preprint2012arXiv

BerkeleyGW: A Massively Parallel Computer Package for the Calculation of the Quasiparticle and Optical Properties of Materials and Nanostructures

BerkeleyGW is a massively parallel computational package for electron excited-state properties that is based on the many-body perturbation theory employing the ab initio GW and GW plus Bethe-Salpeter equation methodology. It can be used in conjunction with many density-functional theory codes for ground-state properties, including PARATEC, PARSEC, Quantum ESPRESSO, OCTOPUS and SIESTA. The package can be used to compute the electronic and optical properties of a wide variety of material systems from bulk semiconductors and metals to nanostructured materials and molecules. The package scales to 10,000's of CPUs and can be used to study systems containing up to 100's of atoms.

preprint2012arXiv

Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach

Ab initio GW calculations are a standard method for computing the spectroscopic properties of many materials. The most computationally expensive part in conventional implementations of the method is the generation and summation over the large number of empty orbitals required to converge the electron self energy. We propose a scheme to reduce the summation over empty states by the use of a modified static-remainder approximation, which is simple to implement and yields accurate self energies for both bulk and molecular systems requiring a small fraction of the typical number of empty orbitals.

preprint2012arXiv

First-principles calculations of the electronic structure of open-shell condensed matter systems

We develop a Green's function approach to quasiparticle excitations of open-shell systems within the GW approximation. It is shown that accurate calculations of the characteristic multiplet structure require a precise knowledge of the self energy and, in particular, its poles. We achieve this by constructing the self energy from appropriately chosen mean-field theories on a fine frequency grid. We apply our method to a two-site Hubbard model, several molecules and the negatively charged nitrogen-vacancy defect in diamond, and obtain good agreement with experiment and other high-level theories.

preprint2012arXiv

Phonon-assisted optical absorption in silicon from first principles

The phonon-assisted interband optical absorption spectrum of silicon is calculated at the quasiparticle level entirely from first principles. We make use of the Wannier interpolation formalism to determine the quasiparticle energies, as well as the optical transition and electron-phonon coupling matrix elements, on fine grids in the Brillouin zone. The calculated spectrum near the onset of indirect absorption is in very good agreement with experimental measurements for a range of temperatures. Moreover, our method can accurately determine the optical absorption spectrum of silicon in the visible range, an important process for optoelectronic and photovoltaic applications that cannot be addressed with simple models. The computational formalism is quite general and can be used to understand the phonon-assisted absorption processes in general.

preprint2012arXiv

Raman spectroscopy study of rotated double-layer graphene: Misorientation-angle dependence of electronic structure

We present a systematic Raman study of unconventionally-stacked double-layer graphene, and find that the spectrum strongly depends on the relative rotation angle between layers. Rotation-dependent trends in the position, width and intensity of graphene 2D and G peaks are experimentally established and accounted for theoretically. Our theoretical analysis reveals that changes in electronic band structure due to the interlayer interaction, such as rotational-angle dependent Van Hove singularities, are responsible for the observed spectra features. Our combined experimental and theoretical study provides a deeper understanding of the electronic band structure of rotated double-layer graphene, and leads to a practical way to identify and analyze rotation angles of misoriented double-layer graphene.

preprint2012arXiv

Resonant Excitation of Graphene K-Phonon and Intra-Landau-Level Excitons in Magneto-Optical Spectroscopy

Precise infrared magnetotransmission experiments have been performed in magnetic fields up to 32 T on a series of multilayer epitaxial graphene samples. We observe changes in the spectral features and broadening of the main cyclotron transition when the incoming photon energy is in resonance with the lowest Landau level separation and the energy of a K point optical phonon. We have developed a theory that explains and quantitatively reproduces the frequency and magnetic field dependence of the phenomenon as the absorption of a photon together with the simultaneous creation of an intervalley, intra-Landau-level exciton, and a K phonon.

preprint2012arXiv

Spin Polarization of Photoelectrons from Topological Insulators

We show that the degree of spin polarization of photoelectrons from the surface states of topological insulators is 100% if fully-polarized light is used as in typical photoemission measurements, and hence can be significantly "higher" than that of the initial state. Further, the spin orientation of these photoelectrons in general can also be very different from that of the initial surface state and is controlled by the photon polarization; a rich set of predicted phenomena have recently been confirmed by spin- and angle-resolved photoemission experiments.

preprint2011arXiv

Anomalous insulator metal transition in boron nitride-graphene hybrid atomic layers

The study of two-dimensional (2D) electronic systems is of great fundamental significance in physics. Atomic layers containing hybridized domains of graphene and hexagonal boron nitride (h-BNC) constitute a new kind of disordered 2D electronic system. Magneto-electric transport measurements performed at low temperature in vapor phase synthesized h-BNC atomic layers show a clear and anomalous transition from an insulating to a metallic behavior upon cooling. The observed insulator to metal transition can be modulated by electron and hole doping and by the application of an external magnetic field. These results supported by ab-initio calculations suggest that this transition in h-BNC has distinctly different characteristics when compared to other 2D electron systems and is the result of the coexistence between two distinct mechanisms, namely, percolation through metallic graphene networks and hopping conduction between edge states on randomly distributed insulating h-BN domains.

preprint2011arXiv

Direct measurement of quantum phases in graphene via photoemission spectroscopy

Quantum phases provide us with important information for understanding the fundamental properties of a system. However, the observation of quantum phases, such as Berry's phase and the sign of the matrix element of the Hamiltonian between two non-equivalent localized orbitals in a tight-binding formalism, has been challenged by the presence of other factors, e.g., dynamic phases and spin/valley degeneracy, and the absence of methodology. Here, we report a new way to directly access these quantum phases, through polarization-dependent angle-resolved photoemission spectroscopy (ARPES), using graphene as a prototypical two-dimensional material. We show that the momentum- and polarization-dependent spectral intensity provides direct measurements of (i) the phase of the band wavefunction and (ii) the sign of matrix elements for non-equivalent orbitals. Upon rotating light polarization by π/2, we found that graphene with a Berry's phase of nπ(n=1 for single- and n=2 for double-layer graphene for Bloch wavefunction in the commonly used form) exhibits the rotation of ARPES intensity by π/n, and that ARPES signals reveal the signs of the matrix elements in both single- and double-layer graphene. The method provides a new technique to directly extract fundamental quantum electronic information on a variety of materials.

preprint2011arXiv

Excitonic effects in the optical properties of SiC sheet and nanotubes

The quasiparticle band structure and optical properties of single-walled zigzag and armchair SiC nanotubes (SiC-NTs) as well as single SiC sheet are investigated by ab initio many-body calculations using the GW and the GW plus Bethe-Salpeter equation (GW+BSE) approaches, respectively. Significant GW quasiparticle corrections of more than 1.0 eV to the Kohn-Sham band gaps from the local density approximation (LDA) calculations are found. The GW self-energy corrections transform the SiC sheet from a indirect LDA band gap to a direct band gap material. Furthermore, the quasiparticle band gaps of SiC-NTs with different chiralities behave very differently as a function of tube diameter, and this can be attributed to the difference in the curvature-induced orbital rehybridization between the different chiral nanotubes. The calculated optical absorption spectra are dominated by discrete exciton peaks due to exciton states with large binding energy up to 2.0 eV in the SiC sheet and SiC-NTs. The formation of strongly bound excitons is attributed to the enhanced electron-hole interaction in these low dimensional systems. Remarkably, the excited electron amplitude of the exciton wavefunction is found to peak on the Si atoms near the hole position (which is on the C site) in the zigzag SiC-NTs, indicating a charge transfer from an anion (hole) to its neighboring cations by photoexcitation. In contrast, this pronounced peak structure disappear in the exciton wavefunction in the armchair SiC-NTs. Furthermore, in the armchair SiC-NTs, the bound exciton wavefunctions are more localized and also strongly cylindrically asymmetric.

preprint2011arXiv

Many-body interactions in quasi-freestanding graphene

The Landau-Fermi liquid picture for quasiparticles assumes that charge carriers are dressed by many-body interactions, forming one of the fundamental theories of solids. Whether this picture still holds for a semimetal like graphene at the neutrality point, i.e., when the chemical potential coincides with the Dirac point energy, is one of the long-standing puzzles in this field. Here we present such a study in quasi-freestanding graphene by using high-resolution angle-resolved photoemission spectroscopy. We see the electron-electron and electron-phonon interactions go through substantial changes when the semimetallic regime is approached, including renormalizations due to strong electron-electron interactions with similarities to marginal Fermi liquid behavior. These findings set a new benchmark in our understanding of many-body physics in graphene and a variety of novel materials with Dirac fermions.

preprint2011arXiv

Quasiparticle Effects in the Bulk and Surface-State Bands of Bi2Se3 and Bi2Te3 Topological Insulators

We investigate the bulk band structures and the surface states of Bi2Se3 and Bi2Te3 topological insulators using first-principles many-body perturbation theory based on the GW approximation. The quasiparticle self-energy corrections introduce significant changes to the bulk band structures, while their effect on the band gaps is opposite in the band-inversion regime compared to the usual situation without band inversion. Parametrized "scissors operators" derived from the bulk studies are then used to investigate the electronic structure of slab models which exhibit topologically protected surface states. The results including self-energy corrections reveal significant shifts of the Dirac points relative to the bulk bands and large gap openings resulting from the interactions between the surface states across the thin slab, both in agreement with experimental data.

preprint2011arXiv

Quasiparticle Excitations and Charge Transition Levels of Oxygen Vacancies in Hafnia

We calculate the quasiparticle defect states and charge transition levels of oxygen vacancies in monoclinic hafnia. The charge transition levels, although they are thermodynamic quantities, can be critically dependent on the band gap owing to localized defect states. These quasiparticle defect level effects are treated using the first principle GW approximation to the self energy. We show that the quality and reliability of the results may be evaluated by calculating the same transition level via two physical paths and that it is important to include the necessary electrostatic corrections in a supercell calculation. Contrary to many previous reports, the oxygen vacancies in monoclinic hafnia are found to be a positive U center, where U is the defect electron addition energy. We identify a physical partitioning of U in terms of an electronic and structural relaxation part.

preprint2011arXiv

Spatially Resolving Spin-split Edge States of Chiral Graphene Nanoribbons

A central question in the field of graphene-related research is how graphene behaves when it is patterned at the nanometer scale with different edge geometries. Perhaps the most fundamental shape relevant to this question is the graphene nanoribbon (GNR), a narrow strip of graphene that can have different chirality depending on the angle at which it is cut. Such GNRs have been predicted to exhibit a wide range of behaviour (depending on their chirality and width) that includes tunable energy gaps and the presence of unique one-dimensional (1D) edge states with unusual magnetic structure. Most GNRs explored experimentally up to now have been characterized via electrical conductivity, leaving the critical relationship between electronic structure and local atomic geometry unclear (especially at edges). Here we present a sub-nm-resolved scanning tunnelling microscopy (STM) and spectroscopy (STS) study of GNRs that allows us to examine how GNR electronic structure depends on the chirality of atomically well-defined GNR edges. The GNRs used here were chemically synthesized via carbon nanotube (CNT) unzipping methods that allow flexible variation of GNR width, length, chirality, and substrate. Our STS measurements reveal the presence of 1D GNR edge states whose spatial characteristics closely match theoretical expectations for GNR's of similar width and chirality. We observe width-dependent splitting in the GNR edge state energy bands, providing compelling evidence of their magnetic nature. These results confirm the novel electronic behaviour predicted for GNRs with atomically clean edges, and thus open the door to a whole new area of applications exploiting the unique magnetoelectronic properties of chiral GNRs.

preprint2011arXiv

Theory of Magnetic Edge States in Chiral Graphene Nanoribbons

Using a model Hamiltonian approach including electron-electron interactions, we systematically investigate the electronic structure and magnetic properties of chiral graphene nanoribbons. We show that the presence of magnetic edge states is an intrinsic feature of smooth graphene nanoribbons with chiral edges, and discover a number of structure-property relations. Specifically, we study the dependence of magnetic moments and edge-state energy splittings on the nanoribbon width and chiral angle as well as the role of environmental screening effects. Our results address a recent experimental observation of signatures of magnetic ordering in chiral graphene nanoribbons and provide an avenue towards tuning their properties via the structural and environmental degrees of freedom.

preprint2011arXiv

Theory of the electronic and transport properties of graphene under a periodic electric or magnetic field

We discuss the novel electronic properties of graphene under an external periodic scalar or vector potential, and the analytical and numerical methods used to investigate them. When graphene is subjected to a one-dimensional periodic scalar potential, owing to the linear dispersion and the chiral (pseudospin) nature of the electronic states, the group velocity of its carriers is renormalized highly anisotropically in such a manner that the velocity is invariant along the periodic direction but is reduced the most along the perpendicular direction. Under a periodic scalar potential, new massless Dirac fermions are generated at the supercell Brillouin zone boundaries. Also, we show that if the strength of the applied scalar potential is sufficiently strong, new zero-energy modes may be generated. With the periodic scalar potential satisfying some special conditions, the energy dispersion near the Dirac point becomes quasi one-dimensional. On the other hand, for graphene under a one-dimensional periodic vector potential (resulting in a periodic magnetic field perpendicular to the graphene plane), the group velocity is reduced isotropically and monotonically with the strength of the potential.

preprint2010arXiv

Angle-resolved photoemission spectra of graphene from first-principles calculations

Angle-resolved photoemission spectroscopy (ARPES) is a powerful experimental technique for directly probing electron dynamics in solids. The energy vs. momentum dispersion relations and the associated spectral broadenings measured by ARPES provide a wealth of information on quantum many-body interaction effects. In particular, ARPES allows studies of the Coulomb interaction among electrons (electron-electron interactions) and the interaction between electrons and lattice vibrations (electron-phonon interactions). Here, we report ab initio simulations of the ARPES spectra of graphene including both electron-electron and electron-phonon interactions on the same footing. Our calculations reproduce some of the key experimental observations related to many-body effects, including the indication of a mismatch between the upper and lower halves of the Dirac cone.

preprint2010arXiv

Calcium-Decorated Graphene-Based Nanostructures for Hydrogen Storage

We report a first-principles study of hydrogen storage media consisting of calcium atoms and graphene-based nanostructures. We find that Ca atoms prefer to be individually adsorbed on the zigzag edge of graphene with a Ca-Ca distance of 10 A without clustering of the Ca atoms, and up to six H2 molecules can bind to a Ca atom with a binding energy of ~0.2 eV/H2. A Ca-decorated zigzag graphene nanoribbon (ZGNR) can reach the gravimetric capacity of ~5 wt % hydrogen. We also consider various edge geometries of the graphene for Ca dispersion.

preprint2010arXiv

Electronic transport in polycrystalline graphene

Most materials in available macroscopic quantities are polycrystalline. Graphene, a recently discovered two-dimensional form of carbon with strong potential for replacing silicon in future electronics, is no exception. There is growing evidence of the polycrystalline nature of graphene samples obtained using various techniques. Grain boundaries, intrinsic topological defects of polycrystalline materials, are expected to dramatically alter the electronic transport in graphene. Here, we develop a theory of charge carrier transmission through grain boundaries composed of a periodic array of dislocations in graphene based on the momentum conservation principle. Depending on the grain boundary structure we find two distinct transport behaviours - either high transparency, or perfect reflection of charge carriers over remarkably large energy ranges. First-principles quantum transport calculations are used to verify and further investigate this striking behaviour. Our study sheds light on the transport properties of large-area graphene samples. Furthermore, purposeful engineering of periodic grain boundaries with tunable transport gaps would allow for controlling charge currents without the need of introducing bulk band gaps in otherwise semimetallic graphene. The proposed approach can be regarded as a means towards building practical graphene electronics.

preprint2010arXiv

EPW: A program for calculating the electron-phonon coupling using maximally localized Wannier functions

EPW (Electron-Phonon coupling using Wannier functions) is a program written in FORTRAN90 for calculating the electron-phonon coupling in periodic systems using density-functional perturbation theory and maximally-localized Wannier functions. EPW can calculate electron-phonon interaction self-energies, electron-phonon spectral functions, and total as well as mode-resolved electron-phonon coupling strengths. The calculation of the electron-phonon coupling requires a very accurate sampling of electron-phonon scattering processes throughout the Brillouin zone, hence reliable calculations can be prohibitively time-consuming. EPW combines the Kohn-Sham electronic eigenstates and the vibrational eigenmodes provided by the Quantum-ESPRESSO package [1] with the maximally localized Wannier functions provided by the wannier90 package [2] in order to generate electron-phonon matrix elements on arbitrarily dense Brillouin zone grids using a generalized Fourier interpolation. This feature of EPW leads to fast and accurate calculations of the electron-phonon coupling, and enables the study of the electron-phonon coupling in large and complex systems.

preprint2010arXiv

Gate-Controlled Ionization and Screening of Cobalt Adatoms on a Graphene Surface

We describe scanning tunneling spectroscopy (STS) measurements performed on individual cobalt (Co) atoms deposited onto backgated graphene devices. We find that Co adatoms on graphene can be ionized by either the application of a global backgate voltage or by the application of a local electric field from a scanning tunneling microscope (STM) tip. Large screening clouds are observed to form around Co adatoms ionized in this way, and we observe that some intrinsic graphene defects display a similar behavior. Our results provide new insight into charged impurity scattering in graphene, as well as the possibility of using graphene devices as chemical sensors.

preprint2010arXiv

Graphene Dirac fermions in one-dimensional inhomogeneous field profiles: Transforming magnetic to electric field

We show that the low-energy electronic structure of graphene under a one-dimensional inhomogeneous magnetic field can be mapped into that of graphene under an electric field or vice versa. As a direct application of this transformation, we find that the carrier velocity in graphene is isotropically reduced under magnetic fields periodic along one direction with zero average flux. This counterintuitive renormalization has its origin in the pseudospin nature of graphene electronic states and is robust against disorder. In magnetic graphene superlattices with a finite average flux, the Landau level bandwidth at high fields exhibits an unconventional behavior of decreasing with increasing strength of the average magnetic field due to the linear energy dispersion of graphene. As another application of our transformation relation, we show that the transmission probabilities of an electron through a magnetic barrier in graphene can directly be obtained from those through an electrostatic barrier or vice versa.

preprint2010arXiv

Selective Functionalization of Halogens on Zigzag Graphene Nanoribbons: A Route to the Separation of Zigzag Graphene Nanoribbons

Using the ab initio pseudopotential density functional method, we investigate the functionalization of halogen molecules into graphene-based nanostructures with zigzag and armchair edges. We find that halogen molecules adsorb through chemisorption on the zigzag edge carbon atoms with a binding energy of ~1-5 eV and their adsorption on a perfect zigzag edge is preferred, in sharp contrast to physisorption on the armchair edge and elsewhere where they adsorb with a binding energy of ~0.07 eV. We suggest that our findings would be utilized for an approach to the separation of zigzag graphene nanoribbons with regular edges with the change of the solubility of the functionalized nanoribbons.

preprint2010arXiv

Spin Polarization and Transport of Surface States in the Topological Insulators Bi2Se3 and Bi2Te3 from First Principles

We investigate the band dispersion and the spin texture of topologically protected surface states in the bulk topological insulators Bi2Se3 and Bi2Te3 by first-principles methods. Strong spin-orbit entanglement in these materials reduces the spin-polarization of the surface states to ~50% in both cases; this reduction is absent in simple models but of important implications to essentially any spintronic application. We propose a way of controlling the magnitude of spin polarization associated with a charge current in thin films of topological insulators by means of an external electric field. The proposed dual-gate device configuration provides new possibilities for electrical control of spin.

preprint2010arXiv

Topological Defects in Graphene: Dislocations and Grain Boundaries

Topological defects in graphene, dislocations and grain boundaries, are still not well understood despites the considerable number of experimental observations. We introduce a general approach for constructing dislocations in graphene characterized by arbitrary Burgers vectors as well as grain boundaries, covering the whole range of possible misorientation angles. By using ab initio calculations we investigate thermodynamic and electronic properties of these topological defects, finding energetically favorable symmetric large-angle grain boundaries, strong tendency towards out-of-plane deformation in the small-angle regimes, and pronounced effects on the electronic structure. The present results show that dislocations and grain boundaries are important intrinsic defects in graphene which may be used for engineering graphene-based nanomaterials and functional devices.

preprint2010arXiv

Tunable Excitons in Biased Bilayer Graphene

Recent measurements have shown that a continuously tunable bandgap of up to 250 meV can be generated in biased bilayer graphene [Y. Zhang et al., Nature 459, 820 (2009)], opening up pathway for possible graphene-based nanoelectronic and nanophotonic devices operating at room temperature. Here, we show that the optical response of this system is dominated by bound excitons. The main feature of the optical absorbance spectrum is determined by a single symmetric peak arising from excitons, a profile that is markedly different from that of an interband transition picture. Under laboratory conditions, the binding energy of the excitons may be tuned with the external bias going from zero to several tens of meV's. These novel strong excitonic behaviors result from a peculiar, effective ``one-dimensional'' joint density of states and a continuously-tunable bandgap in biased bilayer graphene. Moreover, we show that the electronic structure (level degeneracy, optical selection rules, etc.) of the bound excitons in a biased bilayer graphene is markedly different from that of a two-dimensional hydrogen atom because of the pseudospin physics.

preprint2009arXiv

A Tunable Phonon-Exciton Fano System in Bilayer Graphene

Interference between different possible paths lies at the heart of quantum physics. Such interference between coupled discrete and continuum states of a system can profoundly change its interaction with light as seen in Fano resonance. Here we present a unique many-body Fano system composed of a discrete phonon vibration and continuous electron-hole pair transitions in bilayer graphene. Mediated by the electron-phonon interactions, the excited state is described by new quanta of elementary excitations of hybrid phonon-exciton nature. Infrared absorption of the hybrid states exhibit characteristic Fano lineshapes with parameters renormalized by many-body interactions. Remarkably, the Fano resonance in bilayer graphene is continuously tunable through electrical gating. Further control of the phonon-exciton coupling may be achieved with an optical field exploiting the excited state infrared activity. This tunable phonon-exciton system also offers the intriguing possibility of a 'phonon laser' with stimulated phonon amplification generated by population inversion of band-edge electrons.

preprint2009arXiv

Fermi surfaces and quantum oscillations in underdoped high-Tc superconductors YBa2Cu3O6.5 and YBa2Cu4O8

We study underdoped high-Tc superconductors YBa2Cu3O6.5 and YBa2Cu4O8 using first-principles pseudopotential methods with additional Coulomb interactions at the Cu atoms, and obtain Fermi-surface pocket areas in close agreement with measured Shubnikov-de Haas and de Haas-van Alphen oscillations. With antiferromagnetic order in CuO2 planes, stable in the calculations, small hole pockets are formed near the so-called Fermi-arc positions in the Brillouin zone which reproduce the low-frequency oscillations. A large electron pocket, necessary for the negative Hall coefficient, is also formed in YBa2Cu3O6.5, giving rise to the high-frequency oscillations as well. Effective masses and specific heats are also calculated and compared with measurements. Our results highlight the important role of magnetic order in the electronic structure of underdoped high-Tc superconductors.

preprint2009arXiv

GW approach to Anderson model out of equilibrium: Coulomb blockade and false hysteresis in the I-V characteristics

The Anderson model for a single impurity coupled to two leads is studied using the $GW$ approximation in the strong electron-electron interaction regime as a function of the alignment of the impurity level relative to the chemical potentials in the leads. We employ a non-equilibrium Green's function technique to calculate the electron self-energy, the spin density and the current as a function of bias across the junction. In addition we develop an expression for the change in the expectation value of the energy of the system that results when the impurity is coupled to the leads, including the role of Coulomb interactions through the electron self energy in the region of the junction. The current-voltage characteristics calculated within the GW approximation exhibit Coulomb blockade. Depending on the gate voltage and applied bias, we find that there can be more than one steady-state solution for the system, which may give rise to a hysteresis in the I-V characteristics. We show that the hysteresis is an artifact of the $GW$ approximation and would not survive if quantum fluctuations beyond the $GW$ approximation are included.

preprint2009arXiv

GW method with the self-consistent Sternheimer equation

We propose a novel approach to quasiparticle GW calculations which does not require the computation of unoccupied electronic states. In our approach the screened Coulomb interaction is evaluated by solving self-consistent linear-response Sternheimer equations, and the noninteracting Green's function is evaluated by solving inhomogeneous linear systems. The frequency-dependence of the screened Coulomb interaction is explicitly taken into account. In order to avoid the singularities of the screened Coulomb interaction the calculations are performed along the imaginary axis, and the results are analytically continued to the real axis through Pade' approximants. As a proof of concept we implemented the proposed methodology within the empirical pseudopotential formalism and we validated our implementation using silicon as a test case. We examine the advantages and limitations of our method and describe promising future directions.

preprint2009arXiv

Localization of Metal-Induced Gap States at the Metal-Insulator Interface:Origin of Flux Noise in SQUIDs and Superconducting Qubits

The origin of magnetic flux noise in Superconducting Quantum Interference Devices with a power spectrum scaling as $1/f$ ($f$ is frequency) has been a puzzle for over 20 years. This noise limits the decoherence time of superconducting qubits. A consensus has emerged that the noise arises from fluctuating spins of localized electrons with an areal density of $5\times10^{17}$m$^{-2}$. We show that, in the presence of potential disorder at the metal-insulator interface, some of the metal-induced gap states become localized and produce local moments. A modest level of disorder yields the observed areal density.

preprint2007arXiv

Amine-Gold Linked Single-Molecule Junctions: Experiment and Theory

The measured conductance distribution for single molecule benzenediamine-gold junctions, based on 59,000 individual conductance traces recorded while breaking a gold point contact in solution, has a clear peak at 0.0064 G$_{0}$ with a width of $\pm$ 40%. Conductance calculations based on density functional theory (DFT) for 15 distinct junction geometries show a similar spread. Differences in local structure have a limited influence on conductance because the amine-Au bonding motif is well-defined and flexible. The average calculated conductance (0.046 G$_{0}$) is seven times larger than experiment, suggesting the importance of many-electron corrections beyond DFT.

preprint2007arXiv

Reversible Photomechanical Switching of Individual Engineered Molecules at a Surface

We have observed reversible light-induced mechanical switching for a single organic molecule bound to a metal surface. Scanning tunneling microscopy (STM) was used to image the features of an individual azobenzene molecule on Au(111) before and after reversibly cycling its mechanical structure between trans and cis states using light. Azobenzene molecules were engineered to increase their surface photomechanical activity by attaching varying numbers of tert-butyl (TB) ligands ("legs") to the azobenzene phenyl rings. STM images show that increasing the number of TB legs "lifts" the azobenzene molecules from the substrate, thereby increasing molecular photomechanical activity by decreasing molecule-surface coupling.