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

We present a revision to the visual orbit of the young, directly imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric () orbit at an intermediate inclination ( °), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semimajor axis of au and a period of yr, assuming a mass of 1.75 for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigate the potential of using the absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51 Eri derived from a comparison of the Hipparcos and Gaia catalogs was found to be inconsistent at the 2σ–3σ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs and the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using Gaia astrometry alone for a dynamical mass measurement of the planet by simulating Gaia measurements of the motion of the photocenter of the system over the course of the extended 8 yr mission. We find that such a measurement is only possible (>98% probability) given the most optimistic predictions for the Gaia scan astrometric uncertainties for bright stars and a high mass for the planet (≳3.6 M_Jup).

K2-19b and c were among the first planets discovered by NASA’s K2 mission and together stand in stark contrast with the physical and orbital properties of the solar system planets. The planets are between the size of Uranus and Saturn at  7.0 ± 0.2   and  4.1 ± 0.2 , respectively, and reside a mere 0.1% outside the nominal 3:2 mean-motion resonance. They represent a different outcome of the planet formation process than the solar system, as well as the vast majority of known exoplanets. We measured the physical and orbital properties of these planets using photometry from K2, Spitzer, and ground-based telescopes, along with radial velocities from Keck/HIRES. Through a joint photodynamical model, we found that the planets have moderate eccentricities of e ≈ 0.20 and well-aligned apsides Δϖ ≈ 0°. The planets occupy a strictly nonresonant configuration: the resonant angles circulate rather than librate. This defies the predictions of standard formation pathways that invoke convergent or divergent migration, both of which predict Δϖ ≈ 180° and eccentricities of a few percent or less. We measured masses of M_p,b = 32.4 ± 1.7 and M_p,c =  10.8 ± 0.6 . Our measurements, with 5% fractional uncertainties, are among the most precise of any sub-Jovian exoplanet. Mass and size reflect a planet’s core/envelope structure. Despite having a relatively massive core of  ≈ 15 , K2-19b is envelope-rich, with an envelope mass fraction of roughly 50%. This planet poses a challenge to standard models of core-nucleated accretion, which predict that cores ≳10 will quickly accrete gas and trigger runaway accretion when the envelope mass exceeds that of the core.

Future space telescopes may be able to directly image between approximately 10 and 100 planets with sizes and orbits consistent with habitable surface conditions (“exo-Earth candidates” or EECs), but observers will face difficulty in distinguishing these from the potentially hundreds of nonhabitable “false positives” which will also be detected. To maximize the efficiency of follow-up observations, a prioritization scheme must be developed to determine which planets are most likely to be EECs. In this paper, we present a Bayesian method for estimating the likelihood that any directly imaged extrasolar planet is a true EEC by interpreting the planet’s apparent magnitude and separation in the context of existing exoplanet statistics. As a specific application of this general framework, we use published estimates of the discovery yield of future space-based direct-imaging mission concepts to conduct “mock surveys” in which we compute the likelihood that each detected planet is an EEC. We find that it will be difficult to determine which planets are EECs with >50% confidence using single-band photometry immediately upon their detection. The best way to reduce this ambiguity would be to constrain the orbit of a given planet by revisiting the system multiple times or through a radial velocity precursor survey. Astrometric or radial velocity constraints on the mass of the planet would offer a lesser benefit. Finally, we show that a Bayesian approach to prioritizing targets would improve the follow-up efficiency of a direct imaging survey versus a blind approach using the same data. For example, the prioritized approach could reduce the amount of integration time required for the spectral detection (or rejection) of water absorption in most EECs by a factor of two.

EL CVn is the prototype of a class of eclipsing binaries that consist of an A- or F-type main-sequence star and a hot, low-mass, pre-He white dwarf (pre-He WD), the stripped down remains of the former mass donor. Here we present the first direct detection and characterization of the spectrum of the pre-He WD in EL CVn that was made possible through far-UV spectroscopy with the Hubble Space Telescope Cosmic Origins Spectrograph. These spectra straddle the wavelength range where flux dominance shifts from the pre-He WD to the A star. Radial velocities of both components were measured from the far-UV spectra and new optical spectra from the Apache Point Observatory Astrophysical Research Consortium Echelle Spectrograph. We also obtained fast cadence photometry of the eclipses with the pt5m telescope at the Roque de los Muchachos Observatory. A combined analysis of the velocities and light curve yields the component masses and radii. We applied a Doppler tomography algorithm to reconstruct the individual spectra, and we compared these to models to estimate the effective temperatures. The pre-He WD has low mass (0.176 ± 0.004 M_⊙), is small (0.284 ± 0.003 R_⊙), and is relatively hot (11,800 ± 400 K), and these parameters are approximately consistent with predictions for a star stripped through stable mass transfer. The spectral lines of the pre-He WD show that its atmosphere is H-rich, He-depleted, and metal-poor, probably as the result of elemental diffusion that has occurred since mass transfer ceased.

