Showing 1–46 of 46 results for author: Puska, M

Efficient method for calculating Raman spectra of solids with impurities and alloys and its application to two-dimensional transition metal dichalcogenides

Authors: Arsalan Hashemi, Arkady V. Krasheninnikov, Martti Puska, Hannu-Pekka Komsa

Abstract: Raman spectroscopy is a widely used, powerful, and nondestructive tool for studying the vibrational properties of bulk and low-dimensional materials. Raman spectra can be simulated using first-principles methods, but due to the high computational cost calculations are usually limited only to fairly small unit cells, which makes it difficult to carry out simulations for alloys and defects. Here, we… ▽ More Raman spectroscopy is a widely used, powerful, and nondestructive tool for studying the vibrational properties of bulk and low-dimensional materials. Raman spectra can be simulated using first-principles methods, but due to the high computational cost calculations are usually limited only to fairly small unit cells, which makes it difficult to carry out simulations for alloys and defects. Here, we develop an efficient method for simulating Raman spectra of alloys, benchmark it against full density-functional theory calculations, and apply it to several alloys of two-dimensional transition metal dichalcogenides. In this method, the Raman tensor for the supercell mode is constructed by summing up the Raman tensors of the pristine system weighted by the projections of the supercell vibrational modes to those of the pristine system. This approach is not limited to 2D materials and should be applicable to any crystalline solids with defects and impurities. To efficiently evaluate vibrational modes of very large supercells, we adopt mass approximation, although it is limited to chemically and structurally similar atomic substitutions. To benchmark our method, we first apply it to Mo$_x$W$_{(1-x)}$S$_2$ monolayer in the H-phase, where several experimental reports are available for comparison. Second, we consider Mo$_x$W$_{(1-x)}$Te$_2$ in the T'-phase, which has been proposed to be 2D topological insulator, but where experimental results for the monolayer alloy are still missing. We show that the projection scheme also provides a powerful tool for analyzing the origin of the alloy Raman-active modes in terms of the parent system eigenmodes. Finally, we examine the trends in characteristic Raman signatures for dilute concentrations of impurities in MoS$_2$. △ Less

Submitted 6 February, 2019; originally announced February 2019.

Journal ref: Phys. Rev. Materials 3, 023806 (2019)

Plasmon excitations in mixed metallic nanoarrays

Authors: Kevin Conley, Neha Nayyar, Tuomas P. Rossi, Mikael Kuisma, Volodymyr Turkowski, Martti Puska, Talat S. Rahman

Abstract: We study the plasmonic properties of arrays of atomic chains which comprise noble (Cu, Ag, and Au) and transition (Pd, Pt) metal atoms using time-dependent density-functional theory. We show that the response to the electromagnetic radiation is related to both physics, the geometry-dependent confinement of sp-valence electrons, and chemistry, the energy position of d-electrons in the different ato… ▽ More We study the plasmonic properties of arrays of atomic chains which comprise noble (Cu, Ag, and Au) and transition (Pd, Pt) metal atoms using time-dependent density-functional theory. We show that the response to the electromagnetic radiation is related to both physics, the geometry-dependent confinement of sp-valence electrons, and chemistry, the energy position of d-electrons in the different atomic species and the hybridization between d and sp electrons. As a result it is possible to tune the position of the surface plasmon resonance, split it to several peaks, and eventually achieve broadband absorption of radiation. Mixing the arrays with transition metals can strongly attenuate the plasmonic behaviour. We analyze the origin of these phenomena and show that they arise from rich interactions between single-particle electron-hole and collective electron excitations. The tunability of the plasmonic response of arrays of atomic chains, which can be realized on solid surfaces, opens wide possibilities for their applications. In the present study we obtain guidelines how the desired properties can be achieved. △ Less

Submitted 30 September, 2018; originally announced October 2018.

Structural details of Al/Al2O3 junctions and their role in formation of electron tunnel barriers

Authors: M. Koberidze, M. J. Puska, R. M. Nieminen

Abstract: We present a computational study of the adhesive and structural properties of the Al/Al2O3 interfaces as building blocks of the Metal-Insulator-Metal (MIM) tunnel devices, where electron transport is accomplished via tunnelling mechanism through the sandwiched insulating barrier. The main goal of this paper is to understand, on the atomic scale, the role of the geometrical details in the formation… ▽ More We present a computational study of the adhesive and structural properties of the Al/Al2O3 interfaces as building blocks of the Metal-Insulator-Metal (MIM) tunnel devices, where electron transport is accomplished via tunnelling mechanism through the sandwiched insulating barrier. The main goal of this paper is to understand, on the atomic scale, the role of the geometrical details in the formation of the tunnel barrier profiles. To provide reliable results, we carefully assess the accuracy of the traditional methods used to examine Al/Al2O3 interfaces. These are the most widely employed exchange-correlation functionals, LDA, PBE and PW91, the Universal Binding Energy Relation (UBER) for predicting equilibrium interfacial distances and adhesion energies, and the ideal work of separation as a measure of junction stability. Finally, we perform a detailed analysis of the atomic and interplanar relaxations in each junction. Our results imply that the structural irregularities on the surface of the Al film have a significant contribution to lowering the tunnel barrier height, while interplanar relaxations in the Al film, away from the immediate interface do not have a notable impact on the tunnelling properties. On the other hand, up to 5-7 layers of Al2O3 may be involved in shaping the tunnel barriers. Interplanar relaxations of these layers depend on the geometry of the interface and may result in the net contraction by 13% relative to the corresponding thickness in the bulk oxide. This is a significant amount as the tunnelling probability depends exponentially on the barrier width. △ Less

Submitted 1 December, 2017; originally announced December 2017.

Comments: 13 pages, 5 figures

Journal ref: Phys. Rev. B 97, 195406 (2018)

Electonic transport properties of nitrate-doped carbon nanotube networks

Authors: T. Ketolainen, V. Havu, E. Ö. Jónsson, M. J. Puska

Abstract: The conductivity of carbon nanotube (CNT) networks can be improved markedly by doping with nitric acid. In the present work, CNTs and junctions of CNTs functionalized with NO$_3$ molecules are investigated to understand the microscopic mechanism of nitric acid doping. According to our density functional theory band structure calculations, there is charge transfer from the CNT to adsorbed molecules… ▽ More The conductivity of carbon nanotube (CNT) networks can be improved markedly by doping with nitric acid. In the present work, CNTs and junctions of CNTs functionalized with NO$_3$ molecules are investigated to understand the microscopic mechanism of nitric acid doping. According to our density functional theory band structure calculations, there is charge transfer from the CNT to adsorbed molecules indicating p-type doping. The average doping efficiency of the NO$_3$ molecules is higher if the NO$_3$ molecules form complexes with water molecules. In addition to electron transport along individual CNTs, we have also studied electron transport between different types (metallic, semiconducting) of CNTs. Reflecting the differences in the electronic structures of semiconducting and metallic CNTs, we have found that besides turning semiconducting CNTs metallic, doping further increases electron transport most efficiently along semiconducting CNTs as well as through a junction between them. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 13 pages, 12 figures

Journal ref: Phys. Rev. Applied 9, 034010 (2018)

Kohn-Sham decomposition in real-time time-dependent density-functional theory: An efficient tool for analyzing plasmonic excitations

Authors: Tuomas P. Rossi, Mikael Kuisma, Martti J. Puska, Risto M. Nieminen, Paul Erhart

Abstract: The real-time-propagation formulation of time-dependent density-functional theory (RT-TDDFT) is an efficient method for modeling the optical response of molecules and nanoparticles. Compared to the widely adopted linear-response TDDFT approaches based on, e.g., the Casida equations, RT-TDDFT appears, however, lacking efficient analysis methods. This applies in particular to a decomposition of the… ▽ More The real-time-propagation formulation of time-dependent density-functional theory (RT-TDDFT) is an efficient method for modeling the optical response of molecules and nanoparticles. Compared to the widely adopted linear-response TDDFT approaches based on, e.g., the Casida equations, RT-TDDFT appears, however, lacking efficient analysis methods. This applies in particular to a decomposition of the response in the basis of the underlying single-electron states. In this work, we overcome this limitation by developing an analysis method for obtaining the Kohn-Sham electron-hole decomposition in RT-TDDFT. We demonstrate the equivalence between the developed method and the Casida approach by a benchmark on small benzene derivatives. Then, we use the method for analyzing the plasmonic response of icosahedral silver nanoparticles up to Ag$_{561}$. Based on the analysis, we conclude that in small nanoparticles individual single-electron transitions can split the plasmon into multiple resonances due to strong single-electron-plasmon coupling whereas in larger nanoparticles a distinct plasmon resonance is formed. △ Less

