Table of contents for issue 1, volume 162, The Astronomical Journal

The modern search for extraterrestrial intelligence began with the seminal publications of Cocconi & Morrison and Schwartz & Townes, who proposed searching for narrowband signals in the radio spectrum and optical laser pulses. Over the last six decades, more than 100 dedicated search programs have targeted these wavelengths, all with null results. All of these campaigns searched for classical communications, that is, for a significant number of photons above a noise threshold, with the assumption of a pattern encoded in time and/or frequency space. I argue that future searches should also target quantum communications. They are preferred over classical communications with regard to security and information efficiency, and they would have escaped detection in all previous searches. The measurement of Fock state photons or squeezed light would indicate the artificiality of a signal. I show that quantum coherence is feasible over interstellar distances and explain for the first time how astronomers can search for quantum transmissions sent by ETI to Earth using commercially available telescopes and receiver equipment.

The star RIK 90 is a member of a group of recently discovered young T Tauri stars showing complex periodic modulations of their light output, which have been difficult to explain. As one of the brightest, and having a relatively short period, RIK 90 is of particular interest. Thirty-three hours of contemporaneous BVR_CI_C photometry of the star are presented in this paper. An independent secondary periodicity was removed and the remainder of the data phased with the primary 0.25955 day period and smoothed to reveal the detailed structure of the light curves. The ratios of the amplitude in B to those in VR_CI_C are 1.2, 1.5, and 1.8. It is found that a collection of five cool starspots with temperatures in the range 2230–2550 K provides an excellent match to almost all observations, small-scale features in the I_C-band light curve being the notable exception. Starspot modeling of K2 observations of RIK 90 was also attempted. The morphology of the K2 light curve is entirely different from that discussed in this paper but could be accurately modeled by four dark starspots. In order to put these results into context, 15 other distinct K2 light curves belonging to 10 similarly variable T Tauri stars were also studied. Excellent fits were obtained to some of these with models consisting of large cool starspots. The majority of the remaining light curves could be modeled with a combination of hot and cold starspots.

Modern terrestrial-planet formation models are highly successful at consistently generating planets with masses and orbits analogous to those of Earth and Venus. In stark contrast to classic theoretical predictions and inferred demographics of multiplanet systems of rocky exoplanets, the mass (≳10) and orbital period (≳2) ratios between Venus and Earth and the neighboring Mercury and Mars are not common outcomes in numerically generated systems. While viable solutions to the small-Mars problem are abundant in the literature, Mercury’s peculiar origin remains rather mysterious. In this paper, we investigate the possibility that Mercury formed in a mass-depleted, inner region of the terrestrial disk (a < 0.5 au). This regime is often ignored in terrestrial-planet formation models because of the high computational cost of resolving hundreds of short-period objects over ∼100 Myr timescales. By testing multiple disk profiles and mass distributions, we identify several promising sets of initial conditions that lead to remarkably successful analog systems. In particular, our most successful simulations consider moderate total masses of Mercury-forming material (0.1–0.25 Earth masses). While larger initial masses tend to yield disproportionate Mercury analogs, smaller values often inhibit the planets’ formation as the entire region of material is easily accreted by Venus. Additionally, we find that shallow surface density profiles and larger inventories of small planetesimals moderately improve the likelihood of adequately reproducing Mercury.

Cometary outbursts offer a valuable window into the composition of comet nuclei with their forceful ejection of dust and volatiles in explosive events, revealing the interior components of the comet. Understanding how different types of outbursts influence the dust properties and volatile abundances, to better interpret what signatures can be attributed to primordial composition and what features are the result of processing, is an important task best undertaken with a multi-instrument approach. The European Space Agency Rosetta mission to 67P/Churyumov–Gerasimenko carried a suite of instruments capable of carrying out this task in the near-nucleus coma with unprecedented spatial and spectral resolution. In this work, we discuss two outbursts that occurred 2015 November 7 and were observed by three instruments on board: the Alice ultraviolet spectrograph, the Visual Infrared and Thermal Imaging Spectrometer, and the Optical, Spectroscopic, and Infrared Remote Imaging System. Together, the observations show that mixed gas and dust outbursts can have different spectral signatures representative of their initiating mechanisms, with the first outburst showing indicators of a cliff collapse origin and the second more representative of fresh volatiles being exposed via a deepening fracture. This analysis opens up the possibility of remote spectral classification of cometary outbursts with future work.

The Alice ultraviolet spectrograph on board the Rosetta orbiter provided the first near-nucleus ultraviolet observations of a cometary coma from arrival at comet 67P/Churyumov-Gerasimenko in 2014 August through 2016 September. The characterization of atomic and molecular emissions in the coma revealed the unexpected contribution of dissociative electron impact emission at large heliocentric distances and during some outbursts. This mechanism also proved useful for compositional analysis, and Alice observed many cases that suggested elevated levels of the supervolatile O₂, identifiable in part to their emissions resulting from dissociative electron impact. In this paper, we present the first two-dimensional UV maps constructed from Alice observations of atomic emission from 67P during an increase in cometary activity on 2015 November 7–8. Comparisons to observations of the background coma and an earlier collimated jet are used to describe possible changes to the near-nucleus coma and plasma. To verify the mapping method and place the Alice observations in context, comparisons to images derived from the MIRO and VIRTIS-H instruments are made. The spectra and maps we present show an increase in dissociative electron impact emission and an O₂/H₂O ratio of ∼0.3 for the activity; these characteristics have been previously identified with cometary outbursts seen in Alice data. Further, UV maps following the increases in activity show the spatial extent and emission variation experienced by the near-nucleus coma, informing future UV observations of comets that lack the same spatial resolution.

