Table of contents

Recent progress in high-dispersion spectroscopy has revealed the presence of vaporized heavy metals and ions in the atmosphere of hot Jupiters whose dayside temperature is larger than 2000 K, categorized as ultrahot Jupiters (UHJs). Using the archival data of high-resolution transmission spectroscopy obtained with the Subaru telescope, we searched for neutral metals in HD 149026b, a hot Jupiter cooler than UHJs. By removing stellar and telluric absorption and using a cross-correlation technique, we report a tentative detection of neutral titanium with 4.4σ and a marginal signal of neutral iron with 2.8σ in the atmosphere. This is the first detection of neutral titanium in an exoplanetary atmosphere. In this temperature range, titanium tends to form titanium oxide (TiO). The fact that we did not detect any signal from TiO suggests that the C/O ratio in the atmosphere is higher than the solar value. The detection of metals in the atmosphere of hot Jupiters cooler than UHJs will be useful for understanding the atmospheric structure and formation history of hot Jupiters.

We conduct X-ray spectral fits on 184 likely counterparts to Fermi-LAT 3FGL unassociated sources. Characterization and classification of these sources allows for more complete population studies of the high-energy sky. Most of these X-ray spectra are well fit by an absorbed power-law model, as expected for a population dominated by blazars and pulsars. A small subset of seven X-ray sources have spectra unlike the power law expected from a blazar or pulsar and may be linked to coincident stars or background emission. We develop a multiwavelength machine learning classifier to categorize unassociated sources into pulsars and blazars using gamma-ray and X-ray observations. Training a random forest (RF) procedure with known pulsars and blazars, we achieve a cross-validated classification accuracy of 98.6%. Applying the RF routine to the unassociated sources returned 126 likely blazar candidates (defined as P_bzr ≥ 90%) and five likely pulsar candidates (P_bzr ≤ 10%). Our new X-ray spectral analysis does not drastically alter the RF classifications of these sources compared to previous works, but it builds a more robust classification scheme and highlights the importance of X-ray spectral fitting. Our procedure can be further expanded with UV, visual, or radio spectral parameters or by measuring flux variability.

We present results from Speckle inteferometric observations of 15 visual binaries and one double-line spectroscopic binary, carried out with the HRCam Speckle camera of the SOAR 4.1 m telescope. These systems were observed as a part of an on-going survey to characterize the binary population in the solar vicinity, out to a distance of 250 pc. We obtained orbital elements and mass sums for our sample of visual binaries. The orbits were computed using a Markov Chain Monte Carlo algorithm that delivers maximum likelihood estimates of the parameters, as well as posterior probability density functions that allow us to evaluate their uncertainty. Their periods cover a range from 5 yr to more than 500 yr; and their spectral types go from early A to mid M, implying total system masses from slightly more than 4M_⊙ down to 0.2M_⊙. They are located at distances between approximately 12 and 200 pc, mostly at low Galactic latitude. For the double-line spectroscopic binary YSC8, we present the first combined astrometric/radial-velocity orbit resulting from a self-consistent fit, leading to individual component masses of 0.897 ± 0.027 M_⊙ and 0.857 ± 0.026 M_⊙; and an orbital parallax of 26.61 ± 0.29 mas, which compares very well with the Gaia DR2 trigonometric parallax (26.55 ± 0.27 mas). In combination with published photometry and trigonometric parallaxes, we place our objects on an H-R diagram and discuss their evolutionary status. We also present a thorough analysis of the precision and consistency of the photometry available for them.

Two large, ∼500 M_☉ elliptical ring structures have been identified in the high-mass star-forming region NGC 7538. The origin of these ring structures is unknown, so we investigate the possibility that a runaway O- or B-type star may have originated in or passed through the region and created either or both of the ring structures via stellar wind. In testing this hypothesis, we identify one candidate star, BD +61 2408, that may have formed the northern ring. It is a B3e star with a mass of ∼8 M_☉ and a surface temperature of ∼20,000 K. Its position, motion, timescale, and spectral type are all consistent with the star being a candidate for having formed one of the ring structures in NGC 7538.

We present an analysis of 1524 spectra of Vega spanning 10 yr, in which we search for periodic radial-velocity variations. A signal with a periodicity of 0.676 day and a semi-amplitude of ∼10 m s⁻¹ is consistent with the rotation period measured over much shorter time spans by previous spectroscopic and spectropolarimetric studies, confirming the presence of surface features on this A0 star. The activity signal appears to evolve on long timescales, which may indicate the presence of failed fossil magnetic fields on Vega. TESS data reveal Vega's photometric rotational modulation for the first time, with a total amplitude of only 10 ppm. A comparison of the spectroscopic and photometric amplitudes suggests that the surface features may be dominated by bright plages rather than dark spots. For the shortest orbital periods, transit and radial-velocity injection recovery tests exclude the presence of transiting planets larger than 2 R_⊕ and most non-transiting giant planets. At long periods, we combine our radial velocities with direct imaging from the literature to produce detection limits for Vegan planets and brown dwarfs out to distances of 15 au. Finally, we detect a candidate radial-velocity signal with a period of 2.43 days and a semi-amplitude of 6 m s⁻¹. If caused by an orbiting companion, its minimum mass would be ∼20 M_⊕; because of Vega's pole-on orientation, this would correspond to a Jovian planet if the orbit is aligned with the stellar spin. We discuss the prospects for confirmation of this candidate planet.

We present Galaxy Line Emission & Absorption Modeling (gleam), a Python tool for fitting Gaussian models to emission and absorption lines in large samples of 1D extragalactic spectra. gleam is tailored to work well in batch mode without much human interaction. With gleam, users can uniformly process a variety of spectra, including galaxies and active galactic nuclei, in a wide range of instrument setups and signal-to-noise regimes. gleam also takes advantage of multiprocessing capabilities to process spectra in parallel. With the goal of enabling reproducible workflows for its users, gleam employs a small number of input files, including a central, user-friendly configuration in which fitting constraints can be defined for groups of spectra and overrides can be specified for edge cases. For each spectrum, gleam produces a table containing measurements and error bars for the detected spectral lines and continuum and upper limits for nondetections. For visual inspection and publishing, gleam can also produce plots of the data with fitted lines overlaid. In the present paper, we describe gleam's main features, the necessary inputs, expected outputs, and some example applications, including thorough tests on a large sample of optical/infrared multi-object spectroscopic observations and integral field spectroscopic data. gleam is developed as an open-source project hosted at https://github.com/multiwavelength/gleam and welcomes community contributions.

High-dispersion spectra in the Li λ6708 region have been obtained and analyzed in the old, metal-deficient cluster NGC 2243. From Hydra spectra for 29 astrometric and radial velocity members, we derive rotational velocities, as well as [Fe/H], [Ca/H], [Si/H], and [Ni/H] based on 17, 1, 1, and 3 lines, respectively. Using ROBOSPECT, an automatic equivalent width measurement program, we derive [Fe/H] = −0.54 ± 0.11 (MAD), for an internal precision for the cluster [Fe/H] below 0.03 dex. Given the more restricted line set, comparable values for [Ca/H], [Si/H], and [Ni/H] are −0.48 ± 0.19, −0.44 ± 0.11, and −0.61 ± 0.06, respectively. With E(B − V) = 0.055, appropriate isochrones imply (m − M) = 13.2 ± 0.1 and an age of 3.6 ± 0.2 Gyr. Using available VLT spectra and published Li abundances, we construct an Li sample of over 100 stars extending from the tip of the giant branch to 0.5 mag below the Li dip. The Li dip is well populated and, when combined with results for NGC 6819 and Hyades/Praesepe, implies a mass/metallicity slope of 0.4 M_⊙/dex for the high-mass edge of the Li dip. The A(Li) distribution among giants reflects the degree of Li variation among the turnoff stars above the Li dip, itself a function of stellar mass and metallicity and strongly anticorrelated with a v_rot distribution that dramatically narrows with age. Potential implications of these patterns for the interpretation of Li among dwarf and giant field populations, especially selection biases tied to age and metallicity, are discussed.

