We present Very Long Baseline Interferometry (VLBI) observations of 22 GHz H$_2$O masers in the high-mass star-forming region of \objectname{W75N}, carried out with VLBI Exploration of Radio Astrometry (VERA) for three-epochs in 2007 with an angular resolution of $\sim$ 1 mas. We detected H$_2$O maser emission toward the radio jet in VLA 1 and the expanding shell-like structure in VLA 2. We have made elliptical fits to the VLA 2 H$_2$O maser shell-like structure observed in the different epochs (1999, 2005, and 2007), and found that the shell is still expanding eight years after its discovery. From the difference in the size of the semi-major axes of the fitted ellipses in the epochs 1999 ($\simeq$ 71$\pm$1 mas), 2005 ($\simeq$ 97$\pm$3 mas), and 2007 ($\simeq$ 111$\pm$1 mas), we estimate an average expanding velocity of $\sim$ 5 mas yr$^{-1}$, similar to the proper motions measured in the individual H$_2$O maser features. A kinematic age of $\sim$ 20 yr is derived for this structure. In addition, our VERA observations indicate an increase in the ellipticity of the expanding shell around VLA 2 from epochs 1999 to 2007. In fact, the elliptical fit of the VERA data shows a ratio between the minor and major axes of $\sim$ 0.6, in contrast with a most circular shape for the shell detected in 1999 and 2005 (b/a $\sim$ 0.9). This suggests that we are probably observing the formation of a jet-driven H$_2$O maser structure in VLA2, evolving from a non-collimated pulsed-outflow event during the first stages of evolution of a massive young stellar object (YSO). This may support predictions made earlier by other authors on this issue, consistent with recent magnetohydrodynamical simulations. We discuss possible implications of our results in the study of the first stages of evolution of massive YSOs.

The classic HII region M17 is one of the best studied across the electromagnetic spectrum. We present sensitive, high angular resolution observations made with the Jansky Very Large Array (JVLA) at 4.96 and 8.46 GHz that reveal the presence of 38 compact radio sources, in addition to the well known hypercompact cometary HII region M17 UC1. For this last source we find that its spectral index of value $\sim$1 is due to a gradient in opacity across its face. Of the 38 compact radio sources detected, 19 have stellar counterparts detected in the infrared, optical, or X-rays. Finally, we discuss the nature of the radio emission from the massive binary system CEN 1a and 1b, concluding that both are most probably non-thermal emitters, although the first is strongly time variable and the second is steady.

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.

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.

Aims. The interaction of two isotropic stellar winds is studied in order to calculate the free-free emission from the wind collision region. The effects of the binary separation and the wind momentum ratio on the emission from the wind-wind interaction region are investigated. Methods. We developed a semi-analytical model for calculating the thermal emission from colliding stellar winds. Assuming radiative shocks for the compressed layer, which are expected in close binaries, we obtained the emission measure of the thin shell. Then, we computed the total optical depth along each line of sight to obtain the emission from the whole configuration. Results. Here, we present predictions of the free-free emission at radio frequencies from analytic, radiative shock models in colliding wind binaries. It is shown that the emission from the wind collision region mainly arises from the optically thick region of the compressed layer and scales as ~ D^{4/5}, where D is the binary separation. The predicted flux density from the wind collision region becomes more important as the frequency increases, showing higher spectral indices than the expected 0.6 value from the unshocked winds. We also investigate the emission from short-period WR+O systems calculated with our analytic formulation. In particular, we apply the model to the binary systems WR 98 and WR 113 and compare our results with the observations. Our theoretical results are in good agreement with the observed thermal spectra from these sources.

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 present five epochs of VLBI water maser observations around the massive protostar Cepheus A HW2 with 0.4 mas (0.3 AU) resolution. The main goal of these observations was to follow the evolution of the remarkable water maser linear/arcuate structures found in earlier VLBI observations. Comparing the data of our new epochs of observation with those observed five years before, we find that at "large" scales of > 1" (700 AU) the main regions of maser emission persist, implying that both the surrounding medium and the exciting sources of the masers have been relatively stable during that time span. However, at smaller scales of < 0.1" (70 AU) we see large changes in the maser structures, particularly in the expanding arcuate structures R4 and R5. R4 traces a nearly elliptical patchy ring of ~ 70 mas size (50 AU) with expanding motions of ~ 5 mas/yr (15 km/s). This structure is probably driven by the wind of a still unidentified YSO located at the centre of the ring (~ 0.18" south of HW2). On the other hand, the R5 expanding bubble structure (driven by the wind of a previously identified YSO located ~ 0.6" south of HW2) is currently dissipating in the circumstellar medium and losing its previous degree of symmetry, indicating a very short-lived event. In addition, our results reveal, at scales of ~ 1" (700 AU), the simultaneous presence of a relatively slow (~ 10-70 km/s) wide-angle outflow (opening angle of ~ 102 deg, traced by the masers, and the fast (~ 500~km/s) highly collimated radio jet associated with HW2 (opening angle of ~ 18 deg, previously observed with the VLA. This simultaneous presence of a wide-angle outflow and a highly collimated jet associated with a massive protostar is similar to what is found in some low-mass YSOs. The implications of these results in the study of the formation of high-mass stars are discussed.

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 $\mu_\alpha \cos \delta = +10.9 \pm 3.0$ mas yr$^{-1}$; $\mu_\delta = +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.

