Showing 1–50 of 60 results for author: Esquivel, A

Effects of Radiative Transfer on the Observed Anisotropy in MHD Turbulent Molecular Simulations

Authors: D. Hernández-Padilla, A. Esquivel, A. Lazarian, P. F. Velázquez, J. Cho

Abstract: We study the anisotropy of centroid and integrated intensity maps with synthetic observations. We perform post-process radiative transfer including the optically thick regime that was not covered in Hernández-Padilla et al. (2020). We consider the emission in various CO molecular lines, that range from optically thin to optically thick ($\mathrm{^{12}CO}$, $\mathrm{^{13}CO}$, $\mathrm{C^{18}O}$, a… ▽ More We study the anisotropy of centroid and integrated intensity maps with synthetic observations. We perform post-process radiative transfer including the optically thick regime that was not covered in Hernández-Padilla et al. (2020). We consider the emission in various CO molecular lines, that range from optically thin to optically thick ($\mathrm{^{12}CO}$, $\mathrm{^{13}CO}$, $\mathrm{C^{18}O}$, and $\mathrm{C^{17}O}$). The results for the velocity centroids are similar to those in the optically thin case. For instance, the anisotropy observed can be attributed to the Alfvén mode, which dominates over the slow and fast modes when the line of sight is at a high inclination with respect to the mean magnetic field. A few differences arise in the models with higher opacity, where some dependence on the sonic Mach number becomes evident. In contrast to the optically thin case, maps of integrated intensity become more anisotropic in optically thick lines. In this situation the scales probed are restricted, due to absorption, to smaller scales which are known to be more anisotropic. We discuss how the sonic Mach number can affect the latter results, with highly supersonic cases exhibiting a lower degree of anisotropy. △ Less

Submitted 21 June, 2024; v1 submitted 19 June, 2024; originally announced June 2024.

Detecting stellar flares in photometric data using hidden Markov models

Authors: J. Arturo Esquivel, Yunyi Shen, Vianey Leos-Barajas, Gwendolyn Eadie, Joshua Speagle, Radu V Craiu, Amber Medina, James Davenport

Abstract: We present a hidden Markov model (HMM) for discovering stellar flares in light curve data of stars. HMMs provide a framework to model time series data that are not stationary; they allow for systems to be in different states at different times and consider the probabilities that describe the switching dynamics between states. In the context of stellar flares discovery, we exploit the HMM framework… ▽ More We present a hidden Markov model (HMM) for discovering stellar flares in light curve data of stars. HMMs provide a framework to model time series data that are not stationary; they allow for systems to be in different states at different times and consider the probabilities that describe the switching dynamics between states. In the context of stellar flares discovery, we exploit the HMM framework by allowing the light curve of a star to be in one of three states at any given time step: Quiet, Firing, or Decaying. This three state HMM formulation is designed to enable straightforward identification of stellar flares, their duration, and associated uncertainty. This is crucial for estimating the flare's energy, and is useful for studies of stellar flare energy distributions. We combine our HMM with a celerite model that accounts for quasi periodic stellar oscillations. Through an injection recovery experiment, we demonstrate and evaluate the ability of our method to detect and characterize flares in stellar time series. We also show that the proposed HMM flags fainter and lower energy flares more easily than traditional sigma clipping methods. Lastly, we visually demonstrate that simultaneously conducting detrending and flare detection can mitigate biased estimations arising in multistage modelling approaches. Thus, this method paves a new way to calculating stellar flare energy. We conclude with an example application to one star observed by TESS, showing how the HMM compares with sigma clipping when using real data. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: To be submitted to AAS, comments welcomed

Shedding light on the ejection history of molecular outflows: Multiple velocity modes and precession

Authors: Veronica Lora, Thomas Nony, Alejandro Esquivel, Roberto Galván-Madrid

Abstract: Variable accretion has been well studied in evolved stages of low-mass stars formation. However, the accretion history in the initial phases of star formation is still a seldom studied topic. The outflows and jets emerging from protostellar objects could shed some light onto their accretion history. We consider the recently studied case of W43-MM1, a protocluster containing 46 outflows driven by 2… ▽ More Variable accretion has been well studied in evolved stages of low-mass stars formation. However, the accretion history in the initial phases of star formation is still a seldom studied topic. The outflows and jets emerging from protostellar objects could shed some light onto their accretion history. We consider the recently studied case of W43-MM1, a protocluster containing 46 outflows driven by 27 protostellar cores. The outflow kinematics of the individual cores and associated knots in W43-MM1 indicate episodic protostellar ejection. We took the observed parameters of an individual core system (core #8) and performed 3D hydrodynamic simulations of such system including episodic changes in the velocity of the outflow. The simulations consist of a collimated jet emerging from a core taking into account one- and two-velocity modes in the variation of the ejection velocity of the jet. In addition, we investigated the effect of including the precession of the jet to the one- and two-velocity mode models. From the simulations, we constructed position-velocity diagrams and compared them with the observations. We find that including a second mode in the ejection velocity, as well as the precession, are required to explain the positions of the outflow knots and other position-velocity features observed in the outflow of core #8 in W43-MM1. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 13 pages, 8 figures. Accepted for publication in ApJ

Synthetic X-ray emission from white dwarf accreting planetary material

Authors: S. Estrada-Dorado, V. Lora, J. A. Toalá, A. Esquivel, M. A. Guerrero, R. F. Maldonado, Y. -H. Chu

Abstract: The emission of hard X-rays associated with white dwarfs (WD) can be generated by the presence of a stellar companion either by the companion's coronal emission or by an accretion disk formed by material stripped from the companion. Recent studies have suggested that a Jupiter-like planet can also be donor of material whose accretion onto the WD can generate hard X-rays. We use the {\sc guacho} co… ▽ More The emission of hard X-rays associated with white dwarfs (WD) can be generated by the presence of a stellar companion either by the companion's coronal emission or by an accretion disk formed by material stripped from the companion. Recent studies have suggested that a Jupiter-like planet can also be donor of material whose accretion onto the WD can generate hard X-rays. We use the {\sc guacho} code to reproduce the conditions of this WD-planet scenario. With the example of the hard X-ray WD KPD\,0005+5106, we explore different terminal wind velocities and mass-loss rates of a donor planet for a future network of simulations to investigate the luminosity and the spectral and temporal properties of the hard X-ray emission in WD-planet systems. Our simulations show that the material stripped from the planet forms a disk and accretes onto the WD to reach temperatures high enough to generate hard X-rays as usually seen in X-ray binaries with low-mass companions. For high terminal wind velocities, the planet material does not form a disk, but it rather accretes directly onto the WD surface. The simulations reproduce the X-ray luminosity of another X-ray accreting WD (G\,29$-$38), and only for some times reaches the hard X-ray luminosity of KPD\,0005+5106. The X-ray variability is stochastic and does not reproduce the period of KPD\,0005+5106, suggesting that additional physical processes (e.g., hot spots resulting from magnetic channelling of the accreting material) need to be explored. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: 16 pages, 13 figures

Simulated non-thermal emission from SNR G1.9+0.3

Authors: Marco A. Villagran, D. O. Gómez, P. F. Velázquez, D. M. -A. Meyer, A. Chiotellis, A. C. Raga, A. Esquivel, J. C. Toledo-Roy, K. M. Vargas-Rojas, E. M. Schneiter

Abstract: Supernova remnants are the nebular leftover of defunct stellar environments, resulting from the interaction between a supernova blastwave and the circumstellar medium shaped by the progenitor throughout its life. They display a large variety of non-spherical morphologies such as ears that shine non-thermally. % We have modelled the structure and the non-thermal emission of the supernova remnant G1… ▽ More Supernova remnants are the nebular leftover of defunct stellar environments, resulting from the interaction between a supernova blastwave and the circumstellar medium shaped by the progenitor throughout its life. They display a large variety of non-spherical morphologies such as ears that shine non-thermally. % We have modelled the structure and the non-thermal emission of the supernova remnant G1.9+0.3 through 3D magnetohydrodynamic numerical simulations. We propose that the peculiar ear-shaped morphology of this supernova remnant results from the interaction of the its blast wave with a magnetized circumstellar medium, which was previously asymmetrically shaped by the past stellar wind emanating from the progenitor star or its stellar companion. We created synthetic non-thermal radio and x-ray maps from our simulated remnant structure, which are in qualitative agreement with observations, forming ears on the polar directions. Our synthetic map study explains the discrepancies between the measured non-thermal radio and X-ray surface brightness distributions assuming that the Inverse Compton process produces the observed X-ray emission. △ Less

Submitted 17 October, 2023; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS

Emerging planetary nebulae within 3D spiral patterns

Authors: V. Lora, J. A. Toalá, J. I. González-Carbajal, M. A. Guerrero, A. Esquivel, G. Ramos-Larios

Abstract: We present the first 3D radiation-hydrodynamic simulations of the formation of planetary nebulae (PNe) emerging from 3D spiral patterns. We use the GUACHO code to create 3D spiral structures as a consequence of the distortions on the geometry of the intrinsically isotropic wind of an asymptotic giant branch (AGB) star produced by a companion star in a circular orbit. We found that the orbital peri… ▽ More We present the first 3D radiation-hydrodynamic simulations of the formation of planetary nebulae (PNe) emerging from 3D spiral patterns. We use the GUACHO code to create 3D spiral structures as a consequence of the distortions on the geometry of the intrinsically isotropic wind of an asymptotic giant branch (AGB) star produced by a companion star in a circular orbit. We found that the orbital period of the binary producing the 3D spiral pattern has consequences on the formation and shaping of the PN itself. Stellar systems with longer period create less entwined 3D spirals, producing PNe with rounder inner cavities, and prevent the expansion of jet towards the polar directions. The spiral fitting procedure used in the literature to predict the binary's orbital period may be misleading in the case of proto-PNe and PNe as spiral patterns are diluted by their own thermal expansion down to the average AGB density profile within a few hundred years and are further disrupted by the action of jets. By adopting a phase of jet ejections between the AGB and post-AGB stages, we are able to recover the morphologies of proto-PNe and PNe that exhibit ring-like structures in their halos. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 12 pages, 12 figures, 1 table; Accepted to MNRAS

The sculpting of rectangular and jet-like morphologies in supernova remnants by anisotropic equatorially-confined progenitor stellar winds

Authors: P. F. Velazquez, D. M. -A. Meyer, A. Chiotellis, A. E. Cruz-Alvarez, E. M. Schneiter, J. C. Toledo-Roy, E. M. Reynoso, A. Esquivel

Abstract: Thermonuclear and core-collapse supernova remnants (SNRs) are the nebular leftovers of defunct stars. Their morphology and emission properties provide insights into the evolutionary history of the progenitor star. But while some SNRs are spherical, as expected from a point-like explosion expanding into a roughly uniform medium, many others exhibit complex non-spherical morphologies which are often… ▽ More Thermonuclear and core-collapse supernova remnants (SNRs) are the nebular leftovers of defunct stars. Their morphology and emission properties provide insights into the evolutionary history of the progenitor star. But while some SNRs are spherical, as expected from a point-like explosion expanding into a roughly uniform medium, many others exhibit complex non-spherical morphologies which are often not easily explained. In this work, we use three-dimensional magnetohydrodynamic simulations to show that rectangular and jet-like morphologies can be explained by supernovae (SNe), either type Ia or type II, expanding within anisotropic, bipolar stellar wind bubbles driven by the progenitor star. The stellar wind has an anisotropic density distribution, which channels the SN ejecta differently depending on the anisotropy characteristics. We compute synthetic thermal (X-ray) and non-thermal (synchrotron) emission maps from our numerical simulations to compare with observations. We find rectangular morphologies are generated when the stellar wind has a high mass loss rate and forms a dense, narrow disk at the equatorial region. Instead, a jet-like or ear-like morphology is obtained when the stellar wind develops a wide, dense disk. Stellar winds with low mass-loss rates do not strongly influence the SNR morphology. Finally, our synthetic synchrotron and X-ray maps for the high mass-loss rate case qualitatively agree with the observations of the SNRs G332.5-5.6 and G290.1-0.8. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: Accepted to MNRAS

Observational and numerical characterization of a recurrent arc-shaped front propagating along a coronal fan

Authors: M. V. Sieyra, S. Krishna Prasad, G. Stenborg, E. Khomenko, T. Van Doorsselaere, A. Costa, A. Esquivel, J. M. Riedl

Abstract: Recurrent, arc-shaped intensity disturbances were detected by EUV channels in an active region. The fronts were observed to propagate along a coronal loop bundle rooted in a small area within a sunspot umbra. Previous works have linked these intensity disturbances to slow magnetoacoustic waves that propagate from the lower atmosphere to the corona along the magnetic field. The slow magnetoacoustic… ▽ More Recurrent, arc-shaped intensity disturbances were detected by EUV channels in an active region. The fronts were observed to propagate along a coronal loop bundle rooted in a small area within a sunspot umbra. Previous works have linked these intensity disturbances to slow magnetoacoustic waves that propagate from the lower atmosphere to the corona along the magnetic field. The slow magnetoacoustic waves propagate at the local cusp speed. However, the measured propagation speeds from the intensity images are usually smaller as they are subject to projection effects due to the inclination of the magnetic field with respect to the line-of-sight. Here, we aim to understand the effect of projection by comparing observed speeds with those from a numerical model. Using multi-wavelength data we determine the periods present in the observations at different heights of the solar atmosphere through Fourier analysis. We calculate the plane-of-sky speeds along one of the loops from the cross-correlation time lags obtained as a function of distance along the loop. We perform a 2D ideal MHD simulation of an active region embedded in a stratified atmosphere. We drive slow waves from the photosphere with a 3 minutes periodicity. Synthetic time-distance maps are generated from the forward-modelled intensities in coronal wavelengths and the projected propagation speeds are calculated. The intensity disturbances show a dominant period between [2-3] minutes at different heights of the atmosphere. The apparent propagation speeds calculated for coronal channels exhibit an accelerated pattern with values increasing from 40 to 120 km/s as the distance along the loop rises. The propagation speeds obtained from the synthetic time-distance maps also exhibit accelerated profiles within a similar range of speeds. We conclude that the accelerated propagation in our observations is due to the projection effect. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 12 pages, 12 figures

Journal ref: A&A 667, A21 (2022)

Rectangular core-collapse supernova remnants: application to Puppis A

Authors: D. M. -A. Meyer, P. F. Velazquez, O. Petruk, A. Chiotellis, M. Pohl, F. Camps-Farina, M. Petrov, E. M. Reynoso, J. C. Toledo-Roy, E. M. Schneiter, A. Castellanos-Ramirez, A. Esquivel

Abstract: Core-collapse supernova remnants are the gaseous nebulae of galactic interstellar media (ISM) formed after the explosive death of massive stars. Their morphology and emission properties depend both on the surrounding circumstellar structure shaped by the stellar wind-ISM interaction of the progenitor star and on the local conditions of the ambient medium. In the warm phase of the Galactic plane (n… ▽ More Core-collapse supernova remnants are the gaseous nebulae of galactic interstellar media (ISM) formed after the explosive death of massive stars. Their morphology and emission properties depend both on the surrounding circumstellar structure shaped by the stellar wind-ISM interaction of the progenitor star and on the local conditions of the ambient medium. In the warm phase of the Galactic plane (n = 1/cm3, T = 8000 K), an organised magnetic field of strength 7 microG has profound consequences on the morphology of the wind bubble of massive stars at rest. In this paper we show through 2.5D magneto-hydrodynamical simulations, in the context of a Wolf-Rayet-evolving 35 Mo star, that it affects the development of its supernova remnant. When the supernova remnant reaches its middle age (15 to 20 kyr), it adopts a tubular shape that results from the interaction between the isotropic supernova ejecta and the anisotropic, magnetised, shocked stellar progenitor bubble into which the supernova blast wave expands. Our calculations for non-thermal emission, i.e. radio synchrotron and inverse Compton radiation, reveal that such supernova remnants can, due to projection effects, appear as rectangular objects in certain cases. This mechanism for shaping a supernova remnant is similar to the bipolar and elliptical planetary nebula production by wind-wind interaction in the low-mass regime of stellar evolution. If such a rectangular core-collapse supernova remnant is created, the progenitor star must not have been a runaway star. We propose that such a mechanism is at work in the shaping of the asymmetric core-collapse supernova remnant Puppis A. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: Accepted at MNRAS

Huffing, and puffing, and blowing your house in: Strong stellar winds interaction with a super massive black hole

Authors: V. Lora, A. Raga, J. Canto, A. Esquivel

Abstract: We present analytic and numerical models of a cluster wind flow resulting from the interaction of stellar winds of massive stars, with a super massive black hole (SMBH). We consider the motion of the stars as well as the gravitational force of the SMBH. In the numerical simulations we consider two cases: the first one with the stars is in circular orbits, and the second one with the stars in eccen… ▽ More We present analytic and numerical models of a cluster wind flow resulting from the interaction of stellar winds of massive stars, with a super massive black hole (SMBH). We consider the motion of the stars as well as the gravitational force of the SMBH. In the numerical simulations we consider two cases: the first one with the stars is in circular orbits, and the second one with the stars in eccentric orbits around the SMBH. We found that after the system reaches an equilibrium, the circular and elliptical cases are very similar. We found a very good agreement between the analytical and numerical results, not only from our numerical simulations but also from other high resolution numerical calculations. The analytical models are very interesting, since the properties of such complex systems involving strong winds and a massive compact object, can be rapidly inferred without the need of a numerical calculation. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Comments: 14 pages, 8 figures. Accepted for publication in ApJ

Estimating the propagation of a uniformly accelerated jet

Authors: J. I. Castorena, A. C. Raga, A. Esquivel, A. Rodríguez-González, L. Hernández-Martínez, J. Cantó, F. Clever

Abstract: We study the problem of a Herbig-Haro jet with a uniformly accelerating ejection velocity, travelling into a uniform environment. For the ejection density we consider two cases: a time-independent density, and a time-independent mass loss rate. For these two cases, we obtain analytic solutions for the motion of the jet head using a ram-pressure balance and a center of mass equation of motion. We a… ▽ More We study the problem of a Herbig-Haro jet with a uniformly accelerating ejection velocity, travelling into a uniform environment. For the ejection density we consider two cases: a time-independent density, and a time-independent mass loss rate. For these two cases, we obtain analytic solutions for the motion of the jet head using a ram-pressure balance and a center of mass equation of motion. We also compute axisymmetric numerical simulations of the same flow, and compare the time-dependent positions of the leading working surface shocks with the predictions of the two analytic models. We find that if the jet is over-dense and over-pressured (with respect to the environment) during its evolution, a good agreement is obtained with the analytic models, with the flow initially following the center of mass analytic solution, and (for the constant ejection density case) at later times approaching the ram-pressure balance solution. △ Less

Submitted 24 April, 2021; originally announced April 2021.

Comments: 15 pages, 5 figures

Formation and fate of the born-again planetary nebula HuBi 1

Authors: J. A. Toalá, V. Lora, B. Montoro-Molina, M. A. Guerrero, A. Esquivel

Abstract: We present the first 3D radiation-hydrodynamic simulations on the formation and evolution of born-again planetary nebulae (PNe), with particular emphasis to the case of HuBi1, the inside-out PN. We use the extensively-tested GUACHO code to simulate the formation of HuBi1 adopting mass-loss and stellar wind terminal velocity estimates obtained from observations presented by our group. We found that… ▽ More We present the first 3D radiation-hydrodynamic simulations on the formation and evolution of born-again planetary nebulae (PNe), with particular emphasis to the case of HuBi1, the inside-out PN. We use the extensively-tested GUACHO code to simulate the formation of HuBi1 adopting mass-loss and stellar wind terminal velocity estimates obtained from observations presented by our group. We found that, if the inner shell of HuBi1 was formed by an explosive very late thermal pulse (VLTP) ejecting material with velocities of $\sim$300 km s$^{-1}$, the age of this structure is consistent with that of $\simeq$200 yr derived from multi-epoch narrow-band imaging. Our simulations predict that, as a consequence of the dramatic reduction of the stellar wind velocity and photon ionizing flux during the VLTP, the velocity and pressure structure of the outer H-rich nebula are affected creating turbulent ionized structures surrounding the inner shell. These are indeed detected in Gran Telescopio Canarias MEGARA optical observations. Furthermore, we demonstrate that the current relatively low ionizing photon flux from the central star of HuBi1 is not able to completely ionize the inner shell, which favors previous suggestions that its excitation is dominated by shocks. Our simulations suggest that the kinetic energy of the H-poor ejecta of HuBi1 is at least 30 times that of the clumps and filaments in the evolved born-again PNe A30 and A78, making it a truly unique VLTP event. △ Less

Submitted 28 May, 2021; v1 submitted 21 March, 2021; originally announced March 2021.

Comments: 9 pages, 10 figures, no tables; accepted to MNRAS

H$_2$ mass-velocity relationship from 3D numerical simulations of jet-driven molecular outflows

Authors: Adriano Cerqueira, Bertrand Lefloch, Alejandro Esquivel, Pedro Rivera-Ortiz, Claudio Codella, Cecilia Ceccarelli, Linda Podio

Abstract: Previous numerical studies have shown that in protostellar outflows, the mass-velocity distribution $m(v)$ can be well described by a broken power law $\propto v^{- γ}$. On the other hand, recent observations of a sample of outflows show that the CO intensity-velocity distribution, closely related to $m(v)$, follows an exponential law $\propto \exp(-v/v_0)$. In the present work, we revisit the phy… ▽ More Previous numerical studies have shown that in protostellar outflows, the mass-velocity distribution $m(v)$ can be well described by a broken power law $\propto v^{- γ}$. On the other hand, recent observations of a sample of outflows show that the CO intensity-velocity distribution, closely related to $m(v)$, follows an exponential law $\propto \exp(-v/v_0)$. In the present work, we revisit the physical origin of the mass-velocity relationship $m(v)$ in jet-driven protostellar outflows. We investigate the respective contributions of the different regions of the outflow, from the swept-up ambient gas to the jet. We performed 3D numerical simulations of a protostellar jet propagating into a molecular cloud using the hydrodynamical code Yguazu-a. The code takes into account atomic and ionic species and was modified to include the H$_2$ gas. We find that by excluding the jet contribution, $m(v)$ is satisfyingly fitted with a single exponential law, with $v_0$ well in the range of observational values. The jet contribution results in additional components in the mass-velocity relationship. This empirical mass-velocity relationship is found to be valid locally in the outflow. The exponent $v_0$ is almost constant in time and for a given level of mixing between the ambient medium and the jet material. In general, $v_0$ displays only a weak spatial dependence. A simple modeling of the L1157 outflow successfully reproduces the various components of the observed CO intensity-velocity relationship. Our simulations indicate that these components trace the outflow cavity of swept-up gas and the material entrained along the jet, respectively. The CO intensity-velocity exponential law is naturally explained by the jet-driven outflow model. The entrained material plays an important role in shaping the mass-velocity profile. △ Less

Submitted 14 December, 2020; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: 15 pages, 11 figures. Accepted for publication in A&A. The abstract has been abridged due to the limitations on arXiv. The title has been corrected

Journal ref: A&A 645, A135 (2021)

GJ 436b and the stellar wind interaction: simulations constraints using Ly$α$ and H$α$ transits

Authors: C. Villarreal D'Angelo, A. A. Vidotto, A. Esquivel, G. Hazra, A. Youngblood

Abstract: The GJ 436 planetary system is an extraordinary system. The Neptune-size planet that orbits the M3 dwarf revealed in the Ly$α$ line an extended neutral hydrogen atmosphere. This material fills a comet-like tail that obscures the stellar disc for more than 10 hours after the planetary transit. Here, we carry out a series of 3D radiation hydrodynamic simulations to model the interaction of the stell… ▽ More The GJ 436 planetary system is an extraordinary system. The Neptune-size planet that orbits the M3 dwarf revealed in the Ly$α$ line an extended neutral hydrogen atmosphere. This material fills a comet-like tail that obscures the stellar disc for more than 10 hours after the planetary transit. Here, we carry out a series of 3D radiation hydrodynamic simulations to model the interaction of the stellar wind with the escaping planetary atmosphere. With these models, we seek to reproduce the $\sim56\%$ absorption found in Ly$α$ transits, simultaneously with the lack of absorption in H$α$ transit. Varying the stellar wind strength and the EUV stellar luminosity, we search for a set of parameters that best fit the observational data. Based on Ly$α$ observations, we found a stellar wind velocity at the position of the planet to be around [250-460] km s$^{-1}$ with a temperature of $[3-4]\times10^5$ K. The stellar and planetary mass loss rates are found to be $2\times 10^{-15}$ M$_\odot$ yr$^{-1}$ and $\sim[6-10]\times10^9$ g s$^{-1}$, respectively, for a stellar EUV luminosity of $[0.8-1.6]\times10^{27}$ erg s$^{-1}$. For the parameters explored in our simulations, none of our models present any significant absorption in the H$α$ line in agreement with the observations. △ Less

Submitted 9 December, 2020; originally announced December 2020.

Velocity Centroids Anisotropy and the Signature of different MHD Modes in the Turbulent ISM

Authors: D. Hernández-Padilla, A. Esquivel, A. Lazarian, D. Pogosyan, D. Kandel, J. Cho

Abstract: Magnetic turbulence is anisotropic as the directions of motion are constrained by the magnetic field. Such anisotropy can be observed in velocity centroids obtained from spectroscopic observations. We use magnetohydrodynamics (MHD) simulations to produce synthetic spectroscopic observations (position-position-velocity data) and study the anisotropy in the structure function of velocity centroid ma… ▽ More Magnetic turbulence is anisotropic as the directions of motion are constrained by the magnetic field. Such anisotropy can be observed in velocity centroids obtained from spectroscopic observations. We use magnetohydrodynamics (MHD) simulations to produce synthetic spectroscopic observations (position-position-velocity data) and study the anisotropy in the structure function of velocity centroid maps. We decomposed the velocity in the simulations into Alfvén, slow and fast-modes and studied how each of them contribute to the observed anisotropy. We found that when the angle between the line of sight and the mean magnetic field is large the Alfvén-mode dominates the observed anisotropy, while for smaller angles the anisotropy is not large enough to be used to probe the magnetization of the media, and it is dominated by the slow-mode. Our results are in fair agreement with the theoretical predictions in Kandel et al.(2016,2017). △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 17 pages, 7 figures, accepted for publication in The Astrophysical Journal

Solar slow magneto-acoustic-gravity waves: an erratum correction and a revisited scenario

Authors: E. Zurbriggen, M. V. Sieyra, A. Costa, A. Esquivel, G. Stenborg

Abstract: Slow waves are commonly observed on the entire solar atmosphere. Assuming a thin flux tube approximation, the cut-off periods of slow-mode magneto-acoustic-gravity waves that travel from the photosphere to the corona were obtained in Costa et al. (2018). In that paper, however, a typo in the specific heat coefficient at constant pressure $c_{\mathrm{p}}$ value led to an inconsistency in the cut-of… ▽ More Slow waves are commonly observed on the entire solar atmosphere. Assuming a thin flux tube approximation, the cut-off periods of slow-mode magneto-acoustic-gravity waves that travel from the photosphere to the corona were obtained in Costa et al. (2018). In that paper, however, a typo in the specific heat coefficient at constant pressure $c_{\mathrm{p}}$ value led to an inconsistency in the cut-off calculation, which is only significant at the transition region. Due to the abrupt temperature change in the region, a change of the mean atomic weight (by a factor of approximately two) also occurs, but is often overlooked in analytical models for simplicity purposes. In this paper, we revisit the calculation of the cut-off periods of magneto-acoustic-gravity waves in Costa et al. (2018) by considering an atmosphere in hydrostatic equilibrium with a temperature profile, with the inclusion of the variation of the mean atomic weight and the correction of the inconsistency aforementioned. In addition, we show that the cut-off periods obtained analytically are consistent with the corresponding periods measured in observations of a particular active region. △ Less

Submitted 20 April, 2020; originally announced April 2020.

Comments: 12 pages, 7 figures. MNRAS

Star-planet interaction through spectral lines

Authors: C. Villarreal D'Angelo, A. A. Vidotto, A. Esquivel, M. A. Sgró, T. Koskinen, L. Fossati

Abstract: The growth of spectroscopic observations of exoplanetary systems allows the possibility of testing theoretical models and studying the interaction that exoplanetary atmospheres have with the wind and the energetic photons from the star. In this work, we present a set of numerical 3D simulations of HD 209458b for which spectral lines observations of their evaporative atmosphere are available. The d… ▽ More The growth of spectroscopic observations of exoplanetary systems allows the possibility of testing theoretical models and studying the interaction that exoplanetary atmospheres have with the wind and the energetic photons from the star. In this work, we present a set of numerical 3D simulations of HD 209458b for which spectral lines observations of their evaporative atmosphere are available. The different simulations aim to reproduce different scenarios for the star-planet interaction. With our models, we reconstruct the Ly$α$ line during transit and compare with observations. The results allow us to analyse the shape of the line profile under these different scenarios and the comparison with the observations suggests that HD209458b may have a magnetic field off less than 1 G. We also explore the behaviour of the magnesium lines for models with and without magnetic fields. △ Less

Submitted 8 December, 2019; originally announced December 2019.

Hydrodynamical interaction of stellar and planetary winds: effects of charge exchange and radiation pressure on the observed Ly$α$ absorption

Authors: A. Esquivel, M. Schneiter, C. Villarreal D'Angelo, M. A. Sgró, L. Krapp

Abstract: Lyman $α$ observations of the transiting exoplanet HD 209458b enable the study of exoplanets exospheres exposed to stellar EUV fluxes, as well as the interacting stellar wind properties. In this study we present 3D hydrodynamical models for the stellar-planetary wind interaction including radiation pressure and charge exchange, together with photoionization, recombination and collisional ionizatio… ▽ More Lyman $α$ observations of the transiting exoplanet HD 209458b enable the study of exoplanets exospheres exposed to stellar EUV fluxes, as well as the interacting stellar wind properties. In this study we present 3D hydrodynamical models for the stellar-planetary wind interaction including radiation pressure and charge exchange, together with photoionization, recombination and collisional ionization processes. Our models explore the contribution of the radiation pressure and charge exchange on the Ly$α$ absorption profile in a hydrodynamical framework, and for a single set of stellar wind parameters appropriate for HD 209458. We find that most of the absorption is produced by the material from the planet, with a secondary contribution of neutralized stellar ions by charge exchange. At the same time, the hydrodynamic shock heats up the planetary material, resulting in a broad thermal profile. Meanwhile, the radiation pressure yielded a small velocity shift of the absorbing material. While neither charge exchange nor radiation pressure provide enough neutrals at the velocity needed to explain the observations at $-100~\mathrm{km~s^{-1}}$ individually, we find that the two effects combined with the broad thermal profile are able to explain the observations. △ Less

Submitted 10 July, 2019; originally announced July 2019.

Comments: 12 pages, 6 Figures. Accepted for publication in MNRAS

Magnetised winds and their influence in the escaping upper atmosphere of HD 209458b

Authors: Carolina Villarreal D'Angelo, Alejandro Esquivel, Matías Schneiter, Mario Agustín Sgró

Abstract: Lyman $α$ observations during an exoplanet transit have proved to be very useful to study the interaction between the stellar wind and the planetary atmosphere. They have been extensively used to constrain planetary system parameters that are not directly observed, such as the planetary mass loss rate. In this way, Ly $α$ observations can be a powerful tool to infer the existence of a planetary ma… ▽ More Lyman $α$ observations during an exoplanet transit have proved to be very useful to study the interaction between the stellar wind and the planetary atmosphere. They have been extensively used to constrain planetary system parameters that are not directly observed, such as the planetary mass loss rate. In this way, Ly $α$ observations can be a powerful tool to infer the existence of a planetary magnetic field, since it is expected that the latter will affect the escaping planetary material. To explore the effect that magnetic fields have on the Ly $α$ absorption of HD 209458b, we run a set of 3D MHD simulations including dipolar magnetic fields for the planet and the star. We assume values for the surface magnetic field at the poles of the planet in the range of [0-5] G, and from 1 to 5 G at the poles of the star. Our models also include collisional and photo-ionisation, radiative recombination, and an approximation for the radiation pressure. Our results show that the magnetic field of the planet and the star change the shape of the Ly $α$ absorption profile, since it controls the extent of the planetary magnetosphere and the amount of neutral material inside it. The model that best reproduces the absorption observed in HD 209458b (with canonical values for the stellar wind parameters) corresponds to a dipole planetary field of <~ 1 gauss at the poles. △ Less

Submitted 6 June, 2018; originally announced June 2018.

Origin of the bilateral structure of the supernova remnant G296.5+10

Authors: A. Moranchel-Basurto, P. F. Velázquez, E. Giacani, J. C. Toledo-Roy, E. M. Schneiter, F. De Colle, A. Esquivel

Abstract: In the present work we have modeled the supernova remnant (SNR) G296.5+10, by means of 3D magnetohydrodynamics (MHD) simulations. This remnant belongs to the bilateral SNR group and has an additional striking feature: the rotation measure (RM) in its eastern and western parts are very different. In order to explain both the morphology observed in radio-continuum and the RM, we consider that the re… ▽ More In the present work we have modeled the supernova remnant (SNR) G296.5+10, by means of 3D magnetohydrodynamics (MHD) simulations. This remnant belongs to the bilateral SNR group and has an additional striking feature: the rotation measure (RM) in its eastern and western parts are very different. In order to explain both the morphology observed in radio-continuum and the RM, we consider that the remnant expands into a medium shaped by the superposition of the magnetic field of the progenitor star with a constant Galactic magnetic field. We have also carried out a polarization study from our MHD results, obtaining synthetic maps of the linearly polarized intensity and the Stokes parameters. This study reveals that both the radio morphology and the reported RM for G$296.5+10$ can be explained if the quasi-parallel acceleration mechanism is taking place in the shock front of this remnant. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: 9 pages, 7 figures, accepted for publication in MNRAS

Proper motions of the HH1 jet

Authors: A. C. Raga, B. Reipurth, A. Esquivel, A. Castellanos-Ramirez, P. F. Velazquez, L. Hernandez-Martinez, A. Rodriguez-Gonzalez, J. E. Rechy-Garcia, D. Estrella-Trujillo, J. Bally, D. Gonzalez-Gomez, A. Riera

Abstract: We describe a new method for determining proper motions of extended objects, and a pipeline developed for the application of this method. We then apply this method to an analysis of four epochs of [S~II] HST images of the HH~1 jet (covering a period of $\sim 20$~yr). We determine the proper motions of the knots along the jet, and make a reconstruction of the past ejection velocity time-variabili… ▽ More We describe a new method for determining proper motions of extended objects, and a pipeline developed for the application of this method. We then apply this method to an analysis of four epochs of [S~II] HST images of the HH~1 jet (covering a period of $\sim 20$~yr). We determine the proper motions of the knots along the jet, and make a reconstruction of the past ejection velocity time-variability (assuming ballistic knot motions). This reconstruction shows an "acceleration" of the ejection velocities of the jet knots, with higher velocities at more recent times. This acceleration will result in an eventual merging of the knots in $\sim 450$~yr and at a distance of $\sim 80"$ from the outflow source, close to the present-day position of HH~1. △ Less

Submitted 4 August, 2017; originally announced August 2017.

Comments: 12 pages, 8 figures

MHD simulations of coronal supra-arcade downflows including anisotropic thermal conduction

Authors: E. Zurbriggen, A. Costa, A. Esquivel, M. Schneiter, M. Cécere

Abstract: Coronal supra-arcade downflows (SADs) are observed as dark trails descending towards hot turbulent fan shaped regions. Due to the large temperature values, and gradients in these fan regions the thermal conduction should be very efficient. While several models have been proposed to explain the triggering and the evolution of SADs, none of these scenarios address a systematic consideration of therm… ▽ More Coronal supra-arcade downflows (SADs) are observed as dark trails descending towards hot turbulent fan shaped regions. Due to the large temperature values, and gradients in these fan regions the thermal conduction should be very efficient. While several models have been proposed to explain the triggering and the evolution of SADs, none of these scenarios address a systematic consideration of thermal conduction. Thus, we accomplish this task numerically simulating the evolution of SADs within this framework. That is, SADs are conceived as voided (subdense) cavities formed by non-linear waves triggered by downflowing bursty localized reconnection events in a perturbed hot fan. We generate a properly turbulent fan, obtained by a stirring force that permits control of the energy and vorticity input in the medium where SADs develop. We include anisotropic thermal conduction and consider plasma properties consistent with observations. Our aim is to study if it is possible to prevent SADs to vanish by thermal diffusion. We find that this will be the case, depending on the turbulence parameters. In particular, if the magnetic field lines are able to envelope the voided cavities, thermally isolating them from the hot environment. Velocity shear perturbations that are able to generate instabilities of the Kelvin-Helmholtz type help to produce magnetic islands, extending the life-time of SADs. △ Less

Submitted 16 March, 2017; originally announced March 2017.

Journal ref: ApJ, 832:74, 2016

A 3D MHD simulation of SN 1006: a polarized emission study for the turbulent case

Authors: P. F. Velázquez, E. M. Schneiter, E. M. Reynoso, A. Esquivel, F. De Colle, J. C. Toledo-Roy, D. O. Gómez, M. V. Sieyra, A. Moranchel-Basurto

Abstract: Three dimensional magnetohydrodynamical simulations were carried out in order to perform a new polarization study of the radio emission of the supernova remnant SN 1006. These simulations consider that the remnant expands into a turbulent interstellar medium (including both magnetic field and density perturbations). Based on the referenced-polar angle technique, a statistical study was done on obs… ▽ More Three dimensional magnetohydrodynamical simulations were carried out in order to perform a new polarization study of the radio emission of the supernova remnant SN 1006. These simulations consider that the remnant expands into a turbulent interstellar medium (including both magnetic field and density perturbations). Based on the referenced-polar angle technique, a statistical study was done on observational and numerical magnetic field position-angle distributions. Our results show that a turbulent medium with an adiabatic index of 1.3 can reproduce the polarization properties of the SN 1006 remnant. This statistical study reveals itself as a useful tool for obtaining the orientation of the ambient magnetic field, previous to be swept up by the main supernova remnant shock. △ Less

Submitted 19 January, 2017; originally announced January 2017.

Comments: 8 pages, 7 figures, accepted for publication in MNRAS

The metallic winds in dwarf galaxies

Authors: F. Robles-Valdez, A. Rodríguez-González, L. Hernández-Martínez, A. Esquivel

Abstract: We present results from models of galactic winds driven by energy injected from nuclear (at the galactic center) and non-nuclear starbursts. The total energy of the starburst is provided by very massive young stellar clusters,which can push the galactic interstellar medium and produce an important outflow. Such outflow can be a well, or partially mixed wind, or a highly metallic wind. We have perf… ▽ More We present results from models of galactic winds driven by energy injected from nuclear (at the galactic center) and non-nuclear starbursts. The total energy of the starburst is provided by very massive young stellar clusters,which can push the galactic interstellar medium and produce an important outflow. Such outflow can be a well, or partially mixed wind, or a highly metallic wind. We have performed adiabatic 3D N-Body/Smooth Particle Hydrodynamics simulations of galactic winds using the GADGET-2 code. The numerical models cover a wide range of parameters, varying the galaxy concentration index, gas fraction of the galactic disk, and radial distance of the starburst. We show that an off-center starburst in dwarf galaxies is the most effective mechanism to produce a significant loss of metals (material from the starburst itself). At the same time a non-nuclear starburst produce a high efficiency of metal loss, in spite of having a moderate to low mass loss rate. △ Less

Submitted 8 December, 2016; originally announced December 2016.

Comments: 13 pages, 7 figures, accepted for publication in ApJ

The time evolution of HH~2 from four epochs of HST images

Authors: A. C. Raga, B. Reipurth, P. F. Velázquez, A. Esquivel, J. Bally

Abstract: We have analyzed four epochs of H$α$ and [S~II] HST images of the HH~1/2 outflow (covering a time interval from 1994 to 2014) to determine proper motions and emission line fluxes of the knots of HH~2. We find that our new proper motions agree surprisingly well with the motions measured by Herbig \& Jones (1981), although there is partial evidence for a slight deceleration of the motion of the HH~2… ▽ More We have analyzed four epochs of H$α$ and [S~II] HST images of the HH~1/2 outflow (covering a time interval from 1994 to 2014) to determine proper motions and emission line fluxes of the knots of HH~2. We find that our new proper motions agree surprisingly well with the motions measured by Herbig \& Jones (1981), although there is partial evidence for a slight deceleration of the motion of the HH~2 knots from 1945 to 2014. We also measure the time-variability of the H$α$ intensities and the [S~II]/H$α$ line ratios, and find that knots H and A have the largest intensity variabilities (in $1994\to 2014$). Knot H (which now dominates the HH~2 emission) has strengthened substantially, while keeping an approximately constant [S~II]/H$α$ ratio. Knot A has dramatically faded, and at the same time has had a substantial increase in its [S~II]/H$α$ ratio. Possible interpretations of these results are discussed. △ Less

Submitted 6 October, 2016; originally announced October 2016.

Comments: 17 pages, 10 figures. Accepted for publication in The Astronomical Journal

The time evolution of HH 1 from four epochs of HST images

Authors: A. C. Raga, B. Reipurth, A. Esquivel, J. Bally

Abstract: We present an analysis of four epochs of H$α$ and [S II] $λλ$ 6716/6731 HST images of HH 1. For determining proper motions we explore a new method based on analysis of spatially degraded images obtained convolving the images with wavelet functions of chosen widths. With this procedure we are able to generate maps of proper motion velocities along and across the outflow axis, as well as (angularly… ▽ More We present an analysis of four epochs of H$α$ and [S II] $λλ$ 6716/6731 HST images of HH 1. For determining proper motions we explore a new method based on analysis of spatially degraded images obtained convolving the images with wavelet functions of chosen widths. With this procedure we are able to generate maps of proper motion velocities along and across the outflow axis, as well as (angularly integrated) proper motion velocity distributions. From the four available epochs, we find the time evolution of the velocities, intensities and spatial distribution of the line emission. We find that over the last two decades HH 1 shows a clear acceleration. Also, the H$α$ and [S II] intensities have first dropped, and then recovered in the more recent (2014) images. Finally, we show a comparison between the two available HST epochs of [O III] $λ$ 5007 (1994 and 2014), in which we see a clear drop in the value of the [O III]/H$α$ ratio. △ Less

Submitted 4 March, 2016; originally announced March 2016.

Comments: 28 pages, 11 figures, accepted for publication in The Astronomical Journal

Photo-ionization of planetary winds: case study HD209458b

Authors: E. M. Schneiter, A. Esquivel, C. S. Villarreal D'Angelo, P. F. Velazquez, A. C. Raga, A. Costa

Abstract: Close-in hot Jupiters are exposed to a tremendous photon flux that ionizes the neutral escaping material from the planet leaving an observable imprint that makes them an interesting laboratory for testing theoretical models. In this work we present 3D hydrodynamic simulations with radiation transfer calculations of a close-in exoplanet in a blow-off state. We calculate the Ly-$α$ absorption and co… ▽ More Close-in hot Jupiters are exposed to a tremendous photon flux that ionizes the neutral escaping material from the planet leaving an observable imprint that makes them an interesting laboratory for testing theoretical models. In this work we present 3D hydrodynamic simulations with radiation transfer calculations of a close-in exoplanet in a blow-off state. We calculate the Ly-$α$ absorption and compare it with observations of HD 209458b an previous simplified model results.Our results show that the hydrodynamic interaction together with a proper calculation of the photoionization proccess are able to reproduce the main features of the observed Ly-$α$ absorption, in particular at the blue-shifted wings of the line. We found that the ionizing stellar flux produce an almost linear effect on the amount of absorption in the wake. Varying the planetary mass loss rate and the radiation flux, we were able to reproduce the $10\%$ absorption observed at $-100~\mathrm{km~s^{-1}}$. △ Less

Submitted 8 January, 2016; originally announced January 2016.

Comments: 9 pages, 6 figures

A Study of Radio Polarization in Protostellar Jets

Authors: Mariana Cécere, Pablo F. Velázquez, Anabella T. Araudo, Fabio De Colle, Alejandro Esquivel, Carlos Carrasco-González, Luis F. Rodríguez

Abstract: Synchrotron radiation is commonly observed in connection with shocks of different velocities, ranging from relativistic shocks associated with active galactic nuclei, gamma-ray bursts or microquasars to weakly- or non-relativistic flows as those observed in supernova remnants. Recent observations of synchrotron emission in protostellar jets are important not only because they extend the range over… ▽ More Synchrotron radiation is commonly observed in connection with shocks of different velocities, ranging from relativistic shocks associated with active galactic nuclei, gamma-ray bursts or microquasars to weakly- or non-relativistic flows as those observed in supernova remnants. Recent observations of synchrotron emission in protostellar jets are important not only because they extend the range over which the acceleration process works, but also because they allow to determine the jet and/or interstellar magnetic field structure, thus giving insights on the jet ejection and collimation mechanisms. In this paper, we compute for the first time polarized (synchrotron) and non polarized (thermal X-ray) synthetic emission maps from axisymmetrical simulations of magnetized protostellar jets. We consider models with different jet velocities and variability, as well as toroidal or helical magnetic field. Our simulations show that variable, low-density jets with velocities $\sim$ 1000 km s$^{-1}$ and $\sim$ 10 times lighter than the environment can produce internal knots with significant synchrotron emission, and thermal X-rays in the shocked region of the leading bow shock moving in a dense medium. While models with a purely toroidal magnetic field show a very large degree of polarization, models with helical magnetic field show lower values and a decrease of the degree of polarization, in agreement with observations of protostellar jets. △ Less

Submitted 22 November, 2015; v1 submitted 11 June, 2015; originally announced June 2015.

Comments: 10 pages, 8 figures. Accepted for publication in The Astrophysical Journal

The soft and hard X-rays thermal emission from star cluster winds with a supernova explosion

Authors: A. Castellanos-Ramirez, A. Rodriguez-Gonzalez, A. Esquivel, J. C. Toledo-Roy, J. Olivares, P. F. Velazquez

Abstract: Massive young star clusters contain dozens or hundreds of massive stars that inject mechanical energy in the form of winds and supernova explosions, producing an outflow which expands into their surrounding medium, shocking it and forming structures called superbubbles. The regions of shocked material can have temperatures in excess of 10$^6$ K, and emit mainly in thermal X-rays (soft and hard). T… ▽ More Massive young star clusters contain dozens or hundreds of massive stars that inject mechanical energy in the form of winds and supernova explosions, producing an outflow which expands into their surrounding medium, shocking it and forming structures called superbubbles. The regions of shocked material can have temperatures in excess of 10$^6$ K, and emit mainly in thermal X-rays (soft and hard). This X-ray emission is strongly affected by the action of thermal conduction, as well as by the metallicity of the material injected by the massive stars. We present three-dimensional numerical simulations exploring these two effects, metallicity of the stellar winds and supernova explosions, as well as thermal conduction. △ Less

Submitted 10 April, 2015; originally announced April 2015.

Comments: 14 pages, 12 figures. MNRAS. Accepted

Report number: MN-14-1802-MJ.R2

3D MHD simulation of polarized emission in SN 1006

Authors: E. M. Schneiter, P. F. Velázquez, E. M. Reynoso, A. Esquivel, F. De Colle

Abstract: We use three dimensional magnetohydrodynamic (MHD) simulations to model the supernova remnant SN 1006. From our numerical results, we have carried out a polarization study, obtaining synthetic maps of the polarized intensity, the Stokes parameter $Q$, and the polar-referenced angle, which can be compared with observational results. Synthetic maps were computed considering two possible particle acc… ▽ More We use three dimensional magnetohydrodynamic (MHD) simulations to model the supernova remnant SN 1006. From our numerical results, we have carried out a polarization study, obtaining synthetic maps of the polarized intensity, the Stokes parameter $Q$, and the polar-referenced angle, which can be compared with observational results. Synthetic maps were computed considering two possible particle acceleration mechanisms: quasi-parallel and quasi-perpendicular. The comparison of synthetic maps of the Stokes parameter $Q$ maps with observations proves to be a valuable tool to discern unambiguously which mechanism is taking place in the remnant of SN 1006, giving strong support to the quasi-parallel model. △ Less

Submitted 10 February, 2015; originally announced February 2015.

Comments: 6 pages, 4 figures, accepted by MNRAS

A two-mode planetary nebula luminosity function

Authors: A. Rodríguez-González, L. Hernández-Martínez, A. Esquivel, A. C. Raga, G. Stasińska, M. Peña, D. Mayya

Abstract: We propose a new Planetary Nebula Luminosity Function (PNLF) that includes two populations in the distribution. Our PNLF is a direct extension of the canonical function proposed by Jacoby et al. (1987), in order to avoid problems related with the histogram construction, it is cast in terms of cumulative functions. We are interested in recovering the shape of the faint part of the PNLF in a consist… ▽ More We propose a new Planetary Nebula Luminosity Function (PNLF) that includes two populations in the distribution. Our PNLF is a direct extension of the canonical function proposed by Jacoby et al. (1987), in order to avoid problems related with the histogram construction, it is cast in terms of cumulative functions. We are interested in recovering the shape of the faint part of the PNLF in a consistent manner, for galaxies with and without a dip in their PN luminosity functions. The parameters for the two mode PNLF are obtained with a genetic algorithm, which obtains a best fit to the PNLF varying all of the parameters simultaneously in a broad parameter space. We explore a sample of 9 galaxies with various Hubble types and construct their PNLF. All of the irregular galaxies, except one, are found to be consistent with a two-mode population, while the situation is less clear for ellipticals and spirals.For the case of NGC\, 6822, we show that the two-mode PNLF is consistent with previous studies of the star formation history within that galaxy. Our results support two episodes of star formation, in which the latter is significantly stronger. △ Less

Submitted 4 November, 2014; originally announced November 2014.

Comments: 10 pages, 12 figures (accepted). appears in Astronomy & Astrophysics, 2014

Journal ref: A&A 575, A1 (2015)

Measuring the Alfvenic Nature of the Interstellar Medium: Velocity Anisotropy Revisited

Authors: Blakesley Burkhart, A. Lazarian, I. C. Leao, J. R. de Medeiros, A. Esquivel

Abstract: The dynamics of the interstellar medium (ISM) are strongly affected by turbulence, which shows increased anisotropy in the presence of a magnetic field. We expand upon the Esquivel & Lazarian method to estimate the Alfven Mach number using the structure function anisotropy in velocity centroid data from position-position-velocity maps. We utilize 3D magnetohydrodynamic (MHD) simulations of fully d… ▽ More The dynamics of the interstellar medium (ISM) are strongly affected by turbulence, which shows increased anisotropy in the presence of a magnetic field. We expand upon the Esquivel & Lazarian method to estimate the Alfven Mach number using the structure function anisotropy in velocity centroid data from position-position-velocity maps. We utilize 3D magnetohydrodynamic (MHD) simulations of fully developed turbulence, with a large range of sonic and Alfvenic Mach numbers, to produce synthetic observations of velocity centroids with observational characteristics such as thermal broadening, cloud boundaries, noise, and radiative transfer effects of carbon monoxide. In addition, we investigate how the resulting anisotropy-Alfven Mach number dependency found in Esquivel & Lazarian (2011) might change when taking the second moment of the position-position-velocity cube or when using different expressions to calculate the velocity centroids. We find that the degree of anisotropy is related primarily to the magnetic field strength (i.e. Alfven Mach number) and the line-of-sight orientation, with a secondary effect on sonic Mach number. If the line-of-sight is parallel to up to ~45 deg off of the mean field direction, the velocity centroid anisotropy is not prominent enough to distinguish different Alfvenic regimes. The observed anisotropy is not strongly affected by including radiative transfer, although future studies should include additional tests for opacity effects. These results open up the possibility of studying the magnetic nature of the ISM using statistical methods in addition to existing observational techniques. △ Less

Submitted 20 August, 2014; originally announced August 2014.

Journal ref: 2014, ApJ, 790, 130

Modelling the 3D morphology and proper motions of the planetary nebula NGC 6302

Authors: L. Uscanga, P. F. Velázquez, A. Esquivel, A. C. Raga, P. Boumis, J. Cantó

Abstract: We present 3D hydrodynamical simulations of an isotropic fast wind interacting with a previously ejected toroidally-shaped slow wind in order to model both the observed morphology and the kinematics of the planetary nebula (PN) NGC 6302. This source, also known as the Butterfly nebula, presents one of the most complex morphologies ever observed in PNe. From our numerical simulations, we have obtai… ▽ More We present 3D hydrodynamical simulations of an isotropic fast wind interacting with a previously ejected toroidally-shaped slow wind in order to model both the observed morphology and the kinematics of the planetary nebula (PN) NGC 6302. This source, also known as the Butterfly nebula, presents one of the most complex morphologies ever observed in PNe. From our numerical simulations, we have obtained an intensity map for the H$α$ emission to make a comparison with the Hubble Space Telescope (HST) observations of this object. We have also carried out a proper motion (PM) study from our numerical results, in order to compare with previous observational studies. We have found that the two interacting stellar wind model reproduces well the morphology of NGC 6302, and while the PM in the models are similar to the observations, our results suggest that an acceleration mechanism is needed to explain the Hubble-type expansion found in HST observations. △ Less

Submitted 29 May, 2014; originally announced May 2014.

Comments: 5 pages, 3 figures. Accepted by MNRAS

On the sensitivity of extrasolar mass-loss rate ranges: HD 209458b a case study

Authors: C. S. Villarreal D'Angelo, E. M. Schneiter, A. Costa, P. Velázquez, A. Raga, A. Esquivel

Abstract: We present a 3D hydrodynamic study of the effects that different stellar wind conditions and planetary wind structures have on the calculated Ly-$α$ absorptions produced during the transit of HD 209458b. Considering a range of stellar wind speeds $\sim$[350-800] km s$^{-1}$, coronal temperature $\sim$[3-7] $\times10^{6}$ K and two values of the polytropic index $Γ$ $\sim$[1.01-1.13], while keeping… ▽ More We present a 3D hydrodynamic study of the effects that different stellar wind conditions and planetary wind structures have on the calculated Ly-$α$ absorptions produced during the transit of HD 209458b. Considering a range of stellar wind speeds $\sim$[350-800] km s$^{-1}$, coronal temperature $\sim$[3-7] $\times10^{6}$ K and two values of the polytropic index $Γ$ $\sim$[1.01-1.13], while keeping fixed the stellar mass loss rate, we found a that a $\dot M_p$ range between $\sim$[3-5] $\times 10^{10}$g s$^{-1}$ give account for the observational absorption in Ly-$α$ measured for the planetary system. Also, several models with anisotropic evaporation profiles for the planetary escaping atmosphere were carried out, showing that both, the escape through polar regions and through the night side yields larger absorptions than an isotropic planetary wind. △ Less

Submitted 20 December, 2013; originally announced December 2013.

A wind-shell interaction model for multipolar planetary nebulae

Authors: W. Steffen, N. Koning, A. Esquivel, G. Garcia-Segura, Ma. T. Garcia-Diaz, J. A. Lopez, M. Magnor

Abstract: We explore the formation of multipolar structures in planetary and pre-planetary nebulae from the interaction of a fast post-AGB wind with a highly inhomogeneous and filamentary shell structure assumed to form during the final phase of the high density wind. The simulations were performed with a new hydrodynamics code integrated in the interactive framework of the astrophysical modeling package SH… ▽ More We explore the formation of multipolar structures in planetary and pre-planetary nebulae from the interaction of a fast post-AGB wind with a highly inhomogeneous and filamentary shell structure assumed to form during the final phase of the high density wind. The simulations were performed with a new hydrodynamics code integrated in the interactive framework of the astrophysical modeling package SHAPE. In contrast to conventional astrophysical hydrodynamics software, the new code does not require any programming intervention by the user for setting up or controlling the code. Visualization and analysis of the simulation data has been done in SHAPE without external software. The key conclusion from the simulations is that secondary lobes in planetary nebulae, such as Hubble 5 and K3-17, can be formed through the interaction of a fast low-density wind with a complex high density environment, such as a filamentary circumstellar shell. The more complicated alternative explanation of intermittent collimated outflows that change direction, in many cases may therefore not be necessary. We consider that the wind-shell interaction scenario is more likely since the bow-shock shape expected from a strongly cooling bow-shock from jets is different from that of the observed bubbles. Furthermore, the timescales of the wind-wind interaction suggest that the progenitor star was rather massive. △ Less

Submitted 20 August, 2013; originally announced August 2013.

Comments: 10 pages, 11 figures, 2 color figures, accepted for publication in MNRAS

The photoevaporation of a neutral structure by an EUV+FUV radiation field

Authors: V. Lora, M. J. Vasconcelos, A. C. Raga, A. Esquivel, A. H. Cerqueira

Abstract: The EUV photoionizing radiation and FUV dissociating radiation from newly born stars photoevaporate their parental neutral cloud, leading to the formation of dense clumps that could eventually form additional stars. We study the effects of including a photodissociating FUV flux in models of the fragmentation of a photoevaporating, self-gravitating molecular cloud. We compute 3D simulations of the… ▽ More The EUV photoionizing radiation and FUV dissociating radiation from newly born stars photoevaporate their parental neutral cloud, leading to the formation of dense clumps that could eventually form additional stars. We study the effects of including a photodissociating FUV flux in models of the fragmentation of a photoevaporating, self-gravitating molecular cloud. We compute 3D simulations of the interaction of an inhomogeneous, neutral, self-gravitating cloud with external EUV and FUV radiation fields, and calculate the number of collapsing clumps and their mass. We find that the presence of an outer photodissociation region has an important effect on the formation of dense structures due to the expansion of an HII region. In particular, including a FUV field leads to the earlier formation of a larger number of dense clumps, which might lead to the formation of more stars. △ Less

Submitted 21 June, 2013; originally announced June 2013.

Comments: 10 pages, 10 figures. Accepted for publication in Rev Mex AA

Magnetization of cloud cores and envelopes and other observational consequences of reconnection diffusion

Authors: A. Lazarian, A. Esquivel, R. Crutcher

Abstract: Recent observational results for magnetic fields in molecular clouds reviewed by Crutcher (2012) seem to be inconsistent with the predictions of the ambipolar diffusion theory of star formation. These include the measured decrease in mass to flux ratio between envelopes and cores, the failure to detect any self-gravitating magnetically subcritical clouds, the determination of the flat PDF of the t… ▽ More Recent observational results for magnetic fields in molecular clouds reviewed by Crutcher (2012) seem to be inconsistent with the predictions of the ambipolar diffusion theory of star formation. These include the measured decrease in mass to flux ratio between envelopes and cores, the failure to detect any self-gravitating magnetically subcritical clouds, the determination of the flat PDF of the total magnetic field strengths implying that there are many clouds with very weak magnetic fields, and the observed scaling $B \propto ρ^{2/3}$ that implies gravitational contraction with weak magnetic fields. We consider the problem of magnetic field evolution in turbulent molecular clouds and discuss the process of magnetic field diffusion mediated by magnetic reconnection. For this process that we termed "reconnection diffusion" we provide a simple physical model and explain that this process is inevitable in view of the present day understanding of MHD turbulence. We address the issue of the expected magnetization of cores and envelopes in the process of star formation and show that reconnection diffusion provides an efficient removal of magnetic flux that depends only on the properties of MHD turbulence in the core and the envelope. As a result, the magnetic flux trapped during the collapse in the envelope is being released faster than the flux trapped in the core, resulting in much weaker fields in envelopes than in cores, as observed. We provide simple semi-analytical model calculations which support this conclusion and qualitatively agree with the observational results. We argue that magnetic reconnection provides a solution to the magnetic flux problem of star formation that agrees better with observations than the long-standing ambipolar diffusion paradigm. △ Less

Submitted 2 August, 2012; v1 submitted 20 June, 2012; originally announced June 2012.

Comments: 22 pages, 11 figures, submitted to ApJ, revised version. arXiv admin note: substantial text overlap with arXiv:1111.0694

Gravitational drag on a point mass in hypersonic motion through a gaseous medium

Authors: J. Canto, A. C. Raga, A. Esquivel, F. J. Sanchez-Salcedo

Abstract: We explore a ballistic orbit model to infer the gravitational drag force on an accreting point mass M, such as a black hole, moving at a hypersonic velocity v_{0} through a gaseous environment of density ρ_{0}. The streamlines blend in the flow past the body and transfer momentum to it. The total drag force acting on the body, including the nonlinear contribution of those streamlines with small im… ▽ More We explore a ballistic orbit model to infer the gravitational drag force on an accreting point mass M, such as a black hole, moving at a hypersonic velocity v_{0} through a gaseous environment of density ρ_{0}. The streamlines blend in the flow past the body and transfer momentum to it. The total drag force acting on the body, including the nonlinear contribution of those streamlines with small impact parameter that bend significantly and pass through a shock, can be calculated by imposing conservation of momentum. In this fully analytic approach, the ambiguity in the definition of the lower cut-off distance $r_{\rm min}$ in calculations of the effect of dynamical friction is removed. It turns out that $r_{\rm min}=\sqrt{e}GM/2v_{0}^{2}$. Using spherical surfaces of control of different sizes, we carry out a successful comparison between the predicted drag force and the one obtained from a high resolution, axisymmetric, isothermal flow simulation. We demonstrate that ballistic models are reasonably successful in accounting for both the accretion rate and the gravitational drag. △ Less

Submitted 15 August, 2011; originally announced August 2011.

Comments: 8 pages, 6 figures, accepted to MNRAS

Velocity Anisotropy as a Diagnostic of the Magnetization of the Interstellar Medium and Molecular clouds

Authors: Alejandro Esquivel, Alex Lazarian

Abstract: We use a set of magnetohydrodynamics (MHD) simulations of fully-developed (driven) turbulence to study the anisotropy in the velocity field that is induced by the presence of the magnetic field. In our models we study turbulence characterized by sonic Mach numbers M_s from 0.7 to 7.5, and Alfven Mach numbers M_A from 0.4 to 7.7. These are used to produce synthetic observations (centroid maps) that… ▽ More We use a set of magnetohydrodynamics (MHD) simulations of fully-developed (driven) turbulence to study the anisotropy in the velocity field that is induced by the presence of the magnetic field. In our models we study turbulence characterized by sonic Mach numbers M_s from 0.7 to 7.5, and Alfven Mach numbers M_A from 0.4 to 7.7. These are used to produce synthetic observations (centroid maps) that are analyzed. To study the effect of large scale density fluctuations and of white noise we have modified the density fields and obtained new centroid maps, which are analyzed. We show that restricting the range of scales at which the anisotropy is measured makes the method robust against such fluctuations. We show that the anisotropy in the structure function of the maps reveals the direction of the magnetic field for M_A \lesssim 1.5, regardless of the sonic Mach number. We found that the degree of anisotropy can be used to determine the degree of magnetization (i.e. M_A) for M_A \lesssim 1.5. To do this, one needs an additional measure of the sonic Mach number and an estimate of the LOS magnetic field, both feasible by other techniques, offering a new opportunity to study the magnetization state of the interstellar medium. △ Less

Submitted 2 August, 2011; originally announced August 2011.

Comments: Accepted for publication in The Astrophysical Journal (8 pages, 4 figures, one in color)

The cometary cavity created by an aligned streaming environment/collimated outflow interaction

Authors: D. López-Cámara, A. Esquivel, J. Cantó, A. C. Raga, P. F Velázquez, A. Rodríguez-González

Abstract: We present a "thin shell" model of the interaction of a biconical outflow and a streaming environment (aligned with the direction of the flow), as well as numerical (axisymmetric) simulations of such an interaction. A similar situation, although in a more complex setup, takes place at the head of the cometary structure of Mira. Thus, for most of the numerical simulations we explore parameters cons… ▽ More We present a "thin shell" model of the interaction of a biconical outflow and a streaming environment (aligned with the direction of the flow), as well as numerical (axisymmetric) simulations of such an interaction. A similar situation, although in a more complex setup, takes place at the head of the cometary structure of Mira. Thus, for most of the numerical simulations we explore parameters consistent with the observed bipolar outflow from Mira B. For these parameters, the interaction is non-radiative, so that a rather broad jet/streaming environment interaction region is formed. In spite of this, a reasonable agreement between the thin-shell analytic model and the numerical simulations is obtained. △ Less

Submitted 15 June, 2011; originally announced June 2011.

Comments: 9 pages, 11 figures, 1 table. Accepted for publication in RevMexAA

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

Authors: Ary Rodríguez-González, Pablo F. Velázquez, Margarita Rosado, Alejandro Esquivel, Jorge Reyes-Iturbide, J. Claudio Toledo-Roy

Abstract: 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ú-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… ▽ More 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ú-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. △ Less

Submitted 14 February, 2011; originally announced February 2011.

Comments: 9 pages, 11 figures, ApJ accepted

Mass and metal ejection efficiency in disk galaxies driven by young stellar clusters of nuclear starburst

Authors: A. Rodriguez-Gonzalez, A. Esquivel, A. C. Raga, P. Colin

Abstract: We present results from models of galactic winds driven by energy injected by nuclear starbursts. The total energy of the starburst is provided by young central stellar clusters and parts of the galactic interstellar medium are pushed out as part of the galactic wind (in some cases the galactic wind contains an important part of the metals produced in the new generation of stars). We have performe… ▽ More We present results from models of galactic winds driven by energy injected by nuclear starbursts. The total energy of the starburst is provided by young central stellar clusters and parts of the galactic interstellar medium are pushed out as part of the galactic wind (in some cases the galactic wind contains an important part of the metals produced in the new generation of stars). We have performed adiabatic and radiative 3D N-Body/Smooth Particle Hydrodynamics simulations of galactic winds using the GADGET-2 code. The numerical models cover a wide range of starburst (from $\sim10^2$ to $\sim10^7$ M$_\odot$) and galactic gas masses (from $\sim6\times10^6$ to $\sim10^{11}$ M$_\odot$). The concentrated central starburst regions are an efficient engine for producing of the mass and metal loss in galaxies, and also for driving the metal redistribution in the galaxies. △ Less

Submitted 1 February, 2011; originally announced February 2011.

Comments: 15 pages, 7 figures, (in press)

Journal ref: Revista Mexicana de Astronomía y Astrofísica, 2011

A model of Mira's cometary head/tail entering the Local Bubble

Authors: A. Esquivel, A. C. Raga, J. Canto, A. Rodriguez-Gonzalez, D. Lopez-Camara, P. F. Velazquez, F. De Colle

Abstract: We model the cometary structure around Mira as the interaction of an AGB wind from Mira A, and a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment Mira entered the Local Bubble (low density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well colli… ▽ More We model the cometary structure around Mira as the interaction of an AGB wind from Mira A, and a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment Mira entered the Local Bubble (low density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well collimated for a 100 kyr timescale. The result is a broad-head/narrow-tail structure that resembles the observed morphology of Mira's comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail. △ Less

Submitted 15 October, 2010; originally announced October 2010.

Comments: 12 pages, 8 figures (4 in color). Accepted for publication in The Astrophysical Journal

The angular momentum of condensations within elephant trunks

Authors: V. Lora, A. C. Raga, A. Esquivel

Abstract: The radiation from newly born stars photoevaporates their parental neutral cloud, leading to the formation of dense clumps that will eventually form stars. We present 3D simulations of the interaction of a neutral cloud with an external ionising radiation field, and compute the angular momenta of these collapsing clumps. The angular momenta of these collapsing clumps show that they have preferen… ▽ More The radiation from newly born stars photoevaporates their parental neutral cloud, leading to the formation of dense clumps that will eventually form stars. We present 3D simulations of the interaction of a neutral cloud with an external ionising radiation field, and compute the angular momenta of these collapsing clumps. The angular momenta of these collapsing clumps show that they have preferential orient mostly perpendicular to the direction of the incident ionising photon field. Therefore, the axes of the jet systems that will be eventually ejected (from the star + accretion disk systems that will form) will be oriented approximately perpendicular to the direction to the photoionising source. △ Less

Submitted 23 June, 2009; originally announced June 2009.

Comments: 6 pages, 7 figures, accepted in Astronomy and Astrophysics

An upper limit on the mass of the black hole in Ursa Minor dwarf galaxy

Authors: V. Lora, F. J. Sanchez-Salcedo, A. Raga, A. Esquivel

Abstract: The well-established correlations between the mass of massive black holes (BHs) in the nuclei of most studied galaxies and various global properties of their hosting galaxy lend support to the idea that dwarf galaxies and globular clusters could also host a BH in their centers. Direct kinematic detection of BHs in dwarf spheroidal (dSph) galaxies are seriously hindered by the small number of sta… ▽ More The well-established correlations between the mass of massive black holes (BHs) in the nuclei of most studied galaxies and various global properties of their hosting galaxy lend support to the idea that dwarf galaxies and globular clusters could also host a BH in their centers. Direct kinematic detection of BHs in dwarf spheroidal (dSph) galaxies are seriously hindered by the small number of stars inside the gravitational influence region of the BH. The aim of this Letter is to establish an upper dynamical limit on the mass of the putative BH in the Ursa Minor (UMi) dSph galaxy. We present direct N-body simulations of the tidal disruption of the dynamical fossil observed in UMi, with and without a massive BH. We find that the observed substructure is incompatible with the presence of a massive BH of (2-3)x10^4 Msun within the core of UMi. These limits are consistent with the extrapolation of the M_{BH}-sigma relation to the M_{BH}<10^6 Msun regime. We also show that the BH may be off-center with respect to the center of symmetry of the whole galaxy. △ Less

Submitted 4 June, 2009; originally announced June 2009.

Comments: 6 pages, 3 figures, ApJL (in press)

Journal ref: Astrophys.J.699:L113-L117,2009

Tsallis statistics as a tool for studying interstellar turbulence

Authors: A. Esquivel, A. Lazarian

Abstract: We used magnetohydrodynamic (MHD) simulations of interstellar turbulence to study the probability distribution functions (PDFs) of increments of density, velocity, and magnetic field. We found that the PDFs are well described by a Tsallis distribution, following the same general trends found in solar wind and Electron MHD studies. We found that the PDFs of density are very different in subsonic… ▽ More We used magnetohydrodynamic (MHD) simulations of interstellar turbulence to study the probability distribution functions (PDFs) of increments of density, velocity, and magnetic field. We found that the PDFs are well described by a Tsallis distribution, following the same general trends found in solar wind and Electron MHD studies. We found that the PDFs of density are very different in subsonic and supersonic turbulence. In order to extend this work to ISM observations we studied maps of column density obtained from 3D MHD simulations. From the column density maps we found the parameters that fit to Tsallis distributions and demonstrated that these parameters vary with the sonic and Alfvén Mach numbers of turbulence. This opens avenues for using Tsallis distributions to study the dynamical and perhaps magnetic states of interstellar gas. △ Less

Submitted 9 December, 2009; v1 submitted 15 May, 2009; originally announced May 2009.

Comments: 8 pages, 5 figures, 1 table. Accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.710:125-132,2010

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

Authors: J. Reyes-Iturbide, P. F. Velazquez, M. Rosado, A. Rodriguez-Gonzalez, R. F. Gonzalez, A. Esquivel

Abstract: 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… ▽ More 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. △ Less

Submitted 14 January, 2009; originally announced January 2009.

Comments: 8 pages, 6 figures, accepted for publication in MNRAS

Statistics of Centroids of Velocity

Authors: A. Esquivel, A. Lazarian

Abstract: We review the use of velocity centroids statistics to recover information of interstellar turbulence from observations. Velocity centroids have been used for a long time now to retrieve information about the scaling properties of the turbulent velocity field in the interstellar medium. We show that, while they are useful to study subsonic turbulence, they do not trace the statistics of velocity… ▽ More We review the use of velocity centroids statistics to recover information of interstellar turbulence from observations. Velocity centroids have been used for a long time now to retrieve information about the scaling properties of the turbulent velocity field in the interstellar medium. We show that, while they are useful to study subsonic turbulence, they do not trace the statistics of velocity in supersonic turbulence, because they are highly influenced by fluctuations of density. We show also that for sub-Alfvénic turbulence (both supersonic and subsonic) two-point statistics (e.g. correlation functions or power-spectra) are anisotropic. This anisotropy can be used to determine the direction of the mean magnetic field projected in the plane of the sky. △ Less

Submitted 13 January, 2009; originally announced January 2009.

Comments: To appear in proceedings of "Magnetic Fields in the Universe II"

The dynamics of internal working surfaces in MHD jets

Authors: Fabio De Colle, Alejandro Raga, Alejandro Esquivel

Abstract: The dynamical effects of magnetic fields in models of radiative, Herbig-Haro (HH) jets have been studied in a number of papers. For example, magnetized, radiative jets from variable sources have been studied with axisymmetric and 3D numerical simulations. In this paper, we present an analytic model describing the effect of a toroidal magnetic field on the internal working surfaces that result fr… ▽ More The dynamical effects of magnetic fields in models of radiative, Herbig-Haro (HH) jets have been studied in a number of papers. For example, magnetized, radiative jets from variable sources have been studied with axisymmetric and 3D numerical simulations. In this paper, we present an analytic model describing the effect of a toroidal magnetic field on the internal working surfaces that result from a variability in the ejection velocity. We find that for parameters appropriate for HH jets the forces associated with the magnetic field dominate over the gas pressure force within the working surfaces. Depending on the ram pressure radial cross section of the jet, the magnetic field can produce a strong axial pinch, or, alternatively, a broadening of the internal working surfaces. We check the validity of the analytic model with axisymmetric numerical simulations of variable, magnetized jets. △ Less

Submitted 18 August, 2008; originally announced August 2008.

Comments: 14 pages, 4 figures. ApJ in press

Filaments in Galactic Winds Driven by Young Stellar Clusters

Authors: Ary Rodriguez-Gonzalez, Alejandro Esquivel, Pablo Velazquez, Alejandro Raga, Veronica Melo

Abstract: The starburst galaxy M82 shows a system of H$α$-emitting filaments which extend to each side of the galactic disk. We model these filaments as the result of the interaction between the winds from a distribution of Super Stellar Clusters (SSCs). We first derive the condition necessary for producing a radiative interaction between the cluster winds (a condition which is met by the SSC distribution… ▽ More The starburst galaxy M82 shows a system of H$α$-emitting filaments which extend to each side of the galactic disk. We model these filaments as the result of the interaction between the winds from a distribution of Super Stellar Clusters (SSCs). We first derive the condition necessary for producing a radiative interaction between the cluster winds (a condition which is met by the SSC distribution of M82). We then compute 3D simulations for SSC wind distributions which satisfy the condition for a radiative interaction, and also for distributions which do not satisfy this condition. We find that the highly radiative models, that result from the interaction of high metallicity cluster winds, produce a structure of H$α$ emitting filaments, which qualitatively agrees with the observations of the M82, while the non-radiative SSC wind interaction models do not produce filamentary structures. Therefore, our criterion for radiative interactions (which depends on the mass loss rate and the terminal velocity of the SSC winds, and the mean separation between SSCs) can be used to predict whether or not an observed galaxy should have associated H$α$ emitting filaments. △ Less

Submitted 12 August, 2008; originally announced August 2008.

Comments: 10 pages, 6 Figures. ApJ Accepted, August 7, 2008