Table of contents

We present high angular resolution imaging of the quasar PSO J172.3556+18.7734 at z = 6.82 with the Very Long Baseline Array (VLBA). This source currently holds the record of being the highest redshift radio-loud quasar. These observations reveal a dominant radio source with a flux density of 398.4 ± 61.4 μJy at 1.53 GHz, a deconvolved size of 9.9 × 3.5 mas (52.5 × 18.6 pc), and an intrinsic brightness temperature of (4.7 ± 0.7) × 10⁷ K. A weak unresolved radio extension from the main source is also detected at the ∼3.1σ level. The total flux density recovered with the VLBA at 1.53 GHz is consistent with that measured with the Very Large Array (VLA) at a similar frequency. The quasar is not detected at 4.67 GHz with the VLBA, suggesting a steep spectral index with a limit of ${\alpha }_{4.67}^{1.53}\lt -$1.55. The quasar is also not detected with the VLBA at 7.67 GHz. The overall characteristics of the quasar suggest that it is a very young radio source similar to lower redshift gigahertz peaked spectrum radio sources, with an estimated kinematic age of ∼10³ yr. The VLA observations of this quasar revealed a second radio source in the field 231 away. This radio source, which does not have an optical or IR counterpart, is not detected with the VLBA at any of the observed frequencies. Its nondetection at the lowest observed VLBA frequency suggests that it is resolved out, implying a size larger than ∼017. It is thus likely situated at lower redshift than the quasar.

The evanescent wave coronagraph uses the principle of frustrated total internal reflection (FTIR) to suppress the light coming from the star and study its close environment. Its focal plane mask is composed of a lens and a prism placed in contact with each other to produce the coronagraphic effect. In this paper, we present the experimental results obtained using an upgraded focal plane mask of the Evanescent Wave Coronagraph (EvWaCo). These experimental results are also compared to the theoretical performance of the coronagraph obtained through simulations. Experimentally, we reach a raw contrast equal to a few 10⁻⁴ at a distance equal to 3 λ/D over the full I band (λ_c = 800 nm, Δλ/λ ≈ 20%) and equal to 4 λ/D over the full R band (λ_c = 650 nm, Δλ/λ ≈ 23%) in unpolarized light. However, our simulations show a raw contrast close to 10⁻⁴ over the full I band and R band at the same distance, thus confirming the theoretical achromatic advantage of the coronagraph. We also verify the stability of the mask through a series of contrast measurements over a period of 8 months. Furthermore, we measure the sensitivity of the coronagraph to the lateral and longitudinal misalignment of the focal plane mask and to the lateral misalignment of the Lyot stop.

We present high-resolution near-infrared spectra taken during eight transits of 55 Cancri e, a nearby low-density super-Earth with a short orbital period (<18 hr). While this exoplanet's bulk density indicates a possible atmosphere, one has not been detected definitively. Our analysis relies on the Doppler cross-correlation technique, which takes advantage of the high spectral resolution and broad wavelength coverage of our data, to search for the thousands of absorption features from hydrogen-, carbon-, and nitrogen-rich molecular species in the planetary atmosphere. Although we are unable to detect an atmosphere around 55 Cancri e, we do place strong constraints on the levels of HCN, NH₃, and C₂H₂ that may be present. In particular, at a mean molecular weight of 5 amu, we can rule out the presence of HCN in the atmosphere down to a volume mixing ratio (VMR) of 0.02%, NH₃ down to a VMR of 0.08%, and C₂H₂ down to a VMR of 1.0%. If the mean molecular weight is relaxed to 2 amu, we can rule out the presence of HCN, NH₃, and C₂H₂ down to VMRs of 0.001%, 0.0025%, and 0.08%, respectively. Our results reduce the parameter space of possible atmospheres consistent with the analysis of Hubble Space Telescope/WFC3 observations by Tsiaras et al. and indicate that if 55 Cancri e harbors an atmosphere, it must have a high mean molecular weight or clouds.

We observed the occultation of the star Gaia DR2 4056440205544338944 by (28978) Ixion. The event was observed from two Lowell Observatory sites, using the 4.3 m Lowell Discovery Telescope (LDT), near Happy Jack, AZ, USA, and a 0.32 m telescope co-mounted with the Titan Monitoring telescope on Lowell's Mars Hill campus in Flagstaff, AZ. The LDT chord, at 44.86 s, was roughly 30% longer than the longest predicted possible chord. Under the assumption of a spherical body, Ixion's fitted diameter D = 709.6 ± 0.2 km. The LDT light-curve profile was used to place an upper limit on the surface pressure P < 2 μbar on any possible atmosphere of Ixion. At the distance of Ixion, the occulted star had a fitted projected diameter of 19.25 ± 0.3 km assuming uniform disk illumination, giving a stellar angular diameter of 0.675 ± 0.010 mas. Using the Gaia EDR3 parallax of 0.565 mas, the stellar radius is ${130}_{-17}^{+20}\,{R}_{\odot }$. The measured size is consistent with prior spectral classification of this star as a reddened mid-M giant. This is one of only a modest number of M5 III stars to have a directly measured diameter, and is more distant than most.

Using a suite of numerical calculations, we consider the long-term evolution of circumbinary debris from the Pluto–Charon giant impact. Initially, these solids have large eccentricity and pericenters near Charon's orbit. On timescales of 100–1000 yr, dynamical interactions with Pluto and Charon lead to the ejection of most solids from the system. As the dynamics moves particles away from the barycenter, collisional damping reduces the orbital eccentricity of many particles. These solids populate a circumbinary disk in the Pluto–Charon orbital plane; a large fraction of this material lies within a "satellite zone" that encompasses the orbits of Styx, Nix, Kerberos, and Hydra. Compared to the narrow rings generated from the debris of a collision between a trans-Neptunian object and Charon, disks produced after the giant impact are much more extended and may be a less promising option for producing small circumbinary satellites.

