Table of contents

Literature on the theory of exoplanet atmospheric disequilibrium chemistry is rich, although its observational counterpart has yet to emerge beyond the hints provided by a few targets in dedicated studies. We report results from an uniform data reduction and analysis for a catalog of 62 Hubble Space Telescope exoplanet transit spectra where we assess the atmospheric model preference for disequilibrium chemistry (i.e., water vapor is not the dominant absorption spectral signature) over thermal equilibrium chemistry in a comparative planetology context. Where model preference assessment is possible, we find that disequilibrium occurs in about half of the atmospheres, indicating that disequilibrium processes play an important role in the composition of exoplanet atmospheres. While very hot atmospheres, over 1800 K, prefer equilibrium chemistry, we find a clustering of preference for disequilibrium in the 1200–1800 K temperature range. We suggest that UV-augmented thermochemistry may play a significant role for those atmospheres.

The response of the antenna is a source of uncertainty in measurements with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES). We aim to validate the electromagnetic beam model of the low-band (50–100 MHz) dipole antenna with comparisons between models and against data. We find that simulations of a simplified model of the antenna over an infinite perfectly conducting ground plane are, with one exception, robust to changes in the numerical electromagnetic solver code or algorithm. For simulations of the antenna with the actual finite ground plane and realistic soil properties, we find that two out of three numerical solvers agree well. Applying our analysis pipeline to a simulated drift-scan observation from an early EDGES low-band instrument that had a 10 m × 10 m ground plane, we find residual levels after fitting and removing a five-term foreground model from the simulated data binned in local sidereal time (LST) average about 250 mK with ±40 mK variation between numerical solvers. A similar analysis of the primary 30 m × 30 m sawtooth ground plane reduced the LST-averaged residuals to about 90 mK with ±10 mK between the two viable solvers. More broadly we show that larger ground planes generally perform better than smaller ground planes. Simulated data have a power that is within 4% of real observations, a limitation of net accuracy of the sky and beam models. We observe that residual spectral structures after foreground model fits match qualitatively between simulated data and observations, suggesting that the frequency dependence of the beam is reasonably represented by the models. We find that a soil conductivity of 0.02 S m⁻¹ and relative permittivity of 3.5 yield good agreement between simulated spectra and observations. This is consistent with the soil properties reported by Sutinjo et al. for the Murchison Radio-astronomy Observatory, where EDGES is located.

Among the outer solar system minor planet orbits there is an observed gap in perihelion between roughly 50 and 65 au at eccentricities e ≳ 0.65. Through a suite of observational simulations, we show that the gap arises from two separate populations, the Extreme Trans-Neptunian Objects (ETNOs; perihelia q ≳ 40 au and semimajor axes a ≳ 150 au) and the Inner Oort Cloud objects (IOCs; q ≳ 65 au and a ≳ 250 au), and is very unlikely to result from a realistic single, continuous distribution of objects. We also explore the connection between the perihelion gap and a hypothetical distant giant planet, often referred to as Planet 9 or Planet X, using dynamical simulations. Some simulations containing Planet X produce the ETNOs, the IOCs, and the perihelion gap from a simple Kuiper-Belt-like initial particle distribution over the age of the solar system. The gap forms as particles scattered to high eccentricity by Neptune are captured into secular resonances with Planet X where they cross the gap and oscillate in perihelion and eccentricity over hundreds of kiloyears. Many of these objects reach a minimum perihelia in their oscillation cycle within the IOC region increasing the mean residence time of the IOC region by a factor of approximately five over the gap region. Our findings imply that, in the presence of a massive external perturber, objects within the perihelion gap will be discovered, but that they will be only ∼20% as numerous as the nearby IOC population (65 au ≲ q ≲ 100 au).

Themis family is one of the largest and oldest asteroid populations in the main belt. Water ice may widely exist on the parent body (24) Themis. In this work, we employ the Advanced Thermophysical Model as well as midinfrared measurements from NASA's Wide-field Infrared Survey Explorer to explore thermal parameters of 20 Themis family members. Here we show that the average thermal inertia and geometric albedo are 39.5 ± 26.0 J m⁻² s^−1/2 K⁻¹ and 0.067 ± 0.018, respectively. The family members have a relatively moderate roughness fraction on their surfaces. We find that the relatively low albedos of Themis members are consistent with the typical values of B-type and C-type asteroids. As mentioned, the Themis family bears a very low thermal inertia, which indicates a fine and mature regolith on their surfaces. The resemblance of thermal inertia and geometric albedo of Themis members may reveal their close connection in origin and evolution. In addition, we present the compared results of thermal parameters for several prominent families.

The results of speckle interferometric observations at the 4.1 m Southern Astrophysical Research Telescope in 2020, as well as earlier unpublished data, are given, totaling 1735 measurements of 1288 resolved pairs and nonresolutions of 1177 targets. We resolved for the first time 59 new pairs or subsystems in known binaries, mostly among nearby dwarf stars. This work continues our long-term speckle program. Its main goal is to monitor orbital motion of close binaries, including members of high-order hierarchies and Hipparcos pairs in the solar neighborhood. We also report observations of 892 members of young moving groups and associations, where we resolved 103 new pairs.

Globular clusters can form inside their host galaxies at high redshift when gas densities are higher and gas-rich mergers are common. They can also form inside lower-mass galaxies that have since been accreted and tidally disrupted, leaving their globular cluster complement bound to higher-mass halos. We argue that the age–metallicity–specific orbital energy relation in a galaxy's globular cluster system can be used to identify its origin. Gas-rich mergers should produce tightly bound systems in which metal-rich clusters are younger than metal-poor clusters. Globular clusters formed in massive disks and then scattered into a halo should have no relationship between age and specific orbital energy. Accreted globular clusters should produce weakly bound systems in which age and metallicity are correlated with eachother but inversely correlated with specific orbital energy. We use precise relative ages, self-consistent metallicities, and space-based proper motion-informed orbits to show that the Milky Way's metal-poor globular cluster system lies in a plane in age–metallicity–specific orbital energy space. We find that relatively young or metal-poor globular clusters are weakly bound to the Milky Way, while relatively old or metal-rich globular clusters are tightly bound to the Galaxy. While metal-rich globular clusters may be formed either in situ or ex situ, our results suggest that metal-poor clusters are formed outside of the Milky Way in now-disrupted dwarf galaxies. We predict that this relationship between age, metallicity, and specific orbital energy in a L* galaxy's globular cluster system is a natural outcome of galaxy formation in a ΛCDM universe.

