The massive star-forming region W75N~(B) is thought to host a cluster of massive protostars (VLA~1, VLA~2, and VLA~3) undergoing different evolutionary stages. In this work, we present radio continuum data with the highest sensitivity and angular resolution obtained to date in this region, using the VLA-A and covering a wide range of frequencies (4-48~GHz), which allowed us to study the morphology and the nature of the emission of the different radio continuum sources. We also performed complementary studies with multi-epoch VLA data and ALMA archive data at 1.3 mm wavelength. We find that VLA~1 is driving a thermal radio jet at scales of $\approx$0.1 arcsec ($\approx$130 au), but also shows signs of an incipient hyper-compact HII region at scales of $\lesssim$ 1 arcsec ($\lesssim$ 1300~au). VLA~3 is also driving a thermal radio jet at scales of a few tenths of arcsec (few hundred of au). We conclude that this jet is shock-exciting the radio continuum sources Bc and VLA~4 (obscured HH objects), which show proper motions moving outward from VLA~3 at velocities of $\approx$112--118~km/s. We have also detected three new weak radio continuum sources, two of them associated with millimeter continuum cores observed with ALMA, suggesting that these two sources are also embedded YSOs in this massive star-forming region.

Orion BN/KL is an example of a poorly understood phenomena in star forming regions involving the close encounter of young stellar objects. The explosive structure, the great variety of molecules observed, the energy involved in the event and the mass of the region suggest a contribution in the {chemical diversity} of the local interstellar medium. Nevertheless, the frequency and duration of other events like this have not been determined. In this paper, we explore a recent analytic model that takes into account the interaction of a clump with its molecular environment. We show that the widespread kinematic ages of the Orion fingers -- 500 to 4000 years -- is a consequence of the interaction of the explosion debris with the surrounding medium. This model explains satisfactorily the age discrepancy of the Orion fingers, and infers the initial conditions together with the lifetime of the explosion. Moreover, our model can explain why some CO streamers do not have a H$_2$ finger associated.

We model molecular outflows produced by the time dependent interaction between a stellar wind and a rotating cloud envelope in gravitational collapse, studied by Ulrich. We consider spherical and anisotropic stellar winds. We assume that the bipolar outflow is a thin shocked shell, with axial symmetry around the cloud rotation axis and obtain the mass and momentum fluxes into the shell. We solve numerically a set of partial differential equations in space and time, and obtain the shape of the shell, the mass surface density, the velocity field, and the angular momentum of the material in the shell. We find that there is a critical value of the ratio between the wind and the accretion flow momentum rates $\beta$ that allows the shell to expand. As expected, the elongation of the shells increase with the stellar wind anisotropy. In our models, the rotation velocity of the shell is the order to 0.1 - 0.2 km s$^{-1}$, a factor of 5-10 lower than the values measured in several sources. We compare our models with those of Wilkin and Stahler for early evolutionary times and find that our shells have the same sizes at the pole, although we use different boundary conditions at the equator.

The interaction of a high velocity clump of gas has been described by the plasmon model, which considers balance between ram pressure and the internal stratified structure of the decelerated clump. In this paper we propose an analytical model to describe the mass loss of such a clump due the interaction with the environment, describing its influence on the plasmon dynamics. We carry out comparisons between an analytic model and axisymetric gasdynamic simulations of plasmon evolution. From our simulations we were able to find the values of the friction constants $\alpha$ and $\lambda$. Comparing with the complete analytic model from which we can infer the position and the mass loss of the clump as function of the clump's density and the environment ratio.

We present hydrodynamical simulations of interacting Coronal Mass Ejections in the Interplanetary medium (ICMEs). In these events, two consecutive CMEs are launched from the Sun in similar directions within an interval of time of a few hours. In our numerical model, we assume that the ambient solar wind is characterized by its velocity and mass-loss rate. Then, the CMEs are generated when the flow velocity and mass-loss rate suddenly change, with respect to the ambient solar wind conditions during two intervals of time, which correspond to the durations of the CMEs. After their interaction, a merged region is formed and evolve as a single structure into the interplanetary medium. In this work, we are interested in the general morphology of this merged region, which depends on the initial parameters of the ambient solar wind and each of the CMEs involved. In order to understand this morphology, we have performed a parametric study in which we characterize the effects of the initial parameters variations on the density and velocity profiles at 1 AU, using as reference the well-documented event of July 25th, 2004. Based on this parametrization we were able to reproduce with a high accuracy the observed profiles. Then, we apply the parametrization results to the interaction events of May 23, 2010; August 1, 2010; and November 9, 2012. With this approach and using values for the input parameters within the CME observational errors, our simulated profiles reproduce the main features observed at 1 AU. Even though we do not take into account the magnetic field, our models give a physical insight into the propagation and interaction of ICMEs.