Submitted 8 March, 2017; originally announced March 2017.

Comments: 11 pages, 3 figures

Journal ref: Journal of Chemical Theory and Computation 13, 4779 (2017)

Theory and applications of generalized Pipek--Mezey Wannier functions

Authors: Elvar Ö. Jónsson, Susi Lehtola, Martti Puska, Hannes Jónsson

Abstract: The theory for the generation of Wannier functions within the generalized Pipek--Mezey approach [Lehtola, S.; Jónsson, H. J. Chem. Theory Comput. 2014, 10, 642] is presented and an implementation thereof is described. Results are presented for systems with periodicity in one, two and three dimensions as well as isolated molecules. The generalized Pipek--Mezey Wannier functions (PMWF) are highly lo… ▽ More The theory for the generation of Wannier functions within the generalized Pipek--Mezey approach [Lehtola, S.; Jónsson, H. J. Chem. Theory Comput. 2014, 10, 642] is presented and an implementation thereof is described. Results are presented for systems with periodicity in one, two and three dimensions as well as isolated molecules. The generalized Pipek--Mezey Wannier functions (PMWF) are highly localized orbitals consistent with chemical intuition where a distinction is maintained between σしぐま- and πぱい-orbitals. The PMWF method is compared with the so-called maximally localized Wannier functions (MLWF) that are frequently used for the analysis of condensed matter calculations. Whereas PMWFs maximize the localization criterion of Pipek and Mezey, MLWFs maximize that of Foster and Boys and have the disadvantage of mixing σしぐま- and πぱい-orbitals in many cases. The PMWF orbitals turn out to be as localized as the MLWF orbitals as evidenced by cross-comparison of the values of the PMWF and MLWF objective functions for the two types of orbitals. Our implementation in the atomic simulation environment (ASE) is compatible with various representations of the wave function, including real-space grids, plane waves and linear combinations of atomic orbitals. The projector augmented wave formalism for the representation of atomic core electrons is also supported. Results of calculations with the GPAW software are described here, but our implementation can also use output from other electronic structure software such as ABINIT, NWChem and VASP. △ Less

Submitted 14 February, 2017; v1 submitted 23 August, 2016; originally announced August 2016.

Comments: 17 pages, 5 figures

Journal ref: J. Chem. Theory Comput. 13, 460 (2017)

Effect of interface geometry on electron tunnelling in Al/Al$_2$O$_3$/Al junctions

Authors: M Koberidze, A V Feshchenko, M J Puska, R M Nieminen, J P Pekola

Abstract: We investigate how different interface geometries of an Al/Al$_2$O$_3$ junction, a common component of modern tunnel devices, affect electron transport through the tunnel barrier. We study six distinct Al/Al$_2$O$_3$ interfaces which differ in stacking sequences of the metal and the oxide surface atoms and the oxide termination. To construct model potential barrier profiles for each examined geome… ▽ More We investigate how different interface geometries of an Al/Al$_2$O$_3$ junction, a common component of modern tunnel devices, affect electron transport through the tunnel barrier. We study six distinct Al/Al$_2$O$_3$ interfaces which differ in stacking sequences of the metal and the oxide surface atoms and the oxide termination. To construct model potential barrier profiles for each examined geometry, we rely on first-principles density-functional theory (DFT) calculations for the barrier heights and the shapes of the interface regions as well as on experimental data for the barrier widths. We show that even tiny variations in the atomic arrangement at the interface cause significant changes in the tunnel barrier parameters and, consequently, in electron transport properties. Especially, we obtain that variations in the crucial barrier heights and widths can be as large as 2 eV and 5 Å, respectively. Finally, to gain information about the average properties of the measured junction, we fit the conductance calculated within the WKB approximation to the experimental data and interpret the fit parameters with the help of the DFT results. △ Less

Submitted 27 November, 2015; originally announced November 2015.

Comments: 9 pages, 8 figures

Journal ref: Journal of Physics D: Applied Physics 49, 165303 (2016)

Nanoplasmonics simulations at the basis set limit through completeness-optimized, local numerical basis sets

Authors: Tuomas P. Rossi, Susi Lehtola, Arto Sakko, Martti J. Puska, Risto M. Nieminen

Abstract: We present an approach for generating local numerical basis sets of improving accuracy for first-principles nanoplasmonics simulations within time-dependent density functional theory. The method is demonstrated for copper, silver, and gold nanoparticles that are of experimental interest but computationally demanding due to the semi-core d-electrons that affect their plasmonic response. The basis s… ▽ More We present an approach for generating local numerical basis sets of improving accuracy for first-principles nanoplasmonics simulations within time-dependent density functional theory. The method is demonstrated for copper, silver, and gold nanoparticles that are of experimental interest but computationally demanding due to the semi-core d-electrons that affect their plasmonic response. The basis sets are constructed by augmenting numerical atomic orbital basis sets by truncated Gaussian-type orbitals generated by the completeness-optimization scheme, which is applied to the photoabsorption spectra of homoatomic metal atom dimers. We obtain basis sets of improving accuracy up to the complete basis set limit and demonstrate that the performance of the basis sets transfers to simulations of larger nanoparticles and nanoalloys as well as to calculations with various exchange-correlation functionals. This work promotes the use of the local basis set approach of controllable accuracy in first-principles nanoplasmonics simulations and beyond. △ Less

Submitted 3 September, 2015; originally announced September 2015.

Comments: 11 pages, 6 figures

Journal ref: J. Chem. Phys. 142, 094114 (2015)

Quantized evolution of the plasmonic response in a stretched nanorod

Authors: Tuomas P. Rossi, Asier Zugarramurdi, Martti J. Puska, Risto M. Nieminen

Abstract: Quantum aspects, such as electron tunneling between closely separated metallic nanoparticles, are crucial for understanding the plasmonic response of nanoscale systems. We explore quantum effects on the response of the conductively coupled metallic nanoparticle dimer. This is realized by stretching a nanorod, which leads to the formation of a narrowing atomic contact between the two nanorod ends.… ▽ More Quantum aspects, such as electron tunneling between closely separated metallic nanoparticles, are crucial for understanding the plasmonic response of nanoscale systems. We explore quantum effects on the response of the conductively coupled metallic nanoparticle dimer. This is realized by stretching a nanorod, which leads to the formation of a narrowing atomic contact between the two nanorod ends. Based on first-principles time-dependent density-functional-theory calculations, we find a discontinuous evolution of the plasmonic response as the nanorod is stretched. This is especially pronounced for the intensity of the main charge-transfer plasmon mode. We show the correlation between the observed discontinuities and the discrete nature of the conduction channels supported by the formed atomic-sized junction. △ Less

Submitted 7 December, 2015; v1 submitted 3 September, 2015; originally announced September 2015.

Comments: Main text: 6 pages, 2 figures; Supplemental Material: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 115, 236804 (2015)

Full correlation single-particle positron potentials for a positron and a positronium interacting with atoms

Authors: A. Zubiaga, F. Tuomisto, M. Puska

