Search

Dual Coordinate Descent (DCD) is one of the most popular optimization methods. The parallelization of DCD is difficult, as DCD is sequential in nature. As such, simultaneously running multiple DCD threads on batches of data elements causes result inaccuracy and slow convergence, due to the concurrent updates of multiple coordinates. Some parallelization methods adopt separable approximate functions that depend on the degree of parallelism. Such dependencies result in both poor scalability and slow convergence. To address these challenges, in this paper we present a new parallel primal-dual algorithm for DCD, called PPD. In PPD, the block data distribution is utilized to obtain a new approximate function that is independent from the parallelism. Subsequently we design a novel primal-dual acceleration scheme based on PPD to approach the optimal solution closely and quickly. We demonstrate the advantages of PPD in terms of scalability and efficiency through experiments.

New high-level ab initio quartic force field (QFF) methods are explored which provide spectroscopic data for the electronically excited states of the carbon monoxide, water, and formaldehyde cations, sentinel species for expanded, recent cometary spectral analysis. QFFs based on equation-of-motion ionization potential (EOM-IP) with a complete basis set extrapolation and core correlation corrections provide assignment for the fundamental vibrational frequencies of the A˜2B1 and B˜2A1 states of the formaldehyde cation; only three of these frequencies have experimental assignment available. Rotational constants corresponding to these vibrational excitations are also provided for the first time for all electronically excited states of both of these molecules. EOM-IP-CCSDT/CcC computations support tentative re-assignment of the ν1 and ν3 frequencies of the B˜2B2 state of the water cation to approximately 2409.3 cm−1 and 1785.7 cm−1, respectively, due to significant disagreement between experimental assignment and all levels of theory computed herein, as well as work by previous authors. The EOM-IP-CCSDT/CcC QFF achieves agreement to within 12 cm−1 for the fundamental vibrational frequencies of the electronic ground state of the water cation compared to experimental values and to the high-level theoretical benchmarks for variationally-accessible states. Less costly EOM-IP based approaches are also explored using approximate triples coupled cluster methods, as well as electronically excited state QFFs based on EOM-CC3 and the previous (T)+EOM approach. The novel data, including vibrationally corrected rotational constants for all states studied herein, provided by these computations should be useful in clarifying comet evolution or other remote sensing applications in addition to fundamental spectroscopy.

We present the first high spectral resolution mid-infrared survey in the Orion BN/KL region, covering 7.2 to 28.3 micron. With SOFIA/EXES we target the enigmatic source Orion IRc2. While this is in the most prolifically studied massive star-forming region, longer wavelengths and molecular emission lines dominated previous spectral surveys. The mid-infrared observations in this work access different components and molecular species in unprecedented detail. We unambiguously identify two new kinematic components, both chemically rich with multiple molecular absorption lines. The "blue clump" has vLSR = -7.1 \pm 0.7 km/s and the "red clump" 1.4 \pm 0.5 km/s. While the blue and red clumps have similar temperatures and line widths, molecular species in the blue clump have higher column densities. They are both likely linked to pure rotational H2 emission also covered by this survey. This work provides evidence for the scenario that the blue and red clumps are distinct components unrelated to the classic components in the Orion BN/KL region. Comparison to spectroscopic surveys towards other infrared targets in the region show that the blue clump is clearly extended. We analyze, compare, and present in depth findings on the physical conditions of C2H2, 13CCH2, CH4, CS, H2O, HCN, H13CN, HNC, NH3, and SO2 absorption lines and an H2 emission line associated with the blue and red clumps. We also provide limited analysis of H2O and SiO molecular emission lines towards Orion IRc2 and the atomic forbidden transitions [FeII], [SI], [SIII], and [NeII]., Accepted to ApJ; 44 pages, 14 figures, 13 tables

A highly accurate CO