Empirical trends in stellar X-ray and radio luminosities suggest that low-mass ultracool dwarfs (UCDs) should not produce significant radio emission. Defying these expectations, strong nonthermal emission has been observed in a few UCDs in the 1–10 GHz range, with a variable component often attributed to global aurorae and a steady component attributed to other processes, such as gyrosynchrotron emission. While both auroral and gyrosynchrotron emission peak near the critical frequency, only the latter radiation is expected to extend into millimeter wavelengths. We present Atacama Large Millimetre/Submillimeter Array (ALMA) 97.5 GHz and Very Large Array 33 GHz observations of a small survey of 5 UCDs. LP 349-25, LSR J1835+3259, and NLTT 33370 were detected at 97.5 GHz, while LP 423-31 and LP 415-20 resulted in nondetections at 33 GHz. A significant flare was observed in NLTT 33370, which reached a peak flux of 4880 ± 360 μJy, exceeding the quiescent flux by nearly an order of magnitude and lasting 20 s. These ALMA observations show bright 97.5 GHz emission with spectral indices ranging from α = 0.76 to α = −0.29, suggestive of optically thin gyrosynchrotron emission. If such emission traces magnetic reconnection events, then this could have consequences for both UCD magnetic models and the atmospheric stability of planets in orbit around them. Overall, our results provide confirmation that gyrosynchrotron radiation in radio-loud UCDs can remain detectable into the millimeter regime.

We report the direct imaging discovery of a low-mass companion to the nearby accelerating A star, HIP 109427, with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument coupled with the Microwave Kinetic Inductance Detector Exoplanet Camera (MEC) and CHARIS integral field spectrograph. CHARIS data reduced with reference star point spread function (PSF) subtraction yield 1.1–2.4 μm spectra. MEC reveals the companion in Y and J band at a comparable signal-to-noise ratio using stochastic speckle discrimination, with no PSF subtraction techniques. Combined with complementary follow-up L_p photometry from Keck/NIRC2, the SCExAO data favors a spectral type, effective temperature, and luminosity of M4–M5.5, 3000–3200 K, and ${\mathrm{log}}_{10}(L/{L}_{\odot })=-{2.28}_{-0.04}^{+0.04}$, respectively. Relative astrometry of HIP 109427 B from SCExAO/CHARIS and Keck/NIRC2, and complementary Gaia–Hipparcos absolute astrometry of the primary favor a semimajor axis of 6.55^+3.0_−0.48 au, an eccentricity of ${0.54}_{-0.15}^{+0.28}$, an inclination of ${66.7}_{-14}^{+8.5}$ degrees, and a dynamical mass of ${0.280}_{-0.059}^{+0.18}$M_⊙. This work shows the potential for extreme AO systems to utilize speckle statistics in addition to widely used postprocessing methods to directly image faint companions to nearby stars near the telescope diffraction limit.

We analyze existing measurements of [Fe/H] and [α/Fe] for individual red giant branch (RGB) stars in the Giant Stellar Stream (GSS) of M31 to determine whether spatial abundance gradients are present. These measurements were obtained from low- (R ∼ 3000) and moderate- (R ∼ 6000) resolution Keck/DEIMOS spectroscopy using spectral synthesis techniques as part of the Elemental Abundances in M31 survey. From a sample of 62 RGB stars spanning the GSS at 17, 22, and 33 projected kpc, we measure a [Fe/H] gradient of −0.018 ± 0.003 dex kpc⁻¹ and negligible [α/Fe] gradient with M31-centric radius. We investigate GSS abundance patterns in the outer halo using additional [Fe/H] and [α/Fe] measurements for six RGB stars located along the stream at 45 and 58 projected kpc. These abundances provide tentative evidence that the trends in [Fe/H] and [α/Fe] beyond 40 kpc in the GSS are consistent with those within 33 kpc. We also compare the GSS abundances to 65 RGB stars located along the possibly related Southeast (SE) shelf substructure at 12 and 18 projected kpc. The abundances of the GSS and SE shelf are consistent, supporting a common origin hypothesis, although this interpretation may be complicated by the presence of [Fe/H] gradients in the GSS. We discuss the abundance patterns in the context of photometric studies from the literature and explore implications for the properties of the GSS progenitor, suggesting that the high 〈[α/Fe]〉 of the GSS (+0.40 ± 0.05 dex) favors a major merger scenario for its formation.

The dynamical history of stars influences the formation and evolution of planets significantly. To explore the influence of dynamical history on the planet formation and evolution using observations, we assume stars that experienced significantly different dynamical histories tend to have different relative velocities. Utilizing the accurate Gaia–Kepler Stellar Properties Catalog, we select single main-sequence stars and divide these stars into three groups according to their relative velocities, i.e., high-V, medium-V, and low-V stars. After considering the known biases from Kepler data and adopting prior and posterior correction to minimize the influence of stellar properties on planet occurrence rate, we find that high-V stars have a lower occurrence rate of super-Earths and sub-Neptunes (1–4 R_⊕, P < 100 days) and a higher occurrence rate of sub-Earth (0.5–1 R_⊕, P < 30 days) than low-V stars. Additionally, high-V stars have a lower occurrence rate of hot Jupiter-sized planets (4–20 R_⊕, P < 10 days) and a slightly higher occurrence rate of warm or cold Jupiter-sized planets (4–20 R_⊕, 10 < P < 400 days). After investigating multiplicity and eccentricity, we find that high-V planet hosts prefer a higher fraction of multiplanet systems and lower average eccentricity, which are consistent with the eccentricity–multiplicity dichotomy of Kepler planetary systems. All of these statistical results favor the scenario that high-V stars with large relative velocity may experience fewer gravitational events, while low-V stars may be influenced by stellar clustering significantly.

Here we provide the most comprehensive determinations of the rest-frame UV luminosity function (LF) available to date with the Hubble Space Telescope (HST) at z ∼ 2–9. Essentially all of the noncluster extragalactic legacy fields are utilized, including the Hubble Ultra Deep Field, the Hubble Frontier Fields parallel fields, and all five CANDELS fields, for a total survey area of 1136 arcmin². Our determinations include galaxies at z ∼ 2–3 leveraging the deep HDUV, UVUDF, and ERS WFC3/UVIS observations available over an ∼150 arcmin² area in the GOODS-North and GOODS-South regions. All together, our collective samples include >24,000 sources, >2.3× larger than previous selections with HST. We identify 5766, 6332, 7240, 3449, 1066, 601, 246, and 33 sources at z ∼ 2, 3, 4, 5, 6, 7, 8, and 9, respectively. Combining our results with an earlier z ∼ 10 LF determination by Oesch et al., we quantify the evolution of the UV LF. Our results indicate that there is (1) a smooth flattening of the faint-end slope α from α ∼ −2.4 at z ∼ 10 to α ∼ −1.5 at z ∼ 2, (2) minimal evolution in the characteristic luminosity M* at z ≥ 2.5, and (3) a monotonic increase in the normalization ${\mathrm{log}}_{10}{\phi }^{* }$ from z ∼ 10 to 2, which can be well described by a simple second-order polynomial, consistent with an "accelerated" evolution scenario. We find that each of these trends (from z ∼ 10 to 2.5 at least) can be readily explained on the basis of the evolution of the halo mass function and a simple constant star formation efficiency model.

