Numerical simulations suggest that kink and torus instabilities are two potential contributors to the initiation and prorogation of eruptive events. A magnetic parameter called the decay index (i.e., the coronal magnetic gradient of the overlying fields above the eruptive flux ropes) could play an important role in controlling the kinematics of eruptions. Previous studies have identified a threshold range of the decay index that distinguishes between eruptive and confined configurations. Here we advance the study by investigating if there is a clear correlation between the decay index and coronal mass ejection (CME) speed. Thirty-eight CMEs associated with filament eruptions and/or two-ribbon flares are selected using the Hα data from the Global Hα Network. The filaments and flare ribbons observed in Hα associated with the CMEs help to locate the magnetic polarity inversion line, along which the decay index is calculated based on the potential field extrapolation using Michelson Doppler Imager magnetograms as boundary conditions. The speeds of CMEs are obtained from the LASCO C2 CME catalog available online. We find that the mean decay index increases with CME speed for those CMEs with a speed below 1000 km s ^–1 and stays flat around 2.2 for the CMEs with higher speeds. In addition, we present a case study of a partial filament eruption, in which the decay indices show different values above the erupted/non-erupted part.

We introduce a simple model to explain the time-reversed and stretched residuals in gamma-ray burst (GRB) pulse light curves. In this model an impactor wave in an expanding GRB jet accelerates from subluminal to superluminal velocities, or decelerates from superluminal to subluminal velocities. The impactor wave interacts with the surrounding medium to produce Cerenkov and/or other collisional radiation when traveling faster than the speed of light in this medium, and other mechanisms (such as thermalized Compton or synchrotron shock radiation) when traveling slower than the speed of light. These transitions create both a time-forward and a time-reversed set of light-curve features through the process of relativistic image doubling. The model can account for a variety of unexplained yet observed GRB pulse behaviors, including the amount of stretching observed in time-reversed GRB pulse residuals and the relationship between stretching factor and pulse asymmetry. The model is applicable to all GRB classes since similar pulse behaviors are observed in long/intermediate GRBs, short GRBs, and X-ray flares. The free model parameters are the impactor’s Lorentz factor when moving subluminally, its Lorentz factor when moving superluminally, and the speed of light in the impacted medium.

The Event Horizon Telescope (EHT) has recently delivered the first resolved images of M87*, the supermassive black hole in the center of the M87 galaxy. These images were produced using 230 GHz observations performed in 2017 April. Additional observations are required to investigate the persistence of the primary image feature—a ring with azimuthal brightness asymmetry—and to quantify the image variability on event horizon scales. To address this need, we analyze M87* data collected with prototype EHT arrays in 2009, 2011, 2012, and 2013. While these observations do not contain enough information to produce images, they are sufficient to constrain simple geometric models. We develop a modeling approach based on the framework utilized for the 2017 EHT data analysis and validate our procedures using synthetic data. Applying the same approach to the observational data sets, we find the M87* morphology in 2009–2017 to be consistent with a persistent asymmetric ring of ∼40 μas diameter. The position angle of the peak intensity varies in time. In particular, we find a significant difference between the position angle measured in 2013 and 2017. These variations are in broad agreement with predictions of a subset of general relativistic magnetohydrodynamic simulations. We show that quantifying the variability across multiple observational epochs has the potential to constrain the physical properties of the source, such as the accretion state or the black hole spin.

