Angel Rubio

Bio: Angel Rubio is an academic researcher from Max Planck Society. The author has contributed to research in topics: Density functional theory & Carbon nanotube. The author has an hindex of 110, co-authored 930 publications receiving 52731 citations. Previous affiliations of Angel Rubio include University of California, Berkeley & Universidad Autónoma Metropolitana.

Electronic excitations: density-functional versus many-body Green's-function approaches

Abstract: Electronic excitations lie at the origin of most of the commonly measured spectra. However, the first-principles computation of excited states requires a larger effort than ground-state calculations, which can be very efficiently carried out within density-functional theory. On the other hand, two theoretical and computational tools have come to prominence for the description of electronic excitations. One of them, many-body perturbation theory, is based on a set of Green’s-function equations, starting with a one-electron propagator and considering the electron-hole Green’s function for the response. Key ingredients are the electron’s self-energy S and the electron-hole interaction. A good approximation for S is obtained with Hedin’s GW approach, using density-functional theory as a zero-order solution. First-principles GW calculations for real systems have been successfully carried out since the 1980s. Similarly, the electron-hole interaction is well described by the Bethe-Salpeter equation, via a functional derivative of S. An alternative approach to calculating electronic excitations is the time-dependent density-functional theory (TDDFT), which offers the important practical advantage of a dependence on density rather than on multivariable Green’s functions. This approach leads to a screening equation similar to the Bethe-Salpeter one, but with a two-point, rather than a four-point, interaction kernel. At present, the simple adiabatic local-density approximation has given promising results for finite systems, but has significant deficiencies in the description of absorption spectra in solids, leading to wrong excitation energies, the absence of bound excitonic states, and appreciable distortions of the spectral line shapes. The search for improved TDDFT potentials and kernels is hence a subject of increasing interest. It can be addressed within the framework of many-body perturbation theory: in fact, both the Green’s functions and the TDDFT approaches profit from mutual insight. This review compares the theoretical and practical aspects of the two approaches and their specific numerical implementations, and presents an overview of accomplishments and work in progress.

Theory of graphitic boron nitride nanotubes

Abstract: Based upon the similarities in properties between carbon- and BN-based (BN=boron nitride) materials, we propose that BN-based nanotubes can be stable and study their electronic structure. A simple Slater-Koster tight-binding scheme has been applied. All the BN nanotubes are found to be semiconducting materials. The band gaps are larger than 2 eV for most tubes. Depending on the helicity, the calculated band gap can be direct at [Gamma] or indirect. In general, the larger the diameter of the nanotube the larger the band gap, with a saturation value corresponding to the calculated local-density-approximation band gap of hexagonal BN. The higher ionicity of BN is important in explaining the electronic differences between these tubes and similar carbon nanotubes.

Stability and Band Gap Constancy of Boron Nitride Nanotubes

Abstract: Extensive LDA and quasi-particle calculations have been performed on boron nitride (BN) single-wall and multi-wall nanotubes. Strain energies are found to be smaller for BN nanotubes than for carbon nanotubes of the same radius, owing to a buckling effect which stabilizes the BN tubular structure. For tubes larger than 9.5 A in diameter, the lowest conduction band is predicted to be free-electron-like with electronic charge density localized inside the tube. For these tubes, this band is at constant energy above the top of the valence band. Consequently, contrarily to carbon nanotubes, single- and multi-wall BN nanotubes are constant-band-gap materials, independent of their radius and helicity. In addition, we expect them to exhibit remarkable properties under n-type doping.

Germanene: a novel two-dimensional germanium allotrope akin to graphene and silicene

Abstract: We have grown an atom-thin, ordered, two-dimensional multi-phase film in situ through germanium molecular beam epitaxy using a gold (111) surface as a substrate. Its growth is similar to the formation of silicene layers on silver (111) templates. One of the phases, forming large domains, as observed in scanning tunneling microscopy, shows a clear, nearly flat, honeycomb structure. Thanks to thorough synchrotron radiation core-level spectroscopy measurements and advanced density functional theory calculations we can identify it as a ?3????3 R(30?) germanene layer in conjunction with a ?7????7 R(19.1?) Au(111) supercell, presenting compelling evidence of the synthesis of the germanium-based cousin of graphene on gold.

I and i

Abstract: There is, I think, something ethereal about i —the square root of minus one. I remember first hearing about it at school. It seemed an odd beast at that time—an intruder hovering on the edge of reality. Usually familiarity dulls this sense of the bizarre, but in the case of i it was the reverse: over the years the sense of its surreal nature intensified. It seemed that it was impossible to write mathematics that described the real world in …

A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu

Abstract: The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.

Fast parallel algorithms for short-range molecular dynamics

Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

The electronic properties of graphene

Abstract: This article reviews the basic theoretical aspects of graphene, a one-atom-thick allotrope of carbon, with unusual two-dimensional Dirac-like electronic excitations. The Dirac electrons can be controlled by application of external electric and magnetic fields, or by altering sample geometry and/or topology. The Dirac electrons behave in unusual ways in tunneling, confinement, and the integer quantum Hall effect. The electronic properties of graphene stacks are discussed and vary with stacking order and number of layers. Edge (surface) states in graphene depend on the edge termination (zigzag or armchair) and affect the physical properties of nanoribbons. Different types of disorder modify the Dirac equation leading to unusual spectroscopic and transport properties. The effects of electron-electron and electron-phonon interactions in single layer and multilayer graphene are also presented.