We introduce a new binary detection technique, Binary INformation from Open Clusters using SEDs (binocs), which we show is able to determine reliable stellar multiplicity and masses over a much larger mass range than current approaches. This new technique determines accurate component masses of binary and single systems of the open clusters' main sequence by comparing observed magnitudes from multiple photometric filters to synthetic star spectral energy distributions (SEDs), allowing us to systematically probe the binary population for low-mass stars in clusters for eight well-studied open clusters. We provide new deep, infrared photometric catalogs (1.2–8.0 μm) for the key open clusters NGC 1960 (M36), NGC 2099 (M37), NGC 2420, and NGC 2682 (M67), using observations from NOAO/NEWFIRM and Spitzer/IRAC. Using these deep multiwavelength catalogs, the binocs method is applied to these clusters to determine accurate component masses for unresolved cluster binaries. We explore binary fractions as a function of cluster age, Galactic location, and metallicity.

TOI-216 hosts a pair of warm, large exoplanets discovered by the TESS mission. These planets were found to be in or near the 2:1 resonance, and both of them exhibit transit timing variations (TTVs). Precise characterization of the planets' masses and radii, orbital properties, and resonant behavior can test theories for the origins of planets orbiting close to their stars. Previous characterization of the system using the first six sectors of TESS data suffered from a degeneracy between planet mass and orbital eccentricity. Radial-velocity measurements using HARPS, FEROS, and the Planet Finder Spectrograph break that degeneracy, and an expanded TTV baseline from TESS and an ongoing ground-based transit observing campaign increase the precision of the mass and eccentricity measurements. We determine that TOI-216c is a warm Jupiter, TOI-216b is an eccentric warm Neptune, and that they librate in 2:1 resonance with a moderate libration amplitude of ${60}_{-2}^{+2}$ deg, a small but significant free eccentricity of ${0.0222}_{-0.0003}^{+0.0005}$ for TOI-216b, and a small but significant mutual inclination of 12–39 (95% confidence interval). The libration amplitude, free eccentricity, and mutual inclination imply a disturbance of TOI-216b before or after resonance capture, perhaps by an undetected third planet.

We report the discovery of a sextuply eclipsing sextuple star system from TESS data, TIC 168789840, also known as TYC 7037-89-1, the first known sextuple system consisting of three eclipsing binaries. The target was observed in Sectors 4 and 5 during Cycle 1, with lightcurves extracted from TESS Full Frame Image data. It was also previously observed by the WASP survey and ASAS-SN. The system consists of three gravitationally bound eclipsing binaries in a hierarchical structure of an inner quadruple system with an outer binary subsystem. Follow-up observations from several different observatories were conducted as a means of determining additional parameters. The system was resolved by speckle interferometry with a 042 separation between the inner quadruple and outer binary, inferring an estimated outer period of ∼2 kyr. It was determined that the fainter of the two resolved components is an 8.217 day eclipsing binary, which orbits the inner quadruple that contains two eclipsing binaries with periods of 1.570 days and 1.306 days. Markov Chain Monte Carlo (MCMC) analysis of the stellar parameters has shown that the three binaries of TIC 168789840 are "triplets," as each binary is quite similar to the others in terms of mass, radius, and T_eff. As a consequence of its rare composition, structure, and orientation, this object can provide important new insight into the formation, dynamics, and evolution of multiple star systems. Future observations could reveal if the intermediate and outer orbital planes are all aligned with the planes of the three inner eclipsing binaries.

We present the results of commissioning observations for a new digital beam-forming back end for the Focal plane L-band Array for the Robert C. Byrd Green Bank Telescope (FLAG), a cryogenically cooled Phased Array Feed (PAF) with the lowest measured T_sys/η of any PAF outfitted on a radio telescope to date. We describe the custom software used to apply beam-forming weights to the raw element covariances to create research-quality spectral-line images for the new fine-channel mode, study the stability of the beam weights over time, characterize FLAG's sensitivity over a frequency range of 150 MHz, and compare the measured noise properties and observed distribution of neutral hydrogen emission from several extragalactic and Galactic sources with data obtained with the current single-pixel L-band receiver. These commissioning runs establish FLAG as the preeminent PAF receiver currently available for spectral-line observations on the world's major radio telescopes.

We present high-resolution speckle interferometric imaging observations of TESS exoplanet host stars using the NN-EXPLORE Exoplanet and Stellar Speckle Imager instrument at the 3.5 m WIYN telescope. Eight TESS objects of interest that were originally discovered by Kepler were previously observed using the Differential Speckle Survey Instrument. Speckle observations of 186 TESS stars were carried out, and 45 (24%) likely bound companions were detected. This is approximately the number of companions we would expect to observe given the established 46% binarity rate in exoplanet host stars. For the detected binaries, the distribution of stellar mass ratio is consistent with that of the standard Raghavan distribution and may show a decrease in high-q systems as the binary separation increases. The distribution of binary orbital periods, however, is not consistent with the standard Ragahavan model, and our observations support the premise that exoplanet-hosting stars with binary companions have, in general, wider orbital separations than field binaries. We find that exoplanet-hosting binary star systems show a distribution peaking near 100 au, higher than the 40–50 au peak that is observed for field binaries. This fact led to earlier suggestions that planet formation is suppressed in close binaries.

We present a detailed photometric and spectroscopic analysis of DD CMa, based on published survey photometry and new spectroscopic data. We find an improved orbital period of P_o = 2.0084530(6) days. Our spectra reveal Hβ and Hα absorptions with weak emission shoulders, and we also find a color excess in the Wide-field Infrared Survey Explorer multiband photometry, interpreted as signatures of circumstellar matter. We model the V-band orbital light curve derived from the ASAS and ASAS-SN surveys, assuming a semidetached configuration and using the mass ratio and temperature of the hotter star derived from our spectroscopic analysis. Our model indicates that the system consists of a B2.5 dwarf and a B9 giant of radii 3.2 and 3.7 R_⊙, respectively, orbiting in a circular orbit of radius 6.75 R_⊙. We also found M_c = 1.7 ± 0.1 M_⊙, T_c = 11,350 ± 100 K, and M_h = 6.4 ± 0.1 M_⊙, T_h = 20,000 ± 500 K, for the cooler and hotter star, respectively. We find broad single emission peaks in Hα and Hβ after subtracting the synthetic stellar spectra. Our results are consistent with mass exchange between the stars and suggest the existence of a stream of gas being accreted onto the early B-type star.

Starlight subtraction algorithms based on the method of Karhunen–Loève eigenimages have proved invaluable to exoplanet direct imaging. However, they scale poorly in runtime when paired with differential imaging techniques. In such observations, reference frames and frames from which starlight is to be subtracted are drawn from the same set of data, requiring a new subset of references (and eigenimages) for each frame processed to avoid self-subtraction of the signal of interest. The data rates of extreme adaptive optics instruments are such that the only way to make this computationally feasible has been to downsample the data. We develop a technique that updates a precomputed singular value decomposition of the full data set to remove frames (i.e., a "downdate") without a full recomputation, yielding the modified eigenimages. This not only enables analysis of much larger data volumes in the same amount of time, but also exhibits near-linear scaling in runtime as the number of observations increases. We apply this technique to archival data and investigate its scaling behavior for very large numbers of frames N. The resulting algorithm provides speed improvements of 2.6× (for 200 eigenimages at N = 300) to 140× (at N = 10⁴) with the advantage only increasing as N grows. This algorithm has allowed us to substantially accelerate Karhunen–Loève image projection (KLIP) even for modest N, and will let us quickly explore how KLIP parameters affect exoplanet characterization in large-N data sets.

