Showing 1–50 of 66 results for author: Fendt, C

Thin Accretion disks in GR-MHD simulations

Authors: Indu K. Dihingia, Christian Fendt

Abstract: We review some recent results of general relativistic magnetohydrodynamic (GR-MHD) simulations considering the evolution of geometrically thin disks around a central black hole. Thin disk GR-MHD simulations complement the widely used MAD (Magnetically Arrested Disk) or SANE (Standard And Normal Evolution) approaches of evolving from an initial disk torus. In particular, we discuss the dynamical ev… ▽ More We review some recent results of general relativistic magnetohydrodynamic (GR-MHD) simulations considering the evolution of geometrically thin disks around a central black hole. Thin disk GR-MHD simulations complement the widely used MAD (Magnetically Arrested Disk) or SANE (Standard And Normal Evolution) approaches of evolving from an initial disk torus. In particular, we discuss the dynamical evolution of the disk, its role in the formation of disk winds or jets, the impact of disk resistivity, and its potential role in generating magnetic flux by an internal disk dynamo. The main characteristics of a thin disk in our approach are the Keplerian rotation of the disk material, which allows to launch disk outflows by the Blandford-Payne magneto-centrifugal effect, in addition to the Blandford-Znajek-driven spine jet from the black hole ergosphere. Thus, for this approach, we neglect disk thermodynamics and radiative effects, concentrating predominantly on the dynamical evolution of the system. Resistive MHD further allows the investigation of physical reconnection and also dynamo action. Magnetic reconnection may generate magnetic islands of plasmoids that are ejected from the disk along with the outflow. We also discussed potential applications of thin disk in explaining the decaying phase of an outburst in black hole X-ray binaries (BH-XRBs). Post-processing of radiation using the simulated dynamical data allows to derive spectra or fluxes, e.g., in the X-ray band, and to derive potential variability characteristics. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 37 pages, 15 figures, Invited chapter for the edited book New Frontiers in GRMHD Simulations (Eds. C. Bambi, Y. Mizuno, S. Shashank and F. Yuan, Springer Singapore, expected in 2024)

X-Shaped Radio Galaxies: Probing Jet Evolution, Ambient Medium Dynamics, and Their Intricate Interconnection

Authors: Gourab Giri, Christian Fendt, Kshitij Thorat, Gianluigi Bodo, Paola Rossi

Abstract: This review explores the field of X-shaped radio galaxies (XRGs), a distinctive subset of winged radio sources that are identified by two pairs of jetted lobes which aligned by a significant angle, resulting in an inversion-symmetric structure. These lobes, encompassing active (primary) and passive (secondary) phases, exhibit a diverse range of properties across the multiple frequency bands, posin… ▽ More This review explores the field of X-shaped radio galaxies (XRGs), a distinctive subset of winged radio sources that are identified by two pairs of jetted lobes which aligned by a significant angle, resulting in an inversion-symmetric structure. These lobes, encompassing active (primary) and passive (secondary) phases, exhibit a diverse range of properties across the multiple frequency bands, posing challenges in discerning their formation mechanism. The proposed mechanisms can broadly be categorized into those related either to a triaxial ambient medium, into which the jet propagates, or to a complex, central AGN mechanism, where the jet is generated. The observed characteristics of XRGs as discovered in the most substantial sample to date, challenge the idea that there is universal process at work that produces the individual sources of XRGs. Instead, the observational and numerical results rather imply the absence of an universal model and infer that distinct mechanisms may be at play for the specific sources. By scrutinizing salient and confounding properties, this review intends to propose the potential direction for future research to constrain and constrict individual models applicable to XRGs. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Accepted for publication in "Frontiers in Astronomy and Space Sciences"

Deciphering the morphological origins of X-shaped radio galaxies: Numerical modeling of Back-flow vs. Jet-reorientation

Authors: Gourab Giri, Bhargav Vaidya, Christian Fendt

Abstract: X-shaped Radio Galaxies (XRGs) develop when certain extra-galactic jets deviate from their propagation path. An asymmetric ambient medium (Back-flow model) or complex Active Galactic Nuclei activity (Jet-reorientation model) enforcing the jet direction to deviate may cause such structures. In this context, the present investigation focuses on the modeling of XRGs by performing 3D relativistic magn… ▽ More X-shaped Radio Galaxies (XRGs) develop when certain extra-galactic jets deviate from their propagation path. An asymmetric ambient medium (Back-flow model) or complex Active Galactic Nuclei activity (Jet-reorientation model) enforcing the jet direction to deviate may cause such structures. In this context, the present investigation focuses on the modeling of XRGs by performing 3D relativistic magneto-hydrodynamic simulations. We implement different jet propagation models applying an initially identical jet-ambient medium configuration to understand distinctive features. This study, the first of its kind, demonstrates that all adopted models produce XRGs with notable properties, thereby challenging the notion of a universal model. Jet reorientation naturally explains several contentious properties of XRGs, including wing alignment along the ambient medium's primary axis, development of collimated lobes, and the formation of noticeably longer wings than active lobes. Such XRGs disrupt the cluster medium by generating isotropic shocks and channeling more energy than the Back-flow scenario. Our synthetic thermal X-ray maps of the cluster medium reveal four `clear' elongated cavities associated with the wing-lobe alignment, regardless of projection effects, though affecting their age estimation. We show depth and geometric alignment of the evolved cavities may qualify as promising characteristics of XRGs, which may be used to disentangle different formation scenarios. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series (10 Figs, 2 Tables)

Precession-induced Variability in AGN Jets and OJ 287

Authors: Silke Britzen, Michal Zajaček, Gopal-Krishna, Christian Fendt, Emma Kun, Frédéric Jaron, Aimo Sillanpää, Andreas Eckart

Abstract: The combined study of the flaring of Active Galactic Nuclei (AGN) at radio wavelengths and pc-scale jet kinematics with Very Long Baseline Interferometry (VLBI) has led to the view that i) the observed flares are associated with ejections of synchrotron blobs from the core, and ii) most of the flaring would follow a one-to-one correlation with the component ejection. Recent results have provided m… ▽ More The combined study of the flaring of Active Galactic Nuclei (AGN) at radio wavelengths and pc-scale jet kinematics with Very Long Baseline Interferometry (VLBI) has led to the view that i) the observed flares are associated with ejections of synchrotron blobs from the core, and ii) most of the flaring would follow a one-to-one correlation with the component ejection. Recent results have provided mounting evidence that the quasi-regular component injections into the relativistic jet may not be the only cause of the flux variability. We propose that AGN flux variability and jet morphology changes can both be of deterministic nature, i.e. having a geometric/kinetic origin linked to the time-variable Doppler beaming of the jet emission as its direction changes due to precession (and nutation). The physics of the underlying jet leads to shocks, instabilities, or to ejections of plasmoids. The appearance (morphology, flux, etc.) of the jet can, however, be strongly affected and modulated by precession. We demonstrate this modulating power of precession for OJ 287. For the first time, we show that the spectral state of the Spectral Energy Distribution (SED) can be directly related to the jet's precession phase. We model the SED evolution and reproduce the precession parameters. Further, we apply our precession model to eleven prominent AGN. We show that for OJ 287 precession seems to dominate the long-term variability ($\gtrsim 1\,{\rm yr}$) of the AGN flux, SED spectral state, and jet morphology, while stochastic processes affect the variability on short timescales ($\lesssim 0.2\,{\rm yr}$). △ Less

Submitted 13 July, 2023; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: 48 pages, 26 figures, 14 tables; published in the Astrophysical Journal

Journal ref: The Astrophysical Journal (2023), Volume 951, Issue 2, id.106, 38 pp

Particles in Relativistic MHD Jets. I. Role of Jet Dynamics in Particle Acceleration

Authors: Ravi Pratap Dubey, Christian Fendt, Bhargav Vaidya

Abstract: Relativistic jets from (supermassive) black holes are typically observed in non-thermal emission, caused by highly-relativistic electrons. Here, we study the interrelation between three-dimensional (special) relativistic magnetohydrodynamics, and particle acceleration in these jets. We inject Lagrangian particles into the jet that are accelerated through diffusive shock acceleration and radiate en… ▽ More Relativistic jets from (supermassive) black holes are typically observed in non-thermal emission, caused by highly-relativistic electrons. Here, we study the interrelation between three-dimensional (special) relativistic magnetohydrodynamics, and particle acceleration in these jets. We inject Lagrangian particles into the jet that are accelerated through diffusive shock acceleration and radiate energy via synchrotron and inverse Compton processes. We investigate the impact of different injection nozzles on the jet dynamics, propagation, and the spectral energy distribution of relativistic particles. We consider three different injection nozzles -- injecting steady, variable and precessing jets. These jets evolve with substantially different dynamics, driving different levels of turbulence and shock structures. The steady jet shows a strong, stationary shock feature, resulting from a head-on collision with an inner back-flow along the jet axis - a jet inside a jet. This shock represents a site for highly-efficient particle acceleration for electrons upto a few tens of TeV and should be visible in emission as a jet knot. Overall, we find that the total number of shocks is more essential for particle acceleration than the strength of the shocks. The precessing jet is most efficient in accelerating electrons to high energies reaching even few hundred TeVs, with power-law index ranging from 2.3 to 3.1. We compare different outflow components, such as jet and the entrained material concerning particle acceleration. For the precessing nozzle, particle acceleration in the entrained material is as efficient as in the jet stream. This is due to the higher level of turbulence induced by the precession motion. △ Less

Submitted 13 July, 2023; v1 submitted 19 June, 2023; originally announced June 2023.

Comments: 26 pages, 13 figures, published in the Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 952, Number 1, 2023

Signatures of winds and jets in the environment of supermassive black holes

Authors: Bidisha Bandyopadhyay, Christian Fendt, Dominik R. G. Schleicher, Javier Lagunas, Javier Pedreros, Neil M. Nagar, Felipe Agurto

Abstract: The Event Horizon Telescope Collaboration (EHTC) has presented first - dynamic-range limited - images of the black hole shadows in M87 and Sgr A*. The next generation Event Horizon Telescope (ngEHT) will provide higher sensitivity and higher dynamic range images (and movies) of these two sources plus image at least a dozen others at $\leq$100 gravitational radii resolution. We here perform an expl… ▽ More The Event Horizon Telescope Collaboration (EHTC) has presented first - dynamic-range limited - images of the black hole shadows in M87 and Sgr A*. The next generation Event Horizon Telescope (ngEHT) will provide higher sensitivity and higher dynamic range images (and movies) of these two sources plus image at least a dozen others at $\leq$100 gravitational radii resolution. We here perform an exploratory study of the appearance of winds and jets in such future observations. To do this we use M87 and Sgr A* as reference systems: we do not aim to exactly reproduce them, but rather to determine how their observed images will depend on specific physical assumptions. Even in the case of similar or the same dynamics, the images depend significantly on global parameters such as the black hole mass and the mass accretion rate. Our results provide guidance in the interpretation of future high-resolution images, particularly if a wind or jet is detected. △ Less

