Table of contents

Main-sequence, fully convective M dwarfs in eclipsing binaries are observed to be larger than stellar evolutionary models predict by as much as 10%–15%. A proposed explanation for this discrepancy involves effects from strong magnetic fields, induced by rapid rotation via the dynamo process. Although, a handful of single, slowly rotating M dwarfs with radius measurements from interferometry also appear to be larger than models predict, suggesting that rotation or binarity specifically may not be the sole cause of the discrepancy. We test whether single, rapidly rotating, fully convective stars are also larger than expected by measuring their $R\sin i$ distribution. We combine photometric rotation periods from the literature with rotational broadening ($v\sin i$) measurements reported in this work for a sample of 88 rapidly rotating M dwarf stars. Using a Bayesian framework, we find that stellar evolutionary models underestimate the radii by $10 \% \mbox{--}15{ \% }_{-2.5}^{+3}$, but that at higher masses (0.18 < M < 0.4 M_Sun), the discrepancy is only about 6% and comparable to results from interferometry and eclipsing binaries. At the lowest masses (0.08 < M < 0.18 M_Sun), we find that the discrepancy between observations and theory is 13%–18%, and we argue that the discrepancy is unlikely to be due to effects from age. Furthermore, we find no statistically significant radius discrepancy between our sample and the handful of M dwarfs with interferometric radii. We conclude that neither rotation nor binarity are responsible for the inflated radii of fully convective M dwarfs, and that all fully convective M dwarfs are larger than models predict.

We explore KLIP forward modeling spectral extraction on Gemini Planet Imager coronagraphic data of HR 8799, using PyKLIP, and show algorithm stability with varying KLIP parameters. We report new and re-reduced spectrophotometry of HR 8799 c, d, and e in the H and K bands. We discuss a strategy for choosing optimal KLIP PSF subtraction parameters by injecting simulated sources and recovering them over a range of parameters. The K1/K2 spectra for HR 8799 c and d are similar to previously published results from the same data set. We also present a K-band spectrum of HR 8799 e for the first time and show that our H-band spectra agree well with previously published spectra from the VLT/SPHERE instrument. We show that HR 8799 c and d show significant differences in their H and K spectra, but do not find any conclusive differences between d and e, nor between c and e, likely due to large error bars in the recovered spectrum of e. Compared to M-, L-, and T-type field brown dwarfs, all three planets are most consistent with mid- and late-L spectral types. All objects are consistent with low gravity, but a lack of standard spectra for low gravity limit the ability to fit the best spectral type. We discuss how dedicated modeling efforts can better fit HR 8799 planets' near-IR flux, as well as how differences between the properties of these planets can be further explored.

We present Submillimeter Array and Northern Extended Millimeter Array observations of the host galaxy of FRB 121102 in the CO 3–2 and 1–0 transitions, respectively. We do not detect emission from either transition. We set 3σ upper limits to the CO luminosity L_CO < 2.5 × 10⁷ K km s⁻¹ pc⁻² for CO 3–2 and L_CO < 2.3 × 10⁹ K km s⁻¹ pc⁻² for CO 1–0. For Milky Way–like star formation properties, we set a 3σ upper limit on the H₂ mass of 2.5 × 10⁸M_⊙, slightly less than the predictions for the H₂ mass based on the star formation rate. The true constraint on the H₂ mass may be significantly higher, however, because of the reduction in CO luminosity that is common for low-metallicity dwarf galaxies like the FRB host galaxy. These results demonstrate the challenge of identifying the nature of FRB progenitors through study of the host galaxy molecular gas. We also place a limit of 42 μJy (3σ) on the continuum flux density of the persistent radio source at 97 GHz, consistent with a power-law extrapolation of the low-frequency spectrum, which may arise from an active galactic nucleus or other nonthermal source.

Proxima Centauri has become the subject of intense study since the radial-velocity (RV) discovery by Anglada-Escudé et al. of a planet orbiting this nearby M dwarf every ∼11.2 days. If Proxima Centauri b transits its host star, independent confirmation of its existence is possible, and its mass and radius can be measured in units of the stellar host mass and radius. To date, there have been three independent claims of possible transit-like event detections in light curve observations obtained by the MOST satellite (in 2014–15), the Bright Star Survey Telescope telescope in Antarctica (in 2016), and the Las Campanas Observatory (in 2016). The claimed possible detections are tentative, due in part to the variability intrinsic to the host star, and in the case of the ground-based observations, also due to the limited duration of the light curve observations. Here, we present preliminary results from an extensive photometric monitoring campaign of Proxima Centauri, using telescopes around the globe and spanning from 2006 to 2017, comprising a total of 329 observations. Considering our data that coincide directly and/or phased with the previously published tentative transit detections, we are unable to independently verify those claims. We do, however, verify the previously reported ubiquitous and complex variability of the host star. We discuss possible interpretations of the data in light of the previous claims, and we discuss future analyses of these data that could more definitively verify or refute the presence of transits associated with the RV-discovered planet.

The VLBI USNO 2016A (U16A) solution is part of a work-in-progress effort by USNO toward the preparation of the ICRF3. Most of the astrometric improvement with respect to the ICRF2 is due to the re-observation of the VCS sources. Our objective in this paper is to assess U16A's astrometry. A comparison with ICRF2 shows statistically significant offsets of size 0.1 mas between the two solutions. While Gaia DR1 positions are not precise enough to resolve these offsets, they are found to be significantly closer to U16A than ICRF2. In particular, the trend for typically larger errors for southern sources in VLBI solutions is decreased in U16A. Overall, the VLBI-Gaia offsets are reduced by 21%. The U16A list includes 718 sources not previously included in ICRF2. Twenty of those new sources have statistically significant radio-optical offsets. In two-thirds of the cases, these offsets can be explained from PanSTARRS images.

