Source author record

Marco Miceli

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

20 published item(s)

preprint2024arXiv

The impact of effective matter mixing based on three-dimensional hydrodynamical models on the molecule formation in the ejecta of SN 1987A

To investigate the impact of matter mixing on the formation of molecules in the ejecta of SN 1987A, time-dependent rate equations for chemical reactions are solved for one-zone and one-dimensional ejecta models of SN 1987A. The latter models are based on the one-dimensional profiles obtained by angle-averaging of the three-dimensional hydrodynamical models (Ono et al. 2020), which effectively reflect the 3D matter mixing; the impact is demonstrated, for the first time, based on three-dimensional hydrodynamical models. The distributions of initial seed atoms and radioactive $^{56}$Ni influenced by the mixing could affect the formation of molecules. By comparing the calculations for spherical cases and for several specified directions in the bipolar-like explosions in the three-dimensional hydrodynamical models, the impact is discussed. The decay of $^{56}$Ni, practically $^{56}$Co at later phases, could heat the gas and delay the molecule formation. Additionally, Compton electrons produced by the decay could ionize atoms and molecules and could destruct molecules. Several chemical reactions involved with ions such as H$^+$ and He$^+$ could also destruct molecules. The mixing of $^{56}$Ni plays a non-negligible role in both the formation and destruction of molecules through the processes above. The destructive processes of carbon monoxide and silicon monoxide due to the decay of $^{56}$Ni generally reduce the amounts. However, if the molecule formation is sufficiently delayed under a certain condition, the decay of $^{56}$Ni could locally increase the amounts through a sequence of reactions.

preprint2022arXiv

A spatially resolved study of hard X-ray emission in Kepler's SNR: indications of different regimes of particle acceleration

Synchrotron X-ray emission in young supernova remnants (SNRs) is a powerful diagnostic tool to study the population of high energy electrons accelerated at the shock front and the acceleration process. We performed a spatially resolved spectral analysis of NuSTAR and XMM-Newton observations of the young Kepler's SNR, aiming to study in detail its non-thermal emission in hard X-rays. We selected a set of regions all around the rim of the shell and extracted the corresponding spectra. The spectra were analyzed by adopting a model of synchrotron radiation in the loss-limited regime, to constrain the dependence of the cutoff energy of the synchrotron radiation on the shock velocity. We identify two different regimes of particle acceleration, characterized by different Bohm factors. In the north, where the shock interacts with a dense circumstellar medium (CSM), we found a more efficient acceleration than in the south, where the shock velocity is higher and there are no signs of shock interaction with dense CSM. Our results suggest an enhanced efficiency of the acceleration process in regions where the shock-CSM interaction generates an amplified and turbulent magnetic field. By combining hard X-ray spectra with radio and $γ-$ray observations of Kepler's SNR, we modelled the spectral energy distribution. In the light of our results we propose that the observed $γ-$ray emission is mainly hadronic, and originates in the northern part of the shell.

preprint2022arXiv

Additional evidence for a pulsar wind nebula in the heart of SN 1987A from multi-epoch X-ray data and MHD modeling

Since the day of its explosion, supernova (SN) 1987A has been closely monitored to study its evolution and to detect its central compact relic. In fact, the formation of a neutron star is strongly supported by the detection of neutrinos from the SN. However, besides the detection in the Atacama Large Millimeter/submillimeter Array (ALMA) data of a feature that is compatible with the emission arising from a proto-pulsar wind nebula (PWN), the only hint for the existence of such elusive compact object is provided by the detection of hard emission in NuSTAR data up to ~ 20 keV. We report on the simultaneous analysis of multi-epoch observations of SN 1987A performed with Chandra, XMM-Newton and NuSTAR. We also compare the observations with a state-of-the-art 3D magnetohydrodynamic (MHD) simulation of SN 1987A. A heavily absorbed power-law, consistent with the emission from a PWN embedded in the heart of SN 1987A, is needed to properly describe the high-energy part of the observed spectra. The spectral parameters of the best-fit power-law are in agreement with the previous estimate, and exclude diffusive shock acceleration as a possible mechanism responsible for the observed non-thermal emission. The information extracted from our analysis are used to infer the physical characteristics of the pulsar and the broad-band emission of its nebula, in agreement with the ALMA data. Analysis of the synthetic spectra also show that, in the near future, the main contribution to Fe K emission line will originate in the outermost shocked ejecta of SN 1987A.