The late-stage evolution of the most massive stars such as η Carinae is controlled by the effects of mass loss, which may be dominated by poorly understood eruptive mass ejections. Understanding this population is challenging because no true analogs of η Car have been clearly identified in the Milky Way or other galaxies. We utilize Spitzer IRAC images of seven nearby ( 4 Mpc) galaxies to search for such analogs. We find 34 candidates with a flat or rising mid-IR spectral energy distributions toward longer mid-infrared wavelengths that emit >10 ⁵ L _☉ in the IRAC bands (3.6 to 8.0 μm) and are not known to be background sources. Based on our estimates for the expected number of background sources, we expect that follow-up observations will show that most of these candidates are not dust enshrouded massive stars, with an expectation of only 6 ± 6 surviving candidates. Since we would detect true analogs of η Car for roughly 200 years post-eruption, this implies that the rate of eruptions like η Car is less than the core-collapse supernova rate. It is possible, however, that every M > 40 M _☉ star undergoes such eruptions given our initial results. In Paper II we will characterize the candidates through further analysis and follow-up observations, and there is no barrier to increasing the galaxy sample by an order of magnitude. The primary limitation of the present search is that Spitzer's resolution limits us to the shorter wavelength IRAC bands. With the James Webb Space Telescope, such surveys can be carried out at the far more optimal wavelengths of 10-30 μm, allowing identification of η Car analogs for millennia rather than centuries post-eruption.

The structural isomers ethanol (CH ₃CH ₂OH) and dimethyl ether (CH ₃OCH ₃) were detected in several low-, intermediate-, and high-mass star-forming regions, including Sgr B2, Orion, and W33A, with the relative abundance ratios of ethanol/dimethyl ether varying from about 0.03 to 3.4. Until now, no experimental data regarding the formation mechanisms and branching ratios of these two species in laboratory simulation experiments could be provided. Here, we exploit tunable photoionization reflectron time-of-flight mass spectrometry (PI-ReTOF-MS) to detect and analyze the production of complex organic molecules (COMs) resulting from the exposure of water/methane (H ₂O/CH ₄) ices to energetic electrons. The main goal is to understand the formation mechanisms in star-forming regions of two C ₂H ₆O isomers: ethanol (CH ₃CH ₂OH) and dimethyl ether (CH ₃OCH ₃). The results show that the experimental branching ratios favor the synthesis of ethanol versus dimethyl ether (31 ± 11:1). This finding diverges from the abundances observed toward most star-forming regions, suggesting that production routes on interstellar grains to form dimethyl ether might be missing; alternatively, ethanol can be overproduced in the present simulation experiments, such as via radical–radical recombination pathways involving ethyl and hydroxyl radicals. Finally, the PI-ReTOF-MS data suggest the formation of methylacetylene (C ₃H ₄), ketene (CH ₂CO), propene (C ₃H ₆), vinyl alcohol (CH ₂CHOH), acetaldehyde (CH ₃CHO), and methyl hydroperoxide (CH ₃OOH), in addition to ethane (C ₂H ₆), methanol (CH ₃OH), and CO ₂ detected from infrared spectroscopy. The yield of all the confirmed species is also determined.

Magnetic flux ropes are topological structures consisting of twisted magnetic field lines that globally wrap around an axis. The torus instability model predicts that a magnetic flux rope of major radius R undergoes an eruption when its axis reaches a location where the decay index of the ambient magnetic field B _ex is larger than a critical value. In the current-wire model, the critical value depends on the thickness and time evolution of the current channel. We use magnetohydrodynamic simulations to investigate whether the critical value of the decay index at the onset of the eruption is affected by the magnetic flux rope’s internal current profile and/or by the particular pre-eruptive photospheric dynamics. The evolution of an asymmetric, bipolar active region is driven by applying different classes of photospheric motions. We find that the critical value of the decay index at the onset of the eruption is not significantly affected by either the pre-erupitve photospheric evolution of the active region or the resulting different magnetic flux ropes. As in the case of the current-wire model, we find that there is a “critical range” [1.3–1.5], rather than a “critical value” for the onset of the torus instability. This range is in good agreement with the predictions of the current-wire model, despite the inclusion of line-tying effects and the occurrence of tether-cutting magnetic reconnection.

