Pericenter Passage Research Articles

Magnetic field evolution in tidal disruption events

When a star gets tidally disrupted by a supermassive black hole, its magnetic field is expected to pervade its debris. In this paper, we study this process via smoothed particle magnetohydrodynamical simulations of the disruption and early debris evolution including the stellar magnetic field. As the gas stretches into a stream, we show that the magnetic field evolution is strongly dependent on its orientation with respect to the stretching direction. In particular, an alignment of the field lines with the direction of stretching induces an increase of the magnetic energy. For disruptions happening well within the tidal radius, the star compression causes the magnetic field strength to sharply increase by an order of magnitude at the time of pericentre passage. If the disruption is partial, we find evidence for a dynamo process occurring inside the surviving core due to the formation of vortices. This causes an amplification of the magnetic field strength by a factor of $\sim 10$. However, this value represents a lower limit since it increases with numerical resolution. For an initial field strength of 1 G, the magnetic field never becomes dynamically important. Instead, the disruption of a star with a strong 1 MG magnetic field produces a debris stream within which magnetic pressure becomes similar to gas pressure a few tens of hours after disruption. If the remnant of one or multiple partial disruptions is eventually fully disrupted, its magnetic field could be large enough to magnetically power the relativistic jet detected from Swift J1644+57. Magnetized streams could also be significantly thickened by magnetic pressure when it overcomes the confining effect of self-gravity.

A survey of dual active galactic nuclei in simulations of galaxy mergers: frequency and properties

We investigate the simultaneous triggering of active galactic nuclei (AGN) in merging galaxies, using a large suite of high-resolution hydrodynamical simulations. We compute dual-AGN observability time-scales using bolometric, X-ray, and Eddington-ratio thresholds, confirming that dual activity from supermassive black holes (BHs) is generally higher at late pericentric passages, before a merger remnant has formed, especially at high luminosities. For typical minor and major mergers, dual activity lasts ~20-70 and ~100-160 Myr, respectively. We also explore the effects of X-ray obscuration from gas, finding that the dual-AGN time decreases at most by a factor of ~2, and of contamination from star formation. Using projected separations and velocity differences rather than three-dimensional quantities can decrease the dual-AGN time-scales by up to ~4, and we apply filters which mimic current observational-resolution limitations. In agreement with observations, we find that, for a sample of major and minor mergers hosting at least one AGN, the fraction harbouring dual AGN is ~20-30 and ~1-10 per cent, respectively. We quantify the effects of merger mass ratio (0.1 to 1), geometry (coplanar, prograde, retrograde, and inclined), disc gas fraction, and BH properties, finding that the mass ratio is the most important factor, with the difference between minor and major mergers varying between factors of a few to orders of magnitude, depending on the luminosity and filter used. We also find that a shallow imaging survey will require very high angular resolution, whereas a deep imaging survey will be less resolution-dependent.

Observations of Sagittarius A* during the pericenter passage of the G2 object with MAGIC

Context. We present the results of a multi-year monitoring campaign of the Galactic Center (GC) with the MAGIC telescopes. These observations were primarily motivated by reports that a putative gas cloud (G2) would be passing in close proximity to the super-massive black hole (SMBH), associated with Sagittarius A*, located at the center of our galaxy. This event was expected to give astronomers a unique chance to study the effect of in-falling matter on the broad-band emission of a SMBH. Aims. We search for potential flaring emission of very-high-energy (VHE; $\geq$100 GeV) gamma rays from the direction of the SMBH at the GC due to the passage of the G2 object. Using these data we also study the morphology of this complex region. Methods. We observed the GC region with the MAGIC Imaging Atmospheric Cherenkov Telescopes during the period 2012-2015, collecting 67 hours of good-quality data. In addition to a search for variability in the flux and spectral shape of the GC gamma-ray source, we use a point-source subtraction technique to remove the known gamma-ray emitters located around the GC in order to reveal the TeV morphology of the extended emission inside that region. Results. No effect of the G2 object on the VHE gamma-ray emission from the GC was detected during the 4 year observation campaign. We confirm previous measurements of the VHE spectrum of Sagittarius A*, and do not detect any significant variability of the emission from the source. Furthermore, the known VHE gamma-ray emitter at the location of the supernova remnant G0.9+0.1 was detected, as well as the recently discovered VHE source close to the GG radio Arc.

