Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is known about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they underestimate charge gaps.

Universal zero-bias conductance through a quantum wire side-coupled to a quantum dot

A numerical renormalization-group study of the conductance through a quantum wire containing noninteracting electrons side-coupled to a quantum dot is reported. The temperature and the dot-energy dependence of the conductance are examined in the light of a recently derived linear mapping between the temperature-dependent conductance and the universal function describing the conductance for the symmetric Anderson model of a quantum wire with an embedded quantum dot. Two conduction paths, one traversing the wire, the other a bypass through the quantum dot, are identified. A gate potential applied to the quantum wire is shown to control the current through the bypass. When the potential favors transport through the wire, the conductance in the Kondo regime rises from nearly zero at low temperatures to nearly ballistic at high temperatures. When it favors the dot, the pattern is reversed: the conductance decays from nearly ballistic to nearly zero. When comparable currents flow through the two channels, the conductance is nearly temperature independent in the Kondo regime, and Fano antiresonances in the fixed-temperature plots of the conductance as a function of the dot-energy signal interference between them. Throughout the Kondo regime and, at low temperatures, even in the mixed-valence regime, the numerical data are in excellent agreement with the universal mapping.

Nonempirical hyper-generalized-gradient functionals constructed from the Lieb-Oxford bound

A simple and completely general representation of the exact exchange-correlation functional of density-functional theory is derived from the universal Lieb-Oxford bound, which holds for any Coulomb-interacting system. This representation leads to an alternative point of view on popular hybrid functionals, providing a rationale for why they work and how they can be constructed. A similar representation of the exact correlation functional allows to construct fully nonempirical hyper-generalized-gradient approximations (HGGAs), radically departing from established paradigms of functional construction. Numerical tests of these HGGAs for atomic and molecular correlation energies and molecular atomization energies show that even simple HGGAs match or outperform state-of-the-art correlation functionals currently used in solid-state physics and quantum chemistry.

Electronic structure of In(2)O(3) and Sn-doped In(2)O(3) by hard x-ray photoemission spectroscopy

The valence and core levels of In(2)O(3) and Sn-doped In(2)O(3) have been studied by hard x-ray photoemission spectroscopy (hv = 6000 eV) and by conventional Al K alpha (hv = 1486.6 eV) x-ray photoemission spectroscopy. The experimental spectra are compared with density-functional theory calculations. It is shown that structure deriving from electronic levels with significant In or Sn 5s character is selectively enhanced under 6000 eV excitation. This allows us to infer that conduction band states in Sn-doped samples and states at the bottom of the valence band both contain a pronounced In 5s contribution. The In 3d core line measured at hv = 1486.6 eV for both undoped and Sn-doped In(2)O(3) display an asymmetric lineshape, and may be fitted with two components associated with screened and unscreened final states. The In 3d core line spectra excited at hv = 6000 eV for the Sn-doped samples display pronounced shoulders and demand a fit with two components. The In 3d core line spectrum for the undoped sample can also be fitted with two components, although the relative intensity of the component associated with the screened final state is low, compared to excitation at 1486.6 eV. These results are consistent with a high concentration of carriers confined close to the surface of nominally undoped In(2)O(3). This conclusion is in accord with the fact that a conduction band feature observed for undoped In(2)O(3) in Al K alpha x-ray photoemission is much weaker than expected in hard x-ray photoemission.

Theoretical investigation of atomic and electronic structures of Ga(2)O(3)(ZnO)(6)

Transparent conducting oxides (TCO) are widely used in technological applications ranging from photovoltaics to thin-film transparent field-effect transistors. In this work we report a first-principles investigation, based on density-functional theory, of the atomic and electronic properties of Ga(2)O(3)(ZnO)(6) (GZO(6)), which is a promising candidate to be used as host oxide for wide band gap TCO applications. We identify a low-energy configuration for the coherent distribution of the Ga and Zn atoms in the cation positions within the experimentally reported orthorhombic GZO(6) structure. Four Ga atoms are located in four-fold sites, while the remaining 12 Ga atoms in the unit cell form four shared Ga agglomerates (a motif of four atoms). The Zn atoms are distributed in the remaining cation sites with effective coordination numbers from 3.90 to 4.50. Furthermore, we identify the natural formation of twin-boundaries in GZO(6), which can explain the zigzag modulations observed experimentally by high-resolution transmission electron microscopy in GZO(n) (n=9). Due to the intrinsic twin-boundary formation, polarity inversion in the ZnO tetrahedrons is present which is facilitated by the formation of the Ga agglomerates. Our analysis shows that the formation of fourfold Ga sites and Ga agglomerates are stabilized by the electronic octet rule, while the distribution of Ga atoms and the formation of the twin-boundary help alleviate excess strain. Finally we identify that the electronic properties of GZO(6) are essentially determined by the electronic properties of ZnO, i.e., there are slight changes in the band gap and optical absorption properties.

