BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Oliver Porth

Oliver Porth appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-ph.SRastro-ph.GAastro-phastro-ph.IMphysics.comp-phphysics.plasm-ph

Catalog footprint

What is connected

24works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

GRMHD Simulations of Accreting Neutron Stars with Non-Dipole Fields

NASA's NICER telescope has recently provided evidence for non-dipolar magnetic field structures in rotation-powered millisecond pulsars. These stars are assumed to have gone through a prolonged accretion spin-up phase, begging the question of what accretion flows onto stars with complex magnetic fields would look like. We present results from a suite of GRMHD simulations of accreting neutron stars for dipole, quadrupole, and quadrudipolar stellar field geometries. This is a first step towards simulating realistic hotspot shapes in a general relativistic framework to understand hotspot variability in accreting millisecond pulsars. We find that the location and size of the accretion columns resulting in hotspots changes significantly depending on initial stellar field strength and geometry. We also find that the strongest contributions to the stellar torque are from disk-connected fieldlines and the pulsar wind, leading to spin-down in almost all of the parameter regime explored here. We further analyze angular momentum transport in the accretion disk due to large scale magnetic stresses, turbulent stresses, wind- and compressible effects which we identify with convective motions. The disk collimates the initial open stellar flux forming jets. For dipoles, the disk-magnetosphere interaction can either enhance or reduce jet power compared to the isolated case. However for quadrupoles, the disk always leads to an enhanced net open flux making the jet power comparable to the dipolar case. We discuss our results in the context of observed neutron star jets and provide a viable mechanism to explain radio power both in the low- and high-magnetic field case.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Magnetic reconnection and plasmoid formation in three-dimensional accretion flows around black holes

Magnetic reconnection is thought to be one of the main energy-dissipation mechanisms fueling energy to the plasma in the vicinity of a black hole. Indeed, plasmoids formed through magnetic reconnection may play a key role in $γ$-ray, X-ray and near-infrared flares from the black hole at the center of our galaxy, SgrA*. We report the results of three-dimensional general-relativistic ideal and resistive magnetohydrodynamics simulations modelling magnetic reconnection in accretion flows around astrophysical black holes. As an important difference with similar works, our accretion discs have an initial dipolar magnetic-field configuration with loops of alternating polarity. We show that current sheets are formed and destroyed rapidly in the turbulent environment of black-hole accretion. Plasmoids are formed from current sheets close to the event horizon, in a region of $\sim2-15$ gravitational radii. We further quantify the magnetic dissipation and the process of energy transfer to the plasmoids, reporting the reconnection rate, the relative current density with respect to the local magnetic field, and the size of the plasmoids. We find that plasmoids gain energy through reconnection and heat up to relativistic temperatures, with the largest ones being sufficiently energetic to leave the black hole near the polar regions. During their evolution, plasmoids are stretched and elongated, becoming disrupted when the shear is sufficiently large, although some plasmoids survive as well-distinguished structures at distances of $\sim30-40$ gravitational radii from the black hole. Finally, we find that in some cases the plasmoids acquire a super-Keplerian azimuthal velocity, as suggested by recent observations of flares from Sgr~A*.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

3D magnetised jet break-out from neutron-star binary merger ejecta: afterglow emission from the jet and the ejecta

We perform three-dimensional (3D) general-relativistic magnetohydrodynamic simulations to model the jet break-out from the ejecta expected to be produced in a binary neutron-star merger. The structure of the relativistic outflow from the 3D simulation confirms our previous results from 2D simulations, namely, that a relativistic magnetized outflow breaking out from the merger ejecta exhibits a hollow core of $θ_{\rm core}\approx4^{\circ}$, an opening angle of $θ_{\rm jet}\gtrsim10^{\circ}$, and is accompanied by a wind of ejected matter that will contribute to the kilonova emission. We also compute the non-thermal afterglow emission of the relativistic outflow and fit it to the panchromatic afterglow from GRB170817A, together with the superluminal motion reported from VLBI observations. In this way, we deduce an observer angle of $θ_{\rm obs}= 35.7^{\circ \,\,+1.8}_{\phantom{\circ \,\,}-2.2}$. We further compute the afterglow emission from the ejected matter and constrain the parameter space for a scenario in which the matter responsible for the thermal kilonova emission will also lead to a non-thermal emission yet to be observed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Impact of non-thermal particles on the spectral and structural properties of M87