As part of the Panchromatic Exoplanet Treasury program, we have conducted a spectroscopic study of WASP-79b, an inflated hot Jupiter orbiting an F-type star in Eridanus with a period of 3.66 days. Building on the original WASP and TRAPPIST photometry of Smalley et al., we examine Hubble Space Telescope (HST)/Wide Field Camera 3 (WFC3) (1.125–1.650 μm), Magellan/Low Dispersion Survey Spectrograph (LDSS)-3C (0.6–1 μm) data, and Spitzer data (3.6 and 4.5 μm). Using data from all three instruments, we constrain the water abundance to be −2.20 ≤ log(H₂O) ≤ −1.55. We present these results along with the results of an atmospheric retrieval analysis, which favor inclusion of FeH and H⁻ in the atmospheric model. We also provide an updated ephemeris based on the Smalley, HST/WFC3, LDSS-3C, Spitzer, and Transiting Exoplanet Survey Satellite (TESS) transit times. With the detectable water feature and its occupation of the clear/cloudy transition region of the temperature/gravity phase space, WASP-79b is a target of interest for the approved James Webb Space Telescope (JWST) Director’s Discretionary Early Release Science (ERS) program, with ERS observations planned to be the first to execute in Cycle 1. Transiting exoplanets have been approved for 78.1 hr of data collection, and with the delay in the JWST launch, WASP-79b is now a target for the Panchromatic Transmission program. This program will observe WASP-79b for 42 hr in four different instrument modes, providing substantially more data by which to investigate this hot Jupiter.

Self-consistent model grids of brown dwarfs involve complex physics and chemistry, and are often computed using proprietary computer codes, making it challenging to identify the reasons for discrepancies between model and data as well as between the models produced by different research groups. In the current study, we demonstrate a novel method for analyzing brown dwarf spectra, which combines the use of the Sonora, AMES-cond, and HELIOS model grids with the supervised machine-learning method of the random forest. Besides performing atmospheric retrieval, the random forest enables information content analysis of the three model grids as a natural outcome of the method, both individually on each grid and by comparing the grids against one another, by computing large suites of mock retrievals. Our analysis reveals that the different choices made in modeling the alkali line shapes hinder the use of the alkali lines as gravity indicators. Nevertheless, the spectrum longward of 1.2 μm encodes enough information on the surface gravity to allow its inference from retrieval. Temperature may be accurately and precisely inferred independent of the choice of model grid, but not the surface gravity. We apply random forest retrieval to three objects: the benchmark T7.5 brown dwarf GJ 570D, and Indi Ba (T1.5 brown dwarf) and Bb (T6 brown dwarf), which are part of a binary system and have measured dynamical masses. For GJ 570D, the inferred effective temperature and surface gravity are consistent with previous studies. For Indi Ba and Bb, the inferred surface gravities are broadly consistent with the values informed by the dynamical masses.

We present results for the inflated super-Neptune HATS-8b from the Michigan Optical Planetary Spectra Survey (MOPSS). This program is aimed at creating a database of optical planetary transmission spectra all observed, reduced, and analyzed with a uniform method for the benefit of enabling comparative exoplanet studies. HATS-8b orbits a G dwarf and is a low-density super-Neptune with a radius of 0.873 R_Jup, a mass of 0.138 M_Jup, and a density of 0.259 g cm⁻³. Two transits of HATS-8b were observed in 2017 July and August with the Inamori-Magellan Areal Camera and Spectrograph (IMACS) instrument on the Magellan Baade 6.5 m telescope. We find an enhanced scattering slope on each night that agree within 2.3σ. This slope is stronger than one due only to Rayleigh scattering and cannot be fully explained by unocculted starspots. We explore the impact of condensates on the scattering slope and determine that MnS particulates smaller than 10⁻²μm can explain up to 80% of our measured slope if the planet is warmer than equilibrium, or 50% of the slope at the equilibrium temperature of the planet for a low mean molecular weight atmosphere. The scattering slope that we observe is thus beyond even the most extreme haze case we consider. We suggest further follow up on this target and host star to determine if the temporal variation of the slope is primarily due to stellar or planetary effects, and to better understand what these effects may be.

In this paper we report the current status of a new research program. The primary goal of the “Vanishing and Appearing Sources during a Century of Observations” project is to search for vanishing and appearing sources using existing survey data to find examples of exceptional astrophysical transients. The implications of finding such objects extend from traditional astrophysics fields to the more exotic searches for evidence of technologically advanced civilizations. In this first paper we present new, deeper observations of the tentative candidate discovered by Villarroel et al. in 2016. We then perform the first searches for vanishing objects throughout the sky by comparing 600 million objects from the US Naval Observatory Catalogue (USNO) B1.0 down to a limiting magnitude of ∼20–21 with the recent Pan-STARRS Data Release-1 (DR1) with a limiting magnitude of ∼23.4. We find about 150,000 preliminary candidates that do not have any Pan-STARRS counterpart within a 30″ radius. We show that these objects are redder and have larger proper motions than typical USNO objects. We visually examine the images for a subset of about 24,000 candidates, superseding the 2016 study with a sample 10 times larger. We find about 100 point sources visible in only one epoch in the red band of the USNO, which may be of interest in searches for strong M-dwarf flares, high-redshift supernovae, or other categories of unidentified red transients.

