Jet-like Outflows Research Articles

Post-dynamical inspiral phase of common envelope evolution

During common envelope evolution, an initially weak magnetic field may undergo amplification by interacting with spiral density waves and turbulence generated in the stellar envelope by the inspiralling companion. Using 3D magnetohydrodynamical simulations on adaptively refined spherical grids with excised central regions, we studied the amplification of magnetic fields and their effect on the envelope structure, dynamics, and the orbital evolution of the binary during the post-dynamical inspiral phase. About 95% of magnetic energy amplification arises from magnetic field stretching, folding, and winding due to differential rotation and turbulence while compression against magnetic pressure accounts for the remaining ∼5%. Magnetic energy production peaks at a scale of 3ab, where ab is the semimajor axis of the central binary’s orbit. Because the magnetic energy production declines at large radial scales, the conditions are not favorable for the formation of magnetically collimated bipolar jet-like outflows unless they are generated on small scales near the individual cores, which we did not resolve. Magnetic fields have a negligible impact on binary orbit evolution, mean kinetic energy, and the disk-like morphology of angular momentum transport, but turbulent Maxwell stress can dominate Reynolds stress when accretion onto the central binary is allowed, leading to an α-disk parameter of ≃0.034. Finally, we discovered accretion streams arising from the stabilizing effect of the magnetic tension from the toroidal field about the orbital plane, which prevents overdensities from being destroyed by turbulence and enables them to accumulate mass and eventually migrate toward the binary.

Flares, Jets, and Quasiperiodic Outbursts from Neutron Star Merger Remnants

Using numerical relativity simulations with a subgrid dynamo prescription to generate strong initial magnetic fields, we investigate the possibility of launching a jet-like outflow from the hypermassive neutron star (HMNS) during the early stages of the merger, prior to the remnant’s collapse to a black hole. We demonstrate that buoyant instabilities in the strongly magnetized HMNS can lead to a periodic emission of powerful electromagnetic flares shortly after the merger. These are followed by a collimated mildly relativistic outflow. Both types of outflows feature quasiperiodic kilohertz substructure. These early-time outflows may power precursors to short-duration gamma-ray bursts (sGRBs) or in some cases the entire sGRB. While the overall temporal power spectrum we find broadly agrees with the one recently reported for quasiperiodic oscillations in the sGRB GRB910711, our simulations suggest that the periodic electromagnetic substructure is dominated by magnetohydrodynamic shearing processes rather than correlating with the corresponding postmerger gravitational-wave signal.

Supernova double-peaked light curves from double-nickel distribution

Among supernovae (SNe) of different luminosities, many double-peaked light curves (LCs) have been observed, representing a broad morphological variety. In this work, we investigate which of these can be modelled by assuming a double-peaked distribution of their radioactive material, as originally proposed for SN2005bf. The inner zone corresponds to the regular explosive nucleosynthesis and extends outwards, according to the usual scenario of mixing. The outer 56Ni-rich shell may be related to the effect of jet-like outflows that have interacted with more distant portions of the star before the arrival of the SN shock. As the outer layer is covered by matter that is optically less thick, its energy emerges earlier and generates a first peak of radiation. To investigate this scenario in more detail, we have applied our hydrodynamic code that follows the shock propagation through the progenitor star and takes into account the effect of the γ-ray photons produced by the decay of the radioactive isotopes. We present a simple parametric model for the 56Ni abundance profile and explore the consequences on the LC of individually varying the quantities that define this distribution, setting our focus onto the stripped-envelope progenitors. In this first study, we are interested in the applicability of this model to SNe that have not been classified as superluminous, thus, we have selected our parameter space accordingly. Then, within the same mathematical prescription for the 56Ni -profile, we revisited the modelling process for a series of objects: SN2005bf, PTF2011mnb, SN2019cad, and SN2008D. In some cases, a decrease in the gamma ray opacity is required to fit the late time observations. We also discuss the other cases in which this scenario might be likely to explain the LC morphology. A steep initial decline in the observed bolometric LC within less than few days after the explosion becomes less feasible for this model, because it requires a large abundance of 56Ni near the stellar surface, indicating a strongly inverted distribution. An initial bolometric rise before the two peaks seems more favourable for the double-nickel case, particularly as it can be difficult to explain through other scenarios, unless a combination of power sources is invoked.