We present a new method to measure the vertical aerosol optical depth (VAOD) during clear nights using a wide-field imager—a CCD camera with a photographic lens on an equatorial mount. A series of 30 s exposures taken at different altitudes above the horizon can be used to measure the VAOD with a precision better than 0.008 optical depths within a few minutes. Such a measurement does not produce any light and is thus suitable for use at sites where other astronomical instruments are located. The precision of the VAOD measurement depends on the laboratory calibration of spectral properties of the system and of the response of the camera electronics to varying illumination levels, as well as careful considerations of the details of stellar photometry and modeling of the dependence of the measured stellar flux, the star color, and the position within the field of view. The results obtained with robotic setups at the future sites of the Cerenkov Telescope Array show good internal consistency and agreement with simultaneous measurements from a Sun/Moon photometer located at the same sites.

We report the discovery and characterization of seven transiting exoplanets from the HATNet survey. The planets, which are hot Jupiters and Saturns transiting bright Sun-like stars, include: HAT-P-58b (with mass M_p = 0.37 M_J, radius R_p = 1.33 R_J, and orbital period P = 4.0138 days), HAT-P-59b (M_p = 1.54 M_J, R_p = 1.12 R_J, P = 4.1420 days), HAT-P-60b (M_p = 0.57 M_J, R_p = 1.63 R_J, P = 4.7948 days), HAT-P-61b (M_p = 1.06 M_J, R_p = 0.90 R_J, P = 1.9023 days), HAT-P-62b (M_p = 0.76 M_J, R_p = 1.07 R_J, P = 2.6453 days), HAT-P-63b (M_p = 0.61 M_J, R_p = 1.12 R_J, P = 3.3777 days), and HAT-P-64b (M_p = 0.58 M_J, R_p = 1.70 R_J, P = 4.0072 days). The typical errors on these quantities are 0.06 M_J, 0.03 R_J, and 0.2 s, respectively. We also provide accurate stellar parameters for each of the host stars. With V = 9.710 ± 0.050 mag, HAT-P-60 is an especially bright transiting planet host, and an excellent target for additional follow-up observations. With R_p = 1.703 ± 0.070 R_J, HAT-P-64b is a highly inflated hot Jupiter around a star nearing the end of its main-sequence lifetime, and is among the largest known planets. Five of the seven systems have long-cadence observations by TESS which are included in the analysis. Of particular note is HAT-P-59 (TOI-1826.01) which is within the northern continuous viewing zone of the TESS mission, and HAT-P-60, which is the TESS candidate TOI-1580.01.

The Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effect governs the spin evolution of small asteroids. The axial component of YORP, which alters the rotation rate of the asteroid, is mostly independent of its thermal inertia, while the obliquity component is very sensitive to the thermal model of the asteroid. Here, we develop a semi-analytic theory for the obliquity component of YORP. We integrate an approximate thermal model over the surface of an asteroid, and find an analytic expression for the obliquity component in terms of two YORP coefficients. This approach allows us to investigate the overall evolution of asteroid rotation state, and to generalize the results previously obtained in the case of zero thermal inertia. The proposed theory also explains how a nonzero obliquity component of YORP originates even for a symmetric asteroid, due to its finite thermal inertia. In many cases, this causes equatorial planes of asteroids to align with their orbital planes. The studied nontrivial behavior of YORP as a function of thermal model allows for a new kind of rotational equilibria, which can have important evolutionary consequences for asteroids.

Future direct imaging missions will primarily observe planets that have been previously detected, mostly via the radial velocity (RV) technique, to characterize planetary atmospheres. In the meantime, direct imaging may discover new planets within existing planetary systems that have bright enough reflected flux but insufficient signals for other methods to detect. Here we investigate the parameter space within which planets are unlikely to be detected by RV in the near future due to precision limitations but could be discovered through reflected light with future direct imaging missions. We use the HD 134987 system as a working example, combine RV and direct imaging detection limit curves in the same parameter space through various assumptions, and insert a fictitious planet into the system while ensuring that it lies between the RV and imaging detection limits. Planet validity tested through dynamical simulations and retrieval tests revealed that the planet could indeed be detected by imaging while remaining hidden from RV surveys. Direct imaging retrieval was carried out using starshade simulations for two mission concepts: the Starshade Rendezvous Probe, which could be coupled with the Nancy Grace Roman Space Telescope, and the Habitable Exoplanet Observatory. This method is applicable to all other systems and high-contrast direct imaging instruments and could help inform future imaging observations and data analysis on the discovery of new exoplanets.