Integrated light spectroscopy from galaxies can be used to study the stellar populations that cannot be resolved into individual stars. This analysis relies on stellar population synthesis (SPS) techniques to study the formation history and structure of galaxies. However, the spectral templates available for SPS are limited, especially in the near-infrared (near-IR). We present A-LIST (APOGEE Library of Infrared SSP Templates), a new set of high-resolution, near-IR SSP spectral templates spanning a wide range of ages (2–12 Gyr), metallicities ( − 2.2 < [M/H] < + 0.4) and α abundances ( − 0.2 < [α/M] < + 0.4). This set of SSP templates is the highest resolution (R ∼ 22, 500) available in the near-IR, and the first such based on an empirical stellar library. Our models are generated using spectra of ∼300,000 stars spread across the Milky Way, with a wide range of metallicities and abundances, from the APOGEE survey. We show that our model spectra provide accurate fits to M31 globular cluster spectra taken with APOGEE, with best-fit metallicities agreeing with those of previous estimates to within ∼0.1 dex. We also compare these model spectra to lower-resolution E-MILES models and demonstrate that we recover the ages of these models to within ∼1.5 Gyr. This library is available in https://github.com/aishashok/ALIST-library.

With an expected torrent of data from the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the need for automated identification of noisy and sparse light curves will increase drastically. In this paper, we performed classification of multiband astronomical light curves from the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) data set via boosted neural nets, boosted decision trees, and a voted classifier for 14 astronomical categories. In order to deal with noisy features, we used wavelet decomposition together with feature selection. We also performed a feature ranking method using a neural network. Our method may be considered an alternative to random forests, which is known to favor features with more categories as relevant. We also investigated the class importance with neural nets using a one-versus-all approach which reduces the multiclass problem to a binary class problem.

The old, solar-metallicity open cluster Messier 67 has long been considered a lynchpin in the study and understanding of the structure and evolution of solar-type stars. The same is arguably true for stellar remnants; the white dwarf population of M67 provides crucial observational data for understanding and interpreting white dwarf populations and evolution. In this work, we determine the white dwarf masses and derive their progenitor star masses using high signal-to-noise spectroscopy of warm (≳10,000 K) DA white dwarfs in the cluster. From this, we are able to derive each white dwarf's position on the initial–final mass relation (IFMR), with an average M_WD = 0.60 ± 0.01 M_⊙ and progenitor mass M_i = 1.52 ± 0.04 M_⊙. These values are fully consistent with recently published linear and piecewise linear fits to the semiempirical IFMR and provide a crucial, precise anchor point for the IFMR for solar-metallicity, low-mass stars. The mean mass of M67 white dwarfs is also consistent with the sharp narrow peak in the local field white dwarf mass distribution, indicating that a majority of recently formed field white dwarfs come from stars with progenitor masses of ≈1.5 M_⊙. Our results enable more precise modeling of the Galactic star formation rate encoded in the field white dwarf mass distribution.

Stellar light curves are well known to encode physical stellar properties. Precise, automated, and computationally inexpensive methods to derive physical parameters from light curves are needed to cope with the large influx of these data from space-based missions such as Kepler and TESS. Here we present a new methodology that we call "The Swan," a fast, generalizable, and effective approach for deriving stellar surface gravity ($\mathrm{log}g$) for main-sequence, subgiant, and red giant stars from Kepler light curves using local linear regression on the full frequency content of Kepler long-cadence power spectra. With this inexpensive data-driven approach, we recover $\mathrm{log}g$ to a precision of ∼0.02 dex for 13,822 stars with seismic $\mathrm{log}g$ values between 0.2 and 4.4 dex and ∼0.11 dex for 4646 stars with Gaia-derived $\mathrm{log}g$ values between 2.3 and 4.6 dex. We further develop a signal-to-noise metric and find that granulation is difficult to detect in many cool main-sequence stars (T_eff ≲ 5500 K), in particular K dwarfs. By combining our $\mathrm{log}g$ measurements with Gaia radii, we derive empirical masses for 4646 subgiant and main-sequence stars with a median precision of ∼7%. Finally, we demonstrate that our method can be used to recover $\mathrm{log}g$ to a similar mean absolute deviation precision for a TESS baseline of 27 days. Our methodology can be readily applied to photometric time series observations to infer stellar surface gravities to high precision across evolutionary states.

The detection and characterization of young planetary systems offer a direct path to study the processes that shape planet evolution. We report on the discovery of a sub-Neptune-sized planet orbiting the young star HD 110082 (TOI-1098). Transit events we initially detected during TESS Cycle 1 are validated with time-series photometry from Spitzer. High-contrast imaging and high-resolution, optical spectra are also obtained to characterize the stellar host and confirm the planetary nature of the transits. The host star is a late-F dwarf (M_⋆ = 1.2M_⊙) with a low-mass, M dwarf binary companion (M_⋆ = 0.26M_⊙) separated by nearly one arcminute (∼6200 au). Based on its rapid rotation and Lithium absorption, HD 110082 is young, but is not a member of any known group of young stars (despite proximity to the Octans association). To measure the age of the system, we search for coeval, phase-space neighbors and compile a sample of candidate siblings to compare with the empirical sequences of young clusters and to apply quantitative age-dating techniques. In doing so, we find that HD 110082 resides in a new young stellar association we designate MELANGE-1, with an age of ${250}_{-70}^{+50}$ Myr. Jointly modeling the TESS and Spitzer light curves, we measure a planetary orbital period of 10.1827 days and radius of R_p = 3.2 ± 0.1R_⊕. HD 110082 b's radius falls in the largest 12% of field-age systems with similar host-star mass and orbital period. This finding supports previous studies indicating that young planets have larger radii than their field-age counterparts.

We use HST/STIS optical spectroscopy of 10 M dwarfs in five closely separated binary systems to test models of M-dwarf structure and evolution. Individual dynamical masses ranging from 0.083 to 0.405 M_⊙ for all stars are known from previous work. We first derive temperature, radius, luminosity, surface gravity, and metallicity by fitting the BT-Settl atmospheric models. We verify that our methodology agrees with empirical results from long-baseline optical interferometry for stars of similar spectral types. We then test whether or not evolutionary models can predict those quantities given the stars' known dynamical masses and the conditions of coevality and equal metallicity within each binary system. We apply this test to five different evolutionary model sets: the Dartmouth models, the MESA/MIST models, the models of Baraffe et al., the PARSEC models, and the YaPSI models. We find marginal agreement between evolutionary model predictions and observations, with few cases where the models respect the condition of coevality in a self-consistent manner. We discuss the pros and cons of each family of models and compare their predictive power.

We present early results from the Epoch of Giant Planet Migration program, a precise radial velocity (RV) survey of more than 100 intermediate-age (∼20–200 Myr) G and K dwarfs with the Habitable Zone Planet Finder spectrograph (HPF) at McDonald Observatory's Hobby–Eberly Telescope. The goals of this program are to determine the timescale and dominant physical mechanism of giant planet migration interior to the water ice line of Sun-like stars. Here, we summarize results from the first 14 months of this program, with a focus on our custom RV pipeline for HPF, a measurement of the intrinsic near-infrared RV activity of young Solar analogs, and modeling the underlying population-level distribution of stellar jitter. We demonstrate on-sky stability at the sub-2 m s⁻¹ level for the K2 standard HD 3765 using a least-squares matching method to extract precise RVs. Based on a subsample of 29 stars with at least three RV measurements from our program, we find a median rms level of 34 m s⁻¹. This is nearly a factor of 2 lower than the median rms level in the optical of 60 m s⁻¹ for a comparison sample with similar ages and spectral types as our targets. The observed near-infrared jitter measurements for this subsample are well reproduced with a log-normal parent distribution with μ = 4.15 and σ = 1.02. Finally, by compiling rms values from previous planet search programs, we show that near-infrared jitter for G and K dwarfs generally decays with age in a similar fashion to optical wavelengths, albeit with a shallower slope and lower overall values for ages ≲1 Gyr.