We report the discovery of a lithium-rich giant, HD 150382, a post-red-giant-branch clump star. The atmospheric parameters, the chemical abundances for 17 elements, and the isotopic ratio ¹²C/¹³C were determined using the equivalent width and the spectral synthesis methods. The lithium abundance was determined via spectral synthesis of the Li i resonance doublet at λ 6708 yielding log (Li)_NLTE = 2.55. The excited-state line at λ 6104 also was analyzed. In addition, we have used theoretical evolutionary tracks to determine the mass and luminosity of HD 150382 and compare it with other well-known lithium-rich giants in the literature. Our results show that HD 150382 is a slowly rotating K-giant star having near-solar metallicity. We also show that HD 150382 is not s-process enriched, which is in agreement with its evolutionary status. Finally, based on the abundance results, we study the possible scenarios of the lithium enrichment.

A wealth of Earth-sized exoplanets will be discovered in the coming years, providing a large pool of candidates from which the targets for the search for life beyond the solar system will be chosen. The target selection process will require the leveraging of all available information in order to maximize the robustness of the target list and make the most productive use of follow-up resources. Here, we present the results of a suite of n-body simulations that demonstrate the degree to which the orbital architecture of the solar system impacts the variability of Earth’s orbital elements. By varying the orbit of Jupiter and keeping the initial orbits of the other planets constant, we demonstrate how subtle changes in solar system architecture could alter the Earth’s orbital evolution—a key factor in the Milankovitch cycles that alter the amount and distribution of solar insolation, thereby driving periodic climate change on our planet. The amplitudes and frequencies of Earth’s modern orbital cycles fall in the middle of the range seen in our runs for all parameters considered—neither unusually fast nor slow, neither large nor small. This finding runs counter to the “Rare Earth” hypothesis, which suggests that conditions on Earth are so unusual that life elsewhere is essentially impossible. Our results highlight how dynamical simulations of newly discovered exoplanetary systems could be used as an additional means to assess the potential targets of biosignature searches, and thereby help focus the search for life to the most promising targets.

The binary fractions of open and globular clusters yield powerful constraints on their dynamical state and evolutionary history. We apply publicly available Bayesian analysis tools to a UBVRIJHK_S photometric catalog of the open cluster NGC 188 to detect and characterize photometric binaries along the cluster main sequence. This technique has the advantage of self-consistently handling photometric errors, missing data in various bandpasses, and star-by-star prior constraints on cluster membership. Simulations are used to verify uncertainties and quantify selection biases in our analysis, illustrating that among binaries with mass ratios >0.5, we recover the binary fraction to better than 7% in the mean, with no significant dependence on binary fraction and a mild dependence on assumed mass-ratio distribution. Using our photometric catalog, we recover the majority (65% ± 11%) of spectroscopically identified main-sequence binaries, including eight of the nine with spectroscopically measured mass ratios. Accounting for incompleteness and systematics, we derive a mass-ratio distribution that rises toward lower mass ratios (within our q > 0.5 analysis domain). We observe a raw binary fraction for solar-type main-sequence stars with mass ratios q > 0.5 of 42% ± 4%, independent of the assumed mass-ratio distribution to within its uncertainties, consistent with literature values for old open clusters but significantly higher than the field solar-type binary fraction. We confirm that the binaries identified by our method are more concentrated than single stars, in agreement with previous studies, and we demonstrate that the binary nature of those candidates that remain unidentified spectroscopically is strongly supported by photometry from Gaia DR2.

The SR 24 multistar system hosts both circumprimary and circumsecondary disks, which are strongly misaligned with each other. The circumsecondary disk is circumbinary in nature. Interestingly, both disks are interacting, and they possibly rotate in opposite directions. To investigate the nature of this unique twin disk system, we present 01 resolution near-infrared polarized intensity images of the circumstellar structures around SR 24, obtained with HiCIAO mounted on the Subaru 8.2 m telescope. Both the circumprimary disk and the circumsecondary disk are resolved and have elongated features. While the position angle of the major axis and radius of the near-IR (NIR) polarization disk around SR 24S are 55° and 137 au, respectively, those around SR 24N are 110° and 34 au, respectively. With regard to overall morphology, the circumprimary disk around SR 24S shows strong asymmetry, whereas the circumsecondary disk around SR 24N shows relatively strong symmetry. Our NIR observations confirm the previous claim that the circumprimary and circumsecondary disks are misaligned from each other. Both the circumprimary and circumsecondary disks show similar structures in ¹²CO observations in terms of its size and elongation direction. This consistency is because both NIR and ¹²CO are tracing surface layers of the flared disks. As the radius of the polarization disk around SR 24N is roughly consistent with the size of the outer Roche lobe, it is natural to interpret the polarization disk around SR 24N as a circumbinary disk surrounding the SR 24Nb–Nc system.