To fully characterize the atmospheres, or lack thereof, of terrestrial exoplanets, we must include the high-energy environments provided by their host stars. The nearby mid-M dwarf LHS 3844 hosts a terrestrial world that lacks a substantial atmosphere. We present a time-series UV spectrum of LHS 3844 from 1131 to 3215 Å captured by HST/COS. We detect one flare in the FUV that has an absolute energy of 8.96 ± 0.77 × 10²⁸ erg and an equivalent duration of 355 ± 31 s. We extract the flare and quiescent UV spectra separately. For each spectrum, we estimate the Lyα flux using correlations between UV line strengths. We use Swift-XRT to place an upper limit on the soft X-ray flux and construct a differential emission model to estimate flux that is obscured by the interstellar medium. We compare the differential emission model flux estimates in the XUV to other methods that rely on scaling from the Lyα, Si iv, and N v lines in the UV. The XUV, FUV, and NUV flux of LHS 3844 relative to its bolometric luminosity is log₁₀(L_band/L_Bol) =−3.65, −4.16, and −4.48, respectively, for the quiescent state. These values agree with trends in high-energy flux as a function of stellar effective temperature found by the MUSCLES survey for a sample of early-M dwarfs. Many of the most spectroscopically accessible terrestrial exoplanets orbit inactive mid-to-late M dwarfs like LHS 3844. Measurements of M dwarf high-energy spectra are preferable for exoplanet characterization but are not always possible. The spectrum of LHS 3844 is a useful proxy for the current radiation environment for these worlds.

We report the flare activity of Wolf 359, the fifth closest star to the Sun and a candidate exoplanet-hosting M dwarf. The star was a target of the Kepler/K2 mission and was observed by the EDEN project, a global network of 1–2 m class telescopes for detection and characterization of rocky exoplanets in the habitable zones of late-M dwarfs within 50 light year from the solar system. In the combination of the archived K2 data and our EDEN observations, a total of 872 flares have been detected, 861 with the K2 (860 in the short-cadence and 18 in the long-cadence data, with 17 long-cadence events having short-cadence counterparts) and 11 with EDEN. Wolf 359 has relatively strong flare activity even among flaring M dwarfs, in terms of the flare activity indicator (FA) defined as the integrated flare energy relative to the total stellar bolometric energy, where FA = ∑E_f/∫L_boldt ∼ 8.93 × 10⁻⁵ for the long-cadence flares, whereas for K2 short cadence and EDEN flares, the FA values are somewhat larger, FA ≈ 6.67 × 10⁻⁴ and FA ≈ 5.25 × 10⁻⁴, respectively. Such a level of activity, in accordance with the rotation period (P_rot), suggests the star to be in the saturation phase. The size of the starspots is estimated to be at least 1.87% ± 0.59% of the projected disk area of Wolf 359. We find no correlation of FA with the stellar rotational phase. Our analysis indicates a flare frequency distribution in a power-law form of with α = 2.13 ± 0.14, equivalent to an occurrence rate of flares E_f ≥ 10³¹ erg about once per day and of superflares with E_f ≥ 10³³ erg approximately 10 times per year. These superflares may impact the habitability of system in multiple ways, the details of which are topics for future investigations.

The extensive time span of modern radial velocity surveys has made the discovery of long-period substellar companions more common in recent years; however, measuring the true masses of these objects remains challenging. Astrometry from the Gaia mission is expected to provide mass measurements for many of these long-period companions, but these data are not yet available. However, combining proper-motion data from Gaia DR2 and the earlier Hipparcos mission makes it possible to measure true masses of substellar companions in favorable cases. In this work, we combine radial velocities with Hipparcos–Gaia astrometry to measure the true masses of two recently discovered long-period substellar companion candidates, HD 92987 B and HD 221420 b. In both cases, we find that the true masses are significantly higher than implied by radial velocities alone. A 2087 ± 19 m s⁻¹ astrometric signal reveals that HD 92987 B is not close to its 17 M_J minimum mass but is instead a 0.2562 ± 0.0045 M_⊙ star viewed at a near-polar orbital inclination, whereas the 22.9 ± 2.2 M_J HD 221420 b can be plausibly interpreted as a high-mass “superplanet” or a low-mass brown dwarf. With semimajor axes of ∼10 au, both companions are interesting targets for direct imaging, and HD 221420 b in particular would be a benchmark metal-rich substellar object if it proves possible to directly detect. Our results demonstrate the power of Hipparcos–Gaia astrometry for studying long-period planet and brown dwarf candidates discovered from radial velocity surveys.

High precision CCD observations of six totally eclipsing contact binaries were presented and analyzed. It is found that only one target is an A-type contact binary (V429 Cam), while the others are W-type contact ones. By analyzing the times of light minima, we discovered that two of them exhibit secular period increase while three manifest long-term period decrease. For V1033 Her, a cyclic variation superimposed on the long-term increase was discovered. By comparing the Gaia distances with those calculated by the absolute parameters of 173 contact binaries, we found that the Gaia distance can be applied to estimate the absolute parameters for most contact binaries. The absolute parameters of our six targets were estimated using their Gaia distances. The evolutionary status of contact binaries was studied, and we found that the A- and W-subtype contact binaries may have different formation channels. The relationship between the spectroscopic and photometric mass ratios for 101 contact binaries was presented. It was discovered that the photometric mass ratios are in good agreement with the spectroscopic ones for almost all of the totally eclipsing systems, which is corresponding to the results derived by Pribulla et al. and Terrell & Wilson.

