Browsing by Subject "ATOMS"

Sort by: Order: Results:

Now showing items 1-10 of 10
  • Lehtola, Susi (2020)
    Knowledge of the repulsive behavior of potential energy curves V (R) at R -> 0 is necessary for understanding and modeling irradiation processes of practical interest. V (R) is in principle straightforward to obtain from electronic structure calculations; however, commonly used numerical approaches for electronic structure calculations break down in the strongly repulsive region due to the closeness of the nuclei. In this work, we show by comparison to fully numerical reference values that a recently developed procedure [S. Lehtola, J. Chem. Phys. 151, 241102 (2019)] can be employed to enable accurate linear combination of atomic orbitals calculations of V (R) even at small R by a study of the seven nuclear reactions He-2 (sic) Be, HeNe (sic) Mg, Ne-2 (sic) Ca, HeAr (sic) Ca, MgAr (sic) Zn, Ar-2 (sic) Kr, and NeCa (sic) Zn.
  • Buchs, Gilles; Krasheninnikov, Arkady V.; Ruffieux, Pascal; Groening, Pierangelo; Foster, Adam S.; Nieminen, Risto M.; Groening, Oliver (2007)
  • Heikinheimo, J.; Mizohata, K.; Räisänen, J.; Ahlgren, T.; Jalkanen, P.; Lahtinen, A.; Catarino, N.; Alves, E.; Tuomisto, Filip (2019)
    Reliable and accurate knowledge of the physical properties of elementary point defects is crucial for predictive modeling of the evolution of radiation damage in materials employed in harsh conditions. We have applied positron annihilation spectroscopy to directly detect mono-vacancy defects created in tungsten through particle irradiation at cryogenic temperatures, as well as their recovery kinetics. We find that efficient self-healing of the primary damage takes place through Frenkel pair recombination already at 35 K, in line with an upper bound of 0.1 eV for the migration barrier of self-interstitials. Further self-interstitial migration is observed above 50 K with activation energies in the range of 0.12-0.42 eV through the release of the self-interstitial atoms from impurities and structural defects and following recombination with mono-vacancies. Mono-vacancy migration is activated at around 550 K with a migration barrier of E-m(V) = 1.85 +/- 0.05 eV. (C) 2019 Author(s).
  • Lehtola, Susi (2019)
    We present the implementation of a variational finite element solver in the HelFEM program for benchmark calculations on diatomic systems. A basis set of the form chi nlm mu nu phi=Bn mu Ylm nu phi is used, where (mu, nu, phi) are transformed prolate spheroidal coordinates, B-n(mu) are finite element shape functions, and Ylm are spherical harmonics. The basis set allows for an arbitrary level of accuracy in calculations on diatomic molecules, which can be performed at present with either nonrelativistic Hartree-Fock (HF) or density functional (DF) theory. Hundreds of DFs at the local spin density approximation (LDA), generalized gradient approximation (GGA), and the meta-GGA level can be used through an interface with the Libxc library; meta-GGA and hybrid DFs are not available in other fully numerical diatomic program packages. Finite electric fields are also supported in HelFEM, enabling access to electric properties. We introduce a powerful tool for adaptively choosing the basis set by using the core Hamiltonian as a proxy for its completeness. HelFEM and the novel basis set procedure are demonstrated by reproducing the restricted open-shell HF limit energies of 68 diatomic molecules from the first to the fourth period with excellent agreement with literature values, despite requiring orders of magnitude fewer parameters for the wave function. Then, the electric properties of the BH and N-2 molecules under finite field are studied, again yielding excellent agreement with previous HF limit values for energies, dipole moments, and dipole polarizabilities, again with much more compact wave functions than what were needed for the literature references. Finally, HF, LDA, GGA, and meta-GGA calculations of the atomization energy of N-2 are performed, demonstrating the superb accuracy of the present approach.
  • Herbig, Charlotte; Åhlgren, E. Harriet; Jolie, Wouter; Busse, Carsten; Kotakoski, Jani; Krasheninnikov, Arkady V.; Michely, Thomas (2014)
  • Rauhalahti, Markus; Sundholm, Dage; Johansson, Mikael P. (2021)
    The magnetically induced current density of an intriguing naphthalene-fused heteroporphyrin has been studied, using the quantum-chemical, gauge-including magnetically induced currents (GIMIC) method. The ring-current strengths and current-density pathways for the heteroporphyrin, its Pd complex, and the analogous quinoline-fused heteroporphyrin provide detailed information about their aromatic properties. The three porphyrinoids have similar current-density pathways and are almost as aromatic as free-base porphyrin. Notably, we show that the global ring current makes a branch at three specific points. Thus, the global current is composed of a total of eight pathways that include 22 pi-electrons, with no contributions from 18-electron pathways.
  • Granberg, F.; Byggmastar, J.; Nordlund, K. (2021)
    Tungsten has been chosen as the plasma-facing wall material in fusion reactors, due to its high density and melting point. The wall material will not only be sputtered at the surface, but also damaged deep inside the material by energetic particles. We investigate the high-dose damage production and accumulation by computational means using molecular dynamics. We observe that the choice of interatomic potential drastically affects the evolution. The structure and stability of the obtained defect configurations are validated using a quantum-accurate Gaussian approximation potential. (c) 2021 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license ( )
  • Ryazantsev, Sergey V.; Tyurin, Daniil A.; Feldman, Vladimir I.; Khriachtchev, Leonid (2017)
    We report on the preparation and vibrational characterization of the C2H3 center dot center dot center dot CO2 complex, the first example of a stable intermolecular complex involving vinyl radicals. This complex was prepared in Ar and Kr matrices using UV photolysis of propiolic acid (HC3OOH) and subsequent thermal mobilization of H atoms. This preparation procedure provides vinyl radicals formed exclusively as a complex with CO2, without the presence of either CO2 or C2H3 monomers. The absorption bands corresponding to the nu(5)(C2H3), nu(7)(C2H3), nu(8)(C2H3), nu(2)(CO2), and nu(3)(CO2) modes of the C2H3 center dot center dot center dot CO2 complex were detected experimentally. The calculations at the UCCSD(T)/L2a level of theory predict two structures of the C2H3 center dot center dot center dot CO2 complex with C-s and C-1 symmetries and interaction energies of -1.92 and -5.19 kJ mol(-1). The harmonic vibrational frequencies of these structures were calculated at the same level of theory. The structural assignment of the experimental species is not straightforward because of rather small complexation-induced shifts and matrix-site splitting of the bands (for both complex and monomers). We conclude that the C-1 structure is the most probable candidate for the experimental C2H3 center dot center dot center dot O-2 complex based on the significant splitting of the bending vibration of CO2 and on the energetic and structural considerations. Published by AIP Publishing.
  • Shahi, Chandra; Bhattarai, Puskar; Wagle, Kamal; Santra, Biswajit; Schwalbe, Sebastian; Hahn, Torsten; Kortus, Jens; Jackson, Koblar A.; Peralta, Juan E.; Trepte, Kai; Lehtola, Susi; Nepal, Niraj K.; Myneni, Hemanadhan; Neupane, Bimal; Adhikari, Santosh; Ruzsinszky, Adrienn; Yamamoto, Yoh; Baruah, Tunna; Zope, Rajendra R.; Perdew, John P. (2019)
    Semilocal approximations to the density functional for the exchange-correlation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semilocal approximation makes that approximation exact for all one-electron ground- or excited-state densities and accurate for stretched bonds. When the minimization of the PZ total energy is made over real localized orbitals, the orbital densities can be noded, leading to energy errors in many-electron systems. Minimization over complex localized orbitals yields nodeless orbital densities, which reduce but typically do not eliminate the SIC errors of atomization energies. Other errors of PZ SIC remain, attributable to the loss of the exact constraints and appropriate norms that the semilocal approximations satisfy, suggesting the need for a generalized SIC. These conclusions are supported by calculations for one-electron densities and for many-electron molecules. While PZ SIC raises and improves the energy barriers of standard generalized gradient approximations (GGAs) and meta-GGAs, it reduces and often worsens the atomization energies of molecules. Thus, PZ SIC raises the energy more as the nodality of the valence localized orbitals increases from atoms to molecules to transition states. PZ SIC is applied here, in particular, to the strongly constrained and appropriately normed (SCAN) meta-GGA, for which the correlation part is already self-interaction-free. This property makes SCAN a natural first candidate for a generalized SIC. Published under license by AIP Publishing.