A possible time-delayed brightening of the Sgr A* accretion flow after the pericenter passage of the G2 cloud

Abstract A possibility of time-delayed radio brightenings of Sgr A* triggered by the pericenter passage of the G2 cloud is studied by carrying out global three-dimensional magnetohydrodynamic simulations, taking into account the radiative cooling of the tidal debris of the G2 cloud. Magnetic fields in the accretion flow are strongly perturbed and reorganized after the passage of G2. We have found that the magnetic energy in the accretion flow increased by a factor of 3–4 in 5–10 yr after the passage of G2 through a dynamo mechanism driven by the magneto-rotational instability. Since this B-field amplification enhances the synchrotron emission from the disk and the outflow, the radio and the infrared luminosity of Sgr A* are expected to increase some time, around 2020. The time delay of the radio brightening enables us to determine the rotation axis of the preexisting disk.

Weak lensing and spectroscopic analysis of the nearby dissociative merging galaxy cluster Abell 3376

The galaxy cluster Abell~3376 is a nearby (z=0.046) dissociative merging cluster surrounded by two prominent radio relics and showing an X-ray comet-like morphology. The merger system is comprised of the subclusters A3376W & A3376E. Based on new deep multi-wavelength large-field images and published redshifts, we bring new insights about the history of this merger. Despite the difficulty of applying the weak lensing technique at such low redshift, we successfully recovered the mass distribution in the cluster field. Moreover, with the application of a two-body model, we have addressed the dynamics of these merging system. We have found the individual masses of M_{200}^{W}=3.0_{-1.7}^{+1.3}x10^{14} M_{\odot} and M_{200}^{E}=0.9_{-0.8}^{+0.5}x10^{14} M_{\odot}. The cometary shaped X-ray distribution shows only one peak spatially coincident with both Eastern BCG and the A3376E mass peak whereas the gas content of A3376W seems to be stripped out. Our data allowed us to confirm the existence of a third subcluster located at the North, 1147+-62 kpc apart from the neighbour subcluster A3376E and having a mass M_{200}^{N}=1.4_{-1.0}^{+0.7}x10^{14} M_{\odot}. From our dynamical analysis, we found the merging is taking place very close to the plane of the sky, with the merger axis just 10 deg +-11 deg from it. The application of a two-body analysis code showed that the merging cluster is seen 0.9_{-0.3}^{+0.2} Gyr after the pericentric passage and it is currently going to the point of maximum separation between the subclusters.

Tidal dissipation and evolution of white dwarfs around massive black holes: an eccentric path to tidal disruption

Awhite dwarf (WD)captured into a high-eccentricity orbit around amassive black hole (MBH) may undergo many pericentre passages before tidal disruption. During these passages, the tidal potential of the MBH excites internal oscillations or waves in the WD, and the dissipation of these oscillations can significantly influence the physical properties of the WD prior to its disruption. We calculate the amplitude of the tidally excited gravity (buoyancy) waves in the WD as a function of the pericentre distance and eccentricity for realistic WD models, under the assumption that these outgoing gravity waves are efficiently dissipated in the outer layers of the WD by non-linear effects or radiative damping. We obtain fitting formulae for the tidal energy and angular momentum transfer rates as well as the tidal heating rate.We find that these dynamical tides are much weaker than gravitational radiation in driving the orbital decay of theWD–MBH binary, and they are also inefficient in changing theWD spin during the orbital evolution. Incorporating our computed tidal dissipation rate into a MESA-based WD evolution code, we find that tidal heating can lead to appreciable brightening of theWD and may induce runaway fusion in the hydrogen envelope well before the WD undergoes tidal disruption.

Lidov–Kozai Cycles with Gravitational Radiation: Merging Black Holes in Isolated Triple Systems

Abstract We show that a black-hole binary with an external companion can undergo Lidov–Kozai cycles that cause a close pericenter passage, leading to a rapid merger due to gravitational-wave emission. This scenario occurs most often for systems in which the companion has a mass comparable to the reduced mass of the binary and the companion orbit has a semimajor axis within a factor of ∼10 of the binary semimajor axis. Using a simple population-synthesis model and three-body simulations, we estimate the rate of mergers in triple black-hole systems in the field to be about six per Gpc3 per year in the absence of natal kicks during black-hole formation. This value is within the low end of the 90% credible interval for the total black hole–black hole merger rate inferred from the current LIGO results. There are many uncertainties in these calculations, the largest of which is the unknown distribution of natal kicks. Even modest natal kicks of 40 km s−1 will reduce the merger rate by a factor of 40. A few percent of these systems will have eccentricity greater than 0.999 when they first enter the frequency band detectable by aLIGO (above 10 Hz).