Context. Mathematical optimization can be used as a computational tool to obtain the optimal solution to a given problem in a systematic and efficient way. For example, in twice-differentiable functions and problems with no constraints, the optimization consists of finding the points where the gradient of the objective function is zero and using the Hessian matrix to classify the type of each point. Sometimes, however it is impossible to compute these derivatives and other type of techniques must be employed such as the steepest descent/ascent method and more sophisticated methods such as those based on the evolutionary algorithms. Aims. We present a simple algorithm based on the idea of genetic algorithms (GA) for optimization. We refer to this algorithm as AGA (Asexual Genetic Algorithm) and apply it to two kinds of problems: the maximization of a function where classical methods fail and model fitting in astronomy. For the latter case, we minimize the chi-square function to estimate the parameters in two examples: the orbits of exoplanets by taking a set of radial velocity data, and the spectral energy distribution (SED) observed towards a YSO (Young Stellar Object). Methods. The algorithm AGA may also be called genetic, although it differs from standard genetic algorithms in two main aspects: a) the initial population is not encoded, and b) the new generations are constructed by asexual reproduction. Results. Applying our algorithm in optimizing some complicated functions, we find the global maxima within a few iterations. For model fitting to the orbits of exoplanets and the SED of a YSO, we estimate the parameters and their associated errors.

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.

Neutral clumps immersed in HII regions are frequently found in star formation regions. We investigate here the formation of tail of neutral gas, which are not reached by the direct ionizing flux coming from two massive stars, using both an analytical approximation, that allows us to estimate the shadow geometry behind the clumps for different initial geometric configurations, and three-dimensional numerical simulations. We found a good agreement between both approaches to this theoretically interesting problem. A particularly important application could be the proplyds that are found in the Trapezium cluster in Orion, which are being photoevaporated primordially by the O stars theta-1 Ori C and theta-2 Ori A.

Molecular clouds in the interstellar medium suffer gravitational instabilities that lead to the formation of one or multiple stars. A recently formed star inside a cold cloud communicates its gravitational force to the surrounding environment and soon an accretion flow falling into the star develops. After their formation, all stars soon eject a wind of gas that interacts with the external accretion flow. This interaction produces a shock wave that evolves with time. The work presented in this article formulates a simple prescription for the evolution of this interaction. With the aid of this model we construct a few radio continuum maps of the source.

We report centimeter continuum observations with the VLA and the VLBA as well as mid-infrared observations with COMICS/SUBARU toward the components of the YLW 15 very young binary system, VLA 1 and VLA 2. The centimeter emission of the two components traces partially thick free-free emission, likely due to collimated, ionized winds. VLA 1 is an embedded protostar, undetected in the near-IR, and possibly in the Class 0 to Class I transition and powering a Herbig-Haro outflow. Its mid-IR emission appears slightly resolved with a diameter of ~16 AU, possibly tracing circumstellar material from both the envelope and the disk. VLA 2 is a typical Class I object, unresolved in the mid-IR, and is the responsible of the strong X-ray emission associated with YLW 15. The expected centimeter ''peri-stellar'' emission associated with the X-ray emission is not detected with the VLBA at 6 cm likely due to the high optical depth of the free-free emission. Strikingly, the near to mid-IR properties of YLW 15 suggest that VLA 1 is a more embedded YSO, or alternatively, less luminous than VLA 2, whereas orbital proper motions of this binary system indicate that VLA 1 is more massive than VLA 2. This result is apparently against the expected evolutionary scenario, where one expects that the more massive YSO in a binary system is the more evolved and luminous YSO. Finally, the nearby source YLW 16A is detected with the VLA, its position coincides with reported near-IR and X-ray sources.

