First principles investigation of finite-temperature behavior in small sodium clusters

A systematic and detailed investigation of the finite-temperature behavior of small sodium clusters, Na_n, in the size range of n= 8 to 50 are carried out. The simulations are performed using density-functional molecular-dynamics with ultrasoft pseudopotentials. A number of thermodynamic indicators such as specific heat, caloric curve, root-mean-square bond length fluctuation, deviation energy, etc. are calculated for each of the clusters. Size dependence of these indicators reveals several interesting features. The smallest clusters with n= 8 and 10, do not show any signature of melting transition. With the increase in size, broad peak in the specific heat is developed, which alternately for larger clusters evolves into a sharper one, indicating a solidlike to liquidlike transition. The melting temperatures show irregular pattern similar to experimentally observed one for larger clusters [ M. Schmidt et al., Nature (London) 393, 238 (1998) ]. The present calculations also reveal a remarkable size-sensitive effect in the size range of n= 40 to 55. While Na_40 and Na_55 show well developed peaks in the specific heat curve, Na_50 cluster exhibits a rather broad peak, indicating a poorly-defined melting transition. Such a feature has been experimentally observed for gallium and aluminum clusters [ G. A. Breaux et al., J. Am. Chem. Soc. 126, 8628 (2004); G. A.Breaux et al., Phys. Rev. Lett. 94, 173401 (2005) ].Comment: 8 pages, 11 figure

Kohn-Sham Exchange Potential for a Metallic Surface

