We report the discovery of a unique object of uncertain nature—but quite possibly a disintegrating asteroid or minor planet—orbiting one star of the widely separated binary TIC 400799224. We initially identified the system in data from TESS Sector 10 via an abnormally shaped fading event in the light curve (hereafter “dips”). Follow-up speckle imaging determined that TIC 400799224 is actually two stars of similar brightness at 0.″62 separation, forming a likely bound binary with projected separation of ∼300 au. We cannot yet determine which star in the binary is host to the dips in flux. ASAS-SN and Evryscope archival data show that there is a strong periodicity of the dips at ∼19.77 days, leading us to believe that an occulting object is orbiting the host star, though the duration, depth, and shape of the dips vary substantially. Statistical analysis of the ASAS-SN data shows that the dips only occur sporadically at a detectable threshold in approximately one out of every three to five transits, lending credence to the possibility that the occulter is a sporadically emitted dust cloud. The cloud is also fairly optically thick, blocking up to 37% or 75% of the light from the host star, depending on the true host. Further observations may allow for greater detail to be gleaned as to the origin and composition of the occulter, as well as to a determination of which of the two stars comprising TIC 400799224 is the true host star of the dips.

In data from the Kepler mission, the normal F3V star KIC 8462852 (Boyajian’s star) was observed to exhibit infrequent dips in brightness that have not been satisfactorily explained. A previous paper reported the first results of a search for other similar stars in a limited region of the sky around the Kepler field. This paper expands on that search to cover the entire sky between declinations of +22°and +68°. Fifteen new candidates with low rates of dipping, referred to as “slow dippers” in Paper I, have been identified. The dippers occupy a limited region of the HR diagram and an apparent clustering in space is found. This latter feature suggests that these stars are attractive targets for SETI searches.

Recent analyses have shown that distant orbits within the scattered disk population of the Kuiper Belt exhibit an unexpected clustering in their respective arguments of perihelion. While several hypotheses have been put forward to explain this alignment, to date, a theoretical model that can successfully account for the observations remains elusive. In this work we show that the orbits of distant Kuiper Belt objects (KBOs) cluster not only in argument of perihelion, but also in physical space. We demonstrate that the perihelion positions and orbital planes of the objects are tightly confined and that such a clustering has only a probability of 0.007% to be due to chance, thus requiring a dynamical origin. We find that the observed orbital alignment can be maintained by a distant eccentric planet with mass ≳10 m _⊕ whose orbit lies in approximately the same plane as those of the distant KBOs, but whose perihelion is 180° away from the perihelia of the minor bodies. In addition to accounting for the observed orbital alignment, the existence of such a planet naturally explains the presence of high-perihelion Sedna-like objects, as well as the known collection of high semimajor axis objects with inclinations between 60° and 150° whose origin was previously unclear. Continued analysis of both distant and highly inclined outer solar system objects provides the opportunity for testing our hypothesis as well as further constraining the orbital elements and mass of the distant planet.

Megaconstellations of thousands to tens of thousands of artificial satellites (satcons) are rapidly being developed and launched. These satcons will have negative consequences for observational astronomy research, and are poised to drastically interfere with naked-eye stargazing worldwide should mitigation efforts be unsuccessful. Here we provide predictions for the optical brightnesses and on-sky distributions of several satcons, including Starlink, OneWeb, Kuiper, and StarNet/GW, for a total of 65,000 satellites on their filed or predicted orbits. We develop a simple model of satellite reflectivity, which is calibrated using published Starlink observations. We use this model to estimate the visible magnitudes and on-sky distributions for these satellites as seen from different places on Earth, in different seasons, and different times of night. For latitudes near 50° north and south, satcon satellites make up a few percent of all visible point sources all night long near the summer solstice, as well as near sunrise and sunset on the equinoxes. Altering the satellites’ altitudes only changes the specific impacts of the problem. Without drastic reduction of the reflectivities, or significantly fewer total satellites in orbit, satcons will greatly change the night sky worldwide.

