No Arabic abstract
We present the first simulations of the tidal disruption of stars with realistic structures and compositions by massive black holes (BHs). We build stars in the stellar evolution code MESA and simulate their disruption in the 3D adaptive-mesh hydrodynamics code FLASH, using an extended Helmholtz equation of state and tracking 49 elements. We study the disruption of a 1$M_odot$ star and 3$M_odot$ star at zero-age main sequence (ZAMS), middle-age, and terminal-age main sequence (TAMS). The maximum BH mass for tidal disruption increases by a factor of $sim$2 from stellar radius changes due to MS evolution; this is equivalent to varying BH spin from 0 to 0.75. The shape of the mass fallback rate curves is different from the results for polytropes of Guillochon & Ramirez-Ruiz (2013). The peak timescale $t_{rm peak}$ increases with stellar age, while the peak fallback rate $dot M_{rm peak}$ decreases with age, and these effects diminish with increasing impact parameter $beta$. For a $beta=1$ disruption of a 1$M_odot$ star by a $10^6 M_odot$ BH, from ZAMS to TAMS, $t_{rm peak}$ increases from 30 to 54 days, while $dot M_{rm peak}$ decreases from 0.66 to 0.14 $M_odot$/yr. Compositional anomalies in nitrogen, helium, and carbon can occur before the peak timescale for disruptions of MS stars, which is in contrast to predictions from the frozen-in model. More massive stars can show stronger anomalies at earlier times, meaning that compositional constraints can be key in determining the mass of the disrupted star. The abundance anomalies predicted by these simulations provide a natural explanation for the spectral features and varying line strengths observed in tidal disruption events.
Binary stars that are on close orbits around massive black holes (MBH) such as Sgr A* in the center of the Milky Way are liable to undergo tidal disruption and eject a hypervelocity star. We study the interaction between such a MBH and circular binaries for general binary orientations and penetration depths (i.e. binaries penetrate into the tidal radius around the BH). We show that for very deep penetrators, all binaries are disrupted when the binary rotation axis is roughly oriented toward the BH or it is in the opposite direction. The surviving chance becomes significant when the angle between the binary rotation axis and the BH direction is between pi /4 and 3 pi /4. The surviving chance is as high as $sim$ 20$%$ when the binary rotation axis is perpendicular to the BH direction. The angular dependence is opposite for very shallow penetrators where coplanar prograde orbits have the lowest surviving chance (or equivalently most vulnerable). We provide numerical fits to the disruption probability and energy gain at the the BH encounter as a function of the penetration depth. The latter can be simply rescaled in terms of binary masses, their initial separation and the binary-to-BH mass ratio to evaluate the ejection velocity of a binary members in various systems. We also investigate the disruption of coplanar, eccentric binaries by a MBH. It is shown that for highly eccentric binaries retrograde orbits have a significantly increased disruption probability and ejection velocities compared to the circular binaries.
When a star approaches a black hole closely, it may be pulled apart by gravitational forces in a tidal disruption event (TDE). The flares produced by TDEs are unique tracers of otherwise quiescent supermassive black holes (SMBHs) located at the centre of most galaxies. In particular, the appearance of such flares and the subsequent decay of the light curve are both sensitive to whether the star is partially or totally destroyed by the tidal field. However, the physics of the disruption and the fall-back of the debris are still poorly understood. We are here modelling the hydrodynamical evolution of realistic stars as they approach a SMBH on parabolic orbits, using for the first time the moving-mesh code AREPO, which is particularly well adapted to the problem through its combination of quasi-Lagrangian behaviour, low advection errors, and high accuracy typical of mesh-based techniques. We examine a suite of simulations with different impact parameters, allowing us to determine the critical distance at which the star is totally disrupted, the energy distribution and the fallback rate of the debris, as well as the hydrodynamical evolution of the stellar remnant in the case of a partial disruption. Interestingly, we find that the internal evolution of the remnants core is strongly influenced by persistent vortices excited in the tidal interaction. These should be sites of strong magnetic field amplification, and the associated mixing may profoundly alter the subsequent evolution of the tidally pruned star.
Aims: A strong, hard X-ray flare was discovered (IGR J12580+0134) by INTEGRAL in 2011, and is associated to NGC 4845, a Seyfert 2 galaxy never detected at high-energy previously. To understand what happened we observed this event in the X-ray band on several occasions. Methods: Follow-up observations with XMM-Newton, Swift, and MAXI are presented together with the INTEGRAL data. Long and short term variability are analysed and the event wide band spectral shape modelled. Results: The spectrum of the source can be described with an absorbed (N_H ~ 7x10^22 cm^{-2}) power law (Gamma simeq 2.2), characteristic of an accreting source, plus a soft X-ray excess, likely to be of diffuse nature. The hard X-ray flux increased to maximum in a few weeks and decreased over a year, with the evolution expected for a tidal disruption event. The fast variations observed near the flare maximum allowed us to estimate the mass of the central black hole in NGC 4845 as ~ 3x10^5 Msun. The observed flare corresponds to the disruption of about 10% of an object with a mass of 14-30 Jupiter. The hard X-ray emission should come from a corona forming around the accretion flow close to the black hole. This is the first tidal event where such a corona has been observed.
The question of how supermassive black holes (SMBHs) grow over cosmic time is a major puzzle in high-energy astrophysics. One promising approach to this problem is via the study of tidal disruption flares (TDFs). These are transient events resulting from the disruption of stars by quiescent supermassive black holes at centers of galaxies. A meter-class X-ray observatory with a time resolution $sim$ a millisecond and a spectral resolution of a few eV at KeV energies would be revolutionary as it will facilitate high signal to noise spectral-timing studies of several cosmological TDFs. It would open a new era of astrophysics where SMBHs in TDFs at cosmic distances can be studied in similar detail as current studies of much nearer, stellar-mass black hole binaries. Using Athena X-ray observatory as an example, we highlight two specific aspects of spectral-timing analysis of TDFs. (1) Detection of X-ray quasi-periodic oscillations (QPOs) over a redshift range and using these signal frequencies to constrain the spin evolution of SMBHs, and (2) Time-resolved spectroscopy of outflows/winds to probe super-Eddington accretion. SMBH spin distributions at various redshifts will directly allow us to constrain their primary mode of growth as higher spins are predicted due to spin-up for prolonged accretion-mode growth, while lower spins are expected for growth via mergers due to angular momentum being deposited from random directions. A meter-class X-ray telescope will also be able to characterize relativistic TDFs, viz., SwJ1644+57-like events, out to a redshift greater than 8, i.e., it would facilitate detailed spectral-timing studies of TDFs by the youngest SMBHs in the Universe.
We present results from general relativistic calculations of the tidal disruption of white dwarf stars from near encounters with intermediate mass black holes. We follow the evolution of 0.2 and $0.6 M_odot$ stars on parabolic trajectories that approach $10^3$ - $10^4 M_odot$ black holes as close as a few Schwarzschild radii at periapsis, paying particular attention to the effect tidal disruption has on thermonuclear reactions and the synthesis of intermediate to heavy ion elements. These encounters create diverse thermonuclear environments characteristic of Type I supernovae and capable of producing both intermediate and heavy mass elements in arbitrary ratios, depending on the strength (or proximity) of the interaction. Nuclear ignition is triggered in all of our calculations, even at weak tidal strengths $beta sim 2.6$ and large periapsis radius $R_P sim 28$ Schwarzschild radii. A strong inverse correlation exists between the mass ratio of calcium to iron group elements and tidal strength, with $beta lesssim 5$ producing predominately calcium-rich debris. At these moderate to weak interactions, nucleosynthesis is not especially efficient, limiting the total mass and outflows of calcium group elements to $< 15$% of available nuclear fuel. Iron group elements however continue to be produced in greater quantity and ratio with increasing tidal strength, peaking at $sim 60$% mass conversion efficiency in our closest encounter cases. These events generate short bursts of gravitational waves with characteristic frequencies 0.1-0.7 Hz and strain amplitudes $0.5times10^{-22}$ - $3.5times10^{-22}$ at 10 Mpc source distance.