Alpha Centauri A is the closest solar-type star to the Sun and offers the best opportunity to find and ultimately to characterize an Earth-sized planet located in its habitable zone. Here, we describe initial results from an Atacama Large Millimeter/submillimeter Array (ALMA) program to search for planets in the α Cen AB system using differential astrometry at millimeter wavelengths. Our initial results include new absolute astrometric measurements of the proper motion, orbital motion and parallax of the α Cen system. These lead to an improved knowledge of the physical properties of both α Cen A and B. Our estimates of ALMA’s relative astrometric precision suggest that we will ultimately be sensitive to planets of a few tens of Earth mass in orbits from 1 to 3 au, where stable orbits are thought to exist.

We analyze the very short Einstein timescale (t_E ≃ 7 hr) event KMT-2019-BLG-2073. Making use of the pronounced finite-source effects generated by the clump giant source, we measure the Einstein radius θ_E ≃ 4.8 μas and so infer a mass , where π_rel is the lens-source relative parallax. We find no significant evidence for a host of this planetary-mass object, though one could be present at sufficiently wide separation. If so, it would be detectable after about 10 yr. This is the fourth isolated microlens with a measured Einstein radius θ_E < 10 μas, which we argue is a useful threshold for a “likely free-floating planet (FFP)” candidate. We outline a new approach to constructing a homogeneous sample of giant-star finite-source/point-lens (FSPL) events, within which the subsample of FFP candidates can be statistically analyzed. We illustrate this approach using 2019 KMTNet data and show that there appears to be a large θ_E gap between the two FFP candidates and the 11 other FSPL events. We argue that such sharp features are more identifiable in a sample selected on θ_E compared to the traditional approach of identifying candidates based on short t_E.

Over the course of the last two decades, traditional models of planet formation have been repeatedly challenged by the emerging census of extrasolar planets. Key among them is the orbital architecture problem: while standard models of orbital migration predict resonant orbits for short-period objects, most planets do not appear to lie in orbital resonances. Here, we show that the four-planet system Kepler-221, not previously recognized to have active orbital resonances, has a three-body commensurability relation unique within the Kepler sample. Using a suite of numerical experiments as well as a perturbative analysis, we demonstrate that this system likely began as a resonant chain and proceeded to undergo large-scale divergence away from resonance, under the action of tidal dissipation. Our results further indicate that obliquity tides, driven by a secular spin–orbit resonance and mutual inclination, are an excellent candidate for driving this orbital divergence, and that the high tidal luminosity may also explain the anomalous size of planet b, which lies within the Fulton radius gap.

We show that the perturbation at the peak of the light curve of microlensing event KMT-2019-BLG-0371 is explained by a model with a mass ratio between the host star and planet of q ∼ 0.08. Due to the short event duration (t_E ∼ 6.5 days), the secondary object in this system could potentially be a massive giant planet. A Bayesian analysis shows that the system most likely consists of a host star with a mass and a massive giant planet with a mass . However, the interpretation of the secondary as a planet (i.e., as having M_p < 13M_Jup) rests entirely on the Bayesian analysis. Motivated by this event, we conduct an investigation to determine which constraints meaningfully affect Bayesian analyses for microlensing events. We find that the masses inferred from such a Bayesian analysis are determined almost entirely by the measured value of θ_E and are relatively insensitive to other factors such as the direction of the event (ℓ, b), the lens–source relative proper motion μ_rel, or the specific Galactic model prior.

The most challenging limitation in transit photometry arises from the noises in the photometric signal. In particular, the ground-based telescopes are heavily affected by the noise due to perturbation in the Earth’s atmosphere. Use of telescopes with large apertures can improve the photometric signal-to-noise ratio to a great extent. However, detecting a transit signal out of a noisy light curve of the host star and precisely estimating the transit parameters call for various noise reduction techniques. Here, we present multiband transit photometric follow-up observations of five hot Jupiters e.g., HAT-P-30 b, HAT-P-54 b, WASP-43 b, TrES-3 b, and XO-2 N b, using the 2 m Himalayan Chandra Telescope at the Indian Astronomical Observatory, Hanle, and the 1.3 m J. C. Bhattacharya Telescope at the Vainu Bappu Observatory, Kavalur. Our critical noise treatment approach includes techniques such as wavelet denoising and Gaussian process regression, which effectively reduce both time-correlated and time-uncorrelated noise components from our transit light curves. In addition to these techniques, use of our state-of-the-art model algorithms have allowed us to estimate the physical properties of the target exoplanets with a better accuracy and precision compared to the previous studies.