The occurrence of a planet transiting in front of its host star offers the opportunity to observe the planet's atmosphere filtering starlight. The fraction of occulted stellar flux is roughly proportional to the optically thick area of the planet, the extent of which depends on the opacity of the planet's gaseous envelope at the observed wavelengths. Chemical species, haze, and clouds are now routinely detected in exoplanet atmospheres through rather small features in transmission spectra, i.e., collections of planet-to-star area ratios across multiple spectral bins and/or photometric bands. Technological advances have led to a shrinking of the error bars down to a few tens of parts per million (ppm) per spectral point for the brightest targets. The upcoming James Webb Space Telescope (JWST) is anticipated to deliver transmission spectra with precision down to 10 ppm. The increasing precision of measurements requires a reassessment of the approximations hitherto adopted in astrophysical models, including transit light-curve models. Recently, it has been shown that neglecting the planet's thermal emission can introduce significant biases in the transit depth measured with the JWST/Mid-InfraRed Instrument, integrated between 5 and 12 μm. In this paper, we take a step forward by analyzing the effects of the approximation on transmission spectra over the 0.6–12 μm wavelength range covered by various JWST instruments. We present open-source software to predict the spectral bias, showing that, if not corrected, it may affect the inferred molecular abundances and thermal structure of some exoplanet atmospheres.

Turbulence has the potential for creating gas density enhancements that initiate cloud and star formation (SF), and it can be generated locally by SF. To study the connection between turbulence and SF, we looked for relationships between SF traced by FUV images, and gas turbulence traced by kinetic energy density (KED) and velocity dispersion (v_disp) in the LITTLE THINGS sample of nearby dIrr galaxies. We performed 2D cross-correlations between FUV and KED images, measured cross-correlations in annuli to produce correlation coefficients as a function of radius, and determined the cumulative distribution function of the cross-correlation value. We also plotted on a pixel-by-pixel basis the locally excess KED, v_disp, and H i mass surface density, Σ_HI, as determined from the respective values with the radial profiles subtracted, versus the excess SF rate density Σ_SFR, for all regions with positive excess Σ_SFR. We found that Σ_SFR and KED are poorly correlated. The excess KED associated with SF implies a ∼0.5% efficiency for supernova energy to pump local H i turbulence on the scale of the resolution here, which is a factor of ∼2 too small for all of the turbulence on a galactic scale. The excess v_disp in SF regions is also small, only ∼0.37 km s⁻¹. The local excess in Σ_HI corresponding to an excess in Σ_SFR is consistent with a H i consumption time of ∼1.6 Gyr in the inner parts of the galaxies. The similarity between this timescale and the consumption time for CO implies that CO-dark molecular gas has comparable mass to H i in the inner disks.

We present a catalog of eclipsing binaries in the northern Galactic plane from the Kiso Wide-Field Camera Intensive Survey of the Galactic Plane (KISOGP). We visually identified 7055 eclipsing binaries spread across ∼330 deg², including 4197 W Ursa Majoris/EW-type, 1458 β Lyrae/EB-type, and 1400 Algol/EA-type eclipsing binaries. For all systems, I-band light curves were used to obtain accurate system parameters. We derived the distances and extinction values for the EW-type objects from their period–luminosity relation. We also obtained the structure of the thin disk from the distribution of our sample of eclipsing binary systems, combined with those of high-mass star-forming regions and Cepheid tracers. We found that the thin disk is inhomogeneous in number density as a function of Galactic longitude. Using this new set of distance tracers, we constrain the detailed structure of the thin disk. Finally, we report a global parallax zero-point offset of Δπ = −42.1 ± 1.9 (stat.) ± 12.9 (syst.) μas between our carefully calibrated EW-type eclipsing binary positions and those provided by Gaia Early Data Release 3. Implementation of the officially recommended parallax zero-point correction results in a significantly reduced offset. Additionally, we provide a photometric characterization of our EW-type eclipsing binaries that can be applied to further analyses.

We present infrared photometry of all 36 potential JWST calibrators for which there is archival Spitzer IRAC data. This photometry can then be used to inform the stellar models necessary to provide absolute calibration for all JWST instruments. We describe in detail the steps necessary to measure IRAC photometry from archive retrieval to photometric corrections. To validate our photometry, we examine the distribution of uncertainties from all detections in all four IRAC channels as well as compare the photometry and its uncertainties to those from models, ALLWISE, and the literature. Seventy-five percent of our detections have standard deviations per star of all observations within each channel of less than 3%. The median standard deviations are 1.2%, 1.3%, 1.1%, and 1.9% in [3.6]–[8.0], respectively. We find less than 8% standard deviations in differences of our photometry with ALLWISE and excellent agreement with literature values (less than 3% difference), lending credence to our measured fluxes. JWST is poised to do groundbreaking science, and accurate calibration and cross-calibration with other missions will be part of the underpinnings of that science.

In this paper we present DMC, a model and associated tool for polarimetric imaging of very long baseline interferometry data sets that simultaneously reconstructs the full-Stokes emission structure along with the station-based gain and leakage calibration terms. DMC formulates the imaging problem in terms of posterior exploration, which is achieved using Hamiltonian Monte Carlo sampling. The resulting posterior distribution provides a natural quantification of uncertainty in both the image structure and the data calibration. We run DMC on both synthetic and real data sets, the results of which demonstrate its ability to accurately recover both the image structure and calibration quantities, as well as to assess their corresponding uncertainties. The framework underpinning DMC is flexible, and its specific implementation is under continued development.

We present a comprehensive orbital analysis to the exoplanets β Pictoris b and c that resolves previously reported tensions between the dynamical and evolutionary mass constraints on β Pic b. We use the Markov Chain Monte Carlo orbit code orvara to fit 15 years of radial velocities and relative astrometry (including recent GRAVITY measurements), absolute astrometry from Hipparcos and Gaia, and a single relative radial velocity measurement between β Pic A and b. We measure model-independent masses of ${9.3}_{-2.5}^{+2.6}$M_Jup for β Pic b and 8.3 ± 1.0 M_Jup for β Pic c. These masses are robust to modest changes to the input data selection. We find a well-constrained eccentricity of 0.119 ± 0.008 for β Pic b, and an eccentricity of ${0.21}_{-0.09}^{+0.16}$ for β Pic c, with the two orbital planes aligned to within ∼05. Both planets' masses are within ∼1σ of the predictions of hot-start evolutionary models and exclude cold starts. We validate our approach on N-body synthetic data integrated using REBOUND. We show that orvara can account for three-body effects in the β Pic system down to a level ∼5 times smaller than the GRAVITY uncertainties. Systematics in the masses and orbital parameters from orvara's approximate treatment of multiplanet orbits are a factor of ∼5 smaller than the uncertainties we derive here. Future GRAVITY observations will improve the constraints on β Pic c's mass and (especially) eccentricity, but improved constraints on the mass of β Pic b will likely require years of additional radial velocity monitoring and improved precision from future Gaia data releases.

Motivated by the development of high-dispersion spectrographs in the mid-infrared (MIR) range, we study their application to the atmospheric characterization of nearby nontransiting temperate terrestrial planets around M-type stars. We examine the detectability of CO₂, H₂O, N₂O, and O₃ features in high-resolution planetary thermal emission spectra at 12–18 μm assuming an Earth-like profile and a simplified thermal structure. The molecular line width of such planets can be comparable to or broader than the Doppler shift due to the planetary orbital motion. Given the likely difficulty in knowing the high-resolution MIR spectrum of the host star with sufficient accuracy, we propose observing the target system at two quadrature phases and extracting the differential spectra as the planetary signal. In this case, the signals can be substantially suppressed compared with the case where the host star spectrum is perfectly known, as some parts of the spectral features do not remain in the differential spectra. Despite this self-subtraction, the CO₂ and H₂O features of nearby (≲5 pc) systems with mid-/late-M host stars would be feasible with a 6.5 m class cryogenic space telescope, and orbital inclination could also be constrained for some of them. For CO₂ and N₂O in a 1 bar Earth-like atmosphere, this method would be sensitive when the mixing ratio is 1–10³ ppm. The detectability of molecules except O₃ is not significantly improved when the spectral resolution is higher than ${ \mathcal R }\gtrsim {\rm{10,000}}$, although the constraint on the orbital inclination is improved. This study provides some benchmark cases useful for assessing the value of MIR high-resolution spectroscopy in terms of characterization of potentially habitable planets.