The planetary and lunar ephemerides called DE440 and DE441 have been generated by fitting numerically integrated orbits to ground-based and space-based observations. Compared to the previous general-purpose ephemerides DE430, seven years of new data have been added to compute DE440 and DE441, with improved dynamical models and data calibration. The orbit of Jupiter has improved substantially by fitting to the Juno radio range and Very Long Baseline Array (VLBA) data of the Juno spacecraft. The orbit of Saturn has been improved by radio range and VLBA data of the Cassini spacecraft, with improved estimation of the spacecraft orbit. The orbit of Pluto has been improved from use of stellar occultation data reduced against the Gaia star catalog. The ephemerides DE440 and DE441 are fit to the same data set, but DE441 assumes no damping between the lunar liquid core and the solid mantle, which avoids a divergence when integrated backward in time. Therefore, DE441 is less accurate than DE440 for the current century, but covers a much longer duration of years −13,200 to +17,191, compared to DE440 covering years 1550–2650.

The atmospheres and accretion disks of planetary-mass and substellar companions provide an unprecedented look into planet and moon formation processes, most notably the frequency and lifetime of circumplanetary disks. In our ongoing effort to leverage the extraordinary sensitivity of the Spitzer/Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 μm to study wide planetary-mass and substellar companions near the diffraction limit, we present point-spread function fitting photometry of archival Spitzer/IRAC images for nine stars (G0 to M4+M7) in nearby star-forming regions or stellar associations that host companions at separations of ρ = 1.″17–12.″33. We detect all system primaries in all four IRAC channels and recover eight low-mass companions in at least one IRAC channel for our sample, five of which have not been resolved previously in IRAC images. We measure nonphotospheric [3.6]–[8.0] colors for four of the system companions (DH Tau B, 2M0441 B, SR 12 c, and ROXs 42B b), confirming or discovering the presence of circumstellar or circum(sub)stellar disks. We detect fluxes consistent with photospheric emission for four other companions (AB Pic b, CHXR 73 b, 1RXS J1609 b, and HD 203030 b) that are unlikely to host disks. Combined with past detections of accretion or disk indicators, we determine the global disk frequency of young (<15 Myr) wide companions with masses near the deuterium-burning limit to be 56% ± 12%.

We model the settlement of the Galaxy by space-faring civilizations in order to address issues related to the Fermi Paradox. We are motivated to explore the problem in a way that avoids assumptions about the agency (i.e., questions of intent and motivation) of any exo-civilization seeking to settle other planetary systems. We begin by considering the speed of an advancing settlement front to determine if the Galaxy can become inhabited with space-faring civilizations on timescales shorter than its age. Our models for the front speed include the directed settlement of nearby settleable systems through the launching of probes with a finite velocity and range. We also include the effect of stellar motions on the long-term behavior of the settlement front which adds a diffusive component to its advance. As part of our model we also consider that only a fraction, f, of planets will have conditions amenable to settlement by the space-faring civilization. The results of these models demonstrate that the Milky Way can be readily filled-in with settled stellar systems under conservative assumptions about interstellar spacecraft velocities and launch rates. We then move on to consider the question of the Galactic steady state achieved in terms of the fraction X of settled planets. We do this by considering the effect of finite settlement civilization lifetimes on the steady states. We find a range of parameters for which 0 <  X < 1, i.e., the Galaxy supports a population of interstellar space-faring civilizations even though some settleable systems are uninhabited. In addition we find that statistical fluctuations can produce local overabundances of settleable worlds. These generate long-lived clusters of settled systems immersed in large regions that remain unsettled. Both results point to ways in which Earth might remain unvisited in the midst of an inhabited galaxy. Finally we consider how our results can be combined with the finite horizon for evidence of previous settlements in Earth’s geologic record. Using our steady-state model we constrain the probabilities for an Earth visit by a settling civilization before a given time horizon. These results break the link between Hart’s famous “Fact A” (no interstellar visitors on Earth now) and the conclusion that humans must, therefore, be the only technological civilization in the Galaxy. Explicitly, our solutions admit situations where our current circumstances are consistent with an otherwise settled, steady-state galaxy.