Abstract: In this work we define single-particle potentials for a positron and a positronium atom interacting with light atoms (H, He, Li and Be) by inverting a single-particle Schrödinger equation. For this purpose, we use accurate energies and positron densities obtained from the many-body wavefunction of the corresponding positronic systems. The introduced potentials describe the exact correlations for t… ▽ More In this work we define single-particle potentials for a positron and a positronium atom interacting with light atoms (H, He, Li and Be) by inverting a single-particle Schrödinger equation. For this purpose, we use accurate energies and positron densities obtained from the many-body wavefunction of the corresponding positronic systems. The introduced potentials describe the exact correlations for the calculated systems including the formation of a positronium atom. We show that the scattering lengths and the low-energy s-wave phase shifts from accurate many-body calculations are well accounted for by the introduced potential. We also calculate self-consistent two-component density-functional theory positron potentials and densities for the bound positronic systems within the local density approximation. They are in a very good agreement with the many-body results, provided that the finite-positron-density electron-positron correlation potential is used, and they can also describe systems comprising a positronium atom. We argue that the introduced single-particle positron potentials defined for single molecules are transferable to the condensed phase when the inter-molecular interactions are weak. When this condition is fulfilled, the total positron potential can be constructed in a good approximation as the superposition of the molecular potentials. △ Less

Submitted 13 May, 2014; v1 submitted 29 May, 2013; originally announced May 2013.

Comments: 10 pages, 2 tables, 5 figures

MSC Class: 81V55

Journal ref: Phys. Rev. A 89, 052707 (2014)

Ab-initio transport fingerprints for resonant scattering in graphene

Authors: Karri Saloriutta, Andreas Uppstu, Ari Harju, Martti J. Puska

Abstract: We have recently shown that by using a scaling approach for randomly distributed topological defects in graphene, reliable estimates for transmission properties of macroscopic samples can be calculated based even on single-defect calculations [A. Uppstu et al., Phys. Rev. B 85, 041401 (2012)]. We now extend this approach of energy-dependent scattering cross sections to the case of adsorbates on gr… ▽ More We have recently shown that by using a scaling approach for randomly distributed topological defects in graphene, reliable estimates for transmission properties of macroscopic samples can be calculated based even on single-defect calculations [A. Uppstu et al., Phys. Rev. B 85, 041401 (2012)]. We now extend this approach of energy-dependent scattering cross sections to the case of adsorbates on graphene by studying hydrogen and carbon adatoms as well as epoxide and hydroxyl groups. We show that a qualitative understanding of resonant scattering can be gained through density functional theory results for a single-defect system, providing a transmission "fingerprint" characterizing each adsorbate type. This information can be used to reliably predict the elastic mean free path for moderate defect densities directly using ab-initio methods. We present tight-binding parameters for carbon and epoxide adsorbates, obtained to match the density-functional theory based scattering cross sections. △ Less

Submitted 30 November, 2012; originally announced November 2012.

Journal ref: Phys. Rev. B 86, 235417 (2012)

Matter-positronium interaction: An exact diagonalization study of the He atom - positronium system

Authors: A. Zubiaga, F. Tuomisto, M. J. Puska

Abstract: The many-body system comprising a He nucleus, three electrons, and a positron has been studied using the exact diagonalization technique. The purpose has been to clarify to which extent the system can be considered as a distinguishable positronium (Ps) atom interacting with a He atom and, thereby, to pave the way to a practical atomistic modeling of Ps states and annihilation in matter. The maximu… ▽ More The many-body system comprising a He nucleus, three electrons, and a positron has been studied using the exact diagonalization technique. The purpose has been to clarify to which extent the system can be considered as a distinguishable positronium (Ps) atom interacting with a He atom and, thereby, to pave the way to a practical atomistic modeling of Ps states and annihilation in matter. The maximum value of the distance between the positron and the nucleus is constrained and the Ps atom at different distances from the nucleus is identified from the electron and positron densities, as well as from the electron-positron distance and center-of-mass distributions. The polarization of the Ps atom increases as its distance from the nucleus decreases. A depletion of the He electron density, particularly large at low density values, has been observed. The ortho-Ps pick-off annihilation rate calculated as the overlap of the positron and the free He electron densities has to be corrected for the observed depletion, specially at large pores/voids. △ Less

Submitted 29 February, 2012; originally announced February 2012.

Comments: 18 pages, 8 figures

Journal ref: Phys. Rev. A 85, 052707 (2012)

Electronic transport in graphene-based structures: an effective cross section approach

Authors: Andreas Uppstu, Karri Saloriutta, Ari Harju, Martti Puska, Antti-Pekka Jauho

Abstract: We show that transport in low-dimensional carbon structures with finite concentrations of scatterers can be modeled by utilising scaling theory and effective cross sections. Our reults are based on large scale numerical simulations of carbon nanotubes and graphene nanoribbons, using a tightbinding model with parameters obtained from first principles electronic structure calculations. As shown by a… ▽ More We show that transport in low-dimensional carbon structures with finite concentrations of scatterers can be modeled by utilising scaling theory and effective cross sections. Our reults are based on large scale numerical simulations of carbon nanotubes and graphene nanoribbons, using a tightbinding model with parameters obtained from first principles electronic structure calculations. As shown by a comprehensive statistical analysis, the scattering cross sections can be used to estimate the conductance of a quasi-1D system both in the Ohmic and localized regimes. They can be computed with good accuracy from the transmission functions of single defects, greatly reducing the computational cost and paving the way towards using first principles methods to evaluate the conductance of mesoscopic systems, consisting of millions of atoms. △ Less

Submitted 1 November, 2011; originally announced November 2011.

Comments: Submitted to Physical Review Letters

Journal ref: Phys. Rev. B 85, 041401(R) (2012)

Nonadiabatic Ehrenfest molecular dynamics within the projector augmented-wave method

Authors: Ari Ojanperä, Ville Havu, Lauri Lehtovaara, Martti Puska

Abstract: We have derived equations for nonadiabatic Ehrenfest molecular dynamics which conserve the total energy in the case of time-dependent discretization for electrons. A discretization is time-dependent in all cases where it or part of it depends on the positions of the nuclei, for example, in atomic orbital basis sets, and in the projector augmented-wave (PAW) method, where the augmentation functions… ▽ More We have derived equations for nonadiabatic Ehrenfest molecular dynamics which conserve the total energy in the case of time-dependent discretization for electrons. A discretization is time-dependent in all cases where it or part of it depends on the positions of the nuclei, for example, in atomic orbital basis sets, and in the projector augmented-wave (PAW) method, where the augmentation functions depend on the nuclear positions. We have derived, implemented, and analyzed the energy conserving equations and their most common approximations for a 1D test system where we can achieve numerical results converged to a high accuracy. Based on the observations in 1D, we implement and analyze the Ehrenfest molecular dynamics in 3D using the PAW method and the time-dependent density functional formalism. We demonstrate the applicability of our method by carrying out calculations for small and medium sized molecules in both the adiabatic and the nonadiabatic regime. △ Less

Submitted 28 September, 2011; originally announced September 2011.

Comments: 12 pages, 10 figures

Local Semiconducting Transition in Armchair Carbon Nanotubes: The Effect of Periodic Bi-site Perturbation on Electronic and Transport Properties of Carbon Nanotubes

Authors: M. J. Hashemi, K. Sääskilahti, M. J. Puska

Abstract: In carbon nanotubes, the most abundant defects, caused for example by irradiation or chemisorption treatments, are small perturbing clusters, i.e. bi-site defects, extending over both A and B sites. The relative positions of these perturbing clusters play a crucial role in determining the electronic properties of carbon nanotubes. Using bandstructure and electronic transport calculations, we find… ▽ More In carbon nanotubes, the most abundant defects, caused for example by irradiation or chemisorption treatments, are small perturbing clusters, i.e. bi-site defects, extending over both A and B sites. The relative positions of these perturbing clusters play a crucial role in determining the electronic properties of carbon nanotubes. Using bandstructure and electronic transport calculations, we find out that in the case of armchair metallic nanotubes a band gap opens up when the clusters fulfill a certain periodicity condition. This phenomenon might be used in future nanoelectronic devices in which certain regions of single metallic nanotubes could be turned to semiconducting ones. Although in this work we study specifically the effect of hydrogen adatom clusters, the phenomenon is general for different types of defects. Moreover, we study the influence of the length and randomness of the defected region on the electron transport through it. △ Less

Submitted 3 January, 2011; originally announced January 2011.

Comments: 5 Pages, 5 Figures

Adsorption structures of phenol on the Si(001)-(2 \times 1) surface calculated using density functional theory

Authors: Karen Johnston, Andris Gulans, Tuukka Verho, Martti J. Puska