We investigate the origins of the photometrically very red (VR) and less red (LR) trans-Neptunian objects (TNOs). We first reanalyze the data set of Marsset et al. and find that in addition to the known color–inclination correlation in hot TNOs, a similar trend exists for color–eccentricity. We show that VR TNOs are sharply constrained to eccentricities <0.42 and inclinations <21°, leading to a paucity of VR scattered disk and distant mean motion resonance objects. We then interpret these findings using N-body simulations accounting for Neptune’s outward migration into a massless particles disk and find that these observations are best reproduced with an LR-to-VR color transition line between ∼38 and 42 au in the primordial disk, separating the objects’ formation locations. For an initial surface density profile (Σ ∝ 1/r²), a color transition around 38 au is needed to explain the high abundance of VR plutinos, but it creates too many VR scattered disk objects, while a transition line around 42 au seems to better reproduce the scattered disk colors but creates virtually no VR plutinos. Our simulations furthermore show that the rarity of VR particles at high eccentricity is possibly due to the absence of sweeping higher-order MMRs, and secular resonances, beyond 42 au. Inspecting individual populations, we show that the majority of VR SDOs originate as objects trapped in Neptune’s second- and third-order MMRs. These then evolve due to diffusion, scattering, Kozai–Lidov cycles, and secular resonances into their current orbits. Future unbiased color surveys are crucial to better constrain the TNOs dynamical origins.

We answer the question of whether long-period comet C/1961 T1 (Seki) is the parent body of a meteor shower observed in Earth’s atmosphere. We model discrete parts of a meteoroid stream associated with this comet and follow the orbital evolution of its particles via a numerical integration of their orbits. The gravitational perturbations of major planets and the Poynting–Robertson effect, which influences the motion of the particles, are taken into account. We predict two showers which are related to C/1961 T1. One of these showers was identified with the December ρ-Virginids, #502. This identification confirms an earlier suggested relationship between the shower and the comet. The second predicted shower was identified with the γ-Sagittariids, #657. The radiant areas of both showers are symmetrical in respect to the apex of Earth’s heliocentric motion.

We make use of individual (epoch) detection data from the Pan-STARRS “3π” survey for 2863 optical ICRF3 counterparts in the five wavelength bands g, r, i, z, and y, published as part of the Data Release 2. A dedicated method based on the Functional Principal Component Analysis is developed for these sparse and irregularly sampled data. With certain regularization and normalization constraints, it allows us to obtain uniform and compatible estimates of the variability amplitudes and average magnitudes between the passbands and objects. We find that the starting assumption of affinity of the light curves for a given object at different wavelengths is violated for several percent of the sample. The distributions of rms variability amplitudes are strongly skewed toward small values, peaking at ∼0.1 mag with tails stretching to 2 mag. Statistically, the lowest variability is found for the r band and the largest for the reddest y band. A small “brighter-redder” effect is present, with amplitudes in y greater than amplitudes in g in 57% of the sample. The variability versus redshift dependence shows a strong decline with z toward redshift 3, which we interpret as the time dilation of the dominant time frequencies. The colors of radio-loud ICRF3 quasars are correlated with redshift in a complicated, wavy pattern governed by the emergence of brightest emission lines within the five passbands.

We study dynamical systems that admit action-angle variables at leading order, which are subject to nearly resonant perturbations. If the frequencies characterizing the unperturbed system are not in resonance, the long-term dynamical evolution may be integrated by orbit-averaging over the high-frequency angles, thereby evolving the orbit-averaged effect of the perturbations. It is well known that such integrators may be constructed via a canonical transformation, which eliminates the high-frequency variables from the orbit-averaged quantities. An example of this algorithm in celestial mechanics is the von Zeipel transformation. However, if the perturbations are inside or close to a resonance, i.e., the frequencies of the unperturbed system are commensurate; these canonical transformations are subject to divergences. We introduce a canonical transformation that eliminates the high-frequency phase variables in the Hamiltonian without encountering divergences. This leads to a well-behaved symplectic integrator. We demonstrate the algorithm through two examples: a resonantly perturbed harmonic oscillator and the gravitational three-body problem in mean motion resonance.

The proposition that life can spread from one planetary system to another (interstellar panspermia) has a long history, but this hypothesis is difficult to test through observations. We develop a mathematical model that takes parameters such as the microbial survival lifetime, the stellar velocity dispersion, and the dispersion of ejecta into account in order to assess the prospects for detecting interstellar panspermia. We show that the correlations between pairs of life-bearing planetary systems (embodied in the pair-distribution function from statistics) may serve as an effective diagnostic of interstellar panspermia, provided that the velocity dispersion of ejecta is greater than the stellar dispersion. We provide heuristic estimates of the model parameters for various astrophysical environments and conclude that open clusters and globular clusters appear to represent the best targets for assessing the viability of interstellar panspermia.

To reveal the origins of diffuse Hα emissions observed around the Herbig star MWC 1080, we have performed a high-resolution near-infrared (NIR) spectroscopic observation using the Immersion Grating Infrared Spectrograph. In the NIR H and K bands, we detected various emission lines (six hydrogen Brackett lines, seven H₂ lines, and an [Fe ii] line) and compared their spatial locations with the optical (Hα and [S ii]) and radio (¹³CO and CS) line maps. The shock-induced H₂ and [Fe ii] lines indicate the presence of multiple outflows, consisting of at least three associated young stars in this region. The kinematics of H₂ and [Fe ii] near the northeast (NE) cavity edge supports the idea that the NE main outflow from MWC 1080A is the blueshifted one with a low inclination angle. The H₂ and [Fe ii] lines near the southeast molecular region newly reveal that additional highly blueshifted outflows originate from other young stars. The fluorescent H₂ lines were found to trace photodissociation regions formed on the cylindrical surfaces of the main outflow cavity, which are expanding outward with a velocity of about 10–15 km s⁻¹. For the Hα emission, we identify its components associated with two stellar outflows and two young stars in addition to the dominant component of MWC 1080A scattered by dust. We also report a few faint Hα features located ∼0.4 pc away in the southwest direction from MWC 1080A, which lie near the axes of the NE main outflow and one of the newly identified outflows.