We search for escaping helium from the hot super-Earth 55 Cnc e by taking high-resolution spectra of the 1083 nm line during two transits using Keck/NIRSPEC. We detect no helium absorption down to a 90% upper limit of 250 ppm in excess absorption or 0.27 mÅ in equivalent width. This corresponds to a mass-loss rate of less than ∼10⁹ g s⁻¹ assuming a Parker wind model with a plausible exosphere temperature of 5000–6000 K, although the precise constraint is heavily dependent on model assumptions. We consider both hydrogen- and helium-dominated atmospheric compositions and find similar bounds on the mass-loss rate in both scenarios. Our hydrodynamical models indicate that if a lightweight atmosphere exists on 55 Cnc e, our observations would have easily detected it. Together with the nondetection of Lyα absorption by Ehrenreich et al., our helium nondetection indicates that 55 Cnc e either never accreted a primordial atmosphere in the first place or lost its primordial atmosphere shortly after the dissipation of the gas disk.

We present a CCD UBVI photometric study of poorly studied intermediate-age open cluster SAI 35 (Juchert 20) for the first time. To accomplish this study, we also used LAMOST DR5, Two Micron All Sky Survey, and Gaia EDR3 databases. We identified 214 most probable cluster members with membership probability higher than 50%. The mean proper motion of the cluster is found as ${\mu }_{\alpha }\cos \delta =1.10\pm 0.01$ and μ_δ = −1.66 ± 0.01 mas yr⁻¹. We find the normal interstellar extinction law using the various two-color diagrams. The age, distance, reddening, and radial velocity of the cluster are estimated to be 360 ± 40 Myr, 2.9 ± 0.15 kpc, 0.72 ± 0.05 mag, and −91.62 ± 6.39 km s⁻¹, respectively. The overall mass function slope for main-sequence stars is found to be 1.49 ± 0.16 within the mass range 1.1–3.1 M_⊙, which is in agreement with Salpeter's value within uncertainty. The present study demonstrates that SAI 35 is dynamically relaxed. Galactic orbital parameters are determined using Galactic potential models. We found that this object follows a circular path around the Galactic center.

We measured ³⁵Cl abundances in 52 M giants with metallicities in the range −0.5 < [Fe/H] < 0.12. Abundances and atmospheric parameters were derived using infrared spectra from CSHELL on the NASA Infrared Telescope Facility and from optical echelle spectra. We measured Cl abundances by fitting a H³⁵Cl molecular feature at 3.6985 μm with synthetic spectra. We also measured the abundances of O, Ca, Ti, and Fe using atomic absorption lines. We find that the [Cl/Fe] ratio for our stars agrees with chemical evolution models of Cl, and the [Cl/Ca] ratio is broadly consistent with the solar ratio over our metallicity range. Both indicate that Cl is primarily made in core-collapse supernovae with some contributions from Type Ia supernovae. We suggest that other potential nucleosynthesis processes, such as the ν-process, are not significant producers of Cl. Finally, we also find our Cl abundances are consistent with H ii and planetary nebular abundances at a given oxygen abundance, although there is scatter in the data.

Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multiband image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all of the change points (in time domain), as well as all of the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real data sets show that the proposed method enjoys very promising empirical properties. Applications to two real data sets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it.

We apply the Tremaine–Weinberg method to 19 nearby galaxies using stellar mass surface densities and velocities derived from the PHANGS-MUSE survey, to calculate (primarily bar) pattern speeds (Ω_P). After quality checks, we find that around half (10) of these stellar-mass-based measurements are reliable. For those galaxies, we find good agreement between our results and previously published pattern speeds, and we use rotation curves to calculate major resonance locations (corotation radii and Lindblad resonances). We also compare these stellar-mass-derived pattern speeds with Hα (from MUSE) and CO(J = 2 − 1) emission from the PHANGS-ALMA survey. We find that in the case of these clumpy interstellar medium (ISM) tracers, this method erroneously gives a signal that is simply the angular frequency at a representative radius set by the distribution of these clumps (Ω_clump), and that this Ω_clump is significantly different from Ω_P (∼20% in the case of Hα, and ∼50% in the case of CO). Thus, we conclude that it is inadvisable to use "pattern speeds" derived from ISM kinematics. Finally, we compare our derived pattern speeds and corotation radii, along with bar properties, to the global parameters of these galaxies. Consistent with previous studies, we find that galaxies with a later Hubble type have a larger ratio of corotation radius to bar length, more molecular-gas-rich galaxies have higher Ω_P, and more bulge-dominated galaxies have lower Ω_P. Unlike earlier works, however, there are no clear trends between the bar strength and Ω_P, nor between the total stellar mass surface density and the pattern speed.