We report the discovery of deuterium absorption in the very metal-poor ([Fe/H] = –2.88) damped Lyα system at z _abs = 3.06726 toward the QSO SDSS J1358+6522. On the basis of 13 resolved D I absorption lines and the damping wings of the H I Lyα transition, we have obtained a new, precise measure of the primordial abundance of deuterium. Furthermore, to bolster the present statistics of precision D/H measures, we have reanalyzed all of the known deuterium absorption-line systems that satisfy a set of strict criteria. We have adopted a blind analysis strategy (to remove human bias) and developed a software package that is specifically designed for precision D/H abundance measurements. For this reanalyzed sample of systems, we obtain a weighted mean of (D/H) _p = (2.53 ± 0.04) × 10 ^–5, corresponding to a universal baryon density 100 Ω _{b, 0} h ² = 2.202 ± 0.046 for the standard model of big bang nucleosynthesis (BBN). By combining our measure of (D/H) _p with observations of the cosmic microwave background (CMB), we derive the effective number of light fermion species, N _eff = 3.28 ± 0.28. We therefore rule out the existence of an additional (sterile) neutrino (i.e., N _eff = 4.046) at 99.3% confidence (2.7σ), provided that the values of N _eff and of the baryon-to-photon ratio (η ₁₀) did not change between BBN and recombination. We also place a strong bound on the neutrino degeneracy parameter, independent of the ⁴He primordial mass fraction, Y _P: ξ _D = +0.05 ± 0.13 based only on the CMB+(D/H) _p observations. Combining this value of ξ _D with the current best literature measure of Y _P, we find a 2σ upper bound on the neutrino degeneracy parameter, |ξ| ≤ +0.062.

The Parker Solar Probe (PSP) and Solar Orbiter missions are designed to make groundbreaking observations of the Sun and interplanetary space within this decade. We show that a particularly interesting in situ observation of an interplanetary coronal mass ejection (ICME) by PSP may arise during close solar flybys (<0.1 au). During these times, the same magnetic flux rope inside an ICME could be observed in situ by PSP twice, by impacting its frontal part as well as its leg. Investigating the odds of this situation, we forecast the ICME rate in solar cycle 25 based on two models for the sunspot number (SSN): (1) the forecast of an expert panel in 2019 (maximum SSN = 115), and (2) a prediction by McIntosh et al. (2020, maximum SSN = 232). We link the SSN to the observed ICME rates in solar cycles 23 and 24 with the Richardson and Cane list and our own ICME catalog, and calculate that between one and seven ICMEs will be observed by PSP at heliocentric distances <0.1 au until 2025, including 1 σ uncertainties. We then model the potential flux rope signatures of such a double-crossing event with the semiempirical 3DCORE flux rope model, showing a telltale elevation of the radial magnetic field component B _R , and a sign reversal in the component B _N normal to the solar equator compared to field rotation in the first encounter. This holds considerable promise to determine the structure of CMEs close to their origin in the solar corona.

Variations in the total radiative output of the Sun as well as the detailed spectral irradiance are of interest to terrestrial and solar-stellar atmosphere studies. Recent observations provide measurements of spectral irradiance variations at wavelengths in the range 1100-8650 Å with improved accuracy, and correlative studies give procedures for estimating the spectral irradiance changes from solar activity records using indicators such as those derived from Ca II K and Mg II indices. Here we describe our approach to physical modeling of irradiance variations using seven semiempirical models to represent sunspots, plage, network, and quiet atmosphere. This paper gives methods and details, and some preliminary results of our synthesis of the variations of the entire irradiance spectrum. Our calculation uses object-oriented programming techniques that are very efficient and flexible. We compute at high spectral resolution the intensity as a function of wavelength and position on the disk for each of the structure types corresponding to our models. These calculations include three different approximations for the line source function: one suited for the very strong resonance lines where partial redistribution (PRD) is important, another for the most important nonresonance lines, and another approximation for the many narrow lines that are provided in Kurucz's listings. The image analysis and calculations of the irradiance variation as a function of time will be described in a later paper. This work provides an understanding of the sources of variability arising from solar-activity surface structures. We compute the Lyα irradiance to within 3% of the observed values. The difference between our computations and the Neckel & Labs data is 3% or less in the near-IR wavelengths at 8650 Å, and less than 1% in the red at 6080 Å. Near 4100 Å we overestimate the irradiance by 9%-19% because of opacity sources missing in our calculations. We also compute a solar cycle variability of 49% in the Lyα irradiance, which is very close to observed values. At wavelengths between 4100 Å and 1.6 μm, we obtain spectral irradiance variations ranging from -0.06% to 0.46% in the visible—the higher values correspond to the presence of strong lines. The variability in the IR between 1.3 and 2.2 μm is ~-0.15%.

