We present a study of the evolution of the inner few astronomical units of protoplanetary disks around low-mass stars. We consider nearby stellar groups with ages spanning from 1 to 11 Myr, distributed into four age bins. Combining PANSTARSS photometry with spectral types, we derive the reddening consistently for each star, which we use (1) to measure the excess emission above the photosphere with a new indicator of IR excess and (2) to estimate the mass accretion rate ($\dot{M}$) from the equivalent width of the H$\alpha$ line. Using the observed decay of $\dot{M}$ as a constrain to fix the initial conditions and the viscosity parameter of viscous evolutionary models, we use approximate Bayesian modeling to infer the dust properties that produce the observed decrease of the IR excess with age, in the range between 4.5 and $24\,\mu$m. We calculate an extensive grid of irradiated disk models with a two-layered wall to emulate a curved dust inner edge and obtain the vertical structure consistent with the surface density predicted by viscous evolution. We find that the median dust depletion in the disk upper layers is $\epsilon \sim 3 \times 10^{-3}$ at 1.5 Myr, consistent with previous studies, and it decreases to $\epsilon \sim 3 \times 10^{-4}$ by 7.5 Myr. We include photoevaporation in a simple model of the disk evolution and find that a photoevaporative wind mass-loss rate of $\sim 1 -3 \times 10 ^{-9} \, M_{\odot}yr^{-1}$ agrees with the decrease of the disk fraction with age reasonably well. The models show the inward evolution of the H$_2$O and CO snowlines.

It is known that the millimeter dust thermal emission of protoplanetary disks is affected by scattering, such that for optically thick disks the emission decreases with respect to the pure absorption case and the spectral indices can reach values below 2. The latter can also be obtained with temperature gradients. Using simple analytical models of radiative transfer in thin slabs, we quantify the effect of scattering, vertical temperature gradients, and dust settling on the emission and spectral indices of geometrically thin face-on accretion disks around young stars. We find that in vertically isothermal disks with large albedo ($\omega_{\nu} \gtrsim 0.6$), the emergent intensity can increase at optical depths between $10^{-2}$ and $10^{-1}$. We show that dust settling has important effects on the spectral indices in the optically thick regime, since the disk emission mainly traces small dust grains in the upper layers of the disk. The $\lambda = 870 \ \mu$m emission of these small grains can hide large grains at the disk mid plane when the dust surface density is larger than $\sim$ 3.21 g cm$^{-2}$. Finally, because of the change of the shape of the spectral energy distribution, optically thick disks at 1.3 mm and grains with sizes between 300 $\mu$m $< a_{\rm max} <$ 1 mm have a 7 mm flux $\sim 60$\% higher than the extrapolation from higher millimeter frequencies, assumed when scattering is neglected. This effect could provide an explanation to the excess emission at $\lambda = 7$ mm reported in several disks.

The infrared source known as Orion n was detected in 1980 with observations made with the 3.8-m United Kingdom Infrared Telescope. About two decades later, sensitive observations made with the Very Large Array revealed the presence of a mJy double radio source apparently coincident in position with the infrared source n. The radio source was assumed to be the counterpart of the infrared source. However, over the years it has been concluded that the radio source shows large proper motions to the south while the infrared source n is stationary. Here we reanalyze the proper motions of the radio source adding both older and newer VLA observations than previously used. We confirm the proper motions of the radio source that at present no longer coincides positionally with the infrared source. The solution to this problem is, most probably, that the infrared source n and the radio source are not the same object: the infrared source is a stationary object in the region while the radio counterpart is moving as a result of the explosion that took place in this region some 500 years ago and that expelled large amounts of molecular gas as well as several compact sources. Considering the paper where it was first reported, we refer to this double radio source as Orion MR. In addition, we use these new observations to fully confirm the large proper motions of the sources IRc23 and Zapata 11. Together with sources BN, I, Orion MR, and x, there are at least six compact sources that recede from a point in common in Orion BN/KL. However, IRc23 is peculiar in that its ejection age appears to be only $\sim$300 years. The relatively large number of sources rules out as a possible mechanism the classic three-body scenario since then only two escaping bodies are expected: a tight binary plus the third star involved in the encounter.

We model the ALMA and VLA millimeter radial profiles of the disk around HL Tau to constrain the properties of the dust grains. We adopt the disk evolutionary models of Lynden-Bell \& Pringle and calculate their temperature and density structure and emission. These disks are heated by the internal viscosity and irradiated by the central star and a warm envelope. We consider a dust size distribution $n(a) da \propto a^{-3.5} da $, and vary the maximum grain size in the atmosphere and the midplane, $a_{\rm max}=100\ \mu$m, 1 mm, and 1cm. We also include dust settling and vary the dust-to-gas mass ratio from 1 to 9 times the ISM value. We find that the models that can fit the observed level of emission along the profiles at all wavelengths have an atmosphere with a maximum grain size $a_{\rm max} = 100 \ \mu$m, and a midplane with $a_{\rm max}=1$ cm. The disk substructure, with a deficit of emission in the gaps, can be due to dust properties in these regions that are different from those in the rings. We test an opacity effect (different $a_{\rm max}$) and a dust mass deficit (smaller dust-to-gas mass ratio) in the gaps. We find that the emission profiles are better reproduced by models with a dust deficit in the gaps, although a combined effect is also possible. These models have a global dust-to-gas mass ratio twice the ISM value, needed to reach the level of emission of the 7.8 mm VLA profile.

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.