The recent 230 GHz observations of the Event Horizon Telescope (EHT) are able to image the innermost structure of the M87 and show a ring-like structure which is in agreement with thermal synchrotron emission generated in a torus surrounding a supermassive black hole. However, at lower frequencies M87 is characterised by a large-scale and edge-brightened jet with clear signatures of non-thermal emission. In order to bridge the gap between these scales and to provide a theoretical interpretation of these observations we perform general relativistic magnetohydrodynamic simulations of accretion on to black holes and jet launching. M87 has been the target for multiple observations across the entire electromagnetic spectrum. Among these VLBI observations provide unique details on the collimation profile of the jet down to several gravitational radii. In this work we aim to model the observed broad-band spectrum of M87 from the radio to the NIR regime and at the same time fit the jet structure as observed with Global mm-VLBI at 86 GHz. We use general relativistic magnetohydrodynamics and simulate the accretion of the magnetised plasma onto Kerr-black holes in 3D. The radiative signatures of these simulations are computed taking different electron distribution functions into account and a detailed parameter survey is performed in order to match the observations. The results of our simulations show that magnetically arrested disks around fast spinning black holes ($a_\star\geq0.5$) together with a mixture of thermal and non-thermal particle distributions are able to model simultaneously the broad-band spectrum and the innermost jet structure of M87

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deep Horizon; a machine learning network that recovers accreting black hole parameters

The Event Horizon Telescope recently observed the first shadow of a black hole. Images like this can potentially be used to test or constrain theories of gravity and deepen the understanding in plasma physics at event horizon scales, which requires accurate parameter estimations. In this work, we present Deep Horizon, two convolutional deep neural networks that recover the physical parameters from images of black hole shadows. We investigate the effects of a limited telescope resolution and observations at higher frequencies. We trained two convolutional deep neural networks on a large image library of simulated mock data. The first network is a Bayesian deep neural regression network and is used to recover the viewing angle $i$, and position angle, mass accretion rate $\dot{M}$, electron heating prescription $R_{\rm high}$ and the black hole mass $M_{\rm BH}$. The second network is a classification network that recovers the black hole spin $a$. We find that with the current resolution of the Event Horizon Telescope, it is only possible to accurately recover a limited number of parameters of a static image, namely the mass and mass accretion rate. Since potential future space-based observing missions will operate at frequencies above 230 GHz, we also investigated the applicability of our network at a frequency of 690 GHz. The expected resolution of space-based missions is higher than the current resolution of the Event Horizon Telescope, and we show that Deep Horizon can accurately recover the parameters of simulated observations with a comparable resolution to such missions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How to tell an accreting boson star from a black hole

The capability of the Event Horizon Telescope (EHT) to image the nearest supermassive black hole candidates at horizon-scale resolutions offers a novel means to study gravity in its strongest regimes and to test different models for these objects. Here, we study the observational appearance at 230 GHz of a surfaceless black hole mimicker, namely a non-rotating boson star, in a scenario consistent with the properties of the accretion flow onto Sgr A*. To this end, we perform general relativistic magnetohydrodynamic simulations followed by general relativistic radiative transfer calculations in the boson star space-time. Synthetic reconstructed images considering realistic astronomical observing conditions show that, despite qualitative similarities, the differences in the appearance of a black hole -- either rotating or not -- and a boson star of the type considered here are large enough to be detectable. These differences arise from dynamical effects directly related to the absence of an event horizon, in particular, the accumulation of matter in the form of a small torus or a spheroidal cloud in the interior of the boson star, and the absence of an evacuated high-magnetization funnel in the polar regions. The mechanism behind these effects is general enough to apply to other horizonless and surfaceless black hole mimickers, strengthening confidence in the ability of the EHT to identify such objects via radio observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

MPI-AMRVAC: a parallel, grid-adaptive PDE toolkit

We report on the latest additions to our open-source, block-grid adaptive framework MPI-AMRVAC, which is a general toolkit for especially hyperbolic/parabolic partial differential equations (PDEs). Applications traditionally focused on shock-dominated, magnetized plasma dynamics described by either Newtonian or special relativistic (magneto)hydrodynamics, but its versatile design easily extends to different PDE systems. Here, we demonstrate applications covering any-dimensional scalar to system PDEs, with e.g. Korteweg-de Vries solutions generalizing early findings on soliton behaviour, shallow water applications in round or square pools, hydrodynamic convergence tests as well as challenging computational fluid and plasma dynamics applications. The recent addition of a parallel multigrid solver opens up new avenues where also elliptic constraints or stiff source terms play a central role. This is illustrated here by solving several multi-dimensional reaction-diffusion-type equations. We document the minimal requirements for adding a new physics module governed by any nonlinear PDE system, such that it can directly benefit from the code flexibility in combining various temporal and spatial discretisation schemes. Distributed through GitHub, MPI-AMRVAC can be used to perform 1D, 1.5D, 2D, 2.5D or 3D simulations in Cartesian, cylindrical or spherical coordinate systems, using parallel domain-decomposition, or exploiting fully dynamic block quadtree-octree grids.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the opening angle of magnetised jets from neutron-star mergers: the case of GRB170817A