preprint2022arXiv

The supernova remnant SN 1006 as a Galactic particle accelerator

The origin of cosmic rays is a pivotal open issue of high-energy astrophysics. Supernova remnants are strong candidates to be the Galactic factory of cosmic rays, their blast waves being powerful particle accelerators. However, supernova remnants can power the observed flux of cosmic rays only if they transfer a significant fraction of their kinetic energy to the accelerated particles, but conclusive evidence for such efficient acceleration is still lacking. In this scenario, the shock energy channeled to cosmic rays should induce a higher post-shock density than that predicted by standard shock conditions. Here we show this effect, and probe its dependence on the orientation of the ambient magnetic field, by analyzing deep X-ray observations of the Galactic remnant of SN 1006. By comparing our results with state-of-the-art models, we conclude that SN 1006 is an efficient source of cosmic rays and obtain an observational support for the quasi-parallel acceleration mechanism.

preprint2021arXiv

Indication of a Pulsar Wind Nebula in the hard X-ray emission from SN 1987A

Since the day of its explosion, SN 1987A (SN87A) was closely monitored with the aim to study its evolution and to detect its central compact relic. The detection of neutrinos from the supernova strongly supports the formation of a neutron star (NS). However, the constant and fruitless search for this object has led to different hypotheses on its nature. Up to date, the detection in the ALMA data of a feature somehow compatible with the emission arising from a proto Pulsar Wind Nebula (PWN) is the only hint of the existence of such elusive compact object. Here we tackle this 33-years old issue by analyzing archived observations of SN87A performed Chandra and NuSTAR in different years. We firmly detect nonthermal emission in the $10-20$ kev energy band, due to synchrotron radiation. The possible physical mechanism powering such emission is twofold: diffusive shock acceleration (DSA) or emission arising from an absorbed PWN. By relating a state-of-the-art magneto-hydrodynamic simulation of SN87A to the actual data, we reconstruct the absorption pattern of the PWN embedded in the remnant and surrounded by cold ejecta. We found that, even though the DSA scenario cannot be firmly excluded, the most likely scenario that well explains the data is the PWN emission.

preprint2020arXiv

Unveiling pure-metal ejecta X-ray emission in supernova remnants through their radiative recombination continuum

Spectral analysis of X-ray emission from ejecta in supernova remnants (SNRs) is hampered by the low spectral resolution of CCD cameras, which creates a degeneracy between the best-fit values of abundances and emission measure. The combined contribution of shocked ambient medium and ejecta to the X-ray emission complicates the determination of the ejecta mass and chemical composition, leading to big uncertainties in mass estimates and it can introduce a bias in the comparison between the observed ejecta composition and the yields predicted by explosive nucleosynthesis. We explore the capabilities of present and future spectral instruments with the aim of identifying a spectral feature which may allow us to discriminate between metal-rich and pure-metal plasmas in X-ray spectra of SNRs. We studied the behavior of the most common X-ray emission processes of an optically thin plasma in the high-abundance regime. We investigated spectral features of bremsstrahlung, radiative recombination continua (RRC) and line emission, by exploring a wide range of chemical abundances, temperatures and ionization parameters. We synthesized X-ray spectra from a 3D hydrodynamic (HD) simulation of Cas A, by using the response matrix from the Chandra/ACIS-S CCD detector and that of the XRISM/Resolve X-ray calorimeter. We found that a bright RRC shows up when the plasma is made of pure-metal ejecta, and a high spectral resolution is needed to identify this ejecta signature. We verified the applicability of our novel diagnostic tool and we propose a promising target for the future detection of such spectral feature: the southeastern Fe-rich clump of Cas A. While there is no way to unambiguously reveal pure-metal ejecta emission with CCD detectors, X-ray calorimeters will be able to pinpoint the presence of pure-metal RRC and to recover correctly absolute mass and the chemical composition of the ejecta.