ConspectusWe review the Best Theory + Reliable High-Resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is twofold: (a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH

The semi-empirical molecular rovibrational IR line lists, such as ExoMol, TheoReTs, and Ames, combine the experimental accuracy and theoretical power to reach better than 0.1 cm−1 accuracy for line positions and better than 80–90% agreement for line intensities. The quality of these existing semi-empirical IR lists allows further improvements of intensity and line positions for those unobserved minor isotopologues. This paper presents our new BTRHE (Best Theory + Reliable High-resolution Experiment) strategy implementation. For line intensity, the isotopologue consistency and the patterns of mass dependence in the Ames-296 K SO2 and CO2 IR lists are quantitatively presented along the mass-inverse coordinates. The consistency and patterns are better than those in existing experimental data. The methodology proposed here can be used to identify inconsistencies, outliers, and mistakes in intensities, and help improve Effective Dipole Model (EDM) and molecular IR databases. We call for an experimental study on the 50006 and 60007 bands of CO2 628. For line position predictions, a simple approach combining the variational IR line lists with Effective Hamiltonian (EH) model may refine the effective rotational constants A0/B0/C0 and quartic centrifugal distortion constants of minor isotopologues. The prediction accuracy may be improved by two orders of magnitude, i.e. reaching 0–5 MHz prediction accuracy in the range of J

The quality of Ames-296K SO2 Infrared (IR) line list intensities is first validated by quantitative exploration of several dipole moment surfaces (DMSs) and partition sum convergence. The DMSs are computed with several of Dunning's correlation-consistent basis sets and their vibrational dependence are compared to the empirical model derived from Stark effect experiments reported by D. Patel, D. Margolese, and T.R. Dyke [J. Chem. Phys. 70, 2740 (1979)]. The effective dipole deviations from the DMS adopted in the Ames IR lists is 0.2–0.4% for vibrational states up to 3ν3. The vibrational dependence of the dipole moment is also in good agreement, except for nν1. Partition sum convergence at 296 K is confirmed by new calculations with rotational quantum number J up to 150 and upper state E’ up to 8000 cm−1. The isotopologue consistency of the Ames IR line lists is superior relative to the regular Effective Hamiltonian (EH) models and Effective Dipole Moment (EDM) models. The ν1 + ν2 and ν2 + ν3 intensity consistency check reveals the recently reported experimental intensities need significant improvement or re-analysis. After the accuracy, convergence, and isotopologue consistency have been confirmed, the theoretical Ames-296K intensities are combined with the experimental line positions or EH models that experimental spectroscopists published after 2009. Three high-resolution IR line sets are reported for the 32/33/34S16O2, 32S18O2 and 16O32S18O isotopologues: (1) the “New Lines Sets” include experimentally measured line positions; (2) the “Expanded Line Sets” include possible transitions among new rovibrational levels assigned in experiments and ground state (GS) levels predicted by reliable EH models; (3) the “Ames + MARVEL Sets” include possible transitions among all those levels reported in a recent MARVEL analysis. [Tobias et al, JQSRT 208, 152 (2018)]. Compared to the limited data in High-resolution TRANsmission molecular absorption database (HITRAN), these line sets have significantly improved the data coverage up to 4000 cm−1. Some missing bands can be traced to the unpublished experimental data. The isotopologue consistency of these line sets will help identify the uncertainties and defects in the experimental EH and EDM models. These line sets are good candidates for the next HITRAN update, if line shape parameters are available. The line sets can be downloaded from supplementary files or from the Ames Molecular Spectroscopic Database at http://huang.seti.org .

We present the first mid-infrared (MIR) detections of HNC and H13CN in the interstellar medium, and numerous, resolved HCN rovibrational transitions. Our observations span 12.8 to 22.9 micron towards the hot core Orion IRc2, obtained with the Echelon-Cross-Echelle Spectrograph aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). Exceptional, ~5 km/s, resolution distinguishes individual rovibrational transitions of the HNC and HCN P, Q, and R branches; and the H13CN R branch. This allows direct measurement of the species' excitation temperatures, column densities, and relative abundances. HNC and H13CN exhibit a local standard rest velocity of -7 km/s that may be associated with an outflow from nearby radio source I and an excitation temperature of about 100 K. We resolve two velocity components for HCN, the primary component also being at -7 km/s with temperature 165 K. The hottest component, which had never before been observed, is at 1 km/s with temperature 309 K. This is the closest component to the hot core's centre measured to date. The derived 12C/13C=13 is below expectation for Orion's Galactocentric distance, but the derived HCN/HNC=72 is expected for this extreme environment. Compared to previous sub-mm and mm observations, our SOFIA line survey of this region shows that the resolved MIR molecular transitions are probing a distinct physical component and isolating the chemistry closest to the hot core., 20 pages, 6 figures, submitted to ApJ, comments welcome