We investigate the connection between galactic outflows and star formation using two independent data sets covering a sample of 22 galaxies between 1 ≲ z ≲ 1.5. The Hubble Space Telescope WFC3/G141 grism provides low spectral resolution, high spatial resolution spectroscopy yielding Hα emission-line maps from which we measure the spatial extent and strength of star formation. In the rest-frame near-UV, Keck/DEIMOS observes Fe ii and Mg ii interstellar absorption lines, which provide constraints on the intensity and velocity of the outflows. We compare outflow properties from individual and composite spectra with the star formation rate (SFR) and SFR surface density (Σ_SFR), as well as the stellar mass and specific SFR (sSFR). The Fe ii and Mg ii equivalent widths (EWs) increase with both SFR and Σ_SFR at ≳3σ significance, while the composite spectra show larger Fe ii EWs and outflow velocities in galaxies with higher SFR, Σ_SFR, and sSFR. Absorption-line profiles of the composite spectra further indicate that the differences between subsamples are driven by outflows rather than the interstellar medium. While these results are consistent with those of previous studies, the use of Hα images makes them the most direct test of the relationship between star formation and outflows at z > 1 to date. Future facilities such as the James Webb Space Telescope and the upcoming Extremely Large Telescopes will extend these direct, Hα-based studies to lower masses and SFRs, probing galactic feedback across orders of magnitude in galaxy properties and augmenting the correlations we find here.

We report the detection of an atmosphere on a rocky exoplanet, GJ 1132 b, which is similar to Earth in terms of size and density. The atmospheric transmission spectrum was detected using Hubble WFC3 measurements and shows spectral signatures of aerosol scattering, HCN, and CH₄ in a low mean molecular weight atmosphere. We model the atmospheric loss process and conclude that GJ 1132 b likely lost the original H/He envelope, suggesting that the atmosphere that we detect has been reestablished. We explore the possibility of H₂ mantle degassing, previously identified as a possibility for this planet by theoretical studies, and find that outgassing from ultra-reduced magma could produce the observed atmosphere. In this way we use the observed exoplanet transmission spectrum to gain insights into magma composition for a terrestrial planet. The detection of an atmosphere on this rocky planet raises the possibility that the numerous powerfully irradiated super-Earth planets, believed to be the evaporated cores of sub-Neptunes, may, under favorable circumstances, host detectable atmospheres.

Gaia Early Data Release 3 (EDR3) provides trigonometric parallaxes for 1.5 billion stars, with reduced systematics compared to Gaia Data Release 2 and reported precisions better by up to a factor of 2. New to EDR3 is a tentative model for correcting the parallaxes of magnitude-, position-, and color-dependent systematics for five- and six-parameter astrometric solutions, Z₅ and Z₆. Using a sample of over 2000 first-ascent red giant branch stars with asteroseismic parallaxes, I perform an independent check of the Z₅ model in a Gaia magnitude range of 9 ≲ G ≲ 13 and color range of 1.4 μm⁻¹ ≲ ν_eff ≲ 1.5 μm⁻¹. This analysis therefore bridges the Gaia team's consistency check of Z₅ for G > 13 and indications from independent analysis using Cepheids of a ≈15 μas overcorrection for G < 11. I find overcorrection sets in at G ≲ 10.8, such that Z₅-corrected EDR3 parallaxes are larger than asteroseismic parallaxes by 15 ± 3 μas. For G ≳ 10.8, EDR3 and asteroseismic parallaxes in the Kepler field agree up to a constant consistent with expected spatial variations in EDR3 parallaxes after a linear, color-dependent adjustment. I also infer an average underestimation of EDR3 parallax uncertainties in the sample of 22% ± 6%, consistent with the Gaia team's estimates at similar magnitudes and independent analysis using wide binaries. Finally, I extend the Gaia team's parallax spatial covariance model to brighter magnitudes (G < 13) and smaller scales (down to ≈01), where systematic EDR3 parallax uncertainties are at least ≈3–4 μas.

The Ultraviolet Imaging Telescope (UVIT) on board the AstroSat observatory has imaged the Andromeda galaxy (M31) from 2017 to 2019 in the far- and near-UV (FUV and NUV) with the high spatial resolution of ≃1''. The survey covered the large sky area of M31 with a set of observations (Fields), each 28' in diameter. Field 1 was observed in two epochs with the F148W filter, separated by ≃1133 days (≃3.10 yr). The 6.4 kpc diameter Field 1 (at the distance of M31) includes a substantial part of the inner spiral arms of the galaxy. We identify UVIT sources in both epochs of Field 1 and obtain catalogs of sources that are variable in FUV at >3σ and >5σ confidence level. The fraction of FUV-variable sources is higher for brighter sources, and the fraction is higher in the two main spiral arms compared to other areas. This is evidence that a significant fraction of the FUV variables are associated with hot young stars. Source counterparts are found for 42 of the 86 >5σ FUV variables using existing catalogs. The counterparts include 10 star clusters, 6 H II regions, 5 regular or semiregular variables, 6 other variables, and 6 nova or nova candidates. The UVIT FUV–NUV and FUV–FUV color–magnitude diagrams confirm the association of most of the FUV variables with hot young stars. A catalog of UVIT photometry for the variable sources is presented.

The Active Optics System of the Vera C. Rubin Observatory (Rubin) uses information provided by four wave front sensors to determine deviations between the reconstructed wave front and the ideal wave front. The observed deviations are used to adjust the control parameters of the optical system to maintain image quality across the 35 field of view. The baseline approach from the project is to obtain amplitudes of the Zernike polynomials describing the distorted wave front from out-of-focus images collected by the wave front sensors. These Zernike amplitudes are related via an "influence matrix" to the control parameters necessary to correct the wave front. In this paper, we use deep-learning methods to extract the control parameters directly from the images captured by the wave front sensors. Our neural net model uses anti-aliasing pooling to boost performance, and a domain-specific loss function to aid learning and generalization. The accuracy of the control parameters derived from our model exceeds Rubin requirements even in the presence of full-moon background levels and mis-centering of reference stars. Although the training process is time consuming, model evaluation requires only a few milliseconds. This low latency should allow for the correction of the optical configuration during the readout and slew interval between successive exposures.