Submitted 12 March, 2023; originally announced March 2023.

Comments: 12 pages, 5 figures, Submitted to MNRAS

Forbidden emission lines in protostellar outflows and jets with MUSE

Authors: Lizxandra Flores-Rivera, Mario Flock, Nicolás Kurtovic, Bernd Husemann, Andrea Banzatti, Simon C. Ringqvist, Sebastian Kamann, André Müller, Christian Fendt, Rebeca Garcia Lopez, Gabriel-Dominique Marleau, Thomas Henning, Carlos Carrasco-Gonzalez, Roy van Boekel, Miriam Keppler, Ralf Launhardt, Yuhiko Aoyama

Abstract: Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral fiel… ▽ More Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral field spectrograph of the Multi Unit Spectroscopic Explorer (MUSE), at the Very Large Telescope (VLT). With a resolution of 0.025 X 0.025 arcsec$^{2}$, we aim to derive the position angle of the outflow/jet (PA$_{outflow/jet}$) that is connected with the inner disk. The forbidden emission lines analyzed here have their origin at the inner parts of the protoplanetary disk. From the maximum intensity emission along the outflow/jet in DL Tau, CI Tau, DS Tau, IP Tau, and IM Lup, we were able to reliably measure the PA$_{outflow/jet}$ for most of the identified lines. We found that our estimates agree with PA$_{dust}$ for most of the disks. These estimates depend on the signal-to-noise level and the collimation of the outflow (jet). The outflows/jets in CIDA 9, GO Tau, and GW Lup are too compact for a PA$_{outflow/jet}$ to be estimated. Based on our kinematics analysis, we confirm that DL Tau and CI Tau host a strong outflow/jet with line-of-sight velocities much greater than 100 km s$^{-1}$, whereas DS Tau, IP Tau, and IM Lup velocities are lower and their structures encompass low-velocity components to be more associated with winds. Our estimates for the mass-loss rate, $\dot{M}_{loss}$, range between (1.1-6.5)X10$^{-7}$-10$^{-8}$ $M_{\odot}$ yr$^{-1}$ for the disk-outflow/jet systems analyzed here. The outflow/jet systems analyzed here are aligned within around 1 degree between the inner and outer disk. Further observations are needed to confirm a potential misalignment in IM Lup. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: Accepted on December 19 2022

Truncated accretion discs in black hole X-ray binaries: dynamics and variability signatures

Authors: Indu K. Dihingia, Bhargav Vaidya, Christian Fendt

Abstract: Variable features in black hole X-ray Binaries (BH-XRBs) are observed in different energy ranges and time scales. The physical origin of different spectral states in BH-XRBs and their relations with the underlying accretion disc are still elusive. To investigate the intermediate state of BH-XRBs during outburst, we simulate a truncated accretion disc around a Kerr black hole using a general relati… ▽ More Variable features in black hole X-ray Binaries (BH-XRBs) are observed in different energy ranges and time scales. The physical origin of different spectral states in BH-XRBs and their relations with the underlying accretion disc are still elusive. To investigate the intermediate state of BH-XRBs during outburst, we simulate a truncated accretion disc around a Kerr black hole using a general relativistic magneto-hydrodynamical (GRMHD) framework under axisymmetry with adaptively refined mesh. Additionally, we have also carried out radiative transfer calculations for understanding the implications of disc dynamics on emission. Dynamically, the inner edge of the truncated accretion disc oscillates in a quasi-periodic fashion (QPO). The QPO frequency of oscillations $(ν_{\rm QPO, max})$ increases as the magnetic field strength and magnetic resistivity increase. However, as the truncation radius increases, $ν_{\rm QPO, max}$ decreases. In our simulation models, frequency varies between $7\times(10M_{\odot}/M_{\rm BH})$ Hz $\lesssimν_{\rm QPO, max}\lesssim20 \times (10M_{\odot}/M_{\rm BH})$ Hz, which is in the range of low-frequency QPOs. We further find evidence of transient shocks in the highly accreting stage during oscillation. Such a transient shock acts as an extended hot post-shock corona around the black hole that has an impact on its radiative properties. The radiative transfer calculations show signatures of these oscillations in the form of modulation in the edge-brightened structure of the accretion disc. △ Less

Submitted 18 October, 2022; v1 submitted 23 August, 2022; originally announced August 2022.

Comments: 18 pages, 24 figures, accepted for publication in MNRAS, Arxiv version is same as the published version

Jets from accretion disk dynamos: consistent quenching modes for dynamo and resistivity

Authors: Giancarlo Mattia, Christian Fendt

Abstract: Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the magnetic field, which is a key component of the launching process, is still an open question. Here we address the question of how the magnetic field required for jet launching is generated and maintained by a dynamo process. By carrying out non-ideal MHD simu… ▽ More Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the magnetic field, which is a key component of the launching process, is still an open question. Here we address the question of how the magnetic field required for jet launching is generated and maintained by a dynamo process. By carrying out non-ideal MHD simulations (PLUTO code), we investigate how the feedback of the generated magnetic field on the mean-field dynamo affects the disk and jet properties. We find that a stronger quenching of the dynamo leads to a saturation of the magnetic field at a lower disk magnetization. Nevertheless, we find that, while applying different dynamo feedback models, the overall jet properties remain unaffected. We then investigate a feedback model which encompasses a quenching of the magnetic diffusivity. Our modeling considers a more consistent approach for mean-field dynamo modeling simulations, as the magnetic quenching of turbulence should be considered for both, a turbulent dynamo and turbulent magnetic diffusivity. We find that, after the magnetic field is saturated, the Blandford-Payne mechanism can work efficiently, leading to more collimated jets, that move, however, with slower speed. We find strong intermittent periods of flaring and knot ejection for low Coriolis numbers. In particular, flux ropes are built up and advected towards the inner disk thereby cutting off of the inner disk wind, leading to magnetic field reversals, reconnection and the emergence of intermittent flares. △ Less

Submitted 2 July, 2022; originally announced July 2022.

Comments: 24 pages, 17 figures, accepted for publication in ApJ

Journal ref: 2022 ApJ 935 22

Curved jet motion. I. Orbiting and precessing jets

Authors: Christian Fendt, Melis Yardimci

Abstract: Astrophysical jets are often observed as bent or curved structures. We also know that the different jet sources may be binary in nature, which may lead to a regular, periodic motion of the jet nozzle, an orbital motion or precession. Here, we present the results of 2D (M)HD simulations in order to investigate how a precessing or orbiting jet nozzle affects the propagation of a high-speed jet. We h… ▽ More Astrophysical jets are often observed as bent or curved structures. We also know that the different jet sources may be binary in nature, which may lead to a regular, periodic motion of the jet nozzle, an orbital motion or precession. Here, we present the results of 2D (M)HD simulations in order to investigate how a precessing or orbiting jet nozzle affects the propagation of a high-speed jet. We have performed a parameter study of systems with different precession angles, orbital periods or separations, and different magnetic field strengths. We find that these kinds of nozzles lead to curved jet propagation which is determined by the main parameters that define the jet nozzle. We find C-shaped jets from orbiting nozzles and S-shaped jets from precessing nozzles. Over long time and long distances, the initially curved jet motion bores a broad channel into the ambient gas that is filled with high-speed jet material which lateral motion is damped, however. A strong (longitudinal) magnetic field can damp the jet curvature that is enforced by either precession of orbital motion of the jet sources. We have investigated the force balance across the jet and ambient medium and found that the lateral magnetic pressure and gas pressure gradients are almost balanced, but that a lack of gas pressure on the concave side of the curvature is leading to the lateral motion. Magnetic tension does not play a significant role. Our results are obtained in code units, but we provide scaling relations such that our results may be applied to young stars, micro-quasars, symbiotic stars or AGN. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Comments: 20 pp, 14 figs; accepted by ApJ

The physics of the MHD disk-jet transition in binary systems: jetted spiral walls launched from disk spiral arms

Authors: Somayeh Sheikhnezami, Christian Fendt

Abstract: We present a detailed physical analysis of the jet launching mechanism of a circum-stellar disk that is located in a binary system. Applying 3D resistive MHD simulations, we investigate the local and global properties of the system, such as the angular momentum transport and the accretion and ejection mass fluxes. In comparison to previous works, for the first time, we have considered the full mag… ▽ More We present a detailed physical analysis of the jet launching mechanism of a circum-stellar disk that is located in a binary system. Applying 3D resistive MHD simulations, we investigate the local and global properties of the system, such as the angular momentum transport and the accretion and ejection mass fluxes. In comparison to previous works, for the first time, we have considered the full magnetic torque, the presence of an outflow, thus the angular momentum transport by vertical motion, and the binary torque. We discuss its specific 3D structure, and how it is affected by tidal effects. We find that the spiral structure evolving in the disk is {\em launched into the outflow}. We propose to call this newly discovered structure a {\em jet spiral wall}. These spiral features follow the same time evolution, with the jet spiral somewhat lagging the disk spiral. We find that the vertical transport of angular momentum has a significant role in the total angular momentum budget also in a binary system. The same holds for the magnetic torque, however, the contribution from the $φ$-derivative of magnetic pressure and the $B_φB_r$ stresses are small. The gravity torque arising from the time-dependent 3D Roche potential becomes essential, as it constitutes the fundamental cause for all 3D effects appearing in our disk-jet system. Quantitatively, we find that the disk accretion rate in a binary system increases by $20\%$ compared to a disk around a single star. The disk wind mass flux increases by even 50\%. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: 32 pages,24 figures, Accepted for publication in Astrophysical Journal, comments are welcome

Bipolar planetary nebulae from common envelope evolution of binary stars

Authors: P. A. Ondratschek, F. K. Roepke, F. R. N. Schneider, C. Fendt, C. Sand, S. T. Ohlmann, R. Pakmor, V. Springel

Abstract: Asymmetric shapes and evidence for binary central stars suggest a common-envelope origin for many bipolar planetary nebulae. The bipolar components of the nebulae are observed to expand faster than the rest and the more slowly expanding material has been associated with the bulk of the envelope ejected during the common-envelope phase of a stellar binary system. Common-envelope evolution in genera… ▽ More Asymmetric shapes and evidence for binary central stars suggest a common-envelope origin for many bipolar planetary nebulae. The bipolar components of the nebulae are observed to expand faster than the rest and the more slowly expanding material has been associated with the bulk of the envelope ejected during the common-envelope phase of a stellar binary system. Common-envelope evolution in general remains one of the biggest uncertainties in binary star evolution and the origin of the fast outflow has not been explained satisfactorily. We perform three-dimensional magnetohydrodynamic simulations of common-envelope interaction with the moving-mesh code AREPO. Starting from the plunge-in of the companion into the envelope of an asymptotic giant branch star and covering hundreds of orbits of the binary star system, we are able to follow the evolution to complete envelope ejection. We find that magnetic fields are strongly amplified in two consecutive episodes. First, when the companion spirals in the envelope and, second, when it forms a contact binary with the core of the former giant star. In the second episode, a magnetically-driven, high-velocity outflow of gas is launched self-consistently in our simulations. The outflow is bipolar and the gas is additionally collimated by the ejected common envelope. The resulting structure reproduces typical morphologies and velocities observed in young planetary nebulae. We propose that the magnetic driving mechanism is a universal consequence of common envelope interaction responsible for a substantial fraction of observed planetary nebulae. Such a mechanism likely also exists in the common-envelope phase of other binary stars that lead to the formation of Type Ia supernovae, X-ray binaries and gravitational-wave merger events. △ Less