EXPLODING SATELLITES—THE TIDAL DEBRIS OF THE ULTRA-FAINT DWARF GALAXY HERCULES

ABSTRACT The ultra-faint satellite galaxy Hercules has a strongly elongated and irregular morphology with detections of tidal features up to 1.3 deg (3 kpc) from its center. This suggests that Hercules may be dissolving under the Milky Way’s gravitational influence, and hence could be a tidal stream in formation rather than a bound, dark-matter-dominated satellite. Using Bayesian inference in combination with N-body simulations, we show that Hercules has to be on a very eccentric orbit ( ) within the Milky Way in this scenario. On such an orbit, Hercules “explodes” as a consequence of the last tidal shock at pericenter 0.5 Gyr ago. It is currently decelerating toward the apocenter of its orbit with a velocity of V = 157 km s−1—of which 99% is directed radially outwards. Due to differential orbital plane precession caused by the non-spherical nature of the Galactic potential, its debris fans out nearly perpendicular to its orbit. This explains why Hercules has an elongated shape without showing a distance gradient along its main body: it is in fact a stream that is significantly broader than it is long. In other words, it is moving perpendicular to its apparent major axis. In this scenario, there is a spike in the radial velocity profile created by the dominant debris component that formed through the last pericenter passage. This is similar to kinematic substructure that is observed in the real Hercules. Modeling a satellite on such a highly eccentric orbit is strongly dependent on the form of the Galactic potential. We therefore propose that detailed kinematic investigation of Hercules and other exploding satellite candidates can yield strong constraints on the potential of the Milky Way.

TIDAL ORIGIN OF SPIRAL ARMS IN GALAXIES ORBITING A CLUSTER

ABSTRACT One of the scenarios for the formation of grand-design spiral arms in disky galaxies involves their interactions with a satellite or another galaxy. Here we consider another possibility, where the perturbation is instead due to the potential of a galaxy cluster. Using N-body simulations we investigate the formation and evolution of spiral arms in a Milky-Way-like galaxy orbiting a Virgo-like cluster. The galaxy is placed on a few orbits of different size but similar eccentricity and its evolution are followed for 10 Gyr. The tidally induced, two-armed, approximately logarithmic spiral structure forms on each of them during the pericenter passages. The spiral arms dissipate and wind up with time, to be triggered again at the next pericenter passage. We confirm this transient and recurrent nature of the arms by analyzing the time evolution of the pitch angle and the arm strength. We find that the strongest arms are formed on the tightest orbit; however, they wind up rather quickly and are disturbed by another pericenter passage. The arms on the most extended orbit, which we analyze in more detail, wind up slowly and survive for the longest time. Measurements of the pattern speed of the arms indicate that they are kinematic density waves. We attempt a comparison with observations by selecting grand-design spiral galaxies in the Virgo cluster. Among those, we find nine examples bearing no sign of recent interactions or the presence of companions. For three of them we present close structural analogues among our simulated spiral galaxies.

Tidal Disruption of Milky Way Satellites with Shallow Dark Matter Density Profiles

Dwarf galaxies of the Local Group provide unique possibilities to test current theories of structure formation. Their number and properties have put the broadly accepted cold dark matter model into question, posing a few problems. These problems now seem close to resolution due to the improved treatment of baryonic processes in dwarf galaxy simulations which now predict cored rather than cuspy dark matter profiles in isolated dwarfs with important consequences for their subsequent environmental evolution. Using N-body simulations, we study the evolution of a disky dwarf galaxy with such a shallow dark matter profile on a typical orbit around the Milky Way. The dwarf survives the first pericenter passage but is disrupted after the second due to tidal forces from the host. We discuss the evolution of the dwarf’s properties in time prior to and at the time of disruption. We demonstrate that the dissolution occurs on a rather short timescale as the dwarf expands from a spheroid into a stream with non-zero mean radial velocity. We point out that the properties of the dwarf at the time of disruption may be difficult to distinguish from bound configurations, such as tidally induced bars, both in terms of surface density and line-of-sight kinematics.