Measurements of gas mass in protoplanetary gas disks form the basis for estimating the conditions of planet formation. Among the most important constraints derived from disk diagnostics are the abundances of gas-phase species critical for understanding disk chemistry. Toward this end, we present direct line-of-sight measurements of H₂ and CO, employing UV absorption spectroscopy from Hubble Space Telescope Cosmic Origins Spectrograph to characterize disk composition, molecular excitation temperatures, and spatial distribution in the circumstellar material around the Herbig Ae stars HK Ori and T Ori. We observe strong CO (N(CO) = 10^15.5 cm⁻²; T_rot(CO) = 19 K) and H₂ (N(H₂) = 10^20.34 cm⁻²; T_rot(H₂) = 141 K) absorption toward HK Ori with a CO/H₂ ratio of (≡N(CO)/N(H₂)) = ${1.3}_{-0.7}^{+1.6}$ × 10⁻⁵. These measurements place direct empirical constraints on the CO-to-H₂ conversion factor in the disk around a Herbig Ae star for the first time, although there is uncertainty concerning the exact viewing geometry of the disk. The spectra of T Ori show CO (N(CO) = 10^14.9 cm⁻²; T_rot(CO) = 124 K) absorption. Interestingly, we do not detect any H₂ absorption toward this star (N(H₂) < 10^15.9 cm⁻²). We discuss a potential scenario for the detection of CO without H₂, which deserves further investigation. The low abundance ratio measured around HK Ori suggests significant depletion of CO in the circumstellar gas, which conforms with the handful of other recent CO abundance measurements in protoplanetary disks.

We present Tails, an open-source deep-learning framework for the identification and localization of comets in the image data of the Zwicky Transient Facility (ZTF), a robotic optical time-domain survey currently in operation at the Palomar Observatory in California, USA. Tails employs a custom EfficientDet-based architecture and is capable of finding comets in single images in near real time, rather than requiring multiple epochs as with traditional methods. The system achieves state-of-the-art performance with 99% recall, a 0.01% false-positive rate, and a 1–2 pixel rms error in the predicted position. We report the initial results of the Tails efficiency evaluation in a production setting on the data of the ZTF Twilight survey, including the first AI-assisted discovery of a comet (C/2020 T2) and the recovery of a comet (P/2016 J3 = P/2021 A3).

K2 greatly extended Kepler's ability to find new planets, but it was typically limited to identifying transiting planets with orbital periods below 40 days. While analyzing K2 data through the Exoplanet Explorers project, citizen scientists helped discover one super-Earth and four sub-Neptune sized planets in the relatively bright (V = 12.21, K = 10.3) K2-138 system, all which orbit near 3:2 mean-motion resonances. The K2 light curve showed two additional transit events consistent with a sixth planet. Using Spitzer photometry, we validate the sixth planet's orbital period of 41.966 ± 0.006 days and measure a radius of ${3.44}_{-0.31}^{+0.32}\,{R}_{\oplus }$, solidifying K2-138 as the K2 system with the most currently known planets. There is a sizeable gap between the outer two planets, since the fifth planet in the system, K2-138 f, orbits at 12.76 days. We explore the possibility of additional nontransiting planets in the gap between f and g. Due to the relative brightness of the K2-138 host star, and the near resonance of the inner planets, K2-138 could be a key benchmark system for both radial velocity and transit-timing variation mass measurements, and indeed radial velocity masses for the inner four planets have already been obtained. With its five sub-Neptunes and one super-Earth, the K2-138 system provides a unique test bed for comparative atmospheric studies of warm to temperate planets of similar size, dynamical studies of near-resonant planets, and models of planet formation and migration.

We report the discovery of a likely outbursting Class I young stellar object, associated with the star-forming region NGC 281-W (distance ∼2.8 kpc). The source is currently seen only at infrared wavelengths, appearing in both the Palomar Gattini InfraRed (1.2 μm) and the Near-Earth Object Wide-field Infrared Survey Explorer (3.4 and 4.6 μm) photometric time-domain surveys. Recent near-infrared imaging reveals a new, extended scattered light nebula. Recent near-infrared spectroscopy confirms the similarity of PGIR 20dci to FU Ori–type sources, based on strong molecular absorption in CO, H₂O, and OH; weak absorption in several atomic lines; and a warm wind/outflow as indicated by a P Cygni profile in the He iλ10830 line. This is a rare case of an FU Ori star with a well-measured long-term photometric rise before a sharper outburst, and the second instance of an FU Ori star with a documented two-step brightening in the mid-infrared.

We present the photometric and spectroscopic analysis of four W UMa binaries J015829.5+260333 (hereinafter as J0158), J030505.1+293443 (hereinafter as J0305), J102211.7+310022 (hereinafter as J1022), and KW Psc. The VR_cI_c band photometric observations are carried out with the 1.3 m Devasthal Fast Optical Telescope (DFOT). For low-resolution spectroscopy, we used the 2 m Himalayan Chandra Telescope (HCT) as well as the archival data from the 4 m LAMOST survey. The systems J0158 and J0305 show a period increase rate of 5.26( ± 1.72) × 10⁻⁷ days yr⁻¹ and 1.78( ± 1.52) × 10⁻⁶ days yr⁻¹, respectively. The period of J1022 is found to be decreasing with a rate of 4.22 ( ± 1.67) × 10⁻⁶ days yr⁻¹. The period analysis of KW Psc displays no change in its period. The PHOEBE package is used for the light-curve modeling and basic parameters are evaluated with the help of the GAIA parallax. The asymmetry of light curves is explained with the assumption of cool spots at specific positions on one of the components of the system. On the basis of temperatures, mass ratios, fill-out factors, and periods, the system J1022 is identified as a W-subtype system while the others show some mixed properties. To probe the chromospheric activities in these W UMa binaries, their spectra are compared with the known inactive stars' spectra. The comparison shows emission in H_α, H_β, and Ca II. To understand the evolutionary status of these systems, the components are plotted in mass–radius and mass–luminosity planes with other well characterized binary systems. The secondary components of all the systems are away from ZAMS, which indicates that the secondary is more evolved than the primary component.

It is well known that the chances of success of the search for extraterrestrial intelligence depend on the longevity of technological civilizations or, more broadly, on the duration of the signs of their existence, or technosignatures. Here, we re-examine this general tenet in more detail, and we show that its broader implications have not been given their proper significance. In particular, an often overlooked aspect is that the duration of a technosignature is in principle almost entirely separable from the age of the civilization that produces it. We propose a classification scheme of technosignatures based on their duration, and we use Monte Carlo simulations to show that, given an initial generic distribution of Galactic technosignatures, only the ones with the longest duration are likely to be detected. This tells us, among other things, that looking for a large number of short-lived technosignatures is a weaker observational strategy than focusing the search on a few long-lived ones. It also suggests abandoning any anthropocentric bias in approaching the question of extraterrestrial intelligence. We finally give some ideas of possible pathways that can lead to the establishment of long-lived technosignatures.