Submitted 17 March, 2022; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: 9 pages, 11 figures; accepted for publication by A&A

Journal ref: A&A 660, L8 (2022)

Jets, disc-winds and oscillations in general relativistic, magnetically driven flows around black hole

Authors: Indu K. Dihingia, Bhargav Vaidya, Christian Fendt

Abstract: Relativistic jets and disc-winds are typically observed in BH-XRBs and AGNs. However, many physical details of jet launching and the driving of disc winds from the underlying accretion disc are still not fully understood. In this study, we further investigate the role of the magnetic field strength and structure in launching jets and disc winds. In particular, we explore the connection between jet… ▽ More Relativistic jets and disc-winds are typically observed in BH-XRBs and AGNs. However, many physical details of jet launching and the driving of disc winds from the underlying accretion disc are still not fully understood. In this study, we further investigate the role of the magnetic field strength and structure in launching jets and disc winds. In particular, we explore the connection between jet, wind, and the accretion disc around the central black hole. We perform axisymmetric GRMHD simulations of the accretion-ejection system using adaptive mesh refinement. Essentially, our simulations are initiated with a thin accretion disc in equilibrium. An extensive parametric study by choosing different combinations of magnetic field strength and initial magnetic field inclination is also performed. Our study finds relativistic jets driven by the Blandford \& Znajek (BZ) mechanism and the disc-wind driven by the Blandford \& Payne (BP) mechanism. We also find that plasmoids are formed due to the reconnection events, and these plasmoids advect with disc-winds. As a result, the tension force due to the poloidal magnetic field is enhanced in the inner part of the accretion disc, resulting in disc truncation and oscillation. These oscillations result in flaring activities in the jet mass flow rates. We find simulation runs with a lower value of the plasma-$β$, and lower inclination angle parameters are more prone to the formation of plasmoids and subsequent inner disc oscillations. Our models provide a possible template to understand spectral state transition phenomena in BH-XRBs. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Comments: 21 pages, 20 figures, accepted for publication in MNRAS

A ring accelerator? Unusual jet dynamics in the IceCube candidate PKS 1502+106

Authors: Silke Britzen, Michal Zajaček, Luka Č. Popović, Christian Fendt, Andrea Tramacere, Ilya N. Pashchenko, Frédéric Jaron, Radim Pánis, Leonid Petrov, Margo F. Aller, Hugh D. Aller

Abstract: On 2019/07/30.86853 UT, IceCube detected a high-energy astrophysical neutrino candidate. The Flat Spectrum Radio Quasar PKS 1502+106 is located within the 50 percent uncertainty region of the event. Our analysis of 15 GHz Very Long Baseline Array (VLBA) and astrometric 8 GHz VLBA data, in a time span prior and after the IceCube event, reveals evidence for a radio ring structure which develops with… ▽ More On 2019/07/30.86853 UT, IceCube detected a high-energy astrophysical neutrino candidate. The Flat Spectrum Radio Quasar PKS 1502+106 is located within the 50 percent uncertainty region of the event. Our analysis of 15 GHz Very Long Baseline Array (VLBA) and astrometric 8 GHz VLBA data, in a time span prior and after the IceCube event, reveals evidence for a radio ring structure which develops with time. Several arc-structures evolve perpendicular to the jet ridge line. We find evidence for precession of a curved jet based on kinematic modelling and a periodicity analysis. An outflowing broad line region (BLR) based on the C IV line emission (Sloan Digital Sky Survey, SDSS) is found. We attribute the atypical ring to an interaction of the precessing jet with the outflowing material. We discuss our findings in the context of a spine-sheath scenario where the ring reveals the sheath and its interaction with the surroundings (narrow line region, NLR, clouds). We find that the radio emission is correlated with the $γ$-ray emission, with radio lagging the $γ$-rays. Based on the $γ$-ray variability timescale, we constrain the $γ$-ray emission zone to the BLR (30-200 $r_{\rm g}$) and within the jet launching region. We discuss that the outflowing BLR provides the external radiation field for $γ$-ray production via external Compton scattering. The neutrino is most likely produced by proton-proton interaction in the blazar zone (beyond the BLR), enabled by episodic encounters of the jet with dense clouds, i.e. some molecular cloud in the NLR. △ Less

Submitted 27 February, 2021; originally announced March 2021.

Comments: 35 pages, 33 figures, 3 tables; accepted by the MNRAS Main Journal

Relativistic outflows from a GRMHD mean-field disk dynamo

Authors: Christos Vourellis, Christian Fendt

Abstract: We present simulations of thin accretion disks around black holes investigating a mean-field disk dynamo in our resistive GRMHD code (Vourellis et al. 2019) that is able to produce a large scale magnetic flux. We consider a weak seed field in an initially thin disk, a background (turbulent) magnetic diffusivity and the dynamo action. A standard quenching mechanism is applied to mitigate the otherw… ▽ More We present simulations of thin accretion disks around black holes investigating a mean-field disk dynamo in our resistive GRMHD code (Vourellis et al. 2019) that is able to produce a large scale magnetic flux. We consider a weak seed field in an initially thin disk, a background (turbulent) magnetic diffusivity and the dynamo action. A standard quenching mechanism is applied to mitigate the otherwise exponential increase of the magnetic field. Comparison simulations of an initial Fishbone-Moncrief torus suggest that reconnection may provide another quenching mechanism. The dynamo-generated magnetic flux expands from the disk interior into the disk corona, becomes advected by disk accretion, and fills the axial region of the domain. The dynamo leads to an initially rapid increase in magnetic energy and flux, while for later evolutionary stages the growth stabilizes. Accretion towards the black hole depends strongly on the magnetic field structure that develops. The radial field component supports extraction of angular momentum and thus accretion. It also sets the conditions for launching a disk wind, initially from inner disk area. When a strong field has engulfed the disk, strong winds are launched that are predominantly driven by the pressure gradient of the toroidal field. For rotating black holes we identify a Poynting flux-dominated jet, driven by the Blandford-Znajek mechanism. This axial Poynting flux is advected from the disk and therefore accumulates at the expense of the flux carried by the disk wind, that is itself regenerated by the disk dynamo. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: 26 pages, 16 Figures, accepted for publication in ApJ

General relativistic radiation transport: Implications for VLBI/EHT observations of AGN discs, winds and jets

Authors: Bidisha Bandyopadhyay, Christian Fendt, Dominik R. G. Schleicher, Christos Vourellis

Abstract: In 2019, the Event Horizon Telescope Collaboration (EHTC) has published the first image of a supermassive black hole (SMBH) obtained via the Very Large Baseline Interferometry (VLBI) technique. In the future, it is expected that additional and more sensitive VLBI observations will be pursued for other nearby Active Galactic Nuclei (AGN), and it is therefore important to understand which possible f… ▽ More In 2019, the Event Horizon Telescope Collaboration (EHTC) has published the first image of a supermassive black hole (SMBH) obtained via the Very Large Baseline Interferometry (VLBI) technique. In the future, it is expected that additional and more sensitive VLBI observations will be pursued for other nearby Active Galactic Nuclei (AGN), and it is therefore important to understand which possible features can be expected in such images. In this paper, we post-process General Relativistic Magneto-Hydrodynamical (GR-MHD) simulations which include resistivity, thus providing a self-consistent jet formation model, including resistive mass loading of a wind launched from a disc in Keplerian rotation. The ray-tracing is done using the General Relativistic Ray-Tracing code GRTRANS assuming synchrotron emission. We study the appearance of the black hole environment including the accretion disc, winds and jets under a large range of condition, varying black hole mass, accretion rate, spin, inclination angle, disc parameters and observed frequency. When we adopt M87-like parameters, we show that we can reproduce a ring-like feature (similar as observed by the EHT) for some of our simulations. The latter suggests that such Keplerian disc models thus could be consistent with the observed results. Depending on their masses, accretion rates, spin and the sensitivity of the observation, we note that other SMBHs may show additional features like winds and jets in the observations. △ Less

Submitted 6 September, 2021; v1 submitted 21 February, 2021; originally announced February 2021.

Comments: 22 pages, 19 figures, Version matches with the accepted version for publication to MNRAS

An astronomical institute's perspective on meeting the challenges of the climate crisis

Authors: Knud Jahnke, Christian Fendt, Morgan Fouesneau, Iskren Georgiev, Tom Herbst, Melanie Kaasinen, Diana Kossakowski, Jan Rybizki, Martin Schlecker, Gregor Seidel, Thomas Henning, Laura Kreidberg, Hans-Walter Rix

Abstract: Analysing greenhouse gas emissions of an astronomical institute is a first step in reducing its environmental impact. Here, we break down the emissions of the Max Planck Institute for Astronomy in Heidelberg and propose measures for reductions. Analysing greenhouse gas emissions of an astronomical institute is a first step in reducing its environmental impact. Here, we break down the emissions of the Max Planck Institute for Astronomy in Heidelberg and propose measures for reductions. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: Originated at EAS 2020 conference, sustainability session by https://astronomersforplanet.earth - published in Nature Astronomy, September 2020

Journal ref: Nature Astronomy, Volume 4, p. 812-815 (2020)

MHD accretion-ejection: jets launched by a non-isotropic accretion disk dynamo. II. A dynamo tensor defined by the disk Coriolis number

Authors: Giancarlo Mattia, Christian Fendt

Abstract: Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. We continue our work from Paper I (Mattia & Fendt 2020a), considering a non-scalar accretion disk mean-field $α^2Ω$-dynamo in the context of large scal… ▽ More Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. We continue our work from Paper I (Mattia & Fendt 2020a), considering a non-scalar accretion disk mean-field $α^2Ω$-dynamo in the context of large scale disk-jet simulations. We now investigate a disk dynamo that follows analytical solutions of mean-field dynamo theory, essentially based only on a single parameter, the Coriolis number. We thereby confirm the anisotropy of the dynamo tensor acting in accretion disks, allowing to relate both the resistivity and mean-field dynamo to the disk turbulence. Our new model recovers previous simulations applying a purely radial initial field, while allowing for a more stable evolution for seed fields with a vertical component. We also present correlations between the strength of the disk dynamo coefficients and the dynamical parameters of the jet that is launched, and discuss their implication for observed jet quantities. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: 15 pages, 10 figures, Accepted for publication in ApJ

Journal ref: 2020 ApJ 900 60

MHD accretion-ejection: jets launched by a non-isotropic accretion disk dynamo. I. Validation and application of selected dynamo tensorial components

Authors: Giancarlo Mattia, Christian Fendt

