Arxiv:2104.08286V1 [Astro-Ph.IM] 16 Apr 2021 Ture Data

Research Note version April 20, 2021 Typeset using LATEX twocolumn style in AASTeX631

linemake: an Atomic and Molecular Line List Generator

Vinicius M. Placco ,1 Christopher Sneden ,2 Ian U. Roederer ,3, 4 James E. Lawler ,5 Elizabeth A. Den Hartog ,5 Neda Hejazi ,6 Zachary Maas ,2 and Peter Bernath 7, 8

1Community Science and Data Center/NSF’s NOIRLab, 950 N. Cherry Ave., Tucson, AZ 85719, USA 2Department of Astronomy and McDonald Observatory, The University of Texas, Austin, TX 78712, USA 3Department of Astronomy, University of Michigan, 1085 S. University Ave., Ann Arbor, MI 48109, USA 4Joint Institute for Nuclear Astrophysics – Center for the Evolution of the Elements (JINA-CEE), USA 5Department of Physics, University of Wisconsin-Madison, 1150 University Ave., Madison, WI 53706 6Department of Physics and Astronomy, Georgia State University, Atlanta, GA 30302, USA 7Department of Chemistry and Biochemistry, Old Dominion University, Norfolk, VA 23529, USA 8Department of Physics, Old Dominion University, Norfolk, VA 23529, USA

