The role of carbon impurities on the Si(001)-c(4 x 4) surface reconstruction: Theoretical calculations

In this work we employ the state-of-the-art pseudopotential method, within a generalized gradient approximation to the density functional theory, combined with a recently developed method for the calculation of HREELS spectra to study a series of different proposed models for carbon incorporation on the silicon (001) surface. A fully discussion on the geometry, energetics and specially the comparison between experimental and theoretical STM images and electron energy loss spectra indicate that the Si(100)-c(4 x 4) is probably induced by Si-C surface dinners, in agreement with recent experimental findings. (C) 2009 Elsevier B.V. All rights reserved.

Furan interaction with the Si(001)-(2 x 2) surface: structural, energetics, and vibrational spectra from first-principles

In this work we employ the state of the art pseudopotential method, within a generalized gradient approximation to the density functional theory, to investigate the adsorption process of furan on the silicon (001) surface. A direct comparison of different adsorption structures with x-ray photoelectron spectroscopy (XPS), ultra-violet photoelectron spectroscopy (UPS), high resolution electron energy loss spectroscopy (HREELS), near edge x-ray absorption fine structure (NEXAFS), and high resolution spectroscopy experimental data allows us to identify the [4 + 2] cycloaddition reaction as the most probable adsorbate. In addition, theoretical scanning tunnelling microscopy (STM) images are presented, with a view to contributing to further experimental investigations.

Electronic and optical properties of grossular garnet (Ca(3)Al(2)Si(3)O(12)): An ab initio study

The electronic and optical properties of grossular garnet are investigated using density functional theory (DFT) within generalized gradient approximation (GGA). The calculated lattice parameters are in good agreement with the experiment data. The electronic structure shows that grossular has a direct band gap of 5.22 eV. The dielectric functions, reflective index, extinction coefficient, reflectivity and energy-loss spectrum are calculated. The optical properties of grossular are discussed based on the band structure calculations. The O 2p states and Si 3s play a major role in these optical transitions as initial and final states, respectively. The absorption spectrum is localized in the ultraviolet range between 30 and 250 nm. Finally, we concluded that pure grossular crystal does not absorb radiation in the visible range. (c) 2009 Elsevier B.V. All rights reserved.

Dependence of Response Functions and Orbital Functionals on Occupation Numbers

Explicitly orbital-dependent approximations to the exchange-correlation energy functional of density functional theory typically not only depend on the single-particle Kohn-Sham orbitals but also on their occupation numbers in the ground-state Slater determinant. The variational calculation of the corresponding exchange-correlation potentials with the optimized effective potential (OEP) method therefore also requires a variation of the occupation numbers with respect to a variation in the effective single-particle potential, which is usually not taken into account. Here it is shown under which circumstances this procedure is justified.

Tightened Lieb-Oxford Bound for Systems of Fixed Particle Number

The Lieb-Oxford bound is a constraint upon approximate exchange-correlation functionals. We explore a nonempirical tightening of that bound in both universal and electron number-dependent form. The test functional is PBE. Regarding both atomization energies (slightly worsened) and bond lengths (slightly improved), we find the PBE functional to be remarkably insensitive to the value of the Lieb-Oxford bound. This both rationalizes the use of the original Lieb-Oxford constant in PBE and suggests that enhancement factors more sensitive to sharpened constraints await discovery.

Empirical analysis of the Lieb-Oxford bound in ions and molecules

Universal properties of the Coulomb interaction energy apply to all many-electron systems. Bounds on the exchange-correlation energy, in particular, are important for the construction of improved density functionals. Here we investigate one such universal property-the Lieb-Oxford lower bound-for ionic and molecular systems. In recent work [J Chem Phys 127, 054106 (2007)], we observed that for atoms and electron liquids this bound may be substantially tightened. Calculations for a few ions and molecules suggested the same tendency, but were not conclusive due to the small number of systems considered. Here we extend that analysis to many different families of ions and molecules, and find that for these, too, the bound can be empirically tightened by a similar margin as for atoms and electron liquids. Tightening the Lieb-Oxford bound will have consequences for the performance of various approximate exchange-correlation functionals. (C) 2008 Wiley Periodicals Inc.

Exact Search Directions for Optimization of Linear and Nonlinear Models Based on Generalized Orthonormal Functions

A novel technique for selecting the poles of orthonormal basis functions (OBF) in Volterra models of any order is presented. It is well-known that the usual large number of parameters required to describe the Volterra kernels can be significantly reduced by representing each kernel using an appropriate basis of orthonormal functions. Such a representation results in the so-called OBF Volterra model, which has a Wiener structure consisting of a linear dynamic generated by the orthonormal basis followed by a nonlinear static mapping given by the Volterra polynomial series. Aiming at optimizing the poles that fully parameterize the orthonormal bases, the exact gradients of the outputs of the orthonormal filters with respect to their poles are computed analytically by using a back-propagation-through-time technique. The expressions relative to the Kautz basis and to generalized orthonormal bases of functions (GOBF) are addressed; the ones related to the Laguerre basis follow straightforwardly as a particular case. The main innovation here is that the dynamic nature of the OBF filters is fully considered in the gradient computations. These gradients provide exact search directions for optimizing the poles of a given orthonormal basis. Such search directions can, in turn, be used as part of an optimization procedure to locate the minimum of a cost-function that takes into account the error of estimation of the system output. The Levenberg-Marquardt algorithm is adopted here as the optimization procedure. Unlike previous related work, the proposed approach relies solely on input-output data measured from the system to be modeled, i.e., no information about the Volterra kernels is required. Examples are presented to illustrate the application of this approach to the modeling of dynamic systems, including a real magnetic levitation system with nonlinear oscillatory behavior.

Systematics of nuclear densities, deformations and excitation energies within the context of the generalized rotation-vibration model

We present a large-scale systematics of charge densities, excitation energies and deformation parameters For hundreds of heavy nuclei The systematics is based on a generalized rotation vibration model for the quadrupole and octupole modes and takes into account second-order contributions of the deformations as well as the effects of finite diffuseness values for the nuclear densities. We compare our results with the predictions of classical surface vibrations in the hydrodynamical approximation. (C) 2010 Elsevier B V All rights reserved.

On the generalized eigenvalue method for energies and matrix elements in lattice field theory

We discuss the generalized eigenvalue problem for computing energies and matrix elements in lattice gauge theory, including effective theories such as HQET. It is analyzed how the extracted effective energies and matrix elements converge when the time separations are made large. This suggests a particularly efficient application of the method for which we can prove that corrections vanish asymptotically as exp(-(E(N+1) - E(n))t). The gap E(N+1) - E(n) can be made large by increasing the number N of interpolating fields in the correlation matrix. We also show how excited state matrix elements can be extracted such that contaminations from all other states disappear exponentially in time. As a demonstration we present numerical results for the extraction of ground state and excited B-meson masses and decay constants in static approximation and to order 1/m(b) in HQET.

HQET at order 1/m: II. Spectroscopy in the quenched approximation

Using Heavy Quark Effective Theory with non-perturbatively determined parameters in a quenched lattice calculation, we evaluate the splittings between the ground state and the first two radially excited states of the B(s) system at static order. We also determine the splitting between first excited and ground state, and between the B(s)* and B(s) ground states to order 1/m(b). The Generalized Eigenvalue Problem and the use of all-to-all propagators are important ingredients of our approach.