BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

A. C. Raga

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.HEastro-phastro-ph.IMastro-ph.EPphysics.flu-dynastro-ph.CO

Catalog footprint

What is connected

37works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

37 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical simulations of the 1840s great eruption of $η$ Carinae: I. Revisiting the explosion scenario

In this work, we present new two-dimensional hydrodynamical simulations of the major eruption of $η$ Car in the 1840s, which resulted in the formation of a bipolar nebula that is commonly known as the large Homunculus. In our numerical models, we have included the high-speed component of 10000 km s$^{-1}$, detected in recent observations, that provides direct evidence of an explosive event. Here, we investigate whether such a violent explosion is able to explain both the shape and the dynamical evolution of $η$ Car's nebula. As in our previous work, we have assumed a two-stage scenario for $η$ Car's eruption: a slow outflow phase during a few decades before the eruption followed by the explosive event. From the collision of these outflow phases, the large Homunculus is produced. Our numerical simulations show that such scenario does not resemble some of the observed physical features and the expansion of the nebula. Notwithstanding, we also explore other injection parameters (mass-loss rate and ejection velocity) for these outflow phases. In particular, we find that an explosion with an intermediate-speed of 1000 km s$^{-1}$ is able to reproduce the morphology and the kinematical age of the large Homunculus.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The jet/counter-jet symmetry of the HH 212 outflow

We present Spitzer (IRAC) images observations and a VLT 2.1micron image of the HH 212 outflow. We find that this outflow has a strong symmetry, with jet/counterjet knot pairs with Delta x less than 1 arcsec position offsets. We deduce that the jet/counterjet knots are ejected with time differences Delta tau_0 approx. 6 yr and velocity differences Delta v_0~ 2 km/s. We also analyze the deviations of the knot positions perpendicular to the outflow axis, and interpret them in terms of a binary orbital motion of the outflow source. Through this model, we deduce a ~0.7M_solar mass for the outflow source, and a separation of ~80 AU between the components of the binary (assuming equal masses for the two components). Finally, using the IRAC data and the VLT 2.1micron image we have measured the proper motion velocities, obtaining values from 50 to 170km/s.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Observations and 3D Hydrodynamical models of planetary nebulae with Wolf Rayet type central stars

We present high-resolution, long-slit spectroscopic observations of two planetary nebulae with [WC] central stars located near the galactic bulge, M 1-32 and M 3-15. The observations were obtained with the 2.1-m telescope at the Observatorio Astronómico Nacional, San Pedro Mártir. M 1-32 shows wide wings on the base of its emission lines and M 3-15 has two very faint high-velocity knots. In order to model both PNe we built a three-dimensional model consisting of a jet interacting with an equatorially concentrated slow wind, emulating the presence of a dense torus, using the Yguazú hydrodynamical code. From our hydrodynamical models, we obtained position-velocity (PV) diagrams in the [N II]$λ$6583 line for comparison with the observations. We find that the spectral characteristics of M 1-32 and M 3-15 can be explained with the same physical model -a jet moving inside an AGB wind- using different parameters (physical conditions and position angles of the jet). In agreement with our model and observations, these objects contain a dense torus seeing pole-on and a bipolar jet escaping thorough the poles. Then we propose to classify this kind of objects as spectroscopic bipolar nebulae, although they have been classified morphologically as compact, round, or elliptical nebulae or with "close collimated lobes".

0 citations0 reviewsNo code linkedOpen access
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}}$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gemini-IFU spectroscopy of HH 111

We present new optical observations of the HH 111 Herbig-Haro jet using the Gemini Multi Object Spectrograph in its Integral Field Unit mode. Eight fields of 5" x 3.5" have been positioned along and across the HH 111 jet, covering the spatial region from knot E to L in HH 111 (namely, knots E, F, G, H, J, K and L). We present images and velocity channel maps for the [O I] 6300+6360, Halpha, [N II] 6548+6583 and [S II] 6716+6730 lines, as well as for the [S II]6716/6730 line ratio. We find that the HH 111 jet has an inner region with lower excitation and higher radial velocity, surrounded by a broader region of higher excitation and lower radial velocity. Also, we find higher electron densities at lower radial velocities. These results imply that the HH 111 jet has a fast, axial region with lower velocity shocks surrounded by a lower velocity sheath with higher velocity shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

