Energetics of intrinsic point defects in ZrSiO4_4

Using first principles calculations we have studied the formation energies, electron and hole affinities, and electronic levels of intrinsic point defects in zircon. The atomic structures of charged interstitials, vacancies, Frenkel pairs and anti-site defects are obtained. The limit of high concentration of point defects, relevant for the use of this material in nuclear waste immobilization, was studied with a variable lattice relaxation that can simulate the swelling induced by radiation damage. The limit of low concentration of defects is simulated with larger cells and fixed lattice parameters. Using known band offset values at the interface of zircon with silicon, we analyze the foreseeable effect of the defects on the electronic properties of zircon used as gate in metal-oxide-semiconductor devices.Comment: preprint 16 pages, 4 figures, and 5 table

Band bending and quasi-2DEG in the metallized β\beta-SiC(001) surface

We study the mechanism leading to the metallization of the $\beta$-SiC(001) Si-rich surface induced by hydrogen adsorption. We analyze the effects of band bending and demonstrate the existence of a quasi-2D electron gas, which originates from the donation of electrons from adsorbed hydrogen to bulk conduction states. We also provide a simple model that captures the main features of the results of first-principles calculations, and uncovers the basic physics of the process.Comment: accepted for publication in physica status solidi - Rapid Research Letter

Analysis of the Scanning Tunneling Microscopy Images of the Charge Density Wave Phase in Quasi-one-dimensional Rb0.3MoO3

The experimental STM images for the CDW phase of the blue bronze RbMoO3 have been successfully explained on the basis of first-principles DFT calculations. Although the density of states near the Fermi level strongly concentrates in two of the three types of Mo atoms Mo-II and Mo-III, the STM measurement mostly probes the contribution of the uppermost O atoms of the surface, associated with the Mo-IO6 octahedra. In addition, it is found that the surface concentration of Rb atoms plays a key role in determining the surface nesting vector and hence the periodicity of the CDW modulation. Significant experimental inhomogeneities of the b* surface component of the wavevector of the modulation, probed by STM, are reported. The calculated changes in the surface nesting vector are consistent with the observed experimental inhomogeneities.Comment: 4 pages 5 Figure

Interference effects in one-dimensional moiré crystals

Interference effects in finite sections of one-dimensional moiré crystals are investigated using a Landauer-Büttiker formalism within the tight-binding approximation. We explain interlayer transport in double-wall carbon nanotubes and design a predictive model. Wave function interference is visible at the mesoscale: in the strong coupling regime, as a periodic modulation of quantum conductance and emergent localized states; in the localized-insulating regime, as a suppression of interlayer transport, and oscillations of the density of states. These results could be exploited to design quantum electronic devices. © 2021 The Author

Total energy global optimizations using non orthogonal localized orbitals

An energy functional for orbital based $O(N)$ calculations is proposed, which depends on a number of non orthogonal, localized orbitals larger than the number of occupied states in the system, and on a parameter, the electronic chemical potential, determining the number of electrons. We show that the minimization of the functional with respect to overlapping localized orbitals can be performed so as to attain directly the ground state energy, without being trapped at local minima. The present approach overcomes the multiple minima problem present within the original formulation of orbital based $O(N)$ methods; it therefore makes it possible to perform $O(N)$ calculations for an arbitrary system, without including any information about the system bonding properties in the construction of the input wavefunctions. Furthermore, while retaining the same computational cost as the original approach, our formulation allows one to improve the variational estimate of the ground state energy, and the energy conservation during a molecular dynamics run. Several numerical examples for surfaces, bulk systems and clusters are presented and discussed.Comment: 24 pages, RevTex file, 5 figures available upon reques

New Superhard Phases for 3D C60-based Fullerites

We have explored new possible phases of 3D C60-based fullerites using semiempirical potentials and ab-initio density functional methods. We have found three closely related structures - two body centered orthorhombic and one body centered cubic - having 52, 56 and 60 tetracoordinated atoms per molecule. These 3D polymers result in semiconductors with bulk moduli near 300 GPa, and shear moduli around 240 GPa, which make them good candidates for new low density superhard materials.Comment: To be published in Physical Review Letter

Intrinsic point defects and volume swelling in ZrSiO4 under irradiation

The effects of high concentration of point defects in crystalline ZrSiO4 as originated by exposure to radiation, have been simulated using first principles density functional calculations. Structural relaxation and vibrational studies were performed for a catalogue of intrinsic point defects, with different charge states and concentrations. The experimental evidence of a large anisotropic volume swelling in natural and artificially irradiated samples is used to select the subset of defects that give similar lattice swelling for the concentrations studied, namely interstitials of O and Si, and the anti-site Zr(Si), Calculated vibrational spectra for the interstitials show additional evidence for the presence of high concentrations of some of these defects in irradiated zircon.Comment: 9 pages, 7 (color) figure

Basis Functions for Linear-Scaling First-Principles Calculations

In the framework of a recently reported linear-scaling method for density-functional-pseudopotential calculations, we investigate the use of localized basis functions for such work. We propose a basis set in which each local orbital is represented in terms of an array of `blip functions'' on the points of a grid. We analyze the relation between blip-function basis sets and the plane-wave basis used in standard pseudopotential methods, derive criteria for the approximate equivalence of the two, and describe practical tests of these criteria. Techniques are presented for using blip-function basis sets in linear-scaling calculations, and numerical tests of these techniques are reported for Si crystal using both local and non-local pseudopotentials. We find rapid convergence of the total energy to the values given by standard plane-wave calculations as the radius of the linear-scaling localized orbitals is increased.Comment: revtex file, with two encapsulated postscript figures, uses epsf.sty, submitted to Phys. Rev.

Bonding, Moment Formation, and Magnetic Interactions in Ca14MnBi11 and Ba14MnBi11

The ``14-1-11'' phase compounds based on magnetic Mn ions and typified by Ca14MnBi11 and Ba14MnBi11 show unusual magnetic behavior, but the large number (104) of atoms in the primitive cell has precluded any previous full electronic structure study. Using an efficient, local orbital based method within the local spin density approximation to study the electronic structure, we find a gap between a bonding valence band complex and an antibonding conduction band continuum. The bonding bands lack one electron per formula unit of being filled, making them low carrier density p-type metals. The hole resides in the MnBi4 tetrahedral unit and partially compensates the high spin d^5 Mn moment, leaving a net spin near 4 \mu_B that is consistent with experiment. These manganites are composed of two disjoint but interpenetrating `jungle gym' networks of spin 4/2 MnBi4^{9-} units with ferromagnetic interactions within the same network, and weaker couplings between the networks whose sign and magnitude is sensitive to materials parameters. Ca14MnBi11 is calculated to be ferromagnetic as observed, while for Ba14MnBi11 (which is antiferromagnetic) the ferro- and antiferromagnetic states are calculated to be essentially degenerate. The band structure of the ferromagnetic states is very close to half metallic.Comment: 17 pages, containing 10 postscript figures and 5 tables. Two additional figures (Fig.8 and 11 of the paper) are provided in JPG format in separate files. Submitted to Phys. Rev. B on September 20th 200