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.

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 $\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 $\beta < 10^7$.

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.

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\'enic 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_{\phi}$-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.

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.

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.

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.

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.

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, $\mu_{\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\'en 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 $\mu(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 $\mu_{\rm D} \simeq 0.01$ that separates the regimes of magnetocentrifugally-driven and magnetic pressure-driven jets.

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.

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.

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 $\alpha^2\Omega$ 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.

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.

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.

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.

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 \alpha-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.

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.

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.

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.

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.

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.

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.

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 \alpha . 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 \alpha close to unity. However, disk stability to fragmentation due to thermal effects and gravitational instability would require a lower \alpha value. For \alpha = 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.

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.

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).

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.

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.

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.

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.

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.

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.

We analyze the photometric data obtained by PLANET and OGLE on the caustic-crossing binary-lens microlensing event OGLE-2002-BLG-069. Thanks to the excellent photometric and spectroscopic coverage of the event, we are able to constrain the lens model up to the known ambiguity between close and wide binary lenses. The detection of annual parallax in combination with measurements of extended-source effects allows us to determine the mass, distance and velocity of the lens components for the competing models. While the model involving a close binary lens leads to a Bulge-Disc lens scenario with a lens mass of M=(0.51 +- 0.15) M_sol and distance of D_L=(2.9 +- 0.4) kpc, the wide binary lens solution requires a rather implausible binary black-hole lens (M >=126 M_sol). Furthermore we compare current state-of-the-art numerical and empirical models for the surface brightness profile of the source, a G5III Bulge giant. We find that a linear limb-darkening model for the atmosphere of the source star is consistent with the data whereas a PHOENIX atmosphere model assuming LTE and with no free parameter does not match our observations.

We discuss high-resolution, time-resolved spectra of the caustic exit of the binary microlensing event OGLE 2002-BUL-69 obtained with UVES on the VLT. The source star is a G5III giant in the Galactic Bulge. During such events, the source star is highly magnified, and a strong differential magnification around the caustic resolves its surface. Using an appropriate model stellar atmosphere generated by the NextGEN code we obtained a model light curve for the caustic exit and compared it with a dense set of photometric observations obtained by the PLANET microlensing follow up network. We further compared predicted variations in the H alpha equivalent width with those measured from our spectra. While the model and observations agree in the gross features, there are discrepancies suggesting shortcomings in the model, particularly for the H alpha line core, where we have detected amplified emission from the stellar chromosphere as the source star's trailing limb exited the caustic. This achievement became possible by the provision of the OGLE-III Early Warning System, a network of small telescopes capable of nearly-continuous round-the-clock photometric monitoring, on-line data reduction, daily near-real-time modelling in order to predict caustic crossing parameters, and a fast and efficient response of a 8m-class telescope to a ``Target-Of-Opportunity'' observation request.

We present a detailed numerical study of the dynamics and evolution of ultrarelativistic magnetohydrodynamic jets in the black hole-disk system under extreme magnetization conditions. We find that Lorentz factors of up to 3000 are achieved and derived a modifiedMichel scaling (Gamma ~ sigma) which allows for a wide variation in the flow Lorentz factor. Pending contamination induced by mass-entrainment, the linear Michel scaling links modulations in the ultrarelativistic wind to variations in mass accretion in the disk for a given magnetization. The jet is asymptotically dominated by the toroidal magnetic field allowing for efficient collimation. We discuss our solutions (jets) in the context of Gamma ray bursts and describe the relevant features such as the high variability in the Lorentz factor and how high collimation angles (~ 0-5 degrees), or cylindrical jets, can be achieved. We isolate a jet instability mechanism we refer to as the "bottle-neck" instability which essentially relies on a high magnetization and a recollimation of the magnetic flux surfaces. The instability occurs at large radii where any dissipation of the magnetic energy into radiation would in principle result in an optically thin emission.

Turbulent accretion disks are the most probable origin of jets from young stellar objects. Using an extended version of ZEUS-3D code in the axisymmetry option we investigate the evolution of a disk wind into a collimated jet. Resistive jets propagate slower with increasing magnetic diffusivity, and the degree of collimation decreases. The resistive disk included in the simulations ejects persistent, well collimated jet.

We discuss the possibility to launch outflows from the close vicinity of a protoplanetary core considering a scenario where the protoplanet surrounded by a circum-planetary accretion disk is located in a circum-stellar disk. For the circum-planetary disk accretion rate we assume 6x10^{-5} M_jup/yr, implying peak disk temperatures of about 2000K. The estimated disk ionization degree and Reynolds number allow for a sufficient coupling between the disk matter and the magnetic field. We find that the surface magnetic field strength of the protoplanet is probably not more than 10 G, indicating that the global planetary magnetosphere is dominated by the circum-planetary disk magnetic field below 50 G. The estimated field strength and mass flow rates allow for asymptotic outflow velocities of 60 km/s. The overall outflow geometry will resemble a hollow tube or cone perpendicular the disk. Outflows from circum-planetary disks may be visible in shock excited emission lines along a tube of diameter of the orbital radius and thickness of about 100 protoplanetary radii. Energetically, protoplanetary outflows cannot survive the interaction with a protostellar outflow. Due to the efficient angular momentum removal by the outflow, we expect the protoplanetary outflow to influence the early planet angular momentum evolution.