We present a new ground-based visual transmission spectrum of the hot Jupiter WASP-43b, obtained as part of the ACCESS Survey. The spectrum was derived from four transits observed between 2015 and 2018, with combined wavelength coverage between 5300 and 9000 Å and an average photometric precision of 708 ppm in 230 Å bins. We perform an atmospheric retrieval of our transmission spectrum combined with literature Hubble Space Telescope/WFC3 observations to search for the presence of clouds/hazes as well as Na, K, Hα, and H₂O planetary absorption and stellar spot contamination over a combined spectral range of 5318–16420 Å. We do not detect a statistically significant presence of Na i or K i alkali lines, or Hα in the atmosphere of WASP-43b. We find that the observed transmission spectrum can be best explained by a combination of heterogeneities on the photosphere of the host star and a clear planetary atmosphere with H₂O. This model yields a log evidence of 8.26 ± 0.42 higher than a flat (featureless) spectrum. In particular, the observations marginally favor the presence of large, low-contrast spots over the four ACCESS transit epochs with an average covering fraction and temperature contrast ΔT = 132 K ± 132 K. Within the planet’s atmosphere, we recover a log H₂O volume mixing ratio of , which is consistent with previous H₂O abundance determinations for this planet.

We investigated the jet width profile with distance along the jet in the nearby radio galaxy NGC 1052 at radial distances between ∼300 and 4 × 10⁷ Schwarzschild radii (R_S) from the central engine on both their approaching and receding jet sides. The width of jets was measured in images obtained with the Very Long Baseline Interferometry (VLBI) Space Observatory Programme, the Very Long Baseline Array, and the Very Large Array. The jet width profiles of receding jets are apparently consistent with those of approaching jets throughout the measuring distance ranges, indicating symmetry at least up to the sphere of gravitational influence of the central black hole. The power-law index a of the jet width profile (w_jet ∝ r^a, where w_jet is the jet width, r is the distance from the central engine in the unit of R_S) apparently shows a transition from a ∼ 0 to a ∼ 1, i.e., the cylindrical-to-conical jet structures, at a distance of ∼1 × 10⁴R_S. The cylindrical jet shape at small distances is reminiscent of the innermost jets in 3C 84. Both the central engines of NGC 1052 and 3C 84 are surrounded by dense material, part of which is ionized and causes heavy free–free absorption.

To address the statistics of binary stars in the 8 Myr old Upper Scorpius (USco) star formation region, we conducted a speckle interferometric survey of 614 association members more massive than 0.4  (spectral types earlier than M3V) based on the list of Luhman et al. (2018). We resolved 187 pairs, 55 of which are new discoveries. Also using the published data and the Gaia DR2, a catalog of 250 physical binaries was produced. We carefully estimated detection limits for each target and studied binary statistics in the separation range from 006 to 20″ (9–2800 au), as well as clustering at larger separations. The frequency of companions with mass ratios q > 0.3 in this separation range is 0.33 ± 0.04 and 0.35 ± 0.04 for early M- and solar-type stars, respectively, larger by 1.62 ± 0.22 and 1.39 ± 0.18 times compared to field stars of similar masses. The excess is produced mostly by pairs closer than 100 au. At separations from 100 to 10⁴ au, the separation distribution and companion fraction resemble those of solar-type stars in the field. However, unlike in the field, we see a relative deficit of equal-mass binaries at separations of ∼500 au, compared to smaller and larger separations. The distribution of q depends on the separation, with a preference of larger q and a larger fraction of twins with q > 0.95 at smaller separations. The binary population of USco differs from binaries in the field in several ways and suggests that binary statistics is not universal.

Classification and characterization of variable phenomena and transient phenomena are critical for astrophysics and cosmology. These objects are commonly studied using photometric time series or spectroscopic data. Given that many ongoing and future surveys are conducted in a time domain, and given that adding spectra provides further insights but requires more observational resources, it would be valuable to know which objects we should prioritize to have a spectrum in addition to a time series. We propose a methodology in a probabilistic setting that determines a priori which objects are worth taking a spectrum of to obtain better insights, where we focus on the insight of the type of the object (classification). Objects for which we query their spectrum are reclassified using their full spectral information. We first train two classifiers, one that uses photometric data and another that uses photometric and spectroscopic data together. Then for each photometric object we estimate the probability of each possible spectrum outcome. We combine these models in various probabilistic frameworks (strategies), which are used to guide the selection of follow-up observations. The best strategy depends on the intended use, whether it is obtaining more confidence or accuracy. For a given number of candidate objects (127, equal to 5% of the data set) for taking spectra, we improve the class prediction accuracy by 37% as opposed to 20% of a non-naive (non-random) best-baseline strategy. Our approach provides a general framework for follow-up strategies and can be extended beyond classification to include other forms of follow-ups beyond spectroscopy.