ALMA discovery of a dual dense probably rotating outflow from a massive young stellar object G18.88MME

ABSTRACT We report the discovery of a very dense jet-like fast molecular outflow surrounded by a wide-angle wind in a massive young stellar object (MYSO) G18.88MME (stellar mass ∼8 M⊙) powering an extended green object G18.89−0.47. Four cores MM1–4 are identified in the Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum map (resolution ${\sim} 0.8\, \mathrm { arcsec}$) towards G18.88MME, and are seen at the centre of the emission structure (extent ∼0.3 pc × 0.2 pc) detected in the ALMA map. G18.88MME is embedded in the core MM1 (mass ∼13–18 M⊙), where no radio continuum emission is detected. The molecular outflow centred at MM1 is investigated in the SiO(5–4), HC3N(24–23) and 13CO(2–1) lines. The detection of HC3N in the outflow is rare in MYSOs and indicates its very high density. The position–velocity diagrams display a fast narrow outflow (extent ∼28 000 AU) and a slower wide-angle more extended outflow towards MM1, and both of these components show a transverse velocity gradient indicative of a possible rotation. All these observed features together make G18.88MME as a unique object for studying the unification of the jet-driven and wind-driven scenarios of molecular outflows in MYSOs.

Erratum: Simulation of a compact object with outflows moving through a gaseous background

A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are modified by the presence of outflow from the compact object. We show that the destruction of the wake by an outflow reduces dynamical friction, and reverses its sign when the outflow is strong enough, in good quantitative agreement with analytic calculations. For a strong isotropic outflow, the outcome on scales that we have simulated is a negative dynamical friction, i.e., net acceleration. For jet-like outflows driven by reprocessed accretion, both the rate of accretion and the magnitude of dynamical friction drop for more powerful jets. The accretion rate is strongly intermittent when the jet points to the same direction as the motion of the compact object. The dynamical effects of outflows may be important for the evolution of compact objects during the common envelope phase of binary systems, and for accreting compact objects and massive stars encountering AGN discs.

Exact solution of partial differential equations for the creation of jet-like flows in plasmas and neutral fluids

An exact solution of two fluid ideal classical plasma equations is presented which shows that the one-dimensional jet-like axial outflow and two-dimensional magnetic field are generated simultaneously by the density and temperature gradients of both electrons and ions in cylindrical geometry. Particular profiles of density function ψ=ln n¯ (where n¯ is normalized by some constant density N0) and temperatures Tj (for j=e,i) have to be chosen to obtain an exact solution of the complicated nonlinear partial differential equations. But this is a natural analytical exact solution of the ideal two fluid plasma equations. The basic mechanism presented here explains the creation of plasma jet-like flows along with magnetic fields in astrophysical environments such as young stellar objects, active galactic nuclei, solar spicules, flares, and coronal loops. The theoretical model is also applicable to laser induced plasma where magnetic field and plasma ablation are produced simultaneously. An exact analytical solution of ideal neutral fluid equations is also presented which shows that the jet-like outflows can be generated by the density and temperature gradients in such systems.

Bipolar planetary nebulae from outflow collimation by common envelope evolution

ABSTRACT The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and the dense toroidal-shaped ejecta left over from common envelope (CE) evolution. Here we use the 3D hydrodynamic adaptive mesh refinement (AMR) code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the smoothed particle hydrodynamics (SPH) code phantom serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta pre-PN outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post-CE winds may be aspherical, our models show that collimation via ‘inertial confinement’ will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock-focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC-driven bipolar lobes.

Simulation of a compact object with outflows moving through a gaseous background