Abstract: Several dissociated and two non-dissociated adsorption structures of the phenol molecule on the Si(001)-(2 \times 1) surface are studied using density functional theory with various exchange and correlation functionals. The relaxed structures and adsorption energies are obtained and it is found that the dissociated structures are energetically more favourable than the non-dissociated structures. H… ▽ More Several dissociated and two non-dissociated adsorption structures of the phenol molecule on the Si(001)-(2 \times 1) surface are studied using density functional theory with various exchange and correlation functionals. The relaxed structures and adsorption energies are obtained and it is found that the dissociated structures are energetically more favourable than the non-dissociated structures. However, the ground state energies alone do not determine which structure is obtained experimentally. To elucidate the situation core level shift spectra for Si 2p and C 1s states are simulated and compared with experimentally measured spectra. Several transition barriers were calculated in order to determine which adsorption structures are kinetically accessible. Based on these results we conclude that the molecule undergoes the dissociation of two hydrogen atoms on adsorption. △ Less

Submitted 22 June, 2010; originally announced June 2010.

Comments: 11 pages, 6 figures

Journal ref: Phys. Rev. B 81 235428 (2010)

Tunability of the optical absorption in small silver cluster-polymer hybrid systems

Authors: Laura Koponen, Lasse O. Tunturivuori, Martti J. Puska, Y. Hancock

Abstract: We have calculated the absorption characteristics of different hybrid systems consisting of Ag, Ag2 or Ag3 atomic clusters and poly(methacrylic acid) (PMAA) using the time-dependent density-functional theory. The polymer is found to have an extensive structural-dependency on the spectral patterns of the hybrid systems relative to the bare clusters. The absorption spectrum can be `tuned' to the vis… ▽ More We have calculated the absorption characteristics of different hybrid systems consisting of Ag, Ag2 or Ag3 atomic clusters and poly(methacrylic acid) (PMAA) using the time-dependent density-functional theory. The polymer is found to have an extensive structural-dependency on the spectral patterns of the hybrid systems relative to the bare clusters. The absorption spectrum can be `tuned' to the visible range for hybrid systems with an odd number of electrons per silver cluster, whereas for hybrid systems comprising an even number of electrons, the leading absorption edge can be shifted up to about 4.5 eV. The results give theoretical support to the experimental observations on the absorption in the visible range in metal cluster-polymer hybrid structures. △ Less

Submitted 16 April, 2010; v1 submitted 10 March, 2010; originally announced March 2010.

Comments: Updated layout and minor changes in versions 2 and 3

All-electron density functional theory and time-dependent density functional theory with high-order finite elements

Authors: Lauri Lehtovaara, Ville Havu, Martti Puska

Abstract: We present for static density functional theory and time-dependent density functional theory calculations an all-electron method which employs high-order hierarchical finite element bases. Our mesh generation scheme, in which structured atomic meshes are merged to an unstructured molecular mesh, allows a highly nonuniform discretization of the space. Thus it is possible to represent the core and… ▽ More We present for static density functional theory and time-dependent density functional theory calculations an all-electron method which employs high-order hierarchical finite element bases. Our mesh generation scheme, in which structured atomic meshes are merged to an unstructured molecular mesh, allows a highly nonuniform discretization of the space. Thus it is possible to represent the core and valence states using the same discretization scheme, i.e., no pseudopotentials or similar treatments are required. The nonuniform discretization also allows the use of large simulation cells, and therefore avoids any boundary effects. △ Less

Submitted 6 August, 2009; v1 submitted 13 February, 2009; originally announced February 2009.

Comments: 11 pages, 9 figures; final (=published) version

Journal ref: J. Chem. Phys. 131, 054103 (2009)

Energetics of positron states trapped at vacancies in solids

Authors: I. Makkonen, M. J. Puska

Abstract: We report a computational first-principles study of positron trapping at vacancy defects in metals and semiconductors. The main emphasis is on the energetics of the trapping process including the interplay between the positron state and the defect's ionic structure and on the ensuing annihilation characteristics of the trapped state. For vacancies in covalent semiconductors the ion relaxation is… ▽ More We report a computational first-principles study of positron trapping at vacancy defects in metals and semiconductors. The main emphasis is on the energetics of the trapping process including the interplay between the positron state and the defect's ionic structure and on the ensuing annihilation characteristics of the trapped state. For vacancies in covalent semiconductors the ion relaxation is a crucial part of the positron trapping process enabling the localization of the positron state. However, positron trapping does not strongly affect the characteristic features of the electronic structure, e.g., the ionization levels change only moderately. Also in the case of metal vacancies the positron-induced ion relaxation has a noticeable effect on the calculated positron lifetime and momentum distribution of annihilating electron-positron pairs. △ Less

Submitted 7 July, 2007; v1 submitted 17 April, 2007; originally announced April 2007.

Comments: Submitted to Physical Review B on 17 April 2007. Revised version submitted on 6 July 2007

Journal ref: Physical Review B 76, 054119 (2007).

Analysis of electron-positron momentum spectra of metallic alloys as supported by first-principles calculations

Authors: P. Folegati, I. Makkonen, R. Ferragut, M. J. Puska

Abstract: Electron-positron momentum distributions measured by the coincidence Doppler broadening method can be used in the chemical analysis of the annihilation environment, typically a vacancy-impurity complex in a solid. In the present work, we study possibilities for a quantitative analysis, i.e., for distinguishing the average numbers of different atomic species around the defect. First-principles el… ▽ More Electron-positron momentum distributions measured by the coincidence Doppler broadening method can be used in the chemical analysis of the annihilation environment, typically a vacancy-impurity complex in a solid. In the present work, we study possibilities for a quantitative analysis, i.e., for distinguishing the average numbers of different atomic species around the defect. First-principles electronic structure calculations self-consistently determining electron and positron densities and ion positions are performed for vacancy-solute complexes in Al-Cu, Al-Mg-Cu, and Al-Mg-Cu-Ag alloys. The ensuing simulated coincidence Doppler broadening spectra are compared with measured ones for defect identification. A linear fitting procedure, which uses the spectra for positrons trapped at vacancies in pure constituent metals as components, has previously been employed to find the relative percentages of different atomic species around the vacancy [A. Somoza et al. Phys. Rev. B 65, 094107 (2002)]. We test the reliability of the procedure by the help of first-principles results for vacancy-solute complexes and vacancies in constituent metals. △ Less

Submitted 15 February, 2007; v1 submitted 19 September, 2006; originally announced September 2006.

Comments: Submitted to Physical Review B on September 19 2006. Revised version submitted on November 8 2006. Published on February 14 2007

Journal ref: Physical Review B 75, 054201 (2007) (10 pages)

Bound states in two-dimensional shielded potentials

Authors: I. Nagy, M. J. Puska, N. Zabala

Abstract: We study numerically the existence and character of bound states for positive and negative point charges shielded by the response of a two-dimensional homogeneous electron gas. The problem is related to many physical situations and has recently arisen in experiments for impurities on metal surfaces with Schockley surface states. Mathematical theorems ascertain a bound state for two-dimensional c… ▽ More We study numerically the existence and character of bound states for positive and negative point charges shielded by the response of a two-dimensional homogeneous electron gas. The problem is related to many physical situations and has recently arisen in experiments for impurities on metal surfaces with Schockley surface states. Mathematical theorems ascertain a bound state for two-dimensional circularly symmetric potentials $V(r)$ with $\int_{0}^{\infty}{dr r V(r)} \leq 0$. We find that a shielded potential with $\int_{0}^{\infty}{dr r V(r)} > 0$ may also sustain a bound state. Moreover, on the same footing we study the electron-electron interactions in the two-dimensional electron gas finding a bound state with an energy minimum for a certain electron gas density. △ Less

Submitted 17 March, 2006; originally announced March 2006.

Comments: Revtex-file and 8 eps figures

Three real-space discretization techniques in electronic structure calculations

Authors: T. Torsti, T. Eirola, J. Enkovaara, T. Hakala, P. Havu, V. Havu, T. Höynälänmaa, J. Ignatius, M. Lyly, I. Makkonen, T. T. Rantala, J. Ruokolainen, K. Ruotsalainen, E. Räsänen, H. Saarikoski, M. J. Puska

