3D hydrodynamic numerical models for nebulae around runaway Wolf-Rayet stars

We present 3D hydrodynamical simulations of the circumstellar bubble from a Wolf-Rayet runaway star. In the models two properties were taken into account: (a) the proper motion of the central star through the interstellar medium (ISM) and (b) the evolution of the stellar wind from the red supergiant (RSG) stage to the Wolf-Rayet (WR) stage. From the hydrodynamic results synthetic X-ray maps in the [0.3 − 1.2] keV energy range were computed. These maps show that the bubble morphology is affected by the stellar motion, producing a bow shock in the RSG stage that can explain the limb-brightened morphology observed. Additionally, these synthetic maps show filamentary and clumpy appearance produced by instabilities triggered from the interaction between the WR wind and the RSG shell. It was found that these types of collisions can explain the origin of the X-ray emission observed in the nebulae of Wolf- Rayet stars.Fil: Reyes Iturbide, J.. Universidad Nacional Autónoma de México. Instituto de Ciencias Nucleares; México. Tecnológico de Estudios Superiores de Tianguistenco; MéxicoFil: Velázquez, Pablo F.. Universidad Nacional Autónoma de México. Instituto de Ciencias Nucleares; MéxicoFil: Rosado, M.. Universidad Nacional Autónoma de México. Instituto de Astronomía; MéxicoFil: Schneiter, Ernesto Matías. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; ArgentinaFil: Ramírez Ballinas, I.. Universidad Nacional Autónoma de México. Instituto de Astronomía; Méxic

X-ray emission and dynamics from large diameter superbubbles: The case of N 70 superbubble

