The Eccentric Nature of Eccentric Tidal Disruption Events

Upon entering the tidal sphere of a supermassive black hole, a star is ripped apart by tides and transformed into a stream of debris. The ultimate fate of that debris, and the properties of the bright flare that is produced and observed, depends on a number of parameters, including the energy of the center of mass of the original star. Here we present the results of a set of smoothed particle hydrodynamics simulations in which a 1M ⊙, γ = 5/3 polytrope is disrupted by a 106 M ⊙ supermassive black hole. Each simulation has a pericenter distance of r p = r t (i.e., β ≡ r t/r p = 1 with r t the tidal radius), and we vary the eccentricity e of the stellar orbit from e = 0.8 up to e = 1.20 and study the nature of the fallback of debris onto the black hole and the long-term fate of the unbound material. For simulations with eccentricities e ≲ 0.98, the fallback curve has a distinct, three-peak structure that is induced by self-gravity. For simulations with eccentricities e ≳ 1.06, the core of the disrupted star reforms following its initial disruption. Our results have implications for, e.g., tidal disruption events produced by supermassive black hole binaries.

Stars Crushed by Black Holes. I. On the Energy Distribution of Stellar Debris in Tidal Disruption Events

The distribution of orbital energies imparted into stellar debris following the close encounter of a star with a supermassive black hole is the principal factor in determining the rate of return of debris to the black hole, and thus in determining the properties of the resulting lightcurves from such events. We present simulations of tidal disruption events for a range of β ≡ r t/r p where r p is the pericenter distance and r t the tidal radius. We perform these simulations at different spatial resolutions to determine the numerical convergence of our models. We compare simulations in which the heating due to shocks is included or excluded from the dynamics. For β ≲ 8, the simulation results are well-converged at sufficiently moderate-to-high spatial resolution, while for β ≳ 8, the breadth of the energy distribution can be grossly exaggerated by insufficient spatial resolution. We find that shock heating plays a non-negligible role only for β ≳ 4, and that typically the effect of shock heating is mild. We show that self-gravity can modify the energy distribution over time after the debris has receded to large distances for all β. Primarily, our results show that across a range of impact parameters, while the shape of the energy distribution varies with β, the width of the energy spread imparted to the bulk of the debris is closely matched to the canonical spread, Δ E = GM • R ⋆ / r t 2 , for the range of β we have simulated.

The JCMT BISTRO Survey: Evidence for Pinched Magnetic Fields in Quiescent Filaments of NGC 1333

We investigate the internal 3D magnetic structure of dense interstellar filaments within NGC 1333 using polarization data at 850 μm from the B-fields In STar-forming Region Observations survey at the James Clerk Maxwell Telescope. Theoretical models predict that the magnetic field lines in a filament will tend to be dragged radially inward (i.e., pinched) toward the central axis due to the filament’s self-gravity. We study the cross-sectional profiles of the total intensity (I) and polarized intensity (PI) of dust emission in four segments of filaments unaffected by local star formation that are expected to retain a pristine magnetic field structure. We find that the filaments’ FWHMs in PI are not the same as those in I, with two segments being appreciably narrower in PI (FWHM ratio ≃0.7–0.8) and one segment being wider (FWHM ratio ≃1.3). The filament profiles of the polarization fraction (P) do not show a minimum at the spine of the filament, which is not in line with an anticorrelation between P and I normally seen in molecular clouds and protostellar cores. Dust grain alignment variation with density cannot reproduce the observed P distribution. We demonstrate numerically that the I and PI cross-sectional profiles of filaments in magnetohydrostatic equilibrium will have differing relative widths depending on the viewing angle. The observed variations of FWHM ratios in NGC 1333 are therefore consistent with models of pinched magnetic field structures inside filaments, especially if they are magnetically near-critical or supercritical.

Partial, Zombie, and Full Tidal Disruption of Stars by Supermassive Black Holes