In 2018, Solar Cycle 24 entered into a solar minimum phase. During this period, 11 million zenithal night sky brightness (NSB) data were collected at different dark sites around the planet, including astronomical observatories and natural protected areas, with identical broadband Telescope Encoder and Sky Sensor photometers (based on the Unihedron Sky Quality Meter TLS237 sensor). A detailed observational review of the multiple effects that contribute to the NSB measurement has been conducted with optimal filters designed to avoid brightening effects by the Sun, the Moon, clouds, and other astronomical sources (the Galaxy and zodiacal light). The natural NSB has been calculated from the percentiles for 44 different photometers by applying these new filters. The pristine night sky was measured to change with an amplitude of 0.1 mag/arcsec² in all the photometers, which is suggested to be due to NSB variations on scales of up to months and to be compatible with semiannual oscillations. We report the systematic observation of short-time variations in NSB on the vast majority of the nights and find these to be related to airglow events forming above the mesosphere.

Cross-correlation analyses of high-resolution spectroscopic data have recently shown great success in directly detecting planetary signals and enabling the characterization of their atmospheres. One such technique aims to observe a system at multiple epochs and combine the measured planetary radial velocities from each epoch into a measurement of the planetary Keplerian orbital velocity K_p, constituting a direct detection of the planetary signal. Recent work has shown that in few-epoch (∼5) data sets, unintended structure can arise at a high level, obscuring the planetary detection. In this work, we look to simulations to examine whether there are ways to reduce this structured noise in few-epoch data sets by careful planning of observations. The choice of observation date allows observers to select the primary (stellar) velocity through a set systemic velocity and chosen barycentric velocity and the planetary orbital phase so we focus on the effects of these two parameters. We find that epochs taken when the primary velocity is near zero and the stellar lines remain relatively fixed to the telluric rest-frame greatly reduce the level of structured noise and allow for much stronger planetary detections, on average more than twice the significance of detections made with epochs using randomly selected primary velocities. Following these results, we recommend that observers looking to build up high-resolution multi-epoch data sets target nights when their system has a near-zero primary velocity.

Trans-Neptunian Objects (TNOs) in the scattered disk with 50 ≲ a ≲ 100 au are thought to cluster near Neptune’s n:1 resonances (e.g., 3:1, 4:1, and so on). While these objects spend lengthy periods of time at large heliocentric distances, if their perihelia remain less than around 40 au, their dynamical evolution is still largely coupled to Neptune’s. Conversely, around a dozen extreme TNOs with a ≳ 250 au and detached perihelia seem to exist in a regime where they are too distant to be affected by the giant planets and too close for their dynamics to be governed by external forces. Recent work suggests that the apparent alignment of these orbits in physical space is a signature of gravitational shepherding by a distant massive planet. In this paper, we investigate the evolution of TNOs in each of Neptune’s n:1 resonances between the 3:1 and 14:1. We conclude that both resonant and nonresonant objects beyond the 12:1 near ∼157 au are removed rather efficiently via perturbations from the hypothetical Planet Nine. Additionally, we uncover a population of simulated TNOs with a ≲ 100 au, 40 ≲ q ≲ 45 au, and low inclinations that experience episodes of resonant interactions with both Neptune and Planet Nine. Finally, we simulate the evolution of observed objects with a > 100 au and identify several TNOs that are potentially locked in n:1 resonances with Neptune, including the most distant known resonant candidates, 2014 JW₈₀ and 2014 OS₃₉₄,which appear to be in the 10:1 and 11:1 resonances, respectively. Our results suggest that the detection of similar remote objects might provide a useful constraint on hypotheses invoking the existence of additional distant planets.

Gaps in protoplanetary disks have long been hailed as signposts of planet formation. However, a direct link between exoplanets and disks remains hard to identify. We present a large sample study of ALMA disk surveys of nearby star-forming regions to disentangle this connection. All disks are classified as either structured (transition, ring, extended) or nonstructured (compact) disks. Although low-resolution observations may not identify large-scale substructure, we assume that an extended disk must contain substructure from a dust evolution argument. A comparison across ages reveals that structured disks retain high dust masses up to at least 10 Myr, whereas the dust mass of compact, nonstructured disks decreases over time. This can be understood if the dust mass evolves primarily by radial drift, unless drift is prevented by pressure bumps. We identify a stellar mass dependence of the fraction of structured disks. We propose a scenario linking this dependence with that of giant exoplanet occurrence rates. We show that there are enough exoplanets to account for the observed disk structures if transitional disks are created by exoplanets more massive than Jupiter and ring disks by exoplanets more massive than Neptune, under the assumption that most of those planets eventually migrate inwards. On the other hand, the known anticorrelation between transiting super-Earths and stellar mass implies those planets must form in the disks without observed structure, consistent with formation through pebble accretion in drift-dominated disks. These findings support an evolutionary scenario where the early formation of giant planets determines the disk’s dust evolution and its observational appearance.