We report long-slit spectroscopy of the early major merger galaxies NGC 3395/NGC 3396. The spectra are consistent with those for star-forming galaxies, but there is some indication of LINER-like active galactic nucleus activity in the center of NGC 3396. The total star formation rate in the regions observed is 2.83 M_⊙ yr⁻¹, consistent with estimates for the entire galaxies. The highest abundances are in the centers of the galaxies, with the abundances decreasing with distance. There is a correlation between high abundance and high ionization parameter, both of which can be attributed to the presence of massive stars. Modeling with SB 99 indicates the star-forming regions are younger than 10 Myr. There are 1000–2000 WNL stars in the system, along with several thousand O stars, consistent with the ages of the star-forming regions. The highest electron densities are found in young regions with high star formation rates. The electron temperatures are higher than results for non-interacting galaxies, which is probably due to shock waves produced by the galaxy–galaxy interaction, the outflow of gas from massive stars, and/or collisions between gas clouds in the galaxies. There is star formation in the bridge of material between the galaxies. These regions are among the youngest in the system and have low abundances, suggesting the gas was pulled from the outer parts of the galaxies. X-ray point sources, probably high-mass X-ray binaries, are associated with several star-forming regions.

The upcoming launch of the James Webb Space Telescope (JWST) will dramatically increase our understanding of exoplanets, particularly through direct imaging. Microlensing and radial velocity surveys indicate that some M dwarfs host long-period giant planets. Some of these planets will likely be just a few parsecs away and a few astronomical units from their host stars, a parameter space that cannot be probed by existing high-contrast imagers. We studied whether the coronagraphs on the Mid-infrared Instrument on JWST can detect Jovian-type planets around nearby M dwarfs. For a sample of 27 very nearby M dwarfs, we simulated a sample of Saturn–Jupiter-mass planets with three atmospheric configurations and three orbital separations, observed in three different filters. We found that the f1550c 15.5 μm filter is best suited for detecting Jupiter-like planets. Jupiter-like planets with patchy cloud cover, 2 au from their star, are detectable at 15.5 μm around 14 stars in our sample, while Jupiters with clearer atmospheres are detectable around all stars in the sample. Saturns were most detectable at 10.65 and 11.4 μm (f1065c and f1140c filters), but only with cloud-free atmospheres and within 3 pc (six stars). Surveying all 27 stars would take <170 hr of JWST integration time, or just a few hours for a shorter survey of the most favorable targets. There is one potentially detectable known planet in our sample: GJ 832 b. Observations aimed at detecting this planet should occur in 2024–2026, when the planet is maximally separated from the star.

The Transiting Exoplanet Survey Satellite (TESS) is finding transiting planet candidates around bright, nearby stars across the entire sky. The large field of view, however, results in low spatial resolution; therefore, multiple stars contribute to almost every TESS light curve. High angular resolution imaging can detect the previously unknown companions to planetary candidate hosts that dilute the transit depths, lead to host star ambiguity, and, in some cases, are the source of false-positive transit signals. We use speckle imaging on the Southern Astrophysical Research (SOAR) telescope to search for companions to 542 TESS planet candidate hosts in the southern sky. We provide correction factors for the 117 systems with resolved companions due to photometric contamination. The contamination in TESS due to close binaries is similar to that found in surveys of Kepler planet candidates. For the solar-type population, we find a deep deficit of close binary systems with projected stellar separations less than 100 au among planet candidate hosts (44 observed binaries compared to 124 expected based on field binary statistics). The close binary suppression among TESS planet candidate hosts is similar to that seen for the more distant Kepler population. We also find a large surplus of TESS planet candidates in wide binary systems detected in both SOAR and Gaia DR2 (119 observed binaries compared to 77 expected). These wide binaries almost exclusively host giant planets, however, suggesting that orbital migration caused by perturbations from the stellar companion may lead to planet–planet scattering and suppress the population of small planets in wide binaries. Both trends are also apparent in the M dwarf planet candidate hosts.

Massive bodies undergo orbital eccentricity oscillations when embedded in an axisymmetric disk of smaller mass orbits. These eccentricity oscillations are driven by secular torques that seek to equalize the apsidal precession rates of all the orbits in the disk. We investigate this mechanism within the context of detached objects in the outer solar system, but we find that the oscillation timescale is too long for it to be dynamically important. It could however be interesting for phenomenon a bit farther from home, namely, feeding supermassive black holes and polluting the surfaces of white dwarf stars.

We present 10 years of new photometric monitoring of the light variability of five evolved stars with strong mid-infrared emission from surrounding dust. Three are known carbon-rich proto-planetary nebulae (PPNe) with F−G spectral types; the nature of the other two was previously unknown. For the three PPNe, we determine or refine the pulsation periods of IRAS 04296+3429 (71 days), 06530−0213 (80 days), and 23304+6147 (84 days). A secondary period was found for each, with a period ratio P₂/P₁ of 0.9. The light variations are small, 0.1–0.2 mag. These are similar to values found in other PPNe. The other two are found to be giant stars. IRAS 09296+1159 pulsates with a period of only 47 days but reaches pulsational light variations of 0.5 mag. Supplemental spectroscopy reveals the spectrum of a CH carbon star. IRAS 08359−1644 is a G1 III star that does not display pulsational variability; rather, it shows nonperiodic decreases of brightness of up to 0.5 mag over this 10 year interval. These drops in brightness are reminiscent of the light curves of R Corona Borealis variables, but with much smaller decreases in brightness and are likely due to transient dust obscuration. Its spectral energy distribution is very similar to that of the unusual oxygen-rich giant star HDE 233517, which possesses mid-infrared hydrocarbon emission features. These two non-PPNe turn out to be members of the rare group of giant stars with large mid-infrared excesses due to dust, objects which presumably have interesting evolutionary histories.