Using gas clouds to probe the accretion flow around Sgr A*: G2's delayed pericentre passage

We study the dynamical evolution of the putative gas clouds G1 and G2 recently discovered in the Galactic center. Following earlier studies suggesting that these two clouds are part of a larger gas streamer, we combine their orbits into a single trajectory. Since the gas clouds experience a drag force from background gas, this trajectory is not exactly Keplerian. By assuming the G1 and G2 clouds trace this trajectory, we fit for the drag force they experience and thus extract information about the accretion flow at a distance of thousands of Schwarzschild radii from the black hole. This range of radii is important for theories of black hole accretion, but is currently unconstrained by observations. In this paper we extend our previous work by accounting for radial forces due to possible inflow or outflow of the background gas. Such radial forces drive precession in the orbital plane, allowing a slightly better fit to the G1 and G2 data. This precession delays the pericenter passage of G2 by 4-5 months relative to estimates derived from a Keplerian orbital fit; if it proves possible to identify the pericenter time observationally, this enables an immediate test of whether G1 and G2 are gas clouds part of a larger gas streamer. If G2 is indeed a gas cloud, its closest approach likely occurred in late summer 2014, after many of the observing campaigns monitoring G2's anticipated pericenter passage ended. We discuss how this affects interpretation of the G2 observations.

THE ORBITAL STRUCTURE OF A TIDALLY INDUCED BAR

ABSTRACT Orbits are the key building blocks of any density distribution, and their study helps us understand the kinematical structure and the evolution of galaxies. Here, we investigate orbits in a tidally induced bar of a dwarf galaxy, using an N-body simulation of an initially disky dwarf galaxy orbiting a Milky Way-like host. After the first pericenter passage, a tidally induced bar forms in the stellar component of the dwarf. The bar evolution is different than in isolated galaxies and our analysis focuses on the period before it buckles. We study the orbits in terms of their dominant frequencies, which we calculate in a Cartesian coordinate frame rotating with the bar. Apart from the well-known x1 orbits, we find many other types, mostly with boxy shapes of various degree of elongation. Some of them are also near-periodic, admitting frequency ratios of 4/3, 3/2, and 5/3. The box orbits have various degrees of vertical thickness but only a relatively small fraction of those have banana (i.e., smile/frown) or infinity-symbol shapes in the edge-on view. In the very center we also find orbits known from the potential of triaxial ellipsoids. The elongation of the orbits grows with distance from the center of the bar in agreement with the variation of the shape of the density distribution. Our classification of orbits leads to the conclusion that more than 80% of them have boxy shapes, while only 8% have shapes of classical x1 orbits.

Satellite quenching time-scales in clusters from projected phase space measurements matched to simulated orbits

We measure the star formation quenching efficiency and timescale in cluster environments. Our method uses N-body simulations to estimate the probability distribution of possible orbits for a sample of observed SDSS galaxies in and around clusters based on their position and velocity offsets from their host cluster. We study the relationship between their star formation rates and their likely orbital histories via a simple model in which star formation is quenched once a delay time after infall has elapsed. Our orbit library method is designed to isolate the environmental effect on the star formation rate due to a galaxy's present-day host cluster from `pre-processing' in previous group hosts. We find that quenching of satellite galaxies of all stellar masses in our sample ($10^{9}-10^{11.5}\,{\rm M}_\odot$) by massive ($> 10^{13}\,{\rm M}_\odot$) clusters is essentially $100$ per cent efficient. Our fits show that all galaxies quench on their first infall, approximately at or within a Gyr of their first pericentric passage. There is little variation in the onset of quenching from galaxy-to-galaxy: the spread in this time is at most $\sim 2$ Gyr at fixed $M_*$. Higher mass satellites quench earlier, with very little dependence on host cluster mass in the range probed by our sample.

TIDALLY INDUCED BARS OF GALAXIES IN CLUSTERS