The morphology, dynamics and thermal X-ray emission of the superbubble N70 is studied by means of 3D hydrodynamical simulations, carried out with the {\sc{yguaz\'u-a}} code. We have considered different scenarios: the superbubble being the product of a single supernova remnant, of the stellar winds from an OB association, or the result of the joint action of stellar winds and a supernova event. Our results show that, in spite that all scenarios produce bubbles with the observed physical size, only those where the bubble is driven by stellar winds and a SN event are successful to explain the general morphology, dynamics and the X-ray luminosity of N70. Our models predict temperatures in excess of $10^8 \mathrm{K}$ at the interior of the superbubble, however the density is too low and the emission in thermal X-ray above $2 \mathrm{keV}$ is too faint to be detected.Comment: 9 pages, 11 figures, ApJ accepte

3D numerical model of the Omega Nebula (M17): simulated thermal X-ray emission

We present 3D hydrodynamical simulations of the superbubble M17, also known as the Omega nebula, carried out with the adaptive grid code yguazu'-a, which includes radiative cooling. The superbubble is modelled considering the winds of 11 individual stars from the open cluster inside the nebula (NGC 6618), for which there are estimates of the mass loss rates and terminal velocities based on their spectral types. These stars are located inside a dense interstellar medium, and they are bounded by two dense molecular clouds. We carried out three numerical models of this scenario, considering different line of sight positions of the stars (the position in the plane of the sky is known, thus fixed). Synthetic thermal X-ray emission maps are calculated from the numerical models and compared with ROSAT observations of this astrophysical object. Our models reproduce successfully both the observed X-ray morphology and the total X-ray luminosity, without taking into account thermal conduction effects.Comment: 8 pages, 6 figures, accepted for publication in MNRA

An extensive study of dynamical friction in dwarf galaxies: the role of stars, dark matter, halo profiles and MOND

We investigate the in-spiraling timescales of globular clusters in dwarf spheroidal (dSph) and dwarf elliptical (dE) galaxies, due to dynamical friction. We address the problem of these timescales having been variously estimated in the literature as much shorter than a Hubble time. Using self-consistent two-component (dark matter and stars) models, we explore mechanisms which may yield extended dynamical friction timescales in such systems in order to explain why dwarf galaxies often show globular cluster systems. As a general rule, dark matter and stars both give a comparable contribution to the dynamical drag. By exploring various possibilities for their gravitational make-up, it is shown that these studies help constrain the parameters of the dark matter haloes in these galaxies, as well as to test alternatives to dark matter. Under the assumption of a dark haloes having a constant density core, dynamical friction timescales are naturally extended upwards of a Hubble time. Cuspy dark haloes yield timescales $\lesssim$ 4.5 Gyr, for any dark halo parameters in accordance with observations of stellar line-of-sight velocity dispersion in dwarf spheroidal galaxies. We find that under the hypothesis of MOND dynamics, due to the enhanced dynamical drag of the stars, the dynamical friction timescales would be extremely short. Taking the well-measured structural parameters of the Fornax dSph and its globular cluster system as a case study, we conclude that requiring dynamical friction timescales comparable to the Hubble time strongly favours dark haloes with a central core.Comment: 18 pages, four figures, final version, accepted in MNRA

Bias adjustment and ensemble recalibration methods for seasonal forecasting: a comprehensive intercomparison using the C3S dataset

This work presents a comprehensive intercomparison of diferent alternatives for the calibration of seasonal forecasts, ranging from simple bias adjustment (BA)-e.g. quantile mapping-to more sophisticated ensemble recalibration (RC) methods- e.g. non-homogeneous Gaussian regression, which build on the temporal correspondence between the climate model and the corresponding observations to generate reliable predictions. To be as critical as possible, we validate the raw model and the calibrated forecasts in terms of a number of metrics which take into account diferent aspects of forecast quality (association, accuracy, discrimination and reliability). We focus on one-month lead forecasts of precipitation and temperature from four state-of-the-art seasonal forecasting systems, three of them included in the Copernicus Climate Change Service dataset (ECMWF-SEAS5, UK Met Ofce-GloSea5 and Météo France-System5) for boreal winter and summer over two illustrative regions with diferent skill characteristics (Europe and Southeast Asia). Our results indicate that both BA and RC methods efectively correct the large raw model biases, which is of paramount importance for users, particularly when directly using the climate model outputs to run impact models, or when computing climate indices depending on absolute values/thresholds. However, except for particular regions and/or seasons (typically with high skill), there is only marginal added value-with respect to the raw model outputs-beyond this bias removal. For those cases, RC methods can outperform BA ones, mostly due to an improvement in reliability. Finally, we also show that whereas an increase in the number of members only modestly afects the results obtained from calibration, longer hindcast periods lead to improved forecast quality, particularly for RC methods.This work has been funded by the C3S activity on Evaluation and Quality Control for seasonal forecasts. JMG was partially supported by the project MULTI-SDM (CGL2015-66583-R, MINECO/FEDER). FJDR was partially funded by the H2020 EUCP project (GA 776613)

Discovery of a Be/X-ray pulsar binary and associated supernova remnant in the Wing of the SMC

We report on a new Be/X-ray pulsar binary located in the Wing of the Small Magellanic Cloud (SMC). The strong pulsed X-ray source was discovered with the Chandra and XMM-Newton X-ray observatories. The X-ray pulse period of 1062 s is consistently determined from both Chandra and XMM-Newton observations, revealing one of the slowest rotating X-ray pulsars known in the SMC. The optical counterpart of the X-ray source is the emission-line star 2dFS 3831. Its B0-0.5(III)e+ spectral type is determined from VLT-FLAMES and 2dF optical spectroscopy, establishing the system as a Be/X-ray binary (Be-XRB). The hard X-ray spectrum is well fitted by a power-law with additional thermal and blackbody components, the latter reminiscent of persistent Be-XRBs. This system is the first evidence of a recent supernova in the low density surroundings of NGC 602. We detect a shell nebula around 2dFS 3831 in H-alpha and [O III] images and conclude that it is most likely a supernova remnant. If it is linked to the supernova explosion that created this new X-ray pulsar, its kinematic age of (2-4)x10^4 yr provides a constraint on the age of the pulsar.Comment: 5 pages, 5 figures, accepted for publication in MNRAS Letter

Numerical heat conduction in hydrodynamical models of colliding hypersonic flows

Hydrodynamical models of colliding hypersonic flows are presented which explore the dependence of the resulting dynamics and the characteristics of the derived X-ray emission on numerical conduction and viscosity. For the purpose of our investigation we present models of colliding flow with plane-parallel and cylindrical divergence. Numerical conduction causes erroneous heating of gas across the contact discontinuity which has implications for the rate at which the gas cools. We find that the dynamics of the shocked gas and the resulting X-ray emission are strongly dependent on the contrast in the density and temperature either side of the contact discontinuity, these effects being strongest where the postshock gas of one flow behaves quasi-adiabatically while the postshock gas of the other flow is strongly radiative. Introducing additional numerical viscosity into the simulations has the effect of damping the growth of instabilities, which in some cases act to increase the volume of shocked gas and can re-heat gas via sub-shocks as it flows downstream. The resulting reduction in the surface area between adjacent flows, and therefore of the amount of numerical conduction, leads to a commensurate reduction in spurious X-ray emission, though the dynamics of the collision are compromised. The simulation resolution also affects the degree of numerical conduction. A finer resolution better resolves the interfaces of high density and temperature contrast and although numerical conduction still exists the volume of affected gas is considerably reduced. However, since it is not always practical to increase the resolution, it is imperative that the degree of numerical conduction is understood so that inaccurate interpretations can be avoided. This work has implications for the dynamics and emission from astrophysical phenomena which involve high Mach number shocks.Comment: 14 pages, 10 figures, accepted for publication in MNRA