We present the results of a study of the prospect of detecting habitable Trojan planets in the Kepler Habitable Zone circumbinary planetary systems (Kepler-16, -47, -453, -1647, and -1661). We integrated the orbits of 10,000 separate N-body systems (N = 4, 6), each with a one Earth-mass body in a randomly selected orbit near the L₄ and L₅ Lagrangian points of the host HZ circumbinary planet. We find that stable Trojan planets are restricted to a narrow range of semimajor axes in all five systems and limited to small eccentricities in Kepler-16, -47, and -1661. To assess the prospect of the detection of these habitable Trojan planets, we calculated the amplitudes of the variations they cause in the transit timing of their host bodies. Results show that the mean amplitudes of the transit timing variations (TTVs) correlate with the mass of the transiting planet and range from 70 minutes for Kepler-16b to 390 minutes for Kepler-47c. Our analysis indicates that the TTVs of the circumbinary planets caused by these Trojan bodies fall within the detectable range of timing precision obtained from the Kepler telescope's long-cadence data. The latter points to Kepler data as a viable source to search for habitable Trojan planets.

We present the discovery of rapid photometric variability in three ultra-cool dwarfs from long-duration monitoring with the Spitzer Space Telescope. The T7, L3.5, and L8 dwarfs have the shortest photometric periods known to date: ${1.080}_{-0.005}^{+0.004}$ hr, ${1.14}_{-0.01}^{+0.03}$ hr, and ${1.23}_{-0.01}^{+0.01}$ hr, respectively. We confirm the rapid rotation through moderate-resolution infrared spectroscopy, which reveals projected rotational velocities between 79 and 104 km s⁻¹. We compare the near-infrared spectra to photospheric models to determine the objects' fundamental parameters and radial velocities. We find that the equatorial rotational velocities for all three objects are ≳100 km s⁻¹. The three L and T dwarfs reported here are the most rapidly spinning and likely the most oblate field ultra-cool dwarfs known to date. Correspondingly, all three are excellent candidates for seeking auroral radio emission and net optical/infrared polarization. As of this writing, 78 L-, T-, and Y-dwarf rotation periods have now been measured. The clustering of the shortest rotation periods near 1 hr suggests that brown dwarfs are unlikely to spin much faster.

SW Sextantis systems are nova-like cataclysmic variables that have unusual spectroscopic properties, which are thought to be caused by an accretion geometry having part of the mass flux trajectory out of the orbital plane. Accretion onto a magnetic white dwarf is one of the proposed scenarios for these systems. To verify this possibility, we analyzed photometric and polarimetric time-series data for a sample of six SW Sex stars. We report possible modulated circular polarization in BO Cet, SW Sex, and UU Aqr with periods of 11.1, 41.2, and 25.7 minutes, respectively, and less significant periodicities for V380 Oph at 22 minutes and V442 Oph at 19.4 minutes. We confirm previous results that LS Peg shows variable circular polarization. However, we determine a period of 18.8 minutes, which is different from the earlier reported value. We interpret these periods as the spin periods of the white dwarfs. Our polarimetric results indicate that 15% of the SW Sex systems have direct evidence of magnetic accretion. We also discuss SW Sex objects within the perspective of being magnetic systems, considering the latest findings about the demography, formation, and evolution of cataclysmic variables.

We study the radio and optical properties of the brightest group galaxies (BGGs) in a sample of galaxy groups from the Sloan Digital Sky Survey (SDSS) DR7. The luminosity difference between the BGG and the second-ranked galaxy in the group (known as the luminosity, or magnitude, gap) has been used as a probe for the level of galaxy interaction for the BGG within the group. We study the properties of BGGs with magnitude gaps in the range of 0–2.7 mag, in order to investigate any relation between the luminosity gap and the radio properties of the BGG. In order to eliminate selection biases, we ensure that all variations in stellar mass are accounted for. We then confirm that, at fixed stellar mass, there are no significant variations in the optical properties of the BGGs over the full range of luminosity gaps studied. We compare these optical results with the Evolution and Assembly of GaLaxies and their Environments (EAGLE) hydrodynamical simulations and find broad consistency with the observational data. Using EAGLE we also confirm that no trends begin to arise in the simulated data at luminosity gaps beyond our observational limits. Finally, we find that, at a fixed stellar mass, the fraction of BGGs that are radio-loud also shows no trend as a function of luminosity gap. We examine how the BGG offset from the center of the group may affect the radio results and find no significant trend for the fraction of radio-loud BGGs with a magnitude gap in either the BGG samples with greater or less than 100 kpc offset from the center of the group.

We have developed a method that maps large astronomical images onto a two-dimensional map and clusters them. A combination of various state-of-the-art machine-learning algorithms is used to develop a fully unsupervised image-quality assessment and clustering system. Our pipeline consists of a data pre-processing step where individual image objects are identified in a large astronomical image and converted to smaller pixel images. This data is then fed to a deep convolutional auto-encoder jointly trained with a self-organizing map (SOM). This part can be used as a recommendation system. The resulting output is eventually mapped onto a two-dimensional grid using a second, deep, SOM. We use data taken from ground-based telescopes and, as a case study, compare the system's ability and performance with the results obtained by supervised methods presented by Teimoorinia et al. The availability of target labels in this data allowed for a comprehensive performance comparison between our unsupervised and supervised methods. In addition to image-quality assessments performed in this project, our method can have various other applications. For example, it can help experts label images in a considerably shorter time with minimum human intervention. It can also be used as a content-based recommendation system capable of filtering images based on the desired content.

Next-generation space observatories will conduct the first systematic surveys of terrestrial exoplanet atmospheres and search for evidence of life beyond Earth. While in-depth observations of the nearest habitable worlds may yield enticing results, there are fundamental questions about planetary habitability and evolution that can only be answered through population-level studies of dozens to hundreds of terrestrial planets. To determine the requirements for next-generation observatories to address these questions, we have developed Bioverse. Bioverse combines existing knowledge of exoplanet statistics with a survey simulation and hypothesis testing framework to determine whether proposed space-based direct imaging and transit-spectroscopy surveys will be capable of detecting various hypothetical statistical relationships between the properties of terrestrial exoplanets. Following a description of the code, we apply Bioverse to determine whether an ambitious direct imaging or transit survey would be able to determine the extent of the circumstellar habitable zone and study the evolution of Earth-like planets. Given recent evidence that Earth-sized habitable zone planets are likely much rarer than previously believed, we find that space missions with large search volumes will be necessary to study the population of terrestrial and habitable worlds. Moving forward, Bioverse provides a methodology for performing trade studies of future observatory concepts to maximize their ability to address population-level questions, including and beyond the specific examples explored here.