We report a detailed light-curve analysis of the Kepler target Kepler Input Catalog (KIC) 12602250. The results show that KIC 12602250 is a low-amplitude radial double-mode δ-Scuti star with amplitude modulation. The Fourier analysis of the long cadence data (i.e., Q0—Q17, spanning 1471 days) reveals that the variations of the light curve are dominated by the strongest mode with frequency F0 = 11.6141 d⁻¹, suggesting that KIC 12602250 is a δ-Scuti star. The other independent mode F1 = 14.9741 d⁻¹ is newly detected. The amplitude of the light variations of KIC 12602250 is ∼0.06 mag, which indicates that this is a low-amplitude δ-Scuti star; but the ratio of F0/F1 is estimated as 0.7756, which is typical of HADS, and a slow amplitude growth is detected in F1 and f₃, which could be due to stellar evolution, suggesting that KIC 12602250 could be a post-main-sequence δ Scuti that is crossing the instability strip for the first time.

One of the fundamental properties of an intermediate polar is the dynamical nature of the accretion flow as it encounters the white dwarf's (WD's) magnetosphere. Many works have presumed a dichotomy between disk-fed accretion, in which the WD accretes from a Keplerian disk, and stream-fed accretion, in which the matter stream from the donor star directly impacts the WD's magnetosphere without forming a disk. However, there is also a third, poorly understood regime in which the accretion flow consists of a torus of diamagnetic blobs that encircles the WD. This mode of accretion is expected to exist at mass-transfer rates below those observed during disk-fed accretion, but above those observed during pure stream-fed accretion. We invoke the diamagnetic-blob regime to explain the exceptional Transiting Exoplanet Survey Satellite light curve of the intermediate polar TX Col, which transitioned into and out of states of enhanced accretion during Cycles 1 and 3. Power-spectral analysis reveals that the accretion was principally stream fed. However, when the mass-transfer rate spiked, large-amplitude quasi-periodic oscillations (QPOs) abruptly appeared and dominated the light curve for weeks. The QPOs have two striking properties: they appear in a stream-fed geometry at elevated accretion rates, and they occur preferentially within a well-defined range of frequencies (∼10–25 cycle day⁻¹). We propose that during episodes of enhanced accretion, a torus of diamagnetic blobs forms near the binary's circularization radius and that the QPOs are beats between the white dwarf's spin frequency and unstable blob orbits within the WD's magnetosphere. We discuss how such a torus could be a critical step in producing an accretion disk in a formerly diskless system.

Measuring the obliquity distribution of stars hosting warm Jupiters may help us to understand the formation of close-orbiting gas giants. Few such measurements have been performed due to practical difficulties in scheduling observations of the relatively infrequent and long-duration transits of warm Jupiters. Here, we report a measurement of the Rossiter–McLaughlin effect for K2-232 b, a warm Jupiter on an 11.17 day orbit with an eccentricity of 0.26. The data were obtained with the Automated Planet Finder during two separate transits. The planet's orbit appears to be well aligned with the spin axis of the host star, with a projected spin–orbit angle of λ = −111 ± 66. Combined with the other available data, we find that high obliquities are almost exclusively associated with planets that either have an orbital separation greater than 10 stellar radii or orbit stars with effective temperatures hotter than 6000 K. This pattern suggests that the obliquities of the closest-orbiting giant planets around cooler stars have been damped by tidal effects.

Recently, it has been shown that rocky planets orbiting neutron stars can be habitable under plausible circumstances. If a distant, point-like source of visible light, such as a Sun-like main-sequence star or the gravitationally-lensed accretion disk of a supermassive black hole is present, possible temporal variations ${\rm{\Delta }}{\varepsilon }_{{\rm{p}}}\left(t\right)$ of the planet's axial tilt ε_p to the ecliptic plane should be included in the overall habitability budget since the obliquity determines the insolation at a given latitude on a body' s surface. I point out that, for rather generic initial spin–orbit initial configurations, general relativistic and classical spin variations induced by the post-Newtonian de Sitter and Lense–Thirring components of the field of the host neutron star and by its pull to the planetary oblateness ${J}_{2}^{{\rm{p}}}$ may induce huge and very fast variations of ε_p that would likely have an impact on the habitability of such worlds. In particular, for a planet's distance of, say, 0.005 au from a 1.4 M_⊙ neutron star corresponding to an orbital period P_b = 0.109 day, obliquity shifts Δε_p as large as ${\varepsilon }_{{\rm{p}}}^{\max }-{\varepsilon }_{{\rm{p}}}^{\min }\simeq 50^\circ \mbox{--}100^\circ$ over characteristic timescales as short as 10 days (${J}_{2}^{{\rm{p}}}$) to 3 Myr (Lense–Thirring) may occur for arbitrary orientations of the orbital and spin angular momenta L, S_ns, S_p of the planet-neutron star system. In view of this feature of their spins, I dub such hypothetical planets as "nethotrons."

We monitored a 3 deg² area toward Serpens Main in the Pan-STARRS1 r, i, and z bands from 2016 April to September. Light curves of more than 11,000 stars in each band were obtained, and 143 variables have been identified. Among those, 119 variables are new discoveries, while 24 variables were previously known. We present variability classes and periods of 99 stars. Of these, 81 are located in the upper giant branch, displaying long periods, while the remaining 18 variables are pre-main-sequence objects with short periods. We also identify eight eclipsing binary systems, including the known binary V0623 Ser, and derive their physical parameters. According to a clustering analysis of Gaia DR2 stars in the observed field, there are 10 variable members in Serpens Main, where six members have been classified as young stellar objects in previous studies. Here we provide color–magnitude and color–color diagrams for these variables. The color variability of most variables in the color–magnitude diagrams produces the expected displacements, while the movements of cluster members point in different directions; this behavior may be associated with accretion spots or circumstellar disks.

The results of 3611 intensified CCD observations of double stars, made with the 26 inch refractor of the U.S. Naval Observatory, are presented. Each observation of a system represents a combination of over 2000 short-exposure images. These observations are averaged into 1857 mean relative positions and range in separation from 0341 to 128644, with a median separation of 6533. Two systems have their orbits improved and one has its first determination. This is the 24th in this series of papers and covers the period 2018 May 24 through 2019 October 28.