We present Submillimeter Array (SMA) 1.35 mm subarcsecond angular resolution observations toward the LkH{\alpha} 234 intermediate-mass star-forming region. The dust emission arises from a filamentary structure of $\sim$5 arcsec ($\sim$4500 au) enclosing VLA 1-3 and MM 1, perpendicular to the different outflows detected in the region. The most evolved objects are located at the southeastern edge of the dust filamentary structure and the youngest ones at the northeastern edge. The circumstellar structures around VLA 1, VLA 3, and MM 1 have radii between $\sim$200 and $\sim$375 au and masses in the $\sim$0.08-0.3 M$_{\odot}$ range. The 1.35 mm emission of VLA 2 arises from an unresolved (r$< 135$ au) circumstellar disk with a mass of $\sim$0.02 M$_{\odot}$. This source is powering a compact ($\sim$4000 au), low radial velocity ($\sim$7 km s$^{-1}$) SiO bipolar outflow, close to the plane of the sky. We conclude that this outflow is the "large-scale" counterpart of the short-lived, episodic, bipolar outflow observed through H$_2$O masers at much smaller scales ($\sim $180 au), and that has been created by the accumulation of the ejection of several episodic collimated events of material. The circumstellar gas around VLA 2 and VLA 3 is hot ($\sim$130 K) and exhibits velocity gradients that could trace rotation. There is a bridge of warm and dense molecular gas connecting VLA 2 and VLA 3. We discuss the possibility that this bridge could trace a stream of gas between VLA 3 and VLA 2, increasing the accretion rate onto VLA 2 to explain why this source has an important outflow activity.

The current paradigm of star formation through accretion disks, and magnetohydrodynamically driven gas ejections, predicts the development of collimated outflows, rather than expansion without any preferential direction. We present radio continuum observations of the massive protostar W75N(B)-VLA 2, showing that it is a thermal, collimated ionized wind and that it has evolved in 18 years from a compact source into an elongated one. This is consistent with the evolution of the associated expanding water-vapor maser shell, which changed from a nearly circular morphology, tracing an almost isotropic outflow, to an elliptical one outlining collimated motions. We model this behavior in terms of an episodic, short-lived, originally isotropic, ionized wind whose morphology evolves as it moves within a toroidal density stratification.

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.

W75N(B) is a massive star-forming region that contains three radio continuum sources (VLA 1, VLA 2, and VLA 3), which are thought to be three massive young stellar objects at three different evolutionary stages. VLA 1 is the most evolved and VLA 2 the least evolved source. The 22 GHz H2O masers associated with VLA 1 and VLA 2 have been mapped at several epochs over eight years. While the H2O masers in VLA 1 show a persistent linear distribution along a radio jet, those in VLA 2 are distributed around an expanding shell. Furthermore, H2O maser polarimetric measurements revealed magnetic fields aligned with the two structures. Using new polarimetric observations of H2O masers, we aim to confirm the elliptical expansion of the shell-like structure around VLA 2 and, at the same time, to determine if the magnetic fields around the two sources have changed. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz H2O masers towards the massive star-forming region W75N(B). The H2O maser distribution around VLA 1 is unchanged from that previously observed. We made an elliptical fit of the H2O masers around VLA 2. We find that the shell-like structure is still expanding along the direction parallel to the thermal radio jet of VLA 1. While the magnetic field around VLA 1 has not changed in the past 7 years, the magnetic field around VLA 2 has changed its orientation according to the new direction of the major-axis of the shell-like structure and it is now aligned with the magnetic field in VLA 1.