We present spectral and photometric observations of 10 Type Ia supernovae (SNe Ia) in the redshift range 0.16 ≤ z ≤ 0.62. The luminosity distances of these objects are determined by methods that employ relations between SN Ia luminosity and light curve shape. Combined with previous data from our High- z Supernova Search Team and recent results by Riess et al., this expanded set of 16 high-redshift supernovae and a set of 34 nearby supernovae are used to place constraints on the following cosmological parameters: the Hubble constant ( H ₀), the mass density (Ω _M ), the cosmological constant (i.e., the vacuum energy density, Ω _Λ), the deceleration parameter ( q ₀), and the dynamical age of the universe ( t ₀). The distances of the high-redshift SNe Ia are, on average, 10%–15% farther than expected in a low mass density (Ω _M = 0.2) universe without a cosmological constant. Different light curve fitting methods, SN Ia subsamples, and prior constraints unanimously favor eternally expanding models with positive cosmological constant (i.e., Ω _Λ > 0) and a current acceleration of the expansion (i.e., q ₀ < 0). With no prior constraint on mass density other than Ω _M ≥ 0, the spectroscopically confirmed SNe Ia are statistically consistent with q ₀ < 0 at the 2.8 σ and 3.9 σ confidence levels, and with Ω _Λ > 0 at the 3.0 σ and 4.0 σ confidence levels, for two different fitting methods, respectively. Fixing a "minimal" mass density, Ω _M = 0.2, results in the weakest detection, Ω _Λ > 0 at the 3.0 σ confidence level from one of the two methods. For a flat universe prior (Ω _M + Ω _Λ = 1), the spectroscopically confirmed SNe Ia require Ω _Λ > 0 at 7 σ and 9 σ formal statistical significance for the two different fitting methods. A universe closed by ordinary matter (i.e., Ω _M = 1) is formally ruled out at the 7 σ to 8 σ confidence level for the two different fitting methods. We estimate the dynamical age of the universe to be 14.2 ± 1.7 Gyr including systematic uncertainties in the current Cepheid distance scale. We estimate the likely effect of several sources of systematic error, including progenitor and metallicity evolution, extinction, sample selection bias, local perturbations in the expansion rate, gravitational lensing, and sample contamination. Presently, none of these effects appear to reconcile the data with Ω _Λ = 0 and q ₀ ≥ 0.

We present exoplanet occurrence rates estimated with approximate Bayesian computation for planets with radii between 0.5 and 16 R _⊕ and orbital periods between 0.78 and 400 days orbiting FGK dwarf stars. We base our results on an independent planet catalog compiled from our search of all ∼200,000 stars observed over the Kepler mission, with precise planetary radii supplemented by Gaia DR2-incorporated stellar radii. We take into account detection and vetting efficiency, planet radius uncertainty, and reliability against transit-like noise signals in the data. By analyzing our FGK occurrence rates as well as those computed after separating F-, G-, and K-type stars, we explore dependencies on stellar effective temperature, planet radius, and orbital period. We reveal new characteristics of the photoevaporation-driven “radius gap” between ∼1.5 and 2 R _⊕, indicating that the bimodal distribution previously revealed for P < 100 days exists only over a much narrower range of orbital periods, above which sub-Neptunes dominate and below which super-Earths dominate. Finally, we provide several estimates of the “eta-Earth” value—the frequency of potentially habitable, rocky planets orbiting Sun-like stars. For planets with sizes 0.75–1.5 R _⊕ orbiting in a conservatively defined habitable zone (0.99–1.70 au) around G-type stars, we place an upper limit (84.1th percentile) of <0.18 planets per star.

The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package, as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of interoperable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy Project.

We perform a study of stellar flares for the 24,809 stars observed with 2 minute cadence during the first two months of the TESS mission. Flares may erode exoplanets’ atmospheres and impact their habitability, but might also trigger the genesis of life around small stars. TESS provides a new sample of bright dwarf stars in our galactic neighborhood, collecting data for thousands of M dwarfs that might host habitable exoplanets. Here, we use an automated search for flares accompanied by visual inspection. Then, our public allesfitter code robustly selects the appropriate model for potentially complex flares via Bayesian evidence. We identify 1228 flaring stars, 673 of which are M dwarfs. Among 8695 flares in total, the largest superflare increased the stellar brightness by a factor of 16.1. Bolometric flare energies range from 10 ^31.0 to 10 ^36.9 erg, with a median of 10 ^33.1 erg. Furthermore, we study the flare rate and energy as a function of stellar type and rotation period. We solidify past findings that fast rotating M dwarfs are the most likely to flare and that their flare amplitude is independent of the rotation period. Finally, we link our results to criteria for prebiotic chemistry, atmospheric loss through coronal mass ejections, and ozone sterilization. Four of our flaring M dwarfs host exoplanet candidates alerted on by TESS, for which we discuss how these effects can impact life. With upcoming TESS data releases, our flare analysis can be expanded to almost all bright small stars, aiding in defining criteria for exoplanet habitability.