ABSTRACT Using N-body simulations, we study the formation and evolution of tidally induced bars in disky galaxies in clusters. Our progenitor is a massive, late-type galaxy similar to the Milky Way, composed of an exponential disk and a Navarro-Frenk-White dark matter halo. We place the galaxy on four different orbits in a Virgo-like cluster and evolve it for 10 Gyr. As a reference case, we also evolve the same model in isolation. Tidally induced bars form on all orbits soon after the first pericenter passage and survive until the end of the evolution. They appear earlier, are stronger and longer, and have lower pattern speeds for tighter orbits. Only for the tightest orbit are the properties of the bar controlled by the orientation of the tidal torque from the cluster at pericenter. The mechanism behind the formation of the bars is the angular momentum transfer from the galaxy stellar component to its halo. All of the bars undergo extended periods of buckling instability that occur earlier and lead to more pronounced boxy/peanut shapes when the tidal forces are stronger. Using all simulation outputs of galaxies at different evolutionary stages, we construct a toy model of the galaxy population in the cluster and measure the average bar strength and bar fraction as a function of clustercentric radius. Both are found to be mildly decreasing functions of radius. We conclude that tidal forces can trigger bar formation in cluster cores, but not in the outskirts, and thus can cause larger concentrations of barred galaxies toward the cluster center.

Seeding the Galactic Centre gas stream: gravitational instabilities set the initial conditions for the formation of protocluster clouds

Abstract Star formation within the Central Molecular Zone (CMZ) may be intimately linked to the orbital dynamics of the gas. Recent models suggest that star formation within the dust ridge molecular clouds (from G0.253+0.016 to Sgr B2) follows an evolutionary time sequence, triggered by tidal compression during their preceding pericentre passage. Given that these clouds are the most likely precursors to a generation of massive stars and extreme star clusters, this scenario would have profound implications for constraining the time-evolution of star formation. In this Letter, we search for the initial conditions of the protocluster clouds, focusing on the kinematics of gas situated upstream from pericentre. We observe a highly regular corrugated velocity field in {l, vLSR} space, with amplitude and wavelength A = 3.7 ± 0.1 km s−1 and λvel,i = 22.5 ± 0.1 pc, respectively. The extremes in velocity correlate with a series of massive (∼104 M⊙) and compact (Req ∼ 2 pc), quasi-regularly spaced (∼8 pc), molecular clouds. The corrugation wavelength and cloud separation closely agree with the predicted Toomre (∼17 pc) and Jeans (∼6 pc) lengths, respectively. We conclude that gravitational instabilities are driving the condensation of molecular clouds within the Galactic Centre gas stream. Furthermore, we speculate these seeds are the historical analogue of the dust-ridge molecular clouds, representing the initial conditions of star and cluster formation in the CMZ.

MORPHOLOGY AND MOLECULAR GAS FRACTIONS OF LOCAL LUMINOUS INFRARED GALAXIES AS A FUNCTION OF INFRARED LUMINOSITY AND MERGER STAGE

ABSTRACT We present a new, detailed analysis of the morphologies and molecular gas fractions (MGFs) for a complete sample of 65 local luminous infrared galaxies from Great Observatories All-Sky Luminous Infrared Galaxies (LIRG) Survey using high resolution I-band images from The Hubble Space Telescope, the University of Hawaii 2.2 m Telescope and the Pan-STARRS1 Survey. Our classification scheme includes single undisturbed galaxies, minor mergers, and major mergers, with the latter divided into five distinct stages from pre-first pericenter passage to final nuclear coalescence. We find that major mergers of molecular gas-rich spirals clearly play a major role for all sources with however, below this luminosity threshold, minor mergers and secular processes dominate. Additionally, galaxies do not reach until late in the merger process when both disks are near final coalescence. The mean MGF ( ) for non-interacting and early-stage major merger LIRGs is 18 ± 2%, which increases to 33 ± 3%, for intermediate stage major merger LIRGs, consistent with the hypothesis that, during the early-mid stages of major mergers, most of the initial large reservoir of atomic gas (HI) at large galactocentric radii is swept inward where it is converted into molecular gas (H2).

Detection of polarized continuum emission of the Dusty S-cluster Object (DSO/G2)