We report the discovery of two planetary systems around comoving stars: TOI-2076  (TIC 27491137) and TOI-1807 (TIC 180695581). TOI-2076  is a nearby (41.9 pc) multiplanetary system orbiting a young (204 ± 50 Myr), bright (K = 7.115 in TIC v8.1) start. TOI-1807 hosts a single transiting planet and is similarly nearby (42.58 pc), similarly young (180 ± 40 Myr ), and bright. Both targets exhibit significant, periodic variability due to starspots, characteristic of their young ages. Using photometric data collected by TESS  we identify three transiting planets around TOI-2076  with radii of R_b = 3.3 ± 0.04 R_⊕, R_c = 4.4 ± 0.05 R_⊕, and R_d = 4.1 ± 0.07 R_⊕. Planet TOI-2076b  has a period of P_b = 10.356 days. For both TOI-2076c and d, TESS  observed only two transits, separated by a 2 yr interval in which no data were collected, preventing a unique period determination. A range of long periods (<17 days) are consistent with the data. We identify a short-period planet around TOI-1807 with a radius of R_b = 1.8 ± 0.04 R_⊕ and a period of P_b = 0.549 days. Their close proximity, and bright, cool host stars, and young ages make these planets excellent candidates for follow up. TOI-1807b  is one of the best-known small (R < 2 ${R}_{\oplus }$) planets for characterization via eclipse spectroscopy and phase curves with JWST. TOI-1807b  is the youngest ultra-short-period planet discovered to date, providing valuable constraints on formation timescales of short-period planets. Given the rarity of young planets, particularly in multiple-planet systems, these planets present an unprecedented opportunity to study and compare exoplanet formation, and young planet atmospheres, at a crucial transition age for formation theory.

Transit surveys have revealed a significant population of compact multiplanet systems, containing several sub-Neptune–mass planets on close-in, tightly-packed orbits. These systems are thought to have formed through a final phase of giant impacts, which would tend to leave systems close to the edge of stability. Here, we assess this hypothesis, comparing observed eccentricities in systems exhibiting transit-timing variations versus the maximum eccentricities compatible with long-term stability. We use the machine-learning classifier SPOCK (Tamayo et al.) to rapidly classify the stability of numerous initial configurations and hence determine these stability limits. While previous studies have argued that multiplanet systems are often maximally packed, in the sense that they could not host any additional planets, we find that the existing planets in these systems have measured eccentricities below the limits allowed by stability by a factor of 2–10. We compare these results against predictions from the giant-impact theory of planet formation, derived from both N-body integrations and theoretical expectations that, in the absence of dissipation, the orbits of such planets should be distributed uniformly throughout the phase space volume allowed by stability. We find that the observed systems have systematically lower eccentricities than this scenario predicts, with a median eccentricity about four times lower than predicted. This suggests that, if these systems formed through giant impacts, then some dissipation must occur to damp their eccentricities. This may occur through interactions with the natal gas disk or a leftover population of planetesimals, or over longer timescales through the coupling of tidal and secular processes.

Stellar kinematics is a powerful tool for understanding the formation process of stellar associations. Here, we present a kinematic study of the young stellar population in the Rosette nebula using recent Gaia data and high-resolution spectra. We first isolate member candidates using the published mid-infrared photometric data and the list of X-ray sources. A total of 403 stars with similar parallaxes and proper motions are finally selected as members. The spatial distribution of the members shows that this star-forming region is highly substructured. The young open cluster NGC 2244 in the center of the nebula has a pattern of radial expansion and rotation. We discuss its implication on the cluster formation, e.g., monolithic cold collapse or hierarchical assembly. On the other hand, we also investigate three groups located around the border of the H ii bubble. The western group seems to be spatially correlated with the adjacent gas structure, but their kinematics is not associated with that of the gas. The southern group does not show any systematic motion relative to NGC 2244. These two groups might be spontaneously formed in filaments of a turbulent cloud. The eastern group is spatially and kinematically associated with the gas pillar receding away from NGC 2244. This group might be formed by feedback from massive stars in NGC 2244. Our results suggest that the stellar population in the Rosette Nebula may form through three different processes: the expansion of stellar clusters, hierarchical star formation in turbulent clouds, and feedback-driven star formation.

The disintegrating planet candidate K2-22 b shows periodic and stochastic transits best explained by an escaping debris cloud. However, the mechanism that creates the debris cloud is unknown. The grain size of the debris as well as its sublimation rate can be helpful in understanding the environment that disintegrates the planet. Here, we present simultaneous photometry with the g band at 0.48 μm and K_S band at 2.1 μm using the Large Binocular Telescope. During an event with very low dust activity, we put a new upper limit on the size of the planet of 0.71 R_⊕ or 4500 km. We also detected a medium depth transit that can be used to constrain the dust particle sizes. We find that the median particle size must be larger than about 0.5–1.0 μm, depending on the composition of the debris. This leads to a high mass-loss rate of about 3 × 10⁸ kg s⁻¹, which is consistent with hydrodynamic escape models. If they are produced by some alternate mechanism such as explosive volcanism, it would require extraordinary geological activity. Combining our upper limits on the planet size with the high mass-loss rate, we find a lifetime of the planet of less than 370 Myr. This drops to just 21 Myr when adopting the 0.02 M_⊕ mass predicted from hydrodynamical models.

The presence of a stellar companion can place constraints on occurrence and orbital evolution of satellites orbiting exoplanets, i.e., exomoons. In this work we revise earlier orbital stability limits for retrograde orbits in the case of a three-body system consisting of a star, planet, and satellite. The revised limit reads ${a}_{\mathrm{sat}}^{\mathrm{crit}}\approx 0.668(1-1.236{e}_{{\rm{p}}})$ for e_p ≤ 0.8 in units of the Hill Radius and represents the lower critical orbit as a function of the planetary eccentricity e_p. A similar formula is determined for exomoons hosted by planets in binary star systems, where e_p is replaced with the components of free and forced eccentricity from secular orbit evolution theory. By exploring the dynamics of putative exomoons in α Centauri AB we find that the outer stability limit can be much less than half the Hill Radius due to oscillations in the planetary orbital eccentricity caused by the gravitational interaction with the binary star. We show, furthermore, how the resulting truncation of the outer stability limit can affect the outward tidal migration and potential observability of exomoons through transit-timing variations (TTVs). Typical TTV (rms) amplitudes induced by exomoons in binary systems are ≲10 minutes and appear more likely for planets orbiting the less massive stellar component.

MOA-2006-BLG-074 was selected as one of the most promising planetary candidates in a retrospective analysis of the MOA collaboration: its asymmetric high-magnification peak can be perfectly explained by a source passing across a central caustic deformed by a small planet. However, after a detailed analysis of the residuals, we have realized that a single lens and a source orbiting with a faint companion provides a more satisfactory explanation for all the observed deviations from a Paczynski curve and the only physically acceptable interpretation. Indeed the orbital motion of the source is constrained enough to allow a very good characterization of the binary source from the microlensing light curve. The case of MOA-2006-BLG-074 suggests that the so-called xallarap effect must be taken seriously in any attempts to obtain accurate planetary demographics from microlensing surveys.