The young (50–400 Myr) A3V star β Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot (∼1600 K), warm (∼600 K), and cold (∼120 K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Telescope Interferometer (LBTI) as part of its exozodiacal dust survey (HOSTS). The measured excess is 0.47% ± 0.050% within the central 1.5 au, rising to 0.81% ± 0.026% within 4.5 au, outside the habitable zone of β Leo. This dust level is 50 ± 10 times greater than in the solar system's zodiacal cloud. Poynting–Robertson drag on the cold dust detected by Spitzer, and Herschel underpredicts the dust present in the habitable zone of β Leo, suggesting an additional delivery mechanism (e.g., comets) or an additional belt at ∼5.5 au. A model of these dust components is provided that implies the absence of planets more than a few Saturn masses between ∼5 au and the outer belt at ∼40 au. We also observationally constrain giant planets with the LBTI imaging channel at 3.8 μm wavelength. Assuming an age of 50 Myr, any planet in the system between approximately 5–50 au must be less than a few Jupiter masses, consistent with our dust model. Taken together, these observations showcase the deep contrasts and detection capabilities attainable by the LBTI for both warm exozodiacal dust and giant exoplanets in or near the habitable zone of nearby stars.

Crater chronologies are a fundamental tool to assess the relative and absolute ages of planetary surfaces when direct radiometric dating is not available. Martian crater chronologies are derived from lunar crater spatial densities on terrains with known radiometric ages, and thus they critically depend on the Moon-to-Mars extrapolation. This extrapolation requires knowledge of the time evolution of the impact flux, including contributions from various impactor populations, factors that are not trivially connected to the dynamical evolution of the early Solar System. In this paper, we will present a new Martian crater chronology based on current dynamical models, and consider the main sources of uncertainties (e.g., impactor size–frequency distribution; dynamical models with late and early instabilities, etc.). The resulting "envelope" of Martian crater chronologies significantly differs from previous chronologies. The new Martian crater chronology is discussed using two interesting applications: Jezero crater's dark terrain (relevant to the NASA Mars 2020 mission) and the southern heavily cratered highlands. Our results indicate that Jezero's dark terrain may have formed ∼3.1 Ga, i.e., up to 0.5 Gyr older than previously thought. In addition, available crater chronologies (including our own) overestimate the number of craters larger than 150 km on the southern highlands, suggesting either that large craters have been efficiently erased over Martian history or that dynamical models need further refinement. Further, our chronology constrains the age of Isidis basin to be 4.05–4.2 Ga and that of the Borealis basin to be 4.35–4.40 Ga; these are predictions that can be tested with future sample and return missions.

We study the development of activity in the incoming long-period comet C/2017 K2 over the heliocentric distance range 9 ≲ r_H ≲ 16 au. The comet continues to be characterized by a coma of submillimeter-sized and larger particles ejected at low velocity. In a fixed co-moving volume around the nucleus we find that the scattering cross section of the coma, C, is related to the heliocentric distance by a power law, $C\propto {r}_{H}^{-s}$, with heliocentric index s = 1.14 ± 0.05. This dependence is significantly weaker than the r_H⁻² variation of the insolation as a result of two effects. These are, first, the heliocentric dependence of the dust velocity and, second, a lag effect due to very slow-moving particles ejected long before the observations were taken. A Monte Carlo model of the photometry shows that dust production beginning at r_H ∼ 35 au is needed to match the measured heliocentric index, with only a slight dependence on the particle size distribution. Mass-loss rates in dust at 10 au are of order 10³ kg s⁻¹, while loss rates in gas may be much smaller, depending on the unknown dust to gas ratio. Consequently, the ratio of the nongravitational acceleration to the local solar gravity, α', may, depending on the nucleus size, attain values of ∼10⁻⁷ ≲ α' ≲ 10⁻⁵, comparable to values found in short-period comets at much smaller distances. Nongravitational acceleration in C/2017 K2 and similarly distant comets, while presently unmeasured, may limit the accuracy with which we can infer the properties of the Oort cloud from the orbits of long-period comets.

Estimating stellar ages is important for advancing our understanding of stellar and exoplanet evolution and investigating the history of the Milky Way. However, ages for low-mass stars are hard to infer as they evolve slowly on the main sequence. In addition, empirical dating methods are difficult to calibrate for low-mass stars as they are faint. In this work, we calculate ages for Kepler F, G, and crucially K and M dwarfs, using their rotation and kinematic properties. We apply the simple assumption that the velocity dispersion of stars increases over time and adopt an age–velocity-dispersion relation (AVR) to estimate average stellar ages for groupings of coeval stars. We calculate the vertical velocity dispersion of stars in bins of absolute magnitude, temperature, rotation period, and Rossby number and then convert velocity dispersion to kinematic age via an AVR. Using this method, we estimate gyro-kinematic ages for 29,949 Kepler stars with measured rotation periods. We are able to estimate ages for clusters and asteroseismic stars with an rms of 1.22 Gyr and 0.26 Gyr respectively. With our Astraea machine-learning algorithm, which predicts rotation periods, we suggest a new selection criterion (a weight of 0.15) to increase the size of the McQuillan et al. catalog of Kepler rotation periods by up to 25%. Using predicted rotation periods, we estimated gyro-kinematic ages for stars without measured rotation periods and found promising results by comparing 12 detailed age–element abundance trends with literature values.

We present and analyze 120 spectroscopic binary and triple cluster members of the old (4 Gyr) open cluster M67 (NGC 2682). As a cornerstone of stellar astrophysics, M67 is a key cluster in the WIYN Open Cluster Study (WOCS); radial-velocity (RV) observations of M67 are ongoing and extend back over 45 yr, incorporating data from seven different telescopes, and allowing us to detect binaries with orbital periods ≲10⁴ days. Our sample contains 1296 stars (604 cluster members) with magnitudes of 10 ≤ V ≤ 16.5 (about 1.3–0.7 M_⊙), from the giants down to ∼4 mag below the main-sequence turnoff, and extends in radius to 30' (7.4 pc at a distance of 850 pc, or ∼7 core radii). This paper focuses primarily on the main-sequence binaries, but orbital solutions are also presented for red giants, yellow giants, and sub-subgiants. Out to our period detection limit and within our magnitude and spatial domain, we find a global main-sequence incompleteness-corrected binary fraction of 34% ± 3%, which rises to 70% ± 17% in the cluster center. We derive a tidal circularization period of ${P}_{\mathrm{circ}}={11.0}_{-1.0}^{+1.1}\,\mathrm{days}$. We also analyze the incompleteness-corrected distributions of binary orbital elements and masses. The period distribution rises toward longer periods. The eccentricity distribution, beyond P_circ, is consistent with a uniform distribution. The mass-ratio distribution is also consistent with a uniform distribution. Overall, these M67 binaries are closely consistent with similar binaries in the galactic field, as well as with the old (7 Gyr) open cluster NGC 188. WOCS. 83.

The galaxy NGC 7020 displays an exotic hexagonal ringlike central structure with conspicuous ansae located at two opposite vertices and a tenuous external ring populated by H ii regions. Inside and around the hexagonal structure, Hα emission is also present at the inner disk. To characterize the population of the H ii regions, as well as their ionizing clusters, we imaged NGC 7020 with narrowband Hα and nearby continuum filters attached to GMOS-S installed on the Gemini South telescope. We found 202 H ii regions or complexes of H ii regions evenly distributed between the outer ring and the central disk The nucleus and ansae also present Hα emission. The equivalent width of the Hα line (W_Hα) is systematically greater at the regions of the outer ring relative to those of the inner disk. We discuss the influence of the metallicity gradient of the disk and the upper limit of the masses of the initial mass function on W_Hα, and we conclude that the data are still consistent with the occurrence of a younger burst of H ii region formation in the outer ring. The central regions present more massive clusters, M ≥ 10⁶M_⊙, than those of the outer ring (M ≤ 10⁶M_⊙). Three clusters within 5'' of the nucleus present masses higher than 10⁸M_⊙. The presence of diffuse Hα emission in the inner 5'' suggests gas flows in the nuclear region.

We announce the second data release (DR2) of the NOIRLab Source Catalog (NSC), using 412,116 public images from CTIO-4 m+DECam, the KPNO-4 m+Mosaic3, and the Bok-2.3 m+90Prime. NSC DR2 contains over 3.9 billion unique objects, 68 billion individual source measurements, covers ≈35,000 square degrees of the sky, has depths of ≈23 mag in most broadband filters with ≈1%–2% photometric precision, and astrometric accuracy of ≈7 mas. Approximately 1.9 billion objects within ≈30,000 square degrees of sky have photometry in three or more bands. There are several improvements over NSC DR1. DR2 includes 156,662 (61%) more exposures extending over 2 more years than in DR1. The southern photometric zero-points in griz are more accurate by using the Skymapper DR1 and ATLAS-Ref2 catalogs, and improved extinction corrections were used for high-extinction regions. In addition, the astrometric accuracy is improved by taking advantage of Gaia DR2 proper motions when calibrating the astrometry of individual images. This improves the NSC proper motions to ∼2.5 mas yr⁻¹ (precision) and ∼0.2 mas yr⁻¹ (accuracy). The combination of sources into unique objects is performed using a DBSCAN algorithm and mean parameters per object (such as mean magnitudes, proper motion, etc.) are calculated more robustly with outlier rejection. Finally, eight multi-band photometric variability indices are calculated for each object and variable objects are flagged (23 million objects). NSC DR2 will be useful for exploring solar system objects, stellar streams, dwarf satellite galaxies, quasi-stellar objects, variable stars, high proper-motion stars, and transients. Several examples of these science use cases are presented. The NSC DR2 catalog is publicly available via the NOIRLab's Astro Data Lab science platform.

The pulsation periods of RR Lyrae stars usually vary with time, and they are often used as probes to study the mechanism behind the variation. After the early discovery that the pulsation period of the RR Lyrae star AX UMa decreased rapidly, in further research, we made multiband photometric observations of this star using the Sino-Thai 70 cm telescope and the 60 cm telescope at Yunnan Observatories, and collected its light-curve data from several photometry sky surveys. The O–C diagram confirmed that AX UMa has a rapid period decrease with a rate of −7.752 ± 0.005 days Myr⁻¹, which indicates that it is the fastest-period decreasing ab-type RR Lyrae star in the Galactic field. Moreover, the O – C residuals contain additional periodic variations. We suppose that the variation with a long period is probably caused by the light-travel time effect as the star orbits in a binary system. The calculation shows that the lower mass limit of the companion is about 1 M_⊙. Combined with the full amplitudes and color indexes, we suggested that the companion is probably a hot subdwarf star. We compared the light curves of AX UMa and those of another binary evolution pulsator, OGLE-BLG-RRLYR-02792, and found that the former shows the characteristics of ab-type RR Lyrae stars, while the latter is more like an extreme long-period c-type RR Lyrae star. However, the absence of a bump in the light curves implies that the mass loss has occurred in the outer atmosphere of AX UMa. The special features of AX UMa make it worth more attention and further observations.

We present the discovery and characterization of five hot and warm Jupiters—TOI-628 b (TIC 281408474; HD 288842), TOI-640 b (TIC 147977348), TOI-1333 b (TIC 395171208, BD+47 3521A), TOI-1478 b (TIC 409794137), and TOI-1601 b (TIC 139375960)—based on data from NASA's Transiting Exoplanet Survey Satellite (TESS). The five planets were identified from the full-frame images and were confirmed through a series of photometric and spectroscopic follow-up observations by the TESS Follow-up Observing Program Working Group. The planets are all Jovian size (R_P = 1.01–1.77 R_J) and have masses that range from 0.85 to 6.33 M_J. The host stars of these systems have F and G spectral types (5595 ≤ T_eff ≤ 6460 K) and are all relatively bright (9.5 < V < 10.8, 8.2 < K < 9.3), making them well suited for future detailed characterization efforts. Three of the systems in our sample (TOI-640 b, TOI-1333 b, and TOI-1601 b) orbit subgiant host stars ($\mathrm{log}$g < 4.1). TOI-640 b is one of only three known hot Jupiters to have a highly inflated radius (R_P > 1.7 R_J, possibly a result of its host star's evolution) and resides on an orbit with a period longer than 5 days. TOI-628 b is the most massive, hot Jupiter discovered to date by TESS with a measured mass of ${6.31}_{-0.30}^{+0.28}$M_J and a statistically significant, nonzero orbital eccentricity of e = ${0.074}_{-0.022}^{+0.021}$. This planet would not have had enough time to circularize through tidal forces from our analysis, suggesting that it might be remnant eccentricity from its migration. The longest-period planet in this sample, TOI-1478 b (P = 10.18 days), is a warm Jupiter in a circular orbit around a near-solar analog. NASA's TESS mission is continuing to increase the sample of well-characterized hot and warm Jupiters, complementing its primary mission goals.