The landing site of China’s Chang’e-4 (CE-4) probe is located on the mare basalts on the floor of the Von Kármán crater on the lunar far side. The Von Kármán crater is inside the ancient and highly cratered South Pole–Aitken basin, which has experienced complex emplacement sequences of both near and distant ejecta materials. These issues complicated the interpretation of the CE-4 surface in situ measurements of the visible and near-infrared spectrometer and the lunar penetrating radar (LPR) onboard Yutu-2 rover. To evaluate the sources and amounts of all principle foreign materials at the CE-4 landing site, we thoroughly examine the ejecta delivered by crater-forming events that occurred later than the formation of the mare basalts at the CE-4 landing site. We found a total of 16 craters that may have delivered ejecta thicker than 10 cm level superposed on the mare basalts at the CE-4 landing site. Crater Finsen, Von Kármán L, Von Kármán L′, and Maksutov are the top four major foreign material sources, and each of them contributed ejecta thicker than 1 m. Our surveys confirm that the ejecta from Finsen crater are the most dominant foreign materials in the uppermost few meters at the CE-4 landing site and the total impact ejecta deposited upon the mare basalts at the landing site is estimated to be thinner than 30 m. We found that the estimations from Pike’s model are the most consistent with the Yutu-2 LPR observations.

Simulations of exoplanet albedo profiles are key to planning and interpreting future direct imaging observations. In this paper we demonstrate the use of the Planetary Spectrum Generator (PSG) to produce simulations of reflected light exoplanet spectra. We use PSG to examine multiple issues relevant to all models of directly imaged exoplanet spectra and to produce sample spectra of the bright, nearby exoplanet υ Andromedae d (υ And d), a potential direct imaging target for next-generation facilities. We introduce a new, fast, and accurate subsampling technique that enables calculations of disk-integrated spectra one order of magnitude faster than Chebyshev–Gauss sampling for moderate- to high-resolution sampling. Using this method and a first-principles-derived atmosphere for υ And d, we simulate phase-dependent spectra for a variety of different potential atmospheric configurations. The simulated spectra for υ And d include versions with different haze and cloud properties. Based on our combined analysis of this planet’s orbital parameters, phase- and illumination-appropriate model spectra, and realistic instrument noise parameters, we find that υ And d is a potentially favorable direct imaging and spectroscopy target for the Coronagraph Instrument (CGI) on the Nancy Grace Roman Space Telescope. When a noise model corresponding to the Roman CGI SPC spectroscopy mode is included, PSG predicts the time required to reach a signal-to-noise ratio of 10 of the simulated spectra in both the central wavelength bin of the Roman CGI SPC spectroscopy mode (R = 50 spectrum) and of the Band 1 HLC imaging mode is approximately 400 and less than 40 hr, respectively. We also discuss potential pathways to extricating information about the planet and its atmosphere with future observations and find that Roman observations may be able to bound the interior temperature of the planet.

We observed the compact group of galaxies NGC 6845 with GMOS-S imaging and long-slit spectroscopy. The group itself is comprised of two spiral galaxies and two S0 galaxies, and shows signs of gravitational interaction among them. While the spiral galaxies have been already well studied, the two S0 galaxies have received little attention. We investigate the morphology and stellar kinematics of the two S0 galaxies in detail in order to derive the value for their mass and subsequently incorporate them into the mass model for the group. It was found that the internal kinematics of the largest S0 (NGC 6845 C) is mostly regular, while the smallest S0 (NGC 6845 D) shows some hint of interaction, most likely with NGC 6845 C and NGC 6845 A. The new determination of the dynamical mass of the group implies that the group mass has to be 10 times the combined mass of the individual members according to their internal kinematics. The results point to an early stage of the group interaction and also to the possibility that it has still not reached virial equilibrium. From the galaxy evolution viewpoint, this group represents an example of the case of the transformation of Spirals into S0 galaxies in dense environments.

Rotation is a critical physical process operating in the formation of Am stars. There is a strong correlation between low-velocity rotation and chemical peculiarity. However, the existence of many non-chemical-peculiar slow rotators challenges the understanding of Am stars. The purpose of our work is to search for low-velocity rotating non-chemical-peculiar A-type stars and Am stars and to make a comparative analysis. In this paper, we pick out a sample from the LAMOST–Kepler project, including 21 Am stars, 125 non-chemical-peculiar slow rotators, and 53 non-chemical-peculiar fast rotators. We calculate the rotational frequencies through a periodic change of light curves caused by inhomogeneous stellar surfaces and then obtain the rotational velocities. For slow rotators, the age of Am stars is statistically younger than that of non-chemical-peculiar stars in the same temperature bin. In the comparison of the period, the average amplitude, and stellar mass of Am and non-chemical-peculiar stars, we discover that there is no difference in the photometric variability between Am and non-chemical-peculiar stars, which implies similar inhomogeneities on the surfaces. The average amplitude of non-chemical-peculiar stars has a downward trend with the increase of effective temperature and stellar mass, which is consistent with the theoretical prediction caused by weak dynamo-generated magnetic fields in A-type stars. In addition, we confirm four non-chemical-peculiar stars that have flares by checking field-of-view images, pixel images, and pixel-level light curves.