We present long-duration numerical simulations of the tidal disruption of stars modeled with accurate stellar structures and spanning a range of pericenter distances, corresponding to cases where the stars are partially and completely disrupted. We substantiate the prediction that the late-time power-law index of the fallback rate n ∞ ≃ −5/3 for full disruptions, while for partial disruptions—in which the central part of the star survives the tidal encounter intact—we show that n ∞ ≃ −9/4. For the subset of simulations where the pericenter distance is close to that which delineates full from partial disruption, we find that a stellar core can reform after the star has been completely destroyed; for these events the energy of the zombie core is slightly positive, which results in late-time evolution from n ≃ −9/4 to n ≃ −5/3. We find that self-gravity can generate an n(t) that deviates from n ∞ by a small but significant amount for several years post-disruption. In one specific case with the stellar pericenter near the critical value, we find that self-gravity also drives the recollapse of the central regions of the debris stream into a collection of several cores while the rest of the stream remains relatively smooth. We also show that it is possible for the surviving stellar core in a partial disruption to acquire a circumstellar disk that is shed from the rapidly rotating core. Finally, we provide a novel analytical fitting function for the fallback rates that may also be useful in a range of contexts beyond tidal disruption events.

Dynamically Driven Inflow onto the Galactic Center and its Effect upon Molecular Clouds

The Galactic bar plays a critical role in the evolution of the Milky Way’s Central Molecular Zone (CMZ), driving gas toward the Galactic Center via gas flows known as dust lanes. To explore the interaction between the CMZ and the dust lanes, we run hydrodynamic simulations in arepo, modeling the potential of the Milky Way’s bar in the absence of gas self-gravity and star formation physics, and we study the flows of mass using Monte Carlo tracer particles. We estimate the efficiency of the inflow via the dust lanes, finding that only about a third (30% ± 12%) of the dust lanes’ mass initially accretes onto the CMZ, while the rest overshoots and accretes later. Given observational estimates of the amount of gas within the Milky Way’s dust lanes, this suggests that the true total inflow rate onto the CMZ is 0.8 ± 0.6 M ⊙ yr−1. Clouds in this simulated CMZ have sudden peaks in their average density near the apocenter, where they undergo violent collisions with inflowing material. While these clouds tend to counter-rotate due to shear, co-rotating clouds occasionally occur due to the injection of momentum from collisions with inflowing material (∼52% are strongly counter-rotating, and ∼7% are strongly co-rotating of the 44 cloud sample). We investigate the formation and evolution of these clouds, finding that they are fed by many discrete inflow events, providing a consistent source of gas to CMZ clouds even as they collapse and form stars.

Resolving the Fastest Ejecta from Binary Neutron Star Mergers: Implications for Electromagnetic Counterparts

We examine the effect of spatial resolution on initial mass ejection in grid-based hydrodynamic simulations of binary neutron star mergers. The subset of the dynamical ejecta with velocities greater than ∼0.6c can generate an ultraviolet precursor to the kilonova on approximately hour timescales and contribute to a years long nonthermal afterglow. Previous work has found differing amounts of this fast ejecta, by one to two orders of magnitude, when using particle-based or grid-based hydrodynamic methods. Here, we carry out a numerical experiment that models the merger as an axisymmetric collision in a corotating frame, accounting for Newtonian self-gravity, inertial forces, and gravitational wave losses. The lower computational cost allows us to reach spatial resolutions as high as 4 m, or ∼3 × 10−4 of the stellar radius. We find that fast ejecta production converges to within 10% for a cell size of 20 m. This suggests that fast ejecta quantities found in existing grid-based merger simulations are unlikely to increase to the level needed to match particle-based results upon further resolution increases. The resulting neutron-powered precursors are in principle detectable out to distances ≲200 Mpc with upcoming facilities.We also find that head-on collisions at the freefall speed, relevant for eccentric mergers, yield fast and slow ejecta quantities of order 10−2 M ⊙, with a kilonova signature distinct from that of quasi-circular mergers.

D-dimensional self-gravitating lattice gas in general relativity