Direct imaging is likely the best way to characterize the atmospheres of Earth-sized exoplanets in the habitable zone of Sun-like stars. Previously, Stark et al. estimated the Earth twin yield of future direct imaging missions, such as LUVOIR and HabEx. We extend this analysis to other types of planets, which will act as false positives for Earth twins. We define an Earth twin as any exoplanet within half an e-folding of 1 au in semimajor axis and 1 ${R}_{\oplus }$ in planetary radius, orbiting a G-dwarf. Using Monte Carlo analyses, we quantify the biases and planetary false-positive rates of Earth searches. That is, given a pale dot at the correct projected separation and brightness to be a candidate Earth, what are the odds that it is, in fact, an Earth twin? Our notional telescope has a diameter of 10 m, an inner working angle of 3λ/D, and an outer working angle of 10λ/D (62 mas and 206 mas at 1.0 μm). With no precursor knowledge and one visit per star, 77% of detected candidate Earths are actually un-Earths; their mean radius is 2.3 ${R}_{\oplus }$, a sub-Neptune. The odds improve if we image every planet at its optimal orbital phase, either by relying on precursor knowledge, or by performing multi-epoch direct imaging. In such a targeted search, 47% of detected Earth twin candidates are false positives, and they have a mean radius of 1.7 ${R}_{\oplus }$. The false-positive rate is insensitive to stellar spectral type and the assumption of circular orbits.

The unique inner-belt asteroid 311P/PANSTARRS (formerly P/2013 P5) is notable for its sporadic, comet-like ejection of dust in nine distinct epochs spread over ∼250 days in 2013. This curious behavior has been interpreted as the product of localized, equatorward landsliding from the surface of an asteroid rotating at the brink of instability. We obtained new Hubble Space Telescope observations to directly measure the nucleus and to search for evidence of its rapid rotation. We find a nucleus with mid-light absolute magnitude H_V = 19.14 ± 0.02, corresponding to an equal-area circle with radius 190 ± 30 m (assuming geometric albedo p_V = 0.29). However, instead of providing photometric evidence for rapid nucleus rotation, our data set a lower limit to the light-curve period, P ≥ 5.4 hr. The dominant feature of the light curve is a V-shaped minimum, ∼0.3 mag deep, which is suggestive of an eclipsing binary. Under this interpretation, the time-series data are consistent with a secondary/primary mass ratio, m_s/m_p ∼ 1:6, a ratio of separation/primary radius, r/r_p ∼ 4 and an orbit period ∼0.8 days. These properties lie within the range of other asteroid binaries that are thought to be formed by rotational breakup. While the light-curve period is long, centripetal dust ejection is still possible if one or both components rotate rapidly (≲2 hr) and have small light-curve variation because of azimuthal symmetry. Indeed, radar observations of asteroids in critical rotation reveal "muffin-shaped" morphologies, which are closely azimuthally symmetric and which show minimal light curves. Our data are consistent with 311P being a close binary in which one or both components rotates near the centripetal limit. The mass loss in 2013 suggests that breakup occurred recently and could even be on-going. A search for fragments that might have been recently ejected beyond the Hill sphere reveals none larger than effective radius r_e ∼ 10 m.

CRTS J035905.9+175034 is the first eclipsing SU UMa system for which a superoutburst has been observed by Kepler in the short-cadence mode. The light curve contains one superoutburst, eight normal outbursts (including a precursor to the superoutburst), and several minioutbursts that are present before—but not after—the superoutburst. The superoutburst began with a precursor normal outburst, and shortly after the peak of the precursor, the system developed large-amplitude superhumps that achieved their maximum amplitude after just three superhump cycles. The period excess of the initial superhump period relative to the orbital period implies a mass ratio of 0.281 ± 0.015, placing it marginally above most theoretical predictions of the highest-possible mass ratio for superhump formation. In addition, our analysis of the variations in eclipse width and depth, as well as the hot spot amplitudes, generally provides substantiation of the thermal-tidal instability model. The K2 data, in conjunction with our ground-based time-resolved spectroscopy and photometry from 2014 to 2016, allows us to determine many of the fundamental parameters of this system.

The low column density gas at the outskirts of galaxies as traced by the 21 cm hydrogen line emission (H i) represents the interface between galaxies and the intergalactic medium, i.e., where galaxies are believed to get their supply of gas to fuel future episodes of star formation. Photoionization models predict a break in the radial profiles of H i at a column density of ∼5 × 10¹⁹ cm⁻² due to the lack of self-shielding against extragalactic ionizing photons. To investigate the prevalence of such breaks in galactic disks and to characterize what determines the potential edge of the H i disks, we study the azimuthally averaged H i column density profiles of 17 nearby galaxies from the H i Nearby Galaxy Survey and supplemented in two cases with published Hydrogen Accretion in LOcal GAlaxieS data. To detect potential faint H i emission that would otherwise be undetected using conventional moment map analysis, we line up individual profiles to the same reference velocity and average them azimuthally to derive stacked radial profiles. To do so, we use model velocity fields created from a simple extrapolation of the rotation curves to align the profiles in velocity at radii beyond the extent probed with the sensitivity of traditional integrated H i maps. With this method, we improve our sensitivity to outer-disk H i emission by up to an order of magnitude. Except for a few disturbed galaxies, none show evidence of a sudden change in the slope of the H i radial profiles: the alleged signature of ionization by the extragalactic background.

NGC 1232 is a face-on spiral galaxy and a great laboratory for the study of star formation due to its proximity. We obtained high spatial resolution Hα images of this galaxy, with adaptive optics, using the SAM instrument at the SOAR telescope, and used these images to study its H ii regions. These observations allowed us to produce the most complete H ii region catalog for it to date, with a total of 976 sources. This doubles the number of H ii regions previously found for this object. We used these data to construct the H ii luminosity function, and obtained a power-law index lower than the typical values found for Sc galaxies. This shallower slope is related to the presence of a significant number of high-luminosity H ii regions (log L > 39 dex). We also constructed the size distribution function, verifying that, as for most galaxies, NGC 1232 follows an exponential law. We also used the Hα luminosity to calculate the star formation rate. An extremely interesting fact about this galaxy is that X-ray diffuse observations suggest that NGC 1232 recently suffered a collision with a dwarf galaxy. We found an absence of star formation around the region where the X-ray emission is more intense, which we interpret as a star formation quenching due to the collision. Along with that, we found an excess of star-forming regions in the northeast part of the galaxy, where the X-ray emission is less intense.

