Can model Hamiltonians describe the electron–electron interaction in π-conjugated systems?: PAH and graphene

Model Hamiltonians have been, and still are, a valuable tool for investigating the electronic structure of systems for which mean field theories work poorly. This review will concentrate on the application of Pariser–Parr–Pople (PPP) and Hubbard Hamiltonians to investigate some relevant properties of polycyclic aromatic hydrocarbons (PAH) and graphene. When presenting these two Hamiltonians we will resort to second quantisation which, although not the way chosen in its original proposal of the former, is much clearer. We will not attempt to be comprehensive, but rather our objective will be to try to provide the reader with information on what kinds of problems they will encounter and what tools they will need to solve them. One of the key issues concerning model Hamiltonians that will be treated in detail is the choice of model parameters. Although model Hamiltonians reduce the complexity of the original Hamiltonian, they cannot be solved in most cases exactly. So, we shall first consider the Hartree–Fock approximation, still the only tool for handling large systems, besides density functional theory (DFT) approaches. We proceed by discussing to what extent one may exactly solve model Hamiltonians and the Lanczos approach. We shall describe the configuration interaction (CI) method, a common technology in quantum chemistry but one rarely used to solve model Hamiltonians. In particular, we propose a variant of the Lanczos method, inspired by CI, that has the novelty of using as the seed of the Lanczos process a mean field (Hartree–Fock) determinant (the method will be named LCI). Two questions of interest related to model Hamiltonians will be discussed: (i) when including long-range interactions, how crucial is including in the Hamiltonian the electronic charge that compensates ion charges? (ii) Is it possible to reduce a Hamiltonian incorporating Coulomb interactions (PPP) to an 'effective' Hamiltonian including only on-site interactions (Hubbard)? The performance of CI will be checked on small molecules. The electronic structure of azulene and fused azulene will be used to illustrate several aspects of the method. As regards graphene, several questions will be considered: (i) paramagnetic versus antiferromagnetic solutions, (ii) forbidden gap versus dot size, (iii) graphene nano-ribbons, and (iv) optical properties.

Polarized interacting exciton gas in quantum wells and bulk semiconductors

We develop a theory to calculate exciton binding energies of both two- and three-dimensional spin polarized exciton gases within a mean field approach. Our method allows the analysis of recent experiments showing the importance of the polarization and intensity of the excitation light on the exciton luminescence of GaAs quantum wells. We study the breaking of the spin degeneracy observed at high exciton density (5×1010 cm2). Energy level splitting between spin +1 and spin -1 is shown to be due to many-body interexcitonic exchange while the spin relaxation time is controlled by intraexciton exchange. © 1996 The American Physical Society.

Theory of ferromagnetism in planar heterostructures of (Mn,III)-V semiconductors

A density-functional theory of ferromagnetism in heterostructures of compound semiconductors doped with magnetic impurities is presented. The variable functions in the density-functional theory are the charge and spin densities of the itinerant carriers and the charge and localized spins of the impurities. The theory is applied to study the Curie temperature of planar heterostructures of III-V semiconductors doped with manganese atoms. The mean-field, virtual-crystal and effective-mass approximations are adopted to calculate the electronic structure, including the spin-orbit interaction, and the magnetic susceptibilities, leading to the Curie temperature. By means of these results, we attempt to understand the observed dependence of the Curie temperature of planar δ-doped ferromagnetic structures on variation of their properties. We predict a large increase of the Curie temperature by additional confinement of the holes in a δ-doped layer of Mn by a quantum well.

Spin separation in digital ferromagnetic heterostructures

In a study of the ferromagnetic phase of a multilayer digital ferromagnetic semiconductor in the mean-field and effective-mass approximations, we find the exchange interaction to have the dominant energy scale of the problem, effectively controlling the spatial distribution of the carrier spins in the digital ferromagnetic heterostructures. In the ferromagnetic phase, the majority-spin and minority-spin carriers tend to be in different regions of the space (spin separation). Hence, the charge distribution of carriers also changes noticeably from the ferromagnetic to the paramagnetic phase. An example of a design to exploit these phenomena is given here.

Prediction of hidden multiferroic order in graphene zigzag ribbons

An electronic phase with coexisting magnetic and ferroelectric order is predicted for graphene ribbons with zigzag edges. The electronic structure of the system is described with a mean-field Hubbard model that yields results very similar to those of density functional calculations. Without further approximations, the mean-field theory is recasted in terms of a BCS wave function for electron-hole pairs in the edge bands. The BCS coherence present in each spin channel is related to spin-resolved electric polarization. Although the total electric polarization vanishes, due to an internal phase locking of the BCS state, strong magnetoelectric effects are expected in this system. The formulation naturally accounts for the two gaps in the quasiparticle spectrun, Δ0 and Δ1, and relates them to the intraband and interband self-energies.

Vacancy-induced magnetism in graphene and graphene ribbons