Integral field spectroscopy can map astronomical objects spatially and spectroscopically. Due to instrumental and atmospheric effects, it is common for integral field instruments to yield a sampling of the sky image that is both irregular and wavelength dependent. Most subsequent analysis procedures require a regular, wavelength-independent sampling (for example a fixed rectangular grid), and thus an initial step of fundamental importance is to resample the data onto a new grid. The best possible resampling would produce a well-sampled image, with a resolution equal to that imposed by the intrinsic spatial resolution of the instrument, telescope, and atmosphere, and with no statistical correlations between neighboring pixels. A standard method in the field to produce a regular set of samples from an irregular set of samples is Shepard’s method, but Shepard’s method typically yields images with a degraded resolution and large statistical correlations between pixels. Here we introduce a new method, which improves on Shepard’s method in both these respects. We apply this method to data from the Mapping Nearby Galaxies at Apache Point Observatory survey, part of Sloan Digital Sky Survey IV, demonstrating a full width at half maximum close to that of the intrinsic spatial resolution (and ∼16% better than Shepard’s method) and low statistical correlations between pixels. These results nearly achieve the ideal resampling. This method can have broader applications to other integral field data sets and to other astronomical data sets (such as dithered images) with irregular sampling.

Accurately modeling effects from stellar activity is a key step in detecting radial velocity (RV) signals of low-mass and long-period exoplanets. RVs from stellar activity are dominated by magnetic active regions that move in and out of sight as the star rotates, producing signals with timescales related to the stellar rotation period. Methods to characterize RV periodograms assume that peaks from magnetic active regions will typically occur at the stellar rotation period or a related harmonic. However, with surface features unevenly spaced and evolving over time, signals from magnetic activity are not perfectly periodic, and the effectiveness of characterizing them with sine curves is unconfirmed. With a series of simulations, we perform the first test of common assumptions about signals from magnetic active regions in RV periodograms. We simulate RVs with quasi-periodic signals that account for evolution and migration of magnetic surface features. As test cases, we apply our analysis to two exoplanet hosts, Kepler-20 and K2-131. Simulating observing schedules and uncertainties of real RV surveys, we find that magnetic active regions commonly produce maximum periodogram peaks at spurious periods unrelated to the stellar rotation period: 81% and 72% of peaks, respectively, for K2-131 and Kepler-20. These unexpected peaks can potentially lead to inaccuracies in derived planet masses. We also find that these spurious peaks can sometimes survive multiple seasons of observation, imitating signals typically attributed to exoplanet companions.

We present the first high-resolution spectra for the eclipsing binary V404 Lyr showing γ Dor pulsations, which we use to study its absolute properties. By fitting models to the disentangling spectrum of the primary star, we found that it has an effective temperature of T_eff,1 = 7330 ± 150 K and a rotational velocity of v₁ sin i = 148 ± 18 km s⁻¹. The simultaneous analysis of our double-lined radial velocities and the pulsation-subtracted Kepler data gives us accurate stellar and system parameters of V404 Lyr. The masses, radii, and luminosities are M₁ = 2.17 ± 0.06 M_⊙, R₁ = 1.91 ± 0.02 R_⊙, and L₁ = 9.4 ± 0.8 L_⊙ for the primary, and M₂ = 1.42 ± 0.04 M_⊙, R₂ = 1.79 ± 0.02 R_⊙, and L₂ = 2.9 ± 0.2 L_⊙ for the secondary. The tertiary component orbiting the eclipsing pair has a mass of M_3b = 0.71 ± 0.15 M_⊙ in an orbit of P_3b = 642 ± 3 days, e_3b = 0.21 ± 0.04, and a_3b = 509 ± 2 R_⊙. The third light of l₃ = 4.1% ± 0.2% could be partly attributable to the K-type circumbinary object. By applying a multiple frequency analysis to the eclipse-subtracted light residuals, we detected 45 frequencies with signal-to-noise amplitude ratios larger than 4.0. Identified as independent pulsation modes, seven frequencies (f₁ − f₆, f₉), their new pulsation constants, and the location in the Hertzsprung–Russell diagram indicate that the pulsating primary is a γ Dor-type variable star.

The intranight trajectories of asteroids can be approximated by straight lines, and so are their intranight detections. Therefore, the Hough transform, a line detecting algorithm, can be used to connect the line-up detections to find asteroids. We applied this algorithm to a high-cadence Pan-STARRS 1 (PS1) observation, which was originally designed to collect asteroid light curves for rotation period measurements. The algorithm recovered most of the known asteroids in the observing fields and, moreover, discovered 3574 new asteroids with magnitude mainly of 21.5 < w_p1 < 22.5 mag. This magnitude range is equivalent to subkilometer main-belt asteroids (MBAs), which usually lack rotation period measurements due to their faintness. Using the light curves of the 3574 new asteroids, we obtained 122 reliable rotation periods, of which 13 are super-fast rotators (SRFs; i.e., rotation period of <2 hr). The required cohesion to survive these SFRs range from tens to thousands of Pa, a value consistent with the known SFRs and the regolith on the Moon and Mars. The higher chance of discovering SFRs here suggests that subkilometer MBAs probably harbor more SFRs.