In this paper, we present the first results from a CARMA high-resolution ¹²CO(1-0), ¹³CO(1-0), and C¹⁸O(1-0) molecular line survey of the North America and Pelican (NAP) Nebulae. CARMA observations have been combined with single-dish data from the Purple Mountain 13.7 m telescope, to add short spacings and to produce high-dynamic-range images. We find that the molecular gas is predominantly shaped by the W80 H ii bubble, driven by an O star. Several bright rims noted in the observation are probably remnant molecular clouds, heated and stripped by the massive star. Matching these rims in molecular lines and optical images, we construct a model of the three-dimensional structure of the NAP complex. Two groups of molecular clumps/filaments are on the near side of the bubble: one is being pushed toward us, whereas the other is moving toward the bubble. Another group is on the far side of the bubble, and moving away. The young stellar objects in the Gulf region reside in three different clusters, each hosted by a cloud from one of the three molecular clump groups. Although all gas content in the NAP is impacted by feedback from the central O star, some regions show no signs of star formation, while other areas clearly exhibit star formation activity. Additional molecular gas being carved by feedback includes cometary structures in the Pelican Head region, and the boomerang features at the boundary of the Gulf region. The results show that the NAP complex is an ideal place for the study of feedback effects on star formation.

Radial–velocity (RV) planet searches are often polluted by signals caused by gas motion at the star's surface. Stellar activity can mimic or mask changes in the RVs caused by orbiting planets, resulting in false positives or missed detections. Here we use Gaussian process regression to disentangle the contradictory reports of planets versus rotation artifacts from Kapteyn's star. To model rotation, we use joint quasiperiodic kernels for the RV and Hα signals, requiring that their periods and correlation timescales be the same. We find that the rotation period of Kapteyn's star is 125 days, while the characteristic active-region lifetime is 694 days. Adding a planet to the RV model produces a best-fit orbital period of 100 yr, or 10 times the observing time baseline, indicating that the observed RVs are best explained by star rotation only. We also find no significant periodic signals in residual RV data sets constructed by subtracting off realizations of the best-fit rotation model and conclude that both previously reported "planets" are artifacts of the star's rotation and activity. Our results highlight the pitfalls of using sinusoids to model quasiperiodic rotation signals.

The Kepler Mission revolutionized exoplanet science and stellar astrophysics by obtaining highly precise photometry of over 200,000 stars over 4 yr. A critical piece of information to exploit Kepler data is its selection function, since all targets had to be selected from a sample of half a million stars on the Kepler CCDs using limited information. Here we use Gaia DR2 to reconstruct the Kepler selection function and explore possible biases with respect to evolutionary state, stellar multiplicity, and kinematics. We find that the Kepler target selection is nearly complete for stars brighter than Kp < 14 mag and was effective at selecting main-sequence stars, with the fraction of observed stars decreasing from 95% to 60% between 14 < Kp < 16 mag. We find that the observed fraction for subgiant stars is only 10% lower, confirming that a significant number of subgiants selected for observation were believed to be main-sequence stars. Conversely we find a strong selection bias against low-luminosity red giant stars (R ≈ 3–5R_⊙, T_eff ≈ 5500 K), dropping from 90% at Kp = 14 mag to below 30% at Kp = 16 mag, confirming that the target selection was efficient at distinguishing dwarfs from giants. We compare the Gaia Re-normalized Unit Weight Error (RUWE) values of the observed and nonobserved main-sequence stars and find a difference in elevated (>1.2) RUWE values at ∼σ significance, suggesting that the Kepler target selection shows some bias against either close or wide binaries. We furthermore use the Gaia proper motions to show that the Kepler selection function was unbiased with respect to kinematics.

We develop a Hamiltonian analytical theory for the rotation of a Poincaré Earth model (rigid mantle and liquid core) at the second order with respect to the lunisolar potential and moving ecliptic term. Since the Andoyer variables considered in the first-order solution present virtual singularities, i.e., vanishing divisors, we introduce a set of nonsingular complex canonical variables. This choice allows for applying the Hori canonical perturbation method in a standard way. We derive analytical expressions for the first- and second-order solutions of the precession and nutation of the angular momentum axis (Poisson terms). Contrary to first-order theories, there is a part of the Poisson terms that does depend on the Earth's structure. The resulting numerical amplitudes, not incorporated in the International Astronomical Union nutation standard, are not negligible considering current accuracies. They are at the microarcsecond level for a few terms, with a very significant contribution in obliquity of about 40 μas for the nutation argument with period −6798.38 days. The structure-dependent amplitudes present a large amplification with respect to the rigid model due to the fluid core resonance. The features of such resonance, however, are different from those found in first-order solutions. The most prominent is that it does not depend directly on the second-order nutation argument but rather on the combination of first-order arguments generating it. It entails that some first-order approaches, like those based on the transfer function, cannot be applied to obtain the second-order contributions.

In the search for life in the cosmos, NASA's Transiting Exoplanet Survey Satellite (TESS) mission has already monitored about 74% of the sky for transiting extrasolar planets, including potentially habitable worlds. However, TESS only observed a fraction of the stars long enough to be able to find planets like Earth. We use the primary mission data—the first two years of observations—and identify 4239 stars within 210 pc that TESS observed long enough to see three transits of an exoplanet that receives similar irradiation to Earth: 738 of these stars are located within 30 pc. We provide reliable stellar parameters from the TESS Input Catalog that incorporates Gaia DR2 and also calculate the transit depth and radial velocity semiamplitude for an Earth-analog planet. Of the 4239 stars in the Revised TESS HZ Catalog, 9 are known exoplanet hosts—GJ1061, GJ1132, GJ3512, GJ685, Kepler-42, LHS1815, L98-59, RRCae, and TOI700—around which TESS could identify additional Earth-like planetary companions. Thirty-seven additional stars host yet unconfirmed TESS Objects of Interest: three of these orbit in the habitable zone—TOI203, TOI715, and TOI2298. For a subset of 614 of the 4239 stars, TESS has observed the star long enough to be able to observe planets throughout the full temperate, habitable zone out to the equivalent of Mars' orbit. Thus, the Revised TESS Habitable Zone Catalog provides a tool for observers to prioritize stars for follow-up observation to discover life in the cosmos. These stars are the best path toward the discovery of habitable planets using the TESS mission data.