We present an analytical model for the cross section of a turbulent, radiative jet or wake. This model is appropriate for modeling HH jets, or "wakes" left behind by "astrophysical bullets". Even though the model is very simple, it has the benign property of only having four free parameters (the outer radius of the beam, the axial velocity, the velocity at the edge of the beam, and the turbulent velocity width), which can be derived by fitting the radial velocity and line width cross sections of an observed outflow. We illustrate how to do such fits using previously published spectroscopic data of the HH 110 jet.

We have obtained a Halpha position-velocity cube from Fabry-Perot interferometric observations of the HH 110 flow. We analyze the results in terms of anisotropic wavelet transforms, from which we derive the spatial distribution of the knots as well as their characteristic sizes (along and across the outflow axis). We then study the spatial behaviour of the line width and the central radial velocity. The results are interpreted in terms of a simple ``mean flow+turbulent eddy'' jet/wake model. We find that most of the observed kinematics appear to be a direct result of the mean flow, on which are superposed low amplitude (35 km/s) turbulent velocities.

(Abridged) We consider the expansion of an initially self-gravitating, static, singular, isothermal cloud core. For t>0, the gas is ionized and heated to a higher uniform temperature by the formation of a luminous, but massless, star in its center. The approximation that the mass and gravity of the central star is negligible for the subsequent motion of the HII region holds for distances much greater than about 100 AU and for the massive cloud cores that give rise to high-mass stars. If the initial ionization and heating is approximated to occur instantaneously at t=0, then the subsequent flow (for r >> 100 AU) caused by the resulting imbalance between self-gravity and thermal pressure is self-similar. Because of the steep density profile, pressure gradients produce a shock front that travels into the cloud, accelerating the gas to supersonic velocities in what has been called the ``champagne phase.'' We then study the self-similar solutions of the expansion of HII regions embedded in molecular clouds characterized by more general power-law density distributions with exponent 3/2 < n < 3. In these cases, the shock velocity is an increasing function of the exponent n, and diverges as n tends to 3. We show that this happens because the model includes an origin, where the pressure driving the shock diverges because the enclosed heated mass is infinite. Our results imply that the continued photoevaporation of massive reservoirs of neutral gas (e.g., surrounding disks and/or globules) nearby to the embedded ionizing source is required in order to maintain over a significant timescale the emission measure observed in champagne flows.

We present multi-frequency VLA continuum observations towards 8 star forming regions with molecular and optical outflows: L1489, HH 68-69, HH 94-95, NGC 2264D, L1681B, L778, MWC 1080 and V645 Cyg. We detect three thermal radio jets, L1489, YLW 16A in L1681B and NGC 2264D VLA 7, associated with molecular and/or HH outflows. The L1489 and NGC 2264D VLA 7 thermal radio jets appear elongated in the direction of the larger scale outflow. We report the first tentative detection of a non-thermal radio jet, L778 VLA 5, associated with a low mass Class I protostar and powering a molecular outflow. For HH68-69, HH 94-95 and the molecular outflow in NGC 2264D we could not identify a candidate of the exciting source of these outflows. The radio emission associated with V645 Cyg is quite extended, ~0.1 pc, and time variable. We detect three radio sources in the MWC 1080 that could be associated with YSOs.

We explain the Ha emission of the cometary knots in the Helix Nebula (NGC 7293) with an analytical model that describes the emission of the head of the globules as a photoevaporated flow produced by the incident ionizing radiation of the central star.We compare these models with the Ha emission obtained from the HST archival images of the Helix Nebula. From a comparison of the Ha emission with the predictions of the analytical model we obtain a rate of ionizing photons from the central star of about 5e45 s^-1, which is consistent with estimates based on the total Hb flux of the nebula. We also model the tails of the cometary knots as a photoevaporated wind from a neutral shadow region produced by the diffuse ionizing photon field of the nebula. A comparison with the HST images allows us to obtain a direct determination of the value of the diffuse ionizing flux. We compare the ratio of diffuse to direct stellar flux as a function of radius inside an HII region with those obtained from the observational data through the analytical tail and head wind model. The agreement of this model with the values determined from the observations of the knots is excellent.