ABSTRACT A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi–Hoyle–Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are modified by the presence of outflow from the compact object. We show that the destruction of the wake by an outflow reduces dynamical friction, and reverses its sign when the outflow is strong enough, in good quantitative agreement with analytic calculations. For a strong isotropic outflow, the outcome on scales that we have simulated is a negative dynamical friction, i.e. net acceleration. For jet-like outflows driven by reprocessed accretion, both the rate of accretion and the magnitude of dynamical friction drop for more powerful jets. The accretion rate is strongly intermittent when the jet points to the same direction as the motion of the compact object. The dynamical effects of outflows may be important for the evolution of compact objects during the common envelope phase of binary systems, and for accreting compact objects and massive stars encountering active galactic nucleus discs.

Fine-scale Explosive Energy Release at Sites of Prospective Magnetic Flux Cancellation in the Core of the Solar Active Region Observed by Hi-C 2.1, IRIS, and SDO

Abstract The second Hi-C flight (Hi-C 2.1) provided unprecedentedly high spatial and temporal resolution (∼250 km, 4.4 s) coronal EUV images of Fe ix/x emission at 172 Å of AR 12712 on 2018 May 29, during 18:56:21–19:01:56 UT. Three morphologically different types (I: dot-like; II: loop-like; III: surge/jet-like) of fine-scale sudden-brightening events (tiny microflares) are seen within and at the ends of an arch filament system in the core of the AR. Although type Is (not reported before) resemble IRIS bombs (in size, and brightness with respect to surroundings), our dot-like events are apparently much hotter and shorter in span (70 s). We complement the 5 minute duration Hi-C 2.1 data with SDO/HMI magnetograms, SDO/AIA EUV images, and IRIS UV spectra and slit-jaw images to examine, at the sites of these events, brightenings and flows in the transition region and corona and evolution of magnetic flux in the photosphere. Most, if not all, of the events are seated at sites of opposite-polarity magnetic flux convergence (sometimes driven by adjacent flux emergence), implying likely flux cancellation at the microflare’s polarity inversion line. In the IRIS spectra and images, we find confirming evidence of field-aligned outflow from brightenings at the ends of loops of the arch filament system. In types I and II the explosion is confined, while in type III the explosion is ejective and drives jet-like outflow. The light curves from Hi-C, AIA, and IRIS peak nearly simultaneously for many of these events, and none of the events display a systematic cooling sequence as seen in typical coronal flares, suggesting that these tiny brightening events have chromospheric/transition region origin.

Black hole magnetosphere with small-scale flux tubes – II. Stability and dynamics

ABSTRACT In some Seyfert galaxies, the hard X-rays that produce fluorescent emission lines are thought to be generated in a hot corona that is compact and located at only a few gravitational radii above the supermassive black hole. We consider the possibility that this X-ray source may be powered by small-scale magnetic flux tubes attached to the accretion disc near the black hole. We use three-dimensional, time-dependent, special relativistic, force-free simulations in a simplified setting to study the dynamics of such flux tubes as they get continuously twisted by the central compact star/black hole. We find that the dynamical evolution of the flux tubes connecting the central compact object and the accretion disc is strongly influenced by the confinement of the surrounding field. Although differential rotation between the central object and the disc tends to inflate the flux tubes, strong confinement from surrounding field quenches the formation of a jet-like outflow, as the inflated flux tube becomes kink unstable and dissipates most of the extracted rotational energy relatively close to the central object. Such a process may be able to heat up the plasma and produce strong X-ray emission. We estimate the energy dissipation rate and discuss its astrophysical implications.

High-velocity Bullets from V Hydrae, an Asymptotic Giant Branch Star in Transition: Ejection History and Spatio-kinematic Modeling