A line of sight toward the Galactic Center (GC) offers the largest number of potentially habitable systems of any direction in the sky. The Breakthrough Listen program is undertaking the most sensitive and deepest targeted SETI surveys toward the GC. Here, we outline our observing strategies with Robert C. Byrd Green Bank Telescope (GBT) and Parkes telescope to conduct 600 hr of deep observations across 0.7–93 GHz. We report preliminary results from our survey for extraterrestrial intelligence (ETI) beacons across 1–8 GHz with 7.0 and 11.2 hr of observations with Parkes and GBT, respectively. With our narrowband drifting signal search, we were able to place meaningful constraints on ETI transmitters across 1–4 GHz and 3.9–8 GHz with EIRP limits of ≥4 × 10¹⁸ W among 60 million stars and ≥5 × 10¹⁷ W among half a million stars, respectively. For the first time, we were able to constrain the existence of artificially dispersed transient signals across 3.9–8 GHz with EIRP ≥1 × 10¹⁴ W/Hz with a repetition period ≤4.3 hr. We also searched our 11.2 hr of deep observations of the GC and its surrounding region for Fast Radio Burst–like magnetars with the DM up to 5000 pc cm⁻³ with maximum pulse widths up to 90 ms at 6 GHz. We detected several hundred transient bursts from SGR J1745−2900, but did not detect any new transient bursts with the peak luminosity limit across our observed band of ≥10³¹ erg s⁻¹ and burst rate of ≥0.23 burst hr⁻¹. These limits are comparable to bright transient emission seen from other Galactic radio-loud magnetars, constraining their presence at the GC.

We present a new ground-based optical transmission spectrum of the ultrahot Jupiter WASP-103b ( K). Our transmission spectrum is the result of combining five new transits from the ACCESS survey and two new transits from the LRG-BEASTS survey with a reanalysis of three archival Gemini/GMOS transits and one VLT/FORS2 transit. Our combined 11-transit transmission spectrum covers a wavelength range of 3900–9450 Å with a median uncertainty in the transit depth of 148 parts per million, which is less than one atmospheric scale height of the planet. In our retrieval analysis of WASP-103b’s combined optical and infrared transmission spectrum, we find strong evidence for unocculted bright regions (4.3σ) and weak evidence for H₂O (), HCN (), and TiO (), which could be responsible for WASP-103b’s observed temperature inversion. Our optical transmission spectrum shows significant structure that is in excellent agreement with the extensively studied ultrahot Jupiter WASP-121b, for which the presence of VO has been inferred. For WASP-103b, we find that VO can only provide a reasonable fit to the data if its abundance is implausibly high and we do not account for stellar activity. Our results highlight the precision that can be achieved by ground-based observations and the impacts that stellar activity from F-type stars can have on the interpretation of exoplanet transmission spectra.

The Galactic interstellar medium hosts a significant magnetic field, which can be probed through the synchrotron emission produced from its interaction with relativistic electrons. Linearly polarized synchrotron emission is generated throughout the Galaxy and, at longer wavelengths, modified along nearly every path by Faraday rotation in the intervening magneto-ionic medium. Full characterization of the polarized emission requires wideband observations with many frequency channels. We have surveyed polarized radio emission from the Northern sky over the range 1280–1750 MHz, with channel width 236.8 kHz, using the John A. Galt Telescope (diameter 25.6 m) at the Dominion Radio Astrophysical Observatory, as part of the Global Magneto-Ionic Medium Survey. The survey covered 72% of the sky, decl. −30° to +87° at all R.A. The intensity scale was absolutely calibrated, based on the flux density and spectral index of Cygnus A. Polarization angle was calibrated using the extended polarized emission of the Fan Region. Data are presented as brightness temperatures with angular resolution 40′. Sensitivity in Stokes Q and U is 45 mK rms in a 1.18 MHz band. We have applied rotation measure synthesis to the data to obtain a Faraday depth cube of resolution 150 rad m⁻² and sensitivity 3 mK rms of polarized intensity. Features in Faraday depth up to a width of 110 rad m⁻² are represented. The maximum detectable Faraday depth is ±2 × 10⁴ rad m⁻². The survey data are available at the Canadian Astronomy Data Centre.

It is well established that the magnitude of the incident stellar flux is the single most important factor in determining the day–night temperature gradients and atmospheric chemistries of short-period gas giant planets. However, it is likely that other factors, such as planet-to-planet variations in atmospheric metallicity, C/O ratio, and cloud properties, also contribute to the observed diversity of infrared spectra for this population of planets. In this study, we present new 3.6 and 4.5 μm secondary eclipse measurements for five transiting gas giant planets: HAT-P-5b, HAT-P-38b, WASP-7b, WASP-72b, and WASP-127b. We detect eclipses in at least one bandpass for all five planets and confirm circular orbits for all planets except WASP-7b, which shows evidence for a nonzero eccentricity. Building on the work of Garhart et al., we place these new planets into a broader context by comparing them with the sample of all planets with measured Spitzer secondary eclipses. We find that incident flux is the single most important factor for determining the atmospheric chemistry and circulation patterns of short-period gas giant planets. Although we might also expect surface gravity and host star metallicity to play a secondary role, we find no evidence for correlations with either of these two variables.