Sazuna/UAT_keywords huggingface dataset:
Bits Missing: Finding Exotic Pulsars Using bfloat16 on NVIDIA GPUs. The Fourier domain acceleration search (FDAS) is an effective technique for detecting faint binary pulsars in large radio astronomy data sets. This paper quantifies the sensitivity impact of reducing numerical precision in the graphics processing unit (GPU)-accelerated FDAS pipeline of the AstroAccelerate (AA) software package. The prior implementation used IEEE-754 single-precision in the entire binary pulsar detection pipeline, spending a large fraction of the runtime computing GPU-accelerated fast Fourier transforms. AA has been modified to use bfloat16 (and IEEE-754 double-precision to provide a "gold standard" comparison) within the Fourier domain convolution section of the FDAS routine. Approximately 20,000 synthetic pulsar filterbank files representing binary pulsars were generated using SIGPROC with a range of physical parameters. They have been processed using bfloat16, single-precision, and double-precision convolutions. All bfloat16 peaks are within 3% of the predicted signal-to-noise ratio of their corresponding single-precision peaks. Of 14,971 "bright" single-precision fundamental peaks above a power of 44.982 (our experimentally measured highest noise value), 14,602 (97.53%) have a peak in the same acceleration and frequency bin in the bfloat16 output plane, while in the remaining 369 the nearest peak is located in the adjacent acceleration bin. There is no bin drift measured between the single- and double-precision results. The bfloat16 version of FDAS achieves a speedup of approximately 1.6× compared to single-precision. A comparison between AA and the PRESTO software package is presented using observations collected with the GMRT of PSR J1544+4937, a 2.16 ms black widow pulsar in a 2.8 hr compact orbit.

Keywords: Computational astronomy, Astronomy software, Astronomy data analysis, Computational methods, Binary pulsars, GPU computing, Astronomical methods
Study of Variability in Long-term Multiwavelength Optical Lightcurves of Blazar AO 0235+164. We present a long-term and intraday variability study of optical multiwaveband (U BV RI) data from the blazar AO 0235+164 collected by various telescopes for ~44 yr (1975-2019). The blazar was found to be significantly variable over the years in all wave bands with a variation of about 6 mag between its low and active states. The variations in the different wave bands are highly correlated without any time lag. We did not observe any significant trend in color variation with time, but we observed a bluer-when-brighter trend between the B - I color index and the R magnitude. Optical BV R-band spectral energy distributions always show a convex shape. Significant intraday variability was frequently seen in the quasi-simultaneous observations of AO 0235+164 made on 22 nights in the R and V bands by the CASLEO and CAHA telescopes during 1999-2019. We also estimated a central supermassive black hole mass of 7.9 × 107 M ⊙ by analyzing the broad Mg II emission line in AO 0235+164's spectrum. We briefly explore the probable physical scenarios responsible for the observed variability.

Keywords: Active galaxies, BL Lacertae objects, Galaxies, Quasars
Modules for Experiments in Stellar Astrophysics (MESA): Time-dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure. We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthen MESA's implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates to MESA's software infrastructure that enhance source code development and community engagement.

Keywords: Stellar physics, Computational methods, Stellar evolutionary models, Stellar evolution
Sample and Statistical Analysis of the Near-Earth Object Wide-field Infrared Survey Explorer Variability of the 6.7 GHz Methanol Maser Sources. We present a statistical study of infrared variability using the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE) database for a sample consisting of 1085 high-mass young stellar objects (YSOs) related to 6.7 GHz methanol masers. A total of 383 maser sources were identified as NEOWISE variables and classified in two variability behavior classifications: 204 secular (linear, curved, and periodic) and 179 stochastic (burst, drop, and irregular) variables. Statistical analysis of the properties of these variables (e.g., the dust temperature, bolometric luminosity, hydrogen column density, W4 luminosity, and W1-W2 color) has revealed a potential evolutionary sequence among different light-curve types of variables. There is a possible general evolutionary (from less to more evolved) trend between the three variable types from secular to stochastic to nonvariable. For the specific classifications, the evolutionary trend for secular variables is linear to sin to sin+linear, and for stochastic variables it is burst to irregular to drop. These sequences may reflect the evolution of the envelope or accretion disk of high-mass YSOs, from large to small radii due to gravitational collapse. Although no significant variability correlation was established between the 6.7 GHz methanol maser and the W2-band emission based on the data collected so far, a number of candidates were found for further investigating the accretion burst events via future variability monitoring programs of both mid-infrared emission and masers.

