Normalized to: Kissmann, R.
[1]
oai:arXiv.org:1812.02960 [pdf] - 1823014
Proton acceleration in colliding stellar wind binaries
Submitted: 2018-12-07
The interaction between the strong winds in stellar colliding-wind binary
(CWB) systems produces two shock fronts, delimiting the wind collision region
(WCR). There, particles are expected to be accelerated mainly via diffusive
shock acceleration (DSA). We investigate the injection and the acceleration of
protons in typical CWB systems by means of Monte Carlo simulations, with both a
test-particle approach and a non-linear method modelling a shock locally
modified by the backreaction of the accelerated protons. We use
magnetohydrodynamic simulations to determine the background plasma in the WCR
and its vicinity. This allows us to consider particle acceleration at both
shocks, on either side of the WCR, with a realistic large-scale magnetic field.
We highlight the possible effects of particle acceleration on the local shock
profiles at the WCR. We include the effect of magnetic field amplification due
to resonant streaming instability (RSI), and compare results without and with
the backreaction of the accelerated protons. In the latter case we find a lower
flux of the non-thermal proton population, and a considerable magnetic field
amplification. This would significantly increase the synchrotron losses of
relativistic electrons accelerated in CWB systems, lowering the maximal energy
they can reach and strongly reducing the inverse Compton fluxes. As a result,
$\gamma$-rays from CWBs would be predominantly due to the decay of neutral
pions produced in nucleon-nucleon collisions. This might provide a way to
explain why, in the vast majority of cases, CWB systems have not been
identified as $\gamma$-ray sources, while they emit synchrotron radiation.
[2]
oai:arXiv.org:1811.12683 [pdf] - 1791839
The consequence of a new ISRF model of the Milky Way on predictions for
diffuse gamma-ray emission
Submitted: 2018-11-30
We investigate the impact of the recently published ISRF model of the Milky
Way by Popescu et al. on the CR electrons and positrons in the context of
cosmic-ray transport modelling. We also study predictions for diffuse Galactic
gamma radiation and underline the importance of the increased ISRF
intensities.We use the PICARD code for solving the CR transport equation to
obtain predictions based on Galactic cosmic-ray electron fluxes. We show that
the new ISRF yields gamma-ray intensity increases, most particular at the
Galactic Center and in the Galactic Plane. The impact is largest at energies
around 220 GeV.
[3]
oai:arXiv.org:1806.09479 [pdf] - 1709563
The CRONOS Code for Astrophysical Magnetohydrodynamics
Submitted: 2018-06-25, last modified: 2018-07-03
We describe the magnetohydrodynamics (MHD) code CRONOS, which has been used
in astrophysics and space physics studies in recent years. CRONOS has been
designed to be easily adaptable to the problem at hand, where the user can
expand or exchange core modules or add new functionality to the code. This
modularity comes about through its implementation using a C++ class structure.
The core components of the code include solvers for both hydrodynamical (HD)
and MHD problems. These problems are solved on different rectangular grids,
which currently support Cartesian, spherical, and cylindrical coordinates.
CRONOS uses a finite-volume description with different approximate Riemann
solvers that can be chosen at runtime. Here, we describe the implementation of
the code with a view toward its ongoing development. We illustrate the code's
potential through several (M)HD test problems and some astrophysical
applications.
[5]
oai:arXiv.org:1708.05860 [pdf] - 1587279
3D magneto-hydrodynamic models of non-thermal photon emission in the
binary system $\gamma^2$ Velorum
Submitted: 2017-08-19
Recent reports claiming association of the massive star binary system
$\gamma^2 $Velorum (WR 11) with a high-energy $\gamma$-ray source observed by
\textit{Fermi}-LAT contrast the so-far exclusive role of $\eta$ Carinae as the
hitherto only detected $\gamma$-ray emitter in the source class of
particle-accelerating colliding-wind binary (CWB) systems. We offer support to
this claim of association by providing dedicated model predictions for the
nonthermal photon emission spectrum of $\gamma^2$ Velorum.
We use three-dimensional magneto-hydrodynamic modeling (MHD) to investigate
the structure and conditions of the wind-collision region (WCR) of $\gamma^2$
Velorum including the important effect of radiative braking in the stellar
winds. A transport equation is then solved throughout the computational domain
to study the propagation of relativistic electrons and protons. The resulting
distributions of particles are subsequently used to compute nonthermal photon
emission components.
In agreement with observation in X-ray spectroscopy, our simulations yield a
large shock-cone opening angle. We find the nonthermal $\gamma$-ray emission of
$\gamma^2$ Velorum to be of hadronic origin owing to the strong radiation
fields in the binary system which inhibit the acceleration of electrons to
energies sufficiently high for efficient inverse Compton radiation. We also
discuss the strong dependence of a hadronic $\gamma$-ray component on the
energy-dependent diffusion used in the simulations. Of two mass-loss rates for
the WR star found in literature, only the higher one is able to accommodate the
observed $\gamma$-ray spectrum with reasonable values for important simulation
parameters such as the injection ratio of high-energy particles within the WCR.
[6]
oai:arXiv.org:1701.07285 [pdf] - 1535007
Diffuse Gamma Rays in 3D Galactic Cosmic-ray Propagation Models
Submitted: 2017-01-25
The Picard code for the numerical solution of the Galactic cosmic ray
propagation problem allows for high-resolution models that acknowledge the 3D
structure of our Galaxy. Picard was used to determine diffuse gamma-ray
emission of the Galaxy over the energy range from 100 MeV to 100 TeV. We
discuss the impact of a cosmic-ray source distribution aligned with the
Galactic spiral arms for a range of such spiral-arm models. As expected, the
impact on the gamma-ray emission is most distinct in the inverse-Compton
channel, where imprints of the spiral arms are visible and yield predictions
that are no longer symmetric to the rotational axis of the Milkyway. We will
illustrate these differences by a direct comparison to results from previous
axially symmetric Galactic propagation models: we find differences in the
gamma-ray flux both on global scales and on local scales related to the spiral
arm tangents. We compare gamma-ray flux and spectra at on-arm vs. off-arm
projections and characterize the differences to axially symmetric models.
[7]
oai:arXiv.org:1701.07289 [pdf] - 1535010
Combined Magnetohydrodynamic- Monte Carlo Simulations of Proton
Acceleration in Colliding Wind Binaries
Submitted: 2017-01-25
The interaction between the strong winds of the stars in colliding-wind
binary (CWB) systems produces two shock fronts, delimiting the wind collision
region (WCR). There, particles are expected to be accelerated mainly via
diffusive shock acceleration, and to produce $\gamma$-rays, in processes
involving relativistic electrons and/or protons. We investigate the injection
and the acceleration of protons in typical CWB systems by means of Monte Carlo
simulations, with a test-particle approach. We use magnetohydrodynamic
simulations to determine the background conditions in the wind collision
region. This allows us to consider particle acceleration at both shocks, on
either side of the WCR, with a self-consistently determined large-scale
magnetic field, which has an impact on the shape of the WCR, and the topology
of which plays an important role in particle acceleration at collisionless
shocks. Such studies may contribute to improve $\gamma$-ray flux predictions
for CWB systems.
[8]
oai:arXiv.org:1701.07288 [pdf] - 1535009
The Use Case of a New ISRF on Diffuse Gamma-ray Emission Models
Submitted: 2017-01-25
An important contribution to the gamma-ray emission of our Galaxy is the
Galactic Diffuse emission. We present specific developments within the PICARD
code for modeling of Galactic cosmic ray propagation that are relevant for the
computation of the gamma-ray emission to accommodate a new Interstellar
Radiation Field (ISRF) model. We study the differences of the individual ISRF
models on Galactic Diffuse gamma-ray emission predictions in the GeV to the TeV
energy regime. Emphasis was laid on obtaining robust predictions for observable
signatures in the very high energy gamma-ray regime with a special attention to
the energy regime for H.E.S.S. Preliminary longitude-latitude profile and
residuum maps of inverse Compton emissions are shown for axisymmetric and
four-arm spiral Galaxy models.
[9]
oai:arXiv.org:1701.07284 [pdf] - 1535006
MHD Models of Gamma-ray Emission in WR 11
Submitted: 2017-01-25
Recent reports claiming tentative association of the massive star binary
system gamma^2 Velorum (WR 11) with a high-energy gamma-ray source observed by
Fermi-LAT contrast the so-far exclusive role of Eta Carinae as the hitherto
only detected gamma-ray emitter in the source class of particle-accelerating
colliding-wind binary systems. We aim to shed light on this claim of
association by providing dedicated model predictions for the nonthermal photon
emission spectrum of WR 11. We use three-dimensional magneto-hydrodynamic
modeling to trace the structure and conditions of the wind-collision region of
WR 11 throughout its 78.5 day orbit, including the important effect of
radiative braking in the stellar winds. A transport equation is then solved in
the wind-collision region to determine the population of relativistic electrons
and protons which are subsequently used to compute nonthermal photon emission
components. We find that - if WR 11 be indeed confirmed as the responsible
object for the observed gamma-ray emission - its radiation will unavoidably be
of hadronic origin owing to the strong radiation fields in the binary system
which inhibit the acceleration of electrons to energies suffciently high for
observable inverse Compton radiation. Different conditions in wind-collision
region near the apastron and periastron configuration lead to significant
variability on orbital time scales. The bulk of the hadronic gamma-ray emission
originates at a 400 solar radii wide region at the apex.
[11]
oai:arXiv.org:1609.01130 [pdf] - 1510361
Colliding-wind Binaries with strong magnetic fields
Submitted: 2016-09-05
The dynamics of colliding wind binary systems and conditions for efficient
particle acceleration therein have attracted multiple numerical studies in the
recent years. These numerical models seek an explanation of the thermal and
non-thermal emission of these systems as seen by observations. In the
non-thermal regime, radio and X-ray emission is observed for several of these
colliding-wind binaries, while gamma-ray emission has so far only been found in
$\eta$ Carinae and possibly in WR 11. Energetic electrons are deemed
responsible for a large fraction of the observed high-energy photons in these
systems. Only in the gamma-ray regime there might be, depending on the
properties of the stars, a significant contribution of emission from neutral
pion decay. Thus, studying the emission from colliding-wind binaries requires
detailed models of the acceleration and propagation of energetic electrons.
This in turn requires a detailed understanding of the magnetic field, which
will not only affect the energy losses of the electrons but in case of
synchrotron emission also the directional dependence of the emissivity. In this
study we investigate magnetohydrodynamic simulations of different colliding
wind binary systems with magnetic fields that are strong enough to have a
significant effect on the winds. Such strong fields require a detailed
treatment of the near-star wind acceleration zone. We show the implementation
of such simulations and discuss results that demonstrate the effect of the
magnetic field on the structure of the wind collision region.
[12]
oai:arXiv.org:1510.02580 [pdf] - 1290860
Galactic cosmic ray propagation models using Picard
Submitted: 2015-10-09
We present results obtained from our newly developed Galactic cosmic-ray
transport code PICARD, that solves the cosmic-ray transport equation. This code
allows for the computation of cosmic-ray spectra and the resulting gamma-ray
emission. Relying on contemporary numerical solvers allows for efficient
computation of models with deca-parsec resolution. PICARD can handle locally
anisotropic spatial diffusion acknowledging a full diffusion tensor. We used
this framework to investigate the transition from axisymmetric to spiral-arm
cosmic-ray source distributions. Wherever possible we compare model predictions
with constraining observables in cosmic-ray astrophysics.
[14]
oai:arXiv.org:1504.08249 [pdf] - 1095710
Propagation in 3D spiral-arm cosmic-ray source distribution models and
secondary particle production using PICARD
Submitted: 2015-04-30
We study the impact of possible spiral-arm distributions of Galactic
cosmic-ray sources on the flux of various cosmic-ray nuclei throughout our
Galaxy. We investigate model cosmic-ray spectra at the nominal position of the
sun and at different positions within the Galaxy. The modelling is performed
using the recently introduced numerical cosmic ray propagation code
\textsc{Picard}. Assuming non-axisymmetric cosmic ray source distributions
yields new insights on the behaviour of primary versus secondary nuclei.
We find that primary cosmic rays are more strongly confined to the vicinity
of the sources, while the distribution of secondary cosmic rays is much more
homogeneous compared to the primaries. This leads to stronger spatial variation
in secondary to primary ratios when compared to axisymmetric source
distribution models. A good fit to the cosmic-ray data at Earth can be
accomplished in different spiral-arm models, although leading to decisively
different spatial distributions of the cosmic-ray flux. This results in very
different cosmic ray anisotropies, where even a good fit to the data becomes
possible. Consequently, we advocate directions to seek best fit propagation
parameters that take into account the higher complexity introduced by the
spiral-arm structure on the cosmic-ray distribution. We specifically
investigate whether the flux at Earth is representative for a large fraction of
the Galaxy. The variance among possible spiral-arm models allows us to quantify
the spatial variation of the cosmic-ray flux within the Galaxy in presence of
non-axisymmetric source distributions.
[15]
oai:arXiv.org:1504.01858 [pdf] - 1095671
Implementing turbulence transport in the CRONOS framework and
application to the propagation of CMEs
Submitted: 2015-04-08
We present the implementation of turbulence transport equations in addition
to the Reynolds-averaged MHD equations within the Cronos framework. The model
is validated by comparisons with earlier findings before it is extended to be
applicable to regions in the solar wind that are not highly super-Alfv\'enic.
We find that the respective additional terms result in absolute normalized
cross-helicity to decline more slowly, while a proper implementation of the
mixing terms can even lead to increased cross-helicities in the inner
heliosphere. The model extension allows to place the inner boundary of the
simulations closer to the Sun, where we choose its location at 0.1 AU for
future application to the Wang-Sheeley-Arge model. Here, we concentrate on
effects on the turbulence evolution for transient events by injecting a coronal
mass ejection (CME). We find that the steep gradients and shocks associated
with these structures result in enhanced turbulence levels and reduced
cross-helicity. Our results can now be used straightforwardly for studying the
transport of charged energetic particles, where the elements of the diffusion
tensor can now benefit from the self-consistently computed solar wind
turbulence. Furthermore, we find that there is no strong back-reaction of the
turbulence on the large-scale flow so that CME studies concentrating on the
latter need not be extended to include turbulence transport effects.
[16]
oai:arXiv.org:1410.5266 [pdf] - 1222659
Spiral Arms as Cosmic Ray Source Distributions
Submitted: 2014-10-20
There is evidence that the distribution of suspected cosmic ray sources are
associated with the spiral arm structure of galaxies. It is yet not clearly
understood what effect such a cosmic ray source distribution has on the
particle transport in our Galaxy. We use the PICARD code to perform
high-resolution 3D simulations of electrons and protons in galactic propagation
scenarios that include four-arm and two-arm logarithmic spiral cosmic ray
source distributions with and without a central bar structure as well as the
spiral arm configuration of the NE2001 model for the distribution of free
electrons in the Milky Way. Results of these simulation are compared to an
axisymmetric radial source distribution. Also, effects on the cosmic ray flux
and spectra due to different positions of the Earth relative to the spiral
structure are studied. We find that high energy electrons are strongly confined
to their sources and the obtained spectra largely depend on the Earth's
position relative to the spiral arms. Similar finding have been obtained for
low energy protons and electrons albeit at smaller magnitude. We find that even
fractional contributions of a spiral arm component to the total cosmic ray
source distribution influences the spectra on the Earth. This is apparent when
compared to an axisymmetric radial source distribution as well as with respect
to the Earth's position relative to the spiral arm structure. We demonstrate
that the presence of a Galactic bar manifests itself as an overall excess of
low energy electrons at the Earth. Using a spiral arm geometry as a cosmic ray
source distributions offers a genuine new quality of modelling and is used to
explain features in cosmic ray spectra at the Earth that are else-wise
attributed to other propagation effects. We show that realistic cosmic ray
propagation scenarios have to acknowledge non-axisymmetric source
distributions.[Abbreviated]
[17]
oai:arXiv.org:1406.0293 [pdf] - 830040
Cosmic Ray Transport in Heliospheric Magnetic Structures: I. Modeling
Background Solar Wind Using the CRONOS MHD Code
Submitted: 2014-06-02
The transport of energetic particles such as Cosmic Rays is governed by the
properties of the plasma being traversed. While these properties are rather
poorly known for galactic and interstellar plasmas due to the lack of in situ
measurements, the heliospheric plasma environment has been probed by spacecraft
for decades and provides a unique opportunity for testing transport theories.
Of particular interest for the 3D heliospheric transport of energetic particles
are structures such as corotating interaction regions (CIRs), which, due to
strongly enhanced magnetic field strengths, turbulence, and associated shocks,
can act as diffusion barriers on the one hand, but also as accelerators of low
energy CRs on the other hand as well. In a two-fold series of papers we
investigate these effects by modeling inner-heliospheric solar wind conditions
with a numerical magnetohydrodynamic (MHD) setup (this paper), which will serve
as an input to a transport code employing a stochastic differential equation
(SDE) approach (second paper). In this first paper we present results from 3D
MHD simulations with our code CRONOS: for validation purposes we use analytic
boundary conditions and compare with similar work by Pizzo. For a more
realistic modeling of solar wind conditions, boundary conditions derived from
synoptic magnetograms via the Wang-Sheeley-Arge (WSA) model are utilized, where
the potential field modeling is performed with a finite-difference approach
(FDIPS) in contrast to the traditional spherical harmonics expansion often
utilized in the WSA model. Our results are validated by comparing with
multi-spacecraft data for ecliptical (STEREO-A/B) and out-of-ecliptic (Ulysses)
regions.
[18]
oai:arXiv.org:1405.6868 [pdf] - 1209760
Simulating three-dimensional nonthermal high-energy photon emission in
colliding-wind binaries
Submitted: 2014-05-27
Massive stars in binary systems have long been regarded as potential sources
of high-energy gamma rays.The emission is principally thought to arise in the
region where the stellar winds collide and accelerate relativistic particles
which ubsequently emit gamma rays. On the basis of a three-dimensional
distribution function of high-energy particles in the wind collision region -
as obtained by a numerical hydrodynamics and particle transport model - we
present the computation of the three-dimensional nonthermal photon emission for
a given line of sight. Anisotropic inverse Compton emission is modelled using
the target radiation field of both stars. Photons from relativistic
bremsstrahlung and neutral pion decay are computed on the basis of local wind
plasma densities. We also consider photon photon opacity effects due to the
dense radiation fields of the stars. Results are shown for different stellar
separations of a given binary system comprising of a B star and a Wolf-Rayet
star. The influence of orbital orientation with respect to the line of sight is
also studied by using different orbital viewing angles. For the chosen
electron-proton injection ratio of 0.01, we present the ensuing photon emission
in terms of two-dimensional projections maps, spectral energy distributions and
integrated photon flux values in various energy bands. Here, we find a
transition from hadron-dominated to lepton-dominated high-energy emission with
increasing stellar separations. In addition, we confirm findings from previous
analytic modeling that the spectral energy distribution varies significantly
with orbital orientation.
[19]
oai:arXiv.org:1401.4035 [pdf] - 1202792
PICARD: A novel code for the Galactic Cosmic Ray propagation problem
Submitted: 2014-01-16
In this manuscript we present a new approach for the numerical solution of
the Galactic Cosmic Ray propagation problem. We introduce a method using
advanced contemporary numerical algorithms while retaining the general
complexity of other established codes. In this paper we present the underlying
numerical scheme in conjunction with tests showing the correctness of the
scheme. Finally we show the solution of a first example propagation problem
using the new code to show its applicability to Galactic Cosmic Ray
propagation.
[20]
oai:arXiv.org:1401.1323 [pdf] - 779638
High-energy particle transport in 3D hydrodynamic models of
colliding-wind binaries
Submitted: 2014-01-07
Massive stars in binary systems (as WR140, WR147 or $\eta$ Carinae) have long
been regarded as potential sources of high-energy $\gamma$-rays. The emission
is thought to arise in the region where the stellar winds collide and produce
relativistic particles which subsequently might be able to emit $\gamma$-rays.
Detailed numerical hydrodynamic simulations have already offered insight in the
complex dynamics of the wind collision region (WCR), while independent
analytical studies, albeit with simplified descriptions of the WCR, have shed
light on the spectra of charged particles. In this paper, we describe a
combination of these two approaches. We present a 3D-hydrodynamical model for
colliding stellar winds and compute spectral energy distributions of
relativistic particles for the resulting structure of the WCR. The hydrodynamic
part of our model incorporates the line-driven acceleration of the winds,
gravity, orbital motion and the radiative cooling of the shocked plasma. In our
treatment of charged particles we consider diffusive shock acceleration in the
WCR and the subsequent cooling via inverse Compton losses (including
Klein-Nishina effects), bremsstrahlung, collisions and other energy loss
mechanisms.
[21]
oai:arXiv.org:1308.2829 [pdf] - 706019
A new 3D transport and radiation code for galactic cosmic rays
Submitted: 2013-08-13
We show the necessity for a new approach towards comprehensive and consistent
simulations of the propagation of galactic cosmic rays. Our developments are
optimised for addressing the spatially 3-dimensional inhomogeneous diffusion
problem and utilise contemporary numerical methods. We aim to address the
transport problem in a full 3-dimensional environment. For that, we test the
transition from 2D to 3D simulation results within an existing propagation
code. We present sub-kpc scale simulations that allow the investigation of
small- scale structures regarding different model conditions such as variety
regarding non-axisymmetric cosmic ray source distributions. These results are
discussed critically and motivate our development of a new transport code for
galactic cosmic rays. The capabilities of this code are outlined.
[22]
oai:arXiv.org:1303.4260 [pdf] - 640575
MHD Simulation of the Inner-Heliospheric Magnetic Field
Submitted: 2013-03-18
Maps of the radial magnetic field at a heliocentric distance of ten solar
radii are used as boundary conditions in the MHD code CRONOS to simulate a 3D
inner-heliospheric solar wind emanating from the rotating Sun out to 1 AU. The
input data for the magnetic field are the result of solar surface flux
transport modelling using observational data of sunspot groups coupled with a
current sheet source surface model. Amongst several advancements, this allows
for higher angular resolution than that of comparable observational data from
synoptic magnetograms. The required initial conditions for the other MHD
quantities are obtained following an empirical approach using an inverse
relation between flux tube expansion and radial solar wind speed. The
computations are performed for representative solar minimum and maximum
conditions, and the corresponding state of the solar wind up to the Earths
orbit is obtained. After a successful comparison of the latter with
observational data, they can be used to drive outer-heliospheric models.
[23]
oai:arXiv.org:1203.0933 [pdf] - 1117062
Diffusion of energetic particles in turbulent MHD plasmas
Submitted: 2012-03-05
In this paper we investigate the transport of energetic particles in
turbulent plasmas. A numerical approach is used to simulate the effect of the
background plasma on the motion of energetic protons. The background plasma is
in a dynamically turbulent state found from numerical MHD simulations, where we
use parameters typical for the heliosphere. The implications for the transport
parameters (i.e. pitch-angle diffusion coefficients and mean free path) are
calculated and deviations from the quasi-linear theory are discussed.
[24]
oai:arXiv.org:1111.5497 [pdf] - 1091873
Global Structure of Magnetorotationally Turbulent Protoplanetary Discs
Submitted: 2011-11-23
The aim of the present paper is to investigate the spatial structure of a
protoplanetary disc whose dynamics is governed by magnetorotational turbulence.
We perform a series of local 3D chemo-radiative MHD simulations located at
different radii of a disc which is twice as massive as the standard minimum
mass solar nebula of Hayashi (1981). The ionisation state of the disc is
calculated by including collisional ionisation, stellar X-rays, cosmic rays and
the decay of radionuclides as ionisation sources, and by solving a simplified
chemical network which includes the effect of the absorption of free charges by
{\mu}m-sized dust grains. In the region where the ionisation is too low to
assure good coupling between matter and magnetic fields, a non-turbulent
central "dead zone" forms, which ranges approximately from a distance of 2 AU
to 4 AU from the central star. The approach taken in the present work allows
for the first time to derive the global spatial structure of a protoplanetary
disc from a set of physically realistic numerical simulations.
[25]
oai:arXiv.org:1109.1930 [pdf] - 1083952
Investigation of the recombination of the retarded shell of "born-again"
CSPNe by time-dependent radiative transfer models
Submitted: 2011-09-09
A standard planetary nebula stays more than 10 000 years in the state of a
photoionized nebula. As long as the timescales of the most important ionizing
processes are much smaller, the ionization state can be characterized by a
static photoionization model and simulated with codes like CLOUDY (Ferland et
al. 1998). When the star exhibits a late Helium flash, however, its ionizing
flux stops within a very short period. The star then re-appears from itsopaque
shell after a few years (or centuries) as a cold giant star without any hard
ionizing photons. Describing the physics of such behavior requires a fully
time-dependent radiative transfer model. Pollacco (1999), Kerber et al. (1999)
and Lechner & Kimeswenger (2004) used data of the old nebulae around V605 Aql
and V4334 Sgr to derive a model of the pre-outburst state of the CSPN in a
static model. Their argument was the long recombination time scale for such
thin media. With regard to these models Schoenberner (2008) critically raised
the question whether a significant change in the ionization state (and thus the
spectrum) has to be expected after a time of up to 80 years, and whether static
models are applicable at all.
[26]
oai:arXiv.org:1007.4747 [pdf] - 1033942
Vertical Structure and Turbulent Saturation Level in Fully Radiative
Protoplanetary Disc Models
Submitted: 2010-07-27
We investigate a massive ($\varSigma \sim 10000 g cm^{-2}$ at 1 AU)
protoplanetary disc model by means of 3D radiation magnetohydrodynamics
simulations. The vertical structure of the disc is determined self-consistently
by a balance between turbulent heating caused by the MRI and radiative cooling.
Concerning the vertical structure, two different regions can be distinguished:
A gas-pressure dominated, optically thick midplane region where most of the
dissipation takes place, and a magnetically dominated, optically thin corona
which is dominated by strong shocks. At the location of the photosphere, the
turbulence is supersonic ($M \sim 2$), which is consistent with previous
results obtained from the fitting of spectra of YSOs. It is known that the
turbulent saturation level in simulations of MRI-induced turbulence does depend
on numerical factors such as the numerical resolution and the box size.
However, by performing a suite of runs at different resolutions (using up to
64x128x512 grid cells) and with varying box sizes (with up to 16 pressure
scaleheights in the vertical direction), we find that both the saturation
levels and the heating rates show a clear trend to converge once a sufficient
resolution in the vertical direction has been achieved.
[27]
oai:arXiv.org:0901.4495 [pdf] - 315553
Local Turbulence Simulations for the Multiphase ISM
Submitted: 2009-01-28
In this paper we show results of numerical simulations for the turbulence in
the interstellar medium. These results were obtained using a Riemann
solver-free numerical scheme for high-Mach number hyperbolic equations. Here we
especially concentrate on the physical properties of the ISM. That is, we do
not present turbulence simulations trimmed to be applicable to the interstellar
medium. The simulations are rather based on physical estimates for the relevant
parameters of the interstellar gas.
Applying our code to simulate the turbulent plasma motion within a typical
interstellar molecular cloud, we investigate the influence of different
equations of state (isothermal and adiabatic) on the statistical properties of
the resulting turbulent structures. We find slightly different density power
spectra and dispersion maps, while both cases yield qualitatively similar
dissipative structures, and exhibit a departure from the classical Kolmogorov
case towards a scaling described by the She-Leveque model.
Solving the full energy equation with realistic heating/cooling terms
appropriate for the diffuse interstellar gas, we are able to reproduce a
realistic two-phase distribution of cold and warm plasma. When extracting maps
of polarised intensity from our simulation data, we find encouraging similarity
to actual observations. Finally, we compare the actual magnetic field strength
of our simulations to its value inferred from the rotation measure. We find
these to be systematically different by a factor of about 1.5, thus
highlighting the often underestimated influence of varying line-of-sight
particle densities on the magnetic field strength derived from observed
rotation measures.
[28]
oai:arXiv.org:0901.4029 [pdf] - 1001364
Growth of the MRI in Accretion Discs -- the Influence of Radiation
Transport
Submitted: 2009-01-26
In this paper we investigate the influence of radiative transport on the
growth of the magnetorotational instability (MRI) in accretion discs. The
analysis is performed by use of analytical and numerical means. We provide a
general dispersion relation together with the corresponding eigenfunctions
describing the growth rates of small disturbances on a homogeneous background
shear flow. The dispersion relation includes compressibility and radiative
effects in the flux-limited diffusion approximation. By introducing an
effective speed of sound, all the effects of radiation transport can be
subsumed into one single parameter. It can be shown that the growth rates of
the vertical modes -- which are the fastest growing ones -- are reduced by
radiative transport. For the case of non-vertical modes, the growth rates may
instead be enhanced. We quantify the effects of compressibility and radiative
diffusion on the growth rates for the gas-pressure dominated case. The
analytical discussion is supplemented by numerical simulations, which are also
used for a first investigation of the non-linear stage of the MRI in
gas-pressure dominated accretion discs with radiation transport included.