The evolution of an initially stellar dipole type magnetosphere interacting with an accretion disk is investigated using numerical ideal MHD simulations. The simulations follow several 1000 Keplerian periods of the inner disk. The initial dipole type magnetic field develops into a spherically radial outflow pattern with two main components both evolving into a quasi-stationary final state. A neutral field line divides both components, along which small plasmoids are ejected in irregular time intervals. The half opening angle of the stellar wind cone varies from 30 to 55 degrees depending on the ratio of the mass flow rates of disk wind and stellar wind. The maximum speed of the outflow is about the Keplerian speed at the inner disk radius. An axial jet forms during the first decades of rotation. However, this feature does not survive on the very long time scale and a pressure driven low velocity flow along the axis evolves.

We investigate the evolution of a disk wind into a collimated jet under the influence of magnetic diffusivity. Using the ZEUS-3D code in the axisymmetry option we solve the time-dependent resistive MHD equations for a model setup of a central star surrounded by an accretion disk. The disk is taken as a time-independent boundary condition for the mass flow rate and the magnetic flux distribution. We find that the diffusive jets propagate slower into the ambient medium, most probably due to the lower mass flow rate in axial direction. Close to the star we find that a quasi stationary state evolves after several hundreds (weak diffusion) or thousands (strong diffusion) of disk rotations. Magnetic diffusivity affects the protostellar jet structure de-collimating it. We explain these effects in the framework of the Lorentz force.

Highly relativistic jets are most probably driven by strong magnetic fields and are launched from the accretion disk surrounding a central black hole. In this paper we review some of our recent results considering the two-dimensional magnetic field structure and the dynamics of collimating relativistic jets. Applying the stationarity assumption enables us to calculate global solutions of the relevant MHD equations. We discuss solutions of the Grad-Shafranov equation obtained in a general relativistic context applying the 3+1 formalism for Kerr geometry. These solutions extend from the inner light surface around the Kerr black hole to the asymptotic regime of a cylindrically collimated jet with a finite radius. In a further step, we include differential rotation of the foot points of the field lines. Considering the general relativistic MHD wind equation, we investigate the dynamics of the collimating jet, in particular the effects of Kerr metric on the acceleration. Temperature and density follow a power law. Having at hand a relativistic MHD jet solution, we calculate the thermal optically thin X-ray spectrum for the innermost hot part of the jet. Doppler shift and boosting is taken into account. For microquasars we obtain a jet X-ray luminosity 10^33 erg/s. Iron emission lines are clearly visible.

High velocity highly collimated beams of plasma -- the jets -- are known as a general phenomenon among astrophysical sources of different energy and spatial scale. The common MHD scenario of astrophysical jet formation allows to describe the different jet sources by a scaling of the jet magnetization. In this paper we apply the MHD jet theory for GRB parameters and find indeed that main features as ultra-relativistic motion and strong collimation can be explained. Others, like GRB time-variability and the bi-modal GRB duration are (yet) beyond a MHD description.

We investigate the evolution of a disk wind into a collimated jet under the influence of magnetic diffusivity, assuming that the turbulent pattern in the disk will also enter the disk corona and the jet. Using the ZEUS-3D code in the axisymmetry option we solve the time-dependent resistive MHD equations for a model setup of a central star surrounded by an accretion disk. We find that the diffusive jets propagate slower into the ambient medium. Close to the star we find that a quasi stationary state evolves after several hundred (weak diffusion) or thousand (strong diffusion) disk rotations. Magnetic diffusivity affects the protostellar jet structure as follows. The jet poloidal magnetic field becomes de-collimated. The jet velocity increases with increasing diffusivity, while the degree of collimation for the hydrodynamic flow remains more or less the same. We suggest that the mass flux is a proper tracer for the degree of jet collimation and find indications of a critical value for the magnetic diffusivity above which the jet collimation is only weak.

Highly relativistic jets are most probably driven by strong magnetic fields and launched from the accretion disk surrounding a central black hole. Applying the jet flow parameters (velocity, density, temperature) calculated from the magnetohydrodynamic (MHD) equations, we derive the thermal X-ray luminosity along the inner jet flow in the energy range 0.2-10.1 keV. Here, we concentrate on the case of Galactic microquasars emitting highly relativistic jets. For a 5 solar masses central object and a jet mass flow rate of 10^{-8} solar masses/year we obtain a jet X-ray luminosity L_x ~ 10^{33} erg\sec. Emission lines of Fe XXV and Fe XXVI are clearly visible. Relativistic effects such as Doppler shift and Doppler boosting were considered for different inclinations of the jet axis. Due to the chosen geometry of the MHD jet the inner X-ray emitting part is not yet collimated. Therefore, depending on the viewing angle, the Doppler boosting does not play a major role in the total spectra.