We present Advanced Camera for Surveys, NICMOS, and Keck adaptive-optics-assisted photometry of 20 Type Ia supernovae (SNe Ia) from the Hubble Space Telescope ( HST) Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Of these SNe Ia, 14 pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Of our new SNe Ia, 10 are beyond redshift z = 1, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zero point at the count rates appropriate for very distant SNe Ia. Adding these SNe improves the best combined constraint on dark-energy density, ρ _DE( z), at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a flat ΛCDM universe, we find Ω _Λ = 0.729 ± 0.014 (68% confidence level (CL) including systematic errors). For a flat wCDM model, we measure a constant dark-energy equation-of-state parameter w = –1.013 ^+0.068 _–0.073 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w ₀ and w _a . Further tightening the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z > 1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on board HST. The updated supernova Union2.1 compilation of 580 SNe is available at http://supernova.lbl.gov/Union.

We present an improved determination of the Hubble constant from Hubble Space Telescope (HST) observations of 70 long-period Cepheids in the Large Magellanic Cloud (LMC). These were obtained with the same WFC3 photometric system used to measure extragalactic Cepheids in the hosts of SNe Ia. Gyroscopic control of HST was employed to reduce overheads while collecting a large sample of widely separated Cepheids. The Cepheid period–luminosity relation provides a zero-point-independent link with 0.4% precision between the new 1.2% geometric distance to the LMC from detached eclipsing binaries (DEBs) measured by Pietrzyński et al. and the luminosity of SNe Ia. Measurements and analysis of the LMC Cepheids were completed prior to knowledge of the new DEB LMC distance. Combined with a refined calibration of the count-rate linearity of WFC3-IR with 0.1% precision, these three improved elements together reduce the overall uncertainty in the geometric calibration of the Cepheid distance ladder based on the LMC from 2.5% to 1.3%. Using only the LMC DEBs to calibrate the ladder, we find H ₀ = 74.22 ± 1.82 km s ⁻¹ Mpc ⁻¹ including systematic uncertainties, 3% higher than before for this particular anchor. Combining the LMC DEBs, masers in NGC 4258, and Milky Way parallaxes yields our best estimate: H ₀ = 74.03 ± 1.42 km s ⁻¹ Mpc ⁻¹, including systematics, an uncertainty of 1.91%–15% lower than our best previous result. Removing any one of these anchors changes H ₀ by less than 0.7%. The difference between H ₀ measured locally and the value inferred from Planck CMB and ΛCDM is 6.6 ± 1.5 km s ⁻¹ Mpc ⁻¹ or 4.4 σ ( P = 99.999% for Gaussian errors) in significance, raising the discrepancy beyond a plausible level of chance. We summarize independent tests showing that this discrepancy is not attributable to an error in any one source or measurement, increasing the odds that it results from a cosmological feature beyond ΛCDM.

We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the solar system, exploring the transient optical sky, and mapping the Milky Way. LSST will be a large, wide-field ground-based system designed to obtain repeated images covering the sky visible from Cerro Pachón in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg ² field of view, a 3.2-gigapixel camera, and six filters ( ugrizy) covering the wavelength range 320–1050 nm. The project is in the construction phase and will begin regular survey operations by 2022. About 90% of the observing time will be devoted to a deep-wide-fast survey mode that will uniformly observe a 18,000 deg ² region about 800 times (summed over all six bands) during the anticipated 10 yr of operations and will yield a co-added map to r ∼ 27.5. These data will result in databases including about 32 trillion observations of 20 billion galaxies and a similar number of stars, and they will serve the majority of the primary science programs. The remaining 10% of the observing time will be allocated to special projects such as Very Deep and Very Fast time domain surveys, whose details are currently under discussion. We illustrate how the LSST science drivers led to these choices of system parameters, and we describe the expected data products and their characteristics.