Abstract: A characteristic feature of the state-of-the-art of real-space methods in electronic structure calculations is the diversity of the techniques used in the discretization of the relevant partial differential equations. In this context, the main approaches include finite-difference methods, various types of finite-elements and wavelets. This paper reports on the results of several code development… ▽ More A characteristic feature of the state-of-the-art of real-space methods in electronic structure calculations is the diversity of the techniques used in the discretization of the relevant partial differential equations. In this context, the main approaches include finite-difference methods, various types of finite-elements and wavelets. This paper reports on the results of several code development projects that approach problems related to the electronic structure using these three different discretization methods. We review the ideas behind these methods, give examples of their applications, and discuss their similarities and differences. △ Less

Submitted 24 January, 2006; v1 submitted 10 January, 2006; originally announced January 2006.

Comments: 39 pages, 10 figures, accepted to a special issue of "physica status solidi (b) - basic solid state physics" devoted to the CECAM workshop "State of the art developments and perspectives of real-space electronic structure techniques in condensed matter and molecular physics". v2: Minor stylistic and typographical changes, partly inspired by referee comments

Giant vortices in rotating electron droplets

Authors: E. Rasanen, H. Saarikoski, Y. Yu, A. Harju, M. J. Puska, S. M. Reimann

Abstract: We predict the formation of giant vortices in quasi-two-dimensional quantum dots at high magnetic fields, i.e., in rapidly rotating electron droplets. Our numerical results of quantum dots confined by a flat, anharmonic potential show ground states where vortices are accumulated in the center of the dot, thereby leading to large cores in the electron and current densities. The phenomenon is anal… ▽ More We predict the formation of giant vortices in quasi-two-dimensional quantum dots at high magnetic fields, i.e., in rapidly rotating electron droplets. Our numerical results of quantum dots confined by a flat, anharmonic potential show ground states where vortices are accumulated in the center of the dot, thereby leading to large cores in the electron and current densities. The phenomenon is analogous to what was recently found in rotating Bose-Einstein condensates. The giant-vortex states leave measurable signatures in the ground-state energetics. The conditions for the giant-vortex formation as well as the internal structure of the vortex cores are discussed. △ Less

Submitted 7 May, 2006; v1 submitted 26 September, 2005; originally announced September 2005.

Comments: to appear in Phys. Rev. B

Journal ref: Phys. Rev. B 73, 235324 (2006)

Modeling the momentum distributions of annihilating electron-positron pairs in solids

Authors: I. Makkonen, M. Hakala, M. J. Puska

Abstract: Measuring the Doppler broadening of the positron annihilation radiation or the angular correlation between the two annihilation gamma quanta reflects the momentum distribution of electrons seen by positrons in the material.Vacancy-type defects in solids localize positrons and the measured spectra are sensitive to the detailed chemical and geometric environments of the defects. However, the measu… ▽ More Measuring the Doppler broadening of the positron annihilation radiation or the angular correlation between the two annihilation gamma quanta reflects the momentum distribution of electrons seen by positrons in the material.Vacancy-type defects in solids localize positrons and the measured spectra are sensitive to the detailed chemical and geometric environments of the defects. However, the measured information is indirect and when using it in defect identification comparisons with theoretically predicted spectra is indispensable. In this article we present a computational scheme for calculating momentum distributions of electron-positron pairs annihilating in solids. Valence electron states and their interaction with ion cores are described using the all-electron projector augmented-wave method, and atomic orbitals are used to describe the core states. We apply our numerical scheme to selected systems and compare three different enhancement (electron-positron correlation) schemes previously used in the calculation of momentum distributions of annihilating electron-positron pairs within the density-functional theory. We show that the use of a state-dependent enhancement scheme leads to better results than a position-dependent enhancement factor in the case of ratios of Doppler spectra between different systems. Further, we demonstrate the applicability of our scheme for studying vacancy-type defects in metals and semiconductors. Especially we study the effect of forces due to a positron localized at a vacancy-type defect on the ionic relaxations. △ Less

Submitted 14 November, 2005; v1 submitted 1 September, 2005; originally announced September 2005.

Comments: Submitted to Physical Review B on September 1 2005. Revised manuscript submitted on November 14 2005

Journal ref: Physical Review B 73, 035103 (2006)

Spin dependent electron transport through a magnetic resonant tunneling diode

Authors: Paula Havu, Noora Tuomisto, Riikka Vaananen, Martti J. Puska, Risto M. Nieminen

Abstract: Electron transport properties in nanostructures can be modeled, for example, by using the semiclassical Wigner formalism or the quantum mechanical Green's functions formalism. We compare the performance and the results of these methods in the case of magnetic resonant-tunneling diodes. We have implemented the two methods within the self-consistent spin-density-functional theory. Our numerical im… ▽ More Electron transport properties in nanostructures can be modeled, for example, by using the semiclassical Wigner formalism or the quantum mechanical Green's functions formalism. We compare the performance and the results of these methods in the case of magnetic resonant-tunneling diodes. We have implemented the two methods within the self-consistent spin-density-functional theory. Our numerical implementation of the Wigner formalism is based on the finite-difference scheme whereas for the Green's function formalism the finite-element method is used. As a specific application, we consider the device studied by Slobodskyy et all. [Phys. Rev. Lett. 90, 246601 (2003)] and analyze their experimental results. The Wigner and Green's functions formalisms give similar electron densities and potentials but, surprisingly, the former method requires much more computer resources in order to obtain numerically accurate results for currents. Both of the formalisms can successfully be used to model magnetic resonant tunneling diode structures. △ Less

Submitted 24 November, 2004; originally announced November 2004.

Comments: 13 pages and 12 figures

Calculation of valence electron momentum densities using the projector augmented-wave method

Authors: I. Makkonen, M. Hakala, M. J. Puska

Abstract: We present valence electron Compton profiles calculated within the density-functional theory using the all-electron full-potential projector augmented-wave method (PAW). Our results for covalent (Si), metallic (Li, Al) and hydrogen-bonded ((H_2O)_2) systems agree well with experiments and computational results obtained with other band-structure and basis set schemes. The PAW basis set describes… ▽ More We present valence electron Compton profiles calculated within the density-functional theory using the all-electron full-potential projector augmented-wave method (PAW). Our results for covalent (Si), metallic (Li, Al) and hydrogen-bonded ((H_2O)_2) systems agree well with experiments and computational results obtained with other band-structure and basis set schemes. The PAW basis set describes the high-momentum Fourier components of the valence wave functions accurately when compared with other basis set schemes and previous all-electron calculations. △ Less

Submitted 13 December, 2004; v1 submitted 17 September, 2004; originally announced September 2004.

Comments: Submitted to Journal of Physics and Chemistry of Solids on September 17 2004. Revised version submitted on December 13 2004

Journal ref: Journal of Physics and Chemistry of Solids 66, 1128 (2005).

Stability of vortex structures in quantum dots

Authors: H. Saarikoski, S. M. Reimann, E. Rasanen, A. Harju, M. Puska

Abstract: We study the stability and structure of vortices emerging in two-dimensional quantum dots in high magnetic fields. Our results obtained with exact diagonalization and density-functional calculations show that vortex structures can be found in various confining potentials. In non-symmetric external potentials we find off-electron vortices that are localized giving rise to charge deficiency or hol… ▽ More We study the stability and structure of vortices emerging in two-dimensional quantum dots in high magnetic fields. Our results obtained with exact diagonalization and density-functional calculations show that vortex structures can be found in various confining potentials. In non-symmetric external potentials we find off-electron vortices that are localized giving rise to charge deficiency or holes in the electron density with rotating currents around them. We discuss the role of quantum fluctuations and show that vortex formation is observable in the energetics of the system. Our findings suggest that vortices can be used to characterize the solutions in high magnetic fields, giving insight into the underlying internal structure of the electronic wave function. △ Less

Submitted 20 August, 2004; originally announced August 2004.

Comments: 8 pages, 7 figures

Self-consistent study of electron confinement to metallic thin films on solid surfaces