We have investigated magnetically driven superluminal jets originating from rotating black holes. The stationary, general relativistic, magnetohydrodynamic wind equation along collimating magnetic flux surfaces has been solved numerically. In particular, we apply our results to the Galactic superluminal sources GRS 1915+105 and GRO 1655-40. Our results show that the observed speed of more than 0.9 c can be achieved in general by magnetohydrodynamic acceleration. The velocity distribution along the magnetic field has a saturating profile. Further, we find that highly relativistic speeds can be reached also for jets not emerging from a region close to the black hole, if the flow magnetization is sufficiently large. The plasma temperature rapidly decreases from about 10**10 K at the foot point of the jet to about 10**6 K at a distance of 5000 gravitational radii from the source. Temperature and the mass density follow a power law distribution with the radius.

The evolution of an initially stellar dipole type magnetosphere interacting with an accretion disk is investigated numerically using the ideal MHD ZEUS-3D code in the 2D-axisymmetry option. The simulations last several thousands of Keplerian periods of the inner disk. A Keplerian disk is assumed as a boundary condition prescribing a mass inflow into the corona. Additionally, a stellar wind from a rotating central star is prescribed. Our major result is that the initially dipole type field develops into a spherically radial outflow pattern with two main components, a disk wind and a stellar wind component. These components evolve into a quasi- stationary final state. The poloidal field lines follow a conical distribution. As a consequence of the initial dipole, the field direction in the stellar wind is opposite to that in the disk wind. The maximum speed of the outflow is about the Keplerian speed at the inner disk radius. With the chosen mass flow rates and field strength we see almost no indication for a flow self-collimation.

The magnetic flux distribution is determined by the solution of the Grad-Shafranov equation. With differential rotation, i.e. the variation of the iso-rotation parameter, the shape of the light surface must be calculated in an iterative way. For the first time, we have calculated the force-free magnetic structure of truly two-dimensional, relativistic jets, anchored in a differentially rotating disk. Such an approach allows for a direct connection between parameters of the central source (mass, rotation) and the extension of the radio jet. We present an analytical estimate for the jet opening angle along the asymptotic branches of the light surface. In general, differential rotation of the iso-rotation parameter leads to an increase of the jet opening angle. Comparison to the M87 jet shows agreement in the collimation distance. We derive a light cylinder radius of the M87 jet of 50 Schwarzschild radii.

The influence of strong, large-scale magnetic fields on the structure and temperature distribution in white dwarf atmospheres is investigated. Magnetic fields may provide an additional component of pressure support, thus possibly inflating the atmosphere compared to the non-magnetic case. Since the magnetic forces are not isotropic, atmospheric properties may significantly deviate from spherical symmetry. In this paper the magnetohydrostatic equilibrium is calculated numerically in the radial direction for either for small deviations from different assumptions for the poloidal current distribution. We generally find indication that the scale height of the magnetic white dwarf atmosphere enlarges with magnetic field strength and/or poloidal current strength. This is in qualitative agreement with recent spectropolarimetric observations of Grw$+10\degr8247$. Quantitatively, we find for e.g. a mean surface poloidal field strength of 100 MG and a toroidal field strength of 2-10 MG an increase of scale height by a factor of 10. This is indicating that already a small deviation from the initial force-free dipolar magnetic field may lead to observable effects. We further propose the method of finite elements for the solution of the two-dimensional magnetohydrostatic equilibrium including radiation transport in the diffusive approximation. We present and discuss preliminary solutions, again indicating on an expansion of the magnetized atmosphere.

We present, for the first time, the structure of the axisymmetric force-free magnetosphere of an aligned rotating magnetic dipole, in the case in which there exists a sufficiently large charge density (whose origin we do not question) to satisfy the ideal MHD condition, ${\bf E\cdot B}=0$, everywhere. The unique distribution of electric current along the open magnetic field lines which is required for the solution to be continuous and smooth is obtained numerically. With the geometry of the field lines thus determined we compute the dynamics of the associated MHD wind. The main result is that the relativistic outflow contained in the magnetosphere is not accelerated to the extremely relativistic energies required for the flow to generate gamma rays. We expect that our solution will be useful as the starting point for detailed studies of pulsar magnetospheres under more general conditions, namely when either the force-free and/or the ideal MHD condition ${\bf E\cdot B}=0$ are not valid in the entire magnetosphere. Based on our solution, we consider that the most likely positions of such an occurrence are the polar cap, the crossings of the zero space charge surface by open field lines, and the return current boundary, but not the light cylinder.