The observations of GW170817/GRB170817A have confirmed that the coalescence of a neutron-star binary is the progenitor of a short gamma-ray burst. In the standard picture of a short gamma-ray burst, a collimated highly relativistic outflow is launched after merger and it successfully breaks out from the surrounding ejected matter. Using initial conditions inspired from numerical-relativity binary neutron-star merger simulations, we have performed general-relativistic hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations in which the jet is launched and propagates self-consistently. The complete set of simulations suggests that: (i) MHD jets have an intrinsic energy and velocity polar structure with a ``hollow core'' subtending an angle $θ_{\rm core}\approx4^{\circ}-5^{\circ}$ and an opening angle of $θ_{\rm jet}\gtrsim10^{\circ}$; (ii) MHD jets eject significant amounts of matter and two orders of magnitude more than HD jets; (iii) the energy stratification in MHD jets naturally yields the power-law energy scaling $E(>Γβ)\propto(Γβ)^{-4.5}$; (iv) MHD jets provide fits to the afterglow data from GRB170817A that are comparatively better than those of the HD jets and without free parameters; (v) finally, both of the best-fit HD/MHD models suggest an observation angle $θ_{\rm obs} \simeq 21^{\circ}$ for GRB170817A.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Plasmoid formation in global GRMHD simulations and AGN flares

One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Modelling the polarised emission from black holes on event horizon-scales

Upcoming VLBI observations will resolve nearby supermassive black holes, most notably Sagittarius A* and M87, on event horizon-scales. Recent observations of Sagittarius A* with the Event Horizon Telescope have revealed horizon-scale structure. Accordingly, the detection and measurement of the back hole ``shadow" is expected to enable the existence of astrophysical black holes to be verified directly. Although the theoretical description of the shadow is straightforward, its observational appearance is largely determined by the properties of the surrounding accretion flow, which is highly turbulent. We introduce a new polarised general-relativistic radiative transfer code, \texttt{BHOSS}, which accurately solves the equations of polarised radiative transfer in arbitrary strong-gravity environments, providing physically-realistic images of astrophysical black holes on event horizon-scales, as well as also providing insight into the fundamental properties and nature of the surrounding accretion flow environment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solar prominences: 'double, double ... boil and bubble'

Observations revealed rich dynamics within prominences, the cool 10,000 K, macroscopic (sizes of order 100 Mm) "clouds" in the million degree solar corona. Even quiescent prominences are continuously perturbed by hot, rising bubbles. Since prominence matter is hundredfold denser than coronal plasma, this bubbling is related to Rayleigh-Taylor instabilities. Here we report on true macroscopic simulations well into this bubbling phase, adopting a magnetohydrodynamic description from chromospheric layers up to 30 Mm height. Our virtual prominences rapidly establish fully non-linear (magneto)convective motions where hot bubbles interplay with falling pillars, with dynamical details including upwelling pillars forming within bubbles. Our simulations show impacting Rayleigh-Taylor fingers reflecting on transition region plasma, ensuring that cool, dense chromospheric material gets mixed with prominence matter up to very large heights. This offers an explanation for the return mass cycle mystery for prominence material. Synthetic views at extreme ultraviolet wavelengths show remarkable agreement with observations, with clear indications of shear-flow induced fragmentations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The inner knot of the Crab nebula

We model the inner knot of the Crab Nebula as a synchrotron emission coming from the non-spherical MHD termination shock of relativistic pulsar wind. The post-shock flow is mildly relativistic; as a result the Doppler-beaming has a strong impact on the shock appearance. The model can reproduce the knot location, size, elongation, brightness distribution, luminosity and polarization provided the effective magnetization of the section of the pulsar wind producing the knot is low, $σ\leq 1$. In the striped wind model, this implies that the striped zone is rather wide, with the magnetic inclination angle of the Crab pulsar $\ge 45^\circ$; this agrees with the previous model-dependent estimate based on the gamma-ray emission of the pulsar. We conclude that the tiny knot is indeed a bright spot on the surface of a quasi-stationary magnetic relativistic shock and that this shock is a site of efficient particle acceleration. On the other hand, the deduced low magnetization of the knot plasma implies that this is an unlikely site for the Crab's gamma-ray flares, if they are related to the fast relativistic magnetic reconnection events.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interacting tilt and kink instabilities in repelling current channels