preprint2019arXiv

Deep XMM-Newton Observations Reveal the Origin of Recombining Plasma in the Supernova Remnant W44

Recent X-ray studies revealed over-ionized recombining plasmas (RPs) in a dozen mixed-morphology (MM) supernova remnants (SNRs). However, the physical process of the over-ionization has not been fully understood yet. Here we report on spatially resolved spectroscopy of X-ray emission from W44, one of the over-ionized MM-SNRs, using XMM-Newton data from deep observations, aiming to clarify the physical origin of the over-ionization. We find that combination of low electron temperature and low recombination timescale is achieved in the region interacting with dense molecular clouds. Moreover, a clear anti-correlation between the electron temperature and the recombining timescale is obtained from each of the regions with and without the molecular clouds. The results are well explained if the plasma was over-ionized by rapid cooling through thermal conduction with the dense clouds hit by the blast wave of W44. Given that a few other over-ionized SNRs show evidence for adiabatic expansion as the major driver of the rapid cooling, our new result indicates that both processes can contribute to over-ionization in SNRs, with the dominant channel depending on the evolutionary stage.

preprint2019arXiv

Matter Mixing in Aspherical Core-collapse Supernovae: Three-dimensional Simulations with Single Star and Binary Merger Progenitor Models for SN 1987A

We perform three-dimensional hydrodynamic simulations of aspherical core-collapse supernovae focusing on the matter mixing in SN 1987A. The impacts of four progenitor (pre-supernova) models and parameterized aspherical explosions are investigated. The four pre-supernova models include a blue supergiant (BSG) model based on a slow merger scenario developed recently for the progenitor of SN 1987A (Urushibata et al. 2018). The others are a BSG model based on a single star evolution and two red supergiant (RSG) models. Among the investigated explosion (simulation) models, a model with the binary merger progenitor model and with an asymmetric bipolar-like explosion, which invokes a jetlike explosion, best reproduces constraints on the mass of high velocity $^{56}$Ni, as inferred from the observed [Fe II] line profiles. The advantage of the binary merger progenitor model for the matter mixing is the flat and less extended $ρ\,r^3$ profile of the C+O core and the helium layer, which may be characterized by the small helium core mass. From the best explosion model, the direction of the bipolar explosion axis (the strongest explosion direction), the neutron star (NS) kick velocity, and its direction are predicted. Other related implications and future prospects are also given.

preprint2016arXiv

Modeling the shock-cloud interaction in SN 1006: unveiling the origin of nonthermal X-ray and gamma-ray emission

The supernova remnant SN 1006 is a source of high-energy particles and its southwestern limb is interacting with a dense ambient cloud, thus being a promising region for gamma-ray hadronic emission. We aim at describing the physics and the nonthermal emission associated with the shock-cloud interaction to derive the physical parameters of the cloud (poorly constrained by the data analysis), to ascertain the origin of the observed spatial variations in the spectral properties of the X-ray synchrotron emission, and to predict spectral and morphological features of the resulting gamma-ray emission. We performed 3-D magnetohydrodynamic simulations modeling the evolution of SN 1006 and its interaction with the ambient cloud, and explored different model setups. By applying the REMLIGHT code on the model results, we synthesized the synchrotron X-ray emission, and compared it with actual observations, to constrain the parameters of the model. We also synthesized the leptonic and hadronic gamma-ray emission from the models, deriving constraints on the energy content of the hadrons accelerated at the southwestern limb. We found that the impact of the SN 1006 shock front with a uniform cloud with density 0.5 cm^-3 can explain the observed morphology, the azimuthal variations of the cutoff frequency of the X-ray synchrotron emission, and the shock proper motion in the interaction region. Our results show that the current upper limit for the total hadronic energy in the southwestern limb is 2.5e49 erg.

preprint2016arXiv

XMM-Newton Large Program on SN1006 - II: Thermal Emission