Polycyclic aromatic hydrocarbons (PAHs) have been shown to be ubiquitous in a large variety of distinct astrophysical environments and are therefore of great interest to astronomers. The majority of these findings are based on theoretically predicted spectra, which make use of scaled DFT harmonic frequencies for band positions and the double harmonic approximation for intensities. However, these approximations have been shown to fail at predicting high-resolution gas-phase infrared spectra accurately, especially in the CH-stretching region (2950–3150 cm−1, 3 μm). This is particularly worrying for the subset of hydrogenated or methylated PAHs to which astronomers attribute the observed non-aromatic features that appear in the CH-stretching region of spectral observations of the interstellar medium (ISM). In our previous work, we presented the anharmonic theoretical spectra of three linear PAHs and five non-linear PAHs, demonstrating the importance of including anharmonicities into theoretical calculations. In this work we extend these techniques to two methylated PAHs (9-methylanthracene, and 9,10-dimethylanthracene) and four hydrogenated PAHs (9,10-dihydroanthracene, 9,10-dihydrophenanthrene, 1,2,3,4-tetrahydronaphthalene, and 1,2,3,6,7,8-hexahydropyrene) in order to better understand the aliphatic IR features of substituted PAHs. The theoretical spectra are compared with the spectra obtained under matrix isolation low-temperature conditions for the full vibrational fundamental range and under high-resolution, low-temperature gas-phase conditions for the CH-stretching region. Excellent agreement is observed between the theoretical and high-resolution experimental spectra with a deviation of 0.00% ± 0.17%, and changes to the spectra of PAHs upon methylation and hydrogenated are tracked accurately and explained.

Total internal partition sums (TIPS) are reported for 166 isotopologues of 51 molecules important in planetary atmospheres. Molecules 1 to 50 are taken from the HITRAN2016 list, and, in some cases, additional isotopologues are considered for some of the molecules. Molecules 51–53 are C3H4, CH3, and CS2, respectively. TIPS are not reported for the O atom and CF4; thus, while there are 53 species in the list, data are reported for 51 molecules. The TIPS are determined by various methods from 1 K to a Tmax that ensures the TIPS reported have converged. These data are provided with HITRAN2016 and a new version of the TIPS code is available in both FORTRAN and python languages.

A new 626-based Ames-2 PES refinement and Ames-2016 line lists for 13 CO2 isotopologues are reported. A consistent σRMS = ±0.02 cm−1 is established for hundreds of isotopologue band origins using the Ames-2 PES. Ames-2016 line lists are computed at 296 K, 1000 K and 4000 K using the Ames-2 PES and the same DMS-N2 dipole surface used previously, with J up to 150, E′ up to 24,000 cm−1 or 18,000 cm−1 and appropriate intensity cutoffs. The lists are compared to the CDSD-296, CDSD-4000 databases, UCL line lists, and a few recent highly accurate CO2 intensity measurements. Both agreements and discrepancies are discussed. Compared to the old Ames CO2 lists, the Ames-2016 line lists have line position deviations reduced by 50% or more, which consequently leads to more reliable intensities. The line shape parameters in the Ames-2016 line lists are predicted using the newly assigned conventional vibrational polyad quantum numbers for rovibrational levels below 12,000 cm−1 so the quality of the line shape parameters is similar to that of CDSD or HITRAN. This study further proves that a semi-empirically refined PES (Ames-1 and Ames-2) coupled with a high quality ab initio DMS (DMS-N2 and UCL) may generate IR predictions with consistent accuracy and is thus helpful in the analysis of laboratory spectra and simulations of various isotopologues. The Ames-2016 lists based on DMS-N2 have reached the ∼1% intensity prediction accuracy level for the recent 626 30013-00001 and 20013-00001 bands, but further quantification and improvements require sub-percent or sub-half-percent accurate experimental intensities. The inter-isotopologue consistency of the intensity prediction accuracies should have reached better than 1–3% for regular bands not affected by resonances. Since the Effective Dipole Models (EDM) in CDSD and HITRAN have 1–20% or even larger uncertainties, we show that the Ames lists can provide better alternative IR data for many hard-to-determine isotopologue bands. Comparison at 4000 K suggests that the Ames-4000 K 12C16O2 line list is reliable and consistent within the current cutoffs of J ≤ 150 and E′ ≤ 24,000 cm−1, but intensity contributions involving higher energy levels should not be omitted and future computations need to be converged up to at least 32,000 cm−1 or higher. The remaining issues are discussed regarding the source of energy level discrepancies, intensity underestimations by ∼50% for some weak bands, etc. and also future work.