We present Keck/NIRC2 adaptive optics imaging of planetary microlensing event MOA-2007-BLG-400 that resolves the lens star system from the source. We find that the MOA-2007-BLG-400L planetary system consists of a 1.71 ± 0.27M_Jup planet orbiting a 0.69 ± 0.04M_⊙ K-dwarf host star at a distance of 6.89 ± 0.77 kpc from the Sun. So, this planetary system probably resides in the Galactic bulge. The planet–host star projected separation is only weakly constrained due to the close-wide light-curve degeneracy; the 2σ projected separation ranges are 0.6–1.0 au and 4.7–7.7 au for close and wide solutions, respectively. This host mass is at the top end of the range of masses predicted by a standard Bayesian analysis. Our Keck follow-up program has now measured lens-source separations for six planetary microlensing events, and five of these six events have host star masses above the median prediction under the assumption that assumes that all stars have an equal chance of hosting planets detectable by microlensing. This suggests that more massive stars may be more likely to host planets of a fixed mass ratio that orbit near or beyond the snow line. These results also indicate the importance of host star mass measurements for exoplanets found by microlensing. The microlensing survey imaging data from NASA's Nancy Grace Roman Space Telescope (formerly WFIRST) mission will be doing mass measurements like this for a huge number of planetary events.

Barnard's star is among the most studied stars given its proximity to the Sun. It is often considered the radial velocity (RV) standard for fully convective stars due to its RV stability and equatorial decl. Recently, an $M\sin i=3.3{M}_{\oplus }$ super-Earth planet candidate with a 233 day orbital period was announced by Ribas et al. New observations from the near-infrared Habitable-zone Planet Finder (HPF) Doppler spectrometer do not show this planetary signal. We ran a suite of experiments on both the original data and a combined original + HPF data set. These experiments include model comparisons, periodogram analyses, and sampling sensitivity, all of which show the signal at the proposed period of 233 days is transitory in nature. The power in the signal is largely contained within 211 RVs that were taken within a 1000 day span of observing. Our preferred model of the system is one that features stellar activity without a planet. We propose that the candidate planetary signal is an alias of the 145 day rotation period. This result highlights the challenge of analyzing long-term, quasi-periodic activity signals over multiyear and multi-instrument observing campaigns.

We report the discovery of TOI-1444b, a 1.4 R_⊕ super-Earth on a 0.47 day orbit around a Sun-like star discovered by TESS. Precise radial velocities from Keck/HIRES confirmed the planet and constrained the mass to be 3.87 ± 0.71M_⊕. The RV data set also indicates a possible nontransiting, 16 day planet (11.8 ± 2.9M_⊕). We report a tentative detection of phase-curve variation and a secondary eclipse of TOI-1444b in the TESS bandpass. TOI-1444b joins the growing sample of 17 ultra-short-period planets (USPs) with well-measured masses and sizes, most of which are compatible with an Earth-like composition. We take this opportunity to examine the expanding sample of ultra-short-period planets (<2R_⊕) and contrast them with the newly discovered sub-day ultrahot Neptunes (>3R_⊕, >2000F_⊕ TOI-849 b, LTT9779 b, and K2-100). We find that (1) USPs have predominately Earth-like compositions with inferred iron core mass fractions of 0.32 ± 0.04 and have masses below the threshold of runaway accretion (∼10M_⊕), while ultrahot Neptunes are above the threshold and have H/He or other volatile envelopes. (2) USPs are almost always found in multi-planet systems consistent with a secular interaction formation scenario; ultrahot Neptunes (P_orb ≲1 day) tend to be "lonely," similar to longer-period hot Neptunes (P_orb1–10 days) and hot Jupiters. (3) USPs occur around solar-metallicity stars while hot Neptunes prefer higher metallicity hosts. (4) In all these respects, ultrahot Neptunes show more resemblance to hot Jupiters than the smaller USP planets, although ultrahot Neptunes are rarer than both USPs and hot Jupiters by 1–2 orders of magnitude.

In this work, we aimed to derive the gri-band period–luminosity (PL) and period–luminosity–color (PLC) relations for late-type contact binaries, for the first time, located in globular clusters, using the homogeneous light curves collected by the Zwicky Transient Factory (ZTF). We started with 79 contact binaries in 15 globular clusters, and retained 30 contact binaries in 10 globular clusters that have adequate numbers of data points in the ZTF light curves and are unaffected by blending. Magnitudes at mean and maximum light of these contact binaries were determined using a fourth-order Fourier expansion, while extinction corrections were done using the Bayerstar2019 3D reddening map together with adopting the homogeneous distances to their host globular clusters. After removing early-type and "anomaly" contact binaries, our derived gri-band PL and period–Wesenheit (PW) relations exhibited a much larger dispersion with large errors on the fitted coefficients. Nevertheless, the gr-band PL and PW relations based on this small sample of contact binaries in globular clusters were consistent with those based on a larger sample of nearby contact binaries. Good agreements of the PL and PW relations suggested both samples of contact binaries in the local Solar neighborhood and in the distant globular clusters can be combined and used to derive and calibrate the PL, PW, and PLC relations. The final derived gr-band PL, PW, and PLC relations were much improved over those based on the limited sample of contact binaries in the globular clusters.

We carry out a detailed photometric and kinematic study of the poorly studied sparse open clusters SAI 44 and SAI 45 using ground-based BVR_cI_c data supplemented by archival data from Gaia eDR3 and Pan-STARRS. The stellar memberships are determined using a statistical method based on Gaia eDR3 kinematic data, and we found 204 members in SAI 44 while only 74 members are identified in SAI 45. The average distances to SAI 44 and SAI 45 are calculated to be 3670 ± 184 and 1668 ± 47 pc. The logarithmic age of the clusters are determined to be 8.82 ± 0.10 and 9.07 ± 0.10 yr for SAI 44 and SAI 45, respectively. The color–magnitude diagram of SAI 45 hosts an extended main-sequence turnoff (eMSTO). The apparent age spread is found to be similar to the apparent age spread predicted on the basis of the age spread and cluster age relation predicted by rotation models. This indicates that eMSTO is a stellar evolution rather than star formation phenomenon in SAI 45. We conclude that eMSTO in SAI 45 is mainly caused by the different rotation rates of stars as the SYCLIST synthetic population with different rotation rates was able to reproduce the observed eMSTO, and stars in the red part of the eMSTO were preferentially concentrated in the inner region, which again hints at different rotations being the reason for the extension in the upper MS. This finding supports the theory attributing the origin of eMSTO to the different rotations of eMSTO stars. The mass function slopes are obtained as −2.24 ± 0.66 and −2.58 ± 3.20 in the mass rages 2.426–0.990 ${M}_{\odot }$ and 2.167–1.202 ${M}_{\odot }$ for SAI 44 and SAI 45, respectively. SAI 44 exhibits the signature of mass segregation while we found weak evidence of mass segregation in SAI 45 possibly due to tidal stripping. The dynamical relaxation times of these clusters indicate that both clusters are in a dynamically relaxed state. Using the AD-diagram method, the apex coordinates are found to be ($69\buildrel{\circ}\over{.} 79\pm 0\buildrel{\circ}\over{.} 11$, −$30\buildrel{\circ}\over{.} 82\pm 0\buildrel{\circ}\over{.} 15$) for SAI 44 and (−$56\buildrel{\circ}\over{.} 22\pm 0\buildrel{\circ}\over{.} 13$, −$56\buildrel{\circ}\over{.} 62\pm 0\buildrel{\circ}\over{.} 13$) for SAI 45. The average space velocity components of the clusters SAI 44 and SAI 45 are calculated in units of km s⁻¹ as (−15.14 ± 3.90, −19.43 ± 4.41, −20.85 ± 4.57) and (28.13 ± 5.30, −9.78 ± 3.13, −19.59 ± 4.43), respectively.