Based on the XMM-Newton large program on SN1006 and our newly developed spatially resolved spectroscopy tools (Paper~I), we study the thermal emission from ISM and ejecta of SN1006 by analyzing the spectra extracted from 583 tessellated regions dominated by thermal emission. With some key improvements in spectral analysis as compared to Paper~I, we obtain much better spectral fitting results with less residuals. The spatial distributions of the thermal and ionization states of the ISM and ejecta show different features, which are consistent with a scenario that the ISM (ejecta) is heated and ionized by the forward (reverse) shock propagating outward (inward). Different elements have different spatial distributions and origins, with Ne mostly from the ISM, Si and S from the ejecta, and O and Mg from both ISM and ejecta. Fe L-shell lines are only detected in a small shell-like region SE to the center of SN1006, indicating that most of the Fe-rich ejecta has not yet or just recently been reached by the reverse shock. The overall ejecta abundance patterns for most of the heavy elements, except for Fe and sometimes S, are consistent with typical Type~Ia SN products. The NW half of the SNR interior probably represents a region with turbulently mixed ISM and ejecta, so has enhanced emission from O, Mg, Si, S, lower ejecta temperature, and a large diversity of ionization age. In addition to the asymmetric ISM distribution, an asymmetric explosion of the progenitor star is also needed to explain the asymmetric ejecta distribution.

preprint2015arXiv

An X-ray characterization of the central region of the SNR G332.5-5.6

We present an X-ray analysis of the central region of supernova remnant (SNR) G332.5-5.6 through an exhaustive analysis of XMM-Netwon observations with complementary infrared observations. We characterize and discuss the origin of the observed X-ray morphology, which presents a peculiar plane edge over the west side of the central region. The morphology and spectral properties of the X-ray supernova remnant were studied using a single full frame XMM-Newton observation in the 0.3 to 10.0 keV energy band. Archival infrared WISE observations at 8, 12 and 24 μm were also used to investigate the properties of the source and its surroundings at different wavelengths. The results show that the extended X-ray emission is predominantly soft (0.3-1.2 keV) and peaks around 0.5 keV, which shows that it is an extremely soft SNR. X-ray emission correlates very well with central regions of bright radio emission. On the west side the radio/X-ray emission displays a plane-like feature with a terminal wall where strong infrared emission is detected. Our spatially resolved X-ray spectral analysis confirms that the emission is dominated by weak atomic emission lines of N, O, Ne, and Fe, all of them undetected in previous X-ray studies. These characteristics suggest that the X-ray emission is originated in an optically thin thermal plasma, whose radiation is well fitted by a non-equilibrium ionization collisional plasma (VNEI) X-ray emission model. Our study favors a scenario where G332.5-5.6 is expanding in a medium with an abrupt density change (the wall), likely a dense infrared emitting region of dust on the western side of the source.

preprint2015arXiv

X-ray emission from stellar jets by collision against high-density molecular clouds: an application to HH 248

We investigate the plausibility of detecting X-ray emission from a stellar jet that impacts against a dense molecular cloud. This scenario may be usual for classical T Tauri stars with jets in dense star-forming complexes. We first model the impact of a jet against a dense cloud by 2D axisymmetric hydrodynamic simulations, exploring different configurations of the ambient environment. Then, we compare our results with XMM-Newton observations of the Herbig-Haro object HH 248, where extended X-ray emission aligned with the optical knots is detected at the edge of the nearby IC 434 cloud. Our simulations show that a jet can produce plasma with temperatures up to 10 MK, consistent with production of X-ray emission, after impacting a dense cloud. We find that jets denser than the ambient medium but less dense than the cloud produce detectable X-ray emission only at the impact onto the cloud. From the exploration of the model parameter space, we constrain the physical conditions (jet density and velocity, cloud density) that reproduce well the intrinsic luminosity and emission measure of the X-ray source possibly associated with HH 248. Thus, we suggest that the extended X-ray source close to HH 248 corresponds to the jet impacting on a dense cloud.

preprint2015arXiv

XMM-Newton Large Program on SN1006 - I: Methods and Initial Results of Spatially-Resolved Spectroscopy