A few hours after its encounter with Pluto, the Alice UV instrument on the New Horizons spacecraft observed the simultaneous stellar occultation and appulse of two UV-bright stars, HD 43153 and HD 42545, respectively. As during the solar occultation that occurred immediately prior, the transmission of starlight through Pluto’s atmosphere was sensitive to absorption by N₂, CH₄, C₂H₆, C₂H₂, C₂H₄, and haze. The line of sight to each star passed over different areas of Pluto from those probed by the solar occultation, providing insight into the degree of spatial and diurnal variability of atmospheric composition on Pluto. In this work, we report the profile of CH₄ derived from both the stellar occultation and appulse observation, as well as the profiles of C₂H₆, C₂H₂, C₂H₄, and haze for the stellar occultation only. We compare these profiles to the solar occultation results and discuss implications for atmospheric variability on Pluto.

Kuiper Belt Object (385446) Manwë–Thorondor is a multiobject system with mutual events predicted to occur from 2014 to 2019. To detect the events, we observed the system at 4 epochs (UT 2016 August 25 and 26, 2017 July 22 and 25, 2017 November 9, and 2018 October 6) in g, r, and VR bands using the 4 m SOAR and the 8.1 m Gemini South telescopes at Cerro Pachón, Chile, and Lowell Observatory’s 4.3 m Discovery Channel Telescope at Happy Jack, Arizona. These dates overlap the uncertainty range (±0.5 day) for four inferior events (Thorondor eclipsing Manwë). We clearly observe variability for the unresolved system with a double-peaked period 11.88190 ± 0.00005 hr and ∼0.5 mag amplitude together with much longer-term variability. Using a multicomponent model, we simultaneously fit our observations and earlier photometry measured separately for Manwë and Thorondor with the Hubble Space Telescope. Our fit suggests Manwë is bilobed, close to the “barbell” shape expected for a strengthless body with density ∼0.8 g cm⁻³ in hydrostatic equilibrium. For Manwë, we thereby derive maximum width to length ratio ∼0.30, surface area equivalent to a sphere of diameter 190 km, geometric albedo 0.06, mass 1.4 × 10¹⁸ kg, and spin axis oriented ∼75° from Earth’s line of sight. Changes in Thorondor’s brightness by ∼0.6 mag with a ∼300 day period may account for the system’s long-term variability. Mutual events with unexpectedly shallow depth and short duration may account for residuals to the fit. The system is complex, providing a challenging puzzle for future modeling efforts.

Binary stars are a very important source of astronomical information. Those short-period binaries observed by means of both spectroscopy and interferometry can provide, through their corresponding orbits, precise values regarding individual masses as well as orbital parallaxes. For this reason, it is fundamental to attempt to optically resolve double-lined spectroscopic binaries. In the present article, we used the elaboration of a specific algorithm in order to determine the necessary telescope aperture to resolve a concrete spectroscopic binary. To that end, using available photometric and spectroscopic information, we established a 3D model for each spectroscopic binary with an orbit from which we deduced the most probable maximum and minimum separations between the components. In this manner, using the obtained calibrations in our study, we also deduced different physical parameters of the components of each system. The methodology presented here was tested satisfactorily with a list of spectroscopic binaries that had been optically resolved and that have both the spectroscopic and visual orbits.

We analyze the full Kepler phase curve of KOI-964, a binary system consisting of a hot white dwarf on an eclipsing orbit around an A-type host star. Using all 18 quarters of long-cadence photometry, we carry out a joint light-curve fit and obtain improved phase-curve amplitudes, occultation depths, orbital parameters, and transit ephemeris over the previous results of Carter et al. A periodogram of the residuals from the phase-curve fit reveals an additional stellar variability signal from the host star with a characteristic period of 0.620276 ± 0.000011 days and a full amplitude of 24 ± 2 ppm. We also present new Keck/HIRES radial velocity observations, which we use to measure the orbit and obtain a mass ratio of q = 0.106 ± 0.012. Combining this measurement with the results of a stellar isochrone analysis, we find that the masses of the host star and white dwarf companion are 2.23 ± 0.12 M_⊙ and , respectively. The effective temperatures of the two components are K and 15,080 ± 400 K, respectively, and we determine the age of the system to be Gyr. We use the measured system properties to compute predicted phase-curve amplitudes and find that while the measured Doppler-boosting and mutual illumination components agree well with theory, the ellipsoidal distortion amplitude is significantly underestimated. We detail possible explanations for this discrepancy, including interactions between the dynamical tide of the host star and the tidal bulge and possible nonsynchronous rotation of the host star.