Many X-ray bright active galactic nuclei (AGNs) are predicted to follow an extended stage of obscured black hole growth. In support of this picture we examine the X-ray undetected AGNs in the COSMOS field and compare their host galaxies with X-ray bright AGNs. We examine galaxies with M_* > 10^9.5M_⊙ for the presence of AGNs at redshifts z = 0.5–3. We select AGNs in the infrared using Spitzer and Herschel detections and use color selection techniques to select AGNs within strongly star-forming hosts. We stack Chandra X-ray data of galaxies with an infrared (IR) detection but lacking an X-ray detection to obtain soft and hard fluxes, allowing us to measure the energetics of these AGNs. We find a clear correlation between X-ray luminosity and IR AGN luminosity in the stacked galaxies. We also find that X-ray undetected AGNs all lie on the main sequence—the tight correlation between the star formation rate and M_* that holds for the majority of galaxies, regardless of mass or redshift. This work demonstrates that there is a higher population of obscured AGNs than previously thought.

We present a photometric and spectroscopic study of HD 50526, an ellipsoidal binary member of the group Double Periodic Variable stars. Performing data mining in photometric surveys and conducting new spectroscopic observations with several spectrographs during 2008–2015, we obtained orbital and stellar parameters of the system. The radial velocities were analyzed with the genetic PIKAIA algorithm, whereas Doppler tomography maps for the Hα and Hβ lines were constructed with the Total Variation Minimization code. An optimized simplex algorithm was used to solve the inverse problem adjusting the light curve with the best stellar parameters for the system. We find an orbital period of 6701 ± 0001 and a long photometric cycle of 191 ± 2 days. We detected the spectral features of the coldest star and modeled it with a $\mathrm{log}g=2.79\pm 0.02\,\mathrm{dex}$ giant of mass 1.13 ± 0.02 M_⊙ and effective temperature 10500 ± 125 K. In addition, we determine a mass ratio q = 0.206 ± 0.033 and that the hot star is a B-type dwarf of mass 5.48 ± 0.02 M_⊙. The V-band orbital light curve can be modeled including the presence of an accretion disk around the hotter star. This fills the Roche lobe of the hotter star and has a radius 14.74 ± 0.02 R_⊙ and the temperature at the outer edge is 9400 K. Two bright spots located in the disk account for the global morphology of the light curve. The Doppler tomography maps of Hα and Hβ reveal complex structures of mass fluxes in the system.

We present a novel method of classifying Type Ia supernovae using convolutional neural networks, a neural network framework typically used for image recognition. Our model is trained on photometric information only, eliminating the need for accurate redshift data. Photometric data is preprocessed via 2D Gaussian process regression into two-dimensional images created from flux values at each location in wavelength-time space. These "flux heatmaps" of each supernova detection, along with "uncertainty heatmaps" of the Gaussian process uncertainty, constitute the data set for our model. This preprocessing step not only smooths over irregular sampling rates between filters but also allows SCONE to be independent of the filter set on which it was trained. Our model has achieved impressive performance without redshift on the in-distribution SNIa classification problem: 99.73 ± 0.26% test accuracy with no over/underfitting on a subset of supernovae from PLAsTiCC's unblinded test data set. We have also achieved 98.18 ± 0.3% test accuracy performing six-way classification of supernovae by type. The out-of-distribution performance does not fully match the in-distribution results, suggesting that the detailed characteristics of the training sample in comparison to the test sample have a big impact on the performance. We discuss the implication and directions for future work. All of the data processing and model code developed for this paper can be found in the SCONE software package located at github.com/helenqu/scone.

Using the Karl G. Jansky Very Large Array (VLA), we have conducted a survey for 22 GHz, 6_1,6–5_2,3 H₂O masers toward the Serpens South region. The masers were also observed with the Very Long Baseline Array following the VLA detections. We detect for the first time H₂O masers in the Serpens South region that are found to be associated to three Class 0–Class I objects, including the two brightest protostars in the Serpens South cluster, known as CARMA-6 and CARMA-7. We also detect H₂O masers associated to a source with no outflow or jet features. We suggest that this source is most probably a background asymptotic giant branch star projected in the direction of Serpens South. The spatial distribution of the emission spots suggest that the masers in the three Class 0–Class I objects emerge very close to the protostars and are likely excited in shocks driven by the interaction between a protostellar jet and the circumstellar material. Based on the comparison of the distributions of bolometric luminosity of sources hosting 22 GHz H₂O masers and 162 young stellar objects covered by our observations, we identify a limit of L_Bol ≈ 10L_⊙ for a source to host water masers. However, the maser emission shows strong variability in both intensity and velocity spread, and therefore masers associated to lower-luminosity sources may have been missed by our observations. We also report 11 new sources with radio continuum emission at 22 GHz.

Kepler planets (including super-Earths and sub-Neptunes, from 1–4 Earth radii) are likely formed before the gaseous protoplanetary disks have dissipated, as are the Jovian planets. If the metal content in these disks resembles that in the host stars, one might expect Kepler planets to occur more frequently, and to be more massive, around metal-rich stars. Contrary to these expectations, we find that the radii of Kepler planets (a proxy for mass) are independent of host metallicity. Previous claims that larger planets prefer more metal-rich stars can be adequately explained by the combined facts that more massive stars tend to host bigger planets, and that more massive stars are also more metal-rich in the Kepler sample. We interpret this independence as that the mass of a Kepler planet is not determined by the availability of solids, but is instead regulated by an as yet unknown process. Moreover, we find that the occurrence rates of Kepler planets rise only weakly with stellar metallicity, a trend that is further flattened when the influence of close stellar binaries is accounted for. We explain this weak dependence, in contrast to the strong dependence exhibited by Jovian planets, using a phenomenological model, wherein the masses of protoplanetary disks have a much larger spread than the spread in stellar metallicity, and wherein the formation of Jovian planets requires disks that contain some 5 times more solids than that needed to form Kepler planets. This model predicts that stars more metal-poor than half-solar should rarely host any Kepler planets.