Authors: E. Ogando, N. Zabala, E. V. Chulkov, M. J. Puska

Abstract: We present a method for density-functional modeling of metallic overlayers grown on metallic supports. It offers a tool to study nanostructures and combines the power of self-consistent pseudopotential calculations with the simplicity of a one-dimensional approach. The model is applied to Pb layers grown on the Cu(111) surface. More specifically, Pb is modeled as stabilized jellium and the Cu(11… ▽ More We present a method for density-functional modeling of metallic overlayers grown on metallic supports. It offers a tool to study nanostructures and combines the power of self-consistent pseudopotential calculations with the simplicity of a one-dimensional approach. The model is applied to Pb layers grown on the Cu(111) surface. More specifically, Pb is modeled as stabilized jellium and the Cu(111) substrate is represented by a one-dimensional pseudopotential that reproduces experimental positions of both the Cu Fermi level and the energy gap of the band structure projected along the (111) direction. The model is used to study the quantum well states in the Pb overlayer. Their analysis gives the strength of the electron confinement barriers at the interface and at the surface facing the vacuum. Our results and analysis support the interpretation of the quantum well state spectra measured by the scanning tunneling spectroscopy. △ Less

Submitted 6 July, 2004; originally announced July 2004.

Comments: 12 pages, 8 figures and 1 table. Submitted to PRB

Electron transport through quantum wires and point contacts

Authors: Paula Havu, Martti Puska, Risto Nieminen, Ville Havu

Abstract: We have studied quantum wires using the Green's function technique and the density-functional theory, calculating the electronic structure and the conductance. All the numerics are implemented using the finite-element method with a high-order polynomial basis. For short wires, i.e. quantum point contacts, the zero-bias conductance shows, as a function of the gate voltage and at a finite temperat… ▽ More We have studied quantum wires using the Green's function technique and the density-functional theory, calculating the electronic structure and the conductance. All the numerics are implemented using the finite-element method with a high-order polynomial basis. For short wires, i.e. quantum point contacts, the zero-bias conductance shows, as a function of the gate voltage and at a finite temperature, a plateau at around 0.7G_0. (G_0 = 2e^2/h is the quantum conductance). The behavior, which is caused in our mean-field model by spontaneous spin polarization in the constriction, is reminiscent of the so-called 0.7-anomaly observed in experiments. In our model the temperature and the wire length affect the conductance-gate voltage curves in the same way as in the measured data. △ Less

Submitted 7 June, 2004; originally announced June 2004.

Comments: 8 pages

Vortex formation in quantum dots in high magnetic fields

Authors: H. Saarikoski, A. Harju, M. J. Puska, R. M. Nieminen

Abstract: We study electronic structures of two-dimensional quantum dots in high magnetic fields using the density-functional theory (DFT) and the exact diagonalization (ED). With increasing magnetic field, beyond the formation of the totally spin-polarized maximum density droplet (MDD) state, the DFT gives the ground-state total angular momentum as a continuous function with well-defined plateaus. The pl… ▽ More We study electronic structures of two-dimensional quantum dots in high magnetic fields using the density-functional theory (DFT) and the exact diagonalization (ED). With increasing magnetic field, beyond the formation of the totally spin-polarized maximum density droplet (MDD) state, the DFT gives the ground-state total angular momentum as a continuous function with well-defined plateaus. The plateaus agree well with the magic angular momenta of the ED calculation. By constructing a conditional wave function from the Kohn-Sham states we show that vortices enter the quantum dot one-by-one at the transition to the state with the adjacent magic angular momentum. Vortices are also observed outside the high-density region of the quantum dot. These findings are compared to the ED results and we report a significant agreement. We study also interpretations and limitations of the density functional approach in these calculations. △ Less

Submitted 29 April, 2004; originally announced April 2004.

Comments: 4 pages, 5 figures, QD2004 conference paper

Journal ref: Physica E 26, 317 (2005)

Impurity effects in quantum dots: Towards quantitative modeling

Authors: E. Rasanen, J. Konemann, R. J. Haug, M. J. Puska, R. M. Nieminen

Abstract: We have studied the single-electron transport spectrum of a quantum dot in GaAs/AlGaAs resonant tunneling device. The measured spectrum has irregularities indicating a broken circular symmetry. We model the system with an external potential consisting of a parabolic confinement and a negatively charged Coulombic impurity placed in the vicinity of the quantum dot. The model leads to a good agreem… ▽ More We have studied the single-electron transport spectrum of a quantum dot in GaAs/AlGaAs resonant tunneling device. The measured spectrum has irregularities indicating a broken circular symmetry. We model the system with an external potential consisting of a parabolic confinement and a negatively charged Coulombic impurity placed in the vicinity of the quantum dot. The model leads to a good agreement between the calculated single-electron eigenenergies and the experimental spectrum. Furthermore, we use the spin-density-functional theory to study the energies and angular momenta when the system contains many interacting electrons. In the high magnetic field regime the increasing electron number is shown to reduce the distortion induced by the impurity. △ Less

Submitted 23 April, 2004; originally announced April 2004.

Comments: 6 pages, 8 figures, submitted to PRB

Journal ref: Phys. Rev. B 70, 115308 (2004)

Vortex Clusters in Quantum Dots

Authors: H. Saarikoski, A. Harju, M. J. Puska, R. M. Nieminen

Abstract: We study electronic structures of two-dimensional quantum dots in strong magnetic fields using mean-field density-functional theory and exact diagonalization. Our numerically accurate mean-field solutions show a reconstruction of the uniform-density electron droplet when the magnetic field flux quanta enter one by one the dot in stronger fields. These quanta correspond to repelling vortices form… ▽ More We study electronic structures of two-dimensional quantum dots in strong magnetic fields using mean-field density-functional theory and exact diagonalization. Our numerically accurate mean-field solutions show a reconstruction of the uniform-density electron droplet when the magnetic field flux quanta enter one by one the dot in stronger fields. These quanta correspond to repelling vortices forming polygonal clusters inside the dot. We find similar structures in the exact treatment of the problem using a new type of operator for the analysis. We discuss important differences and limitations of the methods used. △ Less

Submitted 24 May, 2004; v1 submitted 20 February, 2004; originally announced February 2004.

Comments: 4 pages, 4 figures

Journal ref: Phys. Rev. Lett. 93, 116802 (2004)

Rectangular quantum dots in high magnetic fields

Authors: E. Räsänen, A. Harju, M. J. Puska, R. M. Nieminen

Abstract: We use density-functional methods to study the effects of an external magnetic field on two-dimensional quantum dots with a rectangular hard-wall confining potential. The increasing magnetic field leads to spin polarization and formation of a highly inhomogeneous maximum-density droplet at the predicted magnetic field strength. At higher fields, we find an oscillating behavior in the electron de… ▽ More We use density-functional methods to study the effects of an external magnetic field on two-dimensional quantum dots with a rectangular hard-wall confining potential. The increasing magnetic field leads to spin polarization and formation of a highly inhomogeneous maximum-density droplet at the predicted magnetic field strength. At higher fields, we find an oscillating behavior in the electron density and in the magnetization of the dot. We identify a rich variety of phenomena behind the periodicity and analyze the complicated many-electron dynamics, which is shown to be highly dependent on the shape of the quantum dot. △ Less

Submitted 28 November, 2003; originally announced November 2003.

Comments: 6 pages, 6 figures, submitted to Phys. Rev. B

Journal ref: Phys. Rev. B 69, 165309 (2004)

Quantum size effects in Pb islands on Cu(111): Electronic-structure calculations

Authors: E. Ogando, N. Zabala, E. V. Chulkov, M. J. Puska

Abstract: The appearance of "magic" heights of Pb islands grown on Cu(111) is studied by self-consistent electronic structure calculations. The Cu(111) substrate is modeled with a one-dimensional pseudopotential reproducing the essential features, i.e. the band gap and the work function, of the Cu band structure in the [111] direction. Pb islands are presented as stabilized jellium overlayers. The experim… ▽ More The appearance of "magic" heights of Pb islands grown on Cu(111) is studied by self-consistent electronic structure calculations. The Cu(111) substrate is modeled with a one-dimensional pseudopotential reproducing the essential features, i.e. the band gap and the work function, of the Cu band structure in the [111] direction. Pb islands are presented as stabilized jellium overlayers. The experimental eigenenergies of the quantum well states confined in the Pb overlayer are well reproduced. The total energy oscillates as a continuous function of the overlayer thickness reflecting the electronic shell structure. The energies for completed Pb monolayers show a modulated oscillatory pattern reminiscent of the super-shell structure of clusters and nanowires. The energy minima correlate remarkably well with the measured most probable heights of Pb islands. The proper modeling of the substrate is crucial to set the quantitative agreement. △ Less

Submitted 22 October, 2003; originally announced October 2003.

Comments: 4 pages, 4 figures. Submitted

Broken Symmetry in Density-Functional Theory: Analysis and Cure

Authors: A. Harju, E. Rasanen, H. Saarikoski, M. J. Puska, R. M. Nieminen, K. Niemela

Abstract: We present a detailed analysis of the broken-symmetry mean-field solutions using a four-electron rectangular quantum dot as a model system. Comparisons of the density-functional theory predictions with the exact ones show that the symmetry breaking results from the single-configuration wave function used in the mean-field approach. As a general cure we present a scheme that systematically incorp… ▽ More We present a detailed analysis of the broken-symmetry mean-field solutions using a four-electron rectangular quantum dot as a model system. Comparisons of the density-functional theory predictions with the exact ones show that the symmetry breaking results from the single-configuration wave function used in the mean-field approach. As a general cure we present a scheme that systematically incorporates several configurations into the density-functional theory and restores the symmetry. This cure is easily applicable to any density-functional approach. △ Less

Submitted 17 October, 2003; originally announced October 2003.

Comments: 4 pages, 4 figures, submitted to PRL

Non-Equilibrium Electron Transport in Two-Dimensional Nano-Structures Modeled by Green's Functions and the Finite-Element Method

Authors: Paula Havu, Ville Havu, Martti Puska, Risto Nieminen

Abstract: We use the effective-mass approximation and the density-functional theory with the local-density approximation for modeling two-dimensional nano-structures connected phase-coherently to two infinite leads. Using the non-equilibrium Green's function method the electron density and the current are calculated under a bias voltage. The problem of solving for the Green's functions numerically is form… ▽ More We use the effective-mass approximation and the density-functional theory with the local-density approximation for modeling two-dimensional nano-structures connected phase-coherently to two infinite leads. Using the non-equilibrium Green's function method the electron density and the current are calculated under a bias voltage. The problem of solving for the Green's functions numerically is formulated using the finite-element method (FEM). The Green's functions have non-reflecting open boundary conditions to take care of the infinite size of the system. We show how these boundary conditions are formulated in the FEM. The scheme is tested by calculating transmission probabilities for simple model potentials. The potential of the scheme is demonstrated by determining non-linear current-voltage behaviors of resonant tunneling structures. △ Less

Submitted 6 August, 2003; originally announced August 2003.

Comments: 13 pages,15 figures

Electronic structure of rectangular quantum dots

Authors: E. Räsänen, H. Saarikoski, V. N. Stavrou, A. Harju, M. J. Puska, R. M. Nieminen

Abstract: We study the ground state properties of rectangular quantum dots by using the spin-density-functional theory and quantum Monte Carlo methods. The dot geometry is determined by an infinite hard-wall potential to enable comparison to manufactured, rectangular-shaped quantum dots. We show that the electronic structure is very sensitive to the deformation, and at realistic sizes the non-interacting… ▽ More We study the ground state properties of rectangular quantum dots by using the spin-density-functional theory and quantum Monte Carlo methods. The dot geometry is determined by an infinite hard-wall potential to enable comparison to manufactured, rectangular-shaped quantum dots. We show that the electronic structure is very sensitive to the deformation, and at realistic sizes the non-interacting picture determines the general behavior. However, close to the degenerate points where Hund's rule applies, we find spin-density-wave-like solutions bracketing the partially polarized states. In the quasi-one-dimensional limit we find permanent charge-density waves, and at a sufficiently large deformation or low density, there are strongly localized stable states with a broken spin-symmetry. △ Less

Submitted 20 February, 2003; originally announced February 2003.

Comments: 8 pages, 9 figures, submitted to PRB

Journal ref: Phys. Rev. B 67, 235307 (2003)

Testing of two-dimensional local approximations in the current-spin and spin-density-functional theories

Authors: H. Saarikoski, E. Räsänen, S. Siljamäki, A. Harju, M. J. Puska, R. M. Nieminen

Abstract: We study a model quantum dot system in an external magnetic field by using both the spin-density-functional theory and the current-spin-density-functional theory. The theories are used with local approximations for the spin-density and the vorticity. The reliabilities of different parametrizations for the exchange-correlation functionals are tested by comparing the ensuing energetics with quantu… ▽ More We study a model quantum dot system in an external magnetic field by using both the spin-density-functional theory and the current-spin-density-functional theory. The theories are used with local approximations for the spin-density and the vorticity. The reliabilities of different parametrizations for the exchange-correlation functionals are tested by comparing the ensuing energetics with quantum Monte Carlo results. The limit where the vorticity dependence should be used in the exchange-correlation functionals is discussed. △ Less

Submitted 7 January, 2003; originally announced January 2003.

Comments: 18 pages, 3 figures, 5 eps-figure files

Journal ref: Phys. Rev. B 67, 205327 (2003)

Model study of adsorbed metallic quantum dots: Na on Cu(111)

Authors: T. Torsti, V. Lindberg, M. J. Puska, B. Hellsing

Abstract: We model electronic properties of the second monolayer Na adatom islands (quantum dots) on the Cu(111) surface covered homogeneously by the first Na monolayer. An axially-symmetric three-dimensional jellium model, taking into account the effects due to the first Na monolayer and the Cu substrate, has been developed. The electronic structure is solved within the local-density approximation of the… ▽ More We model electronic properties of the second monolayer Na adatom islands (quantum dots) on the Cu(111) surface covered homogeneously by the first Na monolayer. An axially-symmetric three-dimensional jellium model, taking into account the effects due to the first Na monolayer and the Cu substrate, has been developed. The electronic structure is solved within the local-density approximation of the density-functional theory using a real-space multigrid method. The model enables the study of systems consisting of thousands of Na-atoms. The results for the local density of states are compared with differential conductance ($dI/dV$) spectra and constant current topographs from Scanning Tunneling Microscopy. △ Less

Submitted 4 September, 2002; originally announced September 2002.

Comments: 10 pages, 8 figures. For better quality figures, download http://www.fyslab.hut.fi/~tto/cylart1.pdf

Electronic resonance states in metallic nanowires during the breaking process simulated with the ultimate jellium model

Authors: E. Ogando, T. Torsti, N. Zabala, M. J. Puska

Abstract: We investigate the elongation and breaking process of metallic nanowires using the ultimate jellium model in self-consistent density-functional calculations of the electron structure. In this model the positive background charge deforms to follow the electron density and the energy minimization determines the shape of the system. However, we restrict the shape of the wires by assuming rotational… ▽ More We investigate the elongation and breaking process of metallic nanowires using the ultimate jellium model in self-consistent density-functional calculations of the electron structure. In this model the positive background charge deforms to follow the electron density and the energy minimization determines the shape of the system. However, we restrict the shape of the wires by assuming rotational invariance about the wire axis. First we study the stability of infinite wires and show that the quantum mechanical shell-structure stabilizes the uniform cylindrical geometry at given magic radii. Next, we focus on finite nanowires supported by leads modeled by freezing the shape of a uniform wire outside the constriction volume. We calculate the conductance during the elongation process using the adiabatic approximation and the WKB transmission formula. We also observe the correlated oscillations of the elongation force. In different stages of the elongation process two kinds of electronic structures appear: one with extended states throughout the wire and one with an atom-cluster like unit in the constriction and with well localized states. We discuss the origin of these structures. △ Less

Submitted 3 September, 2002; originally announced September 2002.

Comments: 11 pages, 8 figures

Wigner molecules in polygonal quantum dots: A density functional study

Authors: E. Rasanen, H. Saarikoski, M. J. Puska, R. M. Nieminen

Abstract: We investigate the properties of many-electron systems in two-dimensional polygonal (triangle, square, pentagon, hexagon) potential wells by using the density functional theory. The development of the ground state electronic structure as a function of the dot size is of particular interest. First we show that in the case of two electrons, the Wigner molecule formation agrees with the previous ex… ▽ More We investigate the properties of many-electron systems in two-dimensional polygonal (triangle, square, pentagon, hexagon) potential wells by using the density functional theory. The development of the ground state electronic structure as a function of the dot size is of particular interest. First we show that in the case of two electrons, the Wigner molecule formation agrees with the previous exact diagonalization studies. Then we present in detail how the spin symmetry breaks in polygonal geometries as the spin density functional theory is applied. In several cases with more than two electrons, we find a transition to the crystallized state, yielding coincidence with the number of density maxima and the electron number. We show that this transition density, which agrees reasonably well with previous estimations, is rather insensitive to both the shape of the dot and the electron number. △ Less

Submitted 13 August, 2002; originally announced August 2002.

Comments: 8 pages, 11 figures

Journal ref: Phys. Rev. B 67, 035326 (2003)

MIKA: a multigrid-based program package for electronic structure calculations

Authors: T. Torsti, M. Heiskanen, M. J. Puska, R. M. Nieminen

Abstract: A general real-space multigrid algorithm MIKA (Multigrid Instead of the K-spAce) for the self-consistent solution of the Kohn-Sham equations appearing in the state-of-the-art electronic-structure calculations is described. The most important part of the method is the multigrid solver for the Schrödinger equation. Our choice is the Rayleigh quotient multigrid method (RQMG), which applies directly… ▽ More A general real-space multigrid algorithm MIKA (Multigrid Instead of the K-spAce) for the self-consistent solution of the Kohn-Sham equations appearing in the state-of-the-art electronic-structure calculations is described. The most important part of the method is the multigrid solver for the Schrödinger equation. Our choice is the Rayleigh quotient multigrid method (RQMG), which applies directly to the minimization of the Rayleigh quotient on the finest level. Very coarse correction grids can be used, because there is in principle no need to represent the states on the coarse levels. The RQMG method is generalized for the simultaneous solution of all the states of the system using a penalty functional to keep the states orthogonal. Special care has been taken to optimize the iterations towards the self-consistency and to run the code in parallel computer architectures. The scheme has been implemented in multiple geometries. We show examples from electronic structure calculations employing nonlocal pseudopotentials and/or the jellium model. The RQMG solver is also applied for the calculation of positron states in solids. △ Less

Submitted 2 May, 2002; originally announced May 2002.

Comments: To appear in a special issue of Int J. Quant. Chem. devoted to the conference proceedings of 9th International Conference on the Applications of the Density Functional Theory in Chemistry and Physics

Journal ref: Int. J. Quantum Chem. 91, 171 (2003)

Conductance oscillations in metallic nanocontacts

Authors: P. Havu, T. Torsti, M. J. Puska, R. M. Nieminen

Abstract: We examine the conductance properties of a chain of Na atoms between two metallic leads in the limit of low bias. Resonant states corresponding to the conductance channel and the local charge neutrality condition cause conductance oscillations as a function of the number of atoms in the chain. Moreover, the geometrical shape of the contact leads influences the conductivity by giving rise to addi… ▽ More We examine the conductance properties of a chain of Na atoms between two metallic leads in the limit of low bias. Resonant states corresponding to the conductance channel and the local charge neutrality condition cause conductance oscillations as a function of the number of atoms in the chain. Moreover, the geometrical shape of the contact leads influences the conductivity by giving rise to additional oscillations as a function of the lead opening angle. △ Less

Submitted 2 May, 2002; originally announced May 2002.

Comments: 5 pages, 7 figures

Spontaneous magnetization of aluminum nanowires deposited on the NaCl(100) surface

Authors: A. Ayuela, H. Räbiger, M. J. Puska, R. M. Nieminen

Abstract: We investigate electronic structures of Al quantum wires, both unsupported and supported on the (100) NaCl surface, using the density-functional theory. We confirm that unsupported nanowires, constrained to be linear, show magnetization when elongated beyond the equilibrium length. Allowing ions to relax, the wires deform to zig-zag structures with lower magnetization but no dimerization occurs.… ▽ More We investigate electronic structures of Al quantum wires, both unsupported and supported on the (100) NaCl surface, using the density-functional theory. We confirm that unsupported nanowires, constrained to be linear, show magnetization when elongated beyond the equilibrium length. Allowing ions to relax, the wires deform to zig-zag structures with lower magnetization but no dimerization occurs. When an Al wire is deposited on the NaCl surface, a zig-zag geometry emerges again. The magnetization changes moderately from that for the corresponding unsupported wire. We analyse the findings using electron band structures and simple model wires. △ Less

Submitted 19 March, 2001; originally announced March 2001.

Comments: submitted to PHys. Rev. B

A novel multigrid method for electronic structure calculations

Authors: M. Heiskanen, T. Torsti, M. J. Puska, R. M. Nieminen

Abstract: A general real-space multigrid algorithm for the self-consistent solution of the Kohn-Sham equations appearing in the state-of-the-art electronic-structure calculations is described. The most important part of the method is the multigrid solver for the Schroedinger equation. Our choice is the Rayleigh quotient multigrid method (RQMG), which applies directly to the minimization of the Rayleigh qu… ▽ More A general real-space multigrid algorithm for the self-consistent solution of the Kohn-Sham equations appearing in the state-of-the-art electronic-structure calculations is described. The most important part of the method is the multigrid solver for the Schroedinger equation. Our choice is the Rayleigh quotient multigrid method (RQMG), which applies directly to the minimization of the Rayleigh quotient on the finest level. Very coarse correction grids can be used, because there is no need to be able to represent the states on the coarse levels. The RQMG method is generalized for the simultaneous solution of all the states of the system using a penalty functional to keep the states orthogonal. The performance of the scheme is demonstrated by applying it in a few molecular and solid-state systems described by non-local norm-conserving pseudopotentials. △ Less

Submitted 16 January, 2001; originally announced January 2001.

Comments: 9 pages, 3 figures

Journal ref: Phys. Rev. B 63, 245106 (2001)

First-Principles Calculations of Positron Annihilation in Solids

Authors: B. Barbiellini, M. Hakala, R. M. Nieminen, M. J. Puska

Abstract: We present first-principles approaches based on density functional theory for calculating positron states and annihilation characteristics in condensed matter. The treatment of the electron-positron correlation effects (the enhancement of the electron density at the positron with respect to mean-field density) is shown to play a crucial role when calculating the annihilation rates. A generalized… ▽ More We present first-principles approaches based on density functional theory for calculating positron states and annihilation characteristics in condensed matter. The treatment of the electron-positron correlation effects (the enhancement of the electron density at the positron with respect to mean-field density) is shown to play a crucial role when calculating the annihilation rates. A generalized gradient approximation (GGA) takes the strong inhomogeneities of the electron density in the ion core region into account and reproduces well the experimental total annihilation rates (inverses of the positron lifetimes) by suppressing the rates given by a local density approximation (LDA). The GGA combined with an electron-state-dependent enhancement scheme gives a good description for the momentum distributions of the annihilating positron-electron pairs reproducing accurately the trends observed in the angular correlation (ACAR) or Doppler broadening measurements of the annihilation radiation. The combination of the present positron lifetime and momentum density calculations with the corresponding measurements yields a unique tool for defect identification. Especially, the investigation of various vacancy-type defects in semiconductors able to trap positrons will be an important field for these methods. We will show that the identification of vacancy-impurity complexes in highly n-Type Si and the study of the SiO$_2$/Si interface are particularly interesting applications. △ Less

Submitted 15 December, 1999; v1 submitted 14 December, 1999; originally announced December 1999.

Comments: 6 pages, 3 figures, MRS style (enclosed), to appear in MRS 1999 Fall Meeting Proc. 579; added acknowledgments

Journal ref: MRS Proc. Fall 1999 Vol. 579 (2000) p. 249