Based on our newly developed methods and the XMM-Newton large program of SN1006, we extract and analyze the spectra from 3596 tessellated regions of this SNR each with 0.3-8 keV counts $>10^4$. For the first time, we map out multiple physical parameters, such as the temperature ($kT$), electron density ($n_e$), ionization parameter ($n_et$), ionization age ($t_{ion}$), metal abundances, as well as the radio-to-X-ray slope ($α$) and cutoff frequency ($ν_{cutoff}$) of the synchrotron emission. We construct probability distribution functions of $kT$ and $n_et$, and model them with several Gaussians, in order to characterize the average thermal and ionization states of such an extended source. We construct equivalent width (EW) maps based on continuum interpolation with the spectral model of each regions. We then compare the EW maps of OVII, OVIII, OVII K$δ-ζ$, Ne, Mg, SiXIII, SiXIV, and S lines constructed with this method to those constructed with linear interpolation. We further extract spectra from larger regions to confirm the features revealed by parameter and EW maps, which are often not directly detectable on X-ray intensity images. For example, O abundance is consistent with solar across the SNR, except for a low-abundance hole in the center. This "O Hole" has enhanced OVII K$δ-ζ$ and Fe emissions, indicating recently reverse shocked ejecta, but also has the highest $n_et$, indicating forward shocked ISM. Therefore, a multi-temperature model is needed to decompose these components. The asymmetric metal distributions suggest there is either an asymmetric explosion of the SN or an asymmetric distribution of the ISM.

preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

preprint2012arXiv

The northwestern ejecta knot in SN 1006

Aims: We want to probe the physics of fast collision-less shocks in supernova remnants. In particular, we are interested in the non-equilibration of temperatures and particle acceleration. Specifically, we aim to measure the oxygen temperature with regards to the electron temperature. In addition, we search for synchrotron emission in the northwestern thermal rim. Methods: This study is part of a dedicated deep observational project of SN 1006 using XMM-Newton, which provides us with currently the best resolution spectra of the bright northwestern oxygen knot. We aim to use the reflection grating spectrometer to measure the thermal broadening of the O vii line triplet by convolving the emission profile of the remnant with the response matrix. Results: The line broadening was measured to be σ_e = 2.4 \pm 0.3 eV, corresponding to an oxygen temperature of 275$^{+72}_{-63}$ keV. From the EPIC spectra we obtain an electron temperature of 1.35 \pm 0.10 keV. The difference in temperature between the species provides further evidence of non-equilibration of temperatures in a shock. In addition, we find evidence for a bow shock that emits X-ray synchrotron radiation, which is at odds with the general idea that due to the magnetic field orientation only in the NE and SW region X-ray synchrotron radiation should be emitted. We find an unusual Hα and X-ray synchrotron geometry, in that the Hα emission peaks downstream of the synchrotron emission. This may be an indication for a peculiar Hα shock, in which the density is lower and neutral fraction are higher than in other supernova remnants, resulting in a peak in Hα emission further downstream of the shock.

preprint2011arXiv

Overionization in X-ray spectra: a new paradigm for Mixed-Morphology SNRs

Mixed-morphology SNRs are characterized by a shell-like radio emission, a centrally peaked X-ray morphology, and by interaction with molecular clouds. Many models have been proposed to explain these peculiar remnants, but their physical origin is still unclear. The recent discovery of over-ionized (i. e. recombining) ejecta in 3 mixed-morphology SNRs has dramatically challenged all the previous models and opened up new, unexpected scenarios. I review the main properties of these remnants and their peculiar X-ray spectral properties. I also discuss the hydrodynamic model developed to explain the presence of over-ionized ejecta in W49B and present a list of open issues that still need to be clarified.

preprint2011arXiv

Suzaku Detection of Diffuse Hard X-Ray Emission outside Vela X

Vela X is a large, 3x2 degrees, radio-emitting pulsar wind nebula (PWN) powered by the Vela pulsar in the Vela supernova remnant. Using four Suzaku/XIS observations pointed just outside Vela X, we find hard X-ray emission extending throughout the fields of view. The hard X-ray spectra are well represented by a power-law. The photon index is measured to be constant at Gamma~2.4, similar to that of the southern outer part of Vela X. The power-law flux decreases with increasing distance from the pulsar. These properties lead us to propose that the hard X-ray emission is associated with the Vela PWN. The larger X-ray extension found in this work strongly suggests that distinct populations relativistic electrons form the X-ray PWN and Vela X, as was recently inferred from multiwavelength spectral modeling of Vela X.