Abstract The carbon star V Hydrae (V Hya) provides new insight into the nature of the launching mechanism of jet-like outflows that are believed to be the cause of the poorly understood transition phase of asymptotic giant branch stars into aspherical planetary nebulae. V Hya has been shown to periodically eject collimated gas blobs at high velocities (“bullets”). By analyzing data from Hubble Space Telescope/Space Telescope Imaging Spectrograph 2D spectra, obtained at six epochs spaced over a decade that show four successively ejected bullets with a spacing of 8.5 yr, we have created kinematic models of the dynamical evolution of a specific bullet (#1) for the first three observed epochs (2002, 2003, 2004) using a 3D spatio-kinematic code, SHAPE. Using these models, we fit the observed morphology, line-of-sight velocity, proper motion, and intensity for the extended, gaseous bullet as a function of time over a period of 2 yr, in order to constrain its 3D movement and the evolution of its physical properties over this period. Our results suggest that although bullet #1’s motion is predominantly ballistic, there are small but significant changes in the position angle and inclination angle of the long (symmetry) axis of the bullet that tilt it progressively toward the symmetry axis of the bipolar molecular nebula around V Hya. In contrast, bullet #3 shows strong acceleration soon after ejection. We discuss the possibilities that bullet acceleration is caused by a nonradial magnetic field and/or by hydrodynamic interaction with the ambient gas through which the bullet is traveling.

Planetary Nebulae Shaped by Common Envelope Evolution

The morphologies of planetary nebula have long been believed to be due to wind shaping processes in which a “fast wind” from the central star impacts a previously ejected envelope. It is assumed that asymmetries existing in the “slow wind” envelope would lead to inertial confinement, shaping the resulting interacting wind flow. We present new results demonstrating the effectiveness of Common Envelope Evolution (CEE) at producing aspherical envelopes which, when impinged upon by a spherical fast stellar wind, produce highly bipolar, jet-like outflows. We have run two simple cases using the output of a single PHANTOM SPH CEE simulation. Our work uses the Adaptive Mesh Refinement code AstroBEAR to track the interaction of the fast wind and CEE ejecta allows us to follow the morphological evolution of the outflow lobes at high resolution in 3-D. Our two models bracket low and high momentum output fast winds. We find the interaction leads to highly collimated bipolar outflows. In addition, the bipolar morphology depends on the fast wind momentum injection rate. With this dependence comes the initiation of significant symmetry breaking between the top and bottom bipolar lobes. Our simulations, though simplified, confirm the long-standing belief that CEE can plan a major role in PPN and PN shaping. These simulations are intended as an initial exploration of the post-CE/PPN flow patterns that can be expected from central source outflows and CE ejecta.

ALMA Observations of the Water Fountain Pre-Planetary Nebula IRAS 16342-3814: High-Velocity Bipolar Jets and an Expanding Torus.

We have mapped 12CO J=3-2 and other molecular lines from the "water-fountain" bipolar pre-planetary nebula (PPN) IRAS 16342-3814 with [Formula: see text] resolution using ALMA. We find (i) two very high-speed knotty, jet-like molecular outflows, (ii) a central high-density (> few × 106 cm-3), expanding torus of diameter 1300 AU, and (iii) the circumstellar envelope of the progenitor AGB, generated by a sudden, very large increase in the mass-loss rate to > 3.5 × 10-4M⊙ yr-1 in the past ~455 yr. Strong continuum emission at 0.89 mm from a central source (690 mJy), if due to thermally-emitting dust, implies a substantial mass (0.017 M⊙) of very large (~mm-sized) grains. The measured expansion ages of the above structural components imply that the torus (age~160 yr) and the younger high-velocity outflow (age~110 yr) were formed soon after the sharp increase in the AGB mass-loss rate. Assuming a binary model for the jets in IRAS 16342, the high momentum rate for the dominant jet-outflow in IRAS 16342 implies a high minimum accretion rate, ruling out standard Bondi-Hoyle-Lyttleton wind accretion and wind Roche lobe overflow (RLOF) models with white-dwarf or main-sequence companions. Most likely, enhanced RLOF from the primary or accretion modes operating within common envelope evolution are needed.

Further evidence for a quasar-driven jet impacting its neighbour galaxy: The saga of HE0450-2958 continues

HE0450-2958, an interacting quasar-starburst galaxy pair at $z$ = 0.285, is one of the best known examples of strong star formation activity in the presence of a quasar-driven jet. We present new multi-band JVLA-imaging covering 1 to 6 GHz and reaching an angular resolution of up to $0{}_{{}^.}^{"}6$ (a 6-fold improvement over existing radio data). We confirm the previous detection of a spatially extended radio component around the quasar indicating that there is on-going star formation activity in the quasar host galaxy. For the first time, we directly detect a jet-like bipolar outflow from the quasar aligned with its companion star-forming galaxy (SFG) and several blobs of ionized gas in its vicinity identified in previous studies. Within the companion SFG we find evidence for a flattening of the synchrotron spectral index towards the point of intersection with the jet axis, further suggesting that the outflow may actually be impacting its interstellar medium (ISM). We discuss two possible mechanisms that could have triggered the starburst in the companion SFG: a wet-dry merger with the quasar and jet-induced star formation. While triggering through interaction-driven gas dynamics cannot be excluded with current data, our new observations make HE0450-2958 a strong candidate for jet-induced star formation, and one of the rare links between local systems (like Minkowski's Object or Centaurus A) and the high-z regime where radio-optical alignments suggest that this phenomenon could be more common.

Erratum and Addendum: Smoothed particle magnetohydrodynamic simulations of protostellar outflows with misaligned magnetic field and rotation axes

We have developed a modified form of the equations of smoothed particle magnetohydrodynamics which are stable in the presence of very steep density gradients. Using this formalism, we have performed simulations of the collapse of magnetised molecular cloud cores to form protostars and drive outflows. Our stable formalism allows for smaller sink particles (< 5 AU) than used previously and the investigation of the effect of varying the angle, {\theta}, between the initial field axis and the rotation axis. The nature of the outflows depends strongly on this angle: jet-like outflows are not produced at all when {\theta} > 30{\deg}, and a collimated outflow is not sustained when {\theta} > 10{\deg}. No substantial outflows of any kind are produced when {\theta} > 60{\deg}. This may place constraints on the geometry of the magnetic field in molecular clouds where bipolar outflows are seen.

DISCOVERY OF AN EXTREMELY WIDE-ANGLE BIPOLAR OUTFLOW IN AFGL 5142

ABSTRACT Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle (∼180°) bipolar outflow (“EWBO”) in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3–2) and HCO+ (3–2) lines. This bipolar outflow is along a north-west to south-east direction with a line of sight flow velocity of about 3 km s−1 and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high-velocity jet-like outflows and the “EWBO” are ∼24.5 M ⊙ and ∼1.7 × 10−3 M ⊙ yr−1, respectively. Global collapse of the clump is revealed by the “blue profile” in the HCO+ (1–0) line. A hierarchical network of filaments was identified in NH3 (1, 1) emission. Clear velocity gradients of the order of 10 km s−1 pc−1 are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is ∼3.1 × 10−3 M ⊙ yr−1, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

COORDINATED X-RAY, ULTRAVIOLET, OPTICAL, AND RADIO OBSERVATIONS OF THE PSR J1023+0038 SYSTEM IN A LOW-MASS X-RAY BINARY STATE

The PSR J1023+0038 binary system hosts a neutron star and a low-mass, main-sequence-like star. It switches on year timescales between states as an eclipsing radio millisecond pulsar and a low-mass X-ray binary (LMXB). We present a multi-wavelength observational campaign of PSR J1023+0038 in its most recent LMXB state. Two long XMM-Newton observations reveal that the system spends ~70% of the time in a ≈3 × 10^(33) erg s^(−1) X-ray luminosity mode, which, as shown in Archibald et al., exhibits coherent X-ray pulsations. This emission is interspersed with frequent lower flux mode intervals with ≈5 x 10^(32) erg s^(−1) and sporadic flares reaching up to ≈10^(34) erg s^(−1), with neither mode showing significant X-ray pulsations. The switches between the three flux modes occur on timescales of order 10 s. In the UV and optical, we observe occasional intense flares coincident with those observed in X-rays. Our radio timing observations reveal no pulsations at the pulsar period during any of the three X-ray modes, presumably due to complete quenching of the radio emission mechanism by the accretion flow. Radio imaging detects highly variable, flat-spectrum continuum radiation from PSR J1023+0038, consistent with an origin in a weak jet-like outflow. Our concurrent X-ray and radio continuum data sets do not exhibit any correlated behavior. The observational evidence we present bears qualitative resemblance to the behavior predicted by some existing propeller and trapped disk accretion models although none can account for key aspects of the rich phenomenology of this system.

High-velocity Wind from IRS 1 in the NGC 2071IR

AbstractWe present the results of 1.3 and 3.6 cm radio continuum emission toward the NGC 2071IR star-forming region, carried out with the VLA in its A configuration. We detect continuum emission toward the infrared sources IRS 1 and IRS 3 at both wavelengths. In particular, IRS 1 breaks up into three continuum peaks (IRS 1E, 1C, and 1W), aligned in the east-west direction, being IRS 1 the central source. The morphology of the condensation IRS 1W is very interesting, which has an elongated structure and shows a significant curvature towards the north. We suggest that this morphology could be explained as the impact of a high-velocity wind or jetlike outflow from IRS 1 on a close companion or other obstruction, which also explains the strong water maser emission observed toward IRS 1W.

Explaining the flat-spectrum radio core Sgr A* with GRMHD simulations of jets

AbstractThe supermassive black hole in the center of the Milky Way, Sgr A*, displays a nearly flat radio spectrum which is typical for jets in Active Galactic Nuclei. Indeed, time dependent, magnetized models of radiatively inefficient accretion flows, which are commonly used to explain the millimeter, near-infrared, and X-ray emission of Sgr A* also often produce jet-like outflows. However, the emission from these models so far has failed to reproduce the flat radio spectrum. We show that current GRMHD simulations can naturally reproduce the flat spectrum, when using a two-temperature plasma in the disk and a constant electron temperature plasma in the jet. This assumption is consistent with current state-of-the art simulations, in which the electron temperature evolution is not explicitly modeled. Stronger magnetization and stronger shearing seen in the jet sheath could possibly explain the difference in electron heating between jet and disk. The model images and spectra are consistent with the radio sizes and spectrum of Sgr A*.

Kinetic or thermal AGN feedback in simulations of isolated and merging disc galaxies calibrated by the M-σ relation

(Abridged) We investigate two modes of coupling the feedback energy from a central AGN to the neighboring gas in galaxy simulations: kinetic - velocity boost, and thermal - heating. We formulate kinetic feedback models for energy-driven wind (EDW) and momentum-driven wind (MDW), using two free parameters: feedback efficiency epsilon_f, and AGN wind velocity v_w. A novel numerical algorithm is implemented in the SPH code GADGET-3, to prevent the expansion of a hole in the gas distribution around the BH. We perform simulations of isolated evolution and merger of disk galaxies, of Milky-Way mass as well as lower and higher masses. We find that in the isolated galaxy BH kinetic feedback generates intermittent bipolar jet-like gas outflows. We infer that current prescriptions for BH subgrid physics in galaxy simulations can grow the BH to observed values even in an isolated disk galaxy. The BH growth is enhanced in a galaxy merger. Comparing the [M_BH - sigma_star] relation obtained in our simulations with observational data, we conclude that it is possible to find parameter sets for a fit in all the models, except for the case with MDW feedback in a galaxy merger, in which the BH is always too massive. The BH thermal feedback implementation of Springel, Di Matteo & Hernquist (2005) within the multiphase star-formation model is found to have negligible impact on gas properties; and the effect claimed in all previous studies is attributed to gas depletion around the BH by the creation of an artificial hole. The BH mass accretion rate in our simulations exhibit heavy fluctuations. The star formation rate is quenched with feedback by removal of gas. The CGM gas at galactocentric distances (20 - 100)/h kpc are found to give the best metallicity observational diagnostic to distinguish between BH models.