ABSTRACT In this research note, we present linemake, an open-source atomic and molecular line list generator. Rather than a replacement for a number of well-established atomic and molecular spectral databases, linemake aims to be a lightweight, easy-to-use tool to generate formatted and curated lists suitable for spectral synthesis work. We encourage users of linemake to understand the sources of their transition data and cite them as appropriate in published work. We provide the code, line database, and an exten- sive list of literature references in a GitHub repository (https://github.com/vmplacco/linemake), which will be updated regularly as new data become available.

Keywords: Spectroscopy (1558), Atomic physics (2063), Laboratory astrophysics (2004), Molecular physics (2058), Spectral line lists (2082)

1. INTRODUCTION all Fe- and lanthanide-group elements, and recently a Stellar and Galactic chemical evolution can be illu- study of Ca I initiates work on lighter elements. minated by abundance analyses, but such studies must With increased interest in low-temperature stars strive for abundance accuracy. This effort depends on (spectral types mid-K and cooler) it has become at- many factors, but none is more crucial than access to tractive to include molecules such as TiO and H2O in trustworthy basic atomic and molecular transition data. linemake. These molecules have substantial numbers of Several excellent large transition databases exist, the transitions available in large databases but are not easy most well-known being VALD1. Here we introduce the to incorporate in synthetic spectrum line lists. They are utility code linemake, created to generate synthetic included in linemake with choices that keep the num- spectrum input line lists by merging Kurucz2 (Kurucz ber of transitions to a manageable size. linemake does 2011) atomic/molecular line compendium information not try to compete with the VALD resource or other more with updated and very accurate transition probabili- comprehensive line compendia. Our effort is to produce ties, hyperfine structure (HFS), and isotopic substruc- a carefully controlled line list of easily referenced sources arXiv:2104.08286v1 [astro-ph.IM] 16 Apr 2021 ture data. The atomic data included in linemake are from laboratory physics. In this document, we outline primarily those that have been published by the Univer- linemake and point to the sources of its data. sity of Wisconsin atomic physics group (e.g. Lawler et al. 2. ABOUT linemake 2009, among others) and the molecular data are from the Old Dominion University molecular physics group linemake produces synthesis line lists compatible with 3 (Bernath 2020). The Wisconsin data now include nearly those needed by the line analysis code MOOG (Sneden 1973). It starts with the Kurucz compendium and then substitutes or supplements these lists with atomic data 1 http://vald.astro.uu.se/ 2 http://kurucz.harvard.edu/linelists.html 3 https://www.as.utexas.edu/~chris/moog.html 2

1 2 I H He 1.00794 4.002602 Hydrogen Helium 3 4 5 6 7 8 9 10 II Li Be B CNOF Ne 6.941 9.012182 10.811 12.0107 14.0067 15.9994 18.9984032 20.1797 Lithium Beryllium Boron Carbon Nitrogen Oxygen Fluorine Neon 11 12 13 14 15 16 17 18 III Na Mg Al Si P S Cl Ar 22.98976928 24.3050 26.9815386 28.0855 30.973762 32.065 35.453 39.948 Sodium Magnesium Aluminum Silicon Phosphorus Sulfur Chlorine Argon 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 IV K Ca Sc Ti V Cr Mn Fe Co Ni Cu Zn Ga Ge As Se Br Kr 39.0983 40.078 44.955912 47.867 50.9415 51.9961 54.938045 55.845 58.933195 58.6934 63.546 65.38 69.723 72.64 74.92160 78.96 79.904 83.798 Potassium Calcium Scandium Titanium Vanadium Chromium Manganese Iron Cobalt Nickel Copper Zinc Gallium Germanium Arsenic Selenium Bromine Krypton 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 V Rb Sr Y Zr Nb Mo TcRu Rh Pd Ag Cd In SnSb Te I Xe 85.4678 87.62 88.90585 91.224 92.90638 95.96 [98] 101.07 102.90550 106.42 107.8682 112.411 114.818 118.710 121.760 127.60 126.90447 131.293 Rubidium Strontium Yttrium Zirconium Niobium Molybdenum Technetium Ruthenium Rhodium Palladium Silver Cadmium Indium Tin Antimony Tellurium Iodine Xenon 55 56 72 73 74 75 76 77 78 79 80 8182 83 84 85 86 VI Cs Ba LA Hf Ta W Re Os Ir Pt Au Hg TlPb Bi Po At Rn 132.9054519 137.327 178.49 180.94788 183.84 186.207 190.23 192.217 195.084 196.966569 200.59 204.3833 207.2 208.98040 [209] [210] [222] Cesium Barium Hafnium Tantalum Tungsten Rhenium Osmium Iridium Platinum Gold Mercury Thallium Lead Bismuth Polonium Astatine Radon 87 88 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 VII Fr Ra AC Rf Db Sg Bh Hs Mt Ds Rg Cn Nh Fl Mc Lv Ts Og [223] [226] [267] [268] [271] [272] [270] [276] [281] [280] [285] [284] [289] [288] [293] [294] [294] Francium Radium Rutherfordium Dubnium Seaborgium Bohrium Hassium Meitnerium Darmstadtium Roentgenium Copernicium Nihomium Flerovium Moscovium Livermorium Tennessine Oganesson

57 58 59 6061 62 63 64 65 66 67 68 69 70 71 La La Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb Lu 138.90547 140.116 140.90765 144.242 [145] 150.36 151.964 157.25 158.92535 162.500 164.93032 167.259 168.93421 173.054 174.9668 Lanthanum Cerium Praseodymium Neodymium Promethium Samarium Europium Gadolinium Terbium Dysprosium Holmium Erbium Thulium Ytterbium Lutetium 8990 91 92 93 94 95 96 97 98 99 100 101 102 103 Ac AcTh Pa U Np Pu Am Cm Bk Cf Es Fm Md No Lr [227] 232.03806 231.03588 238.02891 [237] [244] [243] [247] [247] [251] [252] [257] [258] [259] [262] Actinium Thorium ProtactiniumUranium Neptunium Plutonium Americium Curium Berkelium Californium Einsteinium Fermium Mendelevium Nobelium Lawrencium 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

not available neutral and ion ionized species neutral species molecular species

Figure 1. Elements with curated transitions currently available in linemake. The most up-to-date version of this periodic table is available at the GitHub repository. from the Wisconsin group and molecular data from the sented by multiple isotopes in its transitions, linemake Old Dominion group. linemake is written in standard presents the lists with isotopic identification but no as- FORTRAN, and it can be compiled and executed on all sumptions about isotopic ratios. For a given transition, UNIX/Linux-based operating systems. the total gf value for each isotope is the same. When linemake is primarily aimed for spectroscopic studies applying spectrum synthesis codes like MOOG to the line of stars cooler than B spectral type. Some limitations lists generated by linemake, users must set the desired in the available line lists have been imposed in order to isotopic fractions in parameter files. Secondly, linemake keep the resulting synthesis lists to sensible sizes. The catalogs some Fe-group hyperfine substructure patterns ionization states available are neutral and first ion. The from the Kurucz database when no recent laboratory maximum lower excitation energy of any line is 7.5 eV, studies are available. These can be included in output except for some transitions of the light elements H, C, files at the user’s option, but they will probably add N, O, Mg, Al, Si, P, and S; higher excitation transitions many extra transitions into the synthesis line lists, and are available for these elements. Additionally, for Fe II, the pedigrees of these are not guaranteed. Users are 8.5 eV is the maximum lower energy. Obvious warnings cautioned to be careful about including them and are should be given about the output line lists: we believe encouraged to test them with stellar syntheses. that they are correct, but there is no substitute for users having a close look to assure themselves of the quality 3. BRIEF DESCRIPTION OF THE DATABASE of these lists. If a user is uncertain about an output Here we briefly discuss the atomic and molecular data line list, the individual files for different species can be sources. The detailed description for each category, in- examined easily. These species files are also useful in cluding literature references and decisions made to max- simply making laboratory-based reliable transition lists imize the utility of line lists for high-resolution spectro- for individual (mostly atomic) species. scopic studies, are given in the README.md file of the A few “local rules” apply when using the code. First GitHub repository. Figure1 summarizes the transition and most importantly, in almost all cases in which an availability in the linemake database, and its most up- atomic (and more often molecular) species is repre- to-date version can also be found in the README.md file. 3

3.1. Atomic Species: Fe-group Elements 3.3. Atomic Species: Other Elements The Fe group is defined here as those elements with Light elements (Z ≤ 20) have often not had recent ex- 21 ≤ Z ≤ 30. Significant HFS usually is a feature of tensive laboratory transition probability studies. Some transitions of odd-Z elements. If the README.md de- progress, both in laboratory and theoretical, has been scription indicates that a species includes HFS, then line made, and a few selected elements are included in the lists have full HFS substructure patterns when they are linemake database. available from laboratory studies. Typically, these are known for many of the lines with laboratory transition probabilities. linemake manages these data internally, 3.4. Molecular (mostly diatomic) Species but users should examine the output synthesis lists to Line data for molecular species have improved con- understand whether HFS patterns have been included siderably recently. Motivated by an interest in M- in any transition of interest. Information about isotopic type stars, which are often involved in exoplanet ra- substructure is scarce. This is largely an issue for even- dial velocity or atmospheric transmission studies, we Z elements, as odd-Z elements have few (often just one) have included molecules that only appear in stellar spec- naturally-occurring isotopes. Moreover, for Fe-group el- tra of very cool stars. The data sources for many ements, usually one isotope dominates the solar system molecules are laboratory studies from Peter Bernath’s abundance, and other isotopes can safely be ignored. group (MoLLIST5). For others we have translated into Exceptions are Ti and Ni, and users should account for MOOG format the line information from two prominent isotopic effects on transitions in the red and infrared molecular physics groups: HITRAN6 and EXOMOL7. spectral regions.

3.2. Atomic Species: Neutron-Capture Elements 4. CONCLUSIONS Neutron-capture elements are defined as those with In this note, we have briefly described the main fea- Z > 30. For lanthanide elements and Hf (57 ≤ Z ≤ 72), tures of the linemake code. This discussion is not the Wisconsin group results dominate, but for lighter intended to be exhaustive, but we hope it may serve and heavier elements other sources are included if they as a helpful guideline for potential users. The GitHub have been subject to recent assessment for reliability. In database will continue to be updated with the most re- this case, NIST refers to the National Institute of Stan- cent laboratory and theoretical work, and input from dards and Technology’s Atomic Spectra Database.4 the user community is certainly welcome.

REFERENCES

Bernath, P. F. 2020, Journal of Quantitative Spectroscopy Lawler, J. E., Sneden, C., Cowan, J. J., Ivans, I. I., & Den and Radiative Transfer, 240, 106687, Hartog, E. A. 2009, ApJS, 182, 51, doi: https://doi.org/10.1016/j.jqsrt.2019.106687 doi: 10.1088/0067-0049/182/1/51 Kurucz, R. L. 2011, Canadian Journal of Physics, 89, 417, Sneden, C. A. 1973, PhD thesis, The University of Texas at doi: 10.1139/p10-104 Austin.

4 https://physics.nist.gov/PhysRefData/ASD/lines form.html 5 http://bernath.uwaterloo.ca/molecularlists.php 6 https://hitran.org/lbl/ 7 http://www.exomol.com/data/molecules/