Interlayer vacancy defects in AA-stacked bilayer graphene: density functional theory predictions.

AA-stacked graphite and closely related structures, where carbon atoms are located in registry in adjacent graphene layers, are a feature of graphitic systems including twisted and folded bilayer graphene, and turbostratic graphite. We present the results of ab initio density functional theory calculations performed to investigate the complexes that are formed from the binding of vacancy defects across neighbouring layers in AA-stacked bilayers. As with AB stacking, the carbon atoms surrounding lattice vacancies can form interlayer structures with sp 2 bonding that are lower in energy than in-plane reconstructions. The sp 2 interlayer bonding of adjacent multivacancy defects in registry creates a type of stable sp 2 bonded 'wormhole' or tunnel defect between the layers. We also identify a new class of 'mezzanine' structure characterised by sp 3 interlayer bonding, resembling a prismatic vacancy loop. The V 6 hexavacancy variant, where six sp 3 carbon atoms sit midway between two carbon layers and bond to both, is substantially more stable than any other vacancy aggregate in AA-stacked layers. Our focus is on vacancy generation and aggregation in the absence of extreme temperatures or intense beams.

The Stone-Wales transformation: from fullerenes to graphite, from radiation damage to heat capacity.

The Stone-Wales (SW) transformation, or carbon-bond rotation, has been fundamental to understanding fullerene growth and stability, and ab initio calculations show it to be a high-energy process. The nature and topology of the fullerene energy landscape shows how the Ih-C60 must be the final product, if SW transformations are fast enough, and various mechanisms for their catalysis have been proposed. We review SW transformations in fullerenes and then discuss the analogous transformation in graphite, where they form the Dienes defect, originally posited to be a transition state in the direct exchange of a bonded atom pair. On the basis of density functional theory calculations in the local density approximation, we propose that non-equilibrium concentrations of the Dienes defect arising from displacing radiation are rapidly healed by point defects and that equilibrium concentrations of Dienes defects are responsible for the divergent ultra-high-temperature heat capacity of graphite.This article is part of the themed issue 'Fullerenes: past, present and future, celebrating the 30th anniversary of Buckminster Fullerene'.

On the validity of empirical potentials for simulating radiation damage in graphite: a benchmark.

In this work, the ability of methods based on empirical potentials to simulate the effects of radiation damage in graphite is examined by comparing results for point defects, found using ab initio calculations based on density functional theory (DFT), with those given by two state of the art potentials: the Environment-Dependent Interatomic Potential (EDIP) and the Adaptive Intermolecular Reactive Empirical Bond Order potential (AIREBO). Formation energies for the interstitial, the vacancy and the Stone-Wales (5775) defect are all reasonably close to DFT values. Both EDIP and AIREBO can thus be suitable for the prompt defects in a cascade, for example. Both potentials suffer from arefacts. One is the pinch defect, where two α-atoms adopt a fourfold-coordinated sp(3) configuration, that forms a cross-link between neighbouring graphene sheets. Another, for AIREBO only, is that its ground state vacancy structure is close to the transition state found by DFT for migration. The EDIP fails to reproduce the ground state self-interstitial structure given by DFT, but has nearly the same formation energy. Also, for both potentials, the energy barriers that control diffusion and the evolution of a damage cascade, are not well reproduced. In particular the EDIP gives a barrier to removal of the Stone-Wales defect as 0.9 eV against DFT's 4.5 eV. The suite of defect structures used is provided as supplementary information as a benchmark set for future potentials.

Extended interplanar linking in graphite formed from vacancy aggregates.

The mechanical and electrical properties of graphite and related materials such as multilayer graphene depend strongly on the presence of defects in the lattice structure, particularly those which create links between adjacent planes. We present findings which suggest the existence of a new type of defect in the graphite or graphene structure which connects adjacent planes through continuous hexagonal sp2 bonding alone and can form through the aggregation of individual vacancy defects. The energetics and kinetics of the formation of this type of defect are investigated with atomistic density functional theory calculations. The resultant structures are then employed to simulate high resolution transmission electron microscopy images, which are compared to recent experimental images of electron irradiation damaged graphite.

The contribution made by lattice vacancies to the Wigner effect in radiation-damaged graphite.

Models for radiation damage in graphite are reviewed and compared, leading to a re-examination of the contribution made by vacancies to annealing processes. A method based on density functional theory, using large supercells with orthorhombic and hexagonal symmetry, is employed to calculate the properties and behaviour of lattice vacancies and displacement defects. It is concluded that annihilation of intimate Frenkel defects marks the onset of recovery in electrical resistivity, which occurs when the temperature exceeds about 160 K. The migration of isolated monovacancies is estimated to have an activation energy of E(a) ≈ 1.1 eV. Coalescence into divacancy defects occurs in several stages, with different barriers at each stage, depending on the path. The formation of pairs ultimately yields up to 8.9 eV energy, which is nearly 1.0 eV more than the formation energy for an isolated monovacancy. Processes resulting in vacancy coalescence and annihilation appear to be responsible for the main Wigner energy release peak in radiation-damaged graphite, occurring at about 475 K.

Muonium chemistry and spin dynamics in sulphur, modelling interstitial hydrogen.

The nature of the elusive muonium centre in sulphur is re-examined in the light of new data on its level-crossing resonance and spin-lattice relaxation. The aim is to provide a model for the solid-state chemistry of interstitial hydrogen in this element, which is as yet unknown, as well as to solve one of the longest standing puzzles in µSR spectroscopy, namely the surprisingly strong depolarization of muons mimicking ion-implanted protons in this innocuous non-magnetic material. The paramagnetic muonium (and by inference hydrogen) centre is confirmed to have the character of a molecular radical, but with huge anisotropy at cryogenic temperatures and a striking shift of the resonance at ordinary temperatures, the hyperfine parameters appearing to collapse and vanish towards the melting point. New density-functional supercell calculations identify a number of possible structures for the defect centre, including a novel form of bond-centred muonium in a closed-ring S(8)Mu complex. Simulations of the spin dynamics and fits to the spectra suggest a dynamical equilibrium or chemical exchange between several configurations, with occupancy of the bond-centre site falling from unity at low cryogenic temperatures to zero near the melting point.

First-principles simulations of boron diffusion in graphite.

Boron strongly modifies electronic and diffusion properties of graphite. We report the first ab initio study of boron interaction with the point defects in graphite, which includes structures, thermodynamics, and diffusion. A number of possible diffusion mechanisms of boron in graphite are suggested. We conclude that boron diffuses in graphite by a kick-out mechanism. This mechanism explains the common activation energy, but large magnitude difference, for the rate of boron diffusion parallel and perpendicular to the basal plane.

Metastable Frenkel pair defect in graphite: source of Wigner energy?

The atomic processes associated with energy storage and release in irradiated graphite have long been subject to untested speculation. We examine structures and recombination routes for interstitial-vacancy (I-V) pairs in graphite. Interaction results in the formation of a new metastable defect (an intimate I-V pair) or a Stone-Wales defect. The intimate I-V pair, although 2.9 eV more stable than its isolated constituents, still has a formation energy of 10.8 eV. The barrier to recombination to perfect graphite is calculated to be 1.3 eV, consistent with the experimental first Wigner energy release peak at 1.38 eV. We expect similar defects to form in carbon nanostructures such as nanotubes, nested fullerenes, and onions under irradiation.