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.

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$\alpha$ 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.

Lyman $\alpha$ 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$\alpha$ 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.

Lyman $\alpha$ 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 $\alpha$ 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 $\alpha$ 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 $\alpha$ 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.

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.

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.

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.

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.

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.

We have analyzed four epochs of H$\alpha$ 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$\alpha$ intensities and the [S~II]/H$\alpha$ 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$\alpha$ ratio. Knot A has dramatically faded, and at the same time has had a substantial increase in its [S~II]/H$\alpha$ ratio. Possible interpretations of these results are discussed.

We present an analysis of four epochs of H$\alpha$ and [S II] $\lambda\lambda$ 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$\alpha$ 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] $\lambda$ 5007 (1994 and 2014), in which we see a clear drop in the value of the [O III]/H$\alpha$ ratio.

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

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.

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.

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.

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.

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.

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$\alpha$ 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.

We present a 3D hydrodynamic study of the effects that different stellar wind conditions and planetary wind structures have on the calculated Ly-$\alpha$ 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 $\Gamma$ $\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-$\alpha$ 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.

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.

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.

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 \rho^{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.

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 \rho_{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.

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.

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.

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

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.

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.

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\'en Mach numbers of turbulence. This opens avenues for using Tsallis distributions to study the dynamical and perhaps magnetic states of interstellar gas.

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.

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.

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.

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\'enic 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.

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.

The starburst galaxy M82 shows a system of H$\alpha$-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$\alpha$ 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$\alpha$ emitting filaments.

We explore the dynamics of a ``cluster wind'' flow in the regime in which the shocks resulting from the interaction of winds from nearby stars are radiative. We first show that for a cluster with T Tauri stars and/or Herbig Ae/Be stars, the wind interactions are indeed likely to be radiative. We then compute a set of four, three dimensional, radiative simulations of a cluster of 75 young stars, exploring the effects of varying the wind parameters and the density of the initial ISM that permeates the volume of the cluster. These simulations show that the ISM is compressed by the action of the winds into a structure of dense knots and filaments. The structures that are produced resemble in a qualitative way the observations of the IRAS 18511+0146 of Vig et al. (2007).

We present a 3D numerical simulation of the recently discovered cometary structure produced as Mira travels through the galactic ISM. In our simulation, we consider that Mira ejects a steady, latitude-dependent wind, which interacts with a homogeneous, streaming environment. The axisymmetry of the problem is broken by the lack of alignment between the direction of the relative motion of the environment and the polar axis of the latitude-dependent wind. With this model, we are able to produce a cometary head with a ``double bow shock'' which agrees well with the structure of the head of Mira's comet. We therefore conclude that a time-dependence in the ejected wind is not required for reproducing the observed double bow shock.

The bipolar morphology of the planetary nebula (PN) K 3-35 observed in radio-continuum images was modelled with 3D hydrodynamic simulations with the adaptive grid code yguazu-a. We find that the observed morphology of this PN can be reproduced considering a precessing jet evolving in a dense AGB circumstellar medium, given by a mass loss rate \dot{M}_{csm}=5x10^{-5}M_{\odot}/yr and a terminal velocity v_{w}=10 km/s. Synthetic thermal radio-continuum maps were generated from numerical results for several frequencies. Comparing the maps and the total fluxes obtained from the simulations with the observational results, we find that a model of precessing dense jets, where each jet injects material into the surrounding CSM at a rate \dot{M}_j=2.8x10^{-4} {M_{\odot}/yr (equivalent to a density of 8x10^{4} {cm}^{-3}, a velocity of 1500 km/s, a precession period of 100 yr, and a semi-aperture precession angle of 20 degrees agrees well with the observations.

We present analytic and numerical models of the `cluster wind' resulting from the multiple interactions of the winds ejected by the stars of a dense cluster of massive stars. We consider the case in which the distribution of stars (i.e., the number of stars per unit volume) within the cluster is spherically symmetric, has a power-law radial dependence, and drops discontinuously to zero at the outer radius of the cluster. We carry out comparisons between an analytic model (in which the stars are considered in terms of a spatially continuous injection of mass and energy) and 3D gasdynamic simulations (in which we include 100 stars with identical winds, located in 3D space by statistically sampling the stellar distribution function). From the analytic model, we find that for stellar distributions with steep enough radial dependencies the cluster wind flow develops a very high central density and a non-zero central velocity, and for steeper dependencies it becomes fully supersonic throughout the volume of the cluster (these properties are partially reproduced by the 3D numerical simulations). Therefore, the wind solutions obtained for stratified clusters can differ dramatically from the case of a homogeneous stellar distribution (which produces a cluster wind with zero central velocity, and a fully subsonic flow within the cluster radius). Finally, from our numerical simulations we compute predictions of X-ray emission maps and luminosities, which can be directly compared with observations of cluster wind flows.

We present models of turbulent mixing at the boundaries between hot (T~10^{6-7} K) and warm material (T~10^4 K) in the interstellar medium, using a three-dimensional magnetohydrodynamical code, with radiative cooling. The source of turbulence in our simulations is a Kelvin-Helmholtz instability, produced by shear between the two media. We found, that because the growth rate of the large scale modes in the instability is rather slow, it takes a significant amount of time (~1 Myr) for turbulence to produce effective mixing. We find that the total column densities of the highly ionized species (C IV, N V, and O VI) per interface (assuming ionization equilibrium) are similar to previous steady-state non-equilibrium ionization models, but grow slowly from log N ~10^{11} to a few 10^{12} cm^{-2} as the interface evolves. However, the column density ratios can differ significantly from previous estimates, with an order of magnitude variation in N(C IV)/N(O VI) as the mixing develops.

The investigation of the statistical properties of maps of line centroids has been used for almost 50 years, but there is still no general agreement on their interpretation. We try to quantify which properties of underlying turbulent velocity fields can be derived from centroid velocity maps, and we test conditions under which the scaling behaviour of the centroid velocities matches the scaling of the three-dimensional velocity field. Using fractal cloud models we study systematically the relation between three-dimensional density and velocity fields and the statistical properties of the produced line centroid maps. We put special attention to cases with large density fluctuations resembling supersonic interstellar turbulence. Starting from the Delta-variance analysis we derive a new tool to compute the scaling behaviour of the three-dimensional velocity field from observed intensity and centroid velocity maps. We provide two criteria to decide whether the information from the centroid velocities directly reflects the properties of the underlying velocity field. Applying these criteria allows to understand the different results found so far in the literature on the interpretation of the statistics of velocity centroids. A new iteration scheme can be used to derive the three-dimensional velocity scaling from centroid velocity maps for arbitrary density and velocity fields, but it requires an accurate knowledge of the average density of the considered interstellar cloud.

We use the results of magnetohydrodynamic (MHD) simulations to emulate spectroscopic observations and use maps of centroids to study their statistics. In order to assess under which circumstances the scaling properties of the velocity field can be retrieved from velocity centroids, we compare two point statistics (structure functions and power-spectra) of velocity centroids with those of the underlying velocity field and analytic predictions presented in a previous paper (Lazarian & Esquivel 2003). We tested a criterion for recovering velocity spectral index from velocity centroids derived in our previous work, and propose an approximation of the early criterion using only the variances of ``unnormalized'' velocity centroids and column density maps. It was found that both criteria are necessary, however not sufficient to determine if the centroids recover velocity statistics. Both criteria are well fulfilled for subsonic turbulence. We find that for supersonic turbulence with sonic Mach numbers > 2.5 centroids fail to trace the spectral index of velocity. Asymptotically, however, we claim that recovery of velocity statistics is always possible provided that the density spectrum is steep and the observed inertial range is sufficiently extended. In addition, we show that velocity centroids are useful for anisotropy studies and determining the direction of magnetic field, even if the turbulence is highly supersonic, but only if it is sub-Alfvenic. This provides a tool for mapping the magnetic field direction, and testing whether the turbulence is sub-Alfvenic or super-Alfvenic.

We address the problem of studying interstellar turbulence using spectral line data. We find a criterion when the velocity centroids may provide trustworthy velocity statistics. To enhance the scope of centroids applications, we construct a measure that we term ``modified velocity centroids'' (MVCs) and derive an analytical solution that relates the 2D spectra of the modified centroids with the underlying 3D velocity spectrum. We test our results using synthetic maps constructed with data obtained through simulations of compressible magnetohydrodynamical (MHD) turbulence. We show that the modified velocity centroids (MVCs) are complementary to the the Velocity Channel Analysis (VCA) technique. Employed together, they make determining of the velocity spectral index more reliable and for wider variety of astrophysical situations.

In a previous work Lazarian and Pogosyan suggested a technique to extract velocity and density statistics, of interstellar turbulence, by means of analysing statistics of spectral line data cubes. In this paper we test that technique, by studying the effect of correlation between velocity and density fields, providing a systematic analysis of the uncertainties arising from the numerics, and exploring the effect of a linear shear. We make use of both compressible MHD simulations and synthetic data to emulate spectroscopic observations and test the technique. With the same synthetic spectroscopic data, we also studied anisotropies of the two point statistics and related those anisotropies with the magnetic field direction. This presents a new technique for magnetic field studies. The results show that the velocity and density spectral indices measured are consistent with the analytical predictions. We identified the dominant source of error with the limited number of data points along a given line of sight. We decrease this type of noise by increasing the number of points and by introducing Gaussian smoothing. We argue that in real observations the number of emitting elements is essentially infinite and that source of noise vanishes.

HI data cubes are sources of unique information on interstellar turbulence. Doppler shifts due to supersonic motions contain information on turbulent velocity field which is otherwise difficult to obtain. However, the problem of separation of velocity and density fluctuations within HI data cubes is far from being trivial. Analytical description of the emissivity statistics of channel maps (velocity slices) in Lazarian & Pogosyan (2000) showed that the relative contribution of the density and velocity fluctuations depends on the thickness of the velocity slice. In particular, power-law assymptotics of the emissivity fluctuations change when the dispersion of the velocity at the scale under study becomes of the order of the velocity slice thickness (integrated width of the channel map). These results are the foundations of the Velocity-Channel Analysis (VCA) technique which allows to determine velocity and density statistics using 21-cm data cubes. The VCA has been successfully tested using data cubes obtained via compressible magnetohydrodynamic simulations and applied to Galactic and Magellanic Clouds data. As a tool it has become much more sophisticated recently when effects of absorption were accounted for. The systematic studies of vast 21-cm data sets to correlate the variations in the turbulence statistics with the astrophysical activity is on the agenda. This should allow to determine the interstellar energy injection mechanisms. Going beyond the VCA, we discuss other tools, namely, genus and anisotropy analysis. The first characterises the topology of HI, while the second provides magnetic field directions. We show a few applications of these new tools to HI data and MHD simulations.