Relativistic calculations of atomic structure
Fricke, Burkhard
1984-01-01
A review of relativistic atomic structure calculations is given with a emphasis on the Multiconfigurational-Dirac-Fock method. Its problems and deficiencies are discussed together with the contributions which go beyond the Dirac-Fock procedure.
Atomic Structure Calculations for Neutral Oxygen
Alonizan, Norah; Qindeel, Rabia; Ben Nessib, Nabil
2016-01-01
Energy levels and oscillator strengths for neutral oxygen have been calculated using the Cowan (CW), SUPERSTRUCTURE (SS), and AUTOSTRUCTURE (AS) atomic structure codes. The results obtained with these atomic codes have been compared with MCHF calculations and experimental values from the National Institute of Standards and Technology (NIST) database.
METHODS OF CALCULATING THE ELECTRONIC AND ATOMIC STRUCTURES OF INTERFACES
Sutton, A
1985-01-01
Methods of calculating the electronic and atomic structures of interfaces are described. An introduction to pseudopotentials and LCAO methods is given. Methods of calculating the electronic structure of an interface with a given atomic structure are considered. The feasibility of total energy calculations, in which the atomic and electronic structures are calculated simultaneously, is discussed.
Atomic structure calculations of Mo XV-XL
Energy levels and oscillator strengths were calculated for Mo XV - Mo XL. The computer program for atomic structure calculation, developed by Dr. Robert D. Cowan, Los Alamos National Laboratory, was used in the present work. The scaled energy parameters were empirically determined from the observed spectral data. We present wavelengths and transition probabilities of Mo XV-XL. Energy levels and spectral patterns are presented in figures that are useful for the identification of spectral lines. (author)
Atomic structure calculations on the CRAY X-MP
Atomic structure calculations require both radial and angular integrations, where the latter are often based on Racah algebra. With relatively minor modifications, good performance is obtained on vector machines for radial integrations. Angular integrations, however, present the bottleneck. In this paper some recent improvements in the algorithms for angular integrations are described, as well as some multitasking experiments on the CRAY X-MP and CRAY 2. These show that the workload can easily be distributed evenly among available processors with dynamic scheduling
Multi-million atom electronic structure calculations for quantum dots
Usman, Muhammad
Quantum dots grown by self-assembly process are typically constructed by 50,000 to 5,000,000 structural atoms which confine a small, countable number of extra electrons or holes in a space that is comparable in size to the electron wavelength. Under such conditions quantum dots can be interpreted as artificial atoms with the potential to be custom tailored to new functionality. In the past decade or so, these nanostructures have attracted significant experimental and theoretical attention in the field of nanoscience. The new and tunable optical and electrical properties of these artificial atoms have been proposed in a variety of different fields, for example in communication and computing systems, medical and quantum computing applications. Predictive and quantitative modeling and simulation of these structures can help to narrow down the vast design space to a range that is experimentally affordable and move this part of nanoscience to nano-Technology. Modeling of such quantum dots pose a formidable challenge to theoretical physicists because: (1) Strain originating from the lattice mismatch of the materials penetrates deep inside the buffer surrounding the quantum dots and require large scale (multi-million atom) simulations to correctly capture its effect on the electronic structure, (2) The interface roughness, the alloy randomness, and the atomistic granularity require the calculation of electronic structure at the atomistic scale. Most of the current or past theoretical calculations are based on continuum approach such as effective mass approximation or k.p modeling capturing either no or one of the above mentioned effects, thus missing some of the essential physics. The Objectives of this thesis are: (1) to model and simulate the experimental quantum dot topologies at the atomistic scale; (2) to theoretically explore the essential physics i.e. long range strain, linear and quadratic piezoelectricity, interband optical transition strengths, quantum confined
Atomic structure calculations using the relativistic random phase approximation
A brief review is given for the relativistic random phase approximation (RRPA) applied to atomic transition problems. Selected examples of RRPA calculations on discrete excitations and photoionization are given to illustrate the need of relativistic many-body theories in dealing with atomic processes where both relativity and correlation are important
A new program for calculating matrix elements in atomic structure
The solution of many problems concerning the electronic structure of atoms requires the evaluation of the matrix elements of the Hamiltonian operator, including the electrostatic interaction. These matrix elements may be expressed as weighted sums of radial integrals. The program we describe in this paper evaluates the coefficients of the Slater integrals and, if these are given, computes all the matrix elements for a given set of configurations. This program has nearly the same purposes as Hibbert's program and is also based on the Racach techniques. The main difference between this algorithm and the cited one is the method used to calculate the recoupling coefficients. While Hibbert's programs use Burke's algorithm to calculate these coefficients, in our program they are computed using the graphical techniques developed by Jucys et al. According to this method, that we describe in another paper the formulae needed to calculate the recoupling coefficients are previously derived and simplified (as a first step of the program). The use of this method may considerably reduce the running time, specially in the case of large configuration interaction matrices. (orig.)
Ab initio calculations and modelling of atomic cluster structure
Solov'yov, Ilia; Lyalin, Andrey G.; Solov'yov, Andrey V.;
2004-01-01
framework for modelling the fusion process of noble gas clusters is presented. We report the striking correspondence of the peaks in the experimentally measured abundance mass spectra with the peaks in the size-dependence of the second derivative of the binding energy per atom calculated for the chain...
Atomic structure calculations for F-like tungsten
Sunny, Aggarwal
2014-09-01
Energy levels, wavefunction compositions and lifetimes have been computed for all levels of 1s22s22p5, 1s22s2p6, 1s22s22p43s, 1s22s22p43p, and 1s22s22p43d configurations in highly charged F-like tungsten ion. The multiconfigurational Dirac—Fock method (MCDF) is adopted to generate the wavefunctions. We have also presented the transition wavelengths, oscillator strengths, transition probabilities, and line strengths for the electric dipole (E1) and magnetic quadrupole (M2) transition from the 1s22s22p5 ground configuration. We have performed parallel calculations with the flexible atomic code (FAC) for comparing the atomic data. The reliability of present data is assessed by comparison with other theoretical and experimental data available in the literature. Good agreement is found between our results and those obtained using different approaches confirm the quality of our results. Additionally, we have predicted some new atomic data for F-like W that were not available so far and may be important for plasma diagnostic analysis in fusion plasma.
Hyun-Kyung Chung; Per Jönsson; Alexander Kramida
2013-01-01
Atomic structure and transition probabilities are fundamental physical data required in many fields of science and technology. Atomic physics codes are freely available to other community users to generate atomic data for their interest, but the quality of these data is rarely verified. This special issue addresses estimation of uncertainties in atomic structure and transition probability calculations, and discusses methods and strategies to assess and ensure the quality of theoretical atomic...
Advances in Atomic Structure Calculations%原子结构计算的进展
Charlotte Froese Fischer
2007-01-01
Correlation and relativistic effects are both needed for accurate atomic structure calculations of energy levels and their atomic properties. For transition probabilities of radiative transitions between low-lying levels of an atom or ion, accurate wave functions for the outer region of are required. For lighter atoms, relativistic effects can be included through the Breit-Pauli approximation. This paper outlines the advances in the treatment of correlation and describes the current state of Breit-Pauli calculations for complex systems.
The structural, electronic and magnetic properties of free standing Au-Pd bimetallic atomic chain is studied using ab-initio method. It is found that electronic and magnetic properties of chains depend on position of atoms and number of atoms. Spin polarization factor for different atomic configuration of atomic chain is calculated predicting a half metallic behavior. It suggests a total spin polarised transport in these chains
Amani Tahat; Mahmoud Abu-Allaban; Safeia Hamasha
2011-01-01
In this study, a new atomic physics program (HTAC) is introduced and tested. It is a utility program designed to automate the computation of various atomic structure and spectral data. It is the first comprehensive code that enables performing atomic calculations based on three advanced theories: the fully relativistic configuration interactions approach, the multi-reference many body perturbation theory and the R-Matrix method. It has been designed to generate tabulated atomic data files tha...
A Hartree-Fock program for atomic structure calculations
The Hartree-Fock equations for a general open shell atom are described. The matrix equations that result when the single particle orbitals are written in terms of a linear combination of analytic basis functions are derived. Attention is paid to the complexities that occur when open shells are present. The specifics of a working FORTRAN program which is available for public use are described. The program has the flexibility to handle either Slater-type orbitals or Gaussian-type orbitals. It can be obtained over the internet at http://lacebark.ntu.edu.au/j_mitroy/research/atomic.htm Copyright (1999) CSIRO Australia
Dzuba, V. A.; Flambaum, V. V.
A brief review of the search for variation of the fine structure constant in quasar absorption spectra is presented. Special consideration is given to the role of atomic calculations in the analysis of the observed data. A range of methods which allow to perform calculations for atoms or ions with different electron structure and which cover practically all periodic table of elements is discussed. Critical compilation of the results of the calculations as well as a review of the most recent results of the analysis are presented.
Dzuba, V A
2008-01-01
A brief review of the search for variation of the fine structure constant in quasar absorption spectra is presented. Special consideration is given to the role of atomic calculations in the analysis of the observed data. A range of methods which allow to perform calculations for atoms or ions with different electron structure and which cover practically all periodic table of elements is discussed. Critical compilation of the results of the calculations as well as a review of the most recent results of the analysis are presented.
Theoretical calculations on the atomic and electronic structure of β-SiC(110) surface
无
2002-01-01
We present a theoretical calculation of the atomic and electronic structure of β-SiC and its non-polar (110) surface using the full potential linear augmented plane wave (FPLAPW) approach. The calculated lattice constant and bulk modulus of β-SiC crystal are in excellent agreement with experimental data. The atomic and electronic structure of β-SiC(110) surface has been calculated by employing the slab and supercell model. It is found that the surface is characterized by a top-layer bond-length-contracting rotation relaxation in which the Si-surface atom moves closer towards the substrate while the C-surface atom moves outward. This relaxation is analogous to that of Ⅲ-Ⅴ semi-conductor surface. The driving mechanism for this atomic rearrangement is that the Si atom tends to a planar sp2-like bonding situation with its three N neighbors and the N atom tends to a p3-like bonding with its three Si neighbors. Furthermore, surface relaxation induces the change from metallic to semiconducting characterization.
Ansari, Reza; Ajori, Shahram; Malakpour, Sina
2016-04-01
The considerable demand for novel materials with specific properties has motivated the researchers to synthesize supramolecular nanostructures through different methods. Porous graphene is the first two-dimensional hydrocarbon synthesized quite recently. This investigation is aimed at studying the mechanical properties of atom-decorated (functionalized) porous graphene by employing density functional theory (DFT) calculation within both local density approximations (LDA) and generalized gradient approximations (GGA). The atoms are selected from period 3 of periodic table as well as Li and O atom from period 2. The results reveal that metallic atoms and noble gases are adsorbed physically on porous graphene and nonmetallic ones form chemical bonds with carbon atom in porous graphene structure. Also, it is shown that, in general, atom decoration reduces the values of mechanical properties such as Young's, bulk and shear moduli as well as Poisson's ratio, and this reduction is more considerable in the case of nonmetallic atoms (chemical adsorption), especially oxygen atoms, as compared to metallic atoms and noble gases (physical adsorption).
Per Jönsson; Hyun-Kyung Chung
2013-01-01
There exist several codes in the atomic physics community to generate atomic structure and transition probabilities freely and readily distributed to researchers outside atomic physics community, in plasma, astrophysical or nuclear physics communities. Users take these atomic physics codes to generate the necessary atomic data or modify the codes for their own applications. However, there has been very little effort to validate and verify the data sets generated by non-expert users. [...
Atomic structure calculations and identification of EUV and SXR spectral lines in Sr XXX
Goyal, Arun; Khatri, Indu; Aggarwal, Sunny; Singh, A. K.; Mohan, Man
2015-08-01
We report an extensive theoretical study of atomic data for Sr XXX in a wide range with L-shell electron excitations to the M-shell. We have calculated energy levels, wave-function compositions and lifetimes for lowest 113 fine structure levels and wavelengths of an extreme Ultraviolet (EUV) and soft X-ray (SXR) transitions. We have employed multi-configuration Dirac Fock method (MCDF) approach within the framework of Dirac-Coulomb Hamiltonian including quantum electrodynamics (QED) and Breit corrections. We have also presented the radiative data for electric and magnetic dipole (E1, M1) and quadrupole (E2, M2) transitions from the ground state. We have made comparisons with available energy levels compiled by NIST and achieve good agreement. But due to inadequate data in the literature, analogous relativistic distorted wave calculations have also been performed using flexible atomic code (FAC) to assess the reliability and accuracy of our results. Additionally, we have provided new atomic data for Sr XXX which is not published elsewhere in the literature and we believe that our results may be beneficial in fusion plasma research and astrophysical investigations and applications.
Atomic structure calculations and identification of EUV and SXR spectral lines in Sr XXX
We report an extensive theoretical study of atomic data for Sr XXX in a wide range with L-shell electron excitations to the M-shell. We have calculated energy levels, wave-function compositions and lifetimes for lowest 113 fine structure levels and wavelengths of an extreme Ultraviolet (EUV) and soft X-ray (SXR) transitions. We have employed multi-configuration Dirac Fock method (MCDF) approach within the framework of Dirac–Coulomb Hamiltonian including quantum electrodynamics (QED) and Breit corrections. We have also presented the radiative data for electric and magnetic dipole (E1, M1) and quadrupole (E2, M2) transitions from the ground state. We have made comparisons with available energy levels compiled by NIST and achieve good agreement. But due to inadequate data in the literature, analogous relativistic distorted wave calculations have also been performed using flexible atomic code (FAC) to assess the reliability and accuracy of our results. Additionally, we have provided new atomic data for Sr XXX which is not published elsewhere in the literature and we believe that our results may be beneficial in fusion plasma research and astrophysical investigations and applications. - Highlights: • 113 Lowest levels for Sr XXX are calculated. • Extreme Ultraviolet (EUV) and soft-X ray (SXR) spectral lines are identified. • Wavelengths of EUV and SXR spectral lines are reported. • E1, E2, M1 and M2 transition rates, oscillator strengths and lines strengths for lowest 113 levels are presented. • Lifetimes for lowest 113 fine structure levels are provided
GRASP92: a package for large-scale relativistic atomic structure calculations
Parpia, F. A.; Froese Fischer, C.; Grant, I. P.
2006-12-01
of CSFs sharing the same quantum numbers is determined using the configuration-interaction (CI) procedure that results upon varying the expansion coefficients to determine the extremum of a variational functional. Radial functions may be determined by numerically solving the multiconfiguration Dirac-Fock (MCDF) equations that result upon varying the orbital radial functions or some subset thereof so as to obtain an extremum of the variational functional. Radial wavefunctions may also be determined using a screened hydrogenic or Thomas-Fermi model, although these schemes generally provide initial estimates for MCDF self-consistent-field (SCF) calculations. Transition properties for pairs of ASFs are computed from matrix elements of multipole operators of the electromagnetic field. All matrix elements of CSFs are evaluated using the Racah algebra. Reasons for the new version: During recent studies using the general relativistic atomic structure package (GRASP92), several errors were found, some of which might have been present already in the earlier GRASP92 version (program ABJN_v1_0, Comput. Phys. Comm. 55 (1989) 425). These errors were reported and discussed by Froese Fischer, Gaigalas, and Ralchenko in a separate publication [C. Froese Fischer, G. Gaigalas, Y. Ralchenko, Comput. Phys. Comm. 175 (2006) 738-744. [7
Elmar Träbert
2014-03-01
Full Text Available The interpretation of atomic observations by theory and the testing of computational predictions by experiment are interactive processes. It is necessary to gain experience with “the other side” before claims of achievement can be validated and judged. The discussion covers some general problems in the field as well as many specific examples, mostly organized by isoelectronic sequence, of what level of accuracy recently has been reached or which atomic structure or level lifetime problem needs more attention.
Elmar Träbert
2014-01-01
The interpretation of atomic observations by theory and the testing of computational predictions by experiment are interactive processes. It is necessary to gain experience with “the other side” before claims of achievement can be validated and judged. The discussion covers some general problems in the field as well as many specific examples, mostly organized by isoelectronic sequence, of what level of accuracy recently has been reached or which atomic structure or level lifetime problem need...
This report presents details of a new method to enable the computation of collision strengths for complex ions which is adapted from long established optimisation techniques previously applied to the calculation of atomic structures and oscillator strengths. The procedure involves the adjustment of Slater parameters so that they determine improved energy levels and eigenvectors. They provide a basis for collision strength calculations in ions where ab initio computations break down or result in reducible errors. This application is demonstrated through modifications of the DISTORTED WAVE collision code and SUPERSTRUCTURE atomic-structure code which interface via a transformation code JAJOM which processes their output. (author)
Corsetti, Fabiano
2014-01-01
The implementation of the orbital minimization method (OMM) for solving the self-consistent Kohn-Sham (KS) problem for electronic structure calculations in a basis of non-orthogonal numerical atomic orbitals of finite-range is reported. We explore the possibilities for using the OMM as an exact cubic-scaling solver for the KS problem, and compare its performance with that of explicit diagonalization in realistic systems. We analyze the efficiency of the method depending on the choice of line search algorithm and on two free parameters, the scale of the kinetic energy preconditioning and the eigenspectrum shift. The results of several timing tests are then discussed, showing that the OMM can achieve a noticeable speedup with respect to diagonalization even for minimal basis sets for which the number of occupied eigenstates represents a significant fraction of the total basis size (>15%). We investigate the hard and soft parallel scaling of the method on multiple cores, finding a performance equal to or better ...
Uranium dioxide UO2 is the standard nuclear fuel used in pressurized water reactors. During in-reactor operation, the fission of uranium atoms yields a wide variety of fission products (FP) which create numerous point defects while slowing down in the material. Point defects and FP govern in turn the evolution of the fuel physical properties under irradiation. In this study, we use electronic structure calculations in order to better understand the fuel behavior under irradiation. In particular, we investigate point defect behavior, as well as the stability of three volatile FP: iodine, krypton and xenon. In order to take into account the strong correlations of uranium 5f electrons in UO2, we use the DFT+U approximation, based on the density functional theory. This approximation, however, creates numerous metastable states which trap the system and induce discrepancies in the results reported in the literature. To solve this issue and to ensure the ground state is systematically approached as much as possible, we use a method based on electronic occupancy control of the correlated orbitals. We show that the DFT+U approximation, when used with electronic occupancy control, can describe accurately point defect and fission product behavior in UO2 and provide quantitative information regarding point defect transport properties in the oxide fuel. (author)
Calculations of effective atomic number
Kaliman, Z. [Department of Physics, Faculty of Arts and Sciences, Omladinska 14, Rijeka (Croatia); Orlic, N. [Department of Physics, Faculty of Arts and Sciences, Omladinska 14, Rijeka (Croatia)], E-mail: norlic@ffri.hr; Jelovica, I. [Department of Physics, Faculty of Arts and Sciences, Omladinska 14, Rijeka (Croatia)
2007-09-21
We present and discuss effective atomic number (Z{sub eff}) obtained by different methods of calculations. There is no unique relation between the computed values. This observation led us to the conclusion that any Z{sub eff} is valid only for given process. We illustrate calculations for different subshells of atom Z=72 and for M3 subshell of several other atoms.
Comparative semi-empirical and ab initio atomic structure calculations in Yb-like tungsten W4+
In this paper, we report on extensive calculations of radiative data in Yb-like tungsten ion using several independent atomic structure methods, i.e. the relativistic Hartree–Fock approach, the flexible atomic code and the multiconfiguration Dirac–Fock method. This multi-platform approach allowed us to check the consistency of our results. Advantages and shortcomings of semi-empirical and ab initio methods for atomic structure calculations in such a complex heavy ion are also discussed in detail. A new set of transition probabilities and oscillator strengths is reported for electric dipole lines together with magnetic dipole and electric quadrupole lines in this ion of interest for fusion plasma diagnostics. (paper)
A robust and general Schrödinger and Dirac solver for atomic structure calculations
Čertík, O.; Pask, J.E.; Vackář, Jiří
2013-01-01
Roč. 184, č. 7 (2013), s. 1777-1791. ISSN 0010-4655 R&D Projects: GA MŠk(CZ) LC06040; GA ČR GA101/09/1630 Institutional support: RVO:68378271 Keywords : atom * electronic structure * Dirac equation * density -functional theory Subject RIV: BE - Theoretical Physics Impact factor: 2.407, year: 2013 http://www.sciencedirect.com/science/article/pii/S0010465513000714
Atomic and electronic structure of hydrogen on ZnO (1bar 100) surface: ab initio hybrid calculations
Usseinov, A. B.; Kotomin, E. A.; Zhukovskii, Yu F.; Purans, J.; Sorokin, A. V.; Akilbekov, A. T.
2013-12-01
Hydrogen atoms unavoidably incorporated into ZnO during growth of bulk samples and thin films considerably affect their electrical conductivity. The results of first principles hybrid LCAO calculations are discussed for hydrogen atoms in the bulk and on the non-polar ZnO (1bar 100) surface. The incorporation energy, the atomic relaxation, the electronic density redistribution and the electronic structure modifications are compared for the surface adsorption and bulk interstitial H positions. It is shown that hydrogen has a strong binding with the surface O ions (2.7 eV) whereas its incorporation into bulk is energetically unfavorable. Surface hydrogen atoms are very shallow donors, thus, contributing to the electronic conductivity.
Atomic and Electronic Structures of C_60+BN Nanopeapods from ab initio Pseudopotential Calculations
Trave, Andrea; Ribeiro, Filipe; Louie, Steven G.; Cohen, Marvin L.
2004-03-01
Nanopeapods are structures of nanometric size consisting of an external carbon nanotube encapsulating a chain or complex array of fullerenes. Recent calculations and experiments have proven that nanopeapods can be obtained assembling fullerenes within boron nitride nanotubes, creating novel materials of possible interest for electronic transport applications. To improve the understanding of the properties of these composite systems, as compared to empty nanotubes and carbon nanopeapods, ab-initio total energy calculations have been performed within the pseudopotential Density Functional Theory in local density approximation. Results of these calculations on the energetics and geometrical deformations involved in the encapsulation will be presented, followed by a discussion of the consequences on the electronic structures of these systems, with particular focus on aspects relevant to electronic transport phenomena. This work is supported by NFS (Grant DMR00-87088) and DOE (Contract DE-AC03-76SF00098), using computational resources at NERSC and NPACI.
Landau, Arie; Kaprálová-Žďánská, Petra Ruth; Moiseyev, Nimrod
2015-01-01
Complex eigenvalues, resonances, play an important role in large variety of fields in physics and chemistry. For example, in cold molecular collision experiments and electron scattering experiments, autoionizing and pre-dissociative metastable resonances are generated. However, the computation of complex resonance eigenvalues is difficult, since it requires severe modifications of standard electronic structure codes and methods. Here we show how resonance eigenvalues, positions and widths, can be calculated using the standard, widely used, electronic-structure packages. Our method enables the calculations of the complex resonance eigenvalues by using analytical continuation procedures (such as Pad\\'{e}). The key point in our approach is the existence of narrow analytical passages from the real axis to the complex energy plane. In fact, the existence of these analytical passages relies on using finite basis sets. These passages become narrower as the basis set becomes more complete, whereas in the exact limit,...
An accurate and efficient method is described for the evaluation of electrostatic contributions in LCAO electronic structure calculations. The charge density rho(r) is decomposed into rho/sup(1)(r), a component whose rapid variation near any nucleus reproduces that of rho(r) to a very good approximation, and a remainder density deltarho(r)equivalentrho(r)-rho/sup(1)(r), which is thereby guaranteed to be slowly varying in space. The power of the decomposition resides in the fact that rho/sup(1)(r) can be expressed exactly as a sum of one-center densities, without the use of any fit procedure. Because rho/sup(1)(r) is a sum of one-center multipolar densities, the Hartree potential is a function with a simple one-dimensional integral representation, and its matrix elements can be obtained by performing one-dimensional integrals over it. Since deltarho(r) is spatially slowly varying, the Hartree potential to which it corresponds and the matrix elements of this potential can accurately be evaluated on a relatively coarse coordinate space mesh, using fast Fourier transforms. The method is illustrated via molecular structure calculations for N2 and NH3. The calculations are accurate to a few percent when the required integrals over deltarho(r) and deltaV(r) are performed on a mesh of spacing 0.4 a.u. The N--N bond length and stretch frequency are found to equal 2.10 a.u. and 2.3 x 103 cm-1, respectively. The equilibrium N--H bond length and H--N--H angle are calculated to be 1.93 a.u. and 1050, respectively, while the NH3 inversion barrier turns out to equal 0.25 eV. These results are in good agreement with earlier calculations
The multiconfiguration Dirac-Hartree-Fock model has been employed to calculate the expectation values for the hyperfine splittings of the 5d96s22D3/2 and 5d96s22D5/2 levels of atomic gold. One-, two-, and three-body electron correlation effects involving all 79 electrons have been included in a systematic manner. The approximation employed in this study is equivalent to a complete-active-space approach. Calculated electric field gradients, together with experimental values of the electric quadrupole hyperfine-structure constants, allow us to extract a nuclear electric quadrupole moment Q(197Au)=521.5(5.0) mb.
Landau, Arie; Haritan, Idan; Kaprálová-Žd'ánská, Petra Ruth; Moiseyev, Nimrod
2016-05-19
Complex eigenvalues, resonances, play an important role in a large variety of fields in physics and chemistry. For example, in cold molecular collision experiments and electron scattering experiments, autoionizing and predissociative metastable resonances are generated. However, the computation of complex resonance requires modifications of standard electronic structure codes and methods, which are not always straightforward, in addition, application of complex codes requires more computational efforts. Here we show how resonance eigenvalues, positions and widths, can be calculated using the standard, widely used, electronic-structure packages. Our method enables the calculations of the complex resonance eigenvalues by using analytical continuation procedures (such as Padé). The key point in our approach is the existence of narrow analytical passages from the real axis to the complex energy plane. In fact, the existence of these analytical passages relies on using finite basis sets. These passages become narrower as the basis set becomes more complete, whereas in the exact limit, these passages to the complex plane are closed. As illustrative numerical examples we calculated the autoionization Feshbach resonances of helium, hydrogen anion, and hydrogen molecule. We show that our results are in an excellent agreement with the results obtained by other theoretical methods and with available experimental results. PMID:26677725
Gidofalvi, Gergely
2014-01-01
Molecule-optimized basis sets, based on approximate natural orbitals, are developed for accelerating the convergence of quantum calculations with strongly correlated (multi-referenced) electrons. We use a low-cost approximate solution of the anti-Hermitian contracted Schr{\\"o}dinger equation (ACSE) for the one- and two-electron reduced density matrices (RDMs) to generate an approximate set of natural orbitals for strongly correlated quantum systems. The natural-orbital basis set is truncated to generate a molecule-optimized basis set whose rank matches that of a standard correlation-consistent basis set optimized for the atoms. We show that basis-set truncation by approximate natural orbitals can be viewed as a one-electron unitary transformation of the Hamiltonian operator and suggest an extension of approximate natural-orbital truncations through two-electron unitary transformations of the Hamiltonian operator, such as those employed in the solution of the ACSE. The molecule-optimized basis set from the ACS...
Calculation of the fine-structure parameters for np 5 n' f configurations of rare-gas atoms
Anisimova, G. P.; Kapel'Kina, E. L.; Semenov, R. I.
2000-12-01
In the single-configuration approximation, fine-structure parameters are calculated semiempirically in the intermediate coupling scheme for the configurations 2p 5 nf( n=4 6) of NeI, 3 p 5 nf( n=4 7) of ArI, and 4 p 54 f of KrI. With the fine-structure parameters obtained, the coefficients of expansion of the wave functions in basis functions of the LS-coupling scheme and the gyromagnetic ratios are calculated. To the authors’ knowledge, analogous data are absent in the literature. The correctness of the fine-structure parameters obtained is confirmed by values of the fine-splitting constant, which is well known for other configurations of rare-gas atoms with an almost filled p shell.
dftatom: A robust and general Schrödinger and Dirac solver for atomic structure calculations
Čertík, Ondřej; Pask, John E.; Vackář, Jiří
2013-07-01
Classification: 2.1. External routines: Numpy (http://www.numpy.org/) and Cython (http://cython.org/) Nature of problem: Solution of the Schrödinger, Dirac, and Kohn-Sham equations of Density Functional Theory for isolated atoms. Solution method: Radial integrations are carried out using a combination of asymptotic forms, Runge-Kutta, and implicit Adams methods. Eigenfunctions are determined by a combination of bisection and perturbation methods. An outward Poisson integration is employed to increase accuracy in the core region. Self-consistent field equations are solved by adaptive linear mixing. Restrictions: Spherical symmetry Unusual features: Radial integrators work for general potentials and meshes. No restriction to Coulombic or self-consistent potentials; no restriction to uniform or exponential meshes. Outward Poisson integration. Fallback to bisection for robustness. Running time: For uranium, non-relativistic density functional calculation execution time is around 0.6 s for 10-6 a.u. accuracy in total energy on an Intel Core i7 1.46 GHz processor.
Present investigations deal with diagnostic X-ray qualities as applied in medical diagnostics. The track structure computer program PARTRAC simulating the coupled photon - electron transport and interactions in complex geometries has been used. Differential cross sections for electrons in gold with primary energies in the range of 100 eV to 100 keV have been derived on the basis of Seltzer model in case of ionizations and excitations with new Hartree-Fock input data. The elastic cross section data were taken from Fink et al. and Riley et al. These cross sections have been inserted into the existing track structure computer program to provide detailed simulations of electron interactions and to describe above mentioned interface effects. The tracks were simulated in a target volume of gold surrounded by water vapour to represent cell tissue and irradiated by 30 keV and 100 keV photons. The simulated geometry was chosen as in experiments using exoelectron emission. The tracks were calculated for comparison in pure water vapour under the same irradiation conditions. (author)
Atomic partial charges on CH3NH3PbI3 from first-principles electronic structure calculations
Madjet, Mohamed E.; El-Mellouhi, Fedwa; Carignano, Marcelo A.; Berdiyorov, Golibjon R.
2016-04-01
We calculated the partial charges in methylammonium (MA) lead-iodide perovskite CH3NH3PbI3 in its different crystalline phases using different first-principles electronic charge partitioning approaches, including the Bader, ChelpG, and density-derived electrostatic and chemical (DDEC) schemes. Among the three charge partitioning methods, the DDEC approach provides chemically intuitive and reliable atomic charges for this material, which consists of a mixture of transition metals, halide ions, and organic molecules. The DDEC charges are also found to be robust against the use of hybrid functionals and/or upon inclusion of spin-orbit coupling or dispersive interactions. We calculated explicitly the atomic charges with a special focus on the dipole moment of the MA molecules within the perovskite structure. The value of the dipole moment of the MA is reduced with respect to the isolated molecule due to charge redistribution involving the inorganic cage. DDEC charges and dipole moment of the organic part remain nearly unchanged upon its rotation within the octahedral cavities. Our findings will be of both fundamental and practical importance, as the accurate and consistent determination of the atomic charges is important in order to understand the average equilibrium distribution of the electrons and to help in the development of force fields for larger scale atomistic simulations to describe static, dynamic, and thermodynamic properties of the material.
Bierón, Jacek; Indelicato, Paul; Jönsson, Per; Pyykkö, Pekka
2009-01-01
The multiconfiguration Dirac-Hartree-Fock (MCDHF) model has been employed to calculate the expectation values for the hyperfine splittings of the 5d96s2 2D3/2 and 5d96s2 2D5/2 levels of atomic gold. One-, two-, and three-body electron correlation effects involving all 79 electrons have been included in a systematic manner. The approximation employed in this study is equivalent to a Complete Active Space (CAS) approach. Calculated electric field gradients, together with experimental values of the electric quadrupole hyperfine structure constants, allow us to extract a nuclear electric quadrupole moment Q(197Au)=521.5(5.0) mb.
Full text: We have briefly reported the activities in the study of atomic structure and dynamics calculations using the GRASP (General purpose Atomic Structure Program) family codes and also of some activities in NIFS. We have introduced the following items: 1. Analysis of Visible M1 Lines in Tungsten Ions, 2. Collisional-radiative model for W ions, 3. Code development for single electron capture by H nucleus from metal surface, 3. Kα radiation from low charge chlorine heated by an ion beam for plasma diagnostics, 4. Code Availability. And we have summarized the talk. The visible lights emitted from highly charged tungsten ions are of special interest; they are mainly realized by magnetic dipole transitions between the fine structure levels of ions and therefore they suffer less self-absorption by surrounding plasmas providing us with a great advantage for the diagnostics of core plasmas. We have studied the source of these visible line emissions in terms of accurate non-empirical atomic structure calculations and of population kinematics analyses. The GRASP package for multi-configuration Dirac-Fock atomic structure calculation have provided us with the transition energies that are accurate enough to discriminate the states and the charges of tungsten ions. We tried, further, to reproduce the emission spectra by modifying a population kinetics code that has been used for the analysis of boron like ion plasmas. We have assigned new visible emission lines. We have shown that the use of the GRASP family of codes is quite effective for the spectroscopy of the visible line emissions of tungsten ions and for the diagnostics of the MCF plasmas. We have extended a CR model calculation code to include the case of tungsten ions plasmas; the atomic data have been calculated using HULLAC code. We have applied the code to the analysis of W35+ to W37+ plasmas. The code is now still under development. We have now under the development of a code for single electron capture by H
Relativistic calculations for many electron atoms
Many improvements have now been introduced in ab-initio methods for relativistic atomic structure calculations. After a short description of the different methods, we review the various contributions to energy levels and compare the most recent theoretical and experimental results for few electron heavy ions
This book is written to teach atom structure in very easy way. It is divided into nine chapters, which indicates what is the components of matter? when we divide matter continuously, it becomes atom, what did atom look like? particles comprised of matter is not only atom, discover of particles comprised of atom, symbol of element, various radiation, form alchemy to nuclear transmutation, shape of atom is evolving. It also has various pictures in each chapters to explain easily.
Lim, Suh Yeon
2006-08-15
This book is written to teach atom structure in very easy way. It is divided into nine chapters, which indicates what is the components of matter? when we divide matter continuously, it becomes atom, what did atom look like? particles comprised of matter is not only atom, discover of particles comprised of atom, symbol of element, various radiation, form alchemy to nuclear transmutation, shape of atom is evolving. It also has various pictures in each chapters to explain easily.
Atomic Structure Theory Lectures on Atomic Physics
Johnson, Walter R
2007-01-01
Atomic Structure Theory is a textbook for students with a background in quantum mechanics. The text is designed to give hands-on experience with atomic structure calculations. Material covered includes angular momentum methods, the central field Schrödinger and Dirac equations, Hartree-Fock and Dirac-Hartree-Fock equations, multiplet structure, hyperfine structure, the isotope shift, dipole and multipole transitions, basic many-body perturbation theory, configuration interaction, and correlation corrections to matrix elements. Numerical methods for solving the Schrödinger and Dirac eigenvalue problems and the (Dirac)-Hartree-Fock equations are given as well. B-spline basis sets are used to carry out sums arising in higher-order many-body calculations. Illustrative problems are provided, together with solutions. FORTRAN programs implementing the numerical methods in the text are included.
Bierón, Jacek; Froese Fischer, Charlotte; Indelicato, Paul; Jönsson, Per; Pyykkö, Pekka
2009-01-01
The multiconfiguration Dirac-Hartree-Fock (MCDHF) model has been employed to calculate the expectation values for the hyperfine splittings of the 5d96s2 2D3/2 and 5d96s2 2D5/2 levels of atomic gold. One-, two-, and three-body electron correlation effects involving all 79 electrons have been included in a systematic manner. The approximation employed in this study is equivalent to a Complete Active Space (CAS) approach. Calculated electric field gradients, together with experimental values of ...
Atomic physics: computer calculations and theoretical analysis
Drukarev, E. G.
2004-01-01
It is demonstrated, how the theoretical analysis preceding the numerical calculations helps to calculate the energy of the ground state of helium atom, and enables to avoid qualitative errors in the calculations of the characteristics of the double photoionization.
Atomic and molecular structure
This book is a textbook for an introductory course of atomic physics for students of chemistry. After an introduction to the mathematical and physical foundations the quantum mechanical theory of atoms is described starting from simple examples of quantum mechanics. Then the atomic structure and the chemical bending are extensively discussed. This book is also suited for physicists who are especially interested in the atomic structure and the theory of chemical reactions. (HSI)
Electronic structure interpolation via atomic orbitals
Chen, Mohan; Guo, G-C; HE, LIXIN
2010-01-01
We present an efficient scheme for accurate electronic structure interpolations based on the systematically improvable optimized atomic orbitals. The atomic orbitals are generated by minimizing the spillage value between the atomic basis calculations and the converged plane wave basis calculations on some coarse $k$-point grid. They are then used to calculate the band structure of the full Brillouin zone using the linear combination of atomic orbitals (LCAO) algorithms. We find that usually 1...
Calculated Atomic Volumes of the Actinide Metals
Skriver, H.; Andersen, O. K.; Johansson, B.
1979-01-01
The equilibrium atomic volume is calculated for the actinide metals. It is possible to account for the localization of the 5f electrons taking place in americium.......The equilibrium atomic volume is calculated for the actinide metals. It is possible to account for the localization of the 5f electrons taking place in americium....
Lee, Chi-Cheng; Fleurence, Antoine; Yamada-Takamura, Yukiko; Ozaki, Taisuke; Friedlein, Rainer
2014-01-01
So far, it represents a challenging task to reproduce angle-resolved photoelectron (ARPES) spectra of epitaxial silicene by first-principles calculations. Here, we report on the resolution of the previously controversial issue related to the structural configuration of silicene on the ZrB$_2$(0001) surface and its band structure. In particular, by representing the band structure in a large Brillouin zone associated with a single Si atom, it is found that the imaginary part of the one-particle...
Calculations of atomic multiplets across the periodic table
Zhang, Qian
2014-01-01
The goal of this thesis is to develop a simulation tool for the calculation and visualization of the multiplet structure of all atoms across the periodic table. The starting point are self- consistent calculations in the spherical-potential approximation. For the resulting atomic levels we calculate the ab-initio Slater parameters that define the electron-electron repulsion term of the many-body Hamiltonian. We then construct the eigen-states of the Hamiltonian on an open shell by constructin...
Electronic structure interpolation via atomic orbitals
Chen Mohan; Guo, G-C; He Lixin, E-mail: helx@ustc.edu.cn [Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei, 230026 (China)
2011-08-17
We present an efficient scheme for accurate electronic structure interpolation based on systematically improvable optimized atomic orbitals. The atomic orbitals are generated by minimizing the spillage value between the atomic basis calculations and the converged plane wave basis calculations on some coarse k-point grid. They are then used to calculate the band structure of the full Brillouin zone using the linear combination of atomic orbitals algorithms. We find that usually 16-25 orbitals per atom can give an accuracy of about 10 meV compared to the full ab initio calculations, and the accuracy can be systematically improved by using more atomic orbitals. The scheme is easy to implement and robust, and works equally well for metallic systems and systems with complicated band structures. Furthermore, the atomic orbitals have much better transferability than Shirley's basis and Wannier functions, which is very useful for perturbation calculations.
Our relativistic configuration-interaction (RCI) methodology has been extended to multireference cases, and improved to permit the construction of angular-momentum functions of arbitrary size, and to minimize the number of vectors needed with each configuration. We report RCI calculations on the fine (fs) and hyperfine (hfs) structure for the (d+s)3 J=0.5 and 1.5 levels of Zr II. The average fs error is 0.075 eV, and 17% for hfs, when compared to available experiment. These results indicate that it is possible to correctly position all levels of (d+s)n configurations in the transition-metal atoms
Bork, Nicolai Christian; Eurenius, K. E. J.; Rossmeisl, Jan; Knee, C. S.; Vegge, Tejs
2012-01-01
A combined density functional theory and Fourier transform infrared spectroscopy study of the structure and specific site preference of protons and hydrides in the pyrochlore Sm1.92Ca0.08Sn2O7-δ is presented. Two protonic sites of particular high stability are identified, both located on O(1......) oxygen atoms closely associated with a Ca dopant. Further, the unexpected presence of Ho hydride defects in undoped, oxygen deficient Sm2Sn2O7 is reported. Finally, the stretching frequencies and relative intensities for these and other sites are calculated. The main features of the Fourier transform...
Emissivity: A Program for Atomic Emissivity Calculations
Sochi, Taha
2009-01-01
In this article we report the release of a new program for calculating the emissivity of atomic transitions. The program, which can be obtained with its documentation from our website www.scienceware.net, passed various rigorous tests and was used by the author to generate theoretical data and analyze observational data. It is particularly useful for investigating atomic transition lines in astronomical context as the program is capable of generating a huge amount of theoretical data and comp...
Systematic Calculations of Total Atomic Binding Energies
We have calculated total atomic binding energies of 3- to 91-electron ions of all atoms with Z=3 to 118, in the Dirac-Fock model, for applications to atomic mass determination from highly-charged ions. In this process we have determined the ground-state configuration of many ions for which it was not known. We also provide total electronic correlation including Breit correlation for iso-electronic series of beryllium, neon, magnesium and argon, using the multiconfiguration Dirac-Fock approach.
Languages for structural calculations
The differences between human and computing languages are recalled. It is argued that they are to some extent structured in antagonistic ways. Languages in structural calculation, in the past, present, and future, are considered. The contribution of artificial intelligence is stressed
Theory and calculation of the atomic photoeffect
Sabbatucci, Lorenzo; Salvat, Francesc
2016-04-01
The so-called elementary theory of the atomic photoeffect is presented in a form that is suited for practical numerical calculation of subshell cross sections and angular distributions of emitted photoelectrons. Atomic states are described within the independent-electron approximation, with bound and free one-electron orbitals that are solutions of the Dirac equation with the Dirac-Hartree-Fock-Slater self-consistent potential of the ground-state configuration. Detailed derivations are given of subshell cross sections for both excitation to discrete bound levels and ionization. In the case of ionization, the cross section differential in the direction of the photoelectron is obtained for partially polarized photons, with the polarization specified by means of the Stokes parameters. The theoretical formulas have been implemented in a computer program named PHOTACS that calculates tables of excitation and ionization cross sections for any element and subshell. Numerical calculations are practicable for excitations to final states with the principal quantum number up to about 20 and for ionization by photons with energy up to about 2 MeV. Elaborate extrapolation schemes for determining the subshell cross section for excitation to bound levels with larger principal quantum numbers and for ionization by photons with higher energies are described. The effect of the finite width of atomic energy levels is accounted for by convolving the calculated subshell cross-section with a Lorentzian profile.
Atomic structure and electron correlations
Synchrotron experiments combined with theoretical calculations have already given much information on atomic structure and the effects of electron correlations, and this combination of theory and experiment is expected to yield much new information in coming years. In the calculations of photoabsorption cross sections, it is almost always necessary to include electron correlations in both initial and final states to obtain good agreement with experiment. The main theoretical approaches which include effects of electron correlations have been R-matrix theory, random phase approximation with exchange (RPAE), relativistic random phase approximation with exchange, and many-body perturbation theory
Calculations of electron screening in muonic atoms
The electron screening in mounic atoms (O, Al, Fe, In, Ho, Au, Th) has been calculated for p3/2, d5/2 and f7/2 levels with nμ=3/2, d5/2 and f7/2 muons up to nμ=30. Screening corrections are also given for electron configurations with holes in the K and L3 shell. (orig.)
The effects of the 4f shell of electrons and the relativity of valence electrons are compared. The effect of 4f shell (lanthanide contraction) is estimated from the numerical Hartree-Fock (HF) calculations of pseudo-atoms corresponding to Hf, Re, Au, Hg, Tl, Pb and Bi without 4f electrons and with atomic numbers reduced by 14. The relativistic effect estimated from the numerical Dirac-Hartree-Fock (DHF) calculations of those atoms is comparable in the magnitude with that of the 4f shell of electrons. Both are larger for 6s than for 5d or 6p electrons. The various relativistic effects on valence electrons are discussed in detail to determine the proper level of the approximation for the valence electron calculations of systems with heavy elements. An effective core potential system has been developed for heavy atoms in which relativistic effects are included in the effective potentials
Lee, Y.S.
1977-11-01
The effects of the 4f shell of electrons and the relativity of valence electrons are compared. The effect of 4f shell (lanthanide contraction) is estimated from the numerical Hartree-Fock (HF) calculations of pseudo-atoms corresponding to Hf, Re, Au, Hg, Tl, Pb and Bi without 4f electrons and with atomic numbers reduced by 14. The relativistic effect estimated from the numerical Dirac-Hartree-Fock (DHF) calculations of those atoms is comparable in the magnitude with that of the 4f shell of electrons. Both are larger for 6s than for 5d or 6p electrons. The various relativistic effects on valence electrons are discussed in detail to determine the proper level of the approximation for the valence electron calculations of systems with heavy elements. An effective core potential system has been developed for heavy atoms in which relativistic effects are included in the effective potentials.
Million atom DFT calculations using coarse graining and petascale computing
Nicholson, Don; Odbadrakh, Kh.; Samolyuk, G. D.; Stoller, R. E.; Zhang, X. G.; Stocks, G. M.
2014-03-01
Researchers performing classical Molecular Dynamics (MD) on defect structures often find it necessary to use millions of atoms in their models. It would be useful to perform density functional calculations on these large configurations in order to observe electron-based properties such as local charge and spin and the Helmann-Feynman forces on the atoms. The great number of atoms usually requires that a subset be ``carved'' from the configuration and terminated in a less that satisfactory manner, e.g. free space or inappropriate periodic boundary conditions. Coarse graining based on the Locally Self-consistent Multiple Scattering method (LSMS) and petascale computing can circumvent this problem by treating the whole system but dividing the atoms into two groups. In Coarse Grained LSMS (CG-LSMS) one group of atoms has its charge and scattering determined prescriptively based on neighboring atoms while the remaining group of atoms have their charge and scattering determined according to DFT as implemented in the LSMS. The method will be demonstrated for a one-million-atom model of a displacement cascade in Fe for which 24,130 atoms are treated with full DFT and the remaining atoms are treated prescriptively. Work supported as part of Center for Defect Physics, an Energy Frontier Research Center funded by the U.S. DOE, Office of Science, Basic Energy Sciences, used Oak Ridge Leadership Computing Facility, Oak Ridge National Lab, of DOE Office of Science.
BOOK REVIEW: Computational Atomic Structure
Post, Douglass E.
1998-02-01
The primary purpose of `Computational Atomic Structure' is to give a potential user of the Multi-Configuration Hartree-Fock (MCHF) Atomic Structure Package an outline of the physics and computational methods in the package, guidance on how to use the package, and information on how to interpret and use the computational results. The book is successful in all three aspects. In addition, the book provides a good overview and review of the physics of atomic structure that would be useful to the plasma physicist interested in refreshing his knowledge of atomic structure and quantum mechanics. While most of the subjects are covered in greater detail in other sources, the book is reasonably self-contained, and, in most cases, the reader can understand the basic material without recourse to other sources. The MCHF package is the standard package for computing atomic structure and wavefunctions for single or multielectron ions and atoms. It is available from a number of ftp sites. When the code was originally written in FORTRAN 77, it could only be run on large mainframes. With the advances in computer technology, the suite of codes can now be compiled and run on present day workstations and personal computers and is thus available for use by any physicist, even those with extremely modest computing resources. Sample calculations in interactive mode are included in the book to illustrate the input needed for the code, what types of results and information the code can produce, and whether the user has installed the code correctly. The user can also specify the calculational level, from simple Hartree-Fock to multiconfiguration Hartree-Fock. The MCHF method begins by finding approximate wavefunctions for the bound states of an atomic system. This involves minimizing the energy of the bound state using a variational technique. Once the wavefunctions have been determined, other atomic properties, such as the transition rates, can be determined. The book begins with an
AJAC: Atomic data calculation tool in Python
In this work, new features and extensions of a currently used online atomic database management system are reported. A multiplatform flexible computation package is added to the present system, to allow the calculation of various atomic radiative and collisional processes, based on simplifying the use of some existing atomic codes adopted from the literature. The interaction between users and data is facilitated by a rather extensive Python graphical user interface working online and could be installed in personal computers of different classes. In particular, this study gives an overview of the use of one model of the package models (i.e., electron impact collisional excitation model). The accuracy of computing capability of the electron impact collisional excitation in the adopted model, which follows the distorted wave approximation approach, is enhanced by implementing the Dirac R-matrix approximation approach. The validity and utility of this approach are presented through a comparison of the current computed results with earlier available theoretical and experimental results. Finally, the source code is made available under the general public license and being distributed freely in the hope that it will be useful to a wide community of laboratory and astrophysical plasma diagnostics. (interdisciplinary physics and related areas of science and technology)
Updated Atomic Data and Calculations for X-ray Spectroscopy
Foster, A R; Smith, R K; Brickhouse, N S
2012-01-01
We describe the latest release of AtomDB, version 2.0.2, a database of atomic data and a plasma modeling code with a focus on X-ray astronomy. This release includes several major updates to the fundamental atomic structure and process data held within AtomDB, incorporating new ionization balance data, state-selective recombination data, and updated collisional excitation data for many ions, including the iron L-shell ions from Fe$^{+16}$ to Fe$^{+23}$ and all of the hydrogen- and helium-like sequences. We also describe some of the effects that these changes have on calculated emission and diagnostic line ratios, such as changes in the temperature implied by the He-like G-ratios of up to a factor of 2.
Large-scale atomic calculations using variational methods
Joensson, Per
1995-01-01
Atomic properties, such as radiative lifetimes, hyperfine structures and isotope shift, have been studied both theoretically and experimentally. Computer programs which calculate these properties from multiconfiguration Hartree-Fock (MCHF) and configuration interaction (CI) wave functions have been developed and tested. To study relativistic effects, a program which calculates hyperfine structures from multiconfiguration Dirac-Fock (MCDF) wave functions has also been written. A new method of dealing with radial non-orthogonalities in transition matrix elements has been investigated. This method allows two separate orbital sets to be used for the initial and final states, respectively. It is shown that, once the usual orthogonality restrictions have been overcome, systematic MCHF calculations are able to predict oscillator strengths in light atoms with high accuracy. In connection with recent high-power laser experiments, time-dependent calculations of the atomic response to intense laser fields have been performed. Using the frozen-core approximation, where the atom is modeled as an active electron moving in the average field of the core electrons and the nucleus, the active electron has been propagated in time under the influence of the laser field. Radiative lifetimes and hyperfine structures of excited states in sodium and silver have been experimentally determined using time-resolved laser spectroscopy. By recording the fluorescence light decay following laser excitation in the vacuum ultraviolet spectral region, the radiative lifetimes and hyperfine structures of the 7p{sup 2}P states in silver have been measured. The delayed-coincidence technique has been used to make very accurate measurements of the radiative lifetimes and hyperfine structures of the lowest 2P states in sodium and silver. 77 refs, 2 figs, 14 tabs.
Large-scale atomic calculations using variational methods
Atomic properties, such as radiative lifetimes, hyperfine structures and isotope shift, have been studied both theoretically and experimentally. Computer programs which calculate these properties from multiconfiguration Hartree-Fock (MCHF) and configuration interaction (CI) wave functions have been developed and tested. To study relativistic effects, a program which calculates hyperfine structures from multiconfiguration Dirac-Fock (MCDF) wave functions has also been written. A new method of dealing with radial non-orthogonalities in transition matrix elements has been investigated. This method allows two separate orbital sets to be used for the initial and final states, respectively. It is shown that, once the usual orthogonality restrictions have been overcome, systematic MCHF calculations are able to predict oscillator strengths in light atoms with high accuracy. In connection with recent high-power laser experiments, time-dependent calculations of the atomic response to intense laser fields have been performed. Using the frozen-core approximation, where the atom is modeled as an active electron moving in the average field of the core electrons and the nucleus, the active electron has been propagated in time under the influence of the laser field. Radiative lifetimes and hyperfine structures of excited states in sodium and silver have been experimentally determined using time-resolved laser spectroscopy. By recording the fluorescence light decay following laser excitation in the vacuum ultraviolet spectral region, the radiative lifetimes and hyperfine structures of the 7p2P states in silver have been measured. The delayed-coincidence technique has been used to make very accurate measurements of the radiative lifetimes and hyperfine structures of the lowest 2P states in sodium and silver. 77 refs, 2 figs, 14 tabs
Atomistic electronic structure calculations are performed to study the coherent inter-dot couplings of the electronic states in a single InGaAs quantum dot molecule. The experimentally observed excitonic spectrum by Krenner et al (2005) Phys. Rev. Lett. 94 057402 is quantitatively reproduced, and the correct energy states are identified based on a previously validated atomistic tight binding model. The extended devices are represented explicitly in space with 15-million-atom structures. An excited state spectroscopy technique is applied where the externally applied electric field is swept to probe the ladder of the electronic energy levels (electron or hole) of one quantum dot through anti-crossings with the energy levels of the other quantum dot in a two-quantum-dot molecule. This technique can be used to estimate the spatial electron-hole spacing inside the quantum dot molecule as well as to reverse engineer quantum dot geometry parameters such as the quantum dot separation. Crystal-deformation-induced piezoelectric effects have been discussed in the literature as minor perturbations lifting degeneracies of the electron excited (P and D) states, thus affecting polarization alignment of wavefunction lobes for III-V heterostructures such as single InAs/GaAs quantum dots. In contrast, this work demonstrates the crucial importance of piezoelectricity to resolve the symmetries and energies of the excited states through matching the experimentally measured spectrum in an InGaAs quantum dot molecule under the influence of an electric field. Both linear and quadratic piezoelectric effects are studied for the first time for a quantum dot molecule and demonstrated to be indeed important. The net piezoelectric contribution is found to be critical in determining the correct energy spectrum, which is in contrast to recent studies reporting vanishing net piezoelectric contributions.
Structure and theoretical calculations of clay minerals
Structural and spectroscopic methods are combined to determine the full structure, including hydrogen atom positions, of dickite, which is a member of the kaolin group. Using the structural information obtained, quantum chemical calculations are performed on these kaolin group minerals. Special emphasis is laid on the relationship between the experimentally derived structure and theory. Finally, the application of quantum chemical methods to study clay minerals at several levels of approximation is reviewed
Fisher Information and Atomic Structure
Chatzisavvas, K Ch; Panos, C P; Moustakidis, Ch C
2013-01-01
We present a comparative study of several information and statistical complexity measures in order to examine a possible correlation with certain experimental properties of atomic structure. Comparisons are also carryed out quantitatively using Pearson correlation coefficient. In particular, we show that Fisher information in momentum space is very sensitive to shell effects, and is directly associated with some of the most characteristic atomic properties, such as atomic radius, ionization energy, electronegativity, and atomic dipole polarizability. Finally we present a relation that emerges between Fisher information and the second moment of the probability distribution in momentum space i.e. an energy functional of interest in (e,2e) experiments.
We present the results of calculations of surface relaxations, rumplings, energetics, optical band gaps, and charge distribution for the SrZrO3 and PbZrO3 (001) and (011) surfaces using the ab initio code CRYSTAL and a hybrid description of exchange and correlation. We consider both SrO(PbO) and ZrO2 terminations of the (001) surface and Sr(Pb), ZrO, and O terminations of the polar SrZrO3 and PbZrO3 (011) surfaces. On the (001) surfaces, we find that all upper and third layer atoms relax inward, while outward relaxations of all atoms in the second layer are found with the sole exception of the SrO-terminated SrZrO3 (001) surface second layer O atom. Between all (001) and (011) surfaces the largest relaxations, more than 15% of the bulk lattice constant, are for the Sr- and Pb-terminated SrZrO3 and PbZrO3 (011) surface upper layer Sr and Pb atoms. Our calculated surface rumpling for the SrO- and PbO-terminated SrZrO3 and PbZrO3 (001) surfaces (6.77 and 3.32% of a0) are by a factor of ten larger than the surface rumpling for the ZrO2-terminated (001) surfaces (-0.72 and 0.38% of a0, respectively). In contrast to the surface rumpling, the (001) surface energies are comparable and in the energy range from 0.93 eV/cell for the ZrO2-terminated PbZrO3 surface to 1.24 eV/cell for the ZrO2-terminated SrZrO3 surface. In contrast to the (001) surface, different terminations of the SrZrO3 and PbZrO3 (011) surfaces lead to very different surface energies ranging from 1.74 eV/cell for the Pb-terminated PbZrO3 (011) surface to 3.61 eV/cell for the ZrO-terminated SrZrO3 (011) surface. All our calculated (011) surface energies are considerably larger than the (001) surface energies. Our calculated optical band gap for the SrZrO3 bulk, 5.31 eV, is in fair agreement with the experimental value of 5.6 eV. All our calculated optical band gaps for the SrZrO3 and PbZrO3 (001) and (011) surfaces are reduced with respect to the bulk. We predict a considerable increase in the Zr-O chemical
Analytical relativistic self-consistent-field calculations for atoms
A new second-order representation of the Dirac equation is presented. This representation which is exact for a hydrogen atom is applied to approximate analytical self-consistent-field calculations for atoms. Results are given for the rare-gas atoms from helium to radon and for lead. The results compare favorably with numerical Dirac-Hartree-Fock solutions
Atomic physics processes in radial transport calculations
These lectures were intended as preparation for detailed discussions of the role of atomic and molecular physics in confinement research at the 1982 NATO Advanced Study Institute. They begin with a description of the major approaches to magnetic confinement: tandem (ambipolar) mirrors with their associated auxiliary barriers, tokamaks, and stellarators. The leading alternatives, the ELMO Bumpy Torus and the reversed field pinch, are also treated. The evolution equations for particle, energy, and (where relevant) field diffusion are presented and discussed. This is the context for atomic and molecular processes relevant to confinement
Calculation of Al-Zn diagram from central atoms model
无
1999-01-01
A slightly modified central atoms model was proposed. The probabilities of various clusters with the central atoms and their nearest neighboring shells can be calculated neglecting the assumption of the param eter of energy in the central atoms model in proportion to the number of other atoms i (referred with the central atom). A parameter Pα is proposed in this model, which equals to reciprocal of activity coefficient of a component, therefore, the new model can be understood easily. By this model, the Al-Zn phase diagram and its thermodynamic properties were calculated, the results coincide with the experimental data.
Introduction to Density Functional Theory: Calculations by Hand on the Helium Atom
Baseden, Kyle A.; Tye, Jesse W.
2014-01-01
Density functional theory (DFT) is a type of electronic structure calculation that has rapidly gained popularity. In this article, we provide a step-by-step demonstration of a DFT calculation by hand on the helium atom using Slater's X-Alpha exchange functional on a single Gaussian-type orbital to represent the atomic wave function. This DFT…
Windows based application has been developed for the calculation of atomic hyperfine spectrum of odd isotopes keeping in view of the needs of the atomic spectroscopists. The application can also calculate the hyperfine spectrum of another odd isotope if hyperfine structure constants of one isotope are known. Various features of the developed application are discussed. (author)
Arithmetic aspects of atomic structures
The purpose of this presentation is to give an updated account of some on-going work related to the theory of large atoms in the context of large-Z asymptotics, which began over the last decade or so. The goal of that work is to produce a refined version of the Thomas-Fermi theory that accounts for observed physical features such as electronic orbitals or an atomic shell structure. This presents similarities with issues in quantum chaos. (orig.)
Multilevel domain decomposition for electronic structure calculations
We introduce a new multilevel domain decomposition method (MDD) for electronic structure calculations within semi-empirical and density functional theory (DFT) frameworks. This method iterates between local fine solvers and global coarse solvers, in the spirit of domain decomposition methods. Using this approach, calculations have been successfully performed on several linear polymer chains containing up to 40,000 atoms and 200,000 atomic orbitals. Both the computational cost and the memory requirement scale linearly with the number of atoms. Additional speed-up can easily be obtained by parallelization. We show that this domain decomposition method outperforms the density matrix minimization (DMM) method for poor initial guesses. Our method provides an efficient preconditioner for DMM and other linear scaling methods, variational in nature, such as the orbital minimization (OM) procedure
Universal bosonic tetramers of dimer-atom-atom structure
Deltuva, A.
2012-01-01
Unstable four-boson states having an approximate dimer-atom-atom structure are studied using momentum-space integral equations for the four-particle transition operators. For a given Efimov trimer the universal properties of the lowest associated tetramer are determined. The impact of this tetramer on the atom-trimer and dimer-dimer collisions is analyzed. The reliability of the three-body dimer-atom-atom model is studied.
Theoretical Calculation of Absolute Radii of Atoms and Ions. Part 1. The Atomic Radii
Raka Biswas
2002-02-01
Full Text Available Abstract. A set of theoretical atomic radii corresponding to the principal maximum in the radial distribution function, 4ÃÂ€r2R2 for the outermost orbital has been calculated for the ground state of 103 elements of the periodic table using Slater orbitals. The set of theoretical radii are found to reproduce the periodic law and the Lother MeyerÃ¢Â€Â™s atomic volume curve and reproduce the expected vertical and horizontal trend of variation in atomic size in the periodic table. The d-block and f-block contractions are distinct in the calculated sizes. The computed sizes qualitatively correlate with the absolute size dependent properties like ionization potentials and electronegativity of elements. The radii are used to calculate a number of size dependent periodic physical properties of isolated atoms viz., the diamagnetic part of the atomic susceptibility, atomic polarizability and the chemical hardness. The calculated global hardness and atomic polarizability of a number of atoms are found to be close to the available experimental values and the profiles of the physical properties computed in terms of the theoretical atomic radii exhibit their inherent periodicity. A simple method of computing the absolute size of atoms has been explored and a large body of known material has been brought together to reveal how many different properties correlate with atomic size.
Starting SCF Calculations by Superposition of Atomic Densities
van Lenthe, J.H.; Zwaans, R.; van Dam, H.J.J.; Guest, M.F.
2006-01-01
We describe the procedure to start an SCF calculation of the general type from a sum of atomic electron densities, as implemented in GAMESS-UK. Although the procedure is well-known for closed-shell calculations and was already suggested when the Direct SCF procedure was proposed, the general procedu
SUPERSTRUCTURE - AN ATOMIC STRUCTURE CODE
Eissner, W.
1991-01-01
We summarize the properties of the atomic structure code SUPERSTRUCTURE, which yields bound state energies in LS coupling and intermediate coupling as well as associated radiative data. Other data that can be computed include term coupling coefficients and radiative data with allowance for cascading. Results are given, mainly for members of the Be isoelectronic sequence, to demonstrate the power and range of the code. Other examples deal with "forbidden" transitions in N-like and He-like ions.
Jiang, Jun; Cheng, Yongjun; Bromley, M W J
2014-01-01
Effective oscillator strength distributions are systematically generated and tabulated for the alkali atoms, the alkaline-earth atoms, the alkaline-earth ions, the rare gases and some miscellaneous atoms. These effective distributions are used to compute the dipole, quadrupole and octupole static polarizabilities, and are then applied to the calculation of the dynamic polarizabilities at imaginary frequencies. These polarizabilities can be used to determine the long-range $C_6$, $C_8$ and $C_{10}$ atom-atom interactions for the dimers formed from any of these atoms and ions, and we present tables covering all of these combinations.
A Variational Monte Carlo Approach to Atomic Structure
Davis, Stephen L.
2007-01-01
The practicality and usefulness of variational Monte Carlo calculations to atomic structure are demonstrated. It is found to succeed in quantitatively illustrating electron shielding, effective nuclear charge, l-dependence of the orbital energies, and singlet-tripetenergy splitting and ionization energy trends in atomic structure theory.
Calculating trajectories for atoms in near-resonant lightfields
We review several methods for calculating the time development of the internal state and the external motion of atoms in near-resonant light fields, with emphasis on studying the focussing of atomic beams into microscopic and potentially nanoscopic patterns. Three different approaches are considered: two-level semiclassical, multi-level semiclassical, and the Monte Carlo wavefunction method. The two-level semiclassical technique of McClelland and Scheinfein (1991) and McClelland (1995) is extended to three dimensions, and used to calculate the trajectories of atoms and the imaging properties of a simple lens formed from a near-resonant travelling TEM01 mode laser. The model is then extended to multi-level atoms, where we calculate the density matrix for the internal state of a sample of thermal atoms in a standing wave, and show how cooling processes can be simulated. Finally, we use the Monte Carlo wavefunction method to calculate the internal state of the atom, and compare the results and required computation time to those of the multi-level semiclassical technique. (authors)
Improved density functional calculations for atoms, molecules and surfaces
The non-collinear and collinear descriptions within relativistic density functional theory is described. We present results of both non-collinear and collinear calculations for atoms, diatomic molecules, and some surface simulations. We find that the accuracy of our density functional calculations for the smaller systems is comparable to good quantum chemical calculations, and thus this method provides a sound basis for larger systems where no such comparison is possible. (author)
Update on nuclear structure effects in light muonic atoms
Hernandez, Oscar Javier; Ji, Chen; Bacca, Sonia; Barnea, Nir
2016-01-01
We present calculations of the nuclear structure corrections to the Lamb shift in light muonic atoms, using state-of-the-art nuclear potentials. We outline updated results on finite nucleon size contributions.
Calculations of optical rotation: Influence of molecular structure
Yu Jia
2012-01-01
Full Text Available Ab initio Hartree-Fock (HF method and Density Functional Theory (DFT were used to calculate the optical rotation of 26 chiral compounds. The effects of theory and basis sets used for calculation, solvents influence on the geometry and values of calculated optical rotation were all discussed. The polarizable continuum model, included in the calculation, did not improve the accuracy effectively, but it was superior to γs. Optical rotation of five or sixmembered of cyclic compound has been calculated and 17 pyrrolidine or piperidine derivatives which were calculated by HF and DFT methods gave acceptable predictions. The nitrogen atom affects the calculation results dramatically, and it is necessary in the molecular structure in order to get an accurate computation result. Namely, when the nitrogen atom was substituted by oxygen atom in the ring, the calculation result deteriorated.
Investigation of atomic processes during film growth using semiempirical calculations
Leonardelli, G
2001-01-01
Growth of thin films on solid surfaces is strongly determined by the rates of the individual atomic processes and therefore depends on the energy barriers which must be surmounted during these processes. The diffusion barriers of interlayer diffusion processes are calculated in this work using embedded atom method (EAM) potentials. Great attention is paid to effects of small simulation cells preventing the atoms near the step edge from relaxing completely and thereby modifying the barriers for step descent on steps of the Pt(111) surface. Calculations in this work can also explain experimental data which show Co atoms sitting in special sites like corners and kinks when small amounts of Co are deposited on the Pt(111) surface. The results show why these sites are occupied and why configurations along A-steps are different from those on B-steps. Furthermore, calculations explain the intermixing of adlayer and substrate atoms on fcc(111) surfaces in the vicinity of rough steps occurring when these steps smoothe...
Variable atomic radii for continuum-solvent electrostatics calculation
Zhou, Baojing; Agarwal, Manish; Wong, Chung F.
2008-07-01
We have developed a method to improve the description of solute cavity defined by the interlocking-sphere model for continuum-solvent electrostatics calculations. Many models choose atomic radii from a finite set of atom types or uses an even smaller set developed by Bondi [J. Phys. Chem. 68, 441 (1964)]. The new model presented here allowed each atom to adapt its radius according to its chemical environment. This was achieved by first approximating the electron density of a molecule by a superposition of atom-centered spherical Gaussian functions. The parameters of the Gaussian functions were then determined by optimizing a function that minimized the difference between the properties from the model and those from ab initio quantum calculations. These properties included the electrostatics potential on molecular surface and the electron density within the core of each atom. The size of each atom was then determined by finding the radius at which the electron density associated with the atom fell to a prechosen value. This value was different for different chemical elements and was chosen such that the averaged radius for each chemical element in a training set of molecules matched its Bondi radius. Thus, our model utilized only a few adjustable parameters—the above density cutoff values for different chemical elements—but had the flexibility of allowing every atom to adapt its radius according to its chemical environment. This variable-radii model gave better solvation energy for 31 small neutral molecules than the Bondi radii did, especially for a quantum mechanics/Poisson-Boltzmann approach we developed earlier. The improvement was most significant for molecules with large dipole moment. Future directions for further improvement are also discussed.
Isogeometric analysis in electronic structure calculations
Cimrman, Robert; Kolman, Radek; Tůma, Miroslav; Vackář, Jiří
2016-01-01
In electronic structure calculations, various material properties can be obtained by means of computing the total energy of a system as well as derivatives of the total energy w.r.t. atomic positions. The derivatives, also known as Hellman-Feynman forces, require, because of practical computational reasons, the discretized charge density and wave functions having continuous second derivatives in the whole solution domain. We describe an application of isogeometric analysis (IGA), a spline modification of finite element method (FEM), to achieve the required continuity. The novelty of our approach is in employing the technique of B\\'ezier extraction to add the IGA capabilities to our FEM based code for ab-initio calculations of electronic states of non-periodic systems within the density-functional framework, built upon the open source finite element package SfePy. We compare FEM and IGA in benchmark problems and several numerical results are presented.
Stability and structure of atomic chains on Si(111)
Battaglia, Corsin; Aebi, Philipp; Erwin, Steven C.
2008-01-01
We study the stability and structure of self-assembled atomic chains on Si(111) induced by monovalent, divalent and trivalent adsorbates, using first-principles total-energy calculations and scanning tunneling microscopy. We find that only structures containing exclusively silicon honeycomb or silicon Seiwatz chains are thermodynamically stable, while mixed configurations, with both honeycomb and Seiwatz chains, may be kinetically stable. The stability and structure of these atomic chains can...
Autoionizing States of Atoms Calculated Using Generalized Sturmians
Avery, James Emil; Avery, John Scales
2005-01-01
The generalized Sturmian method is applied to autoionizing states of atoms and ions. If the Goscinskian basis sets allow for a sufficient amount of angular correletion, the calculated energies of doubly-excited (autoionizing) states are found to agree well with the few available experimental...... energies. A large-Z approximation is discussed, and simple formulas are derived which are valid not only for autoionizing states, but for all states of an isoelectronic atomic series. Diagonalization of a small block of the interelectron repulsion matrix yields roots that can be used for a wide range of Z...
High-accuracy coupled cluster calculations of atomic properties
The four-component Fock-space coupled cluster and intermediate Hamiltonian methods are implemented to evaluate atomic properties. The latter include the spectra of nobelium and lawrencium (elements 102 and 103) in the range 20000-30000 cm−1, the polarizabilities of elements 112-114 and 118, required for estimating their adsorption enthalpies on surfaces used to separate them in accelerators, and the nuclear quadrupole moments of some heavy atoms. The calculations on superheavy elements are supported by the very good agreement with experiment obtained for the lighter homologues
[Atomic force field FFsol for calculation of molecular interactions of in water environment].
Pereiaslavets, L B; Finkel'shtein, A V
2010-01-01
Detailed calculations of protein interactions with explicitly considered water takes enormous computer time. The calculation becomes faster if water is considered implicitly (as a continuous media rather than as molecules); however, these calculations are much less precise, unless one uses an additional (and also volumes) computation of the solvent-accessible areas of protein atoms. The aim of our study was to obtain parameters for non-bonded atom-atom interactions for the case when water surrounding is considered implicitly and the solvent-accessible areas are not computed. Since the "in-vacuum" interactions of atoms are obtained from experimental structures of crystals and enthalpies of their sublimation, the "in-water" interactions of atoms must be corrected using solvation free energies of molecules, which can be obtained from the Henry constants. Taken 58 structures of molecular crystals and thermodynamic data on their sublimation and solubility, we obtained parameters for "in-water" attraction and repulsion of atoms typical of protein structures (H, C, N, O, S) in various covalently-bonded states, as well as parameters for electrostatic interactions. All necessary for calculations parameters of covalent interactions have been taken from the ENCAD force field, and partial charges of all atoms of separate molecules of a crystal have been obtained from quantum-mechanical calculations. The sought parameters of the "in-water" van der Waals and electrostatic interactions were optimized so as to achieve the best description of equilibrium crystal structures and their sublimation and solvation at the room temperature. With the optimized parameters, the average error in calculation of the effective cohesion energy of molecules in crystals was less than 10% both in the "in-vacuum" and "in-water" cases. PMID:20586195
Semiempirical studies of atomic structure
The energy level structure, transition probabilities, and general spectroscopic properties of highly-ionized many-electron systems are studied through the combined use of sensitive semiempirical data systematizations, selected precision experimental measurements, and specialized theoretical computations. Measurements are made primarily through the use of fast ion beam excitation methods, which are combined with available data from laser- and tokamak-produced plasmas, astrophysical sources, and conventional light sources. The experimental studies are strengthened through large scale ab initio calculations. Large blocks of data are predictively systematized and parameterized along isoelectric, homologous, isoionic, Rydberg, and yrast series, to provide a comprehensive and reliable data base
Theoretical Calculation of Absolute Radii of Atoms and Ions. Part 1. The Atomic Radii
Raka Biswas; Dulal C. Ghosh
2002-01-01
Abstract. A set of theoretical atomic radii corresponding to the principal maximum in the radial distribution function, 4ÃÂ€r2R2 for the outermost orbital has been calculated for the ground state of 103 elements of the periodic table using Slater orbitals. The set of theoretical radii are found to reproduce the periodic law and the Lother MeyerÃ¢Â€Â™s atomic volume curve and reproduce the expected vertical and horizontal trend of variation in atomic size in the periodic table. The d-block and...
Description of an atomic structure software package
The MCHF method has been shown to be effective for the study of correlation and for croblems where relativistic effects are small. It has been applied to some difficult cases where a perturber is embedded in a Rydberg series and mixing depends critically on correlation. Relativistic corrections may be included through the Breit-Pauli approximation. This approach has been used recently in a study of the Boron sequence where the fine structure splitting obtained from multiconfiguration Dirac-Hartree-Fock calculations was not in good agreement with observation and a correction process was needed. It has also been used to study forbidden transitions in the Carbon, Nitrogen, and Oxygen sequence. Thus the MCHF (or MCHF+BP) method is a versatile method for the study of atomic structure for a large class of problems. The present software package is based on this method
DFT calculations on atom-specific electronic properties of G/SiC(0001)
Kajihara, M.; Suzuki, T.; Shahed, S. M. F.; Komeda, T.; Minamitani, E.; Watanabe, S.
2016-05-01
We investigate the atom-specific interfacial electronic properties of the epitaxial graphene on Si-terminated SiC substrate using density functional theory (DFT) calculation with van der Waals interaction correction, focusing on the dependency of the local electronic state on the chemical environment. The band structure projected on the respective atomic orbitals of the carbon atoms in the buffer layer and uppermost Si atoms demonstrates that the dangling bonds of these atoms form band structures around the Fermi level. The contribution of each atom to the dangling bond states strongly depends on the chemical environment, i.e., the presence/absence of the interlayer Si-C covalent bond. This difference also affects the atom-specific local density of states of the top-layer graphene through its interaction with the substrate/buffer layer. We demonstrate that the bias voltage dependency of the scanning tunneling spectroscopy (STS) mapping image clearly reflects the presence of the dangling bonds of the buffer layer carbon or uppermost Si atom in the substrate, which would enable the detection of the buried dangling bond with an atomic spatial resolution via STS.
Semiempirical studies of atomic structure
The energy level structure, transition probabilities, and general spectroscopic properties of highly ionized many-electron systems are studied through the combined use of sensitive semiempirical data systematizations, selected precision experimental measurements, and specialized theoretical computations. Measurements are made primarily through the use of fast ion beam excitation methods, which are combined with available data from laser-and tokamak-produced plasmas, astrophysical sources, and conventional light sources. The experimental studies are strengthened through large-scale ab initio calculations. Typical examples are the following: lifetime measurements in the neon isoelectronic sequence; multiplexed decay curve measurements of Li-like Si XII; and isoelectronic specification of intershell resonance and intercombination decay rates using measured transition probabilities and spectroscopically determined singlet-mixing amplitudes
Hyperspherical Three-Body Calculation for Exotic Atoms
Ground state energies of atomic three-body systems like negatively charged hydrogen, normal helium, positively charged-lithium, beryllium, carbon, oxygen, neon and negatively charged exotic-muonium and positronium atoms have been calculated adopting hyperspherical harmonics expansion method. Calculation of matrix elements of two body interactions needed in the hyperspherical harmonics expansion method for a three body system is greatly simplified by expanding the bra-and ket-vector states in the hyperspherical harmonics (HH) basis states appropriate for the partition corresponding to the interacting pair. This involves the Raynal-Revai coefficients (RRC), which are the transformation coefficients between the HH bases corresponding to the two partitions. Use of RRC become particularly essential for the numerical solution of three-body Schroedinger equation where the two-body potentials are other than Coulomb or harmonic. However in the present work the technique is used for two electron atoms 1H-(p+e-e-), D-(d+e-e-), Mu-(μ+e-e-), 4He(4He2+e-e-), 6Li(6Li3+e-e-), 10Be( 10Be4+e-e-), 12C(12C6+e-e-), 16O(16O8+e-e-) etc. and the exotic positronium negative ion Ps-(e+e-e-) where the interactions are purely Coulomb. The relative convergence in ground state binding energy with increasing Kmax for 20Ne has been demonstrated as a representative case. The calculated energies at Kmax = 28 using RRC's have been compared with those obtained by a straight forward manner in some representative cases to demonstrate the appropriateness of the use of RRC. The extrapolated energies have also been compared with those found in the literature. The calculated binding energies agree within the computational error. (author)
Recent developments in high-spin calculations in atomic nuclei
A brief introduction to the recent achievements in the high-spin domain in nuclear physics is given. Results of the calculations in highly developed rotational bands in deformed nuclei, as well as the calculations in the structure of the yrast isomers are presented. The calculations fail in two aspects: local minima in the yrast line are not confirmed experimentally, the overall slope of the yrast line in 152Dy is considerably overestimated. The calculations of the yrast line with new Woods-Saxon parameters are now in progress. The parameters are chosen to reproduce the large gap in the levels at proton number Z=64. (M.H.)
Avery, John Scales; Avery, James Emil; Aquilanti, Vincenzo; Caligiana, Andrea
2004-01-01
The generalized Sturmian method for atomic and molecular electronic structure calculations is a direct configuration interaction method in which the configurations are chosen to be isoenergetic solutions of an approximate N-electron Schrödinger equation with a weighted potential, $\\beta_\
Global nuclear-structure calculations
The revival of interest in nuclear ground-state octupole deformations that occurred in the 1980's was stimulated by observations in 1980 of particularly large deviations between calculated and experimental masses in the Ra region, in a global calculation of nuclear ground-state masses. By minimizing the total potential energy with respect to octupole shape degrees of freedom in addition to ε2 and ε4 used originally, a vastly improved agreement between calculated and experimental masses was obtained. To study the global behavior and interrelationships between other nuclear properties, we calculate nuclear ground-state masses, spins, pairing gaps and Β-decay and half-lives and compare the results to experimental qualities. The calculations are based on the macroscopic-microscopic approach, with the microscopic contributions calculated in a folded-Yukawa single-particle potential
The geometry, electronic structure and magnetic property of the hexagonal AlN (h-AlN) sheet doped by 5d atoms (Lu, Hf, Ta, W, Re, Os, Ir, Pt, Au and Hg) are investigated by first-principles calculations based on the density functional theory. The influence of symmetry and symmetry-breaking is also studied. There are two types of local symmetries of the doped systems: C3v and D3h. The symmetry will deviate from exact C3v and D3h for some particular dopants after optimization. The total magnetic moments of the doped systems are 0μB for Lu, Ta and Ir; 1μB for Hf, W, Pt and Hg; 2μB for Re and Au; and 3μB for Os and Al-vacancy. The total densities of state are presented, where impurity energy levels exist. The impurity energy levels and total magnetic moments can be explained by the splitting of 5d orbitals or molecular orbitals under different symmetries. (condensed matter: structural, mechanical, and thermal properties)
LOCAL ATOMIC STRUCTURE OF AMORPHOUS METALS
Egami, T.; Maed, K.; Srolovitz, D.; Vitek, V.
1980-01-01
The local parameters are introduced to describe the local atomic structure of amorphous metals. They define the structural defects which facilitate the explanation of various properties, including the volume change by annealing.
Kinetic-energy density functional: Atoms and shell structure
We present a nonlocal kinetic-energy functional which includes an anisotropic average of the density through a symmetrization procedure. This functional allows a better description of the nonlocal effects of the electron system. The main consequence of the symmetrization is the appearance of a clear shell structure in the atomic density profiles, obtained after the minimization of the total energy. Although previous results with some of the nonlocal kinetic functionals have given incipient structures for heavy atoms, only our functional shows a clear shell structure for most of the atoms. The atomic total energies have a good agreement with the exact calculations. Discussion of the chemical potential and the first ionization potential in atoms is included. The functional is also extended to spin-polarized systems. copyright 1996 The American Physical Society
Atomic probes of surface structure and dynamics
Progress for the period Sept. 15, 1992 to Sept. 14, 1993 is discussed. Semiclassical methods that will allow much faster and more accurate three-dimensional atom--surface scattering calculations, both elastic and inelastic, are being developed. The scattering of He atoms from buckyballs is being investigated as a test problem. Somewhat more detail is given on studies of He atom scattering from defective Pt surfaces. Molecular dynamics simulations of He+ and Ar+ ion sputtering of Pt surfaces are also being done. He atom scattering from Xe overlayers on metal surfaces and the thermalized dissociation of H2 on Cu(110) are being studied. (R.W.R.) 64 refs
A New Pseudospectral Method for Calculations of Hydrogen Atom in Arbitrary External Fields
QIAO Hao-Xue; LI Bai-Wen1
2002-01-01
A new pseudospectral method was introduced to calculate wavefunctions and energy levels of hydrogen atom in arbitrary potential. Some results of hydrogen atom in uniform magnetic fields were presented, high accuracy of results was obtained with simple calculations, and our calculations show very fast convergence. It suggests a new methodfor calculations of hydrogen atom in external fields.
Ab Initio factorized LCAO calculation of the electronic structure of α-SiO2
The authors report on the results of calculations of the electronic structure of α-quartz that were made using first principles, factorized linear combination of atomic orbitals method. Results were obtained for the primitive 9-atom, and orthorhombic 18- and 72-atom unit cells. Application of this method to the calculation of the electronic structure of the neutral oxygen vacancy in α-quartz is discussed and results obtained using a 72-atom unit cell are given
The calculation of photoionization angular distribution parameter β of atomic Na
刘锦超; 蒲丰年; 刘汉奎; 郭建军; 程延松; 杨向东
1999-01-01
The differential cross sections and angular distribution parameter of the photoionization processes 2p~63s→2p~53skl of atomic Na have been calculated by using many-body perturbation theory. In the calculation, the resonant structure of the excitation process 2s→3p has been included. The electron correlation interaction was analyzed by using the effective diagram method. The summation of specific classes of these diagrams is to an infinite order. The results of calculations are compared with experimental data, which are in good agreement with the experiment.
Introduction to Atomic Structure: Demonstrations and Labs.
Ciparick, Joseph D.
1988-01-01
Demonstrates a variety of electrical phenomena to help explain atomic structure. Topics include: establishing electrical properties, electrochemistry, and electrostatic charges. Recommends demonstration equipment needed and an explanation of each. (MVL)
Electronic structure of crystalline uranium nitride: LCAO DFT calculations
The results of electronic structure calculations performed for the first time for crystalline uranium nitride and using a LCAO basis are discussed. For calculations we used the density functional method with the PW91 exchange correlation potential and a variety of relativistic core potentials for the uranium atom. The calculated atomization energy of the crystal agrees well with the experimental data and with the results of calculations with the plane wave basis. It is shown that a chemical bond in crystalline uranium nitride is a metal covalent bond. The metal component of the bond is due to the 5f electrons localized on the uranium atom and having energies near the Fermi level and the bottom of the conduction band. The covalent component of the chemical bond results from an overlap between the uranium 6d and 7s valence orbitals and the nitrogen 2p atomic orbitals. Inclusion of the 5f electrons in the core of the uranium atom introduces relatively minor changes in the calculated binding energy and electron density distribution
Falkenberg, G.; Bunk, O.; Johnson, R.L.; Rodrigues, J.A.; Takeuchi, N.
2002-01-01
Using scanning-tunneling microscopy (STM) and first-principles total-energy calculations, we have determined the atomic geometry of the superstructures formed by the adsorption of up to 0.5 monolayer of indium on Ge(001) and annealing at temperatures above 200 degreesC. A strong interaction between...... indium adatoms and the germanium substrate atoms leads to the formation of two different In-Ge subunits on the Ge(001) surface. In the subsaturation regime separate (nx4) subunits are observed where n can be either 3 or 4 and the STM images resemble those of the Si(001)-(3x4)-In and -Al reconstructions......, Appl. Surf. Sci. 123/124, 104 (1998) for In on Si(001). For the (4x4) subunit, we propose a model that includes the main features of the (3x4) subunit together with additional mixed Ge-In dimers. The atomic positions were optimized using ab initio total-energy calculations. The calculated local...
Calculation of electron-impace excitation and ionization of atoms
Over the last few years it has been the author's goal to develop a open-quotes completeclose quotes electron-atom scattering theory. For a given projectile energy such a theory would be able to correctly predict the major scattering processes. These include elastic, excitation, and ionization cross sections. The convergent close-coupling (CCC) approach is a candidate for such a theory. Hamiltonians in an orthogonal Laguerre basis. The usage of this basis ensures that open-quotes completenessclose quotes is approached as N is increased. The CCC method may be thought of as a more systematic implementation of pseudostate methods. Whereas the success of the close-coupling approach to the calculation of excitation cross sections has been demonstrated for many decades it has rarely been applied to the calculation of ionization processes. By demonstrating the ability to obtain accurate ionization cross sections simultaneously with excitation processes the CCC method appears to have fulfilled these goals. In the talk a general outline of the CCC method will be given and its recent applications to the measurements of electron-impact ionization and excitation of the ground state of helium discussed
Three-dimensional rf structure calculations
The calculation of three-dimensional rf structures is rapidly approaching adolescence, after having been in its infancy for the last four years. This paper will show the kinds of calculations that are currently being performed in the frequency domain and is a companion paper to one in which time-domain calculations are described. 13 refs., 14 figs
Electronic structure of crystalline uranium nitride: LCAO DFT calculations
The results of the first LCAO DFT calculations of cohesive energy, band structure and charge distribution in uranium nitride (UN) crystal are presented and discussed. The calculations are made with the uranium atom relativistic effective core potentials, including 60, 78 and 81 electrons in the core. It is demonstrated that the chemical bonding in UN crystal has a metallic-covalent nature. Three 5f-electrons are localized on the U atom and occupy the states near the Fermi level. The metallic nature of the crystal is due to the f-character of both the valence-band top and the conduction-band bottom. The covalent bonds are formed by the interaction of 7s- and 6d-states of the uranium atom with the 2p-states of the nitrogen atom. It is shown that the inclusion of 5f-electrons in the atomic core introduces small changes in the calculated cohesive energy of UN crystal and electron charge distribution. However, the inclusion of 5s-, 5p-, 5d-electrons in the valence shell allows the better agreement with the calculated and experimental cohesive-energy value. (copyright 2008 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim) (orig.)
Final disposal room structural response calculations
Finite element calculations have been performed to determine the structural response of waste-filled disposal rooms at the WIPP for a period of 10,000 years after emplacement of the waste. The calculations were performed to generate the porosity surface data for the final set of compliance calculations. The most recent reference data for the stratigraphy, waste characterization, gas generation potential, and nonlinear material response have been brought together for this final set of calculations
The equilibrium lattice parameters, formation energy, and diffusion behavior of hydrogen atoms in vanadium hydrides with and without Ti substitution were calculated by ab initio calculations and quantum correction by zero point energy was achieved using phonon vibration calculations. The calculated formation energies indicated that Ti substitution induces instability in the vanadium hydrides and electron density calculations showed that hydrogen has strong electrochemical affinity with Ti. The diffusion behavior was examined by the nudged elastic band (NEB) method to investigate the transition states of the hydrides. It revealed that Ti substitution is shown to reduce the diffusion coefficient and this effect was decreased with increasing temperature. The results of this study are expected to provide useful guidelines for understanding hydrogen absorption and desorption properties of hydrogen storage materials. - Highlights: • Ti substitution expands any crystal structure of vanadium metal or hydride. • H atoms are repulsed by Ti atoms despite their high electrochemical affinity. • Strong electrochemical bonding between Ti and H lowers the formation energy. • Ti substitution reduces hydrogen diffusion by over 90% at room temperature
Kim, Jiwoong, E-mail: jwk@kigam.re.kr; Yoo, Jeong-Hyun, E-mail: yjh0010@naver.com; Cho, Sung-Wook, E-mail: cho@kigam.re.kr
2014-12-15
The equilibrium lattice parameters, formation energy, and diffusion behavior of hydrogen atoms in vanadium hydrides with and without Ti substitution were calculated by ab initio calculations and quantum correction by zero point energy was achieved using phonon vibration calculations. The calculated formation energies indicated that Ti substitution induces instability in the vanadium hydrides and electron density calculations showed that hydrogen has strong electrochemical affinity with Ti. The diffusion behavior was examined by the nudged elastic band (NEB) method to investigate the transition states of the hydrides. It revealed that Ti substitution is shown to reduce the diffusion coefficient and this effect was decreased with increasing temperature. The results of this study are expected to provide useful guidelines for understanding hydrogen absorption and desorption properties of hydrogen storage materials. - Highlights: • Ti substitution expands any crystal structure of vanadium metal or hydride. • H atoms are repulsed by Ti atoms despite their high electrochemical affinity. • Strong electrochemical bonding between Ti and H lowers the formation energy. • Ti substitution reduces hydrogen diffusion by over 90% at room temperature.
2007-01-01
Density functional calculations of electronic structure, total energy, structural distortions, and magnetism for hydrogenated single-layer, bilayer, and multi-layer graphene are performed. It is found that hydrogen-induced magnetism can survives only at very low concentrations of hydrogen (single-atom regime) whereas hydrogen pairs with optimized structure are usually nonmagnetic. Chemisorption energy as a function of hydrogen concentration is calculated, as well as energy barriers for hydrog...
Boukhvalov, D W; Katsnelson, M. I.; Lichtenstein, A. I.
2008-01-01
Density functional calculations of electronic structure, total energy, structural distortions, and magnetism for hydrogenated single-layer, bilayer, and multi-layer graphene are performed. It is found that hydrogen-induced magnetism can survives only at very low concentrations of hydrogen (single-atom regime) whereas hydrogen pairs with optimized structure are usually nonmagnetic. Chemisorption energy as a function of hydrogen concentration is calculated, as well as energy barriers for hydrog...
Progress in numerical calculations of ion-atom collisions
Reading, J.F.; Ford, A.L.; Becker, R.L.
1983-01-01
An ion-atom collision produces a time dependent perturbation of a many fermion system. In this collision, excitation, ionization and charge transfer can occur. The driving mechanism for these processes may be thought of as the potentials seen by individual electrons at any given separation of the projectile and target nuclei. If we think of these potentials as belonging to the target (a nucleus and electrons) and the projectile (another nucleus and electrons) then as detected by an electron the potentials change because: (a) the target and projectile change position, and (b) electrons on the target and projectile change states. Most work in the past fifty years has concentrated on solving the independent particle model (IPM). Cracks are beginning to appear in this model which only allows for type (a) changes in the potential. But in a short review we shall have quite enough to do in understanding the progress made in the last decade on the IPM. This paper is divided into three parts. The first deals with how to reduce the IPM to the single electron model (SEM). The second is on a new method where charge transfer is important. The third confronts some standard models with modern calculations.
Progress in numerical calculations of ion-atom collisions
An ion-atom collision produces a time dependent perturbation of a many fermion system. In this collision, excitation, ionization and charge transfer can occur. The driving mechanism for these processes may be thought of as the potentials seen by individual electrons at any given separation of the projectile and target nuclei. If we think of these potentials as belonging to the target (a nucleus and electrons) and the projectile (another nucleus and electrons) then as detected by an electron the potentials change because: (a) the target and projectile change position, and (b) electrons on the target and projectile change states. Most work in the past fifty years has concentrated on solving the independent particle model (IPM). Cracks are beginning to appear in this model which only allows for type (a) changes in the potential. But in a short review we shall have quite enough to do in understanding the progress made in the last decade on the IPM. This paper is divided into three parts. The first deals with how to reduce the IPM to the single electron model (SEM). The second is on a new method where charge transfer is important. The third confronts some standard models with modern calculations
Monte Carlo algorithm for calculating moments of atomic transition arrays
A new method for calculating moments of atomic transition arrays is described. It is based on the collective vector method described in earlier publications. In this new approach a single collective state vector is generated from a single parent state vector encompassing the entire parent basis. The amplitudes of the basis vectors comprising the parent state vector are randomized. Thus it is representative of the entire parent manifold. We show that a statistical estimate of the transition moments of the array is given by a suitable combination of moments of the matrix element of the E1 (electric dipole) operator between the representative parent state vector and the collective state vector derived from it. Although the detailed characteristics of the dispersion in these statistical results are determined by the detailed characteristics of the Hamiltonian and the model space, we find that overall the dispersion decreases inversely as the square root of the dimension of the model space. This is in keeping with results obtained earlier for the Hamiltonian moments of nuclear systems and means, speaking broadly, the larger the problem the more accurate the method
Atomic vapor spectroscopy in integrated photonic structures
Ritter, Ralf; Pernice, Wolfram; Kübler, Harald; Pfau, Tilman; Löw, Robert
2015-01-01
We investigate an integrated optical chip immersed in atomic vapor providing several waveguide geometries for spectroscopy applications. The narrow-band transmission through a silicon nitride waveguide and interferometer is altered when the guided light is coupled to a vapor of rubidium atoms via the evanescent tail of the waveguide mode. We use grating couplers to couple between the waveguide mode and the radiating wave, which allow for addressing arbitrary coupling positions on the chip surface. The evanescent atom-light interaction can be numerically simulated and shows excellent agreement with our experimental data. This work demonstrates a next step towards miniaturization and integration of alkali atom spectroscopy and provides a platform for further fundamental studies of complex waveguide structures.
Atomic vapor spectroscopy in integrated photonic structures
Ritter, Ralf; Kübler, Harald; Pfau, Tilman; Löw, Robert, E-mail: r.loew@physik.uni-stuttgart.de [5. Physikalisches Institut and Center for Integrated Quantum Science and Technology, Universität Stuttgart, Pfaffenwaldring 57, 70550 Stuttgart (Germany); Gruhler, Nico; Pernice, Wolfram [Institute of Nanotechnology, Karlsruhe Institute of Technology, 76344 Eggenstein-Leopoldshafen (Germany)
2015-07-27
We investigate an integrated optical chip immersed in atomic vapor providing several waveguide geometries for spectroscopy applications. The narrow-band transmission through a silicon nitride waveguide and interferometer is altered when the guided light is coupled to a vapor of rubidium atoms via the evanescent tail of the waveguide mode. We use grating couplers to couple between the waveguide mode and the radiating wave, which allow for addressing arbitrary coupling positions on the chip surface. The evanescent atom-light interaction can be numerically simulated and shows excellent agreement with our experimental data. This work demonstrates a next step towards miniaturization and integration of alkali atom spectroscopy and provides a platform for further fundamental studies of complex waveguide structures.
Study of the embedded atom method of atomistic calculations for metals and alloys
Two projects were completed in the past year. The stability of a series of binary alloys was calculated using the embedded-atom method (EAM) with an analytic form for two-body potentials derived previously. Both disordered alloys and intermetallic compounds with the L10 and L12 structures were studied. The calculated heats of solution of alloys of Cu, Ag, Au, Ni, and Pt were satisfactory, while results for alloys containing Pd were too high. Atomistic calculations using the EAM were also carried out for point defects in hcp metals. By comparison with results in the literature, it was found that many body effects from the EAM significantly alter predicted physical properties of hcp metals. For example, the EAM calculations yield anisotropic vacancy diffusion with greater vacancy mobility in the basal plane, and imply that diffusion will start at a lower fraction of the melting temperature
The binding energy of an interstitial hydrogen atom at various lattice sites around the a0/2[111] screw dislocation core in BCC iron has been determined using the first-principles calculation. The calculation was based on the core structure of a screw dislocation with symmetric displacement field, which was obtained using the large-scale supercell containing 231 atoms and 1 x 1 x 4 k-point samplings. The binding or trapping energy of a hydrogen atom at both the t-site (tetrahedral site) and o-site (octahedral site) near a core is found to be approximately 0.2 eV. (author)
Effective oscillator strength distributions are systematically generated and tabulated for the alkali atoms, the alkaline-earth atoms, the alkaline-earth ions, the rare gases and some miscellaneous atoms. These effective distributions are used to compute the dipole, quadrupole and octupole static polarizabilities, and are then applied to the calculation of the dynamic polarizabilities at imaginary frequencies. These polarizabilities can be used to determine the long-range C6, C8 and C10 atom–atom interactions for the dimers formed from any of these atoms and ions, and we present tables covering all of these combinations
Jiang, Jun, E-mail: phyjiang@yeah.net [Key Laboratory of Atomic and Molecular Physics and Functional Materials of Gansu Province, College of Physics and Electronic Engineering, Northwest Normal University, Lanzhou 730070 (China); School of Engineering, Charles Darwin University, Darwin, Northern Territory, 0909 (Australia); Mitroy, J. [School of Engineering, Charles Darwin University, Darwin, Northern Territory, 0909 (Australia); Cheng, Yongjun, E-mail: cyj83mail@gmail.com [School of Engineering, Charles Darwin University, Darwin, Northern Territory, 0909 (Australia); Academy of Fundamental and Interdisciplinary Science, Harbin Institute of Technology, Harbin 150080 (China); Bromley, M.W.J., E-mail: brom@physics.uq.edu.au [School of Mathematics and Physics, The University of Queensland, Brisbane, Queensland 4075 (Australia)
2015-01-15
Effective oscillator strength distributions are systematically generated and tabulated for the alkali atoms, the alkaline-earth atoms, the alkaline-earth ions, the rare gases and some miscellaneous atoms. These effective distributions are used to compute the dipole, quadrupole and octupole static polarizabilities, and are then applied to the calculation of the dynamic polarizabilities at imaginary frequencies. These polarizabilities can be used to determine the long-range C{sub 6}, C{sub 8} and C{sub 10} atom–atom interactions for the dimers formed from any of these atoms and ions, and we present tables covering all of these combinations.
Electronic Structure and Molecular Dynamics Calculations for KBH4
Papaconstantopoulos, Dimitrios; Shabaev, Andrew; Hoang, Khang; Mehl, Michael; Kioussis, Nicholas
2012-02-01
In the search for hydrogen storage materials, alkali borohydrides MBH4 (M=Li, Na, K) are especially interesting because of their light weight and the high number of hydrogen atoms per metal atom. Electronic structure calculations can give insights into the properties of these complex hydrides and provide understanding of the structural properties and of the bonding of hydrogen. We have performed first-principles density-functional theory (DFT) and tight-binding (TB) calculations for KBH4 in both the high temperature (HT) and low temperature (LT) phases to understand its electronic and structural properties. Our DFT calculations were carried out using the VASP code. The results were then used as a database to develop a tight-binding Hamiltonian using the NRL-TB method. This approach allowed for computationally efficient calculations of phonon frequencies and elastic constants using the static module of the NRL-TB, and also using the molecular dynamics module to calculate mean-square displacements and formation energies of hydrogen vacancies.
Atomic Structure and Properties of Extended Defects in Silicon
Buczko, R.; Chisholm, M.F.; Kaplan, T.; Maiti, A.; Mostoller, M.; Pantelides, S.T.; Pennycook, S.J.
1998-10-15
The Z-contrast technique represents a new approach to high-resolution electron microscopy allowing for the first time incoherent imaging of materials on the atomic scale. The key advantages of the technique, an intrinsically higher resolution limit and directly interpretable, compositionally sensitive imaging, allow a new level of insight into the atomic configurations of extended defects in silicon. This experimental technique has been combined with theoretical calculations (a combination of first principles, tight binding, and classical methods) to extend this level of insight by obtaining the energetic and electronic structure of the defects.
Calculation of tolerances in accelerating structures
A method is suggested for calculating tolerances for similar elements of an accelerating-focusing channel of a charged particle linac the particle dynamics in which is described by linear or non-linear equations. Tolerances for each drift tube of the accelerating structure with modified variable-phase focusing are determined with respect to tolerances for the output parameters of an accelerated beam at preset lengths of drift tubes. The tolerances obtained in supposition of equal effects, equal tolerances and those accounting for the cost of fabrication and assembling of the elements of the structure are compared. The algorithm suggested can also be used for calculating tolerances in structures with hard focusing
Geminal embedding scheme for optimal atomic basis set construction in correlated calculations
We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen
Geminal embedding scheme for optimal atomic basis set construction in correlated calculations.
Sorella, S; Devaux, N; Dagrada, M; Mazzola, G; Casula, M
2015-12-28
We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen. PMID:26723656
Geminal embedding scheme for optimal atomic basis set construction in correlated calculations
Sorella, S., E-mail: sorella@sissa.it [International School for Advanced Studies (SISSA), Via Beirut 2-4, 34014 Trieste, Italy and INFM Democritos National Simulation Center, Trieste (Italy); Devaux, N.; Dagrada, M., E-mail: mario.dagrada@impmc.upmc.fr [Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05 (France); Mazzola, G., E-mail: gmazzola@phys.ethz.ch [Theoretische Physik, ETH Zurich, 8093 Zurich (Switzerland); Casula, M., E-mail: michele.casula@impmc.upmc.fr [CNRS and Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05 (France)
2015-12-28
We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.
First principles calculation of current-induced forces in atomic gold contacts
Brandbyge, Mads; Stokbro, Kurt; Taylor, Jeremy; Mozos, Jose-Luis; Ordejon, Pablo
2002-03-01
We have recently developed an first principles method [1] for calculating the electronic structure, electronic transport, and forces acting on the atoms, for atomic scale systems connected to semi-infinite electrodes and with an applied voltage bias. Our method is based on the density functional theory (DFT) as implemented in the well tested SIESTA program [2]. We fully deal with the atomistic structure of the whole system, treating both the contact and the electrodes on the same footing. The effect of the finite bias (including selfconsistency and the solution of the electrostatic problem) is taken into account using nonequilibrium Green's functions. In this talk we show results for the forces acting on the contact atoms due to the nonequilibrium situation in the electronic subsystem, i.e. in the presence of an electronic current. We concentrate on one atom wide gold contacts/wires connected to bulk gold electrodes. References [1] Our implementation is called TranSIESTA and is described in M. Brandbyge, J. Taylor, K. Stokbro, J-L. Mozos, and P. Ordejon, cond-mat/0110650 [2] D. Sanchez-Portal, P. Ordejon, E. Artacho and J. Soler, Int. J. Quantum Chem. 65, 453 (1997).
Real-time feedback from iterative electronic structure calculations
Vaucher, Alain C; Reiher, Markus
2015-01-01
Real-time feedback from iterative electronic structure calculations requires to mediate between the inherently unpredictable execution times of the iterative algorithm employed and the necessity to provide data in fixed and short time intervals for real-time rendering. We introduce the concept of a mediator as a component able to deal with infrequent and unpredictable reference data to generate reliable feedback. In the context of real-time quantum chemistry, the mediator takes the form of a surrogate potential that has the same local shape as the first-principles potential and can be evaluated efficiently to deliver atomic forces as real-time feedback. The surrogate potential is updated continuously by electronic structure calculations and guarantees to provide a reliable response to the operator for any molecular structure. To demonstrate the application of iterative electronic structure methods in real-time reactivity exploration, we implement self-consistent semi-empirical methods as the data source and a...
Embedded-cluster calculations in a numeric atomic orbital density-functional theory framework
We integrate the all-electron electronic structure code FHI-aims into the general ChemShell package for solid-state embedding quantum and molecular mechanical (QM/MM) calculations. A major undertaking in this integration is the implementation of pseudopotential functionality into FHI-aims to describe cations at the QM/MM boundary through effective core potentials and therewith prevent spurious overpolarization of the electronic density. Based on numeric atomic orbital basis sets, FHI-aims offers particularly efficient access to exact exchange and second order perturbation theory, rendering the established QM/MM setup an ideal tool for hybrid and double-hybrid level density functional theory calculations of solid systems. We illustrate this capability by calculating the reduction potential of Fe in the Fe-substituted ZSM-5 zeolitic framework and the reaction energy profile for (photo-)catalytic water oxidation at TiO2(110)
Boucard, Stéphane
1998-01-01
In this work we present the calculations of transition energies inlithiumlike and in exotics atoms : 1) With the use of new ions beamsource, it is now possible to produce highly charged ions. We thushave studied li-like ions hyperfine structure to obtain information onrelativistic many-body corrections and the magnetic sector of strongfield QED. In heavy ions these are of the order of a few percent ofthe total transition energy. We also evaluated the Bohr-Weisskopfcorrection, which depends on...
Atomic and electronic structure of exfoliated black phosphorus
Black phosphorus, a layered two-dimensional crystal with tunable electronic properties and high hole mobility, is quickly emerging as a promising candidate for future electronic and photonic devices. Although theoretical studies using ab initio calculations have tried to predict its atomic and electronic structure, uncertainty in its fundamental properties due to a lack of clear experimental evidence continues to stymie our full understanding and application of this novel material. In this work, aberration-corrected scanning transmission electron microscopy and ab initio calculations are used to study the crystal structure of few-layer black phosphorus. Directly interpretable annular dark-field images provide a three-dimensional atomic-resolution view of this layered material in which its stacking order and all three lattice parameters can be unambiguously identified. In addition, electron energy-loss spectroscopy (EELS) is used to measure the conduction band density of states of black phosphorus, which agrees well with the results of density functional theory calculations performed for the experimentally determined crystal. Furthermore, experimental EELS measurements of interband transitions and surface plasmon excitations are also consistent with simulated results. Finally, the effects of oxidation on both the atomic and electronic structure of black phosphorus are analyzed to explain observed device degradation. The transformation of black phosphorus into amorphous PO3 or H3PO3 during oxidation may ultimately be responsible for the degradation of devices exposed to atmosphere over time
Atomic and electronic structure of exfoliated black phosphorus
Wu, Ryan J.; Topsakal, Mehmet; Jeong, Jong Seok; Wentzcovitch, Renata M.; Mkhoyan, K. Andre, E-mail: mkhoyan@umn.edu [Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455 (United States); Low, Tony; Robbins, Matthew C.; Haratipour, Nazila; Koester, Steven J. [Department of Electrical and Computer Engineering, University of Minnesota, Minneapolis, Minnesota 55455 (United States)
2015-11-15
Black phosphorus, a layered two-dimensional crystal with tunable electronic properties and high hole mobility, is quickly emerging as a promising candidate for future electronic and photonic devices. Although theoretical studies using ab initio calculations have tried to predict its atomic and electronic structure, uncertainty in its fundamental properties due to a lack of clear experimental evidence continues to stymie our full understanding and application of this novel material. In this work, aberration-corrected scanning transmission electron microscopy and ab initio calculations are used to study the crystal structure of few-layer black phosphorus. Directly interpretable annular dark-field images provide a three-dimensional atomic-resolution view of this layered material in which its stacking order and all three lattice parameters can be unambiguously identified. In addition, electron energy-loss spectroscopy (EELS) is used to measure the conduction band density of states of black phosphorus, which agrees well with the results of density functional theory calculations performed for the experimentally determined crystal. Furthermore, experimental EELS measurements of interband transitions and surface plasmon excitations are also consistent with simulated results. Finally, the effects of oxidation on both the atomic and electronic structure of black phosphorus are analyzed to explain observed device degradation. The transformation of black phosphorus into amorphous PO{sub 3} or H{sub 3}PO{sub 3} during oxidation may ultimately be responsible for the degradation of devices exposed to atmosphere over time.
Bannikov, V. V.; Shein, I. R.; Ivanovskii, A. L.
2011-01-01
First-principles FLAPW-GGA band structure calculations were employed to examine the structural, electronic properties and the chemical bonding picture for four ZrCuSiAs-like Th-based quaternary pnictide oxides ThCuPO, ThCuAsO, ThAgPO, and ThAgAsO. These compounds were found to be semimetals and may be viewed as "intermediate" systems between two main isostructural groups of superconducting and semiconducting 1111 phases. The Th 5f states participate actively in the formation of valence bands ...
Elongation method for electronic structure calculations of random DNA sequences.
Orimoto, Yuuichi; Liu, Kai; Aoki, Yuriko
2015-10-30
We applied ab initio order-N elongation (ELG) method to calculate electronic structures of various deoxyribonucleic acid (DNA) models. We aim to test potential application of the method for building a database of DNA electronic structures. The ELG method mimics polymerization reactions on a computer and meets the requirements for linear scaling computational efficiency and high accuracy, even for huge systems. As a benchmark test, we applied the method for calculations of various types of random sequenced A- and B-type DNA models with and without counterions. In each case, the ELG method maintained high accuracy with small errors in energy on the order of 10(-8) hartree/atom compared with conventional calculations. We demonstrate that the ELG method can provide valuable information such as stabilization energies and local densities of states for each DNA sequence. In addition, we discuss the "restarting" feature of the ELG method for constructing a database that exhaustively covers DNA species. PMID:26337429
Hylleraas-Configuration Interaction calculations on the 1 ^1S ground state of helium atom
Ruiz, Maria Belen
2012-01-01
Hylleraas-Configuration Interaction (Hy-CI) calculations on the ground $1^1$S state of helium atom are presented using s-, p-, d-, and f-Slater orbitals of both real and complex form. Techniques of construction of adapted configurations, optimization of the orbital exponents and structure of the wave function expansion are explored. A new method to evaluate the two-electron kinetic energy integrals occurring in the Hy-CI method has been tested here and compared with other methods. The non-rel...
林圣路; 张秋菊; 赵珂; 宋晓红; 张延惠
2002-01-01
By using the region-splitting consistent and iterative method, we calculate the recurrence spectra of lithium atoms in parallel strong external electric and magnetic fields, and obtain the novel resonance structure in the photoabsorption spectrum above the ionization threshold with a constant scaled electric field at F ＝ 0.036, and a scaled energy at e ＝ 0.58 and e = 0.006, respectively. The results are compared with those of hydrogen obtained by using standard closed orbit theory. It is demonstrated that the core-scattered effects exhibited in combination recurrence play a great role.
Atomic Structure of Ultrathin Gold Nanowires.
Yu, Yi; Cui, Fan; Sun, Jianwei; Yang, Peidong
2016-05-11
Understanding of the atomic structure and stability of nanowires (NWs) is critical for their applications in nanotechnology, especially when the diameter of NWs reduces to ultrathin scale (1-2 nm). Here, using aberration-corrected high-resolution transmission electron microscopy (AC-HRTEM), we report a detailed atomic structure study of the ultrathin Au NWs, which are synthesized using a silane-mediated approach. The NWs contain large amounts of generalized stacking fault defects. These defects evolve upon sustained electron exposure, and simultaneously the NWs undergo necking and breaking. Quantitative strain analysis reveals the key role of strain in the breakdown process. Besides, ligand-like morphology is observed at the surface of the NWs, indicating the possibility of using AC-HRTEM for surface ligand imaging. Moreover, the coalescence dynamic of ultrathin Au NWs is demonstrated by in situ observations. This work provides a comprehensive understanding of the structure of ultrathin metal NWs at atomic-scale and could have important implications for their applications. PMID:27071038
Calculation of a shadow cone formed by scattering of an ion beam from an atom
The shadow cone formed by the scattering of ions from an atom was calculated using Moliere potential. Calculated radii of the shadow cone by replacing the Thomas-Fermi radius, a sub(TF), with 0.85 a sub(TF) show a good agreement with an experiment for 1 keV helium ions on a titanium atom. (author)
Calculations of hyperfine interaction strength and life-times of states in highly ionized atoms, using the GRASP atomic structure package, are reported. The calculations aim at providing calibration for Recoil-in-Vacuum nuclear excited state g-factor measurements. The method is outlined and results compared with experiment. Inclusion of decay of higher electronic states is discussed.
Block Tridiagonal Matrices in Electronic Structure Calculations
Petersen, Dan Erik
This thesis focuses on some of the numerical aspects of the treatment of the electronic structure problem, in particular that of determining the ground state electronic density for the non–equilibrium Green’s function formulation of two–probe systems and the calculation of transmission in the Lan...
Cobalamins uncovered by modern electronic structure calculations
Kepp, Kasper Planeta; Ryde, Ulf
electronic-structure calculations, in particular density functional methods, the understanding of the molecular mechanism of cobalamins has changed dramatically, going from a dominating view of trans-steric strain effects to a much more complex view involving an arsenal of catalytic strategies. Among these...
Multilevel domain decomposition for electronic structure calculations
Barrault, M; Hager, W W; Le Bris, C
2005-01-01
We introduce a new multilevel domain decomposition method (MDD) for electronic structure calculations within semi-empirical and Density Functional Theory (DFT) frameworks. This method iterates between local fine solvers and global coarse solvers, in the spirit of domain decomposition methods.
Charge-correlation effects in calculations of atomic short-range order in metallic alloys
Pinski, F.J. [Department of Physics, University of Cincinnati, Cincinnati, Ohio 45221 (United States); Staunton, J.B. [Department of Physics, University of Warwick, Coventry CV4 7AL (United Kingdom); Johnson, D.D. [Department of Materials Science and Engineering, University of Illinois, Urbana, Illinois 61801 (United States)
1998-06-01
The {open_quotes}local{close_quotes} chemical environment that surrounds an atom directly influences its electronic charge density. These atomic charge correlations play an important role in describing the Coulomb and total energies for random substitutional alloys. Although the electronic structure may be well represented by a single-site theory, such as the coherent potential approximation, the electrostatic energy is not as well represented when these charge correlations are ignored. For metals, including the average effect from the charge correlation coming from only the nearest-neighbor shell has been shown to be sufficient to determine accurately the energy of formation. In this paper, we incorporate such charge correlations into the concentration-wave approach for calculating the atomic short-range order in random (substitutional) alloys. We present changes within the formalism, and apply the resulting equations to equiatomic nickel platinum. By including these effects, we obtain significantly better agreement with experimental data. In fact, particular to NiPt, a consequence of the charge correlation is a screening which cancels much of the electrostatic contribution to the energy and thus to the atomic short-range order, resulting in agreement with a picture originally outlined using only {open_quotes}band-energy{close_quotes} contributions. {copyright} {ital 1998} {ital The American Physical Society}
Multigrid Methods in Electronic Structure Calculations
Briggs, E L; Bernholc, J
1996-01-01
We describe a set of techniques for performing large scale ab initio calculations using multigrid accelerations and a real-space grid as a basis. The multigrid methods provide effective convergence acceleration and preconditioning on all length scales, thereby permitting efficient calculations for ill-conditioned systems with long length scales or high energy cut-offs. We discuss specific implementations of multigrid and real-space algorithms for electronic structure calculations, including an efficient multigrid-accelerated solver for Kohn-Sham equations, compact yet accurate discretization schemes for the Kohn-Sham and Poisson equations, optimized pseudo\\-potentials for real-space calculations, efficacious computation of ionic forces, and a complex-wavefunction implementation for arbitrary sampling of the Brillioun zone. A particular strength of a real-space multigrid approach is its ready adaptability to massively parallel computer architectures, and we present an implementation for the Cray-T3D with essen...
Nuclear structure calculations for astrophysical applications
Here we present calculated results on such diverse properties as nuclear energy levels, ground-state masses and shapes, β-decay properties and fission-barrier heights. Our approach to these calculations is to use a unified theoretical framework within which the above properties can all be studied. The results are obtained in the macroscopic-microscopic approach in which a microscopic nuclear-structure single-particle model with extensions is combined with a macroscopic model, such as the liquid drop model. In this model the total potential energy of the nucleus may be calculated as a function of shape. The maxima and minima in this function correspond to such features as the ground state, fission saddle points and shape-isomeric states. Various transition rate matrix elements are determined from wave-functions calculated in the single-particle model with pairing and other relevant residual interactions taken into account
Electronic Structure of Silicon Nanowires Matrix from Ab Initio Calculations.
Monastyrskii, Liubomyr S; Boyko, Yaroslav V; Sokolovskii, Bogdan S; Potashnyk, Vasylyna Ya
2016-12-01
An investigation of the model of porous silicon in the form of periodic set of silicon nanowires has been carried out. The electronic energy structure was studied using a first-principle band method-the method of pseudopotentials (ultrasoft potentials in the basis of plane waves) and linearized mode of the method of combined pseudopotentials. Due to the use of hybrid exchange-correlation potentials (B3LYP), the quantitative agreement of the calculated value of band gap in the bulk material with experimental data is achieved. The obtained results show that passivation of dangling bonds with hydrogen atoms leads to substantial transformation of electronic energy structure. At complete passivation of the dangling silicon bonds by hydrogen atoms, the band gap value takes the magnitude which substantially exceeds that for bulk silicon. The incomplete passivation gives rise to opposite effect when the band gap value decreases down the semimetallic range. PMID:26768147
Dynamic response of some atoms: Many-body calculations
Tančić Aleksandar R.
2005-01-01
Full Text Available The frequency-dependent polarizability in the Hartree-Fock (HF approximation has been corrected for true correlation effects by means of many-body theory. The polarizability has been computed in the Random Phase Approximation with Exchange (RPAE for He, Ar Xe, Kr, Li, Ca through the second (and some higher order in the correlation potential. With this polarizability as input we obtained the values of some atomic interaction constants.
PAMOP: Petascale Atomic, Molecular and Optical Collision Calculations
McLaughlin, Brendan M.; Ballance, Connor P.; Pindzola, Michael S.; Müller, Alfred
2015-01-01
Petaflop architectures are currently being utilized efficiently to perform large scale computations in Atomic, Molecular and Optical Collisions. We solve the Schr\\"odinger or Dirac equation for the appropriate collision problem using the R-matrix or R-matrix with pseudo-states approach. We briefly outline the parallel methodology used and implemented for the current suite of Breit-Pauli and DARC codes. In this report, various examples are shown from our theoretical results compared with exper...
Generalized atomic displacements in crystallographie structure analysis
An attempt is made to summarize both theoretical and experimental aspects of generalized atomic displacement parameters (ADP's) in crystalline matter. Generalized displacement parameters are used to describe the weakening of Bragg intensities via the anharmonic (static or thermal) Debye-Waller factor (DWF) and its real-space counterpart, the generalized atomic probability density function (p.d.f.). The lattice dynamical base of the harmonic and anharmonic thermal DWF is discussed. It is pointed out that the static DWF frequently contains higher-order terms. The mathematical case for an experimental determination of generalized ADP's is given. The most popular current formulations (one-particle potential and statistical approaches) are reviewed and their individual limitations are discussed. Likewise the demands put on the quality and extent of experimental data are assessed. Some aids to the interpretation of generalized ADP's established by crystallographic least-squares procedures are given and a Monte-Carlo method for the calculation of errors in p.d.f. maps is presented. Finally, some prospects for future work are outlined and a more frequent comparison of theoretical calculations and experimental determinations of generalized ADP's is advocated. (orig.)
Calculating atomic and molecular properties using variational Monte Carlo methods
The authors compute a number of properties for the 1 1S, 21S, and 23S states of helium as well as the ground states of H2 and H/+3 using Variational Monte Carlo. These are in good agreement with previous calculations (where available). Electric-response constants for the ground states of helium, H2 and H+3 are computed as derivatives of the total energy. The method used to calculate these quantities is discussed in detail
Large-scale quantum transport calculations for electronic devices with over ten thousand atoms
Lu, Wenchang; Lu, Yan; Xiao, Zhongcan; Hodak, Miro; Briggs, Emil; Bernholc, Jerry
The non-equilibrium Green's function method (NEGF) has been implemented in our massively parallel DFT software, the real space multigrid (RMG) code suite. Our implementation employs multi-level parallelization strategies and fully utilizes both multi-core CPUs and GPU accelerators. Since the cost of the calculations increases dramatically with the number of orbitals, an optimal basis set is crucial for including a large number of atoms in the ``active device'' part of the simulations. In our implementation, the localized orbitals are separately optimized for each principal layer of the device region, in order to obtain an accurate and optimal basis set. As a large example, we calculated the transmission characteristics of a Si nanowire p-n junction. The nanowire is along (110) direction in order to minimize the number dangling bonds that are saturated by H atoms. Its diameter is 3 nm. The length of 24 nm is necessary because of the long-range screening length in Si. Our calculations clearly show the I-V characteristics of a diode, i.e., the current increases exponentially with forward bias and is near zero with backward bias. Other examples will also be presented, including three-terminal transistors and large sensor structures.
Calculated Structural Phase-Transitions in the Alkaline-Earth Metals
Skriver, Hans Lomholt
1982-01-01
The local-density approximation and the linear muffin-tin orbital method have been used within the atomic-sphere approximation to calculate structural energy differences for all the alkaline earth metals at zero temperature. At ordinary pressure the calculations predict the crystal structure sequ...... sequence hcp→fcc→bcc as a function of atomic number. As a function of pressure they predict the structure sequence fcc→bcc→hcp. The structural transitions and the onset of superconductivity under pressure are correlated with the d occupation number.......The local-density approximation and the linear muffin-tin orbital method have been used within the atomic-sphere approximation to calculate structural energy differences for all the alkaline earth metals at zero temperature. At ordinary pressure the calculations predict the crystal structure...
Fock space relativistic coupled-Cluster calculations of Two-Valence Atoms
Mani, B. K.; Angom, D.
2010-01-01
We have developed an all particle Fock-space relativistic coupled-cluster method for two-valence atomic systems. We then describe a scheme to employ the coupled-cluster wave function to calculate atomic properties. Based on these developments we calculate the excitation energies, magnetic hyperfine constants and electric dipole matrix elements of Sr, Ba and Yb. Further more, we calculate the electric quadrupole HFS constants and the electric dipole matrix elements of Sr$^+$, Ba$^+$ and Yb$^+$...
Calculation of the elastic collision properties of Na and Li atoms at ultracold temperature
Zhang Ji-Cai; Zhang Ying; Du Bing-Ge; Sun Jin-Feng
2009-01-01
This paper firstly reports a theoretical study of elastic scattering properties in a mixture of 23Na and 7Li atoms at cold and ultracold temperatures in detail.Based on the new constructed accurate singlet X1∑+g and the triplet a3∑+u states interatomic potentials for 23Na7Li mixture,it calculates the scattering lengths and the effective ranges by three computational methods,and obtains good agreements.Using the mass scaling method,it also calculates 23Na6Li scattering lengths and s-wave and total elastic cross sections,whose rich resonance structures were found and interpreted in terms of quasibound diatomic levels trapped behind a centrifugal barrier.
PAMOP: Petascale Atomic, Molecular and Optical Collision Calculations
McLaughlin, Brendan M; Pindzola, Michael S; Müller, Alfred
2015-01-01
Petaflop architectures are currently being utilized efficiently to perform large scale computations in Atomic, Molecular and Optical Collisions. We solve the Schr\\"odinger or Dirac equation for the appropriate collision problem using the R-matrix or R-matrix with pseudo-states approach. We briefly outline the parallel methodology used and implemented for the current suite of Breit-Pauli and DARC codes. In this report, various examples are shown from our theoretical results compared with experimental results obtained from Synchrotron Radiation facilities where the Cray architecture at HLRS is playing an integral part in our computational projects.
The variational method in the atomic structure calcularion
The importance and limitations of variational methods on the atomic structure calculations is set into relevance. Comparisons are made to the Perturbation Theory. Ilustrating it, the method is applied to the H-, H+ and H+2 simple atomic structure systems, and the results are analysed with basis on the study of the associated essential eigenvalue spectrum. Hydrogenic functions (where the screening constants are replaced by variational parameters) are combined to construct the wave function with proper symmetry for each one of the systems. This shows the existence of a bound state for H-, but no conclusions can be made for the others, where it may or may not be necessary to use more flexible wave functions, i.e., with greater number of terms and parameters. (author)
Kodera, M.; Higuchi, K.; Narita, A; Higuchi, M
2008-01-01
We present a density-functional scheme for calculating the orbital-dependent exchange-correlation potential using the virial theorem as a sum rule. In order to check the validity of this scheme, atomic-structure calculations only with the exchange potential are performed. The accuracy of this scheme is shown to be comparable to that of the optimized effective potential (OEP) method, while the computational workload is extremely reduced compared to the OEP method.
Lattice location of dopant atoms: An -body model calculation
N K Deepak
2010-03-01
The channelling and scattering yields of 1 MeV -particles in the $\\langle 1 0 0 \\rangle$, $\\langle 1 1 0 \\rangle and $\\langle 1 1 1 \\rangle$ directions of silicon implanted with bismuth and ytterbium have been simulated using -body model. The close encounter yield from dopant atoms in silicon is determined from the flux density, using the Bontemps and Fontenille method. All previous works reported in literature so far have been done with computer programmes using a statistical analytical expression or by a binary collision model or a continuum model. These results at the best gave only the transverse displacement of the lattice site from the concerned channelling direction. Here we applied the superior -body model to study the yield from bismuth in silicon. The finding that bismuth atom occupies a position close to the silicon substitutional site is new. The transverse displacement of the suggested lattice site from the channelling direction is consistent with the experimental results. The above model is also applied to determine the location of ytterbium in silicon. The present values show good agreement with the experimental results.
Calculated Structural Phase-Transitions in the Alkaline-Earth Metals
Skriver, Hans Lomholt
1982-01-01
The local-density approximation and the linear muffin-tin orbital method have been used within the atomic-sphere approximation to calculate structural energy differences for all the alkaline earth metals at zero temperature. At ordinary pressure the calculations predict the crystal structure...
Extracting chemical information from plane wave calculations by a 3D 'fuzzy atoms' analysis
Bakó, I.; Stirling, A.; Seitsonen, A. P.; Mayer, I.
2013-03-01
Bond order and valence indices have been calculated by the method of the three-dimensional 'fuzzy atoms' analysis, using the numerical molecular orbitals obtained from plane wave DFT calculations, i.e., without introducing any external atom-centered functions. Weight functions of both Hirshfeld and Becke types have been applied. The results are rather close to the similar 'fuzzy atoms' ones obtained by using atom-centered basis sets and agree well with the chemical expectations, stressing the power of the genuine chemical concepts.
Atomic Structure of Benzene Which Accounts for Resonance Energy
Heyrovska, Raji
2008-01-01
Benzene is a hexagonal molecule of six carbon atoms, each of which is bound to six hydrogen atoms. The equality of all six CC bond lengths, despite the alternating double and single bonds, and the surplus (resonance) energy, led to the suggestion of two resonanting structures. Here, the new atomic structure shows that the bond length equality is due to three carbon atoms with double bond radii bound to three other carbon atoms with resonance bond radii (as in graphene). Consequently, there ar...
Calculation of electron scattering on atoms and ions
Bray, I.
1995-02-01
This paper reviews the applications of the convergent close-coupling (CCC) method to electron scattering on light atoms and ions. Particular emphasis is given to those areas where other theories have difficulty, e g. total ionization cross sections and the associated spin asymmetries. It begins with the simplest application to the Temkin-Poet model problem of electron-hydrogen scattering, which is used to validate the CCC approach. Subsequently, results are given for electron impact ionization of various initial states of the targets H(1s,2s), He(1{sup 1}S,2{sup 3.1}S), He{sup +}(1s), Li(2s), O{sup 5+}(2s) and Na(3s). 50 refs., 10 figs.
Calculation of electron scattering on atoms and ions
This paper reviews the applications of the convergent close-coupling (CCC) method to electron scattering on light atoms and ions. Particular emphasis is given to those areas where other theories have difficulty, e g. total ionization cross sections and the associated spin asymmetries. It begins with the simplest application to the Temkin-Poet model problem of electron-hydrogen scattering, which is used to validate the CCC approach. Subsequently, results are given for electron impact ionization of various initial states of the targets H(1s,2s), He(11S,23.1S), He+(1s), Li(2s), O5+(2s) and Na(3s). 50 refs., 10 figs
Learning from numerical calculations of ion-atom collisions
Violent collision of two independent many-particle systems, victims, are discussed in the atomic sphere. The asymmetric region where the charge of the projectile Z/sub p/ is less than the target nuclear charge Z/sub n/ is now well understood though interesting details still need to be worked out. Negatively charged projectiles offer a new illustration of Fadeev re-arrangement collisions. Multi-electron coherence effects illustrate the richness of the field but a symmetric (Z/sub p/ approx. Z/sub n/) collision treatment is needed. A new one and a half center expansion method promises a solution to this problem. Future areas of interest are discussed
Prediction of hexagonal Ta2O5 structure by first-principles calculations
The undetermined crystal structure of hexagonal Ta2O5 is predicted by the first-principles ultrasoft pseudopotential calculations. Among the several hypothetical structures, the structure which is energetically favorable and consistent with the available experimental data is determined by minimizing the total energy with respect to the volume, the c/a ratio, and the atomic configurations. The predicted structure has the space group of P6/mmm with two formula units in the unit cell. The calculated lattice constants, a=7.191 Angstrom and c=3.831 Angstrom, are in good agreement with the x-ray-diffraction measurement if a is assumed to be twice as large as the measured value due to the extinction rule. The coordination number of O atom to Ta atom is 8 for one Ta atom and 6 for the other three Ta atoms. copyright 1997 The American Physical Society
Yingfeng Li
2013-05-01
Full Text Available The ways in which carbon atoms coalesce over the steps on copper (111 surface are ascertained by density functional theory (DFT calculations in the context of chemical vapor deposition (CVD growth of graphene. Two strategies, (1 by putting carbon atoms on and under the steps separately and (2 by importing additional carbon atoms between the ones separated by the steps, have been attempted to investigate if an over-step coalescence of carbon atoms could take place. Based on analyses about the optimized configurations and adsorption energies of carbon atoms nearby the steps, as well as the energy evolution curve of the system throughout the geometry optimizations process, we determined the main way in which graphene grows over the steps continuously: the carbon atoms, adsorbed additionally on the locations between the already existing ones which are separated by the steps, link them (these carbon atoms separated by the steps together. The direct over-step coalescence of the carbon atoms separated by the steps is very difficult, although the energy barrier preventing their coalescence can be weakened by importing carbon atoms on and under the steps gradually. Our results imply potential applications in directing the fabrication of graphene with particular structure by controlling the surface topography of copper substrate.
Aggarwal, Kanti M.
2015-11-01
Recently, Goyal et al. [1] reported energies and lifetimes (τ) for the lowest 113 levels of the 2s22p5, 2s2p6, 2s22p43ℓ, 2s2p53ℓ and 2p63ℓ configurations of F-like Sr XXX. For the calculations they adopted the multi-configuration Dirac-Fock (MCDF) and the flexible atomic code (FAC). Additionally, they also listed radiative rates (A- values), oscillator strengths (f- values) and line strengths (S- values) for four types of transitions, namely electric dipole (E1), electric quadrupole (E2), magnetic dipole (M1) and magnetic quadrupole (M2), but only from the ground to the higher excited levels. However, there are two clear anomalies in their reported data. Firstly, the f-values listed from FAC in their Tables 3-6 are larger than from MCDF by a factor of two, for all transitions. This is because they have blindly listed the output from FAC without realising that, unlike MCDF, FAC lists ωf where ω is the statistical weight, and happens to be exactly 2 in the present case. Secondly, their lifetime for level 2 (2s22p51/2 o 2P) is incorrect. This is because the dominant contributing transition for this level is 1-2 M1 for which A=3.25×106 s-1, listed (correctly) in their Table 5, and this leads to τ=3.08×10-7 s, and not 1.54×10-7 s, as listed in their Table 1.
Atomic density functions: atomic physics calculations analyzed with methods from quantum chemistry
Borgoo, Alex; Geerlings, P
2011-01-01
This contribution reviews a selection of findings on atomic density functions and discusses ways for reading chemical information from them. First an expression for the density function for atoms in the multi-configuration Hartree--Fock scheme is established. The spherical harmonic content of the density function and ways to restore the spherical symmetry in a general open-shell case are treated. The evaluation of the density function is illustrated in a few examples. In the second part of the paper, atomic density functions are analyzed using quantum similarity measures. The comparison of atomic density functions is shown to be useful to obtain physical and chemical information. Finally, concepts from information theory are introduced and adopted for the comparison of density functions. In particular, based on the Kullback--Leibler form, a functional is constructed that reveals the periodicity in Mendeleev's table. Finally a quantum similarity measure is constructed, based on the integrand of the Kullback--L...
Effects of NMR spectral resolution on protein structure calculation.
Suhas Tikole
Full Text Available Adequate digital resolution and signal sensitivity are two critical factors for protein structure determinations by solution NMR spectroscopy. The prime objective for obtaining high digital resolution is to resolve peak overlap, especially in NOESY spectra with thousands of signals where the signal analysis needs to be performed on a large scale. Achieving maximum digital resolution is usually limited by the practically available measurement time. We developed a method utilizing non-uniform sampling for balancing digital resolution and signal sensitivity, and performed a large-scale analysis of the effect of the digital resolution on the accuracy of the resulting protein structures. Structure calculations were performed as a function of digital resolution for about 400 proteins with molecular sizes ranging between 5 and 33 kDa. The structural accuracy was assessed by atomic coordinate RMSD values from the reference structures of the proteins. In addition, we monitored also the number of assigned NOESY cross peaks, the average signal sensitivity, and the chemical shift spectral overlap. We show that high resolution is equally important for proteins of every molecular size. The chemical shift spectral overlap depends strongly on the corresponding spectral digital resolution. Thus, knowing the extent of overlap can be a predictor of the resulting structural accuracy. Our results show that for every molecular size a minimal digital resolution, corresponding to the natural linewidth, needs to be achieved for obtaining the highest accuracy possible for the given protein size using state-of-the-art automated NOESY assignment and structure calculation methods.
Protein structure determination using Nuclear Magnetic Resonance (NMR) requires the use of molecular dynamics programs that incorporate both NMR experimental and implicit atomic data. Atomic parameters for each amino acid type are encoded in libraries used by structure calculation programs such as DYANA and AMBER. However, only a few non-standard amino acid library sets are included in these programs or the molecular visualization program MOLMOL. Our laboratory is calculating the phosphorylated and non-phosphorylated states of peptides and proteins using NMR methods. To calculate chemically correct structures, we have extended the available molecular libraries for these programs to include the modified amino acids phosphoserine, phosphothreonine, and phosphotyrosine
Ballance, Connor
2013-05-01
Over the last couple of decades, a number of advanced non-perturbative approaches such as the R-matrix, TDCC and CCC methods have made great strides in terms of improved target representation and investigating fundamental 2-4 electron problems. However, for the electron-impact excitation of near-neutral species or complicated open-shell atomic systems we are forced to make certain compromises in terms of the atomic structure and/or the number of channels included in close-coupling expansion of the subsequent scattering calculation. The availability of modern supercomputing architectures with hundreds of thousands of cores, and the emergence new opportunities through GPU usauge offers one possibility to address some of these issues. To effectively harness this computational power will require significant revision of the existing code structures. I shall discuss some effective strategies within a non-relativistic and relativistic R-matrix framework using the examples detailed below. The goal is to extend existing R-matrix methods from 1-2 thousand close coupled channels to 10,000 channels. With the construction of the ITER experiment in Cadarache, which will have Tungsten plasma-facing components, there is an urgent diagnostic need for the collisional rates for the near-neutral ion stages. In particular, spectroscopic diagnostics of impurity influx require accurate electron-impact excitation and ionisation as well as a good target representation. There have been only a few non-perturbative collisional calculations for this system, and the open-f shell ion stages provide a daunting challenge even for perturbative approaches. I shall present non-perturbative results for for the excitation and ionisation of W3+ and illustrate how these fundamental calculations can be integrated into a meaningful diagnostic for the ITER device. We acknowledge support from DoE fusion.
Analytical calculations of scattering lengths in atomic physics
We describe a method for evaluating analytical long-range contributions to scattering lengths for some potentials used in atomic physics. We assume that an interaction potential between colliding particles consists of two parts. The form of a short-range component, vanishing beyond some distance from the origin (a core radius), need not be given. Instead, we assume that a set of short-range scattering lengths due to that part of the interaction is known. A long-range tail of the potential is chosen to be an inverse power potential, a superposition of two inverse power potentials with suitably chosen exponents or the Lent potential. For these three classes of long-range interactions a radial Schrodinger equation at zero energy may be solved analytically with solutions expressed in terms of the Bessel, Whittaker and Legendre functions, respectively. We utilize this fact and derive exact analytical formulae for the scattering lengths. The expressions depend on the short-range scattering lengths, the core radius and parameters characterizing the long-range part of the interaction. Cases when the long-range potential (or its part) may be treated as a perturbation are also discussed and formulae for scattering lengths linear in strengths of the perturbing potentials are given. It is shown that for some combination of the orbital angular momentum quantum number and an exponent of the leading term of the potential the derived formulae, exact or approximate, take very simple forms and contain only polynomial and trigonometric functions. The expressions obtained in this paper are applicable to scattering of charged particles by neutral targets and to collisions between neutrals. The results are illustrated by accelerating convergence of scattering lengths computed for e--Xe and Cs-Cs systems. (author)
Completed by recent contributions on various topics (atoms and the Brownian motion, the career of Jean Perrin, the evolution of atomic physics since Jean Perrin, relationship between scientific atomism and philosophical atomism), this book is a reprint of a book published at the beginning of the twentieth century in which the author addressed the relationship between atomic theory and chemistry (molecules, atoms, the Avogadro hypothesis, molecule structures, solutes, upper limits of molecular quantities), molecular agitation (molecule velocity, molecule rotation or vibration, molecular free range), the Brownian motion and emulsions (history and general features, statistical equilibrium of emulsions), the laws of the Brownian motion (Einstein's theory, experimental control), fluctuations (the theory of Smoluchowski), light and quanta (black body, extension of quantum theory), the electricity atom, the atom genesis and destruction (transmutations, atom counting)
Hydrogen atom position in hydrated iodide anion from x-ray absorption near edge structure
Hydrogen atom position in the hydrated iodide anion complex is determined from X-ray Absorption Near Edge Structure (XANES) of an aqueous iodide solution at both the K- and L-edges. The spectra are compared with the theoretical ones calculated by using the FEFF method for several model geometries having hydrogen atoms at different positions. Satisfactory agreements are obtained from models with an almost linear alignment of iodine-hydrogen-oxygen atoms, indicating the capability of the XANES analysis when it is combined with the multiple scattering calculations as a method to detect the positions of hydrogen atoms in the first coordination sphere. (author)
Relativistic ab initio calculations for ion-atom collisions
Within the independent particle model we solve the time---dependent single-particle equation using ab initio SCF-DIRAC-FOCK-SLATER wavefunctions as a basis. To reinstate the many-particle aspect of the collision system we use the inclusive probability formalism to answer experimental questions. As an example we show an application to the case of S15+ on Ar where experimental data on the K-K charge transfer are available for a wide range of impact energies from 4.7 to 90 MeV. Our molecular adiabatic calculations and the evaluation using the inclusive probability formalism show good results in the low energy range from 4.7 to 16 MeV impact energy
Isogeometric analysis in electronic structure calculations
Cimrman, R.; Novák, M.; Kolman, Radek; Tůma, Miroslav; Vackář, Jiří
Ostrava: Ústav geoniky AV ČR, 2014 - (Blaheta, R.; Starý, J.; Sysalová, D.). s. 49-49 ISBN 978-80-86407-47-0. [Modelling 2014. 02.06.2014-06.06.2014, Rožnov pod Radhoštěm] R&D Projects: GA ČR(CZ) GAP101/12/2315; GA ČR(CZ) GAP108/11/0853 Institutional support: RVO:61388998 ; RVO:68378271 ; RVO:67985807 Keywords : isogeometric analysis * electronic structure calculations * density functional theory Subject RIV: JJ - Other Materials
Study on structures and electronic properties of NaxV (x = 1-12) atomic clusters
Superatoms, novel entities being studied extensively in recent years, can be stabilized by mixing with transition metal atoms. The aim of this paper is to present some recent theoretical results on the application of quantum calculations for examining the atomic clusters NaxV (x = 1-12) made from the mixing of Nax superatoms with vanadium transition metal atom. Optimized structures of NaxV, NaxV+ and NaxV- are determined by using the TPSSTPSS/DZVP DFT calculations. Characteristics of optimized structures, as point group symmetry, chemical hardness (η), absolute electronegativity (χ), electrophilicity index (ω), fragmentation energy (Ef), secondary energy (∆2E), are calculated. The obtained results point out that among different structures of an atomic cluster, the more negative total energy the more stable structure and the Na8V cluster is the most stable in NaxV (x = 1-12) clusters. (author)
Evolution of atomic structure during nanoparticle formation
Christoffer Tyrsted
2014-05-01
Full Text Available Understanding the mechanism of nanoparticle formation during synthesis is a key prerequisite for the rational design and engineering of desirable materials properties, yet remains elusive due to the difficulty of studying structures at the nanoscale under real conditions. Here, the first comprehensive structural description of the formation of a nanoparticle, yttria-stabilized zirconia (YSZ, all the way from its ionic constituents in solution to the final crystal, is presented. The transformation is a complicated multi-step sequence of atomic reorganizations as the material follows the reaction pathway towards the equilibrium product. Prior to nanoparticle nucleation, reagents reorganize into polymeric species whose structure is incompatible with the final product. Instead of direct nucleation of clusters into the final product lattice, a highly disordered intermediate precipitate forms with a local bonding environment similar to the product yet lacking the correct topology. During maturation, bond reforming occurs by nucleation and growth of distinct domains within the amorphous intermediary. The present study moves beyond kinetic modeling by providing detailed real-time structural insight, and it is demonstrated that YSZ nanoparticle formation and growth is a more complex chemical process than accounted for in conventional models. This level of mechanistic understanding of the nanoparticle formation is the first step towards more rational control over nanoparticle synthesis through control of both solution precursors and reaction intermediaries.
Absence of a stable atomic structure in fluorinated graphene.
Boukhvalov, Danil W
2016-05-21
Based on the results of first-principles calculations we demonstrate that significant distortion of graphene sheets caused by adsorption of fluorine atoms leads to the formation of metastable patterns for which the next step of fluorination is considerably less energetically favorable. Existence of these stable patterns oriented along the armchair direction makes possible the synthesis of various CFx structures. The combination of strong distortion of the nonfluorinated graphene sheet with the doping caused by the polar nature of C-F bonds reduces the energy cost of migration and the energy of migration barriers, making possible the migration of fluorine atoms on the graphene surface as well as transformation of the shapes of fluorinated areas. The decreasing energy cost of migration with increasing fluorine content also leads to increasing numbers of single fluorine adatoms, which could be the source of magnetic moments. PMID:27116897
Electron Affinity Calculations for Atoms: Sensitive Probe of Many-Body Effects
Felfli, Z.; Msezane, A. Z.
2016-05-01
Electron-electron correlations and core-polarization interactions are crucial for the existence and stability of most negative ions. Therefore, they can be used as a sensitive probe of many-body effects in the calculation of the electron affinities (EAs) of atoms. The importance of relativistic effects in the calculation of the EAs of atoms has recently been assessed to be insignificant up to Z of 85. Here we use the complex angular momentum (CAM) methodology wherein is embedded fully the electron-electron correlations, to investigate core-polarization interactions in low-energy electron elastic scattering from the atoms In, Sn, Eu, Au and At through the calculation of their EAs. For the core-polarization interaction we use the rational function approximation of the Thomas-Fermi potential, which can be analytically continued into the complex plane. The EAs are extracted from the large resonance peaks in the CAM calculated low-energy electron-atom scattering total cross sections and compared with those from measurements and sophisticated theoretical methods. It is concluded that when the electron-electron correlations and core polarization interactions (both major many-body effects) are accounted for adequately the importance of relativity on the calculation of the EAs of atoms can be assessed. Even for the high Z (85) At atom relativistic effects are estimated to contribute a maximum of 3.6% to its EA calculation.
Unfolding method for first-principles LCAO electronic structure calculations
Lee, Chi-Cheng; Yamada-Takamura, Yukiko; Ozaki, Taisuke
2013-08-01
Unfolding the band structure of a supercell to a normal cell enables us to investigate how symmetry breakers such as surfaces and impurities perturb the band structure of the normal cell. We generalize the unfolding method, originally developed based on Wannier functions, to the linear combination of atomic orbitals (LCAO) method, and present a general formula to calculate the unfolded spectral weight. The LCAO basis set is ideal for the unfolding method because the basis functions allocated to each atomic species are invariant regardless of the existence of surface and impurity. The unfolded spectral weight is well defined by the property of the LCAO basis functions. In exchange for the property, the non-orthogonality of the LCAO basis functions has to be taken into account. We show how the non-orthogonality can be properly incorporated in the general formula. As an illustration of the method, we calculate the dispersive quantized spectral weight of a ZrB2 slab and show strong spectral broadening in the out-of-plane direction, demonstrating the usefulness of the unfolding method.
Unfolding method for first-principles LCAO electronic structure calculations
Unfolding the band structure of a supercell to a normal cell enables us to investigate how symmetry breakers such as surfaces and impurities perturb the band structure of the normal cell. We generalize the unfolding method, originally developed based on Wannier functions, to the linear combination of atomic orbitals (LCAO) method, and present a general formula to calculate the unfolded spectral weight. The LCAO basis set is ideal for the unfolding method because the basis functions allocated to each atomic species are invariant regardless of the existence of surface and impurity. The unfolded spectral weight is well defined by the property of the LCAO basis functions. In exchange for the property, the non-orthogonality of the LCAO basis functions has to be taken into account. We show how the non-orthogonality can be properly incorporated in the general formula. As an illustration of the method, we calculate the dispersive quantized spectral weight of a ZrB2 slab and show strong spectral broadening in the out-of-plane direction, demonstrating the usefulness of the unfolding method. (paper)
Structure and properties of atomic nanoclusters
Alonso, Julio A
2011-01-01
Atomic clusters are aggregates of atoms containing a few to several thousand atoms. Due to the small size of these pieces of matter, the properties of atomic clusters in general are different from those of the corresponding material in the macroscopic bulk phase. This monograph presents the main developments of atomic clusters and the current status of the field. The book treats different types of clusters with very different properties: clusters in which the atoms or molecules are tied by weak van der Waals interactions, metallic clusters, clusters of ionic materials, and network clusters mad
Radziute, Laima; Jonsson, Per; Biero, Jacek
2013-01-01
The multiconfiguration Dirac-Hartree-Fock (MCDHF) method has been employed to calculate atomic electric dipole moments (EDM) of 225^Ra, 199^Hg, and 171^Yb. For the calculations of the matrix elements we extended the relativistic atomic structure package GRASP2K. The extension includes programs to evaluate matrix elements of (P, T)-odd e-N tensor-pseudotensor and pseudoscalar-scalar interactions, the atomic electric dipole interaction, the nuclear Schiff moment, and the interaction of the electron electric dipole moment with nuclear magnetic moments. The interelectronic interactions were accounted for through valence and core-valence electron correlation effects. The electron shell relaxation was included with separately optimised wave functions of opposite parities.
Structurally uniform and atomically precise carbon nanostructures
Segawa, Yasutomo; Ito, Hideto; Itami, Kenichiro
2016-01-01
Nanometre-sized carbon materials consisting of benzene units oriented in unique geometric patterns, hereafter named nanocarbons, conduct electricity, absorb and emit light, and exhibit interesting magnetic properties. Spherical fullerene C60, cylindrical carbon nanotubes and sheet-like graphene are representative forms of nanocarbons, and theoretical simulations have predicted several exotic 3D nanocarbon structures. At present, synthetic routes to nanocarbons mainly lead to mixtures of molecules with a range of different structures and properties, which cannot be easily separated or refined into pure forms. Some researchers believe that it is impossible to synthesize these materials in a precise manner. Obtaining ‘pure’ nanocarbons is a great challenge in the field of nanocarbon science, and the construction of structurally uniform nanocarbons, ideally as single molecules, is crucial for the development of functional materials in nanotechnology, electronics, optics and biomedical applications. This Review highlights the organic chemistry approach — more specifically, bottom-up construction with atomic precision — that is currently the most promising strategy towards this end.
Inner Space: The Structure of the Atom
Glasstone, Samuel
1972-01-01
The atom is now regarded as the smallest possible particle of an element that retains the identity of that element. The atoms of an element determine the characteristics of that particular element. One of the purposes of this booklet is to explain how the atoms of various elements differ from one another.
Relativistic calculations of quasi-one-electron atoms and ions using Laguerre and Slater spinors
Jiang, Jun; Cheng, Yongjun; Bromley, Michael W J
2016-01-01
A relativistic description of the structure of heavy alkali atoms and alkali-like ions using S-spinors and L-spinors has been developed. The core wavefunction is defined by a Dirac-Fock calculation using an S-spinors basis. The S-spinor basis is then supplemented by a large set of L-spinors for the calculation of the valence wavefunction in a frozen-core model. The numerical stability of the L-spinor approach is demonstrated by computing the energies and decay rates of several low-lying hydrogen eigenstates, along with the polarizabilities of a $Z=60$ hydrogenic ion. The approach is then applied to calculate the dynamic polarizabilities of the $5s$, $4d$ and $5p$ states of Sr$^+$. The magic wavelengths at which the Stark shifts between different pairs of transitions are zero are computed. Determination of the magic wavelengths for the $5s \\to 4d_{\\frac32}$ and $5s \\to 4d_{\\frac52}$ transitions near $417$~nm (near the wavelength for the $5s \\to 5p_j$ transitions) would allow a determination of the oscillator s...
Atomic and electronic structure of divacancies in carbon nanotubes
Berber, Savas; Oshiyama, Atsushi
2008-04-01
We present atomic and electronic structure of divacancies in carbon nanotubes, which is calculated using the density functional theory. Divacancies in carbon nanotubes self-heal by spontaneous reconstructions, which consist of concerted bond formations. Divacancy formation energies EDV , which strongly depend on the divacancy orientation with respect to the tube axis, are in the range of 2.8 4.3 eV for favorable orientations in the nanotubes of 4 9Å diameter, making divacancies more probable than monovacancies in carbon nanotubes. Defect related states lead to a higher density of states around the Fermi level. Semiconducting nanotubes develop midgap levels that may adversely affect the functionality of carbon nanotube based devices. Our spin polarized density functional calculations show that the exchange splitting of defect-related bands in nonsemiconducting defective nanotubes leads to net spin polarizations of ρ↑-ρ↓≤0.5μB per divacancy for some divacancy orientations.
Kevorkyants, Ruslan; Close, David M; Pavanello, Michele
2013-01-01
We present an application of the linear scaling Frozen Density Embedding (FDE) formulation of subsystem DFT to the calculation of isotropic hyperfine coupling constants (hfccs) of atoms belonging to a guanine radical cation embedded in a guanine hydrochloride monohydrate crystal. The model systems considered range from an isolated guanine to a 15,000 atom QM/MM cluster where the QM region is comprised of 36 protonated guanine cations, 36 chlorine anions and 42 water molecules. Our calculations show that the embedding effects of the surrounding crystal cannot be reproduced neither by small model systems nor by a pure QM/MM procedure. Instead, a large QM region is needed to fully capture the complicated nature of the embedding effects in this system. The unprecedented system size for a relativistic all-electron isotropic hfccs calculation can be approached in this work because the local nature of the electronic structure of the organic crystals considered is fully captured by the FDE approach.
Parquet theory in nuclear structure calculations
The thesis concerns a numerical implementation of the Parquet summation of diagrams within Green's functions theory applied to calculations of nuclear systems. The main motivation has been to investigate whether it is possible to develop this approach to a level comparable in accuracy and reliability to other ab initio nuclear structure methods. The Green's functions approach is theoretically well-established in many-body theory, but to our knowledge, no actual application to nuclear systems has been previously published. It has a number of desirable properties, foremost the gently scaling with system size compared to direct diagonalization and the closeness to experimentally accessible quantities. The main drawback is the numerical instabilities due to the pole structure of the one-particle propagator, leading to convergence difficulties. This issue is one of the main focal points of the work presented in this thesis, and strategies to improve the convergence properties are described and investigated. We have applied the method both to a simple model which can be solved by exact diagonalization and to the more realistic 4He system. The results shows that our implementation is close to the exact solution in the simple model as long as the interaction strengths are small. As the number of particles increases, convergence is increasingly hard to obtain. In the 4He case, we obtain results in the vicinity of the results from comparable approaches. The numerical in-stabilities in the current implementation still prevents the desired accuracy and stability necessary to achieve the current benchmark standards. (Author)
The calculation of satellite line structures in highly stripped plasmas
This is the final report of a three-year, Laboratory Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). Recently developed high-resolution x-ray spectrographs have made it possible to measure satellite structures from various plasma sources with great detail. These lines are weak optically thin lines caused by the decay of dielectronic states and generally accompany the resonance lines of H-like and He-like ions. The Los Alamos atomic physics and kinetics codes provide a unique capability for calculating the position and intensities of such lines. These programs have been used to interpret such highly resolved spectral measurements from pulsed power devices and laser produced plasmas. Some of these experiments were performed at the LANL Bright Source and Trident laser facilities. The satellite structures are compared with calculations to diagnose temperatures and densities. The effect of non-thermal electron distributions of electrons on calculated spectra was also considered. Collaborations with Russian scientists have added tremendous value to this research die to their vast experience in x-ray spectroscopy
The calculation of satellite line structures in highly stripped plasmas
Abdallah, J. Jr.; Kilcrease, D.P. [Los Alamos National Lab., NM (United States); Faenov, A.Ya.; Pikuz, T.A. [Multicharged Ion Spectra Data Center, Moscow (Russian Federation)
1998-11-01
This is the final report of a three-year, Laboratory Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). Recently developed high-resolution x-ray spectrographs have made it possible to measure satellite structures from various plasma sources with great detail. These lines are weak optically thin lines caused by the decay of dielectronic states and generally accompany the resonance lines of H-like and He-like ions. The Los Alamos atomic physics and kinetics codes provide a unique capability for calculating the position and intensities of such lines. These programs have been used to interpret such highly resolved spectral measurements from pulsed power devices and laser produced plasmas. Some of these experiments were performed at the LANL Bright Source and Trident laser facilities. The satellite structures are compared with calculations to diagnose temperatures and densities. The effect of non-thermal electron distributions of electrons on calculated spectra was also considered. Collaborations with Russian scientists have added tremendous value to this research die to their vast experience in x-ray spectroscopy.
Structural Features of Boron-Doped Si(113) Surfaces Simulated by ab initio Calculations
LIAO Long-Zhong; LIU Zheng-Hui; ZHANG Zhao-Hui
2008-01-01
Based on ab initio calculations, boron-doped Si(113) surfaces have been simulated and atomic structures of the surfaces have been proposed. It has been determined that surface features of empty and filled states that are separately localized at pentamers and adatoms indicates a low surface density of B atoms, while it is attributed to heavy doping of B atoms at the second layer that pentamers and adatoms are both present in an image of scanning tunnelling microscopy. B doping at the second layer should be balanced by adsorbed B or Si atoms beside the adatoms and inserted B interstitials below the adatoms.
Study of the atomic rare gas behaviour by ab initio calculations
The atomic behaviour of helium and krypton in uranium dioxide has been studied using an ab initio simulation technique. Incorporation energies and solution energies of these two rare gases have been calculated. Krypton atoms are found to be insoluble in this nuclear fuel whatever the trap considered and their presence in the lattice induces swelling when they are located in interstitial sites or in oxygen vacancies. Due to its small atomic size, the predicted helium behaviour is very different. Indeed, helium is found to be soluble in stoichiometric and hyper-stoichiometric uranium dioxide in the presence of uranium vacancies or divacancies constituted by one uranium and one oxygen vacancy. Moreover helium atoms induce a lattice parameter contraction except in interstitial sites where a slight expansion is calculated. (author)
Atomic Structure of Benzene Which Accounts for Resonance Energy
Heyrovska, Raji
2008-01-01
Benzene is a hexagonal molecule of six carbon atoms, each of which is bound to six hydrogen atoms. The equality of all six CC bond lengths, despite the alternating double and single bonds, and the surplus (resonance) energy, led to the suggestion of two resonanting structures. Here, the new atomic structure shows that the bond length equality is due to three carbon atoms with double bond radii bound to three other carbon atoms with resonance bond radii (as in graphene). Consequently, there are two kinds of CH bonds of slightly different lengths. The bond energies account for the resonance energy.
Atomic and magnetic structure of manganites
Complete text of publication follows. Perovskite manganites A1-xA'xMnO3 (A=La, Pr, ..., A' = Ca, Sr,...) have become a hot problem for neutron scattering during last three years. A long story of their neutron diffraction studies, started by a classical work of Wollan and Koehler [1] on La1-xCaxMnO3 and continued by a detailed investigation of Pr1-xCaxMnO3 by Jirak et al. [2], finalised with a burst-like interest in the middle of ninetieths which was provoked by rediscovering of CMR-effect in these compounds. In last three years numerous studies have been performed in all neutron scattering centres, which helped physicists to understand the main features of CMR materials. The most important of them are the strong correlation between transport properties and changes in atomic (charge/orbital ordering) and magnetic (FM/AFM ordering) structure and their dependence on the doping level, the mean A-site canon radius and its variance, as well as cation off-stoichiometry in the compound. Neutron scattering helped also to understand the important role of a local disorder and unusually strong electron-phonon coupling for the properties of manganites. (author) [1] E.O. Wollan and W.C. Koehler, Phys. Rev. 100 (1995) 545. [2] Z. Jirak, S. Krupicka, Z. Simsa, M. Dlouha and S. Vratislav, J. Magn. Magn. Mater. 53 (1985) 153
Computer Simulation of Atoms Nuclei Structure Using Information Coefficients of Proportionality
Labushev, Mikhail M
2012-01-01
The latest research of the proportionality of atomic weights of chemical elements made it possible to obtain 3 x 3 matrices for the calculation of information coefficients of proportionality Ip that can be used for 3D modeling of the structure of atom nucleus. The results of computer simulation show high potential of nucleus structure research for the characterization of their chemical and physical properties.
Analytic methods for the calculation of the electronic structure of solids
Andersen's atomic-sphere approximation has been utilized with approximations based upon linear-combination of atomic orbitals (LCAO) theory to obtain approximate energy-band parameters for solids. Simple analytic expressions for the bandwidth and position of the band center have been derived that require only free-atom wave functions evaluated at the Wigner-Seitz atomic-sphere radius. For convenience, the method has been named the atomic surface method (ASM). The following simple analytic expressions for the band parameters have been derived from the ASM: (i) The bandwidth is equal to the product of h2/m, the gradient of the electron density at the atomic-sphere radius, and the surface area of the sphere; (ii) the average band energy is shifted from the atomic-term-value energy by an amount given by the product of the bandwidth, electron density at the atomic-sphere radius, and atomic-sphere volume. The theory has been applied without adjustable parameters to the transition metals and f-shell metals with use of tabulated Hartree-Fock wave functions and is in reasonable agreement with full band-structure calculations. The same analysis is applied to atomic core states under compression and is also in reasonable agreement with complete band-structure calculations. The 2s and 2p states of Na and Al have been calculated to the point where they merge with the conduction band as free-electron states. These bandwidths and shifts are also written in terms of the atomic term values by using the asymptotic form of the radial wave function. Finally, the LCAO energy bands of Ni are calculated with use of the ASM parameters
Non-perturbative calculations for the multiphoton ionization of hydrogen and lithium atoms
Multiphoton ionization rates for the Hydrogen atom are calculated by direct solution of the time-dependent Schrodinger equation for several intensities at a photon energy of 5.0 eV (KrF laser). Ionization rates for linear polarized light are extracted front the time evolution of the ground state on a 2d cylindrical coordinate lattice, while rates for circular polarized light are extracted from calculations on a 3d Cartesian coordinate lattice. Multiphoton ionization rates for the Lithium atom are calculated in the frozen-core TDHF approximation for a variety of intensities and photon frequencies. The time-dependent equation for the valence HF orbital is solved on a 2d cylindrical coordinate lattice using both fixed and variable grid spacings. The non-perturbative results for both atoms are in sharp contrast to perturbation theory predictions
Spin-splitting calculation for zincblende semiconductors using an atomic bond-orbital model
We develop a 16-band atomic bond-orbital model (16ABOM) to compute the spin splitting induced by bulk inversion asymmetry in zincblende materials. This model is derived from the linear combination of atomic-orbital (LCAO) scheme such that the characteristics of the real atomic orbitals can be preserved to calculate the spin splitting. The Hamiltonian of 16ABOM is based on a similarity transformation performed on the nearest-neighbor LCAO Hamiltonian with a second-order Taylor expansion over k-vector at the Γ point. The spin-splitting energies in bulk zincblende semiconductors, GaAs and InSb, are calculated, and the results agree with the LCAO and first-principles calculations. However, we find that the spin-orbit coupling between bonding and antibonding p-like states, evaluated by the 16ABOM, dominates the spin splitting of the lowest conduction bands in the zincblende materials.
Spin-splitting calculation for zincblende semiconductors using an atomic bond-orbital model
Kao, Hsiu-Fen; Lo, Ikai; Chiang, Jih-Chen; Chen, Chun-Nan; Wang, Wan-Tsang; Hsu, Yu-Chi; Ren, Chung-Yuan; Lee, Meng-En; Wu, Chieh-Lung; Gau, Ming-Hong
2012-10-01
We develop a 16-band atomic bond-orbital model (16ABOM) to compute the spin splitting induced by bulk inversion asymmetry in zincblende materials. This model is derived from the linear combination of atomic-orbital (LCAO) scheme such that the characteristics of the real atomic orbitals can be preserved to calculate the spin splitting. The Hamiltonian of 16ABOM is based on a similarity transformation performed on the nearest-neighbor LCAO Hamiltonian with a second-order Taylor expansion over \\vec{k} at the Γ point. The spin-splitting energies in bulk zincblende semiconductors, GaAs and InSb, are calculated, and the results agree with the LCAO and first-principles calculations. However, we find that the spin-orbit coupling between bonding and antibonding p-like states, evaluated by the 16ABOM, dominates the spin splitting of the lowest conduction bands in the zincblende materials.
A new approach to the structure of the lightest atoms
Full text: Negative ions of the lightest nuclei are often produced in ion sources of modern linear accelerators. Apart from that H- machines are currently used for producing radiopharmaceuticals. The structure of H- atom is missing in many textbooks on quantum mechanics. The H- atom is one of the simplest three-body (or four-body) systems and it is tempting to use the variation method, which leads to the classical formula for the ground state energy of Helium-like atoms: Eo(Z)= -27.2(Z-5/16)2. However, variation methods with 3, 6, and 11 parameters were used in the case of H-. Using the formula, one takes into account the screening of the atomic charge by one electron from the point of view of the other. Nevertheless, the trial function used in classical calculations exhibits the same probability of finding two electrons at the same point as e.g. in two points, where the nucleus in midway between them. This is the reason why one is tempted to modify the classical trial function by the factor s(Z) (r1-r2), in the hope to avoid the necessity of introducing so many variation parameters. In eq.(2) we introduced s(Z) for the new variation parameter and r1 and r2 are position vectors of the two electrons. The expression (2) is zero, when the two electrons coincide. The optimal value of s(1) turns out to be s(1)=0.03037 for Z=1, while for heavier helium-like atoms the optimal s(Z) will be negative. In the case of Z=10, the parameter s(10) reaches almost its asymptotic value Sasmp= -0.0541. The expression for the ground state energy of Helium-like atoms Eo,s(Z) is now, of course, more complicated than eq.(1). In the case of H- we obtain Eo,s(1)=-12.82762 eV, the value slightly smaller than the binding energy of the neutral Hydrogen atom, EH=-13.5985 eV. Since Eo,s(1) is not directly accessible in experiment, we have to start up with the reactions like: H2 + e→H-+H or H2 + e→ H-+H+e. They provide us in fact with the binding energy of H- in the molecule H2, which turns
Atomic and electronic structures of montmorillonite in soft rock
Montmorillonite is a kind of clay mineral which often causes large deformation in soft-rock tunnel engineering and thus brings about safety problems in practice. To deal with these engineering safety problems, the physical and chemical properties of montmorillonite should be studied from basic viewpoints. We study the atomic and electronic structures of montmorillonite by using density-functional theory within the local-density approximation (LDA). The results of calculation show that Al–O bond lengths are longer than Si–O bond lengths. It is found that both the valence band maximum (VBM) and the conduction band minimum (CBM) of montmorillonite are at point Γ, and the calculated direct band gap of montmorillonite is 5.35 eV. We show that the chemical bonding between cations and oxygen anions in montmorillonite is mainly ionic, accompanied as well by a minor covalent component. It is pointed out that the VBM and CBM of montmorillonite consist of oxygen 2p and cation s states, respectively. Our calculated results help to understand the chemical and physical properties of montmorillonite, and are expected to be a guide for solving the problem of large deformation of soft-rock tunnels. (condensed matter: electronic structure, electrical, magnetic, and optical properties)
Energy losses of fast heavy multiply charged structural ions in collisions with complex atoms
Matveev, V. I.; Sidorov, D. B.
2007-07-01
A nonperturbatve theory of energy losses of fast heavy multiply charged structural ions in collisions with neutral complex atoms is elaborated with allowance for simultaneous excitations of ionic and atomic electron shells. Formulas for the effective deceleration that are similar to the well-known Bethe-Bloch formulas are derived. By way of example, the energy lost by partially stripped U q+ ions (10 ≤ q ≤ 70) colliding with argon atoms and also the energy lost by Au, Pb, and Bi ions colliding with various targets are calculated. The results of calculation are compared with experimental data.
The effective atomic numbers of some biomolecules calculated by two methods: A comparative study
Manohara, S.R.; Hanagodimath, S.M.; Gerward, Leif
2009-01-01
The effective atomic numbers Z(eff) of some fatty acids and amino acids have been calculated by two numerical methods, a direct method and an interpolation method, in the energy range of 1 keV-20 MeV. The notion of Z(eff) is given a new meaning by using a modern database of photon interaction cross...... constant and equal to the mean atomic number of the material. Wherever possible, the calculated values of Z(eff) are compared with experimental data....
Band Structure Calculation of MnxCoyFe3-x-yO4
Rosenson, A.; Tailhades, Ph.
1997-01-01
Electronic band structure of MnxCoyFe3-x-yO4 has been calculated in high symmetrical points and lines of the first Brillouin zone within the scope one-electron quasirelativistic pseudopotential approach. Atomic potential form-factors have been calculated in accordance with modified LCAO method. Dependence of energy gap Eg=Eg(x,y) against Mn, Co and Fe concentrations is calculated and presented.
Theoretical development of atomic structure: Past, present and future
Theoretical development of atomic structure is briefly discussed. The role of correlation, relativity, quantum electrodynamic (QED), finite nuclear size (FNS) and parity nonconservation (PNC) in high precision theoretical investigation of properties of atomic and ionic systems is demonstrated. At present, we do not have a comprehensive and practical atomic structure theory which accounts all these physical effects on an equal footing. Suggestions are made for future directions. (author). 108 refs, 5 figs, 9 tabs
The generalized sturmian method for calculating spectra of atoms and ions
Avery, James Emil; Avery, John Scales
2003-01-01
The properties of generalized Sturmian basis sets are reviewed, and functions of this type are used to perform direct configuration interaction calculations on the spectra of atoms and ions. Singlet excited states calculated in this way show good agreement with experimentally measured spectra. When...... Slater–Condon rules must be used. This aspect of the problem is discussed in detail. Finally spectra are calculated in the presence of a strong external electric field. In addition to the expected Stark effect, the calculated spectra exhibit anomalous states. These are shown to be states where one of the...
Determination of atomic cluster structure with cluster fusion algorithm
Obolensky, Oleg I.; Solov'yov, Ilia; Solov'yov, Andrey V.; Greiner, Walter
We report an efficient scheme of global optimization, called cluster fusion algorithm, which has proved its reliability and high efficiency in determination of the structure of various atomic clusters.......We report an efficient scheme of global optimization, called cluster fusion algorithm, which has proved its reliability and high efficiency in determination of the structure of various atomic clusters....
Full text: To understand the plasma-wall interaction on divertor plates, we investigate the interaction of hydrogen atoms and carbon materials used in the high heat flux components by the use of the following simulations. Monte-Carlo (MC) method based on binary collision approximation can calculate the sputtering process of hydrogen atoms on the carbon material quickly. Classical molecular dynamics (MD) method employs multi-body potential models and can treat realistic structures of crystal and molecule. The ab-initio method can calculate electron energy in quantum mechanics, which is regarded as realistic potential for atoms. In the present paper, the interaction of the hydrogen and the carbon material is investigated using the multi-scale (MC, MD and ab-initio) methods. The bombardment of hydrogen atoms onto the carbon material is simulated by the ACAT-code of the MC method, which cannot represent the structure of crystal, and the MD method using modified reactive empirical bond order (REBO) potential, which treats single crystal graphite and amorphous carbon. Consequently, we clarify that the sputtering yield and the reflection rate calculated by the ACAT-code agree with those on the amorphous carbon calculated by the MD. Moreover, there are many kinds of REBO potential for the MD. Adsorption, reflection and penetration rates between a hydrogen atom and a graphene surface are calculated by the MD simulations using the two kinds of potential model. For the incident energy of less than 1 eV, the MD simulation using the modified REBO potential, which is based on Brenner's REBO potential in 2002, shows that reflection is dominant, while the most popular Brenner's REBO potential in 1990 shows that adsorption is dominant. This reflection of the low energy injection is caused by a small potential barrier for the hydrogen atom in the modified REBO potential. The small potential barrier is confirmed by the ab-initio calculations, which are hybrid DFT (B3LYP/cc-pVDZ), ab
The stability and electronic structure of Fe atoms embedded in zigzag graphene nanoribbons
The stability and electronic properties of the Fe atoms embedded in divacancy defects in graphene nanoribbons (GNR) with zigzag-shaped edges have been studied by first-principles calculations. When Fe is positioned in the middle of the ribbon, it has little effect on the edge C atoms, which reserves the flat edges of graphene nanoribbons. On the other hand, when Fe atom is near the edge, structural distortion takes place resulting in tilted-edge structure with low energies. This indicates that the Fe atoms prefer to occupy divacancy sites near the edges. This is also in consistent with the analyses of electronic structures. Meanwhile, our results reveal that embedding Fe atom in the graphene nanoribbons is an effective method to make the GNR possessing metallic properties.
Magnetism, microstructure and First Principles calculations of atomized and annealed Ni{sub 3}Al
García-Escorial, A., E-mail: age@cenim.csic.es [CENIM-CSIC, Avda. Gregorio del Amo, 8, 28040 Madrid (Spain); Crespo, P.; Hernando, A. [Instituto de Magnetismo Aplicado, IMA-UCM, P.O. Box 155, 28230 Madrid (Spain); Lieblich, M. [CENIM-CSIC, Avda. Gregorio del Amo, 8, 28040 Madrid (Spain); Marín, P.; Velasco, V. [Instituto de Magnetismo Aplicado, IMA-UCM, P.O. Box 155, 28230 Madrid (Spain); Ynduráin, F. [Dpto. de Física de la Materia Condensada, UAM, Cantoblanco, 28049 Madrid (Spain)
2014-12-05
Highlights: • The microstructure and order of as-atomized Ni{sub 3}Al powder change with annealing. • The change of the magnetic properties shows the influence of the chemical order. • First Principles calculations show the effect of the density of states to the order. - Abstract: In this work Ni{sub 3}Al powder particles obtained by atomization were characterized magnetically and microstructurally in as-atomized state and after annealing. Upon annealing the X-ray diffraction patterns show a noticeable increase of the signal of the ordered phase γ′-Ni{sub 3}Al, L1{sub 2}, phase and the microstructure evolves from a lamellar and dendrite to a large grain microstructure. The Curie temperature of the as-atomized powder particles is 85 K and decreases after annealing down to 50 K. First Principles calculations were carried out to correlate the experimental observations with local order of Ni and Al atoms and illustrate the importance of the local order in the density of states at the Fermi level, showing how the magnetic moment depends on the Ni and Al atomic position.
Calculations with the quasirelativistic local-spin-density-functional theory for high-Z atoms
The generalized-exchange local-spin-density-functional theory (LSD-GX) with relativistic corrections of the mass velocity and Darwin terms has been used to calculate statistical total energies for the neutral atoms, the positive ions, and the negative ions for high-Z elements. The effect of the correlation and relaxation correction on the statistical total energy is discussed. Comparing the calculated results for the ionization potentials and electron affinities for the atoms (atomic number Z from 37 to 56 and 72 to 80) with experiment, shows that for the atoms rubidium to barium both the LSD-GX and the quasirelativistic LSD-GX, with self-interaction correction, Gopinathan, Whitehead, and Bogdanovic's Fermi-hole parameters [Phys. Rev. A 14, 1 (1976)], and Vosko, Wilk, and Nusair's correlation correction [Can. J. Phys. 58, 1200 (1980)], are very good methods for calculating ionization potentials and electron affinities. For the atoms hafnium to mercury the relativistic effect has to be considered
STARK STRUCTURE OF THE RYDBERG STATES OF ALKALINE-EARTH ATOMS
郅妙婵; 戴长建; 李士本
2001-01-01
The Stark effects of the Rydberg states in the alkaline-earth atoms are studied theoretically. Using a method similar to the treatment of alkali atoms, the properties of the Stark states of Mg, Ca, Sr and Ba atoms in the regions far away from the perturbers are investigated. The Stark maps for Mg (n=16, M=0), Ca (n=10, M=0), Sr (n=12,M=0) and Ba (n=13, |M|=0,1) are presented. Topics such as the general methods of calculation, the treatment of fine structure, and the structure of level anti-crossings are discussed. The comparison between the theoretical and experimental Stark maps is satisfactory.
Using Density Functional Theory (DFT) for the Calculation of Atomization Energies
Bauschlicher, Charles W., Jr.; Partridge, Harry; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
The calculation of atomization energies using density functional theory (DFT), using the B3LYP hybrid functional, is reported. The sensitivity of the atomization energy to basis set is studied and compared with the coupled cluster singles and doubles approach with a perturbational estimate of the triples (CCSD(T)). Merging the B3LYP results with the G2(MP2) approach is also considered. It is found that replacing the geometry optimization and calculation of the zero-point energy by the analogous quantities computed using the B3LYP approach reduces the maximum error in the G2(MP2) approach. In addition to the 55 G2 atomization energies, some results for transition metal containing systems will also be presented.
Post, D; Clark, R E H; Putvinskaya, N
1995-01-01
Reduction of the peak heat loads on the plasma facing components is essential for the success of the next generation of high fusion power tokamaks such as the International Thermonuclear Experimental Reactor (ITER) 1 . Many present concepts for accomplishing this involve the use of atomic processes to transfer the heat from the plasma to the main chamber and divertor chamber walls and much of the experimental and theoretical physics research in the fusion program is directed toward this issue. The results of these experiments and calculations are the result of a complex interplay of many processes. In order to identify the key features of these experiments and calculations and the relative role of the primary atomic processes, simple quasi-analytic models and the latest atomic physics rate coefficients and cross sections have been used to assess the relative roles of central radiation losses through bremsstrahlung, impurity radiation losses from the plasma edge, charge exchange and hydrogen radiation losses f...
Calculation of parity and time invariance violation in the radium atom
Dzuba, V A; Ginges, J S M
1999-01-01
Parity (P) and time (T) invariance violating effects in the Ra atom are strongly enhanced due to close states of opposite parity, the large nuclear charge Z and the collective nature of P,T-odd nuclear moments. We have performed calculations of the atomic electric dipole moments (EDM) produced by the electron EDM and the nuclear magnetic quadrupole and Schiff moments. We have also calculated the effects of parity non-conservation produced by the nuclear anapole moment and the weak charge. Our results show that as a rule the values of these effects are much larger than those considered so far in other atoms (enhancement is up to 10^5 times).
Evaluation and Comparison of the Configuration Interaction Calculations for Complex Atoms
Charlotte Froese Fischer
2014-03-01
Full Text Available Configuration interaction (CI methods are the method of choice for the determination of wave functions for complex atomic systems from which a variety of atomic properties may be computed. When applied to highly ionized atoms, where few, if any, energy levels from observed wavelengths are available, the question arises as to how a calculation may be evaluated. Many different codes are available for such calculations. Agreement between the results from different codes in itself is not a check on accuracy, but may be due to a similarity in the computational procedures. This paper reviews basic theory, which, when applied in a systematic manner, can be the basis for the evaluation of accuracy. Results will be illustrated in the study of 4s24p5 (odd and 4s24p44d (even levels in W39+ and the transitions between them.
Atomic-absorption programmes for the Hewlett Packard H.P. 97S programmable calculator
This report outlines the limitations of the atomic-absorption (AA) programmes supplied with the Hewlett Packard H.P.97S programmable calculator, and proposes ways in which these limitations can be overcome. Three new programmes, for AA on-line analysis, off-line analysis and an automatic system, are described
A full CI treatment of Ne atom - A benchmark calculation performed on the NAS CRAY 2
Bauschlicher, C. W., Jr.; Langhoff, S. R.; Partridge, H.; Taylor, P. R.
1986-01-01
Full CI calculations are performed for Ne atom using Gaussian basis sets of up to triple-zeta plus double polarization quality. The total valence correlation energy through double, triple, quadruple and octuple excitations is compared for eight different basis sets. These results are expected to be an important benchmark for calibrating methods for estimating the importance of higher excitations.
CALCULATED PROPERTIES OF TWO-DIMENSIONAL SPIN-POLARIZED ATOMIC HYDROGEN
Lantto, L.; Nieminen, R.
1980-01-01
Optimal HNC-Jastrow calculations have been carried out for gaseous spin-polarized hydrogen in two space dimensions. Accurate values for the ground state energy, radial distribution function, average exchange energy and momentum distribution are obtained at low atomic densities.
C and N atoms are the most frequent foreign interstitial atoms (FIAs), and often incorporated into the surface layers of steels to enhance their properties by thermochemical treatments. Al, Si, Ti, V, Cr, Mn, Co, Ni, Cu, Nb and Mo are the most common alloying elements in steels, also can be called foreign substitutional atoms (FSAs). The FIA and FSA interactions play an important role in the diffusion of C and N atoms, and the microstructures and mechanical properties of surface modified layers. Ab initio calculations based on the density functional theory are carried out to investigate FIA interactions with FSA in ferromagnetic bcc iron. The FIA–FSA interactions are analyzed systematically from five aspects, including interaction energies, density of states (DOS), bond populations, electron density difference maps and local magnetic moments
Boll, Torben
2012-10-01
In this article the Cu-Au binding energy in Cu3Au is determined by comparing experimental atom probe tomography (APT) results to simulations. The resulting bonding energy is supported by density functional theory calculations. The APT simulations are based on the Müller-Schottky equation, which is modified to include different atomic neighborhoods and their characteristic bonds. The local environment is considered up to the fifth next nearest neighbors. To compare the experimental with simulated APT data, the AtomVicinity algorithm, which provides statistical information about the positions of the neighboring atoms, is applied. The quality of this information is influenced by the field evaporation behavior of the different species, which is connected to the bonding energies. © Microscopy Society of America 2012.
Atomic and electronic structure of surfaces theoretical foundations
Lannoo, Michel
1991-01-01
Surfaces and interfaces play an increasingly important role in today's solid state devices. In this book the reader is introduced, in a didactic manner, to the essential theoretical aspects of the atomic and electronic structure of surfaces and interfaces. The book does not pretend to give a complete overview of contemporary problems and methods. Instead, the authors strive to provide simple but qualitatively useful arguments that apply to a wide variety of cases. The emphasis of the book is on semiconductor surfaces and interfaces but it also includes a thorough treatment of transition metals, a general discussion of phonon dispersion curves, and examples of large computational calculations. The exercises accompanying every chapter will be of great benefit to the student.
Fluid-structure interaction calculations using a linear perturbation method
Aim of the work is to present and validate FSI (Fluid-Structure Interaction) calculations by using a linear perturbation method and commercial Computational Fluid Dynamics (CFD) and structural analysis codes. Star-CD is used for CFD calculations and ABAQUS for structural analysis. The external MpCCI code is used for coupling the CFD and structural analysis codes
Large Scale Electronic Structure Calculations using Quantum Chemistry Methods
Scuseria, Gustavo E.
1998-03-01
This talk will address our recent efforts in developing fast, linear scaling electronic structure methods for large scale applications. Of special importance is our fast multipole method( M. C. Strain, G. E. Scuseria, and M. J. Frisch, Science 271), 51 (1996). (FMM) for achieving linear scaling for the quantum Coulomb problem (GvFMM), the traditional bottleneck in quantum chemistry calculations based on Gaussian orbitals. Fast quadratures(R. E. Stratmann, G. E. Scuseria, and M. J. Frisch, Chem. Phys. Lett. 257), 213 (1996). combined with methods that avoid the Hamiltonian diagonalization( J. M. Millam and G. E. Scuseria, J. Chem. Phys. 106), 5569 (1997) have resulted in density functional theory (DFT) programs that can be applied to systems containing many hundreds of atoms and ---depending on computational resources or level of theory-- to many thousands of atoms.( A. D. Daniels, J. M. Millam and G. E. Scuseria, J. Chem. Phys. 107), 425 (1997). Three solutions for the diagonalization bottleneck will be analyzed and compared: a conjugate gradient density matrix search (CGDMS), a Hamiltonian polynomial expansion of the density matrix, and a pseudo-diagonalization method. Besides DFT, our near-field exchange method( J. C. Burant, G. E. Scuseria, and M. J. Frisch, J. Chem. Phys. 105), 8969 (1996). for linear scaling Hartree-Fock calculations will be discussed. Based on these improved capabilities, we have also developed programs to obtain vibrational frequencies (via analytic energy second derivatives) and excitation energies (through time-dependent DFT) of large molecules like porphyn or C_70. Our GvFMM has been extended to periodic systems( K. N. Kudin and G. E. Scuseria, Chem. Phys. Lett., in press.) and progress towards a Gaussian-based DFT and HF program for polymers and solids will be reported. Last, we will discuss our progress on a Laplace-transformed \\cal O(N^2) second-order pertubation theory (MP2) method.
Mikajlo, E A; Ford, M J
2003-01-01
This paper presents an experimental measurement of the electronic structure of Na sub 2 O in the solid phase using electron momentum spectroscopy and compares the results with ab initio calculations performed within the linear combination of atomic orbitals (LCAO) approximation. While Hartree-Fock (HF) can reproduce elastic properties we find it overestimates splitting of the oxygen valence bands by around 30% and the width of the O 2p band by a factor of 2. Our experimental values are 15.85 +- 0.2 and 0.6 +- 0.2 eV for these two quantities, respectively. Density functional methods are significantly better, with the hybrid functional PBE0 predicting the oxygen bandgap to within the experimental error. PBE0 also gives the best estimate of the Na core level energies. In contrast, HF performs best for the splitting between the oxygen and sodium bands. Our experimental values of 32.85 +- 0.2 and 27.45 +- 0.2 eV for the Na 2p-Na 2s and O 2p-Na 2p splittings agree well with previous measurements. Distribution of el...
Theoretical calculations of electron-impact and radiative processes in atoms
Electron-impact and radiative processes in atoms are investigated with particular attention paid to the effects of electron correlations. Using the optical potential method, the cross section for the elastic scattering of electrons by the neutral argon atom is calculated from 0 to 300 eV. Corrections to the Hartree--Fock cross section are obtained from a many-particle perturbation expansion. The effects of electron correlations are found to be quite significant at low energy. The optical potential results are compared with a polarized orbital calculation, the Born approximation and experiment. The 2s and 2p excitation cross sections for electron scattering on hydrogen are calculated by two similar methods. The distorted wave method is applied and the effect of calculating the outgoing scattered electron in the potential of the initial or final state is investigated. The imaginary part of the optical potential is also calculated in lowest order by the use of many-body diagrams. The subshell photoionization cross sections in argon are calculated using the acceleration, length and velocity forms of the dipole operator. First order electron correlation corrections to the Hartree--Fock approximation are obtained through the use of many-body perturbation theory. Also investigated is the two photon ionization cross section for the neutral argon atom. A double perturbation expansion in the Coulomb correlations and the atom-radiation field interaction is made. Contributions from intermediate states are obtained by direct summation over Hartree--Fock bound and continuum single particle states. The effects of electron correlations and photon radiative corrections are investigated
Atomic displacements in ferroelectric trigonal and orthorhombic boracite structures
Dowty, Eric; Clark, J.R.
1972-01-01
New crystal-structure refinements of Pca21 boracite, Mg3ClB7O13, and R??{lunate}c ericaite, Fe2.4Mg0.6ClB7O13, show that some boron and oxygen atoms are involved in the 'ferro' transitions as well as the metal and halogen atoms. The atomic displacements associated with the polarity changes are as large as 0.6A??. ?? 1972.
Electronic Structure of Helium Atom in a Quantum Dot
Jayanta, K. Saha; Bhattacharyya, S.; T. K., Mukherjee
2016-03-01
Bound and resonance states of helium atom have been investigated inside a quantum dot by using explicitly correlated Hylleraas type basis set within the framework of stabilization method. To be specific, precise energy eigenvalues of bound 1sns (1Se) (n = 1-6) states and the resonance parameters i.e. positions and widths of 1Se states due to 2sns (n = 2-5) and 2pnp (n = 2-5) configurations of confined helium below N = 2 ionization threshold of He+ have been estimated. The two-parameter (Depth and Width) finite oscillator potential is used to represent the confining potential due to the quantum dot. It has been explicitly demonstrated that the electronic structural properties become sensitive functions of the dot size. It is observed from the calculations of ionization potential that the stability of an impurity ion within a quantum dot may be manipulated by varying the confinement parameters. A possibility of controlling the autoionization lifetime of doubly excited states of two-electron ions by tuning the width of the quantum cavity is also discussed here. TKM Gratefully Acknowledges Financial Support under Grant No. 37(3)/14/27/2014-BRNS from the Department of Atomic Energy, BRNS, Government of India. SB Acknowledges Financial Support under Grant No. PSW-160/14-15(ERO) from University Grants Commission, Government of India
Accelerating structure optimization and tolerance calculation
The problem of optimizing the dynamics of a charged particle beam with high volumetric charge density in an accelerating structure with space-homogeneous quadrupole focusing (SHQF) is considered. The mathematical model of interacting particle dynamics is based on equations for beam envelopes. The optimization criterion (functional) is selected from the terms of obtaining a structure with an assigned length, beam energy at the accelerator outlet and the maximum particle capture for acceleration mode. A method for determining tolerance for the structure parameters based on optimization algorithm is proposed for the structure with SHQF. 4 refs.; 2 figs
Wagner, Philipp; Ewels, Christopher P.; Suarez-Martinez, Irene; Guiot, Vincent; Cox, Stephen F. J.; Lord, James S.; Briddon, Patrick R.
2011-01-01
We examine the behaviour of hydrogen ions, atoms and molecules in alpha-boron using density functional calculations. Hydrogen behaves as a negative-U centre, with positive H ions preferring to sit off-center on inter-layer bonds and negative H ions sitting preferably at in-plane sites between three B12 icosahedra. Hydrogen atoms inside B12 icosahedral cages are unstable, drifting off-center and leaving the cage with only a 0.09 eV barrier. While H0 is extremely mobile (diffusion barrier 0.25 ...
Sharma, Pragati; Roy, Sudip; Karimi-Varzaneh, Hossein Ali
2016-02-25
Microsecond atomic-scale molecular dynamics simulation has been employed to calculate the glass-transition temperature (Tg) of cis- and trans-1,4-polybutadiene (PB) and 1,4-polyisoprene (PI). Both all-atomistic and united-atom models have been simulated using force fields, already available in literature. The accuracy of these decade old force fields has been tested by comparing calculated glass-transition temperatures to the corresponding experimental values. Tg depicts the phase transition in elastomers and substantially affects various physical properties of polymers, and hence the reproducibility of Tg becomes very crucial from a thermodynamic point of view. Such validation using Tg also evaluates the ability of these force fields to be used for advanced materials like rubber nanocomposites, where Tg is greatly affected by the presence of fillers. We have calculated Tg for a total of eight systems, featuring all-atom and united-atom models of cis- and trans-PI and -PB, which are the major constituents of natural and synthetic rubber. Tuning and refinement of the force fields has also been done using quantum-chemical calculations to obtain desirable density and Tg. Thus, a set of properly validated force fields, capable of reproducing various macroscopic properties of rubber, has been provided. A novel polymer equilibration protocol, involving potential energy convergence as the equilibration criterion, has been proposed. We demonstrate that not only macroscopic polymer properties like density, thermal expansion coefficient, and Tg but also local structural characteristics like end-to-end distance (R) and radius of gyration (Rg) and mechanical properties like bulk modulus have also been equilibrated using our strategy. Complete decay of end-to-end vector autocorrelation function with time also supports proper equilibration using our strategy. PMID:26836395
Single structured light beam as an atomic cloud splitter
We propose a scheme to split a cloud of cold noninteracting neutral atoms based on their dipole interaction with a single structured light beam which exhibits parabolic cylindrical symmetry. Using semiclassical numerical simulations, we establish a direct relationship between the general properties of the light beam and the relevant geometric and kinematic properties acquired by the atomic cloud as it passes through the beam.