We present a numerical study in resistive magnetohydrodynamics where the initial equilibrium configuration contains adjacent, oppositely directed, parallel current channels. Since oppositely directed current channels repel, the equilibrium is liable to an ideal magnetohydrodynamic tilt instability. This tilt evolution, previously studied in planar settings, involves two magnetic islands or fluxropes, which on Alfvenic timescales undergo a combined rotation and separation. This in turn leads to the creation of (near) singular current layers, posing severe challenges to numerical approaches. Using our open-source grid-adaptive MPI-AMRVAC software, we revisit the planar evolution case in compressible MHD, as well as its extension to 2.5D and full 3D scenarios. As long as the third dimension remains ignorable, pure tilt evolutions result which are hardly affected by out of plane magnetic field components. In all 2.5D runs, our simulations do show secondary tearing type disruptions throughout the near singular current sheets in the far nonlinear saturation regime. In full 3D, both current channels can be liable to additional ideal kink deformations. We discuss the effects of having both tilt and kink instabilities acting simultaneously in the violent, reconnection dominated evolution. In 3D, both the tilt and the kink instabilities can be stabilized by tension forces. As a concrete space plasma application, we argue that interacting tilt-kink instabilities in repelling current channels provide a novel route to initiate solar coronal mass ejections, distinctly different from currently favored pure kink or torus instability routes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rayleigh-Taylor instability in Magnetohydrodynamic Simulations of the Crab Nebula

In this paper we discuss the development of Rayleigh-Taylor filaments in axisymmetric simulations of Pulsar wind nebulae (PWN). High-resolution adaptive mesh refinement magnetohydrodynamic (MHD) simulations are used to resolve the non-linear evolution of the instability. The typical separation of filaments is mediated by the turbulent flow in the nebula and hierarchical growth of the filaments. The strong magnetic dissipation and field-randomization found in recent global three-dimensional simulations of PWN suggests that magnetic tension is not strong enough to suppress the growth of RT filaments, in agreement with the observations of prominent filaments in the Crab nebula. The long-term axisymmetric results presented here confirm this finding.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic AGN jets II. Jet properties and mixing effects for episodic jet activity

Various radio galaxies show signs of having gone through episodic jet outbursts in the past. An example is the class of double-double radio galaxies (DDRGs). However, to follow the evolution of an individual source in real-time is impossible due to the large time scales involved. Numerical studies provide a powerful tool to investigate the temporal behavior of episodic jet outbursts in a (magneto-)hydrodynamical setting. We simulate the injection of two jets from active galactic nuclei (AGN), separated by a short interruption time. Three different jet models are compared. We find that an AGN jet outburst cycle can be divided into four phases. The most prominent phase occurs when the restarted jet is propagating completely inside the hot and inflated cocoon left behind by the initial jet. In that case, the jet-head advance speed of the restarted jet is significantly higher than the initial jet-head. While the head of the initial jet interacts strongly with the ambient medium, the restarted jet propagates almost unimpeded. As a result, the restarted jet maintains a strong radial integrity. Just a very small fraction of the amount of shocked jet material flows back through the cocoon compared to that of the initial jet and much weaker shocks are found at the head of the restarted jet. We find that the features of the restarted jet in this phase closely resemble the observed properties of a typical DDRG.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the in situ condensation process of solar prominences

Prominences in the solar corona are hundredfold cooler and denser than their surroundings, with a total mass of 1.e13 up to 1.e15 g. Here we report on the first comprehensive simulations of three-dimensional, thermally and gravitationally stratified magnetic flux ropes, where in situ condensation to a prominence happens due to radiative losses. After a gradual thermodynamic adjustment, we witness a phase where runaway cooling happens while counter-streaming shearing flows drain off mass along helical field lines. After this drainage, a prominence-like condensation resides in concave upward field regions, and this prominence retains its overall characteristics for more than two hours. While condensing, the prominence establishes a prominence-corona transition region, where magnetic field-aligned thermal conduction is operative during the runaway cooling. The prominence structure represents a force-balanced state in a helical flux rope. The simulated condensation demonstrates a right-bearing barb, as a remnant of the drainage. Synthetic images at extreme ultraviolet wavelengths follow the onset of the condensation, and confirm the appearance of horns and a three-part structure for the stable prominence state, as often seen in erupting prominences. This naturally explains recent Solar Dynamics Observatory views with the Atmospheric Imaging Assembly on prominences in coronal cavities demonstrating horns.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic 3D jet simulations for the X-ray binary SS433