AbstractA peculiar source in the Galactic center known as the Dusty S-cluster Object (DSO/G2) moves on a highly eccentric orbit around the supermassive black hole with the pericenter passage in the spring of 2014. Its nature has been uncertain mainly because of the lack of any information about its intrinsic geometry. For the first time, we use near-infrared polarimetric imaging data to obtain constraints about the geometrical properties of the DSO. We find out that DSO is an intrinsically polarized source, based on the significance analysis of polarization parameters, with the degree of the polarization of ~30% and an alternating polarization angle as it approaches the position of Sgr A*. Since the DSO exhibits a near-infrared excess of Ks-L′ > 3 and remains rather compact in emission-line maps, its main characteristics may be explained with the model of a pre-main-sequence star embedded in a non-spherical dusty envelope.

Temperature Evolution of Molecular Clouds in the Central Molecular Zone

AbstractWe infer the absolute time dependence of kinematic gas temperature along a proposed orbit of molecular clouds in the Central Molecular Zone (CMZ) of the Galactic Center (GC). Ammonia gas temperature maps are one of the results of the “Survey of Water and Ammonia in the Galactic Center” (SWAG, PI: J. Ott); the dynamical model of molecular clouds in the CMZ was taken from Kruijssen et al. (2015). We find that gas temperatures increase as a function of time in both regimes before and after the cloud passes pericenter on its orbit in the GC potential. This is consistent with the recent proposal that pericenter passage triggers gravitational collapse. Other investigated quantities (line width, column density, opacity) show no strong sign of time dependence but are likely dominated by cloud-to-cloud variations.

Simulations of ram-pressure stripping in galaxy-cluster interactions

Observationally, the quenching of star-forming galaxies appears to depend both on their mass and environment. The exact cause of the environmental dependence is still poorly understood, yet semi-analytic models (SAMs) of galaxy formation need to parameterise it to reproduce observations of galaxy properties. In this work, we use hydrodynamical simulations to investigate the quenching of disk galaxies through ram-pressure stripping (RPS) as they fall into galaxy clusters with the goal of characterising the importance of this effect for the reddening of disk galaxies. Our set-up employs a live model of a galaxy cluster that interacts with infalling disk galaxies on different orbits. We use the moving-mesh code AREPO, augmented with a special refinement strategy to yield high resolution around the galaxy on its way through the cluster in a computationally efficient way. Our direct simulations differ substantially from stripping models employed in current SAMs, which in most cases overpredict the mass loss from RPS. Furthermore, after pericentre passage, as soon as ram pressure becomes weaker, gas that remains bound to the galaxy is redistributed to the outer parts, an effect that is not captured in simplified treatements of RPS. Star formation in our model galaxies is quenched mainly because the hot gas halo is stripped, depriving the galaxy of its gas supply. The cold gas disk is only stripped completely in extreme cases, leading to full quenching and significant reddening on timescale of ~200 Myr. On the other hand, galaxies experiencing only mild ram pressure actually show an enhanced star formation rate that is consistent with observations and are quenched on timescales > 1 Gyr. Stripped gas in the wake is mixed efficiently with intracluster gas already a few tens of kpc behind the disk, and this gas is free of residual star formation.

A Magellanic origin of the DES dwarfs

We establish the connection between the Magellanic Clouds (MCs) and the dwarf galaxy candidates discovered in the Dark Energy Survey (DES) by building a dynamical model of the MC satellite populations, based on an extensive suite of tailor-made numerical simulations. Our model takes into account the response of the Galaxy to the MCs infall, the dynamical friction experienced by the MCs and the disruption of the MC satellites by their hosts. The simulation suite samples over the uncertainties in the MC's proper motions, the masses of the MW and the Clouds themselves and allows for flexibility in the intrinsic volume density distribution of the MC satellites. As a result, we can accurately reproduce the DES satellites' observed positions and kinematics. Assuming that Milky Way (MW) dwarfs follow the distribution of subhaloes in $\Lambda$CDM, we further demonstrate that, of 14 observed satellites, the MW halo contributes fewer than 4 (8) of these with 68% (95%) confidence and that 7 (12) DES dwarfs have probabilities greater than 0.7 (0.5) of belonging to the LMC. Marginalising over the entire suite, we constrain the total number of the Magellanic satellites at ~70, the mass of the LMC around $10^{11}M_\odot$ and show that the Clouds have likely endured only one Galactic pericentric passage so far. Finally, we give predictions for the line-of-sight velocities and the proper motions of the satellites discovered in the vicinity of the LMC.