Atomistic origins of the phase transition mechanism in Ge(2)Sb(2)Te(5)

The fast and reversible phase transition mechanism between crystalline and amorphous phases of Ge(2)Sb(2)Te(5) has been in debate for several years. Through employing first-principles density functional theory calculations, we identify a direct structural link between the metastable crystalline and amorphous phases. The phase transition is driven by the displacement of Ge atoms along the rocksalt [111] direction from stable octahedron to high energy unstable tetrahedron sites close to the intrinsic vacancy regions, which generates a high energy intermediate phase between metastable and amorphous phases. Due to the instability of Ge at the tetrahedron sites, the Ge atoms naturally shift away from those sites, giving rise to the formation of local-ordered fourfold motifs and the long-range structural disorder. Intrinsic vacancies, which originate from Sb(2)Te(3), lower the energy barrier for Ge displacements, and hence, their distribution plays an important role in the phase transition. The high energy intermediate configuration can be obtained experimentally by applying an intense laser beam, which overcomes the thermodynamic barrier from the octahedron to tetrahedron sites. The high figure of merit of Ge(2)Sb(2)Te(5) is achieved from the optimal combination of intrinsic vacancies provided by Sb(2)Te(3) and the instability of the tetrahedron sites provided by GeTe.

Quantitative determination of optical transmission through subwavelength slit arrays in Ag films: Role of surface wave interference and local coupling between adjacent slits

Measurement of the transmitted intensity from a coherent monomode light source through a series of subwavelength slit arrays in Ag films, with varying array pitch and number of slits, demonstrates enhancement (suppression) by factors of as much as 6 (9) when normalized to the transmission efficiency of an isolated slit. Pronounced minima in the transmitted intensity are observed at array pitches corresponding to lambda(SPP), 2 lambda(SPP), and 3 lambda(SPP), where lambda(SPP) is the wavelength of the surface plasmon polariton (SPP). The position of these minima arises from destructive interference between incident propagating waves and pi-phase-shifted SPP waves. Increasing the number of slits to four or more does not increase appreciably the per-slit transmission intensity. A simple interference model fits well the measured transmitted intensity profile.

Transition-metal 13-atom clusters assessed with solid and surface-biased functionals

First-principles density-functional theory studies have reported open structures based on the formation of double simple-cubic (DSC) arrangements for Ru(13), Rh(13), Os(13), and Ir(13), which can be considered an unexpected result as those elements crystallize in compact bulk structures such as the face-centered cubic and hexagonal close-packed lattices. In this work, we investigated with the projected augmented wave method the dependence of the lowest-energy structure on the local and semilocal exchange-correlation (xc) energy functionals employed in density-functional theory. We found that the local-density approximation (LDA) and generalized-gradient formulations with different treatment of the electronic inhomogeneities (PBE, PBEsol, and AM05) confirm the DSC configuration as the lowest-energy structure for the studied TM(13) clusters. A good agreement in the relative total energies are obtained even for structures with small energy differences, e. g., 0.10 eV. The employed xc functionals yield the same total magnetic moment for a given structure, i.e., the differences in the bond lengths do not affect the moments, which can be attributed to the atomic character of those clusters. Thus, at least for those systems, the differences among the LDA, PBE, PBEsol, and AM05 functionals are not large enough to yield qualitatively different results. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3577999]

The role of electron localization in the atomic structure of transition-metal 13-atom clusters: the example of Co(13), Rh(13), and Hf(13)

The crystalline structure of transition-metals (TM) has been widely known for several decades, however, our knowledge on the atomic structure of TM clusters is still far from satisfactory, which compromises an atomistic understanding of the reactivity of TM clusters. For example, almost all density functional theory (DFT) calculations for TM clusters have been based on local (local density approximation-LDA) and semilocal (generalized gradient approximation-GGA) exchange-correlation functionals, however, it is well known that plain DFT fails to correct the self-interaction error, which affects the properties of several systems. To improve our basic understanding of the atomic and electronic properties of TM clusters, we report a DFT study within two nonlocal functionals, namely, the hybrid HSE (Heyd, Scuseria, and Ernzerhof) and GGA + U functionals, of the structural and electronic properties of the Co(13), Rh(13), and Hf(13) clusters. For Co(13) and Rh(13), we found that improved exchange-correlation functionals decrease the stability of open structures such as the hexagonal bilayer (HBL) and double simple-cubic (DSC) compared with the compact icosahedron (ICO) structure, however, DFT-GGA, DFT-GGA + U, and DFT-HSE yield very similar results for Hf(13). Thus, our results suggest that the DSC structure obtained by several plain DFT calculations for Rh(13) can be improved by the use of improved functionals. Using the sd hybridization analysis, we found that a strong hybridization favors compact structures, and hence, a correct description of the sd hybridization is crucial for the relative energy stability. For example, the sd hybridization decreases for HBL and DSC and increases for ICO in the case of Co(13) and Rh(13), while for Hf(13), the sd hybridization decreases for all configurations, and hence, it does not affect the relative stability among open and compact configurations.