Here, we present results on the intrinsic collision probabilities, P_I, and range of collision speeds, V_I, as a function of the heliocentric distance, r, in the trans-Neptunian region. The collision speed is one of the parameters that serves as a proxy for a collisional outcome (e.g., disruption and scattering of fragments, or formation of a crater, as both processes are related to the impact energy). We utilize an improved and debiased model of the trans-Neptunian object (TNO) region from the "Outer Solar System Origins Survey" (OSSOS). It provides a well-defined model of TNO orbital distribution, based on multiple opposition observations of more than 1000 bodies. We compute collisional probabilities for the OSSOS models of the main classical, resonant, detached+outer, and scattering TNO populations. The intrinsic collision probabilities and collision speeds are computed using Öpik's approach, as revised and modified by Wetherill for noncircular and inclined orbits. The calculations are carried out for each of the dynamical TNO groups, allowing for inter-population collisions as well as collisions within each TNO population, resulting in 28 combinations in total. Our results indicate that collisions in the trans-Neptunian region are possible over a wide range in (r, V_I) phase space. Although collisions are calculated to happen within r ∼ 20–200 au and V_I ∼ 0.1 km s⁻¹ to as high as V_I ∼ 9 km s⁻¹, most of the collisions are likely to happen at low relative velocities V_I < 1 km s⁻¹ and are dominated by the main classical belt.

The extragalactic γ-rays sky observed by the Fermi Large Area Telescope (LAT) is dominated by blazars. In the fourth release of the Fermi LAT Point Source Catalog (4FGL) are sources showing a multifrequency behavior similar to that of blazars but lacking an optical spectroscopic confirmation of their nature, known as blazar candidates of uncertain type (BCUs). We aim at confirming the blazar nature of BCUs and test if new optical spectroscopic observations can reveal spectral features, allowing us to get a redshift estimate for known BL Lac objects. We also aim to search for and discover changing-look blazars (i.e., blazars that show a different classification at different epochs). We carry out an extensive search for optical spectra available in the Large Sky Area Multi-object Fibre Spectroscopic Telescope (LAMOST) Data Release 5 (DR5) archive. We select sources out of the 4FGL catalog, the list of targets from our follow-up spectroscopic campaign of unidentified or unassociated γ-ray sources, and the multifrequency catalog of blazars: the Roma-BZCAT. We select a total of 392 spectra. We also compare some of the LAMOST spectra with those available in the literature. We classify 20 BCUs confirming their blazar-like nature. Then we obtain 15 new redshift estimates for known blazars. We discover 26 transitional (i.e., changing-look) blazars that changed their classification. Finally, we are able to confirm the blazar-like nature of six BL Lac candidates. All remaining sources analyzed agree with previous classifications. BL Lac objects are certainly the most elusive type of blazars in the γ-ray extragalactic sky.

It is uncertain whether or not low-mass Population III stars ever existed. While limits on the number density of Population III stars with M_* ≈ 0.8 M_⊙ have been derived, using Sloan Digital Sky Survey (SDSS) data, little is known about the occurrence of Population III stars at lower masses. In the absence of reliable parallaxes, the spectra of metal-poor main-sequence (MPMS) stars with M_* ≲ 0.8 M_⊙ can easily be confused with those of cool white dwarfs. To resolve this ambiguity, in this paper we present a classifier that differentiates between MPMS stars and white dwarfs, based on photometry and/or spectroscopy without the use of parallax information. We build and train our classifier using state-of-the-art theoretical spectra, and evaluate it on existing SDSS-based classifications for objects with reliable Gaia DR2 parallaxes. We then apply our classifier to a large catalog of objects with SDSS photometry and spectroscopy to search for MPMS candidates. We discover several previously unknown, extremely metal-poor (EMP) candidate stars, and recover numerous confirmed EMP stars already noted in the literature. We conclude that archival SDSS spectroscopy has already been exhaustively searched for EMP stars. We predict that the lowest-mass stars of primordial composition will have redder optical-to-infrared colors than cool white dwarfs at constant effective temperature, due to surface gravity-dependent collision-induced absorption from molecular hydrogen. We suggest that the application of our classifier to data produced by next-generation spectroscopic surveys will set stronger constraints on the number density of low-mass Population III stars in the Milky Way.

We present the time-dependent properties of a poorly known OH/IR star, IRAS 18278+0931 (hereafter IRAS 18+09), toward the Ophiuchus constellation. We have carried out long-term optical/near-infrared photometric and spectroscopic observations to study the object. From optical R- and I-band light curves, the period of IRAS 18+09 is estimated to be 575 ± 30 days and the variability amplitudes range from ΔR ∼ 4.0 mag to ΔI ∼ 3.5 mag. From the standard period–luminosity relations, the distance (D) to the object, 4.0 ± 1.3 kpc, is estimated. Applying this distance in the radiative transfer model, the spectral energy distribution is constructed from multiwavelength photometric and IRAS-LRS spectral data, which provide the luminosity, optical depth, and gas mass-loss rate of the object to be 9600 ± 500 L_⊙, 9.1 ± 0.6 at 0.55 μm, and 1.0 × 10⁻⁶M_⊙ yr⁻¹, respectively. The current mass of the object is inferred to be in the range 1.0−1.5 M_⊙ assuming solar metallicity. Notably, the temporal variation of atomic and molecular features (e.g., TiO, Na i, Ca i, CO, H₂O) over the pulsation cycle of the OH/IR star illustrates the sensitivity of the spectral features to the dynamical atmosphere as observed in pulsating AGB stars.

Recently, a noticeable number of new star clusters was identified in the outskirts of the Large Magellanic Cloud (LMC) populating the so-called star-cluster age gap, a space of time (∼4–12 Gyr) where the only known star cluster is up-to-date ESO 121-SC 03. We used Survey of the Magellanic Stellar History DR2 data sets, as well as those employed to identify these star-cluster candidates, to produce relatively deep color–magnitude diagrams (CMDs) of 17 out of 20 discovered age-gap star clusters with the aim of investigating them in detail. Our analysis relies on a thorough CMD cleaning procedure of the field-star contamination, which presents variations in its stellar density and astrophysical properties, such as luminosity and effective temperature, around the star-cluster fields. We built star-cluster CMDs from stars with membership probabilities assigned from the cleaning procedure. These CMDs and their respective spatial distribution maps favor the existence of LMC star field density fluctuations rather than age-gap star clusters, although a definitive assessment on them will be possible from further deeper photometry.