We address the electronic structure and magnetic properties of vacancies and voids both in graphene and graphene ribbons. By using a mean-field Hubbard model, we study the appearance of magnetic textures associated with removing a single atom (vacancy) and multiple adjacent atoms (voids) as well as the magnetic interactions between them. A simple set of rules, based on the Lieb theorem, link the atomic structure and the spatial arrangement of the defects to the emerging magnetic order. The total spin S of a given defect depends on its sublattice imbalance, but some defects with S=0 can still have local magnetic moments. The sublattice imbalance also determines whether the defects interact ferromagnetically or antiferromagnetically with one another and the range of these magnetic interactions is studied in some simple cases. We find that in semiconducting armchair ribbons and two-dimensional graphene without global sublattice imbalance, there is a maximum defect density above which local magnetization disappears. Interestingly, the electronic properties of semiconducting graphene ribbons with uncoupled local moments are very similar to those of diluted magnetic semiconductors, presenting giant Zeeman splitting.

Interplay between sublattice and spin symmetry breaking in graphene

We study the effect of sublattice symmetry breaking on the electronic, magnetic, and transport properties of two-dimensional graphene as well as zigzag terminated one- and zero-dimensional graphene nanostructures. The systems are described with the Hubbard model within the collinear mean field approximation. We prove that for the noninteracting bipartite lattice with an unequal number of atoms in each sublattice, in-gap states still exist in the presence of a staggered on-site potential ±Δ/2. We compute the phase diagram of both 2D and 1D graphene with zigzag edges, at half filling, defined by the normalized interaction strength U/t and Δ/t, where t is the first neighbor hopping. In the case of 2D we find that the system is always insulating, and we find the Uc(Δ) curve above which the system goes antiferromagnetic. In 1D we find that the system undergoes a phase transition from nonmagnetic insulator for U

Magnetism in graphene nanoislands

We study the magnetic properties of nanometer-sized graphene structures with triangular and hexagonal shapes terminated by zigzag edges. We discuss how the shape of the island, the imbalance in the number of atoms belonging to the two graphene sublattices, the existence of zero-energy states, and the total and local magnetic moment are intimately related. We consider electronic interactions both in a mean-field approximation of the one-orbital Hubbard model and with density functional calculations. Both descriptions yield values for the ground state total spin S consistent with Lieb’s theorem for bipartite lattices. Triangles have a finite S for all sizes whereas hexagons have S=0 and develop local moments above a critical size of ≈1.5  nm.

Tunnel magnetoresistance in GaMnAs: going beyond Jullière formula

The relation between tunnel magnetoresistance (TMR) and spin polarization is explored for GaMnAs∕GaAlAs∕GaMnAs structures where the carriers experience strong spin–orbit interactions. TMR is calculated using the Landauer approach. The materials are described in the 6 band k⋅p model which includes spin–orbit interaction. Ferromagnetism is described in the virtual crystal mean field approximations. Our results indicate that TMR is a function of spin polarization and barrier thickness. As a result of the stong spin–orbit interactions, TMR also depends on the the angle between current flow direction and the electrode magnetization. These results compromise the validity of Julliere formula.

Transport across two interacting quantum dots: Bulk Kondo, Kondo box, and molecular regimes

We analyze the transport properties of a double quantum dot device with both dots coupled to perfect conducting leads and to a finite chain of N noninteracting sites connecting both of them. The interdot chain strongly influences the transport across the system and the local density of states of the dots. We study the case of a small number of sites, so that Kondo box effects are present, varying the coupling between the dots and the chain. For odd N and small coupling between the interdot chain and the dots, a state with two coexisting Kondo regimes develops: the bulk Kondo due to the quantum dots connected to leads and the one produced by the screening of the quantum dot spins by the spin in the finite chain at the Fermi level. As the coupling to the interdot chain increases, there is a crossover to a molecular Kondo effect, due to the screening of the molecule (formed by the finite chain and the quantum dots) spin by the leads. For even N the two Kondo temperatures regime does not develop and the physics is dominated by the usual competition between Kondo and antiferromagnetism between the quantum dots. We finally study how the transport properties are affected as N is increased. For the study we used exact multiconfigurational Lanczos calculations and finite-U slave-boson mean-field theory at T=0. The results obtained with both methods describe qualitatively and also quantitatively the same physics.

Noncollinear magnetic phases and edge states in graphene quantum Hall bars

Application of a perpendicular magnetic field to charge neutral graphene is expected to result in a variety of broken symmetry phases, including antiferromagnetic, canted, and ferromagnetic. All these phases open a gap in bulk but have very different edge states and noncollinear spin order, recently confirmed experimentally. Here we provide an integrated description of both edge and bulk for the various magnetic phases of graphene Hall bars making use of a noncollinear mean field Hubbard model. Our calculations show that, at the edges, the three types of magnetic order are either enhanced (zigzag) or suppressed (armchair). Interestingly, we find that preformed local moments in zigzag edges interact with the quantum spin Hall like edge states of the ferromagnetic phase and can induce backscattering.