The results of speckle interferometric observations at the 4.1 m Southern Astrophysical Research Telescope in 2016 and 2017 are given, totaling 2483 measurements of 1570 resolved pairs and 609 non-resolutions. We describe briefly recent changes in the instrument and observing method and quantify the accuracy of the pixel scale and position angle calibration. Comments are given on 44 pairs resolved here for the first time. The orbital motion of the newly resolved subsystem BU 83 Aa,Ab roughly agrees with its 36-year astrometric orbit proposed by J. Dommanget. Most Tycho binaries examined here turned out to be spurious.

With the advent of more and deeper sky surveys, the discovery of interstellar small objects entering into the solar system has been finally possible. In 2017 October 19, using observations of the Pan-STARRS survey, a fast moving object, now officially named 1I/2017 U1 ('Oumuamua), was discovered in a heliocentric unbound trajectory, suggesting an interstellar origin. Assessing the provenance of interstellar small objects is key for understanding their distribution, spatial density, and the processes responsible for their ejection from planetary systems. However, their peculiar trajectories place a limit on the number of observations available to determine a precise orbit. As a result, when its position is propagated ∼10⁵–10⁶ years backward in time, small errors in orbital elements become large uncertainties in position in the interstellar space. In this paper we present a general method for assigning probabilities to nearby stars of being the parent system of an observed interstellar object. We describe the method in detail and apply it for assessing the origin of 'Oumuamua. A preliminary list of potential progenitors and their corresponding probabilities is provided. In the future, when further information about the object and/or the nearby stars be refined, the probabilities computed with our method can be updated. We provide all the data and codes we developed for this purpose in the form of an open source C/C++/Python package, iWander, which is publicly available at http://github.com/seap-udea/iWander.

Obliquity variability could play an important role in the climate and habitability of a planet. Orbital modulations caused by planetary companions and the planet's spin axis precession due to the torque from the host star may lead to resonant interactions and cause large-amplitude obliquity variability. Here we consider the spin axis dynamics of Kepler-62f and Kepler-186f, both of which reside in the habitable zone around their host stars. Using N-body simulations and secular numerical integrations, we describe their obliquity evolution for particular realizations of the planetary systems. We then use a generalized analytic framework to characterize regions in parameter space where the obliquity is variable with large amplitude. We find that the locations of variability are fine-tuned over the planetary properties and system architecture in the lower-obliquity regimes (≲40°). As an example, assuming a rotation period of 24 hr, the obliquities of both Kepler-62f and Kepler-186f are stable below ∼40°, whereas the high-obliquity regions (60°–90°) allow moderate variabilities. However, for some other rotation periods of Kepler-62f or Kepler-186f, the lower-obliquity regions could become more variable owing to resonant interactions. Even small deviations from coplanarity (e.g., mutual inclinations ∼3°) could stir peak-to-peak obliquity variations up to ∼20°. Undetected planetary companions and/or the existence of a satellite could also destabilize the low-obliquity regions. In all cases, the high-obliquity region allows for moderate variations, and all obliquities corresponding to retrograde motion (i.e., >90°) are stable.

We present Spitzer Space Telescope time-series photometry at 3.6 and 4.5 μm of 2MASS J11193254−1137466AB and WISEA J114724.10−204021.3, two planetary-mass, late-type (∼L7) brown dwarf members of the ∼10 Myr old TW Hya Association. These observations were taken in order to investigate whether or not a tentative trend of increasing variability amplitude with decreasing surface gravity seen for L3–L5.5 dwarfs extends to later-L spectral types and to explore the angular momentum evolution of low-mass objects. We examine each light curve for variability and find a rotation period of 19.39^+0.33_−0.28 hr and semi-amplitudes of 0.798^+0.081_−0.083% at 3.6 μm and 1.108^+0.093_−0.094% at 4.5 μm for WISEA J114724.10−204021.3. For 2MASS J11193254−1137466AB, we find a single period of 3.02^+0.04_−0.03 hr with semi-amplitudes of 0.230^+0.036_−0.035% at 3.6 μm and 0.453 ± 0.037% at 4.5 μm, which we find is possibly due to the rotation of one component of the binary. Combining our results with 12 other late-type L dwarfs observed with Spitzer from the literature, we find no significant differences between the 3.6 μm amplitudes of low surface gravity and field gravity late-type L brown dwarfs at Spitzer wavelengths, and find tentative evidence (75% confidence) of higher amplitude variability at 4.5 μm for young, late-type Ls. We also find a median rotation period of young brown dwarfs (10–300 Myr) of ∼10 hr, more than twice the value of the median rotation period of field-age brown dwarfs (∼4 hr), a clear signature of brown dwarf rotational evolution.

We derive lower bounds on the orbital distance and inclination of a putative planet beyond the transiting seven planets of TRAPPIST-1, for a range of masses ranging from 0.08 M_Jup to 3.5 M_Jup. While the outer architecture of this system will ultimately be constrained by radial velocity measurements over time, we present dynamical constraints from the remarkably coplanar configuration of the seven transiting planets, which is sensitive to modestly inclined perturbers. We find that the observed configuration is unlikely if a Jovian-mass planet inclined by ≥3° to the transiting planet exists within 0.53 au, exceeding any constraints from transit timing variations (TTV) induced in the known planets from an undetected perturber. Our results will inform RV programs targeting TRAPPIST-1, and for near coplanar outer planets, tighter constraints are anticipated for radial velocity (RV) precisions of ≲140 m s⁻¹. At higher inclinations, putative planets are ruled out to greater orbital distances with orbital periods up to a few years.

A metallicity, chemical composition, and kinematic survey has been conducted for a sample of 340 candidate field red horizontal-branch (RHB) stars. Spectra with high resolution and high signal-to-noise ratio were gathered with the McDonald Observatory 2.7 m Tull and the Hobby–Eberly Telescope echelle spectrographs, and were used to determine effective temperatures, surface gravities, microturbulent velocities, [Fe/H] metallicities, and abundance ratios [X/Fe] for seven α and Fe-group species. The derived temperatures and gravities confirm that at least half of the candidates are true RHB stars, with (average) parameters T_eff ∼ 5000 K and log g ∼ 2.5. From the α abundances alone, the thin and thick Galactic populations are apparent in our sample. Space motions for 90% of the program stars were computed from Hipparcos and Gaia parallaxes and proper motions. Correlations between chemical compositions and Galactic kinematics clearly indicate the existence of both thin-disk and thick-disk RHB stars.

We present a Chandra study of the deeply embedded Serpens South star-forming region, examining cluster structure and disk properties at the earliest stages. In total, 152 X-ray sources are detected. Combined with Spitzer and 2MASS photometry, 66 X-ray sources are reliably matched to an IR counterpart. We identify 21 class I, 6 flat spectrum, 16 class II, and 18 class III young stars; 5 were unclassified. Eighteen sources were variable in X-rays, 8 exhibiting flare-like emission and one source being periodic. The cluster's X-ray luminosity distance was estimated: the best match was to the nearer distance of 260 pc for the front of the Aquila Rift complex. The ratio of N_H to A_K is found to be ∼0.68 × 10²², similar to that measured in other young low-mass regions, but lower than that measured in the interstellar medium and high-mass clusters (∼(1.6–2) × 10²²). We find that the spatial distribution closely follows that of the dense filament from which the stars have formed, with the class II population still strongly associated with the filament. There are four subclusters in the field, with three forming knots in the filament, and a fourth to the west, which may not be associated but may be contributing to the distributed class III population. A high percentage of diskless class IIIs (upper limit 30% of classified X-ray sources) in such a young cluster could indicate that processing of disks is influenced by the cluster environment and is not solely dependent on timescale.

Kiefer et al. reported the detection of infalling Ca ii absorption in HD 172555, a member of the β Pictoris Moving Group (βPMG). We obtained HST Space Telescope Imaging Spectrograph and Cosmic Origins Spectrograph spectroscopy of this star at 2 epochs separated by a week, and we report the discovery of infalling gas in resonant transitions of Si iii and iv, C ii and iv, and neutral atomic oxygen. Variable absorption is seen in the C ii transitions and is optically thick, with covering factors which range between 58% and 68%, similar to features seen in β Pictoris. The O i spectral profile resembles that of C ii, showing a strong low-velocity absorption to +50 km s⁻¹ in the single spectral segment obtained during orbital night, as well as what may be higher-velocity absorption. Studies of the mid-IR spectrum of this system have suggested the presence of silica. The O i absorption differs from that seen in Si iii, suggesting that the neutral atomic oxygen does not originate in SiO dissociation products but in a more volatile parent molecule such as CO.

Two planetary mass objects in the far outer solar system—collectively referred to here as Planet X— have recently been hypothesized to explain the orbital distribution of distant Kuiper Belt Objects. Neither planet is thought to be exceptionally faint, but the sky locations of these putative planets are poorly constrained. Therefore, a wide area survey is needed to detect these possible planets. The Large Synoptic Survey Telescope (LSST) will carry out an unbiased, large area (around 18000 deg²), deep (limiting magnitude of individual frames of 24.5) survey (the "wide-fast-deep (WFD)" survey) of the southern sky beginning in 2022, and it will therefore be an important tool in searching for these hypothesized planets. Here, we explore the effectiveness of LSST as a search platform for these possible planets. Assuming the current baseline cadence (which includes the WFD survey plus additional coverage), we estimate that LSST will confidently detect or rule out the existence of Planet X in 61% of the entire sky. At orbital distances up to ∼75 au, Planet X could simply be found in the normal nightly moving object processing; at larger distances, it will require custom data processing. We also discuss the implications of a nondetection of Planet X in LSST data.

Wide-field surveys for transiting planets, such as the NASA Kepler and TESS missions, are usually conducted without knowing which stars have binary companions. Unresolved and unrecognized binaries give rise to systematic errors in planet occurrence rates, including misclassified planets and mistakes in completeness corrections. The individual errors can have different signs, making it difficult to anticipate the net effect on inferred occurrence rates. Here, we use simplified models of signal-to-noise limited transit surveys to try and clarify the situation. We derive a formula for the apparent occurrence rate density measured by an observer who falsely assumes all stars are single. The formula depends on the binary fraction, the mass function of the secondary stars, and the true occurrence of planets around primaries, secondaries, and single stars. It also takes into account the Malmquist bias by which binaries are over-represented in flux-limited samples. Application of the formula to an idealized Kepler-like survey shows that for planets larger than 2 R_⊕, the net systematic error is of order 5%. In particular, unrecognized binaries are unlikely to be the reason for the apparent discrepancies between hot-Jupiter occurrence rates measured in different surveys. For smaller planets the errors are potentially larger: the occurrence of Earth-sized planets could be overestimated by as much as 50%. We also show that whenever high-resolution imaging reveals a transit host star to be a binary, the planet is usually more likely to orbit the primary star than the secondary star.

We present new observations from 2016 of two Jupiter Trojan asteroids that are targets for the Lucy Discovery mission. The extremely long rotation period of (11351) Leucus is confirmed and refined to a secure value of 445.732 ± 0.021 hr with photometric parameters of H_r = 11.046 ± 0.003 and G_r = 0.58 ± 0.02 in the SDSS r' filter. This leads to a geometric albedo of p_V = 4.7%. The amplitude of the light curve was measured to be 0.61 mag, unchanged from the value of one-fourth of a revolution earlier, suggesting a low obliquity. The first light-curve observations for (15094) Polymele are also presented. This object is revealed to have a much shorter rotation period of 5.8607 ± 0.0005 hr with a very low amplitude of 0.09 mag. Its photometric parameters are H_r = 11.691 ± 0.002 and G_r = 0.22 ± 0.02. These values lead to a refined geometric albedo of p_V = 7.3%. This object is either nearly spherical or was being viewed nearly pole-on in 2016. Further observations are required to fully determine the spin pole orientation and convex-hull shapes.

The Rosetta spacecraft observations revealed that the nucleus of comet 67P/Churyumov–Gerasimenko consists of two similarly sized lobes connected by a narrow neck. Here, we evaluate the possibility that 67P is a collapsed binary. We assume that the progenitor of 67P was a binary and consider various physical mechanisms that could have brought the binary components together, including small-scale impacts and gravitational encounters with planets. We find that 67P could be a primordial body (i.e., not a collisional fragment) if the outer planetesimal disk lasted ≲10 Myr before it was dispersed by migrating Neptune. The probability of binary collapse by impact is ≃30% for tightly bound binaries. Most km-class binaries become collisionally dissolved. Roughly 10% of the surviving binaries later evolve to become contact binaries during the disk dispersal, when bodies suffer gravitational encounters with Neptune. Overall, the processes described in this work do not seem to be efficient enough to explain the large fraction (∼67%) of bi-lobed cometary nuclei inferred from spacecraft imaging.

We conducted time-resolved optical spectroscopy and/or time-series photometry of 15 cataclysmic binaries that were discovered in hard X-ray surveys by the Swift Burst Alert Telescope and the International Gamma-Ray Astrophysics Laboratory, with the goal of measuring their orbital periods and searching for spin periods. Four of the objects in this study are new optical identifications: Swift J0535.2+2830, Swift J2006.4+3645, IGR J21095+4322, and Swift J2116.5+5336. Coherent pulsations are detected from three objects for the first time, Swift J0535.2+2830 (1523 s), 2PBC J1911.4+1412 (747 s), and 1SWXRT J230642.7+550817 (464 s), indicating that they are intermediate polars (IPs). We find two new eclipsing systems in time-series photometry: 2PBC J0658.0−1746, a polar with a period of 2.38 hr, and Swift J2116.5+5336, a disk system that has an eclipse period of 6.56 hr. Exact or approximate spectroscopic orbital periods are found for six additional targets. Of note is the long 4.637-day orbit for Swift J0623.9−0939, which is revealed by the radial velocities of the photospheric absorption lines of the secondary star. We also discover a 12.76 hr orbital period for RX J2015.6+3711, which confirms that the previously detected 2.00 hr X-ray period from this star is the spin period of an IP, as inferred by Coti Zelati et al. These results support the conclusion that hard X-ray selection favors magnetic CVs, with IPs outnumbering polars.

We observed 12 Plutinos over two separated years with the 4.3 m Lowell's Discovery Channel Telescope. Here, we present the first light-curve data for those objects. Three of them (2014 JL₈₀, 2014 JO₈₀, and 2014 JQ₈₀) display a large light-curve amplitude explainable by a single elongated object, but they are most likely caused by a contact binary system due to their light-curve morphology. These potential contact binaries have rotational periods from 6.3 to 34.9 hr and peak-to-peak light-curve variability between 0.6 and 0.8 mag. We present partial light curves, allowing us to constrain the light-curve amplitude and the rotational period of another nine Plutinos. By merging our data with the literature, we estimate that up to ∼40% of the Plutinos could be contact binaries. Interestingly, we found that all of the suspected contact binaries in the 3:2 resonance are small with absolute magnitude H > 6 mag. Based on our sample and the literature, up to ∼50% of the small Plutinos are potential contact binaries.

Observations of clustering among the orbits of the most distant trans-Neptunian objects (TNOs) has inspired interest in the possibility of an undiscovered ninth planet lurking in the outskirts of the solar system. Numerical simulations by a number of authors have demonstrated that, with appropriate choices of planet mass and orbit, such a planet can maintain clustering in the orbital elements of the population of distant TNOs, similar to the observed sample. However, many aspects of the rich underlying dynamical processes induced by such a distant eccentric perturber have not been fully explored. We report the results of our investigation of the dynamics of coplanar test-particles that interact with a massive body on an circular orbit (Neptune) and a massive body on a more distant, highly eccentric orbit (the putative Planet Nine). We find that a detailed examination of our idealized simulations affords tremendous insight into the rich test-particle dynamics that are possible. In particular, we find that chaos and resonance overlap plays an important role in particles' dynamical evolution. We develop a simple mapping model that allows us to understand, in detail, the web of overlapped mean-motion resonances explored by chaotically evolving particles. We also demonstrate that gravitational interactions with Neptune can have profound effects on the orbital evolution of particles. Our results serve as a starting point for a better understanding of the dynamical behavior observed in more complicated simulations that can be used to constrain the mass and orbit of Planet Nine.

The observation that the orbits of long-period Kuiper Belt objects (KBOs) are anomalously clustered in physical space has recently prompted the Planet Nine hypothesis—the proposed existence of a distant and eccentric planetary member of our solar system. Within the framework of this model, a Neptune-like perturber sculpts the orbital distribution of distant KBOs through a complex interplay of resonant and secular effects, such that in addition to perihelion-circulating objects, the surviving orbits get organized into apsidally aligned and anti-aligned configurations with respect to Planet Nine's orbit. In this work, we investigate the role of Kuiper Belt initial conditions on the evolution of the outer solar system using numerical simulations. Intriguingly, we find that the final perihelion distance distribution depends strongly on the primordial state of the system, and we demonstrate that a bimodal structure corresponding to the existence of both aligned and anti-aligned clusters is only reproduced if the initial perihelion distribution is assumed to extend well beyond ∼36 au. The bimodality in the final perihelion distance distribution is due to the existence of permanently stable objects, with the lower perihelion peak corresponding to the anti-aligned orbits and the higher perihelion peak corresponding to the aligned orbits. We identify the mechanisms that enable the persistent stability of these objects and locate the regions of phase space in which they reside. The obtained results contextualize the Planet Nine hypothesis within the broader narrative of solar system formation and offer further insight into the observational search for Planet Nine.

Titan's atmospheric inventory of oxygen compounds (H₂O, CO₂, CO) are thought to result from photochemistry acting on externally supplied oxygen species (O⁺, OH, H₂O). These species potentially originate from two main sources: (1) cryogenic plumes from the active moon Enceladus and (2) micrometeoroid ablation. Enceladus is already suspected to be the major O⁺ source, which is required for CO creation. However, photochemical models also require H₂O and OH influx to reproduce observed quantities of CO₂ and H₂O. Here, we exploit sulphur as a tracer to investigate the oxygen source because it has very different relative abundances in micrometeorites (S/O ∼ 10⁻²) and Enceladus' plumes (S/O ∼ 10⁻⁵). Photochemical models predict most sulphur is converted to CS in the upper atmosphere, so we use Atacama Large Millimeter/submillimeter Array (ALMA) observations at ∼340 GHz to search for CS emission. We determined stringent CS 3σ stratospheric upper limits of 0.0074 ppb (uniform above 100 km) and 0.0256 ppb (uniform above 200 km). These upper limits are not quite stringent enough to distinguish between Enceladus and micrometeorite sources at the 3σ level and a contribution from micrometeorites cannot be ruled out, especially if external flux is toward the lower end of current estimates. Only the high-flux micrometeorite source model of Hickson et al. can be rejected at 3σ. We determined a 3σ stratospheric upper limit for CH₂NH of 0.35 ppb, which suggests cosmic rays may have a smaller influence in the lower stratosphere than predicted by some photochemical models. Disk-averaged C₃H₄ and C₂H₅CN profiles were determined and are consistent with previous ALMA and Cassini/CIRS measurements.

Parallaxes are presented for a sample of 20 nearby dwarf carbon stars. The inferred luminosities cover almost two orders of magnitude. Their absolute magnitudes and tangential velocities confirm prior expectations that some originate in the Galactic disk, although more than half of this sample are halo stars. Three stars are found to be astrometric binaries, and orbital elements are determined; their semimajor axes are 1–3 au, consistent with the size of an AGB mass-transfer donor star.

We develop a new numerical scheme for solving the radiative transfer equation in a spherically symmetric system. This scheme does not rely on any kind of diffusion approximation, and it is accurate for optically thin, thick, and intermediate systems. In the limit of a homogeneously distributed extinction coefficient, our method is very accurate and exceptionally fast. We combine this fast method with a slower but more generally applicable method to describe realistic problems. We perform various test calculations, including a simplified protostellar collapse simulation. We also discuss possible future improvements.

4C+15.05 (also known as NRAO 91, PKS 0202+14, or J0204+15) is a quintessential blazar. It has a luminous, variable radio spectrum, a super-luminal jet, and gamma-ray detections. Arecibo observations with the 700–800 MHz receiver on the 305 m diameter William E. Gordon Telescope detected, serendipitously, H i in absorption against 4C+15.05 while using it as a bandpass calibrator for another object in an H i absorption project. Although the redshift we derive is different from that commonly in use in the literature (nominally z = 0.405), it agrees very well with the value of z = 0.833 determined by Stickel et al. This absorption feature is best fitted by a sum of three Gaussians, which yield an average redshift of z = 0.8336 ± 0.0004, although without corresponding high-resolution imaging it is not possible to say whether the components are parts of outflows or inflows. A total column density of N(H i) = 2.39 ± 0.13 × 10²¹ cm⁻² is derived, relatively high compared to many radio-loud sources. These results are compared to various relationships in the literature.

HAT-P-11 is a mid-K dwarf that hosts one of the first Neptune-sized planets found outside the solar system. The orbit of HAT-P-11b is misaligned with the star's spin—one of the few known cases of a misaligned planet orbiting a star less massive than the Sun. We find an additional planet in the system based on a decade of precision radial velocity (RV) measurements from Keck/High Resolution Echelle Spectrometer. HAT-P-11c is similar to Jupiter in its mass (${M}_{P}\sin i=1.6\pm 0.1$M_J) and orbital period ($P={9.3}_{-0.5}^{+1.0}$ year), but has a much more eccentric orbit (e = 0.60 ± 0.03). In our joint modeling of RV and stellar activity, we found an activity-induced RV signal of ∼7 ${\rm{m}}\,{{\rm{s}}}^{-1}$, consistent with other active K dwarfs, but significantly smaller than the 31 ${\rm{m}}\,{{\rm{s}}}^{-1}$ reflex motion due to HAT-P-11c. We investigated the dynamical coupling between HAT-P-11b and c as a possible explanation for HAT-P-11b's misaligned orbit, finding that planet–planet Kozai interactions cannot tilt planet b's orbit due to general relativistic precession; however, nodal precession operating on million year timescales is a viable mechanism to explain HAT-P-11b's high obliquity. This leaves open the question of why HAT-P-11c may have such a tilted orbit. At a distance of 38 pc, the HAT-P-11 system offers rich opportunities for further exoplanet characterization through astrometry and direct imaging.

We present results from a medium-resolution (R ∼ 2000) spectroscopic follow-up campaign of 1694 bright (V < 13.5), very metal-poor star candidates from the RAdial Velocity Experiment (RAVE). Initial selection of the low-metallicity targets was based on the stellar parameters published in RAVE Data Releases 4 and 5. Follow up was accomplished with the Gemini-N and Gemini-S, the ESO/NTT, the KPNO/Mayall, and the SOAR telescopes. The wavelength coverage for most of the observed spectra allows for the determination of carbon and α-element abundances, which are crucial for considering the nature and frequency of the carbon-enhanced metal-poor (CEMP) stars in this sample. We find that 88% of the observed stars have $[\mathrm{Fe}/{\rm{H}}]$ ≤ −1.0, 61% have $[\mathrm{Fe}/{\rm{H}}]$ ≤ −2.0, and 3% have $[\mathrm{Fe}/{\rm{H}}]$ ≤ −3.0 (with four stars at $[\mathrm{Fe}/{\rm{H}}]$ ≤ −3.5). There are 306 CEMP star candidates in this sample, and we identify 169 CEMP Group I, 131 CEMP Group II, and 6 CEMP Group III stars from the A(C) versus [Fe/H] diagram. Inspection of the $[\alpha /{\rm{C}}]$ abundance ratios reveals that five of the CEMP Group II stars can be classified as "mono-enriched second-generation" stars. Gaia DR1 matches were found for 734 stars, and we show that transverse velocities can be used as a confirmatory selection criteria for low-metallicity candidates. Selected stars from our validated list are being followed-up with high-resolution spectroscopy to reveal their full chemical-abundance patterns for further studies.

K2-18 is a nearby M2.5 dwarf, located at 34 pc and hosting a transiting planet that was first discovered by the K2 mission and later confirmed with Spitzer Space Telescope observations. With a radius of ∼2 R_⊕ and an orbital period of ∼33 days, the planet lies in the temperate zone of its host star and receives stellar irradiation similar to that of Earth. Here we perform radial velocity follow-up observations with the visual channel of CARMENES with the goal of determining the mass and density of the planet. We measure a planetary semi-amplitude of K_b ∼ 3.5 ${\rm{m}}\,{{\rm{s}}}^{-1}$ and a mass of M_b ∼ 9 M_⊕, yielding a bulk density around ${\rho }_{b}\sim 4\,{\rm{g}}\,{\mathrm{cm}}^{-3}$. This indicates a low-mass planet with a composition consistent with a solid core and a volatile-rich envelope. A signal at 9 days was recently reported using radial velocity measurements taken with the HARPS spectrograph. This was interpreted as being due to a second planet. We see a weaker, time- and wavelength-dependent signal in the CARMENES data set and thus favor stellar activity for its origin. K2-18 b joins the growing group of low-mass planets detected in the temperate zone of M dwarfs. The brightness of the host star in the near-infrared makes the system a good target for detailed atmospheric studies with the James Webb Space Telescope.

The brighter-fatter effect has been postulated to arise due to the build up of a transverse electric field, produced as photocharges accumulate in the pixels' potential wells. We investigate the brighter-fatter effect in the Hyper Suprime-Cam by examining flat fields and moments of stars. We observe deviations from the expected linear relation in the photon transfer curve (PTC), luminosity-dependent correlations between pixels in flat-field images, and a luminosity-dependent point-spread function (PSF) in stellar observations. Under the key assumptions of translation invariance and Maxwell's equations in the quasi-static limit, we give a first-principles proof that the effect can be parameterized by a translationally invariant scalar kernel. We describe how this kernel can be estimated from flat fields and discuss how this kernel has been used to remove the brighter-fatter distortions in Hyper Suprime-Cam images. We find that our correction restores the expected linear relation in the PTCs and significantly reduces, but does not completely remove, the luminosity dependence of the PSF over a wide range of magnitudes.

We show that dense OGLE and KMTNet I-band survey data require four bodies (sources plus lenses) to explain the microlensing light curve of OGLE-2015-BLG-1459. However, these can equally well consist of three lenses and one source (3L1S), two lenses and two sources (2L2S), or one lens and three sources (1L3S). In the 3L1S and 2L2S interpretations, the host is a brown dwarf and the dominant companion is a Neptune-class planet, with the third body (in the 3L1S case) being a Mars-class object that could have been a moon of the planet. In the 1L3S solution, the light curve anomalies are explained by a tight (five stellar radii) low-luminosity binary source that is offset from the principal source of the event by $\sim 0.17\,\mathrm{au}$. These degeneracies are resolved in favor of the 1L3S solution by color effects derived from comparison to MOA data, which are taken in a slightly different (R/I) passband. To enable current and future (WFIRST) surveys to routinely characterize exo-moons and distinguish among such exotic systems requires an observing strategy that includes both a cadence faster than 9 minute⁻¹ and observations in a second band on a similar timescale.

We discuss the detection in the Outer Solar System Origins Survey (OSSOS) of two objects in Neptune's distant 9:1 mean motion resonance at semimajor axis a ≈ 130 au. Both objects are securely resonant on 10 Myr timescales, with one securely in the 9:1 resonance's leading asymmetric libration island and the other in either the symmetric or trailing asymmetric island. These objects are the largest semimajor axis objects with secure resonant classifications, and their detection in a carefully characterized survey allows for the first robust resonance population estimate beyond 100 au. The detection of these objects implies a 9:1 resonance population of 1.1 × 10⁴ objects with H_r < 8.66 (D ≳ 100 km) on similar orbits (95% confidence range of ∼(0.4–3) × 10⁴). Integrations over 4 Gyr of an ensemble of clones spanning these objects' orbit-fit uncertainties reveal that they both have median resonance occupation timescales of ∼1 Gyr. These timescales are consistent with the hypothesis that these objects originate in the scattering population but became transiently stuck to Neptune's 9:1 resonance within the last ∼1 Gyr of solar system evolution. Based on simulations of a model of the current scattering population, we estimate the expected resonance sticking population in the 9:1 resonance to be 1000–4500 objects with H_r < 8.66; this is marginally consistent with the OSSOS 9:1 population estimate. We conclude that resonance sticking is a plausible explanation for the observed 9:1 population, but we also discuss the possibility of a primordial 9:1 population, which would have interesting implications for the Kuiper Belt's dynamical history.

We analyze the combined Spitzer and ground-based data for OGLE-2017-BLG-1140 and show that the event was generated by a Jupiter-class $({m}_{p}\simeq 1.6\,{M}_{{\rm{J}}{\rm{u}}{\rm{p}}})$ planet orbiting a mid-late M dwarf $(M\simeq 0.2\,{M}_{\odot })$ that lies ${D}_{\mathrm{LS}}\simeq 1.0\,\mathrm{kpc}$ in the foreground of the microlensed Galactic-bar source star. The planet–host projected separation is ${a}_{\perp }\simeq 1.0\,\mathrm{au}$, i.e., well beyond the snow line. By measuring the source proper motion ${{\boldsymbol{\mu }}}_{s}$ from ongoing long-term OGLE imaging and combining this with the lens-source relative proper motion ${{\boldsymbol{\mu }}}_{\mathrm{rel}}$ derived from the microlensing solution, we show that the lens proper motion ${{\boldsymbol{\mu }}}_{l}={{\boldsymbol{\mu }}}_{\mathrm{rel}}+{{\boldsymbol{\mu }}}_{s}$ is consistent with the lens lying in the Galactic disk, although a bulge lens is not ruled out. We show that while the Spitzer and ground-based data are comparably well fitted by planetary (i.e., binary-lens (2L1S)) and binary-source (1L2S) models, the combination of Spitzer and ground-based data decisively favors the planetary model. This is a new channel to resolve the 2L1S/1L2S degeneracy, which can be difficult to break in some cases.

Asteroids in mean motion resonances with giant planets are common in the solar system, but it was not until recently that several asteroids in retrograde mean motion resonances with Jupiter and Saturn were discovered. A retrograde co-orbital asteroid of Jupiter, 2015 BZ509 is confirmed to be in a long-term stable retrograde 1:1 mean motion resonance with Jupiter, which gives rise to our interests in its unique resonant dynamics. In this paper, we investigate the phase-space structure of the retrograde 1:1 resonance in detail within the framework of the circular restricted three-body problem. We construct a simple integrable approximation for the planar retrograde resonance using canonical contact transformation and numerically employ the averaging procedure in closed form. The phase portrait of the retrograde 1:1 resonance is depicted with the level curves of the averaged Hamiltonian. We thoroughly analyze all possible librations in the co-orbital region and uncover a new apocentric libration for the retrograde 1:1 resonance inside the planet's orbit. We also observe the significant jumps in orbital elements for outer and inner apocentric librations, which are caused by close encounters with the perturber.

We report on the analysis of a microlensing event, OGLE-2014-BLG-1722, that showed two distinct short-term anomalies. The best-fit model to the observed light curves shows that the two anomalies are explained with two planetary mass ratio companions to the primary lens. Although a binary-source model is also able to explain the second anomaly, it is marginally ruled out by 3.1σ. The two-planet model indicates that the first anomaly was caused by planet "b" with a mass ratio of $q=({4.5}_{-0.6}^{+0.7})\times {10}^{-4}$ and projected separation in units of the Einstein radius, s = 0.753 ± 0.004. The second anomaly reveals planet "c" with a mass ratio of ${q}_{2}=({7.0}_{-1.7}^{+2.3})\times {10}^{-4}$ with Δχ² ∼ 170 compared to the single-planet model. Its separation has two degenerated solutions: the separation of planet c is s₂ = 0.84 ± 0.03 and 1.37 ± 0.04 for the close and wide models, respectively. Unfortunately, this event does not show clear finite-source and microlensing parallax effects; thus, we estimated the physical parameters of the lens system from Bayesian analysis. This gives the masses of planets b and c as ${m}_{{\rm{b}}}={56}_{-33}^{+51}\,$ and ${m}_{{\rm{c}}}={85}_{-51}^{+86}\,{M}_{\oplus }$, respectively, and they orbit a late-type star with a mass of ${M}_{\mathrm{host}}\,={0.40}_{-0.24}^{+0.36}\,{M}_{\odot }$ located at ${D}_{{\rm{L}}}={6.4}_{-1.8}^{+1.3}\,\mathrm{kpc}$ from us. The projected distances between the host and planets are ${r}_{\perp ,{\rm{b}}}=1.5\pm 0.6\,\mathrm{au}$ for planet b and ${r}_{\perp ,{\rm{c}}}={1.7}_{-0.6}^{+0.7}\,\mathrm{au}$ and ${r}_{\perp ,{\rm{c}}}={2.7}_{-1.0}^{+1.1}\,\mathrm{au}$ for the close and wide models of planet c. If the two-planet model is true, then this is the third multiple-planet system detected using the microlensing method and the first multiple-planet system detected in low-magnification events, which are dominant in the microlensing survey data. The occurrence rate of multiple cold gas giant systems is estimated using the two such detections and a simple extrapolation of the survey sensitivity of the 6 yr MOA microlensing survey combined with the 4 yr μFUN detection efficiency. It is estimated that 6% ± 2% of stars host two cold giant planets.

We present a new computer program, SRTC++, to solve spatial problems associated with explorations of Saturn's moon Titan. The program implements a three-dimensional structure well-suited to addressing shortcomings arising from plane-parallel radiative transfer approaches. SRTC++'s design uses parallel processing in an object-oriented, compiled computer language (C++) leading to a flexible and fast architecture. We validate SRTC++ using analytical results, semianalytical radiative transfer expressions, and an existing Titan plane-parallel model. SRTC++ complements existing approaches, addressing spatial problems like near-limb and near-terminator geometries, non-Lambertian surface phase functions (including specular reflections), and surface albedo nonuniformity.

We describe the 44 systems discovered to be within 10 pc of the Sun by the RECONS team, primarily via the long-term astrometry program at the CTIO/SMARTS 0.9 m that began in 1999. The systems—including 41 with red dwarf primaries, 2 white dwarfs, and 1 brown dwarf—have trigonometric parallaxes greater than 100 mas, with errors of 0.4–2.4 mas in all but one case. We provide updated astrometric, photometric (VRIJHK magnitudes), spectral type, and multiplicity information here. Among these are 14 systems that are new entries to the 10 pc sample, including the first parallaxes for 9 systems and new values for 5 systems that had previous parallaxes with errors greater than 10 mas or values placing them beyond 10 pc. We also provide new data for 22 systems known to lie within 10 pc and 9 systems reported to be closer than that horizon but for which new parallaxes place them further away, bringing the total to 75 systems. The 44 systems added by RECONS comprise one of every 7 systems known within 10 pc. We illustrate the evolution of the 10 pc sample from the 191 systems known when the final Yale Parallax Catalog was published in 1995 to the 317 systems known today. Even so close to the Sun, additional discoveries of white, red, and brown dwarfs are possible, both as primaries and secondaries, although we estimate that at least 90% of the stellar systems closer than 10 pc have now been identified.

Using an energy balance model with ice sheets, we examine the climate response of an Earth-like planet orbiting a G-dwarf star and experiencing large orbital and obliquity variations. We find that ice caps couple strongly to the orbital forcing, leading to extreme ice ages. In contrast with previous studies, we find that such exo-Milankovitch cycles tend to impair habitability by inducing snowball states within the habitable zone. The large amplitude changes in obliquity and eccentricity cause the ice edge, the lowest-latitude extent of the ice caps, to become unstable and grow to the equator. We apply an analytical theory of the ice edge latitude to show that obliquity is the primary driver of the instability. The thermal inertia of the ice sheets and the spectral energy distribution of the G-dwarf star increase the sensitivity of the model to triggering runaway glaciation. Finally, we apply a machine learning algorithm to demonstrate how this technique can be used to extend the power of climate models. This work illustrates the importance of orbital evolution for habitability in dynamically rich planetary systems. We emphasize that as potentially habitable planets are discovered around G dwarfs, we need to consider orbital dynamics.