Rings are the most frequently revealed substructure in Atacama Large Millimeter/submillimeter Array (ALMA) dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution (∼0.″12) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high-contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between rings and gaps, and the separations of rings and gaps could all be explained if most gaps are opened by low-mass planets (super-Earths and Neptunes) in the condition of low disk turbulence ( α = 10 ⁻⁴). The gap locations are not well correlated with the expected locations of CO and N ₂ ice lines, so condensation fronts are unlikely to be a universal mechanism to create gaps and rings, though they may play a role in some cases.

We present a new and independent determination of the local value of the Hubble constant based on a calibration of the tip of the red giant branch (TRGB) applied to Type Ia supernovae (SNe Ia). We find a value of H ₀ = 69.8 ± 0.8 (±1.1% stat) ± 1.7 (±2.4% sys) km s ⁻¹ Mpc ⁻¹. The TRGB method is both precise and accurate and is parallel to but independent of the Cepheid distance scale. Our value sits midway in the range defined by the current Hubble tension. It agrees at the 1.2 σ level with that of the Planck Collaboration et al. estimate and at the 1.7 σ level with the Hubble Space Telescope ( HST) SHoES measurement of H ₀ based on the Cepheid distance scale. The TRGB distances have been measured using deep HST Advanced Camera for Surveys imaging of galaxy halos. The zero-point of the TRGB calibration is set with a distance modulus to the Large Magellanic Cloud of 18.477 ± 0.004 (stat) ± 0.020 (sys) mag, based on measurement of 20 late-type detached eclipsing binary stars, combined with an HST parallax calibration of a 3.6 μm Cepheid Leavitt law based on Spitzer observations. We anchor the TRGB distances to galaxies that extend our measurement into the Hubble flow using the recently completed Carnegie Supernova Project I ( CSP-I ) sample containing about 100 well-observed SNe Ia . There are several advantages of halo TRGB distance measurements relative to Cepheid variables; these include low halo reddening, minimal effects of crowding or blending of the photometry, only a shallow (calibrated) sensitivity to metallicity in the I band, and no need for multiple epochs of observations or concerns of different slopes with period. In addition, the host masses of our TRGB host-galaxy sample are higher, on average, than those of the Cepheid sample, better matching the range of host-galaxy masses in the CSP-I distant sample and reducing potential systematic effects in the SNe Ia measurements.

We present a systematic numerical relativity study of the mass ejection and the associated electromagnetic transients and nucleosynthesis from binary neutron star (NS) mergers. We find that a few 10 ⁻³ M _⊙ of material is ejected dynamically during the mergers. The amount and the properties of these outflows depend on binary parameters and on the NS equation of state (EOS). A small fraction of these ejecta, typically ∼10 ⁻⁶ M _⊙, is accelerated by shocks formed shortly after merger to velocities larger than 0.6 c and produces bright radio flares on timescales of weeks, months, or years after merger. Their observation could constrain the strength with which the NSs bounce after merger and, consequently, the EOS of matter at extreme densities. The dynamical ejecta robustly produce second and third r-process peak nuclei with relative isotopic abundances close to solar. The production of light r-process elements is instead sensitive to the binary mass ratio and the neutrino radiation treatment. Accretion disks of up to ∼0.2 M _⊙ are formed after merger, depending on the lifetime of the remnant. In most cases, neutrino- and viscously driven winds from these disks dominate the overall outflow. Finally, we generate synthetic kilonova light curves and find that kilonovae depend on the merger outcome and could be used to constrain the NS EOS.