preprint2011arXiv

Unveiling the spatial structure of the overionized plasma in the supernova remnant W49B

W49B is a mixed-morphology supernova remnant with thermal X-ray emission dominated by the ejecta. In this remnant, the presence of overionized plasma has been directly established, with information about its spatial structure. However, the physical origin of the overionized plasma in W49B has not yet been understood. We investigate this intriguing issue through a 2D hydrodynamic model that takes into account, for the first time, the mixing of ejecta with the inhomogeneous circumstellar and interstellar medium, the thermal conduction, the radiative losses from optically thin plasma, and the deviations from equilibrium of ionization induced by plasma dynamics. The model was set up on the basis of the observational results. We found that the thermal conduction plays an important role in the evolution of W49B, in- ducing the evaporation of the circumstellar ring-like cloud (whose presence has been deduced from previous observations) that mingles with the surrounding hot medium, cooling down the shocked plasma, and pushes the ejecta backwards to the center of the remnant, forming there a jet-like structure. During the evolution, a large region of overionized plasma forms within the remnant. The overionized plasma originates from the rapid cooling of the hot plasma originally heated by the shock reflected from the dense ring-like cloud. In particular, we found two different ways for the rapid cooling of plasma to appear: i) the mixing of relatively cold and dense material evaporated from the ring with the hot shocked plasma and ii) the rapid adiabatic expansion of the ejecta. The spatial distribution of the radiative recombination continuum predicted by the numerical model is in good agreement with that observed.

preprint2011arXiv

Widespread nanoflare variability detected with Hinode/XRT in a solar active region

It is generally agreed that small impulsive energy bursts called nanoflares are responsible for at least some of the Sun's hot corona, but whether they are the explanation for most of the multi-million degree plasma has been a matter of ongoing debate. We here present evidence that nanoflares are widespread in an active region observed by the X-Ray Telescope on-board the Hinode mission. The distributions of intensity fluctuations have small but important asymmetries, whether taken from individual pixels, multi-pixel subregions, or the entire active region. Negative fluctuations (corresponding to reduced intensity) are greater in number but weaker in amplitude, so that the median fluctuation is negative compared to a mean of zero. Using MonteCarlo simulations, we show that only part of this asymmetry can be explained by Poisson photon statistics. The remainder is explainable with a tendency for exponentially decreasing intensity, such as would be expected from a cooling plasma produced from a nanoflare. We suggest that nanoflares are a universal heating process within active regions.

preprint2009arXiv

Generation of radiative knots in a randomly pulsed protostellar jet I. Dynamics and energetics

HH objects are characterized by a complex knotty morphology detected mainly along the axis of protostellar jets in a wide range of bands. Evidence of interactions between knots formed in different epochs have been found, suggesting that jets may result from the ejection of plasma blobs from the source. We aim at investigating the physical mechanism leading to the irregular knotty structure observed in jets in different bands and the complex interactions occurring among blobs of plasma ejected from the stellar source. We perform 2D axisymmetric HD simulations of a randomly ejected pulsed jet. The jet consists of a train of blobs which ram with supersonic speed into the ambient medium. The initial random velocity of each blob follows an exponential distribution. We explore the ejection rate parameter to derive constraints on the physical properties of protostellar jets by comparison of model results with observations. Our model takes into account radiative losses and thermal conduction. We find that the mutual interactions of blobs ejected at different epochs and with different speed lead to a variety of plasma components not described by current models. The main features characterizing the random pulsed jet scenario are: single high speed knots, showing a measurable proper motion in nice agreement with observations; irregular chains of knots aligned along the jet axis and possibly interacting with each other; reverse shocks interacting with outgoing knots; oblique shocks produced by the reflection of shocks at the jet cocoon. All these structures concur to determine the morphology of the jet in different bands. We also find that the thermal conduction plays a crucial role in damping out HD instabilities that would develop within the cocoon and that contribute to the jet breaking.