Source author record

A. Esquivel

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2023arXiv

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

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.

preprint2022arXiv

Rectangular core-collapse supernova remnants: application to Puppis A

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.

preprint2021arXiv

Estimating the propagation of a uniformly accelerated jet

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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).

preprint2017arXiv

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

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.

preprint2016arXiv

Photo-ionization of planetary winds: case study HD209458b

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}}$.

preprint2016arXiv

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

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.

preprint2015arXiv

3D MHD simulation of polarized emission in SN 1006

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.

preprint2015arXiv

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

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.

preprint2014arXiv

A two-mode planetary nebula luminosity function

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2013arXiv

A wind-shell interaction model for multipolar planetary nebulae

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

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

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.

preprint2010arXiv

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

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.

preprint2009arXiv

The angular momentum of condensations within elephant trunks

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.

preprint2009arXiv

Tsallis statistics as a tool for studying interstellar turbulence

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.