Due to their ubiquity and very long main-sequence lifetimes, M dwarfs provide an excellent tool to study the formation and chemical enrichment history of our Galaxy. However, owing to their intrinsic faintness, the acquisition of high-resolution, high signal-to-noise spectra of low-mass stars has been limited to small numbers of very nearby stars, mostly from the Galactic disk population. On the other hand, large numbers of low- to medium-resolution spectra of M-type dwarf stars from both the local Galactic disk and halo are available from various surveys. In order to fully exploit these data, we develop a template-fit method using a set of empirically assembled M dwarf/subdwarf classification templates, based on the measurements of the TiO and CaH molecular bands near 7000 Å, which are used to classify M dwarfs/subdwarfs by spectral type and metallicity class. We further present a pipeline to automatically determine the effective temperature T_eff, metallicity [M/H], α-element to iron abundance ratio [α/Fe], and surface gravity log g of M dwarfs/subdwarfs using the latest version of BT-Settl model atmospheres. We apply these methods to a set of low- to medium-resolution spectra of 1544 high proper-motion (μ ≥ 04 yr⁻¹) M dwarfs/subdwarfs, collected at the MDM observatory, Lick Observatory, Kitt-Peak National Observatory, and Cerro-Tololo Interamerican Observatory. Our metallicity estimates appear to be consistent with the expected color–magnitude variation of stars relative to the atmospheric composition, as our sample shows a clear stratification with respect to metallicity in the Hertzsprung–Russel diagram constructed from their Gaia DR2 parallaxes and optical magnitudes. Furthermore, the measured chemical parameters of the two components in 48 binary systems are in good agreement with each other, which suggest a precision of ±0.22 dex in [M/H], ±0.08 dex in [α/Fe], and ±0.16 dex in the combined index [α/Fe] + [M/H]. We find that the relationship between color and spectral subtype depends on metallicity class, as the color G_BP − G_RP is more sensitive to subtype for metal-rich M dwarfs in comparison to metal-poor M subdwarfs. We also demonstrate that effective temperature as a function of spectral subtype has a steeper slope for metal-rich M dwarfs than metal-poor M subdwarfs. There is also a good consistency between “metallicity class,” obtained from the empirical classification templates, and the index [α/Fe] + [M/H] (∼[α/H]), obtained from BT-Settl model fitting, which means that the more easily measured “metallicity class” can be used as a relatively reliable indicator of absolute α-element abundance, [α/H], in low-mass stars. Finally, we examine the distribution of our stars in the [α/Fe] versus [M/H] diagram, which shows evidence of clustering in chemical abundance makeup, suggestive of discrete populations among the local disk and halo stars. We predict that analyses of larger samples of spectra of nearby M-type stars will uncover a complex structure of our Galaxy.

We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light, and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet–disk system. Uniquely, HD 98363 has a wide comoving stellar companion, Wray 15-788, with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one nondetection to one of the brightest disks in the region. With the new results, a total of 15 A/F stars in Sco-Cen have resolved scattered-light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F stars at this important stage in the development of planetary systems.

Transiting planets in nearby young clusters offer the opportunity to study the atmospheres and dynamics of planets during their formative years. To this end, we focused on K2-25b—a close-in (P = 3.48 days), Neptune-sized exoplanet orbiting a M4.5 dwarf in the 650 Myr Hyades cluster. We combined photometric observations of K2-25 covering a total of 44 transits and spanning >2 yr, drawn from a mix of space-based telescopes (Spitzer Space Telescope and K2) and ground-based facilities (Las Cumbres Observatory Global Telescope network and MEarth). The transit photometry spanned 0.6–4.5 μm, which enabled our study of K2-25b’s transmission spectrum. We combined and fit each data set at a common wavelength within a Markov Chain Monte Carlo framework, yielding consistent planet parameters. The resulting transit depths ruled out a solar-composition atmosphere for K2-25b for the range of expected planetary masses and equilibrium temperature at a >4σ confidence level, and are consistent with a flat transmission spectrum. Mass constraints and transit observations at a finer grid of wavelengths (e.g., from the Hubble Space Telescope) are needed to make more definitive statements about the presence of clouds or an atmosphere of high mean molecular weight. Our precise measurements of K2-25b’s transit duration also enabled new constraints on the eccentricity of K2-25's orbit. We find K2-25b’s orbit to be eccentric (e > 0.20) for all reasonable stellar densities and independent of the observation wavelength or instrument. The high eccentricity is suggestive of a complex dynamical history and motivates future searches for additional planets or stellar companions.

We have developed a general statistical procedure for analysis of 2D and 3D finite patterns, which is applied to the data from recently released Gaia-ESA catalog DR2. The 2D analysis clearly confirms our former results on the presence of binaries in the former DR1 catalog. Our main objective is the statistical 3D analysis of DR2. For this, it is essential that the DR2 catalog includes parallaxes and data on the proper motion. The analysis allows us to determine for each pair of stars the probability that it is the binary star. This probability is represented by the function depending on the separation. Furthermore, a combined analysis of the separation with proper motion provides a clear picture of binaries with two components of the motion: parallel and orbital. The result of this analysis is an estimate of the average orbital period and mass of the binary system. The catalog we have created involves 80,560 binary candidates.