StarBench: The D-type expansion of an HII region

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Stellar Wind Origin for the G2 Cloud: Three-Dimensional Numerical Simulations

We present 3D, adaptive mesh refinement simulations of G2, a cloud of gas moving in a highly eccentric orbit towards the galactic center. We assume that G2 originates from a stellar wind interacting with the environment of the Sgr A* black hole. The stellar wind forms a cometary bubble which becomes increasingly elongated as the star approaches periastron. A few months after periastron passage, streams of material begin to accrete on the central black hole with accretion rates $\dot{M} \sim 10^{-8}$ M$_\odot$ yr$^{-1}$. Predicted Br$γ$ emission maps and position-velocity diagrams show an elongated emission resembling recent observations of G2. A large increase in luminosity is predicted by the emission coming from the shocked wind region during periastron passage. The observations, showing a constant Br$γ$ luminosity, remain puzzling, and are explained here assuming that the emission is dominated by the free-wind region. The observed Br$γ$ luminosity ($\sim 8 \times 10^{30}$ erg s$^{-1}$) is reproduced by a model with a $v_w=50$ km s$^{-1}$ wind velocity and a $10^{-7}$ M$_\odot$ yr$^{-1}$ mass loss rate if the emission comes from the shocked wind. A faster and less dense wind reproduces the Br$γ$ luminosity if the emission comes from the inner, free wind region. The extended cometary wind bubble, largely destroyed by the tidal interaction with the black hole, reforms a few years after periastron passage. As a result, the Br$γ$ emission is more compact after periastron passage.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Age, size, and position of H ii regions in the Galaxy. Expansion of ionized gas in turbulent molecular clouds

This work aims at improving the current understanding of the interaction between H ii regions and turbulent molecular clouds. We propose a new method to determine the age of a large sample of OB associations by investigating the development of their associated H ii regions in the surrounding turbulent medium. Using analytical solutions, one-dimensional (1D), and three-dimensional (3D) simulations, we constrained the expansion of the ionized bubble depending on the turbulent level of the parent molecular cloud. A grid of 1D simulations was then computed in order to build isochrone curves for H ii regions in a pressure-size diagram. This grid of models allowed to date large sample of OB associations and was used on the H ii Region Discovery Survey (HRDS). Analytical solutions and numerical simulations showed that the expansion of H ii regions is slowed down by the turbulence up to the point where the pressure of the ionized gas is in a quasi-equilibrium with the turbulent ram pressure. Based on this result, we built a grid of 1D models of the expansion of H ii regions in a profile based on Larson laws. The 3D turbulence is taken into account by an effective 1D temperature profile. The ages estimated by the isochrones of this grid agree well with literature values of well-known regions such as Rosette, RCW 36, RCW 79, and M16. We thus propose that this method can be used to give ages of young OB associations through the Galaxy such as the HRDS survey and also in nearby extra-galactic sources.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An asymmetric jet launching model for the protoplanetary nebula CRL 618

We propose an asymmetrical jet ejection mechanism in order to model the mirror symmetry observed in the lobe distribution of some protoplanetary nebulae (pPNe), such as the pPN CRL 618. 3D hydrodynamical simulations of a precessing jet launched from an orbiting source were carried out including an alternation in the ejections of the two outflow lobes, depending on which side of the precessing accretion disk is hit by the accretion column from a Roche lobe-filling binary companion. Both synthetic optical emission maps and position-velocity (PV) diagrams were obtained from the numerical results with the purpose of carrying out a direct comparison with observations. Depending on the observer's point of view, multipolar morphologies are obtained which exhibit a mirror symmetry at large distances from the central source. The obtained lobe sizes and their spatial distribution are in good agreement with the observed morphology of the pPN CRL 618. We also obtain that the kinematic ages of the fingers are similar to those obtained in the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collisionally excited filaments in HST H$α$ and H$β$ images of HH~1/2