GLOBAL WELL-POSEDNESS AND NON-LINEAR STABILITY OF PERIODIC TRAVELING WAVES FOR A SCHRODINGER-BENJAMIN-ONO SYSTEM

The objective of this paper is two-fold: firstly, we develop a local and global (in time) well-posedness theory for a system describing the motion of two fluids with different densities under capillary-gravity waves in a deep water flow (namely, a Schrodinger-Benjamin-Ono system) for low-regularity initial data in both periodic and continuous cases; secondly, a family of new periodic traveling waves for the Schrodinger-Benjamin-Ono system is given: by fixing a minimal period we obtain, via the implicit function theorem, a smooth branch of periodic solutions bifurcating a Jacobian elliptic function called dnoidal, and, moreover, we prove that all these periodic traveling waves are nonlinearly stable by perturbations with the same wavelength.

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

Reconstruction of core and surface nanoparticles: The example of Pt(55) and Au(55)

Cuboctahedron (CUB) and icosahedron (ICO) model structures are widely used in the study of transition-metal (TM) nanoparticles (NPs), however, it might not provide a reliable description for small TM NPs such as the Pt(55) and Au(55) systems in gas phase. In this work, we combined density-functional theory calculations with atomic configurations generated by the basin hopping Monte Carlo algorithm within the empirical Sutton-Chen embedded atom potential. We identified alternative lower energy configurations compared with the ICO and CUB model structures, e. g., our lowest energy structures are 5.22 eV (Pt(55)) and 2.01 eV (Au(55)) lower than ICO. The energy gain is obtained by the Pt and Au diffusion from the ICO core region to the NP surface, which is driven by surface compression (only 12 atoms) on the ICO core region. Therefore, in the lowest energy configurations, the core size reduces from 13 atoms (ICO, CUB) to about 9 atoms while the NP surface increases from 42 atoms (ICO, CUB) to about 46 atoms. The present mechanism can provide an improved atom-level understanding of small TM NPs reconstructions.

Excitation energies from ground-state density-functionals by means of generator coordinates

The generator-coordinate method is a flexible and powerful reformulation of the variational principle. Here we show that by introducing a generator coordinate in the Kohn-Sham equation of density-functional theory, excitation energies can be obtained from ground-state density functionals. As a viability test, the method is applied to ground-state energies and various types of excited-state energies of atoms and ions from the He and the Li isoelectronic series. Results are compared to a variety of alternative DFT-based approaches to excited states, in particular time-dependent density-functional theory with exact and approximate potentials.

Bulk structures of PtO and PtO(2) from density functional calculations

Platinum plays an important role in catalysis and electrochemistry, and it is known that the direct interaction of oxygen with Pt surfaces can lead to the formation of platinum oxides (PtO(x)), which can affect the reactivity. To contribute to the atomistic understanding of the atomic structure of PtO(x), we report a density functional theory study of the atomic structure of bulk PtO(x) (1 <= x <= 2). From our calculations, we identified a lowest-energy structure (GeS type, space group Pnma) for PtO, which is 0.181 eV lower in energy than the structure suggested by W. J. Moore and L. Pauling [J. Am. Chem. Soc. 63, 1392 (1941)] (PtS type). Furthermore, two atomic structures were identified for PtO(2), which are almost degenerate in energy with the lowest-energy structure reported so far for PtO(2) (CaCl(2) type). Based on our results and analysis, we suggest that Pt and O atoms tend to form octahedron motifs in PtO(x) even at lower O composition by the formation of Pt-Pt bonds.

Abundance of sedentary consumers and sessile organisms along the wave exposure gradient of subtropical rocky shores of the south-west Atlantic

Sedentary consumers play an important role on populations of prey and, hence, their patterns of abundance, distribution and coexistence on shores are important to evaluate their potential influence on ecosystem dynamics. Here, we aimed to describe their spatio-temporal distribution and abundance in relation to wave exposure in the intertidal rocky shores of the south-west Atlantic to provide a basis for further understanding of ecological processes in this system. The abundance and composition of the functional groups of sessile organisms and sedentary consumers were taken by sampling the intertidal of sheltered and moderately exposed shores during a period of one year. The sublittoral fringe of sheltered areas was dominated by macroalgae, while the low midlittoral was dominated by bare rock and barnacles. In contrast, filter-feeding animals prevailed at exposed shores, probably explaining the higher abundance of the predator Stramonita haemastoma at these locations. Limpets were more abundant at the midlittoral zone of all shores while sea urchins were exclusively found at the sublittoral fringe of moderately exposed shores, therefore, adding grazing pressure on these areas. The results showed patterns of coexistence, distribution and abundance of those organisms in this subtropical area, presumably as a result of wave action, competition and prey availability. It also brought insights on the influence of top-down and bottom-up processes in this area.