The spectral analysis and data products in Data Release 16 (DR16; 2019 December) from the high-resolution near-infrared Apache Point Observatory Galactic Evolution Experiment (APOGEE)-2/Sloan Digital Sky Survey (SDSS)-IV survey are described. Compared to the previous APOGEE data release (DR14; 2017 July), APOGEE DR16 includes about 200,000 new stellar spectra, of which 100,000 are from a new southern APOGEE instrument mounted on the 2.5 m du Pont telescope at Las Campanas Observatory in Chile. DR16 includes all data taken up to 2018 August, including data released in previous data releases. All of the data have been re-reduced and re-analyzed using the latest pipelines, resulting in a total of 473,307 spectra of 437,445 stars. Changes to the analysis methods for this release include, but are not limited to, the use of MARCS model atmospheres for calculation of the entire main grid of synthetic spectra used in the analysis, a new method for filling “holes” in the grids due to unconverged model atmospheres, and a new scheme for continuum normalization. Abundances of the neutron-capture element Ce are included for the first time. A new scheme for estimating uncertainties of the derived quantities using stars with multiple observations has been applied, and calibrated values of surface gravities for dwarf stars are now supplied. Compared to DR14, the radial velocities derived for this release more closely match those in the Gaia DR2 database, and a clear improvement in the spectral analysis of the coolest giants can be seen. The reduced spectra as well as the result of the analysis can be downloaded using links provided on the SDSS DR16 web page.

We present calculations of the occurrence rate of small close-in planets around low-mass dwarf stars using the known planet populations from the Kepler and K2 missions. Applying completeness corrections clearly reveals the radius valley in the maximum a posteriori occurrence rates as a function of orbital separation and planet radius. We measure the slope of the valley to be , which bears the opposite sign from that measured around Sun-like stars, thus suggesting that thermally driven atmospheric mass loss may not dominate the evolution of planets in the low stellar mass regime or that we are witnessing the emergence of a separate channel of planet formation. The latter notion is supported by the relative occurrence of rocky to non-rocky planets increasing from 0.5 ± 0.1 around mid-K dwarfs to 8.5 ± 4.6 around mid-M dwarfs. Furthermore, the center of the radius valley at 1.54 ± 0.16 R _⊕ is shown to shift to smaller sizes with decreasing stellar mass, in agreement with physical models of photoevaporation, core-powered mass loss, and gas-poor formation. Although current measurements are insufficient to robustly identify the dominant formation pathway of the radius valley, such inferences may be obtained by the Transiting Exoplanet Survey Satellite with (85,000) mid-to-late M dwarfs observed with 2 minutes cadence. The measurements presented herein also precisely designate the subset of planetary orbital periods and radii that should be targeted in radial velocity surveys to resolve the rocky to non-rocky transition around low-mass stars.

An accurate and precise Kepler Stellar Properties Catalog is essential for the interpretation of the Kepler exoplanet survey results. Previous Kepler Stellar Properties Catalogs have focused on reporting the best-available parameters for each star, but this has required combining data from a variety of heterogeneous sources. We present the Gaia–Kepler Stellar Properties Catalog, a set of stellar properties of 186,301 Kepler stars, homogeneously derived from isochrones and broadband photometry, Gaia Data Release 2 parallaxes, and spectroscopic metallicities, where available. Our photometric effective temperatures, derived from colors, are calibrated on stars with interferometric angular diameters. Median catalog uncertainties are 112 K for , 0.05 dex for , 4% for , 7% for , 13% for , 10% for , and 56% for stellar age. These precise constraints on stellar properties for this sample of stars will allow unprecedented investigations into trends in stellar and exoplanet properties as a function of stellar mass and age. In addition, our homogeneous parameter determinations will permit more accurate calculations of planet occurrence and trends with stellar properties.

