Normalized to: Bandiera, R.
[1]
oai:arXiv.org:1906.02203 [pdf] - 1898124
The dust content of the Crab Nebula
De Looze, I.;
Barlow, M. J.;
Bandiera, R.;
Bevan, A.;
Bietenholz, M. F.;
Chawner, H.;
Gomez, H. L.;
Matsuura, M.;
Priestley, F.;
Wesson, R.
Submitted: 2019-06-05
We have modelled the near-infrared to radio images of the Crab Nebula with a
Bayesian SED model to simultaneously fit its synchrotron, interstellar and
supernova dust emission. We infer an interstellar dust extinction map with an
average $A_{\text{V}}$=1.08$\pm$0.38 mag, consistent with a small contribution
(<22%) to the Crab's overall infrared emission. The Crab's supernova dust mass
is estimated to be between 0.032 and 0.049 M$_{\odot}$ (for amorphous carbon
grains) with an average dust temperature $T_{\text{dust}}$=41$\pm$3K,
corresponding to a dust condensation efficiency of 8-12%. This revised dust
mass is up to an order of magnitude lower than some previous estimates, which
can be attributed to our different interstellar dust corrections, lower SPIRE
flux densities, and higher dust temperature than were used in previous studies.
The dust within the Crab is predominantly found in dense filaments south of the
pulsar, with an average V-band dust extinction of $A_{\text{V}}$=0.20-0.39 mag,
consistent with recent optical dust extinction studies. The modelled
synchrotron power-law spectrum is consistent with a radio spectral index
$\alpha_{\text{radio}}$=0.297$\pm$0.009 and an infrared spectral index
$\alpha_{\text{IR}}$=0.429$\pm$0.021. We have identified a millimetre excess
emission in the Crab's central regions, and argue that it most likely results
from two distinct populations of synchrotron emitting particles. We conclude
that the Crab's efficient dust condensation (8-12%) provides further evidence
for a scenario where supernovae can provide substantial contributions to the
interstellar dust budgets in galaxies.
[2]
oai:arXiv.org:1811.11551 [pdf] - 1795932
Interplay between Physics and Geometry in Balmer filaments: the Case of
SN 1006
Submitted: 2018-11-28
The analysis of Balmer-dominated emission in supernova remnants is
potentially a very powerful way to derive information on the shock structure,
on the physical conditions of the ambient medium and on the cosmic-ray
acceleration efficiency. However, the outcome of models developed in
plane-parallel geometry is usually not easily comparable with the data, since
they often come from regions with rather a complex geometry. We present here a
general scheme to disentangle physical and geometrical effects in the data
interpretation, which is especially powerful when the transition zone of the
shock is spatially resolved and the spectral resolution is high enough to allow
a detailed investigation of spatial changes of the line profile. We then apply
this technique to re-analyze very high quality data of a region along the
northwestern limb of the remnant of SN~1006. We show how some observed
features, previously interpreted only in terms of spatial variations of
physical quantities, naturally arise from geometrical effects. With these
effects under control, we derive new constraints on physical quantities in the
analyzed region, like the ambient density (in the range 0.03-$0.1{\,\rm
cm^{-3}}$), the upstream neutral fraction (more likely in the range 0.01-0.1),
the level of face-on surface brightness variations (with factors up to $\sim
3$) and the typical scale lengths related to such variations ($\ge 0.1{\,\rm
pc}$, corresponding to angular scales $\ge 10{\,\rm arcsec}$).
[4]
oai:arXiv.org:1706.02978 [pdf] - 1584474
Discovery of a Synchrotron Bubble Associated with PSR J1015-5719
Submitted: 2017-06-09
We report the discovery of a synchrotron nebula, G283.1-0.59, associated with
PSR J1015-5719. Radio observations using the Molonglo Observatory Synthesis
Telescope and the Australia Telescope Compact Array at 36, 16, 6, and 3 cm
reveal a complex morphology. The pulsar is embedded in the "head" of the nebula
with fan-shaped diffuse emission. This is connected to a circular bubble of 20"
radius and a collimated tail extending over 1'. Polarization measurements show
a highly ordered magnetic field in the nebula. It wraps around the edge of the
head and shows an azimuthal configuration near the pulsar, then switches
direction quasi-periodically near the bubble and in the tail. Together with the
flat radio spectrum observed, we suggest that this system is most plausibly a
pulsar wind nebula (PWN), with the head as a bow shock that has a low Mach
number and the bubble as a shell expanding in a dense environment. The bubble
could act as a magnetic bottle trapping the relativistic particles. A
comparison with other bow-shock PWNe with higher Mach numbers shows similar
structure and B-field geometry, implying that pulsar velocity may not be the
most critical factor in determining the properties of these systems.
We also derive analytic expressions for the projected standoff distance and
shape of an inclined bow shock. It is found that the projected distance is
always larger than the true distance in three dimensions. On the other hand,
the projected shape is not sensitive to the inclination after rescaling with
the projected standoff distance.
[5]
oai:arXiv.org:1705.05803 [pdf] - 1583479
Radio polarization maps of shell-type SNRs II. Sedov models with
evolution of turbulent magnetic field
Submitted: 2017-05-16
Polarized radio emission has been mapped with great detail in several
Galactic supernova remnants (SNRs), but has not yet been exploited to the
extent it deserves. We have developed a method to model maps of the Stokes
parameters for shell-like SNRs during their Sedov evolution phase. At first,
3-dimensional structure of a SNR has been computed, by modeling the
distribution of the magnetohydrodynamic parameters and of the accelerated
particles. The generation and dissipation of the turbulent component of
magnetic field everywhere in SNR are also considered taking into account its
interaction with accelerated particles. Then, in order to model the emission,
we have used a generalization of the classical synchrotron theory, valid for
the case in which the magnetic field has ordered and disordered components.
Finally, 2-dimensional projected maps have been derived, for different
orientations of SNR and of interstellar magnetic field with respect to the
observer. An important effect to consider is the Faraday rotation of the
polarization planes inside the SNR interior. In this paper we present details
of the model, and describe general properties of the images.
[6]
oai:arXiv.org:1704.06546 [pdf] - 1582760
Modeling the effect of small-scale magnetic turbulence on the X-ray
properties of Pulsar Wind Nebulae
Submitted: 2017-04-21
Pulsar Wind Nebulae (PWNe) constitute an ideal astrophysical environment to
test our current understanding of relativistic plasma processes. It is well
known that magnetic fields play a crucial role in their dynamics and emission
properties. At present, one of the main issues concerns the level of magnetic
turbulence present in these systems, which in the absence of space resolved
X-ray polarization measures cannot be directly constrained. In this work we
investigate, for the first time using simulated synchrotron maps, the effect of
a small scale fluctuating component of the magnetic field on the emission
properties in X-ray. We illustrate how to include the effects of a turbulent
component in standard emission models for PWNe, and which consequences are
expected in terms of net emissivity and depolarization, showing that the X-ray
surface brightness maps can provide already some rough constraints. We then
apply our analysis to the Crab and Vela nebulae and, by comparing our model
with Chandra and Vela data, we found that the typical energies in the turbulent
component of the magnetic field are about 1.5 to 3 times the one in the ordered
field.
[7]
oai:arXiv.org:1703.09183 [pdf] - 1582153
From molecules to Young Stellar Clusters: the star formation cycle
across the M33 disk
Corbelli, Edvige;
Brain, Jonathan;
Bandiera, Rino;
Brouillet, Nathalie;
Combes, Françoise;
Druard, Clement;
Gratier, Pierre;
Mata, Jimmy;
Schuster, Karl;
Xilouris, Manolis;
Palla, Francesco
Submitted: 2017-03-27
To shed light on the time evolution of local star formation episodes in M33,
we study the association between 566 Giant Molecular Clouds (GMCs), identified
through the CO (J=2-1) IRAM-all-disk survey, and 630 Young Stellar Cluster
Candidates (YSCCs), selected via Spitzer-24~$\mu$m emission. The spatial
correlation between YSCCs and GMCs is extremely strong, with a typical
separation of 17~pc, less than half the CO(2--1) beamsize, illustrating the
remarkable physical link between the two populations. GMCs and YSCCs follow the
HI filaments, except in the outermost regions where the survey finds fewer GMCs
than YSCCs, likely due to undetected, low CO-luminosity clouds. The GMCs have
masses between 2$\times 10^4$ and 2$\times 10^6$ M$_\odot$ and are classified
according to different cloud evolutionary stages: inactive clouds are 32$\%$ of
the total, classified clouds with embedded and exposed star formation are
16$\%$ and 52$\%$ of the total respectively. Across the regular southern spiral
arm, inactive clouds are preferentially located in the inner part of the arm,
possibly suggesting a triggering of star formation as the cloud crosses the
arm. Some YSCCs are embedded star-forming sites while the majority have
GALEX-UV and H$\alpha$ counterparts with estimated cluster masses and ages. The
distribution of the non-embedded YSCC ages peaks around 5~Myrs with only a few
being as old as 8--10~Myrs. These age estimates together with the number of
GMCs in the various evolutionary stages lead us to conclude that 14~Myrs is a
typical lifetime of a GMC in M33, prior to cloud dispersal. The inactive and
embedded phases are short, lasting about 4 and 2~Myrs respectively. This
underlines that embedded YSCCs rapidly break out from the clouds and become
partially visible in H$\alpha$ or UV long before cloud dispersal.
[8]
oai:arXiv.org:1610.05217 [pdf] - 1580415
Unveiling the Magnetic Structure of VHE SNRs/PWNe with XIPE, the X-ray
Imaging-Polarimetry Satellite
Submitted: 2016-10-17, last modified: 2016-10-18
The dynamics, energetics and evolution of pulsar wind nebulae (PWNe) and
supernova remnants (SNRs), are strongly affected by their magnetic field
strength and distribution. They are usually strong, extended, sources of
non-thermal X-ray radiation, producing intrinsically polarised radiation. The
energetic wind around pulsars produces a highly-magnetised, structured flow,
often displaying a jet and a torus and different features (i.e. wisps, knots).
This magnetic-dominant wind evolves as it moves away from the pulsar
magnetosphere to the surrounding large-scale nebula, becoming kinetic-dominant.
Basic aspects such how this conversion is produced, or how the jets and torus
are formed, as well as the level of turbulence in the nebula are still unknown.
Likewise, the processes ruling the acceleration of particles in shell-like SNRs
up to 1e15 eV, including the amplification of the magnetic field, are not clear
yet. Imaging polarimetry in this regard is crucial to localise the regions of
shock acceleration and to measure the strength and the orientation of the
magnetic field at these emission sites. X-ray polarimetry with the X-ray
Imaging Polarimetry Explorer (XIPE) will allow the understanding of the
magnetic field structure and intensity on different regions in SNRs and PWNe,
helping to unveil long-standing questions such as i.e. acceleration of cosmic
rays in SNRs or magnetic-to-kinetic energy transfer. SNRs and PWNe also
represent the largest population of Galactic very-high energy gamma-ray
sources, therefore the study of their magnetic distribution with XIPE will
provide fundamental ingredients on the investigation of those sources at very
high energies. We will discuss the physics case related to SNRs and PWNe and
the expectations of the XIPE observations of some of the most prominent SNRs
and PWNe.
[10]
oai:arXiv.org:1603.01199 [pdf] - 1379031
Radio polarization maps of shell-type SNRs I. Effects of a random
magnetic field component, and thin-shell models
Submitted: 2016-03-03
The maps of intensity and polarization of the radio synchrotron emission from
shell-type supernova remnants (SNRs) contain a considerable amount of
information, although of not easy interpretation. With the aim of deriving
constraints on the 3-D spatial distribution of the emissivity, as well as on
the structure of both ordered and random magnetic fields (MFs), we present here
a scheme to model maps of the emission and polarization in SNRs.
We first generalize the classical treatment of the synchrotron emission to
the case in which the MF is composed by an ordered MF plus an isotropic random
component, with arbitrary relative strengths. In the case of a power-law
particle energy distribution, we derive analytic formulae that formally
resemble those for the classical case. We also treat the case of a shock
compression of a fully random upstream field and we predict that the
polarization fraction in this case should be higher than typically measured in
SNRs. We implement the above treatment into a code, which simulates the
observed polarized emission of an emitting shell, taking into account also the
effect of the internal Faraday rotation.
Finally, we show simulated maps for different orientations with respect to
the observer, levels of the turbulent MF component, Faraday rotation levels,
distributions of the emissivity (either barrel-shaped or limited to polar
caps), and geometries for the ordered MF component (either tangential to the
shell, or radial). Their analysis allows us to outline properties useful for
the interpretation of radio intensity and polarization maps.
[12]
oai:arXiv.org:1311.1300 [pdf] - 1180521
Cosmic Ray acceleration and Balmer emission from RCW 86 (G315.4-2.3)
Submitted: 2013-11-06, last modified: 2014-02-12
Context. Observation of Balmer lines from the region around the forward shock
of supernova remnants (SNR) may provide valuable information on the shock
dynamics and the efficiency of particle acceleration at the shock.
Aims. We calculated the Balmer line emission and the shape of the broad
Balmer line for parameter values suitable for SNR RCW 86 (G315.4-2.3) as a
function of the cosmic-ray (CR) acceleration efficiency and of the level of
thermal equilibration between electrons and protons behind the shock. This
calculation aims at using the width of the broad Balmer-line emission to infer
the CR acceleration efficiency in this remnant.
Methods. We used the recently developed nonlinear theory of diffusive
shock-acceleration in the presence of neutrals. The semianalytical approach we
developed includes a description of magnetic field amplification as due to
resonant streaming instability, the dynamical reaction of accelerated particles
and the turbulent magnetic field on the shock, and all channels of interaction
between neutral hydrogen atoms and background ions that are relevant for the
shock dynamics.
Results. We derive the CR acceleration efficiency in the SNR RCW 86 from the
Balmer emission. Since our calculation used recent measurements of the shock
proper motion, the results depend on the assumed distance to Earth. For a
distance of 2 kpc the measured width of the broad Balmer line is compatible
with the absence of CR acceleration. For a distance of 2.5 kpc, which is a
widely used value in current literature, a CR acceleration efficiency of 5-30%
is obtained, depending upon the electron-ion equilibration and the ionization
fraction upstream of the shock. By combining information on Balmer emission
with the measured value of the downstream electron temperature, we constrain
the CR acceleration efficiency to be ~20%.
[13]
oai:arXiv.org:1311.1746 [pdf] - 1180572
Modeling statistical properties of the X-ray emission from aged Pulsar
Wind Nebulae
Submitted: 2013-11-07
The number of known PWNe has recently increased considerably, and the
majority of them are now middle-age objects. Recent studies have shown a clear
correlation of both X-ray luminosity and size with the PWN age, but fail in
providing a thorough explanation of the observed trends. Here I propose a
different approach to these effects, based on the hypothesis that the observed
trends do not simply reproduce the evolution of a "typical" PWN, but are a
combined effect of PWNe evolving under different ambient conditions, the
leading parameter being the ambient medium density. Using a simple analytic
approach, I show that most middle-age PWNe are more likely observable during
the reverberation phase, and I succeed in reproducing trends consistent with
those observed, provided that the evolution of the X-ray emitting electrons
keeps adiabatic over the whole reverberation phase. As a direct consequence, I
show that the X-ray spectra of older PWNe should be harder: also this is
consistent with observations.
[14]
oai:arXiv.org:1310.8496 [pdf] - 1180388
On the MHD modeling of the Crab Nebula radio emission
Submitted: 2013-10-31
In recent years, it has become a well-established paradigm that many aspects
of the physics of Pulsar Wind Nebulae (PWNe) can be fully accounted for within
a relativistic MHD description. Numerical simulations have proven extremely
successful in reproducing the X-ray morphology of the Crab Nebula, down to very
fine detail. Radio emission, instead, is currently one of the most obscure
aspects of the physics of these objects, and one that holds important
information about pulsar properties and their role as antimatter factories.
Here we address the question of radio emission morphology and integrated
spectrum from the Crab Nebula, by using for the first time an axisymmetric
dynamical model with parameters chosen to best reproduce its X-ray morphology.
Based on our findings we discuss constraints on the origin of the radio
emitting particles.
[15]
oai:arXiv.org:1306.6762 [pdf] - 1172366
Cosmic Ray acceleration and Balmer emission from SNR 0509-67.5
Submitted: 2013-06-28, last modified: 2013-07-31
Context: Observation of Balmer lines from the region around the forward shock
of supernova remnants may provide precious information on the shock dynamics
and on the efficiency of particle acceleration at the shock.
Aims: We calculate the Balmer line emission and the shape of the broad Balmer
line for parameter values suitable for SNR 0509-67.5, as a function of the
cosmic ray acceleration efficiency and of the level of thermal equilibration
between electrons and protons behind the shock. This calculation aims at using
the width of the broad Balmer line emission to infer the cosmic ray
acceleration efficiency in this remnant.
Methods: We use the recently developed non-linear theory of diffusive shock
acceleration in the presence of neutrals. The semi-analytical approach that we
developed includes a description of magnetic field amplification as due to
resonant streaming instability, the dynamical reaction of both accelerated
particles and turbulent magnetic field on the shock, and all channels of
interaction between neutral atoms and background plasma that change the shock
dynamics.
Results: We achieve a quantitative assessment of the CR acceleration
efficiency in SNR 0509-67.5 as a function of the shock velocity and different
levels of electron-proton thermalization in the shock region. If the shock
moves faster than ~4500 km/s, one can conclude that particle acceleration must
be taking place with efficiency of several tens of percent. For lower shock
velocity the evidence for particle acceleration becomes less clear because of
the uncertainty in the electron-ion equilibration downstream. We also discuss
the role of future measurements of the narrow Balmer line.
[16]
oai:arXiv.org:1306.6454 [pdf] - 1172329
Broad Balmer line emission and cosmic ray acceleration efficiency in
supernova remnant shocks
Submitted: 2013-06-27, last modified: 2013-07-31
Balmer emission may be a powerful diagnostic tool to test the paradigm of
cosmic ray (CR) acceleration in young supernova remnant (SNR) shocks. The width
of the broad Balmer line is a direct indicator of the downstream plasma
temperature. In case of efficient particle acceleration an appreciable fraction
of the total kinetic energy of the plasma is channeled into CRs, therefore the
downstream temperature decreases and so does the broad Balmer line width. This
width also depends on the level of thermal equilibration between ions and
neutral hydrogen atoms in the downstream. Since in general in young SNR shocks
only a few charge exchange (CE) reactions occur before ionization,
equilibration between ions and neutrals is not reached, and a kinetic
description of the neutrals is required in order to properly compute Balmer
emission.
We provide a method for the calculation of Balmer emission using a
self-consistent description of the shock structure in the presence of neutrals
and CRs. We use a recently developed semi-analytical approach, where neutral
particles, ionized plasma, accelerated particles and magnetic fields are all
coupled together through the mass, momentum and energy flux conservation
equations. The distribution of neutrals is obtained from the full Boltzmann
equation in velocity space, coupled to Maxwellian ions through ionization and
CE processes. The computation is also improved with respect to previous work
thanks to a better approximation for the atomic interaction rates. We find that
for shock speeds >2500km/s the distribution of broad neutrals never approaches
a Maxwellian and its moments differ from those of the ionized component. These
differences reflect into a smaller FWHM than predicted in previous
calculations, where thermalization was assumed. The method presented here
provides a realistic estimate of particle acceleration efficiency in Balmer
dominated shocks.
[18]
oai:arXiv.org:1211.6148 [pdf] - 1158064
Collisionless shocks in a partially ionized medium: III. Efficient
cosmic ray acceleration
Submitted: 2012-11-26, last modified: 2013-03-19
In this paper we present the first formulation of the theory of non-linear
particle acceleration in collisionless shocks in the presence of neutral
hydrogen in the acceleration region. The dynamical reaction of the accelerated
particles, the magnetic field amplification and the magnetic dynamical effects
on the shock are also included. The main new aspect consists however in
accounting for charge exchange and ionization of neutral hydrogen, which
profoundly change the structure of the shock, as discussed in our previous
work. This important dynamical effect of neutrals is mainly associated to the
so-called neutral return flux, namely the return of hot neutrals from the
downstream region to the upstream, where they deposit energy and momentum
through charge exchange and ionization. We also present the self-consistent
calculation of Balmer line emission from the shock region and discuss how to
use measurements of the anomalous width of the different components of the
Balmer line to infer the cosmic ray acceleration efficiency in supernova
remnants showing Balmer emission: the broad Balmer line, which is due to charge
exchange of hydrogen atoms with hot ions downstream of the shock, is shown to
become narrower as a result of the energy drainage into cosmic rays, while the
narrow Balmer line, due to charge exchange in the cosmic-ray--induced
precursor, is shown to become broader. In addition to these two well-known
components, the neutral return flux leads to the formation of a third component
with intermediate width: this too contains information on ongoing processes at
the shock.
[19]
oai:arXiv.org:1210.4296 [pdf] - 1152162
Collisionless shock in a partially ionized medium. II. Balmer emission
Submitted: 2012-10-16
Strong shocks propagating into a partially ionized medium are often
associated with optical Balmer lines. This emission is due to impact excitation
of neutral hydrogen by hot protons and electrons in the shocked gas. The
structure of such Balmer-dominated shocks has been computed in a previous paper
(Blasi et al. 2012), where the distribution function of neutral particles was
derived from the appropriate Boltzmann equation including coupling with ions
and electrons through charge exchange and ionization. This calculation showed
how the presence of neutrals can significantly modify the shock structure
through the formation of a "neutral-induced" precursor ahead of the shock. Here
we follow up on our previous work and investigate the properties of the
resulting Balmer emission, with the aim of using the observed radiation as a
diagnostic tool for shock parameters. Our main focus is on Supernova Remnant
shocks, and we find that, for typical parameters, the H{\alpha} emission
typically has a three-component spectral profile, where: 1) a narrow component
originates from upstream cold hydrogen atoms, 2) a broad component comes from
hydrogen atoms that have undergone charge exchange with shocked protons
downstream of the shock, and 3) an intermediate component is due to hydrogen
atoms that have undergone charge exchange with warm protons in the
neutral-induced precursor. The width of the intermediate line reflects the
temperature in the precursor, while the width of the narrow one is left
unaltered by the precursor. In addition, we show that the profiles of both the
intermediate and broad components generally depart from a thermal distribution.
Finally, we show that a significant amount of Balmer emission can be produced
in the precursor region if efficient electron heating takes place.
[20]
oai:arXiv.org:1202.3080 [pdf] - 1116630
Collisionless shocks in a partially ionized medium: I. Neutral return
flux and its effects on acceleration of test particles
Submitted: 2012-02-14, last modified: 2012-06-19
A collisionless shock may be strongly modified by the presence of neutral
atoms through the processes of charge exchange between ions and neutrals and
ionization of the latter. These two processes lead to exchange of energy and
momentum between charged and neutral particles both upstream and downstream of
the shock. In particular, neutrals that suffer a charge exchange downstream
with shock-heated ions generate high velocity neutrals that have a finite
probability of returning upstream. These neutrals might then deposit heat in
the upstream plasma through ionization and charge exchange, thereby reducing
the fluid Mach number. A consequence of this phenomenon, that we refer to as
the "neutral return flux", is a reduction of the shock compression factor and
the formation of a shock precursor upstream. The scale length of the precursor
is determined by the ionization and charge exchange interaction lengths of fast
neutrals moving towards upstream infinity. In the case of a shock propagating
in the interstellar medium, the effects of ion-neutral interactions are
especially important for shock velocities < 3000 km/s. Such propagation
velocities are common among shocks associated with supernova remnants, the
primary candidate sources for the acceleration of Galactic cosmic rays. We then
investigate the effects of the return flux of neutrals on the spectrum of
test-particles accelerated at the shock. We find that, for shocks slower than
~3000 km/s, the particle energy spectrum steepens appreciably with respect to
the naive expectation for a strong shock, namely E^-2.
[21]
oai:arXiv.org:1105.4373 [pdf] - 514130
Infrared imaging and polarimetric observations of the pulsar wind nebula
in SNR G21.5-0.9
Submitted: 2011-05-22, last modified: 2012-05-18
We present infrared observations of the supernova remnant G21.5-0.9 with the
Very Large Telescope, the Canada-France-Hawaii Telescope and the Spitzer Space
Telescope. Using the VLT/ISAAC camera equipped with a narrow-band [FeII] 1.64um
filter the entire pulsar wind nebula in SNR G21.5-0.9 was imaged. This led to
detection of iron line-emitting material in the shape of a broken ring-like
structure following the nebula's edge. The detected emission is
limb-brightened. We also detect the compact nebula surrounding PSR J1833-1034,
both through imaging with the CFHT/AOB-KIR instrument (K' band) and the IRAC
camera (all bands) and also through polarimetric observations performed with
VLT/ISAAC (Ks band). The emission from the compact nebula is highly polarised
with an average value of the linear polarisation fraction $P_{L}^{avg} \simeq
0.47$, and the swing of the electric vector across the nebula can be observed.
The infrared spectrum of the compact nebula can be described as a power law of
index $\alpha_{IR} = 0.7 \pm 0.3$, and suggests that the spectrum flattens
between the infrared and X-ray bands.
[22]
oai:arXiv.org:1004.3515 [pdf] - 1026424
XMM-Newton and Suzaku detection of an X-ray emitting shell around the
pulsar wind nebula G54.1+0.3
Submitted: 2010-04-20
Recent X-ray observations have proved to be very effective in detecting
previously unknown supernova remnant shells around pulsar wind nebulae (PWNe),
and in these cases the characteristics of the shell provide further clues on
the evolutionary stage of the embedded PWN. However, it is not clear why some
PWNe are still "naked". We carried out an X-ray observational campaign targeted
at the PWN G54.1+0.3, the "close cousin" of the Crab, with the aim to detect
the associated SNR shell. We analyzed an XMM-Newton and Suzaku observations of
G54.1+0.3 and we model out the contribution of dust scattering halo. We
detected an intrinsic faint diffuse X-ray emission surrounding a hard spectrum,
which can be modeled either with a power-law (gamma= 2.9) or with a thermal
plasma model (kT=2.0 keV.). If the shell is thermal, we derive an explosion
energy E=0.5-1.6x10^51 erg, a pre-shock ISM density of 0.2 cm^-3 and an age of
about 2000 yr. Using these results in the MHD model of PWN-SNR evolution, we
obtain an excellent agreement between the predicted and observed location of
the shell and PWN shock.
[23]
oai:arXiv.org:0911.0829 [pdf] - 1018225
A statistical approach to radio emission from shell-type SNRs. I. Basic
ideas, techniques, and first results
Submitted: 2009-11-04
Shell-type supernova remnants (SNRs) exhibit correlations between radio
surface brightness, SNR diameter, and ambient medium density. We investigate
these correlations, to extract useful information about the typical
evolutionary stage of radio SNRs, and to obtain insight into the origin of the
relativistic electrons and magnetic fields responsible for the radio emission.
We propose a scenario, according to which the observed correlations are the
combined effect of SNRs evolving in a wide range of ambient conditions, rather
than the evolutionary track of a "typical" SNR. We then develop a parametric
approach to interpret the statistical data, and apply it to the data sample
previously published by Berkhuijsen, as well as to a sample of SNRs in the
galaxy M33. We find that SNRs cease to emit effectively in radio at a stage
near the end of their Sedov evolution, and that models of synchrotron emission
with constant efficiencies in particle acceleration and magnetic field
amplification do not provide a close match to the data. We discuss the problem
of the cumulative distribution in size, showing that the slope of this
distribution does not relate to the expansion law of SNRs, as usually assumed,
but only to the ambient density distribution. This solves a long-standing
paradox: the almost linear cumulative distribution of SNRs led several authors
to conclude that these SNRs are still in free expansion, which also implies
very low ambient densities. Within this framework, we discuss the case of the
starburst galaxy M82. Statistical properties of SNR samples may be used to shed
light on both the physics of electron acceleration and the evolution of SNRs.
More precise results could be obtained by combining data of several surveys of
SNRs in nearby galaxies.
[24]
oai:arXiv.org:0809.2159 [pdf] - 16253
On the X-ray feature associated with the Guitar Nebula
Submitted: 2008-09-12
Context: A mysterious X-ray nebula, showing a remarkably linear geometry, was
recently discovered close to the Guitar Nebula, the bow-shock nebula associated
with B2224+65, which is the fastest pulsar known. The nature of this X-ray
feature is unknown, and even its association with pulsar B2224+65 is unclear.
Aims: We attempt to develop a self-consistent scenario to explain the complex
phenomenology of this object. Methods: We assume that the highest energy
electrons accelerated at the termination shock escape from the bow shock and
diffuse into the ambient medium, where they emit synchrotron X-rays. The linear
geometry should reflect the plane-parallel geometry of its ambient field.
Results: We estimate the Lorentz factor of the X-ray emitting electrons and the
strength of the magnetic field. The former (~10^8) is close to its maximum
possible value, while the latter, at ~45 uG, is higher than typical
interstellar values and must have been amplified in some way. The magnetic
field must also be turbulent to some degree to trap the electrons sufficiently
for synchrotron X-ray emission to occur effectively. We propose a
self-consistent scenario in which, by some streaming instability, the electrons
themselves generate a turbulent field in which they then diffuse. Some
numerical coincidences are explained, and tests are proposed to verify our
scenario. Conclusions: Electron leaking may be common in the majority of pulsar
bow-shock nebulae, even though the X-ray nebulosity in general is too diffuse
to be detectable.
[25]
oai:arXiv.org:astro-ph/0606072 [pdf] - 82512
Influence of Thermalisation on Electron Injection in Supernova Remnant
Shocks
Submitted: 2006-06-05
Within a test-particle description of the acceleration process in parallel
nonrelativistic shocks, we present an analytic treatment of the electron
injection. We estimate the velocity distribution of the injected electrons as
the product of the post-shock thermal distribution of electrons times the
probability for electrons with a given velocity to be accelerated; the
injection efficiency is then evaluated as the integral of this velocity
distribution. We estimate the probability of a particle to be injected as that
of going back to the upstream region at least once. This is the product of the
probability of returning to the shock from downstream times that of recrossing
the shock from downstream to upstream. The latter probability is expected to be
sensitive to details of the process of electron thermalisation within the
(collisionless) shock, a process that is poorly known. In order to include this
effect, for our treatment we use results of a numerical, fully kinetic study,
by Bykov & Uvarov (1999). According to them, the probability of recrossing
depends on physics of thermalisation through a single free parameter (Gamma),
which can be expressed as a function of the Mach number of the shock, of the
level of electron-ion equilibration, as well as of the spectrum of turbulence.
It becomes apparent, from our analysis, that the injection efficiency is
related to the post-shock electron temperature, and that it results from the
balance between two competing effects: the higher the electron temperature, the
higher the fraction of downstream electrons with enough velocity to return to
the shock and thus to be ready to cross the shock from downstream to upstream;
at the same time, however, the higher the turbulence, which would hinder the
crossing.
[26]
oai:arXiv.org:astro-ph/0507518 [pdf] - 74687
The nature of the X-ray halo of the plerion G21.5-0.9 unveiled by
XMM-Newton and Chandra
Submitted: 2005-07-21
The nature of the radio-quiet X-ray halo around the plerionic SNR G21.5-0.9
is under debate. On the basis of spatial and spectral analysis of a large
Chandra and XMM-Newton dataset of this source, we have developed a
self-consistent scenario which explains all the observational features. We
found that the halo is composed by diffuse extended emission due to dust
scattering of X-rays from the plerion, by a bright limb which traces particle
acceleration in the fast forward shock of the remnant, and by a bright spot
(the ``North Spur'') which may be a knot of ejecta in adiabatic expansion. By
applying a model of interaction between the PWN, the SNR and supernova
environment, we argue that G21.5-0.9 progenitor may be of Type IIP or Ib/Ic,
and that the remnant may be young (200-1000 yr).
[27]
oai:arXiv.org:astro-ph/0504340 [pdf] - 1468766
Can hidden correlations mimic a variable fine structure constant?
Submitted: 2005-04-15
Murphy et al. (2003a, MNRAS, 345, 609) claim to find evidence of cosmological
variations of the fine structure constant $\alpha$ in the spectra of
intervening QSO absorption line systems. We find that this result is affected
by systematic effects. The $\alpha$ values estimated in individual line systems
depend on the set of atomic transitions used and therefore the quoted
dependence on the cosmic age may reflect the fact that different sets of
transitions are used at different redshifts. A correlation between line shifts
and relative optical depths of the atomic transitions is also present. This
correlation is very tight for a high-redshift subsample and accounts for the
anomalous dispersion of the $\alpha$ values found by Murphy et al. (2003a) in
this subsample. The above correlations are consistent with a scenario in which
gravitational redshift, caused by the gravity of the dark halo, contributes to
the shift in frequency of individual components. Gravitational redshift causes
differential spectral shifts of the same order as magnitude of those measured.
In the presence of line misidentification, these shifts can be interpreted in
terms of a variable $\alpha$. In order to verify the gravitational redshift
hypothesis, a direct access to Murphy et al. (2003a) data, or to a large amount
of new high resolution data, is necessary.
[28]
oai:arXiv.org:astro-ph/0405276 [pdf] - 64796
Magnetic Rayleigh-Taylor instability for Pulsar Wind Nebulae in
expanding Supernova Remnants
Submitted: 2004-05-14
We present a numerical investigation of the development of Rayleigh-Taylor
instability at the interface between an expanding Pulsar Wind Nebula and its
surrounding Supernova Remnant. These systems have long been thought to be
naturally subject to this kind of instability, given their expansion behavior
and the density jump at the contact discontinuity. High resolution images of
the Crab Nebula at optical frequencies show the presence of a complex network
of line-emitting filaments protruding inside the synchrotron nebula. These
structures are interpreted as the observational evidence that Rayleigh-Taylor
instability is in fact at work. The development of this instability in the
regime appropriate to describe Supernova Remnant-Pulsar Wind Nebula systems is
non-trivial. The conditions at the interface are likely close to the stability
threshold, and the inclusion of the nebular magnetic field, which might play an
important role in stabilizing the system, is essential to the modeling. If
Rayleigh-Taylor features can grow efficiently a mixing layer in the outer
portion of the nebula might form where most of the supernova material is
confined. When a magnetic field close to equipartition is included we find that
the interface is stable, and that even a weaker magnetic field affects
substantially the growth and shape of the fingers.
[29]
oai:arXiv.org:astro-ph/0404328 [pdf] - 64236
The effects of spin-down on the structure and evolution of pulsar wind
nebulae
Submitted: 2004-04-16
We present high resolution spherically symmetric relativistic
magnetohydrodynamical simulations of the evolution of a pulsar wind nebula
inside the free expanding ejecta of the supernova progenitor. The evolution is
followed starting from a few years after the supernova explosion and up to an
age of the remnant of 1500 years. We consider different values of the pulsar
wind magnetization parameter and also different braking indices for the
spin-down process. We compare the numerical results with those derived through
an approximate semi-analytical approach that allows us to trace the time
evolution of the positions of both the pulsar wind termination shock and the
contact discontinuity between the nebula and the supernova ejecta. We also
discuss, whenever a comparison is possible, to what extent our numerical
results agree with former self-similar models, and how these models could be
adapted to take into account the temporal evolution of the system. The inferred
magnetization of the pulsar wind could be an order of magnitude lower than that
derived from time independent analytic models.
[30]
oai:arXiv.org:astro-ph/0402598 [pdf] - 63112
Analytic Solutions for the Evolution of Radiative Supernova Remnants
Submitted: 2004-02-25
We present the general analytic solution for the evolution of radiative
supernova remnants in a uniform interstellar medium, under thin-shell
approximation. This approximation is shown to be very accurate approach to this
task. For a given set of parameters, our solution closely matches the results
of numerical models, showing a transient in which the deceleration parameter
reaches a maximum value of 0.33, followed by a slow convergence to the
asymptotic value 2/7. Oort (1951) and McKee and Ostriker (1977) analytic
solutions are discussed, as special cases of the general solution we have
found.
[31]
oai:arXiv.org:astro-ph/0305398 [pdf] - 56871
The X-ray Halo of G21.5-0.9
Submitted: 2003-05-21
The emission of the plerion G21.5-0.9 appears more extended in X rays than in
radio. This is an unexpected result because it would imply that short-lived
X-ray electrons may reach distances even larger than radio electrons. Applying
an empirical relationship between dust scattering optical depth and
photoelectric column density, the measured column density leads to a large
optical depth at 1 keV, of about 1. Therefore we investigate the hypothesis
that the detected halo be an effect of dust scattering, re-analyzing an Cal/PV
XMM-Newton observation of G21.5-0.9 and critically examining it in terms of a
dust scattering model. We also present a spectral analysis of a prominent
extended feature in the northern sector of the halo.
[32]
oai:arXiv.org:astro-ph/0211205 [pdf] - 52940
BeppoSAX observation of the composite remnant G327.1-1.1
Submitted: 2002-11-11
We report an X-ray study of the composite supernova remnant G327.1-1.1, with
particular emphasis on its thermal emission. By virtue of a combined spatial
and spectral analysis, we have been able to model the X-ray emission of the
remnant as a sum of two components: a non-thermal component, due to the pulsar
nebula and the pulsar itself, and a thermal component, of which we have
analysed spectrum and morphology, after proper subtraction of the plerion. We
discuss three possible interpretations of the thermal emission of \src: pure
Sedov expansion, expansion through a inhomogeneous medium with evaporation from
ISM clouds, and radiative expansion. On the light of our new data and
interpretation, we have re-derived all the physical parameters of this SNR. In
the framework of Sedov or radiative expansion we derive a longer age than
previously estimated ($1.1\times 10^4$), thus not requiring a high velocity for
the pulsar.
[33]
oai:arXiv.org:astro-ph/0206091 [pdf] - 49722
A catalogue of damped Lyman alpha absorption systems and radio flux
densities of the background quasars
Submitted: 2002-06-05
We present a catalogue of the 322 damped Lyman alpha absorbers taken from the
literature. All damped Lyman alpha absorbers are included, with no selection on
redshift or quasar magnitude. Of these, 123 are candidates and await
confirmation using high resolution spectroscopy. For all 322 objects we
catalogue the radio properties of the background quasars, where known. Around
60 quasars have radio flux densities above 0.1 Jy and approximately half of
these have optical magnitudes brighter than V = 18. This compilation should
prove useful in several areas of extragalactic/cosmological research.
[34]
oai:arXiv.org:astro-ph/0204050 [pdf] - 48595
BeppoSAX observations of the young pulsar in the Kes 75 Supernova
Remnant
Submitted: 2002-04-03
We present the results of BeppoSAX observations of the young X-ray pulsar PSR
J1846-0258, recently discovered at the center of the composite supernova
remnant Kes 75. The pulsar (plus nebula) spectrum can be fitted by an absorbed
power law with photon index alpha_ph=2.16+/-0.15, N_H=(4.7+/-0.8) x 10^22
cm^-2, and unabsorbed flux ~3.9x10^-11 erg cm^-2 s^-1 (2-10 keV). By joining
two observations taken at an interval of two weeks we have been able to obtain
a precise measurement of the spin period (P=324.818968+/-0.000006 ms). This
value, when combined with previous measurements, cannot be fitted by a smooth
frequency evolution with a canonical braking index n=3. With the hypothesis of
no glitches and/or significant timing noise, the braking index would be
n=1.86+/-0.08 and, assuming a short initial period, the pulsar age would be
\~1700 years, closer to that of the supernova remnant than the simple estimate
\tau=P/2\pdot=723 years. Other likely possibilities involve the presence of
glitches and lead to a wide range of acceptable ages. For example, we obtain n
in the range 1.8-2.5, if a glitch occurred near MJD 51500, while for a glitch
between October 1993 and March 1999 we can only get a lower limit of n>1.89.
[35]
oai:arXiv.org:astro-ph/0112459 [pdf] - 46832
The Crab Nebula at 1.3 mm: evidence for a new synchrotron component
Submitted: 2001-12-19, last modified: 2002-03-04
We present the results of 1.3 mm observations of the Crab Nebula, performed
with the MPIfR bolometer arrays at the IRAM 30-m telescope. The maps obtained,
of unprecedented quality at these wavelengths, allow a direct comparison with
high-resolution radio maps. Although the spatial structure of the Crab Nebula
does not change much from radio to millimetre wavelengths, we have detected
significant spatial variations of the spectral index between 20 cm and 1.3 mm.
The main effect is a spectral flattening in the inner region, which can be
hardly explained just in terms of the evolution of a single population of
synchrotron emitting electrons. We propose instead that this is the result of
the emergence of a second synchrotron component, that we have tried to extract
from the data. Shape and size of this component resemble those of the Crab
Nebula in X rays. However, while the more compact structure of the Crab Nebula
in X rays is commonly regarded as an effect of synchrotron downgrading, it
cannot be explained why a similar structure is present also at mm wavelengths,
where the electron lifetimes far exceed the nebular age. Our data, combined
with published upper limits on spatial variations of the radio spectral index,
also imply a low-energy cutoff for the distribution of electrons responsible
for this additional synchrotron component. Although no model has been developed
so far to explain the details of this component, one may verify that the total
number of the electrons responsible for it is in agreement with what predicted
by the classical pulsar-wind models, which otherwise are known to fail in
accounting for the number of radio emitting electrons. We have also detected a
spectral steepening at mm wavelengths in some elongated regions, whose
positions match those of radio synchrotron filaments.
[36]
oai:arXiv.org:astro-ph/0111580 [pdf] - 46346
Implications of the synchrotron emission in plerionic nebulae: a second
component in the Crab Nebula
Submitted: 2001-11-30
High-quality spatially resolved spectra of synchrotron emission from plerions
represent a powerful diagnostic tool for the physical conditions in these
objects. A fundamental question is to what extent spatial variations of the
synchrotron spectrum originate from the evolution of a single population of
injected electrons, and when instead they imply the coexistence of different
populations of electrons. Here I shall discuss observations of the Crab Nebula
at millimetric wavelengths which, compared with a radio map, show evidence for
the emergence of a second synchrotron component, characterized by a flatter
spectrum and undetected in the radio range. This component is confined to the
inner part of the nebula, and cannot originate from synchrotron-dominated
evolution, also because in the Crab Nebula the electrons emitting at millimeter
wavelengths evolve almost adiabatically. Then its discovery indicates that
there are more populations of injected particles. Finally, the detection of a
low-frequency break in the spectra of synchrotron filaments can be interpreted
as an effect of higher magnetic fields in filaments.
[37]
oai:arXiv.org:astro-ph/0111299 [pdf] - 46065
The effects of an ionizing background on the HI column density
distribution in the local Universe
Submitted: 2001-11-15
Using data on the HI column density distribution in the local Universe,
f(N_{HI}), we show how to determine g(N_{H}), the distribution of the total gas
(HI+HII) column density. A simple power law fit to f(N_{HI}) fails due to
bendings in the distributions when N_{HI}<10^{20} cm^{-2} and H is no longer
fully neutral. If an ultraviolet background is responsible for the gas
ionization, and g(N_{H}) is proportional N_{H}^{-alpha}, we find the values of
alpha and of the intensity of the background radiation which are compatible
with the present data. These best fitting values, however, depend upon the
scaling law of the the gas volume densities with N_{H} and cannot be determined
unambiguously. We examine in detail two models: one in which the average gas
volume density decreases steadily with N_H, while in the other it stays
constant at low column densities. The former model leads to a steep power law
fit for g(N_{H}), with alpha of order 3.3 +- 0.4 and requires an ultraviolet
flux larger than what the QSOs alone produce at z=0. For the latter alpha is
1.5 +- 0.1 and a lower ionizing flux is required. The ambiguities about the
modelling and the resulting steep or shallow N_{H} distribution can be resolved
only if new 21-cm observations and QSOs Lyman limit absorbers searches will
provide more data in the HI-HII transition region at low redshifts. Using the
best fit obtained for higher redshift data we outline two possible scenarios
for the evolution of gaseous structures, compatible with the available data at
z=0.
[38]
oai:arXiv.org:astro-ph/0101335 [pdf] - 40473
A Comprehensive Statistical Analysis of the Gas Distribution in
Lyman-limit and Damped Lyman-alpha Absorption Systems
Submitted: 2001-01-18
In this paper we show how to use data on Lyman-limit and Damped Lyman-alpha
absorption systems to derive the hydrogen ionization fractions and the
distribution of the face-on total gas column density. We consider axially
symmetric, randomly oriented absorbers, ionized by an external background
radiation field in order to relate the face-on total gas distribution to that
of the neutral hydrogen observed along the line of sight. We devise a
statistical procedure based on the Maximum Likelihood criterion, that is able
to treat simultaneously data coming from different surveys and statistically
recovers the "true" column densities in the presence of large uncertainties:
this is especially important for Lyman-limit systems which leave an
unmeasurable residual flux at wavelengths shorter than the Lyman break. We make
use of simulated data to look for possible observational biases and extensively
test our procedure. For a large statistical sample of real data in the redshift
range [1.75,3.25] (collected from all published surveys) our Maximum Likelihood
procedure gives a power-law slope for the total hydrogen distribution of -2.7.
All together Lyman-limit systems therefore contain more gas than Damped
Lyman-alpha systems. Analysis of data at other redshifts shows that more
observations are needed to reach a compelling evidence for a cosmological
evolution of the slope of the gas distribution.
[39]
oai:arXiv.org:astro-ph/0011268 [pdf] - 39275
Sharp HI edges at high z: the gas distribution from Damped Lyman-alpha
to Lyman-limit absorption systems
Submitted: 2000-11-14
We derive the distribution of neutral and ionized gas in high redshift clouds
which are optically thick to hydrogen ionizing radiation, using published data
on Lyman-limit and Damped Lyman-alpha absorption systems in the redshift range
1.75 < z < 3.25. We assume that the distribution of the hydrogen total (HI+HII)
column density in the absorbers follows a power law K N_H^{-alpha}, whereas the
observed HI column density distribution deviates from a pure power law as a
result of ionization from a background radiation field. Comparison of the
models and observations give Maximum Likelihood solutions for the exponent
alpha and for X, the value of log(N_H/N_HI) when the Lyman-limit optical depth
is unity: alpha=2.7^{+1.0}_{-0.7} and X=2.75\pm0.35. X is much lower than what
would be obtained for a gaseous distribution in equilibrium under its own
gravity but the ratio of dark matter to gas density is not well constrained
being log(eta_0)=1.1\pm 0.8. An extrapolation of our derived power law
distribution towards systems of lower column density, the Lyman-alpha forest,
favours models with log(eta_0) < 1.1 and alpha=2.7-3.3. With alpha appreciably
larger than 2, Lyman-limit systems contain more gas than Damped Lyman-alpha
systems and Lyman-alpha forest clouds even more. Estimates of the cosmological
gas and dark matter density due to absorbers of different column density around
z=2.5 are also given.
[40]
oai:arXiv.org:astro-ph/9911163 [pdf] - 109291
A model for the spatial distribution of relativistic electrons in the
Crab Nebula
Submitted: 1999-11-10, last modified: 2000-05-23
We extend Pacini & Salvati's model for the synchrotron emission from plerions
to allow more realistic predictions of the surface brightness profiles of the
Crab Nebula at radio, optical and X-ray frequencies. Abandoning the assumption
of a uniform particle distribution, we assume that particles are injected into
the synchrotron nebula only in the vicinity of the pulsar. Their distribution
is then determined, to zeroth order, by MHD propagation in an azimuthal
magnetic field of spatially constant but time-dependent intensity. A good
agreement with observations across the electromagnetic spectrum is obtained if
the finite extension of the particle injection region and particle diffusion
with respect to the azimuthal field lines are taken into account.
[41]
oai:arXiv.org:astro-ph/9911416 [pdf] - 109544
Properties of the Non-Thermal Emission in Plerions
Submitted: 1999-11-22
The synchrotron emission observed in plerions depends on the characteristics
of the magnetic fields and relativistic electrons in these supernova remnants.
Therefore an analysis of the spectral and spatial properties of this emission,
combined with models for the evolution and structure of plerions, would allow
one to investigate the evolution of the synchrotron nebula, the structure of
the magnetic field, the distribution of the relativistic electrons, as well as
to put constraints on the history of the energy input from the associated
neutron star. The identification of a new class of plerionic remnants, with
spectral properties different from the Crab Nebula, has been proposed. The
spectra of these objects typically show a sharp spectral break at very low
frequencies (below 50 GHz), with a steep spectrum right beyond the break. In
order to model the properties of their emission, a non-standard evolution of
the pulsar output seems to be required. X-ray observations of the synchrotron
emission from the Crab Nebula are shown. They are compared with previous data,
and their implications on the structure of the Crab are discussed. Recent
millimetric data of this object are also presented. A spatially resolved
analysis, based on radio, millimetric and X-ray data, will be carried on also
for other plerions.
[42]
oai:arXiv.org:astro-ph/9604085 [pdf] - 94456
The Supernova Remnant G11.2-0.3 and its central Pulsar
Submitted: 1996-04-16
The plerion inside the composite Supernova Remnant G11.2-0.3 appears to be
dominated by the magnetic field to an extent unprecedented among well known
cases. We discuss its evolution as determined by a central pulsar and the
interaction with the surrounding thermal remnant, which in turn interacts with
the ambient medium. We find that a plausible scenario exists, where all the
observations can be reproduced with rather typical values for the parameters of
the system; we also obtain the most likely period for the still undetected
pulsar.
[43]
oai:arXiv.org:astro-ph/9305004 [pdf] - 90761
Modelling Non-Axisymmetric Bow Shocks
Submitted: 1993-05-06
I present an extension of the thin layer approximation to non-axisymmetric
bow shocks. By choosing a suitable set of curvilinear coordinates that matches
the geometry of a generally distorted bow shock surface, I derive the fluid
equations for the flow. Analytical expansions are given for the matter flow
near the bow shock stagnation point. Numerical solutions are also obtained, by
using a revised Runge-Kutta method. This numerical method results to be
computationally stable and fast. When two winds interact, axisymmetric bow
shocks are expected to form only if both winds are spherically symmetric. But
stellar winds often present strong anisotropies: for all these cases a
non-axisymmetric model of bow shock is required.
