Nanocontacts: probing electronic structure under extreme uniaxial strains

Nanometer-sized metallic necks have the unique ability to sustain extreme uniaxial loads (about 20 times greater than the bulk material). We present an experimental and theoretical study of the electronic transport properties under such extreme conditions. Conductance measurements on gold and aluminum necks show a strikingly different behavior: While gold shows the expected conductance decrease with increasing elastic elongation of the neck, aluminum necks behave in the opposite way. We have performed first-principles electronic-structure calculations which reproduce this behavior, showing that it is an intrinsic property of the bulk band structure under high uniaxial strain.

Emergence of half-metallicity in suspended NiO chains: Ab initio electronic structure and quantum transport calculations

Contrary to the antiferromagnetic and insulating character of bulk NiO, one-dimensional chains of this material can become half metallic due to the lower coordination of their atoms. Here we present ab initio electronic structure and quantum transport calculations of ideal infinitely long NiO chains and of more realistic short ones suspended between Ni electrodes. While infinite chains are insulating, short suspended chains are half-metallic minority-spin conductors that displays very large magnetoresistance and a spin-valve behavior controlled by a single atom.

Electronic structure and transport properties of atomic NiO spinvalves

Ab initio quantum transport calculations show that short NiO chains suspended in Ni nanocontacts present a very strong spin-polarization of the conductance.The generalized gradient approximation we use here predicts a similar polarization of the conductance as the one previously computed with non-local exchange, confirming the robustness of the result. Their use as nanoscopic spinvalves is proposed.

Electronic structure of gated graphene and graphene ribbons

We study the electronic structure of gated graphene sheets. We consider both infinite graphene and finite width ribbons. The effect of Coulomb interactions between the electrically injected carriers and the coupling to the external gate are computed self-consistently in the Hartree approximation. We compute the average density of extra carriers n2D, the number of occupied subbands, and the density profiles as a function of the gate potential Vg. We discuss quantum corrections to the classical capacitance and we calculate the threshold Vg above which semiconducting armchair ribbons conduct. We find that the ideal conductance of perfectly transmitting wide ribbons is proportional to the square root of the gate voltage.

Role of potassium orbitals in the metallic behavior of K3picene

Detailed electronic structure calculations of picene clusters doped by potassium modeling the crystalline K3picene structure show that while two electrons are completely transferred from potassium atoms to the lowest-energy unoccupied molecular orbital of pristine picene, the third one remains closely attached to both material components. Multiconfigurational analysis is necessary to show that many structures of almost degenerate total energies compete to define the cluster ground state. Our results prove that the 4s orbital of potassium should be included in any interaction model describing the material. We propose a quarter-filled two-orbital model as the most simple model capable of describing the electronic structure of K-intercalated picene. Precise solutions obtained by a development of the Lanczos method show low-energy electronic excitations involving orbitals located at different positions. Consequently, metallic transport is possible in spite of the clear dominance of interaction over hopping.

Electronic properties of the MoS2-WS2 heterojunction

We study the electronic structure of a heterojunction made of two monolayers of MoS2 and WS2. Our first-principles density functional calculations show that, unlike in the homogeneous bilayers, the heterojunction has an optically active band gap, smaller than the ones of MoS2 and WS2 single layers. We find that the optically active states of the maximum valence and minimum conduction bands are localized on opposite monolayers, and thus the lowest energy electron-holes pairs are spatially separated. Our findings portray the MoS2-WS2 bilayer as a prototypical example for band-gap engineering of atomically thin two-dimensional semiconducting heterostructures.

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.

On the forbidden gap of finite graphene nanoribbons

The electronic structure of isolated finite graphene nanoribbons is investigated by solving, at the Hartree-Fock (HF) level, the Pariser, Parr and Pople (PPP) many-body Hamiltonian. The study is mainly focused on 7-AGNR and 13-AGNR (Armchair Graphene Nano-Ribbons), whose electronic structures have been recently experimentally investigated. Only paramagnetic solutions are considered. The characteristics of the forbidden gap are studied as a function of the ribbon length. For a 7-AGNR, the gap monotonically decreases from a maximum value of ~6.5 eV for short nanoribbons to a very small value of ~0.12 eV for the longer calculated systems. Gap edges are defined by molecular orbitals that are spatially localized near the nanoribbon extremes, that is, near both zig-zag edges. On the other hand, two delocalized orbitals define a much larger gap of about 5 eV. Conductance measurements report a somewhat smaller gap of ~3 eV. The small real gap lies in the middle of the one given by extended states and has been observed by STM and reproduced by DFT calculations. On the other hand, the length dependence of the gap is not monotonous for a 13-AGNR. It decreases initially but sharply increases for lengths beyond 30 Å remaining almost constant thereafter at a value of ~2.1 eV. Two additional states localized at the nanoribbon extremes show up at energies 0.31 eV below the HOMO (Highest Occupied Molecular Orbital) and above the LUMO (Lowest Unoccupied Molecular Orbital). These numbers compare favorably with those recently obtained by means of STS for a 13-AGNR sustained by a gold surface, namely 1.4 eV for the energy gap and 0.4 eV for the position of localized band edges. We show that the important differences between 7- and 13-AGNR should be ascribed to the charge rearrangement near the zig-zag edges obtained in our calculations for ribbons longer than 30 Å, a feature that does not show up for a 7-AGNR no matter its length.