We study dust concentration in axisymmetric gas rings in protoplanetary disks. Given the gas surface density, we derived an analytical total dust surface density by taking into account the differential concentration of all the grain sizes. This model allows us to predict the local dust-to-gas mass ratio and the slope of the particle size distribution, as a function of radius. We test this analytical model comparing it with a 3D magneto-hydrodynamical simulation of dust evolution in an accretion disk. The model is also applied to the disk around HD 169142. By fitting the disk continuum observations simultaneously at $\lambda = 0.87$, 1.3, 3.0 mm, we obtain a global dust-to-gas mass ratio $\epsilon_{\rm global} = 1.05 \times 10^{-2}$ and a viscosity coefficient $\alpha = 1.35 \times 10^{-2}$. This model can be easily implemented in numerical simulations of accretion disks.

We present Very Large Array observations at 33.0 GHz that detect emission coincident with $\epsilon$ Eridani to within $0\rlap.{"}07$ (0.2 AU at the distance of this star), with a positional accuracy of $0\rlap.{"}05$. This result strongly supports the suggestion of previous authors that the quiescent centimeter emission comes from the star and not from a proposed giant exoplanet with a semi-major axis of $\sim1\rlap.{"}0$ (3.4 AU). The centimeter emission is remarkably flat and is consistent with optically thin free-free emission. In particular, it can be modeled as a stellar wind with a mass loss rate of the order of $6.6 \times 10^{-11}~ M_\odot ~yr^{-1}$, which is 3,300 times the solar value, exceeding other estimates of this star's wind. However, interpretation of the emission in terms of other thermal mechanisms like coronal free-free and gyroresonance emission cannot be discarded.

In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $\alpha_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclination angle $i$ and observing frequency $\nu$ from the VLA Ka band ($\sim 37$ GHz) to ALMA Band 8 ($\sim 405$ GHz). We also consider the effects of decreasing the index of the (sub-)mm dust opacity power law $\beta$ from 1.7 to 1. We find that a distribution of disk masses in the range $M_{\rm disk} = 0.01-2~M_\odot$ is needed to reproduce the observed distribution of spectral indices, and that assuming a fixed $\beta =1.7$ gives better results than $\beta=1$. A wide distribution of disk masses is also needed to produce some cases with $\alpha <2$, as reported for some sources in the literature. Such extremely low spectral indices arise naturally when the selected observing frequencies sample the appropriate change in the temperature structure of the optically thick model disk. Our results show that protostellar disk masses could often be underestimated by $> \times10$, and are consistent with recent hydrodynamical simulations. Although we do not rule out the possibility of some grain growth occurring within the short protostellar timescales, we conclude that self-obscuration needs to be taken into account.

We calculate the emission of protoplanetary disks threaded by a poloidal magnetic field and irradiated by the central star. The radial structure of these disks was studied by Shu and collaborators and the vertical structure was studied by Lizano and collaborators. We consider disks around low mass protostars, T Tauri stars, and FU Ori stars with different mass-to-flux ratios $\lambda_{\rm sys}$. We calculate the spectral energy distribution and the antenna temperature profiles at 1 mm and 7 mm convolved with the ALMA and VLA beams. We find that disks with weaker magnetization (high values of $\lambda_{\rm sys}$) emit more than disks with stronger magnetization (low values of $\lambda_{\rm sys}$). This happens because the former are denser, hotter and have larger aspect ratios, receiving more irradiation from the central star. The level of magnetization also affects the optical depth at millimeter wavelengths, being larger for disks with high $\lambda_{\rm sys}$. In general, disks around low mass protostars and T Tauri stars are optically thin at 7 mm while disks around FU Ori are optically thick. A qualitative comparison of the emission of these magnetized disks, including heating by an external envelope, with the observed millimeter antenna temperature profiles of HL Tau indicates that large cm grains are required to increase the optical depth and reproduce the observed 7 mm emission at large radii.

We study the dust concentration and emission in protoplanetary disks vortices. We extend the Lyra-Lin solution for the dust concentration of a single grain size to a power-law distribution of grain sizes $n(a) \propto a^{-p}$. Assuming dust conservation in the disk, we find an analytic dust surface density as a function of the grain radius. We calculate the increase of the dust to gas mass ratio $\epsilon$ and the slope $p$ of the dust size distribution due to grain segregation within the vortex. We apply this model to a numerical simulation of a disk containing a persistent vortex. Due to the accumulation of large grains towards the vortex center, $\epsilon$ increases by a factor of 10 from the background disk value, and $p$ decreases from 3.5 to 3.0. We find the disk emission at millimeter wavelengths corresponding to synthetic observations with ALMA and VLA. The simulated maps at 7 mm and 1 cm show a strong azimuthal asymmetry. This happens because, at these wavelengths, the disk becomes optically thin while the vortex remains optically thick. The large vortex opacity is mainly due to an increase in the dust to gas mass ratio. In addition, the change in the slope of the dust size distribution increases the opacity by a factor of 2. We also show that the inclusion of the dust scattering opacity substantially changes the disks images.

Sensitive, imaging observations of the 1.1 mm dust continuum emission from a 1 deg^2 area collected with the AzTEC bolometer camera on the Large Millimeter Telescope are presented. A catalog of 1545 compact sources is constructed based on a Wiener-optimization filter. These sources are linked to larger clump structures identified in the Bolocam Galactic Plane Survey. Hydrogen column densities are calculated for all sources and mass and mean volume densities are derived for the subset of sources for which kinematic distances can be assigned. The AzTEC sources are localized, high density peaks within the massive clumps of molecular clouds and comprise 5-15% of the clump mass. We examine the role of the gravitational instability in generating these fragments by comparing the mass of embedded AzTEC sources to the Jeans' mass of the parent BGPS object. For sources with distances less than 6 kpc the fragment masses are comparable to the clump Jeans' mass, despite having isothermal Mach numbers between 1.6 and 7.2. AzTEC sources linked to ultra-compact HII regions have mass surface densities greater than the critical value implied by the mass-size relationship of infrared dark clouds with high mass star formation while AzTEC sources associated with Class II methanol masers have mass surface densities greater than 0.7 g cm^{-2} that approaches the proposed threshold required to form massive stars.

The possible occurrence of dead zones in protoplanetary disks subject to the magneto-rotational instability highlights the importance of disk ionization. We present a closed-form theory for the deep-down ionization by X-rays at depths below the disk surface dominated by far-ultraviolet radiation. Simple analytic solutions are given for the major ion classes, electrons, atomic ions, molecular ions and negatively charged grains. In addition to the formation of molecular ions by X-ray ionization of H2 and their destruction by dissociative recombination, several key processes that operate in this region are included, e.g., charge exchange of molecular ions and neutral atoms and destruction of ions by grains. Over much of the inner disk, the vertical decrease in ionization with depth into the disk is described by simple power laws, which can easily be included in more detailed modeling of magnetized disks. The new ionization theory is used to illustrate the non-ideal MHD effects of Ohmic, Hall and Ambipolar diffusion for a magnetic model of a T Tauri star disk using the appropriate Elsasser numbers.

We have extended the time baseline for observations of the proper motions of radio sources in the Orion BN/KL region from 14.7 to 22.5 years. We present improved determinations for the sources BN and I. In addition, we address the proper motions of the double radio source n, that have been questioned in the literature. We confirm that all three sources are moving away at transverse velocities of tens of km s$^{-1}$ from a region in-between them, where they were located about 500 years ago. Source n exhibits a new component that we interpret as due to a one-sided ejection of free-free emitting plasma that took place after 2006.36. We used the highly accurate relative proper motions between sources BN and I to determine that their closest separation took place in the year 1475$\pm$6, when they were within $\sim$100 AU or less from each other in the plane of the sky.

Photoevaporation is probably the main agent for gas dispersal during the last stages of protoplanetary disk evolution. However, the overall mass loss rate in the photoevaporative wind and its driving mechanism are still not well understood. Here we report multi-configuration Very Large Array observations at 0.7, 3, and 5 cm towards the transitional disk of GM Aur. Our radio continuum observations allow us to image and spatially resolve, for the first time, the three main components at work in this stage of the disk evolution: the disk of dust, the ionized radio jet perpendicular to it, and the photoevaporative wind arising from the disk. The mass loss rate inferred from the flux density of the radio jet is consistent with the ratio between ejection and accretion rates found in younger objects, suggesting that transitional disks can power collimated ejections of material apparently following the same physical mechanisms as much younger protostars. Our results indicate that extreme-UV (EUV) radiation is the main ionizing mechanism of the photoevaporative wind traced by the free-free emission. The required low EUV photon luminosity of $\sim6\times10^{40}$ s$^{-1}$ would produce a photoevaporation rate of only $\dot{M}_{\rm w,EUV}\simeq1.3\times10^{-10}~M_{\odot}$ yr$^{-1}$. Therefore, other mechanisms are required to disperse the disk in the timescale imposed by observations.

Epsilon Eridani is a nearby, young Sun-like star that hosts a ring of cool debris analogous to the solar system's Edgeworth-Kuiper belt. Early observations at (sub-)mm wavelengths gave tentative evidence of the presence of inhomogeneities in the ring, which have been ascribed to the effect of a putative low eccentricity planet, orbiting close to the ring. The existence of these structures have been recently challenged by high resolution interferometric millimeter observations. Here we present the deepest single-dish image of Epsilon Eridani at millimeter wavelengths, obtained with the Large Millimeter Telescope Alfonso Serrano (LMT). The main goal of these LMT observations is to confirm (or refute) the presence of non-axisymmetric structure in the disk. The dusty ring is detected for the first time along its full projected elliptical shape. The radial extent of the ring is not spatially resolved and shows no evidence, to within the uncertainties, of dust density enhancements. Additional features of the 1.1 mm map are: (i) the presence of significant flux in the gap between the ring and the star, probably providing the first exo-solar evidence of Poynting-Robertson drag, (ii) an unambiguous detection of emission at the stellar position with a flux significantly above that expected from Epsilon Eridani's photosphere, and (iii) the identification of numerous unresolved sources which could correspond to background dusty star-forming galaxies.

We model the vertical structure of magnetized accretion disks subject to viscous and resistive heating, and irradiation by the central star. We apply our formalism to the radial structure of magnetized accretion disks threaded by a poloidal magnetic field dragged during the process of star formation developed by Shu and coworkers. We consider disks around low mass protostars, T Tauri, and FU Orionis stars. We consider two levels of disk magnetization, $\lambda_{sys} = 4$ (strongly magnetized disks), and $\lambda_{sys} = 12$ (weakly magnetized disks). The rotation rates of strongly magnetized disks have large deviations from Keplerian rotation. In these models, resistive heating dominates the thermal structure for the FU Ori disk. The T Tauri disk is very thin and cold because it is strongly compressed by magnetic pressure; it may be too thin compared with observations. Instead, in the weakly magnetized disks, rotation velocities are close to Keplerian, and resistive heating is always less than 7\% of the viscous heating. In these models, the T Tauri disk has a larger aspect ratio, consistent with that inferred from observations. All the disks have spatially extended hot atmospheres where the irradiation flux is absorbed, although most of the mass ($\sim 90-95$ \%) is in the disk midplane. With the advent of ALMA one expects direct measurements of magnetic fields and their morphology at disk scales. It will then be possible to determine the mass-to-flux ratio of magnetized accretion disks around young stars, an essential parameter for their structure and evolution. Our models contribute to the understanding of the vertical structure and emission of these disks.

We discuss the effects of the magnetic field observed in molecular clouds on the process of star formation, concentrating on the phase of gravitational collapse of low-mass dense cores, cradles of sunlike stars. We summarize recent analytic work and numerical simulations showing that a substantial level of magnetic field diffusion at high densities has to occur in order to form rotationally supported disks. Furthermore, newly formed accretion disks are threaded by the magnetic field dragged from the parent core during the gravitational collapse. These disks are expected to rotate with a sub-Keplerian speed because they are partially supported by magnetic tension against the gravity of the central star. We discuss how sub-Keplerian rotation makes it difficult to eject disk winds and accelerates the process of planet migration. Moreover, magnetic fields modify the Toomre criterion for gravitational instability via two opposing effects: magnetic tension and pressure increase the disk local stability, but sub-Keplerian rotation makes the disk more unstable. In general, magnetized disks are more stable than their nonmagnetic counterparts; thus, they can be more massive and less prone to the formation of giant planets by gravitational instability.

We present $^{12}$CO(2-1) line and 1300 $\mu$m continuum observations made with the Submillimeter Array (SMA) of the young star DG Tau B. We find, in the continuum observations, emission arising from the circumstellar disk surrounding DG Tau B. The $^{12}$CO(2-1) line observations, on the other hand, revealed emission associated with the disk and the asymmetric outflow related with this source. Velocity asymmetries about the flow axis are found over the entire length of the flow. The amplitude of the velocity differences is of the order of 1 -- 2 km s$^{-1}$ over distances of about 300 -- 400 AU. We interpret them as a result of outflow rotation. The sense of the outflow and disk rotation is the same. Infalling gas from a rotating molecular core cannot explain the observed velocity gradient within the flow. Magneto-centrifugal disk winds or photoevaporated disk winds can produce the observed rotational speeds if they are ejected from a keplerian disk at radii of several tens of AU. Nevertheless, these slow winds ejected from large radii are not very massive, and cannot account for the observed linear momentum and angular momentum rates of the molecular flow. Thus, the observed flow is probably entrained material from the parent cloud. DG Tau B is a good laboratory to model in detail the entrainment process and see if it can account for the observed angular momentum.

Although gravitational collapse is supposed to play an essential role in the star formation process, infall motions have been always elusive to detect. So far, only a few observational signatures have been commonly used to claim for the presence of infall. Often these features consist in either "blue-asymmetries" or absorption at red-shifted velocities (e.g., inverse P-Cygni profiles). Both signatures are based only on the shape of the line profile and they do not guarantee by themselves the presence of dominant infall motions. More robust "mapping signatures" can be obtained from images that angularly resolve the infalling gas. Here we present VLA observations of the ammonia inversion transitions (2,2), (3,3), (4,4), (5,5), and (6,6) towards the hot molecular core (HMC) near G31.41+0.31 that show the signatures of protostellar infall theoretically predicted by Anglada et al. (1991). The intensity of the ammonia emission is compact and sharply increases towards the centre in the blue-shifted velocity channel maps, while it shows a more flattened distribution in the red-shifted velocity channels. Additionally, the emission becomes more compact with increasing (relative) velocity for both red and blue-shifted channels. We introduce a new infall signature, the "central blue spot", easily identifiable in the first-order moment maps. We show that rotation produces an additional, independent signature, making the distribution of the emission in the channel maps asymmetric with respect to the central position, but without masking the infall signatures. All these mapping signatures, which are identified here for the first time, are present in the observed ammonia transitions of G31 HMC.

We present new Karl G. Jansky Very Large Array observations of the compact (~ 0.05"), time-variable radio source projected near the center of the ultracompact HII region W3(OH). The analysis of our new data as well as of VLA archival observations confirms the variability of the source on timescales of years and for a given epoch indicates a spectral index of \alpha = 1.3 +- 0.3 (S_\nu \propto \nu^\alpha). This spectral index and the brightness temperature of the source (~6,500 K) suggest that we are most probably detecting partially optically thick free-free radiation. The radio source is probably associated with the ionizing star of W3(OH), but an interpretation in terms of an ionized stellar wind fails because the detected flux densities are orders of magnitude larger than expected. We discuss several scenarios and tentatively propose that the radio emission could arise in a static ionized atmosphere around a fossil photoevaporated disk.

We present sub-arcsecond observations toward the massive star forming region G75.78+0.34. We used the Very Large Array to study the centimeter continuum and H2O and CH3OH maser emission, and the Owens Valley Radio Observatory and Submillimeter Array to study the millimeter continuum and recombination lines (H40alpha and H30alpha). We found radio continuum emission at all wavelengths, coming from three components: (1) a cometary ultracompact (UC) HII region with an electron density 3.7x10^4 cm^{-3}, excited by a B0 type star, and with no associated dust emission; (2) an almost unresolved UCHII region (EAST), located 6 arcsec to the east of the cometary UCHII region, with an electron density 1.3x10^5 cm^{-3}, and associated with a compact dust clump detected at millimeter and mid-infrared wavelengths; and (3) a compact source (CORE), located 2 arcsec to the southwest of the cometary arc, with a flux density increasing with frequency, and embedded in a dust condensation of 30 Msun. The CORE source is resolved into two compact and unresolved sources which can be well-fit by two homogeneous hypercompact HII regions each one photo-ionized by a B0.5 ZAMS star, or by free-free radiation from shock-ionized gas resulting from the interaction of a jet/outflow system with the surrounding environment. The spatial distribution and kinematics of water masers close to the CORE-N and S sources, together with excess emission at 4.5 mum and the detected dust emission, suggest that the CORE source is a massive protostar driving a jet/outflow.

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.

Young stars typically have strong magnetic fields, so that the magnetospheres of newly formed close binaries can interact, dissipate energy, and produce synchrotron radiation. The V773 Tau A binary system, a pair of T Tauri stars with a 51 day orbit, displays such a signature, with peak emission taking place near periastron. This paper proposes that the observed emission arises from the change in energy stored in the composite magnetic field of the system. We model the fields using the leading order (dipole) components and show that this picture is consistent with current observations. In this model, the observed radiation accounts for a fraction of the available energy of interaction between the magnetic fields from the two stars. Assuming antisymmetry, we compute the interaction energy $E_{\rm int}$ as a function of the stellar radii, the stellar magnetic field strengths, the binary semi-major axis, and orbital eccentricity, all of which can be measured independently of the synchrotron radiation. The variability in time and energetics of the synchrotron radiation depend on the details of the annihilation of magnetic fields through reconnection events, which generate electric fields that accelerate charged particles, and how those charged particles, especially fast electrons, are removed from the interaction region. However, the major qualitative features are well described by the background changes in the global magnetic configuration driven by the orbital motion. The theory can be tested by observing a collection of pre-main-sequence binary systems.

This paper considers gravitational perturbations in geometrically thin disks with rotation curves dominated by a central object, but with substantial contributions from magnetic pressure and tension. The treatment is general, but the application is to the circumstellar disks that arise during the gravitational collapse phase of star formation. We find the dispersion relation for spiral density waves in these generalized disks and derive the stability criterion for axisymmetric $(m=0)$ disturbances (the analog of the Toomre parameter $Q_T$) for any radial distribution of the mass-to-flux ratio $\lambda$. The magnetic effects work in two opposing directions: on one hand, magnetic tension and pressure stabilize the disk against gravitational collapse and fragmentation; on the other hand, they also lower the rotation rate making the disk more unstable. For disks around young stars the first effect generally dominates, so that magnetic fields allow disks to be stable for higher surface densities and larger total masses. These results indicate that magnetic fields act to suppress the formation of giant planets through gravitational instability. Finally, even if gravitational instability can form a secondary body, it must lose an enormous amount of magnetic flux in order to become a planet; this latter requirement represents an additional constraint for planet formation via gravitational instability and places a lower limit on the electrical resistivity.

We present high angular and spectral resolution HI 21~cm line observations toward the cometary-shaped compact HII region G213.880-11.837 in the GGD~14 complex.The kinematics and morphology of the photodissociated region, traced by the HI line emission, reveal that the neutral gas is part of an expanding flow. The kinematics of the HI gas along the major axis of G213.880-11.837 shows that the emission is very extended toward the SE direction, reaching LSR radial velocities in the tail of about 14 km/s. The ambient LSR radial velocity of the molecular gas is 11.5 km/s, which suggests a champagne flow of the HI gas. This is the second (after G111.61+0.37) cometary HII/HI region known.

We present sensitive NIR (J, H and K) imaging observations toward four luminous massive star forming regions in the Norma Spiral Arm: G324.201+0.119, G328.307+0.432, G329.337+0.147 and G330.949-0.174. We identify three clusters of young stellar objects (YSO) based on surface density diagnostics. We also find that sources detected only in the H and K-bands and with colors corresponding to spectral types earlier than B2, are likely YSOs. We analyze the spatial distribution of stars of different masses and find signatures in two clusters of primordial mass segregation which can't be explained as due to incompleteness effects. We show that dynamic interactions of cluster members with the dense gas from the parent core can explain the observed mass segregation, indicating that the gas plays an important role in the dynamics of young clusters.

Magnetic fields are dragged in from the interstellar medium during the gravitational collapse that forms star/disk systems. Consideration of mean field magnetohydrodynamics (MHD) in these disks shows that magnetic effects produce subkeplerian rotation curves and truncate the inner disk. This letter explores the ramifications of these predicted disk properties for the migration of extrasolar planets. Subkeplerian flow in gaseous disks drives a new migration mechanism for embedded planets and modifies the gap opening processes for larger planets. This subkeplerian migration mechanism dominates over Type I migration for sufficiently small planets (m_P < 1 M_\earth) and/or close orbits (r < 1 AU). Although the inclusion of subkeplerian torques shortens the total migration time by only a moderate amount, the mass accreted by migrating planetary cores is significantly reduced. Truncation of the inner disk edge (for typical system parameters) naturally explains final planetary orbits with periods P=4 days. Planets with shorter periods P=2 days can be explained by migration during FU-Ori outbursts, when the mass accretion rate is high and the disk edge moves inward. Finally, the midplane density is greatly increased at the inner truncation point of the disk (the X-point); this enhancement, in conjunction with continuing flow of gas and solids through the region, supports the in situ formation of giant planets.

We present a model aimed to reproduce the observed spectral energy distribution (SED) as well as the ammonia line emission of the G31.41+0.31 hot core. The core is modeled as an infalling envelope onto a massive star that is undergoing an intense accretion phase. We assume an envelope with a density and velocity structure resulting from the dynamical collapse of a singular logatropic sphere. The stellar and envelope physical properties are determined by fitting the SED. From these physical conditions, the ammonia line emission is calculated and compared with subarcsecond resolution VLA data of the (4,4) transition. The only free parameter in this line fitting is the ammonia abundance. The observed properties of the NH3(4,4) lines and their spatial distribution can be well reproduced provided it is taken into account the steep increase of the abundance in the hotter (> 100 K), inner regions of the core produced by the sublimation of icy mantles where ammonia molecules are trapped. The model predictions for the (2,2), (4,4), and (5,5) transitions are also in reasonably agreement with the single-dish spectra available in the literature. The best fit is obtained for a model with a star of 25 Msun, a mass accretion rate of 0.003 Msun/yr, and a total luminosity of 200,000 Lsun. The gas-phase ammonia abundance ranges from 2 \times 10^{-8} in the outer region to 3 \times 10^{-6} in the inner region. To our knowledge, this is the first time that the dust and molecular line data of a hot molecular core, including subarcsecond resolution data that spatially resolve the structure of the core, have been simultaneously explained by a physically self-consistent model. This modeling shows that massive protostars are able to excite high excitation ammonia transitions up to the outer edge (30,000 AU) of the large scale envelope.

Strong magnetization makes the disks surrounding young stellar objects rotate at rates that are too sub-Keplerian to enable the thermal launching of disk winds from their surfaces unless the rate of gas diffusion across field lines is dynamically fast. This under-appreciated implication of disk magnetization poses a considerable challenge for disk-wind theory.

We present absolute astrometry of four radio sources in the Becklin-Neugebauer/Kleinman-Low (BN/KL) region, derived from archival data (taken in 1991, 1995, and 2000) as well as from new observations (taken in 2006). All data consist of 3.6 cm continuum emission and were taken with the Very Large Array in its highest angular resolution A configuration. We confirm the large proper motions of the BN object, the radio source I (GMR I) and the radio counterpart of the infrared source n (Orion-n), with values from 15 to 26 km/s. The three sources are receding from a point between them from where they seem to have been ejected about 500 years ago, probably via the disintegration of a multiple stellar system. We present simulations of very compact stellar groups that provide a plausible dynamical scenario for the observations. The radio source Orion-n appeared as a double in the first three epochs, but as single in 2006. We discuss this morphological change. The fourth source in the region, GMR D, shows no statistically significant proper motions. We also present new, accurate relative astrometry between BN and radio source I that restrict possible dynamical scenarios for the region. During the 2006 observations, the radio source GMR A, located about 1' to the NW of the BN/KL region, exhibited an increase in its flux density of a factor of ~3.5 over a timescale of one hour. This rapid variability at cm wavelengths is similar to that previously found during a flare at millimeter wavelengths that took place in 2003.

We consider the formation of OB associations from two perspectives: (a) the fractional gas consumption in star formation,epsilon, per dynamical time scale t_dyn in a galaxy, and (b) the origin of the so-called Kennicutt-Schmidt law that the rate of star formation per unit area is proportional to a power, alpha, of the surface density in HI and H_2 gas when certain thresholds are crossed. The empirical findings that epsilon is approximately 0.01 and alpha is approximately 1.4 or 1.5 have simple explanations if the rate of star formation is magnetically regulated. An empirical test of the ideas resides in an analysis of why giant OB associations are ``strung out like pearls along the arms" of spiral galaxies.

We consider the accretion process in a disk with magnetic fields that are dragged in from the interstellar medium by gravitational collapse. Two diffusive processes are at work in the system: (1) "viscous" torques exerted by turbulent and magnetic stresses, and (2) "resistive" redistribution of mass with respect to the magnetic flux arising from the imperfect conduction of current. In steady state, self-consistency between the two rates of drift requires that a relationship exists between the coefficients of turbulent viscosity and turbulent resistivity. Ignoring any interactions with a stellar magnetosphere, we solve the steady-state equations for a magnetized disk under the gravitational attraction of a mass point and threaded by an amount of magnetic flux consistent with calculations of magnetized gravitational collapse in star formation. Our model mean-field equations have an exact analytical solution that corresponds to magnetically diluted Keplerian rotation about the central mass point. The solution yields the strength of the magnetic field and the surface density as functions of radial position in the disk and their connection with the departure from pure Keplerian rotation in representative cases. We compare the predictions of the theory with the available observations concerning T Tauri stars, FU Orionis stars, and low- and high-mass protostars. Finally, we speculate on the physical causes for high and low states of the accretion disks that surround young stellar objects. One of the more important results of this study is the physical derivation of analytic expressions for the turbulent viscosity and turbulent resistivity.

We study the self-similar collapse of an isothermal magnetized rotating cloud in the ideal magnetohydrodynamic (MHD) regime. In the limit of small distance from the accreting protostar we find an analytic solution that corresponds to free-fall onto a central mass point. The density distribution is not spherically symmetric but depends on the mass loading of magnetic field lines, which can be obtained by matching our inner solution to an outer collapse solution previously computed by Allen, Shu & Li. The concentration of magnetic field trapped by the central mass point under field-freezing, independent on the details of the starting state, creates a split monopole configuration where the magnetic field strength increases as the inverse square of the distance from the center. Under such conditions, the inflow eventually becomes subalfvenic and the outward transfer of angular momentum by magnetic braking very efficient, thus preventing the formation of a centrifugally supported disk. Instead, the azimuthal velocity of the infalling gas decreases to zero at the center, and the gas spirals into the star. Therefore, the dissipation of dynamically important levels of magnetic field is a fundamental requisite for the formation of protoplanetary disks around young stars.

We formulate the problem of magnetic field dissipation during the accretion phase of low-mass star formation, and we carry out the first step of an iterative solution procedure by assuming that the gas is in free-fall along radial field lines. This so-called ``kinematic approximation'' ignores the back reaction of the Lorentz force on the accretion flow. In quasi steady-state, and assuming the resistivity coefficient to be spatially uniform, the problem is analytically soluble in terms of Legendre's polynomials and confluent hypergeometric functions. The dissipation of the magnetic field occurs inside a region of radius inversely proportional to the mass of the central star (the ``Ohm radius''), where the magnetic field becomes asymptotically straight and uniform. In our solution, the magnetic flux problem of star formation is avoided because the magnetic flux dragged in the accreting protostar is always zero. Our results imply that the effective resistivity of the infalling gas must be higher by several orders of magnitude than the microscopic electric resistivity, to avoid conflict with measurements of paleomagnetism in meteorites and with the observed luminosity of regions of low-mass star formation.

We review the phenomenon of ultra-compact H II regions (UCHIIs) as a key phase in the early lives of massive stars. This most visible manifestation of massive star formation begins when the Lyman continuum output from the massive young stellar object becomes sufficient to ionize the surroundings from which it was born. Knowledge of this environment is gained through an understanding of the morphologies of UCHII regions and we examine the latest developments in deep radio and mid-IR imaging. SPITZER data from the GLIMPSE survey are an important new resource in which PAH emission and the ionizing stars can be seen. We review the role played by strong stellar winds from the central stars in sweeping out central cavities and causing the limb-brightened appearance. A range of evidence from velocity structure, proper motions, the molecular environment and recent hydrodynamical modeling indicates that cometary UCHII regions require a combination of champagne flow and bow shock motion. Finally, we discuss the class of hyper-compact H II regions or broad recombination line objects. They are likely to mark the transition soon after the breakout of the Lyman continuum radiation from the young star. Models for these objects are presented, including photo-evaporating disks and ionized accretion flows that are gravitationally trapped. Evolutionary scenarios tracing young massive stars passage through these ionized phases are discussed.

We present absolute astrometry of 35 radio sources in the Orion Trapezium and Becklin-Neugebauer/Kleinman-Low regions, obtained from Very Large Array archival observations collected over a period of 15 years. By averaging the results for all the sources, we estimate the mean absolute proper motion of Orion to be --in Galactic coordinates-- $\mu_\ell \cos b$ = +2.1 $\pm$ 0.2 mas yr$^{-1}$; $\mu_b$ = $-$0.1 $\pm$ 0.2 mas yr$^{-1}$. These values agree remarkably well with those expected from the differential rotation of the Milky Way. Subtraction of this mean motion from the individual measurements allows us to register all proper motions to the rest frame of the Orion nebula, and identify radio sources with large residual velocities. In the KL region, we find three sources in this situation: the BN object, the radio source I, and the radio counterpart of the infrared source n. All three objects appear to be moving away from a common point where they must all have been located about 500 years ago. This suggests that all three sources were originally part of a multiple massive stellar system that recently disintegrated as a result of a close dynamical

We present absolute astrometry of the core of the Orion molecular cloud, made with Very Large Array archive data taken over the last two decades. Our analysis reveals that both the BN object and the radio source I have proper motions: the BN object has a proper motion of $12.6 \pm 0.6$ mas yr$^{-1}$ (corresponding to a velocity of $27 \pm 1$ km s$^{-1}$ at an adopted distance of 450 pc) to the northwest, while the radio source I has a proper motion of $5.6 \pm 0.7$ mas yr$^{-1}$ (corresponding to a velocity of $12 \pm 2$ km s$^{-1}$) to the southeast. The motion of the two sources is nearly antiparallel, diverging from a point in between them, where they were located about 500 years ago. These results suggest that the BN object and the radio source I were part of a multiple young stellar system that disintegrated in the recent past.

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

In order to compare the x-wind with observations, one needs to be able to calculate its thermal and ionization properties. We formulate the physical basis for the streamline-by-streamline integration of the ionization and heat equations of the steady x-wind. In addition to the well-known processes associated with the interaction of stellar and accretion-funnel hot-spot radiation with the wind, we include X-ray heating and ionization, mechanical heating, and a revised calculation of ambipolar diffusion heating. The mechanical heating arises from fluctuations produced by star-disk interactions of the time dependent x-wind that are carried by the wind to large distances where they are dissipated in shocks, MHD waves, and turbulent cascades. We model the time-averaged heating by the scale-free volumetric heating rate, $\Gamma_{\rm mech} = \alpha \rho v^3 s^{-1}$, where $\rho$ and $v$ are the local mass density and wind speed, respectively, $s$ is the distance from the origin, and $\alpha$ is a phenomenological constant. When we consider a partially-revealed but active young stellar object, we find that choosing $\alpha \sim 10^{-3}$ in our numerical calculations produces temperatures and electron fractions that are high enough for the x-wind jet to radiate in the optical forbidden lines at the level and on the spatial scales that are observed. We also discuss a variety of applications of our thermal-chemical calculations that can lead to further observational checks of x-wind theory.

We review the difficulties of the classical fission and fragmentation hypotheses for the formation of binary and multiple stars. A crucial missing ingredient in previous theoretical studies is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore formulate and solve the problem of the equilibria of isopedically magnetized, singular isothermal disks, without the assumption of axial symmetry. Considerable analytical progress can be made if we restrict our attention to models that are scale-free, i.e., that have surface densities that vary inversely with distance from the rotation axis of the system. In agreement with earlier analysis by Syer and Tremaine, we find that lopsided (M=1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into M=2, 3, 4, ... separate bodies. On the basis of our experience in this paper, and the preceding Paper I, we advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

A crucial missing ingredient in previous theoretical studies of fragmentation is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore consider the equilibrium configuration of isopedically magnetized, scale-free, singular isothermal disks, without the assumption of axial symmetry. We find that lopsided (M=1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into separate bodies. We advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

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.

Cloud environment is thought to play a critical role in determining the mechanism of formation of massive stars. In this contribution we review the physical characteristics of the environment around recently formed massive stars. Particular emphasis is given to recent high angular resolution observations which have improved our knowledge of the physical conditions and kinematics of compact regions of ionized gas and of dense and hot molecular cores associated with luminous O and B stars. We will show that this large body of data, gathered during the last decade, has allowed significant progress in the understanding of the physical processes that take place during the formation and early evolution of massive stars.

We investigate the equilibrium properties of self-gravitating magnetized clouds with polytropic equations of state with negative index n. In particular, we consider scale-free isopedic configurations that have constant dimensionless spherical mass-to-flux ratio lambda_r and that may constitute ``pivotal'' states for subsequent dynamical collapse to form groups or clusters of stars. For given Gamma = 1 + 1/n, equilibria with smaller values of lambda_r are more flattened, ranging from spherical configurations with lambda_r = infinity to completely flattened states for lambda_r=1. For a given amount of support provided by the magnetic field as measured by the dimensionless parameter H_0, equilibria with smaller values of Gamma are more flattened. However, logatropic (defined by Gamma = 0) disks do not exist. The only possible scale-free isopedic equilibria with logatropic equation of state are spherical uniformly magnetized clouds.

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