We present new H$α$ and H$β$ images of the HH~1/2 system, and we find that the H$α$/H$β$ ratio has high values in ridges along the leading edges of the HH~1 bow shock and of the brighter condensations of HH~2. These ridges have H$α$/H$β=4\to 6$, which is consistent with collisional excitation from the $n=1$ to the $n=3$ and 4 levels of hydrogen in a gas of temperatures $T=1.5\to 10\times 10^4$~K. This is therefore the first direct proof that the collisional excitation/ionization region of hydrogen right behind Herbig-Haro shock fronts is detected.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interaction between the IGM and a dwarf galaxy

Dwarf Galaxies are the most common objects in the Universe and are believed to contain large amounts of dark matter. There are mainly three morphologic types of dwarf galaxies: dwarf ellipticals, dwarf spheroidals and dwarf irregulars. Dwarf irregular galaxies are particularly interesting in dwarf galaxy evolution, since dwarf spheroidal predecessors could have been very similar to them. Therefore, a mechanism linked to gas-loss in dwarf irregulars should be observed, i.e. ram pressure stripping. In this paper, we study the interaction between the ISM of a dwarf galaxy, and a flowing IGM. We derive the weak-shock, plasmon solution corresponding to the balance between the post-bow shock pressure and the pressure of the stratified ISM (which we assume follows the fixed stratification of a gravitationally dominant dark matter halo). We compare our model with previously published numerical simulations and with the observed shape of the HI cloud around the Ho II and Pegasus dwarf irregular galaxies. We show that such a comparison provides a straightforward way for estimating the Mach number of the impinging flow.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Proper Motions of Young Stellar Outflows in the Mid-Infrared with Spitzer. II. HH 377/Cep E

We have used multiple mid-infrared observations at 4.5 micron obtained with the Infrared Array Camera, of the compact (~1.4 arcmin) young stellar bipolar outflow Cep E to measure the proper motion of its brightest condensations. The images span a period of ~6 yr and have been reprocessed to achieve a higher angular resolution (~0.8 arcsec) than their normal beam (2 arcsec). We found that for a distance of 730 pc, the tangential velocities of the North and South outflow lobes are 62+/-29 and 94+/-6 km/s respectively, and moving away from the central source roughly along the major axis of the flow. A simple 3D hydrodynamical simulation of the H2 gas in a precessing outflow supports this idea. Observations and model confirm that the molecular Hydrogen gas, traced by the pure rotational transitions, moves at highly supersonic velocities without being dissociated. This suggests either a very efficient mechanism to reform H2 molecules along these shocks or the presence of some other mechanism (e.g. strong magnetic field) that shields the H2 gas.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational drag on a point mass in hypersonic motion within a Gaussian disk

We develop an analytical model for the accretion and gravitational drag on a point mass that moves hypersonically in the midplane of a gaseous disk with a Gaussian vertical density stratification. Such a model is of interest for studying the interaction between a planet and a protoplanetary disk, as well as the dynamical decay of massive black holes in galactic nuclei. The model considers that the flow is ballistic, and gives fully analytical expressions for both the accretion rate onto the point mass, and the gravitational drag it suffers. The expressions are further simplified by taking the limits of a thick, and of a thin disk. The results for the thick disk reduce correctly to those for a uniform density environment (Cantó et al. 2011). We find that for a thin disk (small vertical scaleheight compared to the gravitational radius) the accretion rate is proportional to the mass of the moving object and to the surface density of the disk, while the drag force is independent of the velocity of the object. The gravitational deceleration of the hypersonic perturber in a thin disk was found to be independent of its parameters (i.e. mass or velocity) and depends only on the surface mass density of the disk. The predictions of the model are compared to the results of three-dimensional hydrodynamical simulations, with a reasonable agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Proper Motions of Young Stellar Outflows in the Mid-Infrared with Spitzer (IRAC). I. The NGC 1333 region

We use two 4.5micron Spitzer (IRAC) maps of the NGC 1333 region taken over approx. 7 yr interval to determine proper motions of its associated outflows. This is a first, successful attempt at obtaining proper motions of stellars outflow from Spitzer observations. For the outflow formed by the Herbig-Haro objects HH7, 8 and 10, we find proper motions of approx. 9-13 km/s, which are consistent with previously determined optical proper motions of these objects. We determine proper motions for a total of 8 outflows, ranging from approx. 10 to 100 km/s. The derived proper motions show that out of these 8 outflows, 3 have tangential velocities less or equal to 20 km/s. This result shows that a large fraction of the observed outflows have low intrinsic velocities, and that the low proper motions are not merely a projection effect.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spitzer observations of the HH 1/2 system. The discovery of the counterjet