We present a detailed analysis of two well-localized, highly offset short gamma-ray bursts—GRB 070809 and GRB 090515—investigating the kinematic evolution of their progenitors from compact object formation until merger. Calibrating to observations of their most probable host galaxies, we construct semi-analytic galactic models that account for star formation history and galaxy growth over time. We pair detailed kinematic evolution with compact binary population modeling to infer viable post-supernova velocities and inspiral times. By populating binary tracers according to the star formation history of the host and kinematically evolving their post-supernova trajectories through the time-dependent galactic potential, we find that systems matching the observed offsets of the bursts require post-supernova systemic velocities of hundreds of kilometers per second. Marginalizing over uncertainties in the stellar mass–halo mass relation, we find that the second-born neutron star in the GRB 070809 and GRB 090515 progenitor systems received a natal kick of at the and credible levels, respectively. Applying our analysis to the full catalog of localized short gamma-ray bursts will provide unique constraints on their progenitors and help unravel the selection effects inherent to observing transients that are highly offset with respect to their hosts.

In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. Particularly crucial is the question of what sets the peak of the distribution. To investigate this issue, we perform high-resolution numerical simulations with radiative feedback exploring, in particular, the role of the stellar and accretion luminosities. We also perform simulations with a simple effective equation of state (EOS), and we investigate 1000 solar-mass clumps having, respectively, 0.1 and 0.4 pc of initial radii. We found that most runs, both with radiative transfer or an EOS, present similar mass spectra with a peak broadly located around 0.3–0.5 M _⊙ and a power-law-like mass distribution at higher masses. However, when accretion luminosity is accounted for, the resulting mass spectrum of the most compact clump tends to be moderately top-heavy. The effect remains limited for the less compact one, which overall remains colder. Our results support the idea that rather than the radiative stellar feedback, this is the transition from the isothermal to the adiabatic regime, which occurs at a gas density of about 10 ¹⁰ cm ⁻³, that is responsible for setting the peak of the IMF. This stems from (i) the fact that extremely compact clumps for which the accretion luminosity has a significant influence are very rare and (ii) the luminosity problem, which indicates that the effective accretion luminosity is likely weaker than expected.

From soft X-ray emission, the solar flare temperatures are from several MK to dozens of times MK, which are higher than the preflare coronal temperatures. A combination of several heating mechanisms may contribute to the heating problem in solar flare loops. In this paper, we propose an important mechanism of solar flaring loops heating, in which the excited electron acoustic wave (EAW) by flare-accelerated fast electron beams can lead to electron heating via collisionless Landau damping effect produced by wave–particle resonant interaction. Taking account of the return-current effect of fast electron beams, by use of numerical and analytic solutions, the plasma wave instability driven by fast electron beams is investigated in typical solar flare loop plasma parameters. The results show that the EAW is the strongest unstable wave mode rather than other wave modes. The dissipation of EAW via collisionless Landau damping and its application to solar flaring loops heating are discussed in detail.

Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment ( AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds, and supernovae in FLASH. We here describe the algorithms used for each of these processes. We denote this integrated package Torch. We then use this novel numerical method to simulate the formation and early evolution of several examples of open clusters of ∼1000 stars formed from clouds with a mass range of 10 ³ M _⊙ to 10 ⁵ M _⊙. Analyzing the effects of stellar feedback on the gas and stars of the natal clusters, we find that in these examples, the stellar clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H  ii region shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations.

The gravitational pull of a large number of asteroids perturbs a pulsar’s motion to a degree that is detectable through precision timing of millisecond pulsars. The result is a low-frequency, correlated noise process, similar in form to the red timing noise known to affect canonical pulsars, or to the signal expected from a stochastic gravitational-wave background. Motivated by the observed fact that many millisecond pulsars are in binary systems, we describe the ways in which the presence of a binary companion to the pulsar would affect the signal produced by an asteroid belt. The primary effect of the companion is to destabilize the shortest-period orbits, cutting off the high-frequency component of the signal from the asteroid belt. We also discuss the implications of asteroid belts for gravitational-wave search efforts. Compared to the signal from a stochastic gravitational-wave background, asteroid-belt noise has a similar frequency and amplitude, and is similarly independent of radio frequency, but is not correlated between different pulsars, which should allow the two kinds of signal to be distinguished.