Keywords: Massive stars, Young stellar objects, Infrared photometry, Interstellar masers
GECAM Localization of High-energy Transients and the Systematic Error. The Gravitational Wave High-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a pair of microsatellites (i.e., GECAM-A and GECAM-B) dedicated to monitoring gamma-ray transients including the high-energy electromagnetic counterparts of gravitational waves, such as gamma-ray bursts, soft gamma-ray repeaters, solar flares, and terrestrial gamma-ray flashes. Since launch in 2020 December, GECAM-B has detected hundreds of astronomical and terrestrial events. For these bursts, localization is the key for burst identification and classification as well as follow-up observations in multiple wavelengths. Here, we propose a Bayesian localization method with Poisson data with Gaussian background profile likelihood to localize GECAM bursts based on the distribution of burst counts in detectors with different orientations. We demonstrate that this method can work well for all kinds of bursts, especially extremely short ones. In addition, we propose a new method to estimate the systematic error of localization based on a confidence level test, which can overcome some problems of the existing method in the literature. We validate this method by Monte Carlo simulations, and then apply it to a burst sample with accurate location and find that the mean value of the systematic error of GECAM-B localization is ~2.°5. By considering this systematic error, we can obtain a reliable localization probability map for GECAM bursts. Our methods can be applied to other gamma-ray monitors.

Keywords: Observational astronomy, Gamma-ray astronomy, Astronomical methods, Time domain astronomy
The Nature of Blue Stars with Mid-infrared Excesses in the Large Magellanic Cloud. We present low-resolution optical spectra and classifications of 92 blue objects with mid-infrared excesses in the Large Magellanic Cloud. The majority of these objects were selected with the criteria of U - B &lt; 0 and V &lt; 17 from the potential young stellar object (YSO) candidates in Gruendl &amp; Chu (GC09), which were identified based on Spitzer Infrared Array Camera and Multiband Imaging Photometer for Spitzer observations in conjunction with optical photometry from the Magellanic Clouds Photometric Survey. Many of the sample objects have ambiguous classifications. We examined the properties of these 92 objects using low-resolution optical spectra obtained with the SOAR 4.1 m Telescope at Cerro Pachon and the Blanco 4 m Telescope at Cerro Tololo Inter-American Observatory, supplemented by available photometric and imaging observations. We estimated the spectral types, temperatures, and luminosities of these objects from the optical to near-IR spectral energy distributions based on the photometric data, and further examined stellar absorption line features in the optical spectra to verify the spectral types. The interstellar/circumstellar environments, assessed from nebular line imaging observations and nebular lines detected in the stellar spectra, further helped constrain the nature of stars. Among these 92 objects, we confirm 42 stars as YSOs, and the remaining 50 objects as protoplanetary nebulae, post-AGB/RGB stars, blue evolved massive stars, stars with dust in vicinity, or uncertain classifications. Our results show that the photometric criteria in GC09 are generally effective in the initial selection of YSO candidates, and the low-resolution spectroscopy combined with environment assessment can be useful to better constrain the classifications and ameliorate most ambiguities.

Keywords: Stellar classification, Young stellar objects, Large Magellanic Cloud
Modeling the Solar Corona with an Implicit High-order Reconstructed Discontinuous Galerkin Scheme. The present study aims to apply an implicit high-order reconstructed discontinuous Galerkin (DG) scheme (rDG(P 1 P 2)) to simulate the steady-state solar corona. In this scheme, a piecewise quadratic polynomial solution, P 2, is obtained from the underlying piecewise linear DG solution, P 1, by least-squares reconstruction with a weighted essentially nonoscillatory limiter. The reconstructed quadratic polynomial solution is then used for the computation of the fluxes and source terms. In addition, an implicit time integration method with large time steps is considered in this work. The resulting large linear algebraic system of equations from the implicit discretization is solved by the cellwise relaxation implicit scheme which can make full use of the compactness of the DG scheme. The code of the implicit high-order rDG scheme is developed in Fortran language with message passing interface parallelization in Cartesian coordinates. To validate this code, we first test a problem with an exact solution, which confirms the expected third-order accuracy. Then we simulate the solar corona for Carrington rotations 2167, 2183, and 2210, and compare the modeled results with observations. We find that the numerical results basically reproduce the large-scale observed structures of the solar corona, such as coronal holes, helmet streamers, pseudostreamers, and high- and low-speed streams, which demonstrates the capability of the developed scheme.