The origin of warm Jupiters (gas giant planets with periods between 10 and 200 days) is an open question in exoplanet formation and evolution. We investigate a particular migration theory in which a warm Jupiter is coupled to a perturbing companion planet that excites secular eccentricity oscillations in the warm Jupiter, leading to periodic close stellar passages that can tidally shrink and circularize its orbit. If such companions exist in warm Jupiter systems, they are likely to be massive and close-in, making them potentially detectable. We generate a set of warm Jupiter-perturber populations capable of engaging in high-eccentricity tidal migration and calculate the detectability of the perturbers through a variety of observational metrics. We show that a small percentage of these perturbers should be detectable in the Kepler light curves, but most should be detectable with precise radial velocity measurements over a 3 month baseline and Gaia astrometry. We find these results to be robust to the assumptions made for the perturber parameter distributions. If a high-precision radial velocity search for companions to warm Jupiters does not find evidence of a significant number of massive companions over a 3 month baseline, it will suggest that perturber-coupled high-eccentricity migration is not the predominant delivery method for warm Jupiters.

The sizes of small planets are known to be bimodal, with a gap separating planets that have lost their primordial atmospheres (super-Earths) and the ones that retain them (mini-Neptunes). Here, we report evidence for another distinct population at smaller sizes. By focusing on planets orbiting around GK dwarfs inward of 16 days and correcting for observational completeness, we find that the number of super-Earths peaks around 1.4 Earth radii and disappears shortly below this size. Instead, a new population of planets (sub-Earths) appears to dominate at sizes below ∼1 Earth radius, with an occurrence that increases with decreasing size. This pattern is also observed in ultra-short-period planets. The end of super-Earths supports earlier claims that super-Earths and mini-Neptunes, planets that likely form in gaseous protoplanetary disks, have a narrow mass distribution. Sub-Earths, in contrast, can be described by a power-law mass distribution and may be explained by the theory of terrestrial planet formation. We therefore speculate that they are formed well after the gaseous disks have dissipated. The extension of these sub-Earths toward longer orbital periods, currently invisible, may be the true terrestrial analogs. This strongly motivates new searches.

Exoplanet transit-timing variations (TTVs) caused by gravitational forces between planets can be used to determine planetary masses and orbital parameters. Most of the observed TTVs are small and sinusoidal in time, leading to degeneracies between the masses and orbital parameters. Here we report a TTV analysis of Kepler-90g and Kepler-90h, which exhibit large TTVs up to 25 hr. With optimization, we find a unique solution that allows us to constrain all of the orbital parameters. The best-fit masses for Kepler-90g and 90h are ${15.0}_{-0.8}^{+0.9}$M_⊕ (Earth mass) and ${203}_{-5}^{+5}{M}_{\oplus }$, respectively, with Kepler-90g having an unusually low apparent density of 0.15 ± 0.05 g cm⁻³. The uniqueness of orbital parameter solution enables a long-term dynamical integration, which reveals that although their periods are close to 2:3 orbital resonance, they are not locked in resonance, and the configuration is stable over billions of years. The dynamical history of the system suggests that planet interactions are able to raise the eccentricities and break the resonant lock after the initial formation.

The planet–metallicity correlation serves as a potential link between exoplanet systems as we observe them today and the effects of bulk composition on the planet formation process. Many observers have noted a tendency for Jovian planets to form around stars with higher metallicities; however, there is no consensus on a trend for smaller planets. Here, we investigate the planet–metallicity correlation for rocky planets in single- and multi-planet systems around Kepler M-dwarf and late-K-dwarf stars. Due to molecular blanketing and the dim nature of these low-mass stars, it is difficult to make direct elemental abundance measurements via spectroscopy. We instead use a combination of accurate and uniformly measured parallaxes and photometry to obtain relative metallicities and validate this method with a subsample of spectroscopically determined metallicities. We use the Kolmogorov–Smirnov (K-S) test, Mann–Whitney U-test, and Anderson–Darling (AD) test to compare the compact multiple planetary systems with single-transiting planet systems and systems with no detected transiting planets. We find that the compact multiple planetary systems are derived from a statistically more metal-poor population, with a p-value of 0.015 in the K-S test, a p-value of 0.005 in the Mann–Whitney U-test, and a value of 2.574 in the AD test statistic, which exceeds the derived threshold for significance by a factor of 25. We conclude that metallicity plays a significant role in determining the architecture of rocky planet systems. Compact multiples either form more readily, or are more likely to survive on gigayear timescales, around metal-poor stars.

In 1972, Zinn, Newell, & Gibson (ZNG) published a list of 156 candidate UV-bright stars they had found in 27 Galactic globular clusters (GCs), based on photographs in the U and V bands. UV-bright stars lie above the horizontal branch (HB) and blueward of the asymptotic giant branch (AGB) and red giant branch in the clusters' color–magnitude diagrams. They are in rapid evolutionary phases—if they are members and not unrelated bright foreground stars. The ZNG list has inspired numerous follow-up studies, aimed at understanding late stages of stellar evolution. However, the ZNG candidates were presented only in finding charts, and celestial coordinates were not given. Using my own collection of CCD frames in u and V, I have identified all of the ZNG objects, and have assembled their coordinates, parallaxes, and proper motions from the recent Gaia Early Data Release 3 (EDR3). Based on the Gaia astrometry, I have determined which objects are probable cluster members (45% of the sample). For the members, using photometry from EDR3, I have assigned the stars to various evolutionary stages, including luminous post-AGB stars, and stars above the HB. I point out several ZNG stars of special interest that have still, to my knowledge, never been studied in detail. This study is an adjunct to a forthcoming survey of the Galactic GCs in the uBVI photometric system, designed for detection of low-gravity stars with large Balmer discontinuities.

Two dwarf galaxies, WOC2017-07 and PGC 704814, located in the vicinity of the nearby luminous spiral galaxy NGC 253 were observed with the Advanced Camera for Surveys on the Hubble Space Telescope. Their distances of 3.62 ± 0.18 Mpc and 3.66 ± 0.18 Mpc were derived using the tip of the red giant branch method. These distances are consistent with the dwarf galaxies being members of the NGC 253 group. Based on the radial velocities and projected separations of seven assumed dwarf companions, we estimated the total mass of NGC 253 to be (8.1 ± 2.6)10¹¹M_⊙, giving a total-mass-to-K-luminosity ratio of M_orb/L_K = (8.5 ± 2.7)M_⊙/L_⊙. A notable property of NGC 253 is its declined rotation curve. NGC 253 joins four other luminous spiral galaxies in the Local Volume with declined rotation curves (NGC 2683, NGC 2903, NGC 3521, and NGC 5055) that together have the low average total-mass-to-luminosity ratio, M_orb/L_K = (5.5 ± 1.1)M_⊙/L_⊙. This value is only ∼1/5 of the corresponding ratio for the Milky Way and M31.

Young massive clusters and super star clusters (SSCs) represent an extreme mode of star formation. Far-infrared imaging of the Magellanic Clouds has identified one potential embedded SSC, HSO BMHERICC J72.971176-69.391112 (H72.97-69.39 in short), in the southwest outskirts of the Large Magellanic Cloud. We present Gemini Flamingos 2 and GSAOI near-infrared imaging of a 3' × 3' region around H72.97-69.39 in order to characterize the stellar content of the cluster. The stellar content is probed down to 1.5 M_⊙. We find substantial dust extinction across the cluster region, extending up to A_K of 3. Deeply embedded stars are associated with ALMA-detected molecular gas suggesting that star formation is ongoing. The high spatial resolution of the GSAOI data allows identification of the central massive object associated with the ¹³CO ALMA observations and detection of fainter low-mass stars around the H30α ALMA source. The morphology of the molecular gas and the nebulosity from adjacent star formation suggest they have interacted covering a region of several parsecs. The total stellar content in the cluster is estimated from the intermediate- and high-mass stellar content to be at least 10,000 M_⊙, less than R136 with up to 100,000 M_⊙ within 4.7 pc radius, but places it in the regime of an SSC. Based on the extinction determination of individual stars we estimate a molecular gas mass in the vicinity of H72.97-69.39 of 6600 M_⊙, suggesting more star formation can be expected.