We present a Bayesian method to cross-match 5,827,988 high proper-motion Gaia sources (μ > 40 mas yr⁻¹) to various photometric surveys: Two Micron All Sky Survey, AllWISE data release from the Wide-field Infrared Explorer (WISE) mission, Galaxy Evolution Explorer, Radial Velocity Experiment, Sloan Digital Sky Survey, and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). To efficiently associate these objects across catalogs, we develop a technique that compares the multidimensional distribution of all sources in the vicinity of each Gaia star to a reference distribution of random field stars obtained by extracting all sources in a region on the sky displaced 2'. This offset preserves the local field stellar density and magnitude distribution, allowing us to characterize the frequency of chance alignments. The resulting catalog with Bayesian probabilities >95% has a marginally higher match rate than current internal Gaia data release 2 (DR2) matches for most catalogs. However, a significant improvement is found with Pan-STARRS, where ∼99.8% of the sample within the Pan-STARRS footprint is recovered, as compared to a low ∼20.8% in Gaia DR2. Using these results, we train a Gaussian process regressor to calibrate two photometric metallicity relationships. For dwarfs of 3500 < T_eff < 5280 K, we use metallicity values of 4378 stars from the Apache Point Observatory Galactic Evolution Experiment and Hejazi et al. to calibrate the relationship, producing results with a 1σ precision of 0.12 dex and few systematic errors. We then indirectly infer the metallicity of 4018 stars with 2850 < T_eff < 3500 K, which are wide companions of primaries whose metallicities are estimated with our first regressor, to produce a relationship with a 1σ precision of 0.21 dex and significant systematic errors. Additional work is needed to better remove unresolved binaries from this second sample to reduce these systematic errors.

We present the confirmation of the eccentric warm giant planet TOI-201 b, first identified as a candidate in Transiting Exoplanet Survey Satellite photometry (Sectors 1–8, 10–13, and 27–28) and confirmed using ground-based photometry from Next Generation Transit Survey and radial velocities from FEROS, HARPS, CORALIE, and Minerva-Australis. TOI-201 b orbits a young (${0.87}_{-0.49}^{+0.46}\,\mathrm{Gyr}$) and bright (V = 9.07 mag) F-type star with a 52.9781 day period. The planet has a mass of ${0.42}_{-0.03}^{+0.05}\,{M}_{{\rm{J}}}$, a radius of ${1.008}_{-0.015}^{+0.012}\,{R}_{{\rm{J}}}$, and an orbital eccentricity of ${0.28}_{-0.09}^{+0.06};$ it appears to still be undergoing fairly rapid cooling, as expected given the youth of the host star. The star also shows long-term variability in both the radial velocities and several activity indicators, which we attribute to stellar activity. The discovery and characterization of warm giant planets such as TOI-201 b are important for constraining formation and evolution theories for giant planets.

Direct imaging of exoplanets is usually limited by quasi-static speckles. These uncorrected aberrations in a star's point-spread function (PSF) obscure faint companions and limit the sensitivity of high-contrast imaging instruments. Most current approaches to processing differential imaging sequences like angular differential imaging and spectral differential imaging produce a self-calibrating data set that is combined using a linear least-squares solution to minimize the noise. Due to temporal and chromatic evolution of a telescope's PSF, the best correlated reference images are usually the most contaminated by the planet, leading to self-subtraction and reducing the planet throughput. In this paper, we present an algorithm that directly optimizes the nonlinear equation for planet signal-to-noise ratio (S/N). This new algorithm does not require us to reject adjacent reference images and optimally balances noise reduction with self-subtraction. We then show how this algorithm can be applied to multiple images simultaneously for a further reduction in correlated noise, directly maximizing the S/N of the final combined image. Finally, we demonstrate the technique on an illustrative sequence of HR8799 using the new Julia-based Signal to Noise Analysis Pipeline. We show that S/N optimization can provide up to a 5× improvement in contrast close to the star. Applicable to both new and archival data, this technique will allow for the detection of fainter, lower mass, and closer-in companions, or achieve the same sensitivity with less telescope time.

The large time span and precise observational data of natural satellites is of great significance for updating their ephemerides and studying their dynamic characteristics. With the help of the new image-processing methods and the Gaia DR2 catalog, all CCD images of Triton taken with the 1.56 m telescope of Shanghai Astronomical Observatory during 2005–2009 were reanalyzed. The median filtering algorithm is used for image preprocessing to remove the influence of the halo of Neptune, and an upgraded modified moment, called the intensity-square-weighted centroiding method, is applied to determine the centroids of the stars and Triton. A total of 2299 positions of Triton were obtained, including 263 new observed positions and 2036 updated observed positions. Such five-year time span data with high precision will be very helpful to improve the orbit parameters of Triton.

As the earliest stage of planet formation, massive, optically thick, and gas-rich protoplanetary disks provide key insights into the physics of star and planet formation. When viewed edge-on, high-resolution images offer a unique opportunity to study both the radial and vertical structures of these disks and relate this to vertical settling, radial drift, grain growth, and changes in the midplane temperatures. In this work, we present multi-epoch Hubble Space Telescope and Keck scattered light images, and an Atacama Large Millimeter/submillimeter Array 1.3 mm continuum map for the remarkably flat edge-on protoplanetary disk SSTC2DJ163131.2–242627, a young solar-type star in ρ Ophiuchus. We model the 0.8 μm and 1.3 mm images in separate Markov Chain Monte Carlo (MCMC) runs to investigate the geometry and dust properties of the disk using the MCFOST radiative transfer code. In scattered light, we are sensitive to the smaller dust grains in the surface layers of the disk, while the submillimeter dust continuum observations probe larger grains closer to the disk midplane. An MCMC run combining both data sets using a covariance-based log-likelihood estimation was marginally successful, implying insufficient complexity in our disk model. The disk is well characterized by a flared disk model with an exponentially tapered outer edge viewed nearly edge-on, though some degree of dust settling is required to reproduce the vertically thin profile and lack of apparent flaring. A colder than expected disk midplane, evidence for dust settling, and residual radial substructures all point to a more complex radial density profile to be probed with future, higher-resolution observations.