Keywords: Solar magnetic fields, Solar corona, Magnetohydrodynamical simulations
Identifying 46 New Open Cluster Candidates in Gaia EDR3 Using a Hybrid pyUPMASK and Random Forest Method. Open clusters (OCs) are regarded as tracers to understand stellar evolution theory and validate stellar models. In this study, we presented a robust approach to identifying OCs. A hybrid method consisting of pyUPMASK and the random forest (RF) algorithm is first used to remove field stars and determine more reliable members. An identification model based on the RF algorithm built based on 3714 OC samples from Gaia DR2 and EDR3 is then applied to identify OC candidates. The OC candidates are obtained after isochrone fitting, advanced stellar population synthesis model fitting, and visual inspection. Using the proposed approach, we revisited 868 candidates and preliminarily clustered them by the friends-of-friends algorithm in Gaia EDR3. Excluding OCs that have already been reported, we focused on the remaining 300 unknown candidates. From high to low fitting quality, these unrevealed candidates were further classified into Class A (59), Class B (21), and Class C (220). As a result, 46 new reliable OC candidates among Classes A and B are identified after visual inspection.

Keywords: Open star clusters
CLASS Survey Description: Coronal-line Needles in the SDSS Haystack. Coronal lines are a powerful, yet poorly understood, tool to identify and characterize active galactic nuclei. There have been few large-scale surveys of coronal lines in the general galaxy population in the literature so far. Using a novel preselection technique with a flux-to-rms ratio ${ \mathcal F }$ , followed by Markov Chain Monte Carlo fitting, we searched for the full suite of 20 coronal lines in the optical spectra of almost 1 million galaxies from the Sloan Digital Sky Survey Data Release 8. We present a catalog of the emission-line parameters for the resulting 258 galaxies with detections. The Coronal Line Activity Spectroscopic Survey includes line properties, host-galaxy properties, and selection criteria for all galaxies in which at least one line is detected. This comprehensive study reveals that a significant fraction of coronal-line activity is missed in past surveys based on a more limited set of coronal lines; ~60% of our sample do not display the more widely surveyed [Fe X] λ6374. In addition, we discover a strong correlation between coronal-line and Wide-field Infrared Survey Explorer W2 luminosities, suggesting that the mid-infrared flux can be used to predict coronal-line fluxes. For each line we also provide a confidence level that the line is present, generated by a novel neural network, trained on fully simulated data. We find that after training the network to detect individual lines using 100,000 simulated spectra, we achieve an overall true-positive rate of 75.49% and a false-positive rate of only 3.96%.

Keywords: Starburst galaxies, Active galaxies, Dwarf galaxies, Solar coronal lines, Infrared galaxies, Galaxy evolution, Interstellar medium
Prospects for Ray-tracing Light Intensity and Polarization in Models of Accreting Compact Objects Using a GPU. The Event Horizon Telescope (EHT) has recently released high-resolution images of accretion flows onto two supermassive black holes. Our physical understanding of these images depends on the accuracy and precision of numerical models of plasma and radiation around compact objects. The goal of this work is to speed up radiative-transfer simulations used to create mock images of black holes for comparison with the EHT observations. A ray-tracing code for general relativistic and fully polarized radiative transfer through plasma in strong gravity is ported onto a graphics processing unit (GPU). We describe our GPU implementation and carry out speedup tests using models of optically thin advection-dominated accretion flow onto a black hole realized semianalytically and in 3D general relativistic magnetohydrodynamic simulations, low and very high image pixel resolutions, and two different sets of CPU+GPUs. We show that a GPU with high double precision computing capability can significantly reduce the image production computational time, with a speedup factor of up to approximately 1200. The significant speedup facilitates, e.g., dynamic model fitting to the EHT data, including polarimetric data. The method extension may enable studies of emission from plasma with nonthermal particle distribution functions for which accurate approximate synchrotron emissivities are not available. The significant speedup reduces the carbon footprint of the generation of the EHT image libraries by at least an order of magnitude.

Keywords: Black holes, Astronomy software, Computational methods, General relativity, Radiative transfer, GPU computing
Neural Stellar Population Synthesis Emulator for the DESI PROVABGS. The Probabilistic Value-Added Bright Galaxy Survey (PROVABGS) catalog will provide the posterior distributions of physical properties of &gt;10 million DESI Bright Galaxy Survey galaxies. Each posterior distribution will be inferred from joint Bayesian modeling of observed photometry and spectroscopy using Markov Chain Monte Carlo sampling and the Hahn et al. stellar population synthesis (SPS) model. To make this computationally feasible, PROVABGS will use a neural emulator for the SPS model to accelerate the posterior inference. In this work, we present how we construct the emulator using the Alsing et al. approach and verify that it can be used to accurately infer galaxy properties. We confirm that the emulator is in excellent agreement with the original SPS model with ≪1% error and is 100× faster. In addition, we demonstrate that the posteriors of galaxy properties derived using the emulator are also in excellent agreement with those inferred using the original model. The neural emulator presented in this work is essential in bypassing the computational challenge posed in constructing the PROVABGS catalog. Furthermore, it demonstrates the advantages of emulation for scaling sophisticated analyses to millions of galaxies.