Abstract: Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the jet launching magnetic field is one of the open questions for modeling the accretion-ejection process. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. Applying the PLUTO code, we pr… ▽ More Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the jet launching magnetic field is one of the open questions for modeling the accretion-ejection process. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. Applying the PLUTO code, we present the first resistive MHD simulations of jet launching including a non-scalar accretion disk mean-field $α^2Ω$-dynamo in the context of large scale disk-jet simulations. Essentially, we find the $α_φ$-dynamo component determining the amplification of the poloidal magnetic field, which is strictly related to the disk magnetization (and, as a consequence, to the jet speed, mass and collimation), while the $α_R$ and $α_θ$-dynamo components trigger the formation of multiple, anti-aligned magnetic loops in the disk, with strong consequences on the stability and dynamics of the disk-jet system. In particular, such loops trigger the formation of dynamo inefficient zones, which are characterized by a weak magnetic field, and therefore a lower value of the magnetic diffusivity. The jet mass, speed and collimation are strongly affected by the formation of the dynamo inefficient zones. Moreover, the $θ$-component of the $α$-dynamo plays a key role when interacting with a non-radial component of the seed magnetic field. We also present correlations between the strength of the disk toy dynamo coefficients and the dynamical parameters of the jet that is launched. △ Less

Submitted 30 July, 2020; originally announced July 2020.

Comments: 20 pages, 9 figures, Accepted for publication in ApJ

Journal ref: 2020 ApJ 900 59

Birth of convective low-mass to high-mass second Larson cores

Authors: Asmita Bhandare, Rolf Kuiper, Thomas Henning, Christian Fendt, Mario Flock, Gabriel-Dominique Marleau

Abstract: Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the depen… ▽ More Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the dependence of its properties on the initial cloud core mass. We used the PLUTO code to perform high resolution, 1D and 2D RHD collapse simulations. We include self-gravity and use a grey FLD approximation for the radiative transfer. Additionally, we use for the gas EOS density- and temperature-dependent thermodynamic quantities to account for the effects such as dissociation, ionisation, and molecular vibrations and rotations. Properties of the second core are investigated using 1D studies spanning a wide range of initial cloud core masses from 0.5 to 100 $M_{\odot}$. Furthermore, we expand to 2D collapse simulations for a few cases of 1, 5, 10, and 20 $M_{\odot}$. We follow the evolution of the second core for $\geq$ 100 years after its formation, for each of these non-rotating cases. Our results indicate a dependence of several second core properties on the initial cloud core mass. For the first time, due to an unprecedented resolution, our 2D non-rotating collapse studies indicate that convection is generated in the outer layers of the second core, which is formed due to the gravitational collapse of a 1 $M_{\odot}$ cloud core. Additionally, we find large-scale oscillations of the second accretion shock front triggered by the standing accretion shock instability, which has not been seen before in early evolutionary stages of stars. We predict that the physics within the second core would not be significantly influenced by the effects of magnetic fields or an initial cloud rotation. △ Less

Submitted 26 April, 2020; originally announced April 2020.

Comments: 22 pages, 20 figures, and 2 tables. Accepted for publication in A&A. Collapse movie: https://keeper.mpdl.mpg.de/f/f04abdeabdf3472fb56d/

Journal ref: A&A 638, A86 (2020)

Global axisymmetric simulations of photoevaporation and magnetically driven protoplanetary disk winds

Authors: Peter J. Rodenkirch, Hubert Klahr, Christian Fendt, Cornelis P. Dullemond

Abstract: Photoevaporation and magnetically driven winds are two independent mechanisms to remove mass from protoplanetary disks. In addition to accretion, the effect of these two principles acting concurrently could be significant and the transition between those two has not been extensively studied and quantified in the literature yet. In order to contribute to the understanding of disk winds, we present… ▽ More Photoevaporation and magnetically driven winds are two independent mechanisms to remove mass from protoplanetary disks. In addition to accretion, the effect of these two principles acting concurrently could be significant and the transition between those two has not been extensively studied and quantified in the literature yet. In order to contribute to the understanding of disk winds, we present the phenomena emerging in the framework of two-dimensional axisymmetric, non-ideal magnetohydrodynamic simulations including EUV-/ X-ray driven photoevaporation. Of particular interest are the examination of the transition region between photoevaporation and magnetically driven wind, the possibility of emerging magneto-centrifugal wind effects, as well as the morphology of the wind itself depending on the strength of the magnetic field. We use the PLUTO code in a 2.5D axisymmetric configuration with additional treatment of EUV-/ X-ray heating and dynamic ohmic diffusion based on a semi-analytical chemical model. We identify the transition between both outflow types to occur for values of the initial plasma beta $β\geq 10^7$, while magnetically driven winds generally outperform photoevaporation for stronger fields. In our simulations we observe irregular and asymmetric outflows for stronger magnetic fields. In the weak field regime the photoevaporation rates are slightly lowered by perturbations of the gas density in the inner regions of the disk. Overall, our results predict a wind with a lever arm smaller than 1.5, consistent with a hot magneto-thermal wind. Stronger accretion flows are present for values of $β< 10^7$. △ Less

Submitted 15 January, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: Published in A&A 633, A21 (2020)

Journal ref: A&A 633, A21 (2020)

Approaching the black hole by numerical simulations

Authors: Christian Fendt

Abstract: Black holes represent extreme conditions of physical laws. Being predicted about a century ago, they are now accepted as astrophysical reality by most of the scientific community. Only recently more direct evidence of their existence has been found - the detection of gravitational waves from black hole mergers and of the shadow of a supermassive black hole in the center of a galaxy. Astrophysical… ▽ More Black holes represent extreme conditions of physical laws. Being predicted about a century ago, they are now accepted as astrophysical reality by most of the scientific community. Only recently more direct evidence of their existence has been found - the detection of gravitational waves from black hole mergers and of the shadow of a supermassive black hole in the center of a galaxy. Astrophysical black holes are typically embedded in an active environment which is affected by the strong gravity. When the environmental material emits radiation, this radiation may carry imprints of the black hole that is hosting the radiation source. In order to understand the physical processes that take place in the close neighbourhood of astrophysical black holes, numerical methods and simulations play an essential role. This is simply because the dynamical evolution and the radiative interaction are far too complex in order to allow for an analytic solution of the physical equations. A huge progress has been made over the last decade(s) in the numerical code development, as well as in the computer power that is needed to run these codes. This review tries to summarize the basic questions and methods that are involved in the undertaking of investigating the astrophysics of black holes by numerical means. It is intended for a non-expert audience interested in an overview over this broad field. The review comes along without equations and thus without a detailed expert discussion of the underlying physical processes or numerical specifics. Instead, it intends to illustrate the richness of the field and to motivate for further reading. The review puts some emphasis on magnetohydrodynamic simulations, but also touches radiation transfer and merger simulations, in particular pointing out differences in these approaches. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 18 pages, 6 figures, 130 references

Journal ref: 2019, Universe, vol. 5, issue 5, p. 99

GR-MHD disk winds and jets from black holes and resistive accretion disks

Authors: Christos Vourellis, Christian Fendt, Qian Qian, Scott Noble

Abstract: We perform GR-MHD simulations of outflow launching from thin accretion disks. As in the non-relativistic case, resistivity is essential for the mass loading of the disk wind. We implemented resistivity in the ideal GR-MHD code HARM3D, extending previous works (Qian et al. 2017, 2018) for larger physical grids, higher spatial resolution, and longer simulation time. We consider an initially thin, re… ▽ More We perform GR-MHD simulations of outflow launching from thin accretion disks. As in the non-relativistic case, resistivity is essential for the mass loading of the disk wind. We implemented resistivity in the ideal GR-MHD code HARM3D, extending previous works (Qian et al. 2017, 2018) for larger physical grids, higher spatial resolution, and longer simulation time. We consider an initially thin, resistive disk orbiting the black hole, threaded by a large-scale magnetic flux. As the system evolves, outflows are launched from the black hole magnetosphere and the disk surface. We mainly focus on disk outflows, investigating their MHD structure and energy output in comparison with the Poynting-dominated black hole jet. The disk wind encloses two components -- a fast component dominated by the toroidal magnetic field and a slower component dominated by the poloidal field. The disk wind transitions from sub to super-Alfvénic speed, reaching velocities $\simeq 0.1c$. We provide parameter studies varying spin parameter and resistivity level, and measure the respective mass and energy fluxes. A higher spin strengthens the $B_φ$-dominated disk wind along the inner jet. We disentangle a critical resistivity level that leads to a maximum matter and energy output for both, resulting from the interplay between re-connection and diffusion, which in combination govern the magnetic flux and the mass loading. For counter-rotating black holes the outflow structure shows a magnetic field reversal. We estimate the opacity of the inner-most accretion stream and the outflow structure around it. This stream may be critically opaque for a lensed signal, while the axial jet funnel remains optically thin. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 27 pages, 31 figures, accepted for publication in ApJ

First Core Properties: From Low- to High-mass Star Formation

Authors: Asmita Bhandare, Rolf Kuiper, Thomas Henning, Christian Fendt, Gabriel-Dominique Marleau, Anders Kölligan

Abstract: In this study, the main goal is to understand the molecular cloud core collapse through the stages of first and second hydrostatic core formation. We investigate the properties of Larson's first and second cores following the evolution of the molecular cloud core until formation of Larson's cores. We expand these collapse studies for the first time to span a wide range of initial cloud masses from… ▽ More In this study, the main goal is to understand the molecular cloud core collapse through the stages of first and second hydrostatic core formation. We investigate the properties of Larson's first and second cores following the evolution of the molecular cloud core until formation of Larson's cores. We expand these collapse studies for the first time to span a wide range of initial cloud masses from 0.5 to 100 Msun. Understanding the complexity of the numerous physical processes involved in the very early stages of star formation requires detailed thermodynamical modeling in terms of radiation transport and phase transitions. For this we use a realistic gas equation of state via a density and temperature-dependent adiabatic index and mean molecular weight to model the phase transitions. We use a gray treatment of radiative transfer coupled with hydrodynamics to simulate Larson's collapse in spherical symmetry. We reveal a dependence of a variety of first core properties on the initial cloud mass. The first core radius and mass increase from the low-mass to the intermediate-mass regime and decrease from the intermediate-mass to the high-mass regime. The lifetime of first cores strongly decreases towards the intermediate- and high-mass regime. Our studies show the presence of a transition region in the intermediate-mass regime. Low-mass protostars tend to evolve through two distinct stages of formation which are related to the first and second hydrostatic cores. In contrast, in the high-mass star formation regime, the collapsing cloud cores rapidly evolve through the first collapse phase and essentially immediately form Larson's second cores. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: Accepted for publication in A&A. 20 pages, 11 figures and 3 tables

Journal ref: A&A 618, A95 (2018)

Long-term simulation of MHD jet launching from an orbiting star-disk system

Authors: Somayeh Sheikhnezami, Christian Fendt