We present high-resolution ¹²CO and ¹³CO 2–1 ALMA observations, as well as optical and near-infrared spectroscopy, of the highly inclined protoplanetary disk around SSTC2D J163131.2–242627. The spectral type we derive for the source is consistent with a 1.2 M_⊙ star inferred from the ALMA observations. Despite its massive circumstellar disk, we find little to no evidence for ongoing accretion on the star. The CO maps reveal a disk that is unusually compact along the vertical direction, consistent with its appearance in scattered light images. The gas disk extends about twice as far away as both the submillimeter continuum and the optical scattered light. CO is detected from two surface layers separated by a midplane region in which CO emission is suppressed, as expected from freeze-out in the cold midplane. We apply a modified version of the tomographically reconstructed distribution method presented by Dutrey et al. to derive the temperature structure of the disk. We find a temperature in the CO-emitting layers and the midplane of ∼33 K and ∼20 K at R < 200 au, respectively. Outside of R > 200 au, the disk's midplane temperature increases to ∼30 K, with a nearly vertically isothermal profile. The transition in CO temperature coincides with a dramatic reduction in the submicron and submillimeter emission from the disk. We interpret this as interstellar UV radiation providing an additional source of heating to the outer part of the disk.

The absence of planets interior to Mercury continues to puzzle terrestrial-planet formation models, particularly when contrasted with the relatively high derived occurrence rates of short-period planets around Sun-like stars. Recent work proposed that the majority of systems hosting hot super-Earths attain their orbital architectures through an epoch of dynamical instability after forming in quasi-stable, tightly packed configurations. Isotopic evidence seems to suggest that the formation of objects in the super-Earth-mass regime is unlikely to have occurred in the solar system as the terrestrial-forming disk is thought to have been significantly mass deprived starting around 2 Myr after the formation of calcium-aluminum-rich inclusions—a consequence of either Jupiter's growth or an intrinsic disk feature. Nevertheless, terrestrial-planet formation models and high-resolution investigations of planetesimal dynamics in the gas-disk phase occasionally find that quasi-stable protoplanets with mass comparable to that of Mars emerge in the vicinity of Mercury's modern orbit. In this paper, we investigate whether it is possible for a primordial configuration of such objects to be cataclysmically destroyed in a manner that leaves Mercury behind as the sole survivor without disturbing the other terrestrial worlds. We use numerical simulations to show that this scenario is plausible. In many cases, the surviving Mercury analog experiences a series of erosive impacts, thereby boosting its Fe/Si ratio. A caveat of our proposed genesis scenario for Mercury is that Venus typically experiences at least one late giant impact.

The eccentricity of a planet's orbit and the inclination of its orbital plane encode important information about its formation and history. However, exoplanets detected via direct imaging are often only observed over a very small fraction of their period, making it challenging to perform reliable physical inferences given wide, unconstrained posteriors. The aim of this project is to investigate biases (deviation of the median and mode of the posterior from the true values of orbital parameters, and the width and coverage of their credible intervals) in the estimation of orbital parameters of directly imaged exoplanets, particularly their eccentricities, and to define general guidelines to perform better estimations of uncertainty. For this, we constructed various orbits and generated mock data for each spanning ∼0.5% of the orbital period. We used the Orbits For The Impatient algorithm to compute orbit posteriors and compared those to the true values of the orbital parameters. We found that the inclination of the orbital plane is the parameter that most affects our estimations of eccentricity, with orbits that appear near edge on producing eccentricity distributions skewed away from the true values and often bimodal. We also identified a degeneracy between eccentricity and inclination that makes it difficult to distinguish posteriors of face-on, eccentric orbits and edge-on, circular orbits. For the exoplanet-imaging community, we propose practical recommendations, guidelines, and warnings relevant to orbit fitting.