Context. Modern high resolution observations allow to view closer into the objects powering relativistic jets. This is especially the case for SS433, an X-ray binary from which a precessing jet is observed down to the sub-parsec scale. Aims. We want to study full 3D dynamics of relativistic jets associated with AGN or XRB. We study the precessing motion of a jet as a model for the jet associated with the XRB SS433. Our study of the jet dynamics in this system focuses on the sub-parsec scales. We investigate the impact of jet precession and the variation of the Lorentz factor of the injected matter on the general 3D jet dynamics and its energy transfer to the surrounding medium. We realize synthetic radio mapping of the data, to compare our results with observations. Methods. For our study we use the code MPI-AMRVAC with SRHD model of a baryonic jet. We use a AMR scheme and an inner time-dependent boundary prescription to inject the jets. Parameters extracted from observations were used. 3D jet realizations that match the SS433 jet are intercompared. We track the energy content, as deposited in different regions of the domain affected by the jet. Our code also follows a population of particles injected with the jet. This evolving energy spectrum of accelerated electrons, allows to obtain the radio emission from our simulation. Results. we obtain meaningful observations. We find increased energy transfer for a precessing jet compared to a standing jets. We obtain synthetic radio maps for all jets. Conclusions. The synthetic radio map matches best for a model using the canonical kinematic model. Overdense precessing jets experience significant deceleration in their propagation through the ISM, and while the overall jet is of helical shape. This argument show that the kinematic model for SS433 assuming ballistic propagation has to be corrected for this decelretaion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three-Dimensional Magnetohydrodynamic Simulations of the Crab Nebula

In this paper we give a detailed account of the first 3D relativistic magnetohydrodynamic (MHD) simulations of Pulsar Wind Nebulae (PWN), with parameters most suitable for the Crab Nebula. In order to clarify the new features specific to 3D models, reference 2D simulations have been carried out as well. Compared to the previous 2D simulations, we considered pulsar winds with much stronger magnetisation, up to σ=3, and accounted more accurately for the anticipated magnetic dissipation in the striped zone of these winds. While the 3D models preserve the separation of the post termination shock flow into the equatorial and polar components, their relative strength and significance differ. Whereas the highly magnetised 2D models produce highly coherent and well collimated polar jets capable of efficient "drilling" through the supernova shell, in the corresponding 3D models the jets are disrupted by the kink mode current driven instability and "dissolve" into the main body of PWN after propagation of several shock radii. Our results show that contrary to the expectations based on 1D analytical and semi-analytical models, our numerical solutions with highly magnetized pulsar winds still produce termination shocks with radii comparable to those deduced from the observations. We present polarized synchrotron maps constructed from the 3D simulations, showing that the wealth of morphological features observed with Hubble and Chandra can well be reproduced with 3D-MHD models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Evolution of growing black holes in axisymmetric galaxy cores

NBody realizations of axisymmetric collisional galaxy cores (e.g. M32, M33, NGC205, Milky Way) with embedded growing black holes are presented. Stars which approach the disruption sphere are disrupted and accreted to the black hole. We measure the zone of influence of the black hole and disruption rates in relaxation time scales. We show that secular gravitational instabilities dominate the initial core dynamics, while the black hole is small and growing due to consumption of stars. Later, the black hole potential dominates the core, and loss cone theory can be applied. Our simulations show that central rotation in galaxies can not be neglected for relaxed systems, and compare and discuss our results with the standard theory of spherically symmetric systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Jet formation from massive young stars: Magnetohydrodynamics versus radiation pressure

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Synchrotron radiation of self-collimating relativistic MHD jets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Escape from the vicinity of fractal basin boundaries of a star cluster

The dissolution process of star clusters is rather intricate for theory. We investigate it in the context of chaotic dynamics. We use the simple Plummer model for the gravitational field of a star cluster and treat the tidal field of the Galaxy within the tidal approximation. That is, a linear approximation of tidal forces from the Galaxy based on epicyclic theory in a rotating reference frame. The Poincaré surfaces of section reveal the effect of a Coriolis asymmetry. The system is non-hyperbolic which has important consequences for the dynamics. We calculated the basins of escape with respect to the Lagrangian points $L_1$ and $L_2$. The longest escape times have been measured for initial conditions in the vicinity of the fractal basin boundaries. Furthermore, we computed the chaotic saddle for the system and its stable and unstable manifolds. The chaotic saddle is a fractal structure in phase space which has the form of a Cantor set and introduces chaos into the system.

0 citations0 reviewsNo code linkedOpen access