Abstract: We present fully three-dimensional magnetohydrodynamic jet launching simulations of a jet source orbiting in a binary system. We consider a time-dependent binary gravitational potential, thus all tidal forces that are experienced in the non-inertial frame of the jet-launching primary. We investigate systems with different binary separation, different mass ratio, and different inclination between t… ▽ More We present fully three-dimensional magnetohydrodynamic jet launching simulations of a jet source orbiting in a binary system. We consider a time-dependent binary gravitational potential, thus all tidal forces that are experienced in the non-inertial frame of the jet-launching primary. We investigate systems with different binary separation, different mass ratio, and different inclination between the disk plane and the orbital plane. The simulations run over a substantial fraction of the binary orbital period. All simulations show similar local and global non-axisymmetric effects such as local instabilities in the disk and the jet, or global features such as disk spiral arms and warps, or a global re-alignment of the inflow-outflow structure. The disk accretion rate is higher than for axisymmetric simulations, most probably due to the enhanced angular momentum transport by spiral waves. The disk outflow leaves the Roche lobe of the primary and becomes disturbed by tidal effects. While a disk-orbit inclination of $10\degr$ still allows for a persistent outflow, an inclination of $30\degr$ does not, suggesting a critical angle in between. For moderate inclination we find indication for jet precession such that the jet axis starts to follow a circular pattern with an opening cone of $\simeq 8\degr$. Simulations with different mass ratio indicate a change of time scales for the tidal forces to affect the disk-jet system. A large mass ratio (a massive secondary) leads to stronger spiral arms, a higher (average) accretion and a more pronounced jet-counter jet asymmetry. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: 22 pages, 19 Figures, Accepted for publication in APJ

Jet launching in resistive GR-MHD black hole - accretion disk systems

Authors: Qian Qian, Christian Fendt, Christos Vourellis

Abstract: We investigate the launching mechanism of relativistic jets from black hole sources, in particular the strong winds from the surrounding accretion disk. Numerical investigations of the disk wind launching - the simulation of the accretion-ejection transition - have so far almost only been done for non-relativistic systems. From these simulations we know that resistivity, or magnetic diffusivity, p… ▽ More We investigate the launching mechanism of relativistic jets from black hole sources, in particular the strong winds from the surrounding accretion disk. Numerical investigations of the disk wind launching - the simulation of the accretion-ejection transition - have so far almost only been done for non-relativistic systems. From these simulations we know that resistivity, or magnetic diffusivity, plays an important role for the launching process. Here, we extend this treatment to general relativistic magnetohydrodynamics (GR-MHD) applying the resistive GR-MHD code rHARM. Our model setup considers a thin accretion disk threaded by a large-scale open magnetic field. We run a series of simulations with different Kerr parameter, field strength and diffusivity level. Indeed we find strong disk winds with, however, mildly relativistic speed, the latter most probably due to our limited computational domain. Further, we find that magnetic diffusivity lowers the efficiency of accretion and ejection, as it weakens the efficiency of the magnetic lever arm of the disk wind. As major driving force of the disk wind we disentangle the toroidal magnetic field pressure gradient, however,magneto-centrifugal driving may also contribute. Black hole rotation in our simulations suppresses the accretion rate due to an enhanced toroidal magnetic field pressure that seems to be induced by frame-dragging. Comparing the energy fluxes from the Blandford-Znajek-driven central spine and the surrounding disk wind, we find that the total electromagnetic energy flux is dominated by the total matter energy flux of the disk wind (by a factor 20). The kinetic energy flux of the matter outflow is comparatively small and comparable to the Blandford-Znajek electromagnetic energy flux. △ Less

Submitted 25 April, 2018; originally announced April 2018.

Comments: 21 pages, many figures, ApJ accepted

rHARM: Accretion and Ejection in Resistive GR-MHD

Authors: Qian Qian, Christian Fendt, Scott Noble, Matteo Bugli

Abstract: Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity, respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests and preliminary applications. The test simulations of the new code… ▽ More Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity, respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests and preliminary applications. The test simulations of the new code rHARM are compared with an analytic solution of the diffusion equation and a classical shock tube problem. We have further investigated the evolution of the magneto-rotational instability (MRI) in tori around black holes for a range of magnetic diffusivities. We find indication for a critical magnetic diffusivity (for our setup) beyond which no MRI develops in the linear regime and for which accretion of torus material to the black hole is delayed. Preliminary simulations of magnetically diffusive thin accretion disks around Schwarzschild black holes that are threaded by a large-scale poloidal magnetic field show the launching of disk winds with mass fluxes of about 50% of the accretion rate. The disk magnetic diffusivity allows for efficient disk accretion that replenishes the mass reservoir of the inner disk area and thus allows for long-term simulations of wind launching for more than 5000 time units. △ Less

Submitted 14 October, 2016; originally announced October 2016.

Comments: 21 pages, 43 figures, accepted by ApJ

An extensive numerical survey of the correlation between outflow dynamics and accretion disk magnetization

Authors: Deniss Stepanovs, Christian Fendt

Abstract: We investigate the accretion-ejection process of jets from magnetized accretion disks. We apply a novel approach to the jet-launching problem in order to obtain correlations between the physical properties of the jet and the underlying disk. We extend and confirm the previous works of \citet{2009MNRAS.400..820T} and \citet{2010A&A...512A..82M} by scanning a large parameter range for the disk magne… ▽ More We investigate the accretion-ejection process of jets from magnetized accretion disks. We apply a novel approach to the jet-launching problem in order to obtain correlations between the physical properties of the jet and the underlying disk. We extend and confirm the previous works of \citet{2009MNRAS.400..820T} and \citet{2010A&A...512A..82M} by scanning a large parameter range for the disk magnetization, $μ_{\rm D} = 10^{-3.5} ... 10^{-0.7}$. We disentangle the disk magnetization at the foot point of the outflow as the main parameter that governs the properties of the outflow. We show how the four jet integrals known from steady-state MHD are correlated to the disk magnetization at the jet foot point. This agrees with the usual findings of the steady-state theory, however, here we obtain these correlations from time-dependent simulations that include the dynamical evolution of the disk in the treatment. In particular, we obtain robust correlations between the local disk magnetization and (i)the outflow velocity, (ii) the jet mass loading, (iii) jet angular momentum, and (iv) the local mass accretion rate. Essentially we find that strongly magnetized disks launch more energetic and faster jets, and, due to a larger Alfvén lever arm, these jets extract more angular momentum from the underlying disk. These kind of disk-jet systems have, however, a smaller mass loading parameter and a lower mass ejection-to-accretion ratio. The jets are launched at the disk surface where the magnetization is $μ(r,z) \simeq 0.1$. The magnetization rapidly increases vertically providing the energy reservoir for subsequent jet acceleration. We find indication for a critical disk magnetization $μ_{\rm D} \simeq 0.01$ that separates the regimes of magnetocentrifugally-driven and magnetic pressure-driven jets. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: accepted by ApJ

Wobbling and precessing jets from warped disks in binary systems

Authors: Somayeh Sheikhnezami, Christian Fendt

Abstract: We present results of the first ever three-dimensional (3D) magnetohydrodynamic (MHD) simulations of the accretion-ejection structure. We investigate the 3D evolution of jets launched symmetrically from single stars but also jets from warped disks in binary systems. We have applied various model setups and tested them by simulating a stable and bipolar symmetric 3D structure from a single star-dis… ▽ More We present results of the first ever three-dimensional (3D) magnetohydrodynamic (MHD) simulations of the accretion-ejection structure. We investigate the 3D evolution of jets launched symmetrically from single stars but also jets from warped disks in binary systems. We have applied various model setups and tested them by simulating a stable and bipolar symmetric 3D structure from a single star-disk-jet system. Our reference simulation maintains a good axial symmetry and also a bipolar symmetry for more than 600 rotations of the inner disk confirming the quality of our model setup. We have then implemented a 3D gravitational potential (Roche potential) due to a companion star and run a variety of simulations with different binary separations and mass ratios. These simulations show typical 3D deviations from axial symmetry, such as jet inclination outside the Roche lobe or spiral arms forming in the accretion disk. In order to find indication for precession effects, we have also run an exemplary parameter setup, essentially governed by a small binary separation of only $\simeq{200}$ inner disk radii. This simulation shows strong indication that we observe the onset of a jet precession caused by the wobbling of the jet launching disk. We estimate the opening angle of the precession cone defined by the lateral motion of the the jet axis of about 4 degree after about 5000 dynamical time steps. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 21 pages, 15 figures, accepted for publication in APJ, comments are welcome

Disks and Jets - Gravity, Rotation and Magnetic Fields

Authors: John F. Hawley, Christian Fendt, Martin Hardcastle, Elena Nokhrina, Alexander Tchekhovskoy

Abstract: Magnetic fields are fundamental to the dynamics of both accretion disks and the jets that they often drive. We review the basic physics of these phenomena, the past and current efforts to model them numerically with an emphasis on the jet-disk connection, and the observational constraints on the role of magnetic fields in the jets of active galaxies on all scales. Magnetic fields are fundamental to the dynamics of both accretion disks and the jets that they often drive. We review the basic physics of these phenomena, the past and current efforts to model them numerically with an emphasis on the jet-disk connection, and the observational constraints on the role of magnetic fields in the jets of active galaxies on all scales. △ Less

Submitted 11 August, 2015; originally announced August 2015.

Comments: Invited review for Space Science Reviews, 34 pages

Modelling MHD accretion-ejection - episodic ejections of jets triggered by a mean-field disk dynamo

Authors: Deniss Stepanovs, Christian Fendt, Somayeh Sheikhnezami

Abstract: We present MHD simulations exploring the launching, acceleration and collimation of jets and disk winds. The evolution of the disk structure is consistently taken into account. Extending our earlier studies, we now consider the self-generation of the magnetic field by an $α^2Ω$ mean-field dynamo. The disk magnetization remains on a rather low level, that helps to evolve the simulations for… ▽ More We present MHD simulations exploring the launching, acceleration and collimation of jets and disk winds. The evolution of the disk structure is consistently taken into account. Extending our earlier studies, we now consider the self-generation of the magnetic field by an $α^2Ω$ mean-field dynamo. The disk magnetization remains on a rather low level, that helps to evolve the simulations for $T > 10,000$ dynamical time steps on a domain extending 1500 inner disk radii. We find a magnetic field of the inner disk similar to the commonly found open field structure, favoring magneto-centrifugal launching. The outer disk field is highly inclined and predominantly radial. Here, differential rotation induces a strong toroidal component that plays a key role in outflow launching. These outflows from the outer disk are slower, denser, and less collimated. If the dynamo action is not quenched, magnetic flux is continuously generated, diffuses outward the disk, and fills the entire disk. We have invented a toy model triggering a time-dependent mean-field dynamo. The duty cycles of this dynamo lead to episodic ejections on similar timescales. When the dynamo is suppressed as the magnetization falls below a critical value, the generation of the outflows and also accretion is inhibited. The general result is that we can steer episodic ejection and {\em large-scale jet knots} by a {\em disk-intrinsic dynamo} that is time-dependent and regenerates the jet-launching magnetic field. △ Less