We derive analytic, closed-form solutions for the light curve of a planet transiting a star with a limb-darkening profile that is a polynomial function of the stellar elevation, up to an arbitrary integer order. We provide improved analytic expressions for the uniform, linear, and quadratic limb-darkened cases, as well as novel expressions for higher-order integer powers of limb darkening. The formulae are crafted to be numerically stable over the expected range of usage. We additionally present analytic formulae for the partial derivatives of instantaneous flux with respect to the radius ratio, impact parameter, and limb-darkening coefficients. These expressions are rapid to evaluate and compare quite favorably in speed and accuracy to existing transit light-curve codes. We also use these expressions to numerically compute the first partial derivatives of exposure-time-averaged transit light curves with respect to all model parameters. An additional application is modeling eclipsing binary or eclipsing multiple star systems in cases where the stars may be treated as spherically symmetric. We provide code which implements these formulae in C++, Python, IDL, and Julia, with tests and examples of usage ( https://github.com/rodluger/Limbdark.jl).

We present a catalog of unbound stellar pairs, within 100 pc of the Sun, that are undergoing close, hyperbolic, encounters. The data are drawn from the GAIA EDR3 catalog, and the limiting factors are errors in the radial distance and unknown velocities along the line of sight. Such stellar pairs have been suggested to be possible events associated with the migration of technological civilizations between stars. As such, this sample may represent a finite set of targets for a SETI search based on this hypothesis. Our catalog contains a total of 132 close passage events, featuring stars from across the entire main sequence, with 16 pairs featuring at least one main-sequence star of spectral type between K1 and F3. Many of these stars are also in binaries, so that we isolate eight single stars as the most likely candidates to search for an ongoing migration event—HD 87978, HD 92577, HD 50669, HD 44006, HD 80790, LSPM J2126+5338, LSPM J0646+1829 and HD 192486. Among host stars of known planets, the stars GJ 433 and HR 858 are the best candidates.

A challenge in absolute calibration is to relate very bright stars with physical flux measurements to faint ones within range of modern instruments, e.g., those on large ground-based telescopes or the James Webb Space Telescope (JWST). We propose Sirius as the fiducial color standard. It is an A0V star that is slowly rotating and does not have infrared excesses due to either hot dust or a planetary debris disk; it also has a number of accurate (∼1%–2%) absolute flux measurements. We accurately transfer the near-infrared flux from Sirius to BD +60 1753, an unobscured early A-type star (A1V, V ≈ 9.6, E(B – V) ≈ 0.009) that is faint enough to serve as a primary absolute flux calibrator for JWST. Its near-infrared spectral energy distribution and that of Sirius should be virtually identical. We have determined its output relative to that of Sirius in a number of different ways, all of which give consistent results within ∼1%. We also transfer the calibration to GSPC P330-E, a well-calibrated close solar analog (G2V). We have emphasized the 2MASS K_S band, since it represents a large number and long history of measurements, but the theoretical spectra (i.e., from CALSPEC) of these stars can be used to extend this result throughout the near- and mid-infrared.

We report Spitzer 3.6 and 4.5 μm photometry of 11 bright stars relative to Sirius, exploiting the unique optical stability of the Spitzer Space Telescope point-spread function (PSF). Spitzer's extremely stable beryllium optics in its isothermal environment enables precise comparisons in the wings of the PSF from heavily saturated stars. These bright stars stand as the primary sample to improve stellar models, and to transfer the absolute flux calibration of bright standard stars to a sample of fainter standards useful for missions like JWST and for large ground-based telescopes. We demonstrate that better than 1% relative photometry can be achieved using the PSF wing technique in the radial range of 20''–100'' for stars that are fainter than Sirius by 8 mag (from outside the saturated core to a large radius where a high signal-to-noise ratio profile can still be obtained). We test our results by (1) comparing the [3.6]−[4.5] color with that expected between the WISE W1 and W2 bands, (2) comparing with stars where there is accurate K_S photometry, and (3) also comparing with relative fluxes obtained with the DIRBE instrument on COBE. These tests confirm that relative photometry is achieved to better than 1%.

We apply the automated AnomalyFinder algorithm of Paper I to 2018–2019 light curves from the ≃13 deg² covered by the six KMTNet prime fields, with cadences Γ ≥ 2 hr⁻¹. We find a total of 11 planets with mass ratios q < 2 × 10⁻⁴, including 6 newly discovered planets, 1 planet that was reported in Paper I, and recovery of 4 previously discovered planets. One of the new planets, OGLE-2018-BLG-0977Lb, is in a planetary caustic event, while the other five (OGLE-2018-BLG-0506Lb, OGLE-2018-BLG-0516Lb, OGLE-2019-BLG-1492Lb, KMT-2019-BLG-0253, and KMT-2019-BLG-0953) are revealed by a "dip" in the light curve as the source crosses the host-planet axis on the opposite side of the planet. These subtle signals were missed in previous by-eye searches. The planet-host separations (scaled to the Einstein radius), s, and planet-host mass ratios, q, are, respectively, (s, q × 10⁵) = (0.88, 4.1), (0.96 ± 0.10, 8.3), (0.94 ± 0.07, 13), (0.97 ± 0.07, 18), (0.97 ± 0.04, 4.1), and (0.74, 18), where the " ± " indicates a discrete degeneracy. The 11 planets are spread out over the range . Together with the two planets previously reported with q ∼ 10⁻⁵ from the 2018–2019 nonprime KMT fields, this result suggests that planets toward the bottom of this mass-ratio range may be more common than previously believed.

We present a re-analysis of transit depths of KELT-19Ab, WASP-156b, and WASP-121b, including data from the Transiting Exoplanet Survey Satellite (TESS). The large ∼21'' TESS pixels and point-spread function result in significant contamination of the stellar flux by nearby objects. We use Gaia data to fit for and remove this contribution, providing general-purpose software for this correction. We find all three sources have a larger inclination, compared to earlier work. For WASP-121b, we find significantly smaller values (135) of the inclination when using the 30 minute cadence data compared to the 2 minute cadence data. Using simulations, we demonstrate that the radius ratio of exoplanet to star (R_p/R_*) is biased small relative to data taken with a larger sampling interval although oversampling corrections mitigate the bias. This is particularly important for deriving subpercent transit differences between bands. We find the radius ratio of exoplanet to star (R_p/R_*) in the TESS band is 7.5σ smaller than previous work for KELT-19Ab, but consistent to within ∼2σ for WASP-156b and WASP-121b. The difference could be due to specific choices in the analysis, not necessarily due to the presence of atmospheric features. The result for KELT-19Ab possibly favors a haze-dominated atmosphere. We do not find evidence for the ∼0.95 μm water feature contaminating transit depths in the TESS band for these stars but show that with photometric precision of 500 ppm and with a sampling of about 200 observations across the entire transit, this feature could be detectable in a more narrow z-band.

The Gas Pixel Detector (GPD) is an X-ray polarimeter to fly onboard IXPE and other missions. To correctly measure the source polarization, the response of IXPE’s GPDs to unpolarized radiation has to be calibrated and corrected. In this paper, we describe the way such response is measured with laboratory sources and the algorithm to apply such correction to the observations of celestial sources. The latter allows to correct the response to polarization of single photons, therefore allowing great flexibility in all the subsequent analysis. Our correction approach is tested against both monochromatic and nonmonochromatic laboratory sources and with simulations, finding that it correctly retrieves the polarization up to the statistical limits of the planned IXPE observations.

We present mass and radius measurements of K2-79b and K2-222b, two transiting exoplanets orbiting active G-type stars observed with HARPS-N and K2. Their respective 10.99 day and 15.39 day orbital periods fall near periods of signals induced by stellar magnetic activity. The two signals might therefore interfere and lead to an inaccurate estimate of exoplanet mass. We present a method to mitigate these effects when radial velocity (RV) and activity-indicator observations are available over multiple observing seasons and the orbital period of the exoplanet is known. We perform correlation and periodogram analyses on subsets composed of each target's two observing seasons, in addition to the full data sets. For both targets, these analyses reveal an optimal season with little to no interference at the orbital period of the known exoplanet. We make a confident mass detection of each exoplanet by confirming agreement between fits to the full RV set and the optimal season. For K2-79b, we measure a mass of 11.8 ± 3.6 M _⊕ and a radius of 4.09 ± 0.17 R _⊕. For K2-222b, we measure a mass of 8.0 ± 1.8 M _⊕ and a radius of 2.35 ± 0.08 R _⊕. According to model predictions, K2-79b is a highly irradiated Uranus analog and K2-222b hosts significant amounts of water ice. We also present a RV solution for a candidate second companion orbiting K2-222 at 147.5 days.

We present a re-analysis of transit depths of KELT-19Ab, WASP-156b, and WASP-121b, including data from the Transiting Exoplanet Survey Satellite (TESS). The large ∼21″ TESS pixels and point-spread function result in significant contamination of the stellar flux by nearby objects. We use Gaia data to fit for and remove this contribution, providing general-purpose software for this correction. We find all three sources have a larger inclination, compared to earlier work. For WASP-121b, we find significantly smaller values (13.°5) of the inclination when using the 30 minute cadence data compared to the 2 minute cadence data. Using simulations, we demonstrate that the radius ratio of exoplanet to star ( R _p / R _*) is biased small relative to data taken with a larger sampling interval although oversampling corrections mitigate the bias. This is particularly important for deriving subpercent transit differences between bands. We find the radius ratio of exoplanet to star ( R _p / R _*) in the TESS band is 7.5 σ smaller than previous work for KELT-19Ab, but consistent to within ∼2 σ for WASP-156b and WASP-121b. The difference could be due to specific choices in the analysis, not necessarily due to the presence of atmospheric features. The result for KELT-19Ab possibly favors a haze-dominated atmosphere. We do not find evidence for the ∼0.95 μm water feature contaminating transit depths in the TESS band for these stars but show that with photometric precision of 500 ppm and with a sampling of about 200 observations across the entire transit, this feature could be detectable in a more narrow z-band.

A challenge in absolute calibration is to relate very bright stars with physical flux measurements to faint ones within range of modern instruments, e.g., those on large ground-based telescopes or the James Webb Space Telescope (JWST). We propose Sirius as the fiducial color standard. It is an A0V star that is slowly rotating and does not have infrared excesses due to either hot dust or a planetary debris disk; it also has a number of accurate (∼1%–2%) absolute flux measurements. We accurately transfer the near-infrared flux from Sirius to BD +60 1753, an unobscured early A-type star (A1V, V ≈ 9.6, E( B – V) ≈ 0.009) that is faint enough to serve as a primary absolute flux calibrator for JWST. Its near-infrared spectral energy distribution and that of Sirius should be virtually identical. We have determined its output relative to that of Sirius in a number of different ways, all of which give consistent results within ∼1%. We also transfer the calibration to GSPC P330-E, a well-calibrated close solar analog (G2V). We have emphasized the 2MASS K _S band, since it represents a large number and long history of measurements, but the theoretical spectra (i.e., from CALSPEC) of these stars can be used to extend this result throughout the near- and mid-infrared.

We report Spitzer 3.6 and 4.5 μm photometry of 11 bright stars relative to Sirius, exploiting the unique optical stability of the Spitzer Space Telescope point-spread function (PSF). Spitzer's extremely stable beryllium optics in its isothermal environment enables precise comparisons in the wings of the PSF from heavily saturated stars. These bright stars stand as the primary sample to improve stellar models, and to transfer the absolute flux calibration of bright standard stars to a sample of fainter standards useful for missions like JWST and for large ground-based telescopes. We demonstrate that better than 1% relative photometry can be achieved using the PSF wing technique in the radial range of 20″–100″ for stars that are fainter than Sirius by 8 mag (from outside the saturated core to a large radius where a high signal-to-noise ratio profile can still be obtained). We test our results by (1) comparing the [3.6]−[4.5] color with that expected between the WISE W1 and W2 bands, (2) comparing with stars where there is accurate K _S photometry, and (3) also comparing with relative fluxes obtained with the DIRBE instrument on COBE. These tests confirm that relative photometry is achieved to better than 1%.

In data from the Kepler mission, the normal F3V star KIC 8462852 (Boyajian’s star) was observed to exhibit infrequent dips in brightness that have not been satisfactorily explained. A previous paper reported the first results of a search for other similar stars in a limited region of the sky around the Kepler field. This paper expands on that search to cover the entire sky between declinations of +22°and +68°. Fifteen new candidates with low rates of dipping, referred to as “slow dippers” in Paper I, have been identified. The dippers occupy a limited region of the HR diagram and an apparent clustering in space is found. This latter feature suggests that these stars are attractive targets for SETI searches.

Owing to recent advances in radial-velocity instrumentation and observation techniques, the detection of Earth-mass planets around Sun-like stars may soon be primarily limited by intrinsic stellar variability. Several processes contribute to this variability, including starspots, pulsations, and granulation. Although many previous studies have focused on techniques to mitigate signals from pulsations and other types of magnetic activity, granulation noise has to date only been partially addressed by empirically motivated observation strategies and magnetohydrodynamic simulations. To address this deficit, we present the GRanulation And Spectrum Simulator ( GRASS), a new tool designed to create time-series synthetic spectra with granulation-driven variability from spatially and temporally resolved observations of solar absorption lines. In this work, we present GRASS, detail its methodology, and validate its model against disk-integrated solar observations. As a first-of-its-kind empirical model for spectral variability due to granulation in a star with perfectly known center-of-mass radial-velocity behavior, GRASS is an important tool for testing new methods of disentangling granular line-shape changes from true Doppler shifts.

We study the excitation of mutual inclination between planetary orbits by a novel secular-orbital resonance in multi-planet systems perturbed by binary companions, which we call “ivection.” The ivection resonance happens when the nodal precession rate of the planet matches a multiple of the orbital frequency of the binary, and its physical nature is similar to the previously studied evection resonance. Capture into an ivection resonance requires encountering the resonance with slowly increasing nodal precession rate, and it can excite the mutual inclination of the planets without affecting their eccentricities. We discuss the possible outcomes of ivection resonance capture, and we use simulations to illustrate that it is a promising mechanism for producing the mutual inclination in systems where planets have significant mutual inclination but modest eccentricity, such as Kepler-108. We also find an apparent deficit of multi-planet systems that would have a nodal precession period comparable to the binary orbital period, suggesting that ivection resonance may inhibit formation of or destablize multi-planet systems with an external binary companion.

Second-order mean-motion resonances lead to an interesting phenomenon in the sculpting of the period-ratio distribution, due to their shape and width in period-ratio/eccentricity space. As the osculating periods librate in resonance, the time-averaged period ratio approaches the exact commensurability. The width of second-order resonances increases with increasing eccentricity, and thus more eccentric systems have a stronger peak at commensurability when averaged over sufficient time. The libration period is short enough that this time-averaging behavior is expected to appear on the timescale of the Kepler mission. Using N-body integrations of simulated planet pairs near the 5:3 and 3:1 mean-motion resonances, we investigate the eccentricity distribution consistent with the planet pairs observed by Kepler. This analysis, an approach independent from previous studies, shows no statistically significant peak at the 3:1 resonance and a small peak at the 5:3 resonance, placing an upper limit on the Rayleigh scale parameter, σ, of the eccentricity of the observed Kepler planets at σ = 0.245 (3:1) and σ = 0.095 (5:3) at 95% confidence, consistent with previous results from other methods.

Estimating the amount of foreground extinction due to the Milky Way dust along the line of sight is often a first step in determining the luminosity of an object. The amount of Galactic dust inferred by infrared emission maps can be contaminated by infrared light from nearby galaxies. By comparing extinction values at and around the location of nearby galaxies, we compile a list of 95 galaxies that likely contaminate the maps with an excess or improperly subtracted galaxian infrared emission, and tabulate our recommended values for the MW contribution. In addition to M82, which inspired this work, six more sources have an excess visual extinction A _V of at least 0.05 mag greater than our annular values; including M83, NGC 1313, NGC 6822, NGC 918, UGC 11501, and UGC 11797. M33 is shown to be oversubtracted. NGC 88 and the outskirts of NGC 4258 are located in gaps in the Infrared Astronomical Satellite imaging. The recommended dust map values for the LMC, SMC, and M31 may also not be correctly returned by some software packages. Accurate reddening estimates are important for measuring stellar and supernova luminosities in these nearby galaxies.