We present unpublished Spitzer IRAC observations of the HH 1/2 young stellar outflow processed with a high angular resolution deconvolution algorithm that produces sub-arcsecond (approx. 0.6" - 0.8") images. In the resulting mid-infrared images the optically invisible counterjet is detected for the first time. The counterjet is approximately half as bright as the jet at 4.5 micron (the IRAC band that best traces young stellar outflows) and has a length of approx. 10". The NW optical jet itself can be followed back in the mid-IR to the position of the exciting VLA 1 source. An analysis of the IRAC colors indicates that the jet/counterjet emission is dominated by collisionally excited H2 pure rotational lines arising from a medium with a neutral Hydrogen gas density of 1000-2000 per cubic cm and a temperature of 1500 K. The observed jet/counterjet brightness asymmetry is consistent with an intrinsically symmetric outflow with extinction from a dense, circumstellar structure of 6" size (along the outflow axis), and with a mean visual extinction of Av=11 mag.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radio continuum emission from knots in the DG Tau jet

Context: HH 158, the jet from the young star DG Tau, is one of the few sources of its type where jet knots have been detected at optical and X-ray wavelengths. Aims: To search, using Very Large Array observations of this source, radio knots and if detected, compare them with the optical and X-ray knots. To model the emission from the radio knots. Methods: We analyzed archive data and also obtained new Very Large Array observations of this source, as well as an optical image, to measure the present position of the knots. We also modeled the radio emission from the knots in terms of shocks in a jet with intrinsically time-dependent ejection velocities. Results: We detected radio knots in the 1996.98 and 2009.62 VLA data. These radio knots are,within error, coincident with optical knots. We also modeled satisfactorily the observed radio flux densities as shock features from a jet with intrinsic variability. All the observed radio, optical, and X-ray knot positions can be intepreted as four successive knots, ejected with a period of 4.80 years and traveling away from the source with a velocity of 198 km s$^{-1}$ in the plane of the sky. Conclusions: The radio and optical knots are spatially correlated and our model can explain the observed radio flux densities. However, the X-ray knots do not appear to have optical or radio counterparts and their nature remains poorly understood.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The jet/counterjet IR symmetry of HH34 and the size of the jet formation region

We present a new Spitzer IRAC images of the HH 34 outflow. These are the first images that detect both the knots along the southern jet and the northern counterjet (the counterjet knots were only detected previously in a long slit spectrum). This result removes the problem of the apparent coexistence of a large scale symmetry (at distances of up to approx.1 pc) and a complete lack of symmetry close to the source (at distances of approx. 1.e+17 cm) for this outflow. We present a quantitative evaluation of the newly found symmetry between the HH 34 jet and counterjet, and show that the observed degree of symmetry implies that the jet production region has a characteristic size <2.8 AU. This is the strongest constraint yet derived for the size of the region in which HH jets are produced.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The precession of the HH 111 flow in the infrared

We present Spitzer IRAC images of the HH 111 outflow, that show a wealth of condensations/knots in both jet and counterjet. Studying the positional distribution of these knots, we find very suggestive evidence of a mirror symmetric pattern in the jet/counterjet flow. We model this pattern as the result of an orbital motion of the jet source around a binary companion. From a fit of an analytic, ballistic model to the observed path of the HH 111 system, we find that the motion in a binary with two approx. 1 Msolar stars (one of them being the HH 111 source), in a circular orbit with a separation of approx. 186 AU would produce the mirror symmetric pattern seen in the outflow.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A reassessment of the kinematics of PV Cephei based on accurate proper motion measurements

We present two Very Large Array observations of the pre-main-sequence star PV Cephei, taken with a separation of 10.5 years. These data show that the proper motions of this star are $μ_α\cos δ= +10.9 \pm 3.0$ mas yr$^{-1}$; $μ_δ= +0.2 \pm 1.8$ mas yr$^{-1}$, very similar to those -previously known- of HD 200775, the B2Ve star that dominates the illumination of the nearby reflection nebula NGC 7023. This result suggests that PV Cephei is not a rapidly moving run-away star as suggested by previous studies. The large velocity of PV Cephei had been inferred from the systematic eastward displacement of the bisectors of successive pairs of Herbig Haro knots along its flow. These systematic shifts might instead result from an intrinsic dissymmetry in the ejection mechanisms, or from an asymmetric distribution of the circumstellar material.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Side-entrainment in a jet embedded in a sidewind