We report multi-epoch VLBI H$_2$O maser observations towards the compact cluster of YSOs close to the Herbig Be star LkH$\alpha$ 234. This cluster includes LkH$\alpha$ 234 and at least nine more YSOs that are formed within projected distances of $\sim$10 arcsec ($\sim$9,000 au). We detect H$_2$O maser emission towards four of these YSOs. In particular, our VLBI observations (including proper motion measurements) reveal a remarkable very compact ($\sim$0.2 arcsec = $\sim$180 au), bipolar H$_2$O maser outflow emerging from the embedded YSO VLA 2. We estimate a kinematic age of $\sim$40 yr for this bipolar outflow, with expanding velocities of $\sim$20 km s$^{-1}$ and momentum rate $\dot M_w V_w$ $\simeq$ $10^{-4}-10^{-3}$ M$_{\odot}$ yr$^{-1}$ km s$^{-1}$$\times (\Omega$/$4\pi)$, powered by a YSO of a few solar masses. We propose that the outflow is produced by recurrent episodic jet ejections associated with the formation of this YSO. Short-lived episodic ejection events have previously been found towards high-mass YSOs. We show now that this behaviour is also present in intermediate-mass YSOs. These short-lived episodic ejections are probably related to episodic increases in the accretion rate, as observed in low-mass YSOs. We predict the presence of an accretion disk associated with VLA 2. If detected, this would represent one of the few known examples of intermediate-mass stars with a disk-YSO-jet system at scales of a few hundred au.

We present multi-epoch Very Long Baseline Array (VLBA) H$_2$O maser observations toward the massive young stellar objects (YSOs) VLA 2 and VLA 3 in the star-forming region AFGL 2591. Through these observations, we have extended the study of the evolution of the masers towards these objects up to a time span of $\sim$ 10 yrs, measuring their radial velocities and proper motions. The H$_2$O masers in VLA 3, the most massive YSO in AFGL 2591 ($\sim$ 30--40~M$_{\odot}$), are grouped within projected distances of $\lesssim$ 40 mas ($\lesssim$ 130 AU) from VLA 3. In contrast to other H$_2$O masers in AFGL 2591, the masers associated with VLA 3 are significantly blueshifted (up to $\sim$ 30 km s$^{-1}$) with respect to the velocity of the ambient molecular cloud. We find that the H$_2$O maser cluster as a whole, has moved westwards of VLA~3 between the 2001 and 2009 observations, with a proper motion of $\sim$ 1.2 mas yr$^{-1}$ ($\sim$ 20 km s$^{-1}$). We conclude that these masers are tracing blueshifted outflowing material, shock excited at the inner parts of a cavity seen previously in ammonia molecular lines and infrared images, and proposed to be evacuated by the outflow associated with the massive VLA 3 source. The masers in the region of VLA 2 are located at projected distances of $\sim$ 0.7$''$ ($\sim$ 2300 AU) north from this source, with their kinematics suggesting that they are excited by a YSO other than VLA 2. This driving source has not yet been identified.

We present a formalism of the dynamics of internal shocks in relativistic jets where the source has a time-dependent injection velocity and mass-loss rate. The variation of the injection velocity produces a two-shock wave structure, the working surface, that moves along the jet. This new formalism takes into account the fact that momentum conservation is not valid for relativistic flows where the relativistic mass lost by radiation must be taken into account, in contrast to the classic regime. We find analytic solutions for the working surface velocity and radiated energy for the particular case of a step function variability of the injection parameters. We model two cases: a pulse of fast material and a pulse of slow material (with respect to the mean flow). Applying these models to gamma ray burst light curves, one can determine the ratio of the Lorentz factors gamma_2 / gamma_1 and the ratio of the mass-loss rates dot{m_2} / dot{m_1} of the upstream and downstream flows. As an example, we apply this model to the sources GRB 080413B and GRB 070318 and find the values of these ratios. Assuming a Lorentz factor gamma_1=100, we further estimate jet mass-loss rates between dot{m_1} ~ 10^{-5}-1 Msun.yr^{-1}. We also calculate the fraction of the injected mass lost by radiation. For GRB 070318 this fraction is ~7%. In contrast, for GRB 080413B this fraction is larger than 50%; in this case radiation losses clearly affect the dynamics of the internal shocks.