Keywords: Astronomy data modeling, Spectral energy distribution, Galaxies, Neural networks
Determining the Ionization Rates of Interstellar Neutral Species Using Direct-sampling Observations of Their Direct and Indirect Beams. A good understanding of the ionization rates of neutral species in the heliosphere is important for studies of the heliosphere and planetary atmospheres. So far, the intensities of the ionization reactions have been studied based on observations of the contributing phenomena, such as the solar spectral flux in the EUV band and the flux of the solar wind protons, alpha particles, and electrons. The results strongly depend on absolute calibration of these measurements, which, especially for the EUV measurements, is challenging. Here, we propose a novel method of determining the ionization rate of neutral species based on direct sampling of interstellar neutral gas from two locations in space distant to each other. In particular, we suggest performing observations from the vicinity of Earth's orbit and using ratios of fluxes of interstellar neutral He for the direct and indirect orbits of interstellar atoms. We identify the most favorable conditions and observations geometries, suitable for implementation on the forthcoming NASA Interstellar Mapping and Acceleration Probe mission.

Keywords: Solar extreme ultraviolet emission, Interstellar atomic gas, Heliosphere, Solar wind, Interstellar medium, Solar activity, Interstellar medium wind
The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Quasar Survey: Quasar Properties from Data Releases 6 to 9. We report the fourth installment in the series of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) quasar survey, which includes quasars observed between 2017 September and 2021 June. There are in total 13,066 quasars reliably identified, of which 6685 are newly discovered that are not reported in the Sloan Digital Sky Survey (SDSS) DR14 quasar catalog or Million Quasars catalog. Because LAMOST does not provide accurate absolute flux calibration, we recalibrate the spectra with the SDSS/Pan-STARRS1 multiband photometric data. The emission-line properties of Hα, Hβ, Mg II, and C IV and the continuum luminosities are measured by fitting the recalibrated spectra. We also estimate the single-epoch virial black hole masses (M BH) using the derived emission-line and continuum parameters. This is the first time that the emission-line and continuum fluxes were estimated based on LAMOST recalibrated quasar spectra. The catalog and spectra for these quasars are available online. After the 9 yr LAMOST quasar survey, there are in total 56,175 identified quasars, of which 24,127 are newly discovered. The LAMOST quasar survey not only discovers a great number of new quasars but also provides a database for investigating the spectral variability of the quasars observed by both LAMOST and SDSS and finding rare quasars, including changing-look quasars and broad absorption line quasars.

Keywords: Surveys, Catalogs, Quasars
Extended MCDHF Calculations of Energy Levels and Transition Data for N I. Accurate and extensive atomic data are essential for spectroscopic analyses of stellar atmospheres and other astronomical objects. We present energy levels, lifetimes, and transition probabilities for neutral nitrogen, the sixth most abundant element in the cosmos. The calculations employ the fully relativistic multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods, and span the 103 lowest states up to and including 2s22p25s. Our theoretical energies are in excellent agreement with the experimental data, with an average relative difference of 0.07%. In addition, our transition probabilities are in good agreement with available experimental and theoretical data. We further verify the agreement of our data with experimental results via a reanalysis of the solar nitrogen abundance, with the results from the Babushkin and Coulomb gauges consistent to 2% or 0.01 dex. We estimated the uncertainties of the computed transition data based on a statistical analysis of the differences between the transition rates in the Babushkin and Coulomb gauges. Out of the 1701 computed electric dipole transitions in this work, 83 (536) are associated with uncertainties smaller than 5% (10%).

Keywords: Solar abundances, Atomic data
Young Stellar Objects, Accretion Disks, and Their Variability with Rubin Observatory LSST. Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, can be explored by taking advantage of the multiband filters option available in Rubin LSST, in particular the u, g, r, i filters that enable us to discriminate between photospheric stellar properties and accretion signatures. The homogeneity and depth of sky coverage that will be achieved with LSST will provide us with a unique opportunity to characterize the time evolution of disk accretion as a function of age and varying environmental conditions (e.g., field crowdedness, massive neighbors, metallicity) by targeting different star-forming regions. In this contribution to the Rubin LSST Survey Strategy Optimization Focus Issue, we discuss how implementing a dense observing cadence to explore short-term variability in YSOs represents a key complementary effort to the Wide-Fast-Deep observing mode that will be used to survey the sky over the full duration of the main survey (≈10 yr). The combination of these two modes will be vital to investigate the connection between the inner-disk dynamics and longer-term eruptive variability behaviors, such as those observed on EX Lupi-type objects.

Keywords: Light curves, Young star clusters, Classical T Tauri stars, T Tauri stars, Young stellar objects, Low mass stars, Photometry, Circumstellar disks, Stellar accretion disks, Variable stars