Numerical simulations of HH jets never show side-entrainment of environmental material into the jet beam. This is because the bow shock associated with the jet head pushes the surrounding environment into a dense shell, which is never in direct contact with the sides of the jet beam. We present 3D simulations in which a side-streaming motion (representing the motion of the outflow source through the surrounding medium) pushes the post-bow shock shell into direct contact with the jet beam. This is a possible mechanism for modelling well collimated "molecular jets" as an atomic/ionic flow which entrains molecules initially present only in the surrounding environment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical Analysis of Water Masers in Star-Forming Regions: Cepheus A and W75 N

We have done a statistical analysis of Very Long Baseline Array (VLBA) data of water masers in the star-forming regions (SFRs) Cepheus A and W75 N, using correlation functions to study the spatial clustering and Doppler-velocity distribution of these masers. Two-point spatial correlation functions show a characteristic scale size for clusters of water maser spots < or ~1 AU, similar to the values found in other SFRs. This suggests that the scale for water maser excitation tends to be < or ~1 AU. Velocity correlation functions show power-law dependences with indices that can be explained by regular velocity fields, such as expansion and/or rotation. These velocity fields are similar to those indicated by the water maser proper-motion measurements; therefore, the velocity correlation functions appear to reveal the organized motion of water maser spots on scales larger than 1 AU.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Low and high velocity clouds produced by young stellar clusters

Intermediate and high velocity HI clouds rain onto the plane of our Galaxy. They are observed at heights of between 500 and 1500 pc, falling onto the Galactic plane at velocities from 50 to 140 km s$^{-1}$. To explain the origin of these clouds, we present a galactic fountain model, driven by the wind from a super stellar cluster (SSC). We solve the equations for a steady, radiative de Laval nozzle flow. We consider two effects not considered previously in astrophysical nozzle flow models: cooling functions for different metallicities, and the direct action of the galactic gravitational field on the gas flowing along the nozzle. For an adiabatic nozzle flow, the gravity acting directly on the gas within the nozzle "stalls" the nozzle flow for initial wind velocities lower than the escape velocity from the Galaxy. For the same wind velocity, a radiative nozzle flow stalls at lower altitudes above the galactic plane. We find that SSC winds with velocities of $v_w=500 - 800$ km s$^{-1}$ produce nozzles stall at heights of $x_m=1 - 15$ kpc. The stalled nozzle flow then rains back onto the galactic plane at velocities in the range observed in intermediate and high velocity HI clouds. We study a nozzle flow driven by a wind from a SSC close to the Galactic centre. We find that for velocities within the range expected for a SSC wind, we can produce nozzle flows that stall above the galactic plane. These stalled flows produce cool, infalling clouds with velocities similar to those of intermediate and high velocity HI clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Photoionized mixing layer models of the diffuse ionized gas

It is generally believed that O stars, confined near the galactic midplane, are somehow able to photoionize a significant fraction of what is termed the "diffuse ionized gas" (DIG) of spiral galaxies, which can extend up to 1-2 kpc above the galactic midplane. The heating of the DIG remains poorly understood, however, as simple photoionization models do not reproduce the observed line ratio correlations well or the DIG temperature. We present turbulent mixing layer models in which warm photoionized condensations are immersed in a hot supersonic wind. Turbulent dissipation and mixing generate an intermediate region where the gas is accelerated, heated and mixed. The emission spectrum of such layers are compared with observations of Rand (ApJ 462, 712) of the DIG in the edge-on spiral NGC2363. We generate two sequence of models that fit the line ratio correlations between [SII]/H-alpha, [OI]/H-alpha, [NII]/[SII] and [OIII]/H-beta reasonably well. In one sequence of models the hot wind velocity increases while in the other the ionization parameter and layer opacity increases. Despite the success of the mixing layer models, the overall efficiency in reprocessing the stellar UV is much too low, much less than 1%, which compels us to reject the TML model in its present form.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Revisiting 2D Numerical Models for the 19th century outbursts of $η$ Carinae