Using a lattice equation of state combined with the D-dimensional Tolman–Oppenheimer–Volkoff equation and the Friedmann equations, we investigate the possibility of the formation of compact objects as well as the time evolution of the scale factor and the density profile of a self-gravitating material cluster. The numerical results show that in a ([Formula: see text])-dimensional space–time, the mass is independent of the central pressure. Hence, the formation of only compact objects with a finite constant mass similar to the white dwarf is possible. However, in a ([Formula: see text])-dimensional space–time, self-gravity leads to the formation of compact objects with a large gap of mass and the corresponding phase diagram has the same structure as the one for Neutron Star. The results also show that beyond certain critical central pressure, the star is unstable against gravitational collapse, and it may end in a black hole. Analysis of space–times of higher dimensions shows that gravity has the stronger effect in [Formula: see text] dimensions. Numerical solutions of the Friedmann equations show that the effect of the curvature of space–time increases with the increasing temperature, but decreases with the increasing dimensionality beyond [Formula: see text].

Asteroids are young: global-scale reshaping and resurfacing by small-scale impacts

The fraction of the asteroid population that survived since the Solar System formation has experienced numerous collisional, dynamical and thermal events, which have shaped their structures and orbital properties. Small asteroids are often considered to be rubble-pile objects, aggregates held together only by self-gravity or small cohesive forces (1; 2). The artificial impact experiment of JAXA’s Hayabusa2 mission on the surface of asteroid Ryugu (3) created a surprisingly large crater (≈14 m). This unexpected result suggests that at least the near-surface of the asteroid is controlled to a large extent by its rather weak gravity rather than strength. Due to the inability to re-create these impact conditions in laboratory experiments, this observed regime of low-gravity, low-strength cratering remained largely unexplored so far. In addition, the very large times scales involved in the crater growth made it impossible to numerically simulate these impact processes up to now. Here we use a novel approach to model the entire cratering process resulting from impacts on small, weak asteroids, which uses shock physics code calculations directly. We found that small impacts can significantly deform weak asteroids, causing global resurfacing at the same time. We also show that even very low asteroid cohesions can drastically influence the outcome of an impact and that the collisional life-time of the overall asteroid shapes is significantly lower than the traditionally used life-time based on catastrophic disruption events. Consequently, we predict that NASA’s Double Asteroid Redirection Test (DART) impact on Dimorhpos (4; 5) will not lead to a cratering event, as originally anticipated (i.e., 6; 7). Rather, the impact is going to change the global morphology of the asteroid, if the surface cohesion is less than ≈ 10 Pa. Our results, together with the future observations by the ESA’s Hera mission (8) will provide constraints regarding the evolution of the shapes and structures of small asteroids by sub-catastrophic impacts.

Fragmentation of ring galaxies and transformation to clumpy galaxies

We study the fragmentation of collisional ring galaxies (CRGs) using a linear perturbation analysis that computes the physical conditions of gravitational instability, as determined by the balance of self-gravity of the ring against pressure and Coriolis forces. We adopt our formalism to simulations of CRGs and show that the analysis can accurately characterise the stability and onset of fragmentation, although the linear theory appears to under-predict the number of fragments of an unstable CRG by a factor of 2. In addition, since the orthodox ‘density-wave’ model is inapplicable to such self-gravitating rings, we devise a simple approach that describes the rings propagating as material waves. We find that the toy model can predict whether the simulated CRGs fragment or not using information from their pre-collision states. We also apply our instability analysis to a CRG discovered at a high redshift, z = 2.19. We find that a quite high velocity dispersion is required for the stability of the ring, and therefore the CRG should be unstable to ring fragmentation. CRGs are rarely observed at high redshifts, and this may be because CRGs are usually too faint. Since the fragmentation can induce active star formation and make the ring bright enough to observe, the instability could explain this rarity. An unstable CRG fragments into massive clumps retaining the initial disc rotation, and thus it would evolve into a clumpy galaxy with a low surface density in an inter-clump region.