An official website of the United States government
Abstract We examine the geometry of the post–asymptotic giant branch (AGB) star binary AC Her and its circumbinary disk. We show that the observations describe a binary orbit that is perpendicular to the disk with an angular momentum vector that is within 9° of the binary eccentricity vector, meaning that the disk is close to a stable polar alignment. The most likely explanation for the very large inner radius of the dust is a planet within the circumbinary disk. This is therefore both the first reported detection of a polar circumbinary disk around a post-AGB binary and the first evidence of a polar circumbinary planet. We consider the dynamical constraints on the circumbinary disk size and mass. The polar circumbinary disk feeds circumstellar disks with gas on orbits that are highly inclined with respect to the binary orbit plane. The resulting circumstellar disk inclination could be anywhere from coplanar to polar depending upon the competition between the mass accretion and binary torques.
more »
« less
The Santa Barbara Binary−disk Code Comparison
Abstract We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary−disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary−disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a2ΩB/ν) binary orbits, whereνis the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.
more »
« less
- PAR ID:
- 10527826
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 970
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 156
- Size(s):
- Article No. 156
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with a minidisk around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by am= 1 overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the minidisk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the minidisks. We find that minidisks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases we find that most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes, and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.more » « less
-
Abstract Linear analysis of gas flows around orbiting binaries suggests that a centrifugal barrier ought to clear a low-density cavity around the binary and inhibit mass transfer onto it. Modern hydrodynamics simulations have confirmed the low-density cavity, but show that any mass flowing from large scales into the circumbinary disk is eventually transferred onto the binary components. Even though many numerical studies confirm this picture, it is still not understood precisely how gas parcels overcome the centrifugal barrier and ultimately accrete. We present a detailed analysis of the binary accretion process, using an accurate prescription for evolving grid-based hydrodynamics with Lagrangian tracer particles that track the trajectories of individual gas parcels. We find that binary accretion can be described in four phases: (1) gas is viscously transported through the circumbinary disk up to the centrifugal barrier at the cavity wall, (2) the cavity wall is tidally distorted into accretion streams consisting of near-ballistic gas parcels on eccentric orbits, (3) the portion of each stream moving inwards of an accretion horizon radius r ¯ ≃ a —the radius beyond which no material is returned to the cavity wall—becomes bound to a minidisk orbiting an individual binary component, and (4) the minidisk gas accretes onto the binary component through the combined effect of viscous and tidal stresses.more » « less
-
Abstract The detection of GW170817/AT2017gfo inaugurated an era of multimessenger astrophysics, in which gravitational-wave and multiwavelength photon observations complement one another to provide unique insight into astrophysical systems. A broad theoretical consensus exists, in which the photon phenomenology of neutron star mergers largely rests upon the evolution of the small amount of matter left on bound orbits around the black hole or massive neutron star remaining after the merger. Because this accretion disk is far from inflow equilibrium, its subsequent evolution depends very strongly on its initial state, yet very little is known about how this state is determined. Using both snapshot and tracer particle data from a numerical relativity/MHD simulation of an equal-mass neutron star merger that collapses to a black hole, we show how gravitational forces arising in a nonaxisymmetric, dynamical spacetime supplement hydrodynamical effects in shaping the initial structure of the bound debris disk. The work done by hydrodynamical forces is ∼10 times greater than that due to time-dependent gravity. Although gravitational torques prior to remnant relaxation are an order of magnitude larger than hydrodynamical torques, their intrinsic sign symmetry leads to strong cancellation; as a result, hydrodynamical and gravitational torques have a comparable effect. We also show that the debris disk’s initial specific angular momentum distribution is sharply peaked at roughly the specific angular momentum of the merged neutron star’s outer layers, a fewrgc, and identify the regulating mechanism.more » « less
-
Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process (r-process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We presentνh-amr, a novel neutrino-transport general relativistic magnetohydrodynamic (νGRMHD) code, which we use to perform the first 3D globalνGRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, , do not produce significant neutron-rich (Ye < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, (e.g., for more centrally concentrated progenitors), can unbindMej ≲ M⊙of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflowYe; however, the finalr-process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines becomer-process factories.more » « less