The simple aromatic hydrocarbon, cyclopropenylidene (c-C3H2), is a known, naturally-occurring molecule. The question remains as to whether its isoelectronic, cyclic, fellow aromatics of c-N3+, c-CNN, HCNN+, and c-CNC− are as well. Each of these are exciting objects for observation of Titan, and the rotational constants and vibrational frequencies produced here will allow for remote sensing of Titan's atmosphere or other astrophysical or terrestrial sources. None of these four aromatic species are vibrationally strong absorbers/emitters, but the two ions, HCNN+ and c-CNC−, have dipole moments of greater than 3 D and 1 D, respectively, making them good targets for rotational spectroscopic observation. Each of these molecules is shown here to exhibit its own, unique vibrational properties, but the general trends put the vibrational behavior for corresponding fundamental modes within close ranges of one another, even producing nearly the same heavy atom, symmetric stretching frequencies for HCNN+ and c-C3H2 at 1600 cm−1. The c-N3+ cation is confirmed to be fairly unstable and has almost no intensity in its ν2 fundamental. Hence, it will likely remain difficult to characterize experimentally.

Using the latest published, empirically refined potential energy surface (PES) Ames-1 and purely ab initio CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), we have computed infrared line lists for five symmetric isotopologues of sulfur dioxide: 32S16O2 (626), 33S16O2 (636), 34S16O2 (646), 36S16O2 (666), and 32S18O2 (828). The line lists are based on J = 0–80 rovibrational variational calculations with E′ ⩽ 8000 cm−1. The 34S16O2 and 33S16O2 line lists are compared to the experiment-based models in the HIgh-resolution TRANsmission molecular absorption database (HITRAN2012, http://www.cfa.harvard.edu/hitran/ ) and the Cologne Database for Molecular Spectroscopy, CDMS ( http://www.astro.uni-koeln.de/cdms/ ). The accuracy for computed 646 band origins is similar to what has been reported for the main isotopologue, i.e. 0.01–0.03 cm−1 for bands up to 5500 cm−1. For rovibrational transitions, the 646 line position and intensity deviation patterns are much simpler and more self-consistent than those of the main isotopologue 626. The discrepancies are mainly found for higher Ka/J transitions. 626 and 646 exhibit comparable line position and intensity agreement for lower Ka/J transitions. The line position deviations for the 636 purely rotational band are parallel to those of 626 and 646, while its line intensity deviations do not show branching patterns as we found in the 626 and 646 cases. Predictions for the other minor isotopologues are expected to exhibit similar accuracy. These line lists are accurate enough to provide alternatives for missing bands of 626 and the minor isotopologues. It may significantly facilitate the laboratory spectroscopic measurement and analysis, as well as to identify these isotopologues in various astrophysical environments.

Reliable infrared (IR) line lists are reported for the 13 isotopologues of carbon dioxide in HITRAN notation: 626, 636, 628, 627, 828, 727, 827, 638, 637, 737, 838, 738, and 646. Three IR lists are available for each istotopologue: a complete list at 296 K, a reduced-size list at 296 K, plus a reduced-size list at 1000 K. They are denoted Ames-296K, Ames-296K.reduced and Ames-1000K.reduced. With J up to 150, and energy up to 18,000 cm−1 above the zero point energy, these lists are expected to cover the temperature range up to 1500 K. Line shape parameters including temperature dependence are calculated and reported for four temperature ranges: Mars, Earth, Venus, and Hotter (700–2000 K). Comparisons are made against the available transition data in the HITRAN2012 models. Line position accuracy for most transitions up to 10,000–13,000 cm−1 is better than 0.03–0.05 cm−1. Computed transition intensities agree well with most HITRAN data but there exist suspicious exceptions for isotopologues. These line lists will expedite CO2 IR experimental data analysis and provide the scientific community with trustworthy alternatives for unknown IR bands. These line lists may be combined with existing experimental databases to facilitate the analysis of future laboratory experiments or astronomical observations.

We present high-spectral resolution observations from 12.96 - 13.33 microns towards Orion IRc2 using the mid-infrared spectrograph, EXES, on SOFIA. These observations probe the physical and chemical conditions of the Orion hot core, which is sampled by a bright, compact, mid-infrared background continuum source in the region, IRc2. All ten of the rovibrational C2H2 transitions expected in our spectral coverage, are detected with high S/N, yielding continuous coverage of the R-branch lines from J=9-8 to J=18-17, including both ortho and para species. Eight of these rovibrational transitions are newly reported detections. The isotopologue, 13CCH2, is clearly detected with high signal-to-noise. This enabled a direct measurement of the 12C/13C isotopic ratio for the Orion hot core of 14 +/- 1 and an estimated maximum value of 21. We also detected several HCN rovibrational lines. The ortho and para C2H2 ladders are clearly separate and tracing two different temperatures, 226 K and 164 K, respectively, with a non-equilibrium ortho to para ratio (OPR) of 1.7 +\- 0.1. Additionally, the ortho and para V_LSR values differ by about 1.8 +/- 0.2 km/s, while, the mean line widths differ by 0.7 +/- 0.2 km/s, suggesting that these species are not uniformly mixed along the line of sight to IRc2. We propose that the abnormally low C2H2 OPR could be a remnant from an earlier, colder phase, before the density enhancement (now the hot core) was impacted by shocks generated from an explosive event 500 yrs ago., accepted for publication in The Astrophysical Journal

New semi-empirical Infrared (IR) line lists for 12C16O2, Ames-296 K and Ames-1000 K, have been computed using a newly updated ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface (denoted DMS-N2) and an empirically refined potential energy surface (Ames-1). J=0–150 rovibrational levels are computed up to 30,000 cm−1, and related transitions are cut off at 1E−42 cm molecule−1 (296 K) and 1E−36 cm molecule−1 (1000 K). These are the first line lists available to cover reliably the energy region as high as ~20,000 cm−1. Recent experimental data at 1.1 μm has confirmed the predicted intensities for the 50013-00001 and 50014-00001 band transitions have better than 90% agreement. Comparisons are made against the Wattson 750 K line list and HITEMP/HITRAN at 300 K, 500 K, 725 K, 1000 K, 1500 K, 2000 K and 3000 K. The temperature dependence and accuracy of the new Ames-296 K/1000 K line lists are investigated and we claim both line lists are capable of providing reliable opacities up to 18,000–23,000 cm−1, while the highest applicable wavenumber range drops as T rises. We suggest caution is used for T>1000 K simulations. Comparison to recent experiments at 1000 K, 1550 K and 1773 K shows that the Ames-1000 K line list and HITEMP perform similarly in the 3200–3800 cm−1 and 4600–5200 cm−1 ranges. In the 2000–2100 cm−1 range, Ames-1000 K yields better agreement relative to experiment. Existing problems and possible future solutions in the new Ames-296 K/1000 K line lists for line positions and intensities are also discussed.