Double Hybrid Functionals and the Π-System Bond Length Alternation Challenge: Rivaling Accuracy of Post-HF Methods

Predicting accurate bond length alternations (BLAs) in long conjugated oligomers has been a significant challenge for electronic-structure methods for many decades, made particularly important by the close relationships between BLA and the rich optoelectronic properties of π-delocalized systems. Here, we test the accuracy of recently developed, and increasingly popular, double hybrid (DH) functionals, positioned at the top of Jacobs Ladder of DFT methods of increasing sophistication, computational cost, and accuracy, due to incorporation of MP2 correlation energy. Our test systems comprise oligomeric series of polyacetylene, polymethineimine, and polysilaacetylene up to six units long. MP2 calculations reveal a pronounced shift in BLAs between the 6-31G(d) basis set used in many studies of BLA to date and the larger cc-pVTZ basis set, but only modest shifts between cc-pVTZ and aug-cc-pVQZ results. We hence perform new reference CCSD(T)/cc-pVTZ calculations for all three series of oligomers against which we assess the performance of several families of DH functionals based on BLYP, PBE, and TPSS, along with lower-rung relatives including global- and range-separated hybrids. Our results show that DH functionals systematically improve the accuracy of BLAs relative to single hybrid functionals. xDH-PBE0 (N4 scaling using SOS-MP2) emerges as a DH functional rivaling the BLA accuracy of SCS-MP2 (N5 scaling), which was found to offer the best compromise between computational cost and accuracy the last time the BLA accuracy of DFT- and wave function-based methods was systematically investigated. Interestingly, xDH-PBE0 (XYG3), which differs to other DHs in that its MP2 term uses PBE0 (B3LYP) orbitals that are not self-consistent with the DH functional, is an outlier of trends of decreasing average BLA errors with increasing fractions of MP2 correlation and HF exchange.

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.

Coherent transport in graphene nanoconstrictions

We study the effect of a structural nanoconstriction on the coherent transport properties of otherwise ideal zigzag-edged infinitely long graphene ribbons. The electronic structure is calculated with the standard one-orbital tight-binding model and the linear conductance is obtained using the Landauer formula. We find that, since the zero-bias current is carried in the bulk of the ribbon, this is very robust with respect to a variety of constriction geometries and edge defects. In contrast, the curve of zero-bias conductance versus gate voltage departs from the (2n+1)e2∕h staircase of the ideal case as soon as a single atom is removed from the sample. We also find that wedge-shaped constrictions can present nonconducting states fully localized in the constriction close to the Fermi energy. The interest of these localized states in regards to the formation of quantum dots in graphene is discussed.

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.

Mn-doped II-VI quantum dots: artificial molecular magnets

The notion of artificial atom relies on the capability to change the number of carriers one by one in semiconductor quantum dots, and the resulting changes in their electronic structure. Organic molecules with transition metal atoms that have a net magnetic moment and display hysteretic behaviour are known as single molecule magnets (SMM). The fabrication of CdTe quantum dots chemically doped with a controlled number of Mn atoms and with a number of carriers controlled either electrically or optically paves the way towards a new concept in nanomagnetism: the artificial single molecule magnet. Here we study the magnetic properties of a Mn-doped CdTe quantum dot for different charge states and show to what extent they behave like a single molecule magnet.

Graphene single-electron transistor as a spin sensor for magnetic adsorbates

We study single-electron transport through a graphene quantum dot with magnetic adsorbates. We focus on the relation between the spin order of the adsorbates and the linear conductance of the device. The electronic structure of the graphene dot with magnetic adsorbates is modeled through numerical diagonalization of a tight-binding model with an exchange potential. We consider several mechanisms by which the adsorbate magnetic state can influence transport in a single-electron transistor: tuning the addition energy, changing the tunneling rate, and in the case of spin-polarized electrodes, through magnetoresistive effects. Whereas the first mechanism is always present, the others require that the electrode has to have either an energy- or spin-dependent density of states. We find that graphene dots are optimal systems to detect the spin state of a few magnetic centers.

Modeling contact formation between atomic-sized gold tips via molecular dynamics

The formation and rupture of atomic-sized contacts is modelled by means of molecular dynamics simulations. Such nano-contacts are realized in scanning tunnelling microscope and mechanically controlled break junction experiments. These instruments routinely measure the conductance across the nano-sized electrodes as they are brought into contact and separated, permitting conductance traces to be recorded that are plots of conductance versus the distance between the electrodes. One interesting feature of the conductance traces is that for some metals and geometric configurations a jump in the value of the conductance is observed right before contact between the electrodes, a phenomenon known as jump-to-contact. This paper considers, from a computational point of view, the dynamics of contact between two gold nano-electrodes. Repeated indentation of the two surfaces on each other is performed in two crystallographic orientations of face-centred cubic gold, namely (001) and (111). Ultimately, the intention is to identify the structures at the atomic level at the moment of first contact between the surfaces, since the value of the conductance is related to the minimum cross-section in the contact region. Conductance values obtained in this way are determined using first principles electronic transport calculations, with atomic configurations taken from the molecular dynamics simulations serving as input structures.