The behavior of the surface barrier that forms at the metal-vacuum interface is important for several fields of surface science. Within the Density Functional Theory framework, this surface barrier has two non-trivial components: exchange and correlation. Exact results are provided for the exchange component, for a jellium metal-vacuum interface, in a slab geometry. The Kohn-Sham exact-exchange potential $V_{x}(z)$ has been generated by using the Optimized Effective Potential method, through an accurate numerical solution, imposing the correct boundary condition. It has been proved analytically, and confirmed numerically, that $V_{x}(z\to \infty)\to - e^{2}/z$; this conclusion is not affected by the inclusion of correlation effects. Also, the exact-exchange potential develops a shoulder-like structure close to the interface, on the vacuum side. The issue of the classical image potential is discussed.Comment: Phys. Rev. Lett. (to appear

Exchange parameters from approximate self-interaction correction scheme

The approximate atomic self-interaction corrections (ASIC) method to density functional theory is put to the test by calculating the exchange interaction for a number of prototypical materials, critical to local exchange and correlation functionals. ASIC total energy calculations are mapped onto an Heisenberg pair-wise interaction and the exchange constants J are compared to those obtained with other methods. In general the ASIC scheme drastically improves the bandstructure, which for almost all the cases investigated resemble closely available photo-emission data. In contrast the results for the exchange parameters are less satisfactory. Although ASIC performs reasonably well for systems where the magnetism originates from half-filled bands, it suffers from similar problems than those of LDA for other situations. In particular the exchange constants are still overestimated. This reflects a subtle interplay between exchange and correlation energy, not captured by the ASIC.Comment: 10 page

Origin of negative differential resistance in a strongly coupled single molecule-metal junction device

A new mechanism is proposed to explain the origin of negative differential resistance (NDR) in a strongly coupled single molecule-metal junction. A first-principles quantum transport calculation in a Fe-terpyridine linker molecule sandwiched between a pair of gold electrodes is presented. Upon increasing applied bias, it is found that a new phase in the broken symmetry wavefunction of the molecule emerges from the mixing of occupied and unoccupied molecular orbital. As a consequence, a non-linear change in the coupling between molecule and lead is evolved resulting to NDR. This model can be used to explain NDR in other class of metal-molecule junction device.Comment: Submitted for review on Feb 4, 200

The effects of temperature gradient and growth rate on the morphology and fatigue properties of MAR-M246(Hf)

MAR-M246(Hf) is a nickel based superalloy used in the turbopump blades of the Space Shuttle main engines. The effects are considered of temperature gradient (G) and growth rate (R) on the microstructure and fatigue properties of this superalloy. The primary dendrite arm spacings were found to be inversely proportional to both temperature gradient and growth rate. Carbide and gamma - gamma prime morphology trends were related to G/R ratios. Weibull analysis of fatigue results shows the characteristic life to be larger by a factor of 10 for the low gradient/fast rate pairing of G and R, while the reliability (beta) was lower

Novel properties of the Kohn-Sham exchange potential for open systems: application to the two-dimensional electron gas

The properties of the Kohn-Sham (KS) exchange potential for open systems in thermodynamical equilibrium, where the number of particles is non-conserved, are analyzed with the Optimized Effective Potential (OEP) method of Density Functional Theory (DFT) at zero temperature. The quasi two-dimensional electron gas (2DEG) is used as an illustrative example. The main findings are that the KS exchange potential builds a significant barrier-like structure under slight population of the second subband, and that both the asymptotic value of the KS exchange potential and the inter-subband energy jump discontinuously at the one-subband (1S) -> two-subband (2S) transition. The results obtained in this system offer new insights on open problems of semiconductors, such as the band-gap underestimation and the band-gap renormalization by photo-excited carriers.Comment: 7 pages, 3 figures, uses epl.cls(included), accepted for publication in Europhysics Letter

The quenching of compressible edge states around antidots

We provide a systematic quantitative description of the edge state structure around a quantum antidot in the integer quantum Hall regime. The calculations for spinless electrons within the Hartree approximation reveal that the widely used Chklovskii et al. electrostatic description greatly overestimates the widths of the compressible strips; the difference between these approaches diminishes as the size of the antidot increases. By including spin effects within density functional theory in the local spin-density approximation, we demonstrate that the exchange interaction can suppress the formation of compressible strips and lead to a spatial separation between the spin-up and spin-down states. As the magnetic field increases, the outermost compressible strip, related to spin-down states starts to form. However, in striking contrast to quantum wires, the innermost compressible strip (due to spin-up states) never develops for antidots.Comment: submitted to Phys. Rev. Let

Magnetosubband and edge state structure in cleaved-edge overgrown quantum wires

We provide a systematic quantitative description of the structure of edge states and magnetosubband evolution in hard wall quantum wires in the integer quantum Hall regime. Our calculations are based on the self-consistent Green's function technique where the electron- and spin interactions are included within the density functional theory in the local spin density approximation. We analyze the evolution of the magnetosubband structure as magnetic field varies and show that it exhibits different features as compared to the case of a smooth confinement. In particularly, in the hard-wall wire a deep and narrow triangular potential well (of the width of magnetic length $l_B$) is formed in the vicinity of the wire boundary. The wave functions are strongly localized in this well which leads to the increase of the electron density near the edges. Because of the presence of this well, the subbands start to depopulate from the central region of the wire and remain pinned in the well region until they are eventually pushed up by increasing magnetic field. We also demonstrate that the spin polarization of electron density as a function of magnetic field shows a pronounced double-loop pattern that can be related to the successive depopulation of the magnetosubbands. In contrast to the case of a smooth confinement, in hard-wall wires the compressible strips do not form in the vicinity of wire boundaries and spatial spin separation between spin-up and spin-down states near edges is absent.Comment: 9 pages, submitted to Phys. Rev.

Density Functional Theory for the Photoionization Dynamics of Uracil

Photoionization dynamics of the RNA base Uracil is studied in the framework of Density Functional Theory (DFT). The photoionization calculations take advantage of a newly developed parallel version of a multicentric approach to the calculation of the electronic continuum spectrum which uses a set of B-spline radial basis functions and a Kohn-Sham density functional hamiltonian. Both valence and core ionizations are considered. Scattering resonances in selected single-particle ionization channels are classified by the symmetry of the resonant state and the peak energy position in the photoelectron kinetic energy scale; the present results highlight once more the site specificity of core ionization processes. We further suggest that the resonant structures previously characterized in low-energy electron collision experiments are partly shifted below threshold by the photoionization processes. A critical evaluation of the theoretical results providing a guide for future experimental work on similar biosystems

High-density correlation energy expansion of the one-dimensional uniform electron gas

We show that the expression of the high-density (i.e small-$r_s$) correlation energy per electron for the one-dimensional uniform electron gas can be obtained by conventional perturbation theory and is of the form \Ec(r_s) = -\pi^2/360 + 0.00845 r_s + ..., where $r_s$ is the average radius of an electron. Combining these new results with the low-density correlation energy expansion, we propose a local-density approximation correlation functional, which deviates by a maximum of 0.1 millihartree compared to the benchmark DMC calculations.Comment: 7 pages, 2 figures, 3 tables, accepted for publication in J. Chem. Phy