We introduce the Automatic Learning for the Rapid Classification of Events (ALeRCE) broker, an astronomical alert broker designed to provide a rapid and self-consistent classification of large etendue telescope alert streams, such as that provided by the Zwicky Transient Facility (ZTF) and, in the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). ALeRCE is a Chilean-led broker run by an interdisciplinary team of astronomers and engineers working to become intermediaries between survey and follow-up facilities. ALeRCE uses a pipeline that includes the real-time ingestion, aggregation, cross-matching, machine-learning (ML) classification, and visualization of the ZTF alert stream. We use two classifiers: a stamp-based classifier, designed for rapid classification, and a light curve–based classifier, which uses the multiband flux evolution to achieve a more refined classification. We describe in detail our pipeline, data products, tools, and services, which are made public for the community (see https://alerce.science). Since we began operating our real-time ML classification of the ZTF alert stream in early 2019, we have grown a large community of active users around the globe. We describe our results to date, including the real-time processing of 1.5 × 10⁸ alerts, the stamp classification of 3.4 × 10⁷ objects, the light-curve classification of 1.1 × 10⁶ objects, the report of 6162 supernova candidates, and different experiments using LSST-like alert streams. Finally, we discuss the challenges ahead in going from a single stream of alerts such as ZTF to a multistream ecosystem dominated by LSST.

Turbulence is a prevalent phenomenon in the interstellar medium, and in particular, the environment at the centers of galaxies. For example, detailed observations of the Milky Way's Central Molecular Zone (CMZ) revealed that it has a complex and turbulent structure. Turbulence on galactic scales is often modeled using star formation and feedback. However, these effects do not appear to be sufficient for explaining the high-velocity dispersion observed in the CMZ, indicating that additional gas-stirring processes are likely to be operating. Here we introduce a proof-of-concept method to drive turbulence in gas that orbits under the influence of a galactic potential. Instead of relying on a particular physical mechanism, we have adopted a Fourier forcing module and have applied it using a smoothed particle hydrodynamics code. To test our method, we performed simulations of a simplistic model of the CMZ. Our turbulence injection method is capable of balancing the self-gravity of the gas, which allows us to run the simulations for long timescales and thereby follow the evolution of the CMZ. Our results show that turbulence induces a flocculent spiral pattern in our model, analogous to that found in galactic-scale simulations. Furthermore, we find that our turbulence injection method induces inward migration of gas, a result consistent with previous numerical simulations. We submit that this injection method is a promising new tool to simulate turbulence in galactic centers.

Recent discoveries of young exoplanets within their natal disks offer exciting opportunities to study ongoing planet formation. In particular, a planet's mass accretion rate can be constrained by observing the accretion-induced excess emission. So far, planetary accretion is only probed by the Hα line, which is then converted to a total accretion luminosity using correlations derived for stars. However, the majority of the accretion luminosity is expected to emerge from hydrogen continuum emission, and is best measured in the ultraviolet (UV). In this paper, we present HST/WFC3/UVIS F336W (UV) and F656N (Hα) high-contrast imaging observations of PDS 70. Applying a suite of novel observational techniques, we detect the planet PDS 70 b with signal-to-noise ratios of 5.3 and 7.8 in the F336W and F656N bands, respectively. This is the first time that an exoplanet has been directly imaged in the UV. Our observed Hα flux of PDS 70 b is higher by $3.5\sigma$ than the most recent published result. However, the light curve retrieved from our observations does not support greater than 30% variability in the planet's Hα emission in six epochs over a five month timescale. We estimate a mass accretion rate of $1.4\pm 0.2\times {10}^{-8}{M}_{\mathrm{Jup}}\,{\mathrm{yr}}^{-1}$. Hα accounts for 36% of the total accretion luminosity. Such a high proportion of energy released in line emission suggests efficient production of Hα emission in planetary accretion, and motivates using the Hα band for searches of accreting planets. These results demonstrate HST/WFC3/UVIS's excellent high-contrast imaging performance and highlight its potential for planet formation studies.

Wide-field small aperture optical telescopes are widely used in large-scale surveys currently and they have made great contributions in a number of astronomical applications. However, specific challenges arise owing to the defects caused by the optical system, and the image quality and reduction precision are negatively affected. An innovative method is proposed to address these challenges and achieve a high-precision source intensity estimation. In implementation, first a dedicated pipeline is developed to investigate the point-spread function (PSF) components from large amounts of images, using principal component analysis. Then the PSF model that reveals the actual characteristics of the optical system is constructed based on the evaluation. Last the equations for centroid and intensity estimation are constructed and the results are obtained. A trial of observations is performed with a wide-field small aperture telescope, and a large number of raw images, as well as simulated images, are acquired to test the efficiency of our method. The intensity measurement is performed with our method and other common algorithms, including the modified moment, Gaussian fitting, and SExtractor. Based on the comparison it is demonstrated that our proposed method outperforms the others. The results indicate that our method explores the limitations of such a system and additional gains can be achieved in wider applications.

We identify a set of planetary systems observed by Kepler that merit transit-timing variation (TTV) analysis given the orbital periods of transiting planets, the uncertainties for their transit times, and the number of transits observed during the Kepler mission. We confirm the planetary nature of four Kepler Objects of Interest within multicandidate systems. We forward-model each of the planetary systems identified to determine which systems are likely to yield mass constraints that may be significantly improved upon with follow-up transit observations. We find projected TTVs diverge by more than 90 minutes after 6000 days in 27 systems, including 22 planets with orbital periods exceeding 25 days. Such targets would benefit the most from additional transit-timing data. TTV follow-up could push exoplanet characterization to lower masses, at greater orbital periods and at cooler equilibrium temperatures than is currently possible from the Kepler data set alone. Combining TTVs and recently revised stellar parameters, we characterize an ensemble of homogeneously selected planets and identify planets in the Kepler field with large-enough estimated transmission annuli for atmospheric characterization with James Webb Space Telescope.

In this work, we present an analysis of 33,054 M-dwarf stars, located within 100 parsecs, via the Transiting Exoplanet Survey Satellite (TESS) full-frame images (FFIs) of observed sectors 1–5. We present a new pipeline called NEMESIS, developed to extract detrended photometry, and to perform transit searches of single-sector data in TESS FFIs. As many M-dwarfs are faint, and are not observed with a two-minute cadence by TESS, FFI transit surveys can provide an empirical validation of how many planets are missed, using the 30-minute cadence data. In this work, we detect 183 threshold crossing events, and present 29 candidate planets for sectors 1–5, 24 of which are new detections. Our sample contains orbital periods ranging from 1.25 to 6.84 days, and planetary radii from 1.26 to 5.31 R_⊕. With the addition of our new planet candidate detections, along with detections previously observed in sectors 1–5, we calculate an integrated occurrence rate of 2.49 ± 1.58 planets per star, for the period range ∈ [1, 9] days, and planet radius range ∈ [0.5,11] R_⊕. We project an estimated yield of 122 ± 11 transit detections of nearby M-dwarfs. Of our new candidates, 23 have signal-to-noise ratios >7, transmission spectroscopy metrics >38, and emission spectroscopy metrics >10. We present all of our data products for our planet candidates via the Filtergraph data visualization service, located at https://filtergraph.com/NEMESIS.

The B emission-line stars are rapid rotators that were probably spun up by mass and angular momentum accretion through mass transfer in an interacting binary. Mass transfer will strip the donor star of its envelope to create a small and hot subdwarf remnant. Here we report on Hubble Space Telescope/STIS far-ultraviolet spectroscopy of a sample of Be stars that reveals the presence of the hot sdO companion through the calculation of cross-correlation functions of the observed and model spectra. We clearly detect the spectral signature of the sdO star in 10 of the 13 stars in the sample, and the spectral signals indicate that the sdO stars are hot, relatively faint, and slowly rotating as predicted by models. A comparison of their temperatures and radii with evolutionary tracks indicates that the sdO stars occupy the relatively long-lived, He-core burning stage. Only 1 of the 10 detections was a known binary prior to this investigation, which emphasizes the difficulty of finding such Be+sdO binaries through optical spectroscopy. However, these results and others indicate that many Be stars probably host hot subdwarf companions.

Inconsistencies regarding the nature of globular cluster (GC) multiple population radial distributions is a matter for concern given their role in testing or validating cluster dynamical evolution modeling. In this study, we present a reanalysis of eight GC radial distributions using publicly available ground-based ugriz and UBVRI photometry; correcting for a systematic error identified in the literature. We detail the need for including and considering not only Kolmogorov–Smirnov (K-S) probabilities but critical K-S statistic values as well when drawing conclusions from radial distributions, as well as the impact of sample incompleteness. Revised cumulative radial distributions are presented, and the literature of each cluster is reviewed to provide a fuller picture of our results. We find that many multiple populations are not as segregated as once thought, and that there is a pressing need for better understanding of the spatial distributions of multiple populations in GCs.