In this paper we analyze multi-epoch Very Long Baseline Interferometry (VLBI) water maser observations carried out with the Very Long Baseline Array (VLBA) toward the high-mass star-forming region AFGL 2591. We detected maser emission associated with the radio continuum sources VLA 2 and VLA 3. In addition, a water maser cluster, VLA 3-N, was detected ~ 0.5" north of VLA 3. We concentrate the discussion of this paper on the spatio-kinematical distribution of the water masers towards VLA 3-N. The water maser emission toward the region VLA 3-N shows two bow shock-like structures, Northern and Southern, separated from each other by ~ 100 mas (~ 330 AU). The spatial distribution and kinematics of the water masers in this cluster have persisted over a time span of seven years. The Northern bow shock has a somewhat irregular morphology, while the Southern one has a remarkably smooth morphology. We measured the proper motions of 33 water maser features, which have an average proper motion velocity of ~ 1.3 mas/yr (~ 20 km/s). The morphology and the proper motions of this cluster of water masers show systematic expanding motions that could imply one or two different centers of star formation activity. We made a detailed model for the Southern structure, proposing two different kinematic models to explain the 3-dimensional spatio-kinematical distribution of the water masers: (1) a static central source driving the two bow-shock structures; (2) two independent driving sources, one of them exciting the Northern bow-shock structure, and the other one, a young runaway star moving in the local molecular medium exciting and molding the remarkably smoother Southern bow-shock structure. Future observations will be necessary to discriminate between the two scenarios, in particular by identifying the still unseen driving source(s).

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

DG Tau B is a Class I young stellar source that drives the asymmetric HH 159 bipolar jet. At optical wavelengths it is obscured by circumstellar optically-thick material. Using VLA and JVLA observations, we determine for the first time the proper motions of this source and find them to be consistent, within error, with those of the nearby young star DG Tau. We also discuss an ejection event that is evident in the 1994 VLA data. As the optical and molecular outflows, this ejection traced in the radio continuum is markedly asymmetric and was detected only to the NW of the star. We propose that this knot, no longer detectable in the radio, could be observed in future optical images of DG Tau B. The positions of the VLA source and of a nearby infrared object are not coincident and we suggest that the VLA source traces the exciting object, while the infrared source could be a reflection lobe.

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.

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.

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 $\eta$ 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 $\eta$ 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 $\eta$ 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 $\eta$ Car.

We have studied the kinematics traced by the water masers located at the centre of the planetary nebula (PN) K3-35, using data from previous Very Large Array (VLA) observations. An analysis of the spatial distribution and line-of-sight velocities of the maser spots allows us to identify typical patterns of a rotating and expanding ring in the position-velocity diagrams, according to our kinematical model. We find that the distribution of the masers is compatible with tracing a circular ring with a ~0.021 arcsec (~100 AU) radius, observed with an inclination angle with respect to the line of sight of 55 degrees. We derive expansion and rotation velocities of 1.4 and 3.1 km/s, respectively. The orientation of the ring projected on the plane of the sky, at PA 158 degrees, is almost orthogonal to the direction of the innermost region of the jet observed in K3-35, suggesting the presence of a disc or torus that may be related to the collimation of the outflow.

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.

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 have studied the maser emission from a thin, planar, gaseous ring in Keplerian rotation around a central mass observed edge-on. The absorption coefficient within the ring is assumed to follow a power law dependence with the distance from the central mass as, k=k0r^{-q}. We have calculated position-velocity diagrams for the most intense maser features, for different values of the exponent q. We have found that, depending on the value of q, these diagrams can be qualitatively different. The most intense maser emission at a given velocity can either come mainly from regions close to the inner or outer edges of the amplifying ring or from the line perpendicular to the line of sight and passing through the central mass (as is commonly assumed). Particularly, when q>1 the position-velocity diagram is qualitatively similar to the one observed for the water maser emission in the nucleus of the galaxy NGC 4258. In the context of this simple model, we conclude that in this object the absorption coefficient depends on the radius of the amplifying ring as a decreasing function, in order to have significant emission coming from the inner edge of the ring.