We present here new results of two-dimensional hydrodynamical simulations of the eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suffered by the massive star $η$ Car. The two bipolar nebulae commonly known as the Homunculus and the little Homunculus were formed from the interaction of these eruptive events with the underlying stellar wind. As in previous work (Gonzalez et al. 2004a, 2004b), we assume here an interacting, nonspherical multiple-phase wind scenario to explain the shape and the kinematics of both Homunculi, but adopt a more realistic parametrization of the phases of the wind. During the 1890s eruptive event, the outflow speed {\it decreased} for a short period of time. This fact suggests that the little Homunculus is formed when the eruption ends, from the impact of the post-outburst $η$ Car wind (that follows the 1890s event) with the eruptive flow (rather than by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous models; Gonzalez et al. 2004a, 2004b). Our simulations reproduce quite well the shape and the observed expansion speed of the large Homunculus. The little Homunculus (which is embedded within the large Homunculus) becomes Rayleigh-Taylor unstable and develop filamentary structures that resembles the spatial features observed in the polar caps. In addition, we find that the interior cavity between the two Homunculi is partially filled by material that is expelled during the decades following the great eruption. This result may be connected with the observed double-shell structure in the polar lobes of the $η$ Car nebula. Finally, as in previous work, we find the formation of tenuous, equatorial, high-speed features that seem to be related to the observed equatorial skirt of $η$ Car.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Shock excitation of the knots of Hen 3-1475

We present new optical STIS HST spectroscopic observations of the jets of the proto-planetary nebula Hen 3-1475. The excitation conditions of the knots of Hen 3-1475 are derived from the observed optical spectra, confirming that the knots are shock excited. The shocked spectra are qualitatively reproduced by simple ``3/2''D bow shock models. We present a set of bow shock models devoted to planetary nebulae, and discuss the effects of the pre-ionization conditions, the bow shock velocity, the bow shock shape and the chemical abundances on the predicted spectra. To explore the reliability of the ``3/2''D bow shock models, we also compare the observed spectra of other three proto-planetary nebulae (M 1-92, M 2-56 and CRL 618) to the predicted spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The transverse velocity and excitation structure of the HH 110 jet

We present long-slit spectroscopic observations of the HH 110 jet obtained with the 4.2~m William Herschel Telescope. We have obtained for the first time, spectra for slit positions along and across the jet axis (at the position of knots B, C, I, J and P) to search for the observational signatures of entrainment and turbulence by studying the kinematics and the excitation structure. We find that the HH 110 flow accelerates from a velocity of 35 km/s in knot A up to 110 km/s in knot P. We find some systematic trends for the variation of the emission line ratios along the jet. No clear trends for the variation of the radial velocity are seen across the width of the jet beam. The cross sections of the jet show complex radial velocity and line emission structures which differ quite strongly from each other.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Jet cocoons and the formation of Narrow Line Clouds in Seyfert galaxies

We present non-adiabatic hydrodynamic simulations of a supersonic light jet propagating into a fully ionized medium of uniform density on a scale representative of the narrow line region (NLR) in Seyfert galaxies with associated radio jets. In this regime the cooling distance of the swept up gas in the bowshock of the jet is of the same order as the transverse extent of the jet bowshock, as opposed to the more extreme regimes found for more powerful adiabatic large scale jets or the slow galactic jets which have been simulated previously. We calculate the emissivity for the H_alpha line and radio synchrotron emission. We find that the structure of the line emitting cold envelope of the jet cocoon is strongly dependent on the non-stationary dynamics of the jet head as it propagates through the ambient medium. We observe the formation of cloud-like high density regions which we associate with NLR clouds and filaments. We find that some of these clouds might be partially neutral and represent sites of jet induced star formation. The calculated H_alpha flux and the spectral line width are consistent with NLR observations. The simulation of the radio-optical emission with radiative cooling confirms the basic result of the geometric bowshock model developed by Taylor et al. (1989) that the start of noticeable optical line emission can be significantly offset from the hotspot of the radio emission. However, the time-dependent nature of the jet dynamics implies significant differences from their geometric bowshock model.

0 citations0 reviewsNo code linkedOpen access