We construct detailed vertical structure models of irradiated accretion disks around T Tauri stars with interstellar medium dust uniformly mixed with gas. The dependence of the structure and emission properties on mass accretion rate, viscosity parameter, and disk radius is explored using these models. The theoretical spectral energy distributions (SEDs) and images for all inclinations are compared with observations of the entire population of Classical T Tauri stars (CTTS) and Class I objects in Taurus. In particular, we find that the median near-infrared fluxes can be explained within the errors with the most recent values for the median accretion rates for CTTS. We further show that the majority of the Class I sources in Taurus cannot be Class II sources viewed edge-on because they are too luminous and their colors would be consistent with disks seen only in a narrow range of inclinations. Our models appear to be too geometrically thick at large radii, as suggested by: (a) larger far-infrared disk emission than in the typical SEDs of T Tauri stars; (b) wider dark dust lanes in the model images than in the images of HH30 and HK Tau/c; and (c) larger predicted number of stars extincted by edge-on disks than consistent with current surveys. The large thickness of the model is a consequence of the assumption that dust and gas are well-mixed, suggesting that some degree of dust settling may be required to explain the observations.

We discuss the properties of an accretion disk around a star with parameters typical of classical T Tauri stars (CTTS), and with the average accretion rate for these disks. The disk is assumed steady and geometrically thin. The turbulent viscosity coefficient is expressed using the alpha prescription and the main heating mechanisms considered are viscous dissipation and irradiation by the central star. The energy is transported by radiation, turbulent conduction and convection. We find that irradiation from the central star is the main heating agent of the disk, except in the innermost regions, R less than 2 AU. The irradiation increases the temperature of the outer disk relative to the purely viscous case. As a consequence, the outer disk (R larger than 5 AU) becomes less dense, optically thin and almost vertically isothermal, with a temperature distribution T proportional to R^{-1/2}. The decrease in surface density at the outer disk, decreases the disk mass by a factor of 4 respect to a purely viscous case. In addition, irradiation tends to make the outer disk regions stable against gravitational instabilities.

We study the evolution of dense neutral clumps located in the outer parts of planetary nebulae. These clumps will be photo-ionized by the ionizing radiation from the central star and change their structure in the process. The main effect of the ionization process is the setting up of a photo-evaporation flow and a shock running through the clump. Once this shock has moved through the entire clump it starts to accelerate because of the `rocket effect'. This continues until the entire clump has been photo-ionized. We present an analytic model for the shock and accelerating phases and also the results of numerical simulations which include detailed microphysics. We find a good match between the analytic description and the numerical results and use the numerical results to produce some of the clump's observational characteristics at different phases of its evolution. We compare the results with the properties of the fast moving low ionization knots (ansae or FLIERs) seen in a number of planetary nebulae. We find that the models match many of the kinematic and emission properties of FLIERs.

We present numerical simulations of an isothermal turbulent gas undergoing gravitational collapse, aimed at testing for ``logatropic'' behavior of the form $P_t \sim \log \rho$, where $P_t$ is the ``turbulent pressure'' and $\rho$ is the density. To this end, we monitor the evolution of the turbulent velocity dispersion $\sigma$ as the density increases during the collapse. A logatropic behavior would require that $\sigma \propto \rho^{-1/2}$, a result which, however, is not verified in the simulations. Instead, the velocity dispersion increases with density, implying a polytropic behavior of $P_t$. This behavior is found both in purely hydrodynamic as well as hydromagnetic runs. For purely hydrodynamic and rapidly-collapsing magnetic cases, the velocity dispersion increases roughly as $\sigma \propto \rho^{1/2}$, implying $P_t\sim \rho^2$, where $P_t$ is the turbulent pressure. For slowly-collapsing magnetic cases the behavior is close to $\sigma \propto \rho^{1/4}$, which implies $P_t \sim \rho^{3/2}$. We thus suggest that the logatropic ``equation of state'' may represent only the statistically most probable state of an ensemble of clouds in equilibrium between self-gravity and kinetic support, but does not adequately represent the behavior of the ``turbulent pressure'' within a cloud undergoing a dynamic compression due to gravitational collapse. Finally, we discuss the importance of the underlying physical model for the clouds (in equilibrium vs. dynamic) on the results obtained.