The near-Sun comet C/2019 Y4 (ATLAS) is the first member of a long-period comet group observed to disintegrate well before perihelion. Here we present our investigation into this disintegration event using images obtained in a three-day Hubble Space Telescope campaign. We identify two fragment clusters produced by the initial disintegration event, corresponding to fragments C/2019 Y4-A and C/2019 Y4-B identified in ground-based data. These two clusters started with similar integrated brightness but exhibit different evolutionary behavior. C/2019 Y4-A was much shorter-lived compared to C/2019 Y4-B and showed signs of significant mass loss and changes in size distribution throughout the three-day campaign. The cause of the initial fragmentation is undetermined by the limited evidence but crudely compatible with either the spin-up disruption of the nucleus or runaway sublimation of subsurface supervolatile ices, either of which would lead to the release of a large amount of gas as inferred from the significant bluing of the comet observed shortly before its disintegration. Gas can only be produced by the sublimation of volatile ices, which must have survived at least one perihelion passage at a perihelion distance of q = 0.25 au. We speculate that Comet ATLAS is derived from the ice-rich interior of a nonuniform, kilometer-wide progenitor that split during its previous perihelion. This suggests that comets down to a few kilometers in diameter can still possess complex, nonuniform interiors that can protect ices against intense solar heating.

We report the discovery of a new emission-line object, named SPH 4-South = (GAIA EDR3 5616553300192230272), toward the dark cloud LDN 1667. This object came to our attention after inspecting public images that show a faint diffuse nebula a few arcseconds south of SPH 4, an emission-line object previously classified as a T Tauri star. We present high-resolution spectra and analyzed JHK photometry of SPH 4 and SPH 4-South and new narrowband and archival broadband images of these objects. A comparison of the spectra of SPH 4 and SPH 4-South with high-resolution ones of DG Cir and R Mon strongly suggests that SPH 4 and SPH 4-South are Herbig Ae/Be stars. The classification of SPH 4-South is further supported by using a k-NN algorithm to its position in an H−K versus J−H color–color diagram. Both stars are detected in the four WISE bands and the WISE colors allow us to classify SPH 4 as a Class I and SPH 4-South as a Class II source. We also show that the faint nebula is most probably associated with SPH 4-South. Using published results on LDN 1667 and the Gaia Early Data Release 3, we conclude that SPH 4 is a member of LDN 1667. The case of SPH 4-South is not clear because the determination of its distance and proper motion could be affected by the nebulosity around the star, although membership of SPH 4-South to LDN 1667 cannot be ruled out.

We present a method of selecting quasars up to redshift ≈6 with random forests, a supervised machine-learning method, applied to Pan-STARRS1 and WISE data. We find that, thanks to the increasing set of known quasars, we can assemble a training set that enables supervised machine-learning algorithms to become a competitive alternative to other methods up to this redshift. We present a candidate set for the redshift range 4.8–6.3, which includes the region around z = 5.5 where selecting quasars is difficult due to their photometric similarity to red and brown dwarfs. We demonstrate that, under our survey restrictions, we can reach a high completeness (66% ± 7% below redshift 5.6/${83}_{-9}^{+6} \%$ above redshift 5.6) while maintaining a high selection efficiency (${78}_{-8}^{+10} \%$/${94}_{-8}^{+5} \%$). Our selection efficiency is estimated via a novel method based on the different distributions of quasars and contaminants on the sky. The final catalog of 515 candidates includes 225 known quasars. We predict the candidate catalog to contain additional ${148}_{-33}^{+41}$ new quasars below redshift 5.6 and ${45}_{-8}^{+5}$ above, and we make the catalog publicly available. Spectroscopic follow-up observations of 37 candidates led us to discover 20 new high redshift quasars (18 at 4.6 ≤ z ≤ 5.5, 2 z ∼ 5.7). These observations are consistent with our predictions on efficiency. We argue that random forests can lead to higher completeness because our candidate set contains a number of objects that would be rejected by common color cuts, including one of the newly discovered redshift 5.7 quasars.

The present-day envelope of gaseous planets is a relic of how these giant planets originated and evolved. Measuring their elemental composition therefore presents a powerful opportunity to answer long-standing questions regarding planet formation. Obtaining precise observational constraints on the elemental inventory of giant exoplanets has, however, remained challenging owing to the limited simultaneous wavelength coverage of current space-based instruments. Here, we present thermal emission observations of the nontransiting hot Jupiter τ Boo b using the new wide wavelength coverage (0.95–2.50 μm) and high spectral resolution (R = 70,000) CFHT/SPIRou spectrograph. By combining a total of 20 hr of SPIRou data obtained over five nights in a full atmospheric retrieval framework designed for high-resolution data, we constrain the abundances of all the major oxygen- and carbon-bearing molecules and recover a noninverted temperature structure using a new free-shape, nonparametric temperature–pressure profile retrieval approach. We find a volume mixing ratio of log(CO) = −${2.46}_{-0.29}^{+0.25}$ and a highly depleted water abundance of less than 0.0072 times the expected value for a solar composition envelope. Combined with upper limits on the abundances of CH₄, CO₂, HCN, TiO, and C₂H₂, this results in a gas-phase C/H ratio of ${5.85}_{-2.82}^{+4.44}$ × solar, consistent with the value of Jupiter, and an envelope C/O ratio robustly greater than 0.60, even when taking into account the oxygen that may be sequestered out of the gas phase. Combined, the inferred supersolar C/H, O/H, and C/O ratios on τ Boo b support a formation scenario beyond the water snowline in a disk enriched in CO owing to pebble drift.

We report and analyze updated molecular abundances in 20 comets obtained by employing modern data reduction procedures and molecular models. Using box and scatter plots, we examine how the different molecular species are distributed among the comet population, while by means of pie charts, we investigate the relative proportions of these molecular species in each comet. We compare these results with the orbital parameters of the selected comets to identify trends related to their dynamical history. With these analyses, we tentatively identify at least three chemical classes based mainly on relative abundances of CO, CH₃OH, CH₄, C₂H₆, HCN, and NH₃. The combination of relative abundances and orbital parameters is then compared with recent chemical models of planetary system formation. This approach may help in investigating the origins and evolution of the material in cometary nuclei. Among other aspects, we underline the need to increase our sample size, especially for hypervolatiles (i.e., CH₄ and CO) in Jupiter family comets.