Submitted 24 September, 2014; originally announced September 2014.

Comments: 15 pages, 12 figures; accepted in ApJ; comments are welcome

Modelling MHD accretion-ejection - from the launching area to propagation scales

Authors: Deniss Stepanovs, Christian Fendt

Abstract: We present results of axisymmetric magnetohydrodynamic (MHD) simulations investigating the launching of jets and outflows from a magnetically diffusive accretion disk. The time evolution of the disk structure is self-consistently taken into account. In contrast to previous works we have applied {\em spherical} coordinates for the numerical grid, implying substantial benefits concerning the numeric… ▽ More We present results of axisymmetric magnetohydrodynamic (MHD) simulations investigating the launching of jets and outflows from a magnetically diffusive accretion disk. The time evolution of the disk structure is self-consistently taken into account. In contrast to previous works we have applied {\em spherical} coordinates for the numerical grid, implying substantial benefits concerning the numerical resolution and the stability of the simulation. Thanks to the new setup we were able to run simulations for more than 150,000 dynamical times on a domain extending 1500 inner disk radii with a resolution of up to 24 cells per disk height in the inner disk. Depending on the disk magnetization, jet launching occurs in two different but complementary regimes - jets driven predominantly by centrifugal or magnetic forces. These regimes differ in the ejection efficiency concerning mass, energy and angular momentum. We show that it is the {\em actual} disk magnetization and not so much the initial magnetization which describes the disk-jet evolution best. Considering the actual disk magnetization we also find that simulations starting with different initial magnetization evolve in a similar - typical - way as due to advection and diffusion the magnetic flux in the disk evolves in time. Exploring a new, modified diffusivity model we confirm the self-similar structure of the global jet-launching disk, obtaining power laws for the radial profiles of the disk physical variables such as density, magnetic field strength, or accretion velocity. △ Less

Submitted 11 July, 2014; originally announced July 2014.

Comments: 22 pages, 25 figures; submitted to ApJ, revised; suggestions welcome

Formation and collimation of relativistic MHD jets - simulations and radio maps

Authors: Christian Fendt, Oliver Porth, Somayeh Sheikhnezami

Abstract: We present results of magnetohydrodynamic (MHD) simulations of jet formation and propagation, discussing a variety of astrophysical setups. In the first approach we consider simulations of relativistic MHD jet formation, considering jets launched from the surface of a Keplerian disk, demonstrating numerically - for the first time - the self-collimating ability of relativistic MHD jets. We obtain L… ▽ More We present results of magnetohydrodynamic (MHD) simulations of jet formation and propagation, discussing a variety of astrophysical setups. In the first approach we consider simulations of relativistic MHD jet formation, considering jets launched from the surface of a Keplerian disk, demonstrating numerically - for the first time - the self-collimating ability of relativistic MHD jets. We obtain Lorentz factors up to about 10 while acquiring a high degree of collimation of about 1 degree. We then present synchrotron maps calculated from the intrinsic jet structure derived from the MHD jet formation simulation. We finally present (non-relativistic) MHD simulations of jet lauching, treating the transition between accretion and ejection. These setups include a physical magnetic diffusivity which is essential for loading the accretion material onto the outflow. We find relatively high mass fluxes in the outflow, of the order of 20-40 % of the accretion rate. △ Less

Submitted 20 December, 2013; originally announced December 2013.

Comments: 7 pages, 4 figures, conference proceedings HEPRO 4

Bipolar jets launched from accretion disks. II. Formation of symmetric and asymmetric jets and counter jets

Authors: Christian Fendt, Somayeh Sheikhnezami

Abstract: We investigate the jet launching process from accretion disks extending our recent study (paper I) to a truly bipolar setup. We perform axisymmetric MHD simulations of the disk-jet interaction on a computational domain covering both hemispheres, in particular addressing the question of an intrinsically asymmetric origin of jet / counter jet systems. Treating both hemispheres simultaneously, we ove… ▽ More We investigate the jet launching process from accretion disks extending our recent study (paper I) to a truly bipolar setup. We perform axisymmetric MHD simulations of the disk-jet interaction on a computational domain covering both hemispheres, in particular addressing the question of an intrinsically asymmetric origin of jet / counter jet systems. Treating both hemispheres simultaneously, we overcome the equatorial plane symmetry boundary condition used in most previous studies which naturally fosters a symmetric evolution. For the magnetic diffusivity prescription we apply an alpha-parametrisation, considering both, globally models of diffusivity, and local models. We first approve the quality of our numerical setup by generating perfectly symmetric jets, lasting over a 1000s of dynamical time scales. We then disturb the hemispheric symmetry in the disk, and investigate the subsequent evolution of the outflow. The evolution first leads to a substantial disk warping with electric currents intersecting the equatorial plane. We investigate two models, i) a disk with (initially) different thermal scale height in both hemispheres, and ii) a symmetric disk into which a local disturbance is injected in one hemisphere. In both cases the disk asymmetry results in asymmetric outflows with mass fluxes differing by 10-20%. We find up to 30% difference in mass flux between jet and counter jet for this setup, lasting over 1000s of dynamical time scales (i.e. lasting for the whole simulation). In summary, our results suggest that the jet asymmetries in protostellar and extragalactic jets can indeed be generated intrinsically and maintained over long time by disk asymmetries and the standard jet launching mechanism. △ Less

Submitted 6 May, 2013; originally announced May 2013.

Comments: Submitted to ApJ, follow-up of paper I, 2012, comments welcome

Bipolar jets launched from magnetically diffusive accretion disks. I. Ejection efficiency vs field strength and diffusivity

Authors: Somayeh Sheikhnezami, Christian Fendt, Oliver Porth, Bhargav Vaidya, Jamshid Ghanbari

Abstract: We investigate the launching of jets and outflows from magnetically diffusive accretion disks. Using the PLUTO code we solve the time-dependent resistive MHD equations taking into account the disk and jet evolution simultaneously. The main question we address is which kind of disks do launch jets and which kind of disks do not? In particular, we study how the magnitude and distribution of the (tur… ▽ More We investigate the launching of jets and outflows from magnetically diffusive accretion disks. Using the PLUTO code we solve the time-dependent resistive MHD equations taking into account the disk and jet evolution simultaneously. The main question we address is which kind of disks do launch jets and which kind of disks do not? In particular, we study how the magnitude and distribution of the (turbulent) magnetic diffusivity affect mass loading and jet acceleration. We have applied a turbulent magnetic diffusivity based on α-prescription, but have also investigate examples where the scale height of diffusivity is larger than that of the disk gas pressure. We further investigate how the ejection efficiency is governed by the magnetic field strength. Our simulations last for up to 5000 dynamical time scales corresponding to 900 orbital periods of the inner disk. As a general result we observe a continuous and robust outflow launched from the inner part of the disk, expanding into a collimated jet of super fast magneto-sonic speed. For long time scales the disk internal dynamics changes, as due to outflow ejection and disk accretion the disk mass decreases. For magneto-centrifugally driven jets we find that for i) less diffusive disks, ii) a stronger magnetic field, iii) a low poloidal diffusivity, or a iv) lower numerical diffusivity (resolution), the mass loading of the outflow is increased - resulting in more powerful jets with high mass flux. For weak magnetization the (weak) outflow is driven by the magnetic pressure gradient. We further investigate the jet asymptotic velocity and the jet rotational velocity in respect of the different launching scenarios. We find a lower degree of jet collimation than previous studies, most probably due to our revised outflow boundary condition. △ Less

Submitted 25 July, 2012; originally announced July 2012.

Report number: ApJ87055R2

Magnetic field structure of relativistic jets without current sheets

Authors: K. N. Gourgouliatos, Ch. Fendt, E. Clausen-Brown, M. Lyutikov

Abstract: We present an analytical class of equilibrium solutions for the structure of relativistic sheared and rotating magnetized jets that contain no boundary current sheets. We demonstrate the overall dynamical stability of these solutions and, most importantly, a better numerical resistive stability than the commonly employed force-free structures which inevitably require the presence of dissipative su… ▽ More We present an analytical class of equilibrium solutions for the structure of relativistic sheared and rotating magnetized jets that contain no boundary current sheets. We demonstrate the overall dynamical stability of these solutions and, most importantly, a better numerical resistive stability than the commonly employed force-free structures which inevitably require the presence of dissipative surface currents. The jet is volumetrically confined by the external pressure, with no pressure gradient on the surface. We calculate the expected observed properties of such jets. Given the simplicity of these solution we suggest them as useful initial conditions for relativistic jet simulations. △ Less

Submitted 4 October, 2011; originally announced October 2011.

Comments: 13 pages, 13 figures, Accepted by MNRAS

Jet formation from massive young stars: Magnetohydrodynamics versus radiation pressure

Authors: Bhargav Vaidya, Christian Fendt, Henrik Beuther, Oliver Porth

Abstract: Observations indicate that outflows from massive young stars are more collimated during their early evolution compared to later stages. Our paper investigates various physical processes that impacts the outflow dynamics, i.e. its acceleration and collimation. We perform axisymmetric MHD simulations particularly considering the radiation pressure exerted by the star and the disk. We have modified t… ▽ More Observations indicate that outflows from massive young stars are more collimated during their early evolution compared to later stages. Our paper investigates various physical processes that impacts the outflow dynamics, i.e. its acceleration and collimation. We perform axisymmetric MHD simulations particularly considering the radiation pressure exerted by the star and the disk. We have modified the PLUTO code to include radiative forces in the line-driving approximation. We launch the outflow from the innermost disk region (r < 50 AU) by magneto-centrifugal acceleration. In order to disentangle MHD effects from radiative forces, we start the simulation in pure MHD, and later switch on the radiation force. We perform a parameter study considering different stellar masses (thus luminosity), magnetic flux, and line-force strength. For our reference simulation - assuming a 30 Msun star, we find substantial de-collimation of 35 % due to radiation forces. The opening angle increases from 20 deg to 32 deg for stellar masses from 20 Msun to 60 Msun. A small change in the line-force parameter 'alpha' from 0.60 to 0.55 changes the opening angle by ~ 8 deg. We find that it is mainly the stellar radiation which affects the jet dynamics. Unless the disk extends very close to the star, its pressure is too small to have much impact. Essentially, our parameter runs with different stellar mass can be understood as a proxy for the time evolution of the star-outflow system. Thus, we have shown that when the stellar mass (thus luminosity) increases (with age), the outflows become less collimated. △ Less

Submitted 24 August, 2011; originally announced August 2011.

Comments: Accepted for publication in ApJ. 22 pages (emulateapj), 14 figures

Jet rotation driven by MHD shocks in helical magnetic fields

Authors: Christian Fendt