We present high-angular-resolution imaging observations of 517 host stars of TESS exoplanet candidates using the 'Alopeke and Zorro speckle cameras at Gemini North and South. The sample consists mainly of bright F, G, K stars at distances of less than 500 pc. Our speckle observations span angular resolutions of ∼20 mas out to 12, yielding spatial resolutions of <10–500 au for most stars, and our contrast limits can detect companion stars 5–9 mag fainter than the primary at optical wavelengths. We detect 102 close stellar companions and determine the separation, magnitude difference, mass ratio, and estimated orbital period for each system. Our observations of exoplanet host star binaries reveal that they have wider separations than field binaries, with a mean orbital semimajor axis near 100 au. Other imaging studies have suggested this dearth of very closely separated binaries in systems which host exoplanets, but incompleteness at small separations makes it difficult to disentangle unobserved companions from a true lack of companions. With our improved angular resolution and sensitivity, we confirm that this lack of close exoplanet host binaries is indeed real. We also search for a correlation between planetary orbital radii versus binary star separation; but, given the very short orbital periods of the TESS planets, we do not find any clear trend. We do note that in exoplanet systems containing binary host stars, there is an observational bias against detecting Earth-size planet transits due to transit depth dilution caused by the companion star.

Blazars represent the dominant class of associated γ-ray sources detected by the Fermi Large Area Telescope (LAT). However, in the more recent release of the Fourth Fermi-LAT Point Source Catalog (4FGL), ∼25% of the sources associated with lower-energy counterparts show a multifrequency behavior similar to that of blazars, but lacks an optical spectroscopic confirmation of their nature and are therefore classified as Blazar Candidates of Uncertain Type (BCUs). A particularly challenging task in blazar studies is to classify these BCUs and, when possible to estimate their redshifts, in particular for BL Lac objects, characterized by almost featureless optical spectra with only weak emission lines. Continuing our study of blazars with Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) optical spectral data, we perform an extensive search for optical spectra available in the LAMOST Data Release 6 archive. Our aim is confirming the blazar nature of BCUs and to test if new data can allow us to get a redshift estimate for BL Lac objects that lack measurement, as well as to search for and discover changing-look blazars. We selected sources out of the 4FGL catalog, the list of targets from our follow-up spectroscopic campaign of unidentified and/or unassociated γ-ray sources, and the Roma-BZCAT multifrequency catalog of blazars, finding a total of 42 sources with available LAMOST DR6 spectra. We confirmed the blazar-like nature of four blazar candidates . For the remaining 37 sources we confirm their previous classification.

We report an analysis of the planetary microlensing event OGLE-2018-BLG-1185, which was observed by a large number of ground-based telescopes and by the Spitzer Space Telescope. The ground-based light curve indicates a low planet–host star mass ratio of q = (6.9 ± 0.2) × 10⁻⁵, which is near the peak of the wide-orbit exoplanet mass-ratio distribution. We estimate the host star and planet masses with a Bayesian analysis using the measured angular Einstein radius under the assumption that stars of all masses have an equal probability of hosting the planet. The flux variation observed by Spitzer is marginal, but still places a constraint on the microlens parallax. Imposing a conservative constraint that this flux variation should be Δf_Spz < 4 instrumental flux units yields a host mass of ${M}_{\mathrm{host}}={0.37}_{-0.21}^{+0.35}\ {M}_{\odot }$ and a planet mass of ${m}_{{\rm{p}}}={8.4}_{-4.7}^{+7.9}\ {M}_{\oplus }$. A Bayesian analysis including the full parallax constraint from Spitzer suggests smaller host star and planet masses of ${M}_{\mathrm{host}}={0.091}_{-0.018}^{+0.064}\ {M}_{\odot }$ and ${m}_{{\rm{p}}}={2.1}_{-0.4}^{+1.5}\ {M}_{\oplus }$, respectively. Future high-resolution imaging observations with the Hubble Space Telescope or Extremely Large Telescope could distinguish between these two scenarios and help reveal the planetary system properties in more detail.

We describe for the first time in scientific literature the Planetary Ephemeris Program (PEP), an open-source general-purpose astrometric data-analysis program. We discuss, in particular, the implementation of pulsar timing analysis, which was recently upgraded in PEP to handle more options. This implementation was done independently of other pulsar programs, with minor exceptions that we discuss. We illustrate the implementation of this capability by comparing the postfit residuals from the analyses of time-of-arrival observations by both PEP and Tempo2. The comparison shows substantial agreement: 22 ns rms differences for 1065 pulse time-of-arrival measurements for the millisecond pulsar in a binary system, PSR J1909-3744 (pulse period 2.947108 ms; full width half maximum of pulse 43 μs), for epochs in the interval from 2002 December to 2011 February.

Studies of close-in planets orbiting M dwarfs have suggested that the M-dwarf radius valley may be well explained by distinct formation timescales between enveloped terrestrials and rocky planets that form at late times in a gas-depleted environment. This scenario is at odds with the picture that close-in rocky planets form with a primordial gaseous envelope that is subsequently stripped away by some thermally driven mass-loss process. These two physical scenarios make unique predictions of the rocky/enveloped transition's dependence on orbital separation such that studying the compositions of planets within the M-dwarf radius valley may be able to establish the dominant physics. Here, we present the discovery of one such keystone planet: the ultra-short-period planet TOI-1634 b (P = 0.989 days, $F=121{F}_{\oplus }$, ${r}_{p}={1.790}_{-0.081}^{+0.080}$R_⊕) orbiting a nearby M2 dwarf (K_s = 8.7, R_s = 0.450 R_⊙, M_s = 0.502 M_⊙) and whose size and orbital period sit within the M-dwarf radius valley. We confirm the TESS-discovered planet candidate using extensive ground-based follow-up campaigns, including a set of 32 precise radial velocity measurements from HARPS-N. We measure a planetary mass of ${4.91}_{-0.70}^{+0.68}$M_⊕, which makes TOI-1634 b inconsistent with an Earth-like composition at $5.9\sigma$ and thus requires either an extended gaseous envelope, a large volatile-rich layer, or a rocky composition that is not dominated by iron and silicates to explain its mass and radius. The discovery that the bulk composition of TOI-1634 b is inconsistent with that of Earth supports the gas-depleted formation mechanism to explain the emergence of the radius valley around M dwarfs with ${M}_{s}\lesssim 0.5$M_⊙.

The Extragalactic Distance Database (EDD) was created as a repository for high-quality, redshift-independent distances. A key component of EDD is the Color–Magnitude Diagrams/Tip of the Red Giant Branch (CMDs/TRGB) catalog, which provides information on the stellar content of nearby galaxies observed with the Hubble Space Telescope (HST). Here we provide a decadal update to this catalog, which has now doubled in size to over 500 galaxies. We highlight the additions to our data reduction and analysis techniques and provide examples of the science that has been made possible with this large data set. We find the TRGB to be a reliable measure for distance, and we aim to extend its distance coverage with HST to every galaxy within 10 Mpc. In the near future, the combination of the James Webb Space Telescope and the Nancy Grace Roman Space Telescope will dramatically increase the number of targets within our grasp.