Abstract: In this paper we present a detailed numerical investigation of the hypothesis that a rotation of astrophysical jets can be caused by magnetohydrodynamic shocks in a helical magnetic field. Shock compression of the helical magnetic field results in a toroidal Lorentz force component which will accelerate the jet material in toroidal direction. This process transforms magnetic angular momentum (magn… ▽ More In this paper we present a detailed numerical investigation of the hypothesis that a rotation of astrophysical jets can be caused by magnetohydrodynamic shocks in a helical magnetic field. Shock compression of the helical magnetic field results in a toroidal Lorentz force component which will accelerate the jet material in toroidal direction. This process transforms magnetic angular momentum (magnetic stress) carried along the jet into kinetic angular momentum (rotation). The mechanism proposed here only works in a helical magnetic field configuration. We demonstrate the feasibility of this mechanism by axisymmetric MHD simulations in 1.5D and 2.5D using the PLUTO code. In our setup the jet is injected into the ambient gas with zero kinetic angular momentum (no rotation). Different dynamical parameters for jet propagation are applied such as the jet internal Alfven Mach number and fast magnetosonic Mach number, the density contrast of jet to ambient medium, or the external sonic Mach number of the jet. The mechanism we suggest should work for a variety of jet applications, e.g. protostellar or extragalactic jets, and internal jet shocks (jet knots) or external shocks between the jet and ambient gas (entrainment). For typical parameter values for protostellar jets, the numerically derived rotation feature looks consistent with the observations, i.e. rotational velocities of 0.1-1 percent of the jet bulk velocity. △ Less

Submitted 31 May, 2011; originally announced May 2011.

Comments: 20 pages, 15 figures, accepted by ApJ

Synchrotron radiation of self-collimating relativistic MHD jets

Authors: Oliver Porth, Christian Fendt, Zakaria Meliani, Bhargav Vaidya

Abstract: The goal of this paper is to derive signatures of synchrotron radiation from state-of-the-art simulation models of collimating relativistic magnetohydrodynamic (MHD) jets featuring a large-scale helical magnetic field. We perform axisymmetric special relativistic MHD simulations of the jet acceleration region using the PLUTO code. The computational domain extends from the slow magnetosonic launchi… ▽ More The goal of this paper is to derive signatures of synchrotron radiation from state-of-the-art simulation models of collimating relativistic magnetohydrodynamic (MHD) jets featuring a large-scale helical magnetic field. We perform axisymmetric special relativistic MHD simulations of the jet acceleration region using the PLUTO code. The computational domain extends from the slow magnetosonic launching surface of the disk up to 6000^2 Schwarzschild radii allowing to reach highly relativistic Lorentz factors. The Poynting dominated disk wind develops into a jet with Lorentz factors of 8 and is collimated to 1 degree. In addition to the disk jet, we evolve a thermally driven spine jet, emanating from a hypothetical black hole corona. Solving the linearly polarized synchrotron radiation transport within the jet, we derive VLBI radio and (sub-) mm diagnostics such as core shift, polarization structure, intensity maps, spectra and Faraday rotation measure (RM), directly from the Stokes parameters. We also investigate depolarization and the detectability of a lambda^2-law RM depending on beam resolution and observing frequency. We find non-monotonic intrinsic RM profiles which could be detected at a resolution of 100 Schwarzschild radii. In our collimating jet geometry, the strict bi-modality in polarization direction (as predicted by Pariev et al.) can be circumvented. Due to relativistic aberration, asymmetries in the polarization vectors across the jet can hint to the spin direction of the central engine. △ Less

Submitted 21 May, 2011; originally announced May 2011.

Comments: Submitted to ApJ

Probing the Transition Between the Synchrotron and Inverse-compton Spectral Components of 1ES 1959+650

Authors: E. Bottacini, M. Böttcher, P. Schady, A. Rau, X. -L. Zhang, M. Ajello, C. Fendt, J. Greiner

Abstract: 1ES 1959+650 is one of the most remarkable high-peaked BL Lacertae objects (HBL). In 2002, it exhibited a TeV gamma-ray flare without a similar brightening of the synchrotron component at lower energies. We present the results of a multifrequency campaign, triggered by the INTEGRAL IBIS detection of 1ES 1959+650. Our data range from the optical to hard X-ray energies, thus covering the synchrotron… ▽ More 1ES 1959+650 is one of the most remarkable high-peaked BL Lacertae objects (HBL). In 2002, it exhibited a TeV gamma-ray flare without a similar brightening of the synchrotron component at lower energies. We present the results of a multifrequency campaign, triggered by the INTEGRAL IBIS detection of 1ES 1959+650. Our data range from the optical to hard X-ray energies, thus covering the synchrotron and inverse-Compton components simultaneously. We observed the source with INTEGRAL, the Swift X-Ray Telescope, and the UV-Optical Telescope, and nearly simultaneously with a ground-based optical telescope. The steep spectral component at X-ray energies is most likely due to synchrotron emission, while at soft gamma-ray energies the hard spectral index may be interpreted as the onset of the high-energy component of the blazar spectral energy distribution (SED). This is the first clear measurement of a concave X-ray-soft gamma-ray spectrum for an HBL. The SED can be well modeled with a leptonic synchrotron self-Compton model. When the SED is fitted this model requires a very hard electron spectral index of q ~ 1.85, possibly indicating the relevance of second-order Fermi acceleration. △ Less

Submitted 15 October, 2010; originally announced October 2010.

Comments: 5 pages, 2 postscript figures

Journal ref: APJL 719 (2010) 162-166

Acceleration and collimation of relativistic MHD disk winds

Authors: O. Porth, C. Fendt

Abstract: We perform axisymmetric relativistic magnetohydrodynamic (MHD) simulations to investigate the acceleration and collimation of jets and outflows from disks around compact objects. The fiducial disk surface (respectively a slow disk wind) is prescribed as boundary condition for the outflow. We apply this technique for the first time in the context of relativistic jets. The strength of this approac… ▽ More We perform axisymmetric relativistic magnetohydrodynamic (MHD) simulations to investigate the acceleration and collimation of jets and outflows from disks around compact objects. The fiducial disk surface (respectively a slow disk wind) is prescribed as boundary condition for the outflow. We apply this technique for the first time in the context of relativistic jets. The strength of this approach is that it allows us to run a parameter study in order to investigate how the accretion disk conditions govern the outflow formation. Our simulations using the PLUTO code run for 500 inner disk rotations and on a physical grid size of 100x200 inner disk radii. In general, we obtain collimated beams of mildly relativistic speed and mass-weighted half-opening angles of 3-7 degrees. When we increase the outflow Poynting flux by injecting an additional disk toroidal field into the inlet, Lorentz factors up to 6 are reached. These flows gain super-magnetosonic speed and remain Poynting flux dominated. The light surface of the outflow magnetosphere tends to align vertically - implying three relativistically distinct regimes in the flow - an inner sub-relativistic domain close to the jet axis, a (rather narrow) relativistic jet and a surrounding subrelativistic outflow launched from the outer disk surface - similar to the spine-sheath structure currently discussed for asymptotic jet propagation and stability. The outer subrelativistic disk wind is a promising candidate for the X-ray absorption winds that are observed in many radio-quiet AGN. △ Less

Submitted 14 December, 2009; v1 submitted 16 November, 2009; originally announced November 2009.

Comments: 22 pages, 15 figures; accepted for publication in ApJ; incorporates changes according to referee

Journal ref: Astrophys.J.709:1100-1118,2010

Accretion disks around massive stars: Hydrodynamic structure, stability and dust sublimation

Authors: B. Vaidya, C. Fendt, H. Beuther

Abstract: We investigate the structure of accretion disks around massive protostar applying steady state models of thin disks. The thin disk equations are solved with proper opacities for dust and gas taking into account the huge temperature variation along the disk. We explore a wide parameter range concerning stellar mass, accretion rate, and viscosity parameter α. The most essential finding is a very h… ▽ More We investigate the structure of accretion disks around massive protostar applying steady state models of thin disks. The thin disk equations are solved with proper opacities for dust and gas taking into account the huge temperature variation along the disk. We explore a wide parameter range concerning stellar mass, accretion rate, and viscosity parameter α. The most essential finding is a very high temperature of the inner disk. For e.g. a 10 M_sun protostar and an accretion rate of 10^-4 M_sun/yr, the disk midplane temperature may reach almost 10^5 K. The disk luminosity in this case is about 10^4 L_sun and, thus, potentially higher than that of a massive protostar. We motivate our disk model with similarly hot disks around compact stars. We calculate a dust sublimation radius by turbulent disk self-heating of more than 10AU, a radius, which is 3 times larger than caused by stellar irradiation. We discuss implications of this result on the flashlight effect and the consequences for the radiation pressure of the central star. In difference to disks around low mass protostars our models suggest rather high values for the disk turbulence parameter αclose to unity. However, disk stability to fragmentation due to thermal effects and gravitational instability would require a lower αvalue. For α= 0.1 we find stable disks out to 80AU. Essentially, our model allows to compare the outer disk to some of the observed massive protostellar disk sources, and from that, extrapolate on the disk structure close to the star which is yet impossible to observe. △ Less

Submitted 30 June, 2009; originally announced June 2009.

Comments: 23 pages, 8 figures, accepted by ApJ, For full version see http://www.mpia.de/homes/vaidya/

Journal ref: Astrophys.J.702:567-579,2009

Formation of relativistic MHD jets: stationary state solutions & numerical simulations

Authors: Christian Fendt, Elisabetta Memola

Abstract: We discuss numerical results of relativistic magnetohydrodynamic (MHD) jet formation models. We first review some examples of stationary state solutions treating the collimation and acceleration process of relativistic MHD jets. We provide an a posteriori check for the MHD condition in highly magnetized flows, namely the comparison of particle density to Goldreich-Julian density. Using the jet d… ▽ More We discuss numerical results of relativistic magnetohydrodynamic (MHD) jet formation models. We first review some examples of stationary state solutions treating the collimation and acceleration process of relativistic MHD jets. We provide an a posteriori check for the MHD condition in highly magnetized flows, namely the comparison of particle density to Goldreich-Julian density. Using the jet dynamical parameters calculated from the MHD model we show the rest-frame thermal X-ray spectra of the jet, from which we derive the overall spectrum taking into account a variation of Doppler boosting and Doppler shift of emission lines along the outflow. Finally, we present preliminary results of relativistic MHD simulations of jet formation demonstrating the acceleration of a low velocity (0.01c) disk wind to a collimated high velocity (0.8c). △ Less

Submitted 20 November, 2008; originally announced November 2008.

Comments: 6 pp, 5 figs; Invited talk at High Energy Processes in Relativistic Jets, Dublin, 2007, in press

Journal ref: Int.J.Mod.Phys.D17:1677-1686,2008

Magnetohydrodynamic jets from different magnetic field configurations

Authors: Christian Fendt

Abstract: Using axisymmetric MHD simulations we investigate how the overall jet formation is affected by a variation in the disk magnetic flux profile and/or the existence of a central stellar magnetosphere. Our simulations evolve from an initial, hydrostatic equilibrium state in a force-free magnetic field configuration. We find a unique relation between the collimation degree and the disk wind magnetiza… ▽ More Using axisymmetric MHD simulations we investigate how the overall jet formation is affected by a variation in the disk magnetic flux profile and/or the existence of a central stellar magnetosphere. Our simulations evolve from an initial, hydrostatic equilibrium state in a force-free magnetic field configuration. We find a unique relation between the collimation degree and the disk wind magnetization power law exponent. The collimation degree decreases for steeper disk magnetic field profiles. Highly collimated outflows resulting from a flat profile tend to be unsteady. We further consider a magnetic field superposed of a stellar dipole and a disk field in parallel or anti-parallel alignment. Both stellar and disk wind may evolve in a pair of outflows, however, a reasonably strong disk wind component is essential for jet collimation. Strong flares may lead to a sudden change in mass flux by a factor two. We hypothesize that such flares may eventually trigger jet knots. △ Less

Submitted 20 November, 2008; originally announced November 2008.

Comments: 5 pages, 4 figures; proceedings from conference: Protostellar Jets in Context, held in Rhodes, July 7-12, 2008

Formation of protostellar jets as two-component outflows from star-disk magnetospheres

Authors: Christian Fendt

Abstract: Axisymmetric magnetohydrodynamic (MHD) simulations have been applied to investigate the interrelation of a central stellar magnetosphere and stellar wind with a surrounding magnetized disk outflow and how the overall formation of a large scale jet is affected. The initial magnetic field distribution applied is a superposition of two components - a stellar dipole and a surrounding disk magnetic f… ▽ More Axisymmetric magnetohydrodynamic (MHD) simulations have been applied to investigate the interrelation of a central stellar magnetosphere and stellar wind with a surrounding magnetized disk outflow and how the overall formation of a large scale jet is affected. The initial magnetic field distribution applied is a superposition of two components - a stellar dipole and a surrounding disk magnetic field, in both either parallel or anti-parallel alignment. Correspondingly, the mass outflow is launched as stellar wind plus a disk wind. Our simulations evolve from an initial state in hydrostatic equilibrium and an initially force-free magnetic field configuration. Due to initial differential rotation and induction of a strong toroidal magnetic field the stellar dipolar field inflates and is disrupted on large scale. Stellar and disk wind may evolve in a pair of collimated outflows. The existence of a reasonably strong disk wind component is essential for collimation. A disk jet as known from previous numerical studies will become de-collimated by the stellar wind. In some simulations we observe the generation of strong flares triggering a sudden change in the outflow mass loss rate by a factor of two and also a re-distribution in the radial profile of momentum flux and jet velocity across the jet. We discuss the hypothesis that these flares may trigger internal shocks in the asymptotic jets which are observed as knots. △ Less

Submitted 22 October, 2008; originally announced October 2008.

Comments: 35 pages, 14 x 5 figures; accepted by ApJ

Journal ref: Astrophys.J.692:346-363,2009

Can Protostellar Jets Drive Supersonic Turbulence in Molecular Clouds?

Authors: Robi Banerjee, Ralf S. Klessen, Christian Fendt

Abstract: Jets and outflows from young stellar objects are proposed candidates to drive supersonic turbulence in molecular clouds. Here, we present the results from multi-dimensional jet simulations where we investigate in detail the energy and momentum deposition from jets into their surrounding environment and quantify the character of the excited turbulence with velocity probability density functions.… ▽ More Jets and outflows from young stellar objects are proposed candidates to drive supersonic turbulence in molecular clouds. Here, we present the results from multi-dimensional jet simulations where we investigate in detail the energy and momentum deposition from jets into their surrounding environment and quantify the character of the excited turbulence with velocity probability density functions. Our study include jet--clump interaction, transient jets, and magnetised jets. We find that collimated supersonic jets do not excite supersonic motions far from the vicinity of the jet. Supersonic fluctuations are damped quickly and do not spread into the parent cloud. Instead subsonic, non-compressional modes occupy most of the excited volume. This is a generic feature which can not be fully circumvented by overdense jets or magnetic fields. Nevertheless, jets are able to leave strong imprints in their cloud structure and can disrupt dense clumps. Our results question the ability of collimated jets to sustain supersonic turbulence in molecular clouds. △ Less

Submitted 25 June, 2007; originally announced June 2007.

Comments: 33 pages, 18 figures, accepted by ApJ, version with high resolution figures at: http://www.ita.uni-heidelberg.de/~banerjee/publications/jet_paper.pdf

Magnetic interaction of jets and molecular clouds in NGC 4258

Authors: Marita Krause, Christian Fendt, Nikolaus Neininger

Abstract: NGC 4258 is a well known spiral galaxy with a peculiar large scale jet flow detected in the radio and in H-alpha. Due to the special geometry of the galaxy, the jets emerge from the nuclear region through the galactic disk. Also the distribution of molecular gas looks different from that in other spiral galaxies: [12]CO(1-0)emission has only been detected in the center and along the jets and onl… ▽ More NGC 4258 is a well known spiral galaxy with a peculiar large scale jet flow detected in the radio and in H-alpha. Due to the special geometry of the galaxy, the jets emerge from the nuclear region through the galactic disk. Also the distribution of molecular gas looks different from that in other spiral galaxies: [12]CO(1-0)emission has only been detected in the center and along the jets and only up to distances of about 50 arcsec (1.8 kpc) from the nucleus. The reason for the CO concentration along the inner jets in NGC 4258 was not understood and is the motivation for the observations presented here. Using the IRAM interferometer, we mapped the [12]CO(1-0) emission of the central part of NGC 4258 along the nuclear jet direction in the inner 3 kpc. We detected two parallel CO ridges along a position angle of -25 degr with a total length of about 80 arcsec (2.8 kpc), separated by a CO-depleted funnel with a width of about 5 arcsec (175 pc). The H-alpha emission is more extended and broader than the CO emission with its maximum just in between the two CO ridges. In CO we see a peculiar velocity distribution in the iso-velocity map and p-v diagrams. We discuss different scenarios for an interpretation and present a model which can explain the observational results consistently. We propose here that the concentration of CO along the ridges is due to interaction of the rotating gas clouds with the jet's magnetic field by ambipolar diffusion. This magnetic interaction is thought to increase the time the molecular clouds reside near the jet thus leading to the quasi-static CO ridge. △ Less

Submitted 13 March, 2007; originally announced March 2007.

Comments: 13 pages, accepted by A&A

Disk Winds, Jets, and Outflows: Theoretical and Computational Foundations

Authors: Ralph E. Pudritz, Rachid Ouyed, Christian Fendt, Axel Brandenburg

Abstract: We review advances in the theoretical and computational studies of disk winds, jets and outflows including: the connection between accretion and jets, the launch of jets from magnetized disks, the coupled evolution of jets and disks, the interaction of magnetized young stellar objects with their surrounding disks and the relevance to outflows, and finally, the link between jet formation and grav… ▽ More We review advances in the theoretical and computational studies of disk winds, jets and outflows including: the connection between accretion and jets, the launch of jets from magnetized disks, the coupled evolution of jets and disks, the interaction of magnetized young stellar objects with their surrounding disks and the relevance to outflows, and finally, the link between jet formation and gravitational collapse. We also address the predictions that the theory makes about jet kinematics, collimation, and rotation, that have recently been confirmed by high spatial and spectral resolution observations. Disk winds have a universal character that may account for jets and outflows during the formation of massive stars as well as brown dwarfs. △ Less

Submitted 21 March, 2006; originally announced March 2006.

Comments: 18 pages, 5 figures, review to appear in Protostars and Planets V, B. Reipurth, D. Jewitt, and K. Keil (eds.), University of Arizona Press, Tucson, 2006

Cannonballs in the context of Gamma Ray Bursts: Formation sites ?

Authors: Jan E. Staff, Christian Fendt, Rachid Ouyed

Abstract: We investigate possible formation sites of the cannonballs (in the gamma ray bursts context) by calculating their physical parameters, such as density, magnetic field and temperature close to the origin. Our results favor scenarios where the cannonballs form as instabilities (knots) within magnetized jets from hyperaccreting disks. These instabilities would most likely set in beyond the light cy… ▽ More We investigate possible formation sites of the cannonballs (in the gamma ray bursts context) by calculating their physical parameters, such as density, magnetic field and temperature close to the origin. Our results favor scenarios where the cannonballs form as instabilities (knots) within magnetized jets from hyperaccreting disks. These instabilities would most likely set in beyond the light cylinder where flow velocity with Lorentz factors as high as 2000 can be achieved. Our findings challenge the cannonball model of gamma ray bursts if these indeed form inside core-collapse supernovae (SNe) as suggested in the literature; unless hyperaccreting disks and the corresponding jets are common occurrences in core-collapse SNe. △ Less

Submitted 5 December, 2005; originally announced December 2005.

Comments: 10 pages, 12 figures

Collimation of astrophysical jets - the role of the accretion disk magnetic field distribution

Authors: Christian Fendt

Abstract: We have applied axisymmetric MHD simulations to investigate the impact of the accretion disk magnetic flux profile on the jet collimation. Using the ZEUS-3D code modified for magnetic diffusivity, our simulations evolve from an initial hydrostatic equilibrium state in a force-free magnetic field. Considering a power law for the disk poloidal magnetic field profile Bp ~ r^{-mu} and for the disk w… ▽ More We have applied axisymmetric MHD simulations to investigate the impact of the accretion disk magnetic flux profile on the jet collimation. Using the ZEUS-3D code modified for magnetic diffusivity, our simulations evolve from an initial hydrostatic equilibrium state in a force-free magnetic field. Considering a power law for the disk poloidal magnetic field profile Bp ~ r^{-mu} and for the disk wind density profile rho ~ r^{-mu_rho} we performed a systematic study over a wide parameter range mu and mu_rho. We find a degree of collimation (ratio of mass flow rates in axial and lateral direction) decreasing for steeper disk magnetic field profiles (increasing mu). Varying the total magnetic flux doesn't change the degree of jet collimation substantially, it only affects the time scale of outflow evolution and the terminal jet speed. As our major result we find a general relation between the collimation degree with the disk wind magnetization power law exponent. Outflows with high collimation degree resulting from a flat disk magnetic field profile tend to be unsteady, producing axially propagating knots as discussed earlier. Depending slightly on the inflow density profile this unsteady behavior sets in for mu < 0.4. We also performed simulations of jet formation with artificially enhanced decay of the toroidal magnetic field in order to investigate the idea of a purely "poloidal collimation" discussed in the literature. These outflows remain weakly collimated and propagate with lower velocity. Thanks to our large numerical grid size (7x14 AU for protostars), we may apply our results to recently observed hints of jet rotation (DG Tau) indicating a relatively flat disk magnetic field profile, mu ~ 0.5. In general, our results are applicable to both stellar and extragalactic sources of MHD jets. △ Less

Submitted 25 July, 2006; v1 submitted 21 November, 2005; originally announced November 2005.

Comments: accepted by ApJ, high resolution version under www.mpia-hd.mpg.de/homes/fendt/

Journal ref: Astrophys.J.651:272-287,2006