A theory of potential vorticity (PV) mixing in a disordered (tangled) magnetic field is presented. The analysis is in the context of $\beta$-plane MHD, with a special focus on the physics of momentum transport in the stably stratified, quasi-2D solar tachocline. A physical picture of mean PV evolution by vorticity advection and tilting of magnetic fields is proposed. In the case of weak-field perturbations, quasi-linear theory predicts that the Reynolds and magnetic stresses balance as turbulence Alfv\'enizes for a larger mean magnetic field. Jet formation is explored quantitatively in the mean field-resistivity parameter space. However, since even a modest mean magnetic field leads to large magnetic perturbations for large magnetic Reynolds number, the physically relevant case is that of a strong but disordered field. We show that numerical calculations indicate that the Reynolds stress is modified well before Alfv\'enization -- i.e. before fluid and magnetic energies balance. To understand these trends, a double-average model of PV mixing in a stochastic magnetic field is developed. Calculations indicate that mean-square fields strongly modify Reynolds stress phase coherence and also induce a magnetic drag on zonal flows. The physics of transport reduction by tangled fields is elucidated and linked to the related quench of turbulent resistivity. We propose a physical picture of the system as a resisto-elastic medium threaded by a tangled magnetic network. Applications of the theory to momentum transport in the tachocline and other systems are discussed in detail.

The nearby ultra-luminous infrared galaxy (ULIRG) Arp 220 is an excellent laboratory for studies of extreme astrophysical environments. For 20 years, Very Long Baseline Interferometry (VLBI) has been used to monitor a population of compact sources thought to be supernovae (SNe), supernova remnants (SNRs) and possibly active galactic nuclei (AGNs). Using new and archival VLBI data spanning 20 years, we obtain 23 high-resolution radio images of Arp 220 at wavelengths from 18 cm to 2 cm. From model-fitting to the images we obtain estimates of flux densities and sizes of all detected sources. We detect radio continuum emission from 97 compact sources and present flux densities and sizes for all analysed observation epochs. We find evidence for a LD-relation within Arp 220, with larger sources being less luminous. We find a compact source LF $n(L)\propto L^\beta$ with $\beta=-2.19\pm0.15$, similar to SNRs in normal galaxies. Based on simulations we argue that there are many relatively large and weak sources below our detection threshold. The observations can be explained by a mixed population of SNe and SNRs, where the former expand in a dense circumstellar medium (CSM) and the latter interact with the surrounding interstellar medium (ISM). Nine sources are likely luminous, type IIn SNe. This number of luminous SNe correspond to few percent of the total number of SNe in Arp 220 which is consistent with a total SN-rate of 4 yr$^{-1}$ as inferred from the total radio emission given a normal stellar initial mass function (IMF). Based on the fitted luminosity function, we argue that emission from all compact sources, also below our detection threshold, make up at most 20\% of the total radio emission at GHz frequencies.

In this paper we examine the role of weak magnetic fields in breaking Kelvin's circulation theorem and in vortex breakup in two-dimensional magnetohydrodynamics for the physically important case of a low magnetic Prandtl number (low $Pm$) fluid. We consider three canonical inviscid solutions for the purely hydrodynamical problem, namely a Gaussian vortex, a circular vortex patch and an elliptical vortex patch. We examine how magnetic fields lead to an initial loss of circulation $\Gamma$ and attempt to derive scaling laws for the loss of circulation as a function of field strength and diffusion as measured by two non-dimensional parameters. We show that for all cases the loss of circulation depends on the integrated effects of the Lorentz force, with the patch cases leading to significantly greater circulation loss. For the case of the elliptical vortex the loss of circulation depends on the total area swept out by the rotating vortex and so this leads to more efficient circulation loss than for a circular vortex.

We report the detection of an ultra-bright fast radio burst (FRB) from a modest, 3.4-day pilot survey with the Australian Square Kilometre Array Pathfinder. The survey was conducted in a wide-field fly's-eye configuration using the phased-array-feed technology deployed on the array to instantaneously observe an effective area of $160$ deg$^2$, and achieve an exposure totaling $13200$ deg$^2$ hr. We constrain the position of FRB 170107 to a region $8'\times8'$ in size (90% containment) and its fluence to be $58\pm6$ Jy ms. The spectrum of the burst shows a sharp cutoff above $1400$ MHz, which could be either due to scintillation or an intrinsic feature of the burst. This confirms the existence of an ultra-bright ($>20$ Jy ms) population of FRBs.

Continuous plasma coherent emission is maintained by repetitive Langmuir collapse driven by the nonlinear evolution of a strong electron two-stream instability. The Langmuir waves are modulated by solitary waves in the linear stage, and by electrostatic whistler waves in the nonlinear stage. Modulational instability leads to Langmuir collapse and electron heating that fills in cavitons. The high pressure is released via excitation of a short wavelength ion acoustic mode that is damped by electrons and that re-excites small-scale Langmuir waves---this process closes a feedback loop that maintains the continuous coherent emission.

The unexpected energy spectrum of the positron/electron ratio is interpreted astrophysically, with a possible exception of the 100-300 GeV range. The data indicate that this ratio, after a decline between $0.5-8$ GeV, rises steadily with a trend towards saturation at 200-400GeV. These observations (except for the trend) appear to be in conflict with the diffusive shock acceleration (DSA) mechanism, operating in a \emph{single} supernova remnant (SNR) shock. We argue that $e^{+}/e^{-}$ ratio can still be explained by the DSA if positrons are accelerated in a \emph{subset} of SNR shocks which: (i) propagate in clumpy gas media, and (ii) are modified by accelerated CR \emph{protons}. The protons penetrate into the dense gas clumps upstream to produce positrons and, \emph{charge the clumps positively}. The induced electric field expels positrons into the upstream plasma where they are shock-accelerated. Since the shock is modified, these positrons develop a harder spectrum than that of the CR electrons accelerated in other SNRs. Mixing these populations explains the increase in the $e^{+}/e^{-}$ ratio at $E>8$ GeV. It decreases at $E<8$ GeV because of a subshock weakening which also results from the shock modification. Contrary to the expelled positrons, most of the antiprotons, electrons, and heavier nuclei, are left unaccelerated inside the clumps. Scenarios for the 100-300 GeV AMS-02 fraction exceeding the model prediction, including, but not limited to, possible dark matter contribution, are also discussed.

We use observations from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to search for transient radio sources in the field around the intermittent pulsar PSR J1107-5907. The pulsar is thought to switch between an "off" state in which no emission is detectable, a weak state and a strong state. We ran three independent transient detection pipelines on two-minute snapshot images from a 13 hour BETA observation in order to 1) study the emission from the pulsar, 2) search for other transient emission from elsewhere in the image and 3) to compare the results from the different transient detection pipelines. The pulsar was easily detected as a transient source and, over the course of the observations, it switched into the strong state three times giving a typical timescale between the strong emission states of 3.7 hours. After the first switch it remained in the strong state for almost 40 minutes. The other strong states lasted less than 4 minutes. The second state change was confirmed using observations with the Parkes radio telescope. No other transient events were found and we place constraints on the surface density of such events on these timescales. The high sensitivity Parkes observations enabled us to detect individual bright pulses during the weak state and to study the strong state over a wide observing band. We conclude by showing that future transient surveys with ASKAP will have the potential to probe the intermittent pulsar population.

We report the discovery of a new 21-cm HI absorption system using commissioning data from the Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder (ASKAP). Using the 711.5 - 1015.5 MHz band of ASKAP we were able to conduct a blind search for the 21-cm line in a continuous redshift range between $z = 0.4$ and 1.0, which has, until now, remained largely unexplored. The absorption line is detected at $z = 0.44$ towards the GHz-peaked spectrum radio source PKS B1740$-$517 and demonstrates ASKAP's excellent capability for performing a future wide-field survey for HI absorption at these redshifts. Optical spectroscopy and imaging using the Gemini-South telescope indicates that the HI gas is intrinsic to the host galaxy of the radio source. The narrow OIII emission lines show clear double-peaked structure, indicating either large-scale outflow or rotation of the ionized gas. Archival data from the \emph{XMM-Newton} satellite exhibit an absorbed X-ray spectrum that is consistent with a high column density obscuring medium around the active galactic nucleus. The HI absorption profile is complex, with four distinct components ranging in width from 5 to 300 km s$^{-1}$ and fractional depths from 0.2 to 20 per cent. In addition to systemic HI gas, in a circumnuclear disc or ring structure aligned with the radio jet, we find evidence for a possible broad outflow of neutral gas moving at a radial velocity of $v \sim 300$ km s$^{-1}$. We infer that the expanding young radio source ($t_{\rm age} \approx 2500$ yr) is cocooned within a dense medium and may be driving circumnuclear neutral gas in an outflow of $\sim$ 1 $\mathrm{M}_{\odot}$ yr$^{-1}$.

We present HI imaging of the galaxy group IC 1459 carried out with six antennas of the Australian SKA Pathfinder equipped with phased-array feeds. We detect and resolve HI in eleven galaxies down to a column density of $\sim10^{20}$ cm$^{-2}$ inside a ~6 deg$^2$ field and with a resolution of ~1 arcmin on the sky and ~8 km/s in velocity. We present HI images, velocity fields and integrated spectra of all detections, and highlight the discovery of three HI clouds -- two in the proximity of the galaxy IC 5270 and one close to NGC 7418. Each cloud has an HI mass of $10^9$ M$_\odot$ and accounts for ~15% of the HI associated with its host galaxy. Available images at ultraviolet, optical and infrared wavelengths do not reveal any clear stellar counterpart of any of the clouds, suggesting that they are not gas-rich dwarf neighbours of IC 5270 and NGC 7418. Using Parkes data we find evidence of additional extended, low-column-density HI emission around IC 5270, indicating that the clouds are the tip of the iceberg of a larger system of gas surrounding this galaxy. This result adds to the body of evidence on the presence of intra-group gas within the IC 1459 group. Altogether, the HI found outside galaxies in this group amounts to several times $10^9$ M$_\odot$, at least 10% of the HI contained inside galaxies. This suggests a substantial flow of gas in and out of galaxies during the several billion years of the group's evolution.

This paper describes the system architecture of a newly constructed radio telescope - the Boolardy Engineering Test Array, which is a prototype of the Australian Square Kilometre Array Pathfinder telescope. Phased array feed technology is used to form multiple simultaneous beams per antenna, providing astronomers with unprecedented survey speed. The test array described here is a 6-antenna interferometer, fitted with prototype signal processing hardware capable of forming at least 9 dual-polarisation beams simultaneously, allowing several square degrees to be imaged in a single pointed observation. The main purpose of the test array is to develop beamforming and wide-field calibration methods for use with the full telescope, but it will also be capable of limited early science demonstrations.

Using the Very Long Baseline Array and the European VLBI Network, we have observed 22.2 GHz H_2O and 1612 MHz OH masers in the "water fountain" source IRAS 18460-0151. The H_2O maser spectrum has a very wide line-of-sight velocity range (~310 km/s) and consists of three groups of emission features at the blue-shifted (-68 km/s <~ V_LSR <~ -17 km/s) and red-shifted (V_LSR ~= 240 km/s) edges as well as around the systemic velocity (112 km/s <~ V_LSR <~ 133 km/s). The first two H_2O spectral components exhibit a highly-collimated high-velocity bipolar jet on the sky, with an angular separation of ~120 milliarcseconds (mas) (240 AU in linear length) and a three-dimensional flow velocity of ~160 km/s. The flow dynamical age is estimated to be only ~6 yr (at the time of the observation epochs of 2006--2007). Interestingly, the systemic velocity component clearly exhibits a spherically-expanding outflow with a radius of ~36 AU and a flow velocity of ~9 km/s. On the other hand, the OH maser spectrum shows double peaks with a velocity separation of ~25 km/s (V_LSR=$111--116 and 138--141 km/s), as typically seen in circumstellar envelopes of OH/IR stars. The angular offset between the velocity-integrated brightness peaks of the two high-velocity H_2O components is ~25 mas (50 AU). The offset direction and the alignment of the red-shifted maser spots are roughly perpendicular to the axis of the H_2O maser flow. High-accuracy astrometry for the H_2O and OH masers demonstrates that the collimated fast jet and the slowly expanding outflow originate from a single or multiple sources which are located within 15 mas (30 AU). On the other hand, the estimated systemic velocity of the collimated jet (V_sys ~87--113 km/s) has a large uncertainty. This makes it difficult to provide strong constraints on models of the central stellar system of IRAS 18460-0151.

We present the latest and final version of the movie of the SiO masers toward the Mira variable TX Cam. The new version consists of 112 frames (78 successfully reduced epochs) with data covering almost three complete stellar cycles between 24 May 1997 to 25 January 2002, observed with the VLBA. In this paper we examine the global morphology, kinematics and variability of the masering zone. The morphology of the emission is confined in a structure that usually resembles a ring or an ellipse, with occasional deviations due to localised phenomena. The ring appears to be contracting and expanding, although for the first cycle contraction is not observed. The width and outer boundary of the masering zone follow the stellar pulsation. Our data seem to be consistent with a shock being created once per stellar cycle at maximum that propagates with a velocity of ~7 km/s. The difference in velocities along different axes strongly suggests that the outflow in TX Cam is bipolar. The contribution of projection is examined and our results are compared with the latest theoretical model.

Supernova remnants (SNRs), as the major contributors to the galactic cosmic rays (CR), are believed to maintain an average CR spectrum by diffusive shock acceleration (DSA) regardless of the way they release CRs into the interstellar medium (ISM). However, the interaction of the CRs with nearby gas clouds crucially depends on the release mechanism. We call into question two aspects of a popular paradigm of the CR injection into the ISM, according to which they passively and isotropically diffuse in the prescribed magnetic fluctuations as test particles. First, we treat the escaping CR and the Alfven waves excited by them on an equal footing. Second, we adopt field aligned CR escape outside the source, where the waves become weak. An exact analytic self-similar solution for a CR "cloud" released by a dimmed accelerator strongly deviates from the test-particle result. The CR diffusion coefficient $D_{NL}$ is strongly suppressed compared to its background ISM value $D_{ISM}$: $D_{NL}\sim D_{ISM}\exp(-\Pi)<< D_{ISM}$ for sufficiently high field-line-integrated CR partial pressure, $\Pi$. When $\Pi>>1$, the CRs drive Alfven waves efficiently enough to build a transport barrier that strongly reduces the leakage. The solution has a spectral break at $p=p_{br}$, where $p_{br}$ satisfies the following equation $D_{NL}(p_{br})\simeq z^{2}/t$.

Silicon monoxide maser emission has been detected in the circumstellar envelopes of many evolved stars in various vibrationally-excited rotational transitions. It is considered a good tracer of the wind dynamics close to the photosphere of the star. We have investigated the polarization morphology in the circumstellar envelope of an AGB star, R Cas. We mapped the linear and circular polarization of SiO masers in the v=1, J=1-0 transition. The linear polarization is typically a few tens of percent while the circular polarization is a few percent. The fractional polarization tends to be higher for emission of lower total intensity. We found that, in some isolated features the fractional linear polarization appears to exceed 100%. We found the Faraday rotation is not negligible but is ~15 deg., which could produce small scale structure in polarized emission whilst total intensity is smoother and partly resolved out. The polarization angles vary considerably from feature to feature but there is a tendency to favour the directions parallel or perpendicular to the radial direction with respect to the star. In some features, the polarization angle abruptly flips 90 deg. We found that our data are in the regime where the model of Goldreich et al (1973) can be applied and the polarization angle flip is caused when the magnetic field is at close to 55 deg. to the line of sight. The polarization angle configuration is consistent with a radial magnetic field although other configurations are not excluded.

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10 deg < l < 65 deg), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.5" resolution Stokes I map with a root-mean-squared noise level better than 0.4 mJy/beam. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalogue of compact radio-emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimised to detect emission on size scales up to 14" and for unresolved sources the catalogue is more than 90 percent complete at a flux density of 3.9 mJy. We have detected 3,062 sources above a 7-sigma detection limit and present their ensemble properties. The catalogue is highly reliable away from regions containing poorly-sampled extended emission, which comprise less than two percent of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artefacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: regions, planetary nebulae and radio-galaxies. The CORNISH data and catalogue are available online at http://cornish.leeds.ac.uk

We describe the motivation, design and implementation of the CORNISH survey, an arcsecond resolution radio continuum survey of the inner Galactic plane at 5 GHz using the Karl G. Jansky Very Large Array (VLA). It is a blind survey co-ordinated with the northern Spitzer GLIMPSE I region covering 10 deg < l < 65 deg and |b| < 1 deg at similar resolution. We discuss in detail the strategy that we employed to control the shape of the synthesised beam across this survey that covers a wide range of fairly low declinations. Two snapshots separated by 4 hours in hour angle kept the beam elongation to less that 1.5 over 75% of the survey area and less than 2 over 98% of the survey. The prime scientific motivation is to provide an unbiased survey for ultra-compact HII regions to study this key phase in massive star formation. A sensitivity around 2 mJy will allow the automatic distinction between radio loud and quiet mid-IR sources found in the Spitzer surveys. This survey has many legacy applications beyond star formation including evolved stars, active stars and binaries, and extragalactic sources. The CORNISH survey for compact ionized sources complements other Galactic plane surveys that target diffuse and non-thermal sources as well as atomic and molecular phases to build up a complete picture of the ISM in the Galaxy.

We present the results of single-dish and VLBI observations for the water-vapor masers at the nucleus of the Seyfert 2, IC 2560. We monitored velocities of the maser features with the 45-m telescope of the Nobeyama Radio Observatory. Using the data of 1995--2006, the velocity drift rate was detected to be a = +2.57 +/- 0.04 km/s/yr on the average for 6 systemic features. The Very Long Baseline Array (VLBA) with the Very Large Array (VLA) firstly detected a red-shifted and a blue-shifted maser features of IC 2560, in addition to systemic maser features and a continuum component. We propose a maser disk in the nuclear region. The systemic and red-shifted features are emitted from a nearly edge-on disk with the position angle of PA = -46 deg, which is almost perpendicular to the galactic disk. Assuming the Keplerian rotation, the radii of the maser disk are r = 0.087--0.335 pc, and the thickness is 2H <= 0.025 pc. The binding mass is 3.5 x 10^6 Msun at a distance of D = 26 Mpc, and the mean volume density within the inner radius is 1.3 x 10^9 Msun/pc^3, strongly suggesting a massive black hole at the center. A continuum component was detected at the 0.2 pc southwest of the disk center, and considered as a jet ejected from the nucleus, with an angle of 70 deg from the disk. The blue-shifted maser feature is located on the continuum component, being interpreted to be a "jet maser". The distance to IC 2560 is estimated to be D = 31^{+12}_{-14} Mpc from the geometry of the maser disk and the velocity drift rate.

A survey of the Milky Way disk and the Magellanic System at the wavelengths of the 21-cm atomic hydrogen (HI) line and three 18-cm lines of the OH molecule will be carried out with the Australian Square Kilometre Array Pathfinder telescope. The survey will study the distribution of HI emission and absorption with unprecedented angular and velocity resolution, as well as molecular line thermal emission, absorption, and maser lines. The area to be covered includes the Galactic plane (|b|< 10deg) at all declinations south of delta = +40deg, spanning longitudes 167deg through 360deg to 79deg at b=0deg, plus the entire area of the Magellanic Stream and Clouds, a total of 13,020 square degrees. The brightness temperature sensitivity will be very good, typically sigma_T ~ 1 K at resolution 30arcsec and 1 km/s. The survey has a wide spectrum of scientific goals, from studies of galaxy evolution to star formation, with particular contributions to understanding stellar wind kinematics, the thermal phases of the interstellar medium, the interaction between gas in the disk and halo, and the dynamical and thermal states of gas at various positions along the Magellanic Stream.

The proof of cosmic ray (CR) origin in supernova remnants (SNR) must hinge on full consistency of the CR acceleration theory with the observations; direct proof is impossible because of the orbit stochasticity of CR particles. Recent observations of a number of galactic SNR strongly support the SNR-CR connection in general and the Fermi mechanism of CR acceleration, in particular. However, many SNR expand into weakly ionized dense gases, and so a significant revision of the mechanism is required to fit the data. We argue that strong ion-neutral collisions in the remnant surrounding lead to the steepening of the energy spectrum of accelerated particles by \emph{exactly one power}. The spectral break is caused by a partial evanescence of Alfven waves that confine particles to the accelerator. The gamma-ray spectrum generated in collisions of the accelerated protons with the ambient gas is also calculated. Using the recent Fermi spacecraft observation of the SNR W44 as an example, we demonstrate that the parent proton spectrum is a classical test particle power law $\propto E^{-2}$, steepening to $E^{-3}$ at $E_{br}\approx7GeV$.

Context. The two nuclei of the starburst galaxy Arp220 contain multiple compact radio sources previously identified as radio supernovae or supernova remnants. Aims. In order to search for an embedded radio AGN, or other possible exotic objects, we have carried out a program of VLBI monitoring at 6 cm over three epochs each separated by four months. Methods. Combining the new data with existing data at 6 cm and 18 cm (spanning 4 and 12 years respectively) we are able to characterise source flux density variability on a range of time-scales. Additionally we analyse the variability of sources in shape and position. Results. We detect rapid (< 4 months) variability in three sources (W7, W26 and W29). These sources show possible superluminal motion (> 4c) of jet-like features near rapidly varying almost stationary components. These enigmatic sources might be associated with an AGN or a highly beamed microquasar (i.e. microblazar). Other hypotheses include that the apparent variability is intrinsic and is produced by neutron star powered central components within a supernova remnant, by a sequence of several supernovae within super star clusters, or is extrinsic and is produced by Galactic interstellar scintillation of very compact non-varying objects. Conclusions. A microquasar/microblazar origin seems to be the best explanation for the nature of the variable sources in Arp220.

Most of PNe are not spherical. The loss of spherical symmetry occurs somewhere between the AGB and PN phase. The cause of this change of morphology is not yet well known, but magnetic fields are one of the possible agents. Its origin remains to be determined, and potentially requires the presence of a massive companion to the AGB star. Therefore, further detections of the magnetic field around evolved stars (in particular those thought to be part of a binary system) are crucial to improve our understanding of the origin and role of magnetism on evolved stars. One such binaries is the pre-PN OH231.8, around which a magnetic field was detected in the OH maser region of the outer circumstellar envelope. We aim to detect and infer the properties of the magnetic field of this source in the water maser region. We observed the 6_{1,6}-5_{2,3} water maser rotational transition to determine its linear and circular polarization. These emissions are located within the inner regions of OH231.8 (at few tens of AU). We detected 30 water maser features, which occur in two distinct regions that are moving apart with a velocity on the sky of 2.3 mas/year. Taking into account the inclination angle of the source, this corresponds to an average separation velocity of 21 km/s. Based on the velocity gradient of the maser emission, the masers appear to be dragged along the direction of the nebula jet. Linear polarization is present in 3 of the features, and circular polarization was detected in the 2 brightest ones. We found that the strength of the magnetic field is |B_{||}|~45 mG which, when assuming a toroidal magnetic field, implies B~2.5 G on the stellar surface. The morphology of the field is not yet determined, but the high scatter found on the directions of the linear polarization vectors could indicate that the masers occur near the tangent points of a toroidal field.

The much-anticipated proof of cosmic ray (CR) acceleration in supernova remnants (SNR) must hinge on full consistency of acceleration theory with the observations; direct proof is impossible because of the orbit scrambling of CR particles. The PAMELA orbital telescope revealed deviation between helium and proton CR spectra deemed inconsistent with the theory, since the latter does not differentiate between elements of ultrarelativistic rigidity. By considering an initial (injection-) phase of the diffusive shock acceleration (DSA), where elemental similarity does not apply, we demonstrate that the spectral difference is, in fact, a unique signature of the DSA. Collisionless plasma SNR shocks inject more He2+ relative to protons when they are stronger and so produce harder helium spectra. The injection bias is due to Alfven waves driven by the more abundant protons, so the He2+ ions are harder to trap by these waves because of the larger gyroradii. By fitting the p/He ratio to the PAMELA data, we bolster the DSA-case for resolving the century-old mystery of CR origin.

This paper examines the detailed sub-milliarcsecond polarization properties of an individual SiO maser feature displaying a rotation in polarization electric vector position angle of approximately \pi/2 across the feature. Such rotations are a characteristic observational signature of circumstellar SiO masers detected toward a number of late-type, evolved stars. We employ a new calibration method for accurate circular VLBI polarimetry at millimeter wavelengths, to present the detailed Stokes {I,Q,U,V} properties for this feature. We analyze the fractional linear and circular polarization as a function of projected angular distance across the extent of the feature, and compare these measurements against several theoretical models proposed for sharp rotations of electric vector position angle in polarized SiO maser emission. We find that the rotation is most likely caused by the angle \theta between the line of sight and a projected magnetic field crossing the critical Van Vleck angle for maser propagation. The fractional linear polarization profile m_l(\theta) is well-fit by standard models for polarized maser transport, but we find less agreement for the fractional circular polarization profile m_c(\theta).

In shock precursors populated by accelerated cosmic rays (CR), the CR return current instability is believed to significantly enhance the pre-shock perturbations of magnetic field. We have obtained fully-nonlinear exact ideal MHD solutions supported by the CR return current. The solutions occur as localized spikes of circularly polarized Alfven envelopes (solitons, or breathers). As the conventional (undriven) solitons, the obtained magnetic spikes propagate at a speed $C$ proportional to their amplitude, $C=C_{A}B_{{\rm max}}/\sqrt{2}B_{0}$. The sufficiently strong solitons run thus ahead of the main shock and stand in the precursor, being supported by the return current. This property of the nonlinear solutions is strikingly different from the linear theory that predicts non-propagating (that is, convected downstream) circularly polarized waves. The nonlinear solutions may come either in isolated pulses (solitons) or in soliton-trains (cnoidal waves). The morphological similarity of such quasi-periodic soliton chains with recently observed X-ray stripes in Tycho supernova remnant (SNR) is briefly discussed. The magnetic field amplification determined by the suggested saturation process is obtained as a function of decreasing SNR blast wave velocity during its evolution from the ejecta-dominated to the Sedov-Taylor stage.

We present 2 cm and 3.6 cm wavelength very long baseline interferometry images of the compact radio continuum sources in the nearby ultra-luminous infrared galaxy Arp220. Based on their radio spectra and variability properties, we confirm these sources to be a mixture of supernovae (SNe) and supernova remnants (SNRs). Of the 17 detected sources we resolve 7 at both wavelengths. The SNe generally only have upper size limits. In contrast all the SNRs are resolved with diameters {\geq} 0.27 pc. This size limit is consistent with them having just entered their Sedov phase while embedded in an interstellar medium (ISM) of density 10^4 cm^{-3} . These objects lie on the diameter-luminosity correlation for SNRs (and so also on the diameter-surface brightness relation) and extend these correlations to very small sources. The data are consistent with the relation L {\propto} D^{-9/4}. Revised equipartition arguments adjusted to a magnetic field to relativistic particle energy density ratio of 1% combined with a reasonable synchrotron-emitting volume filling factor of 10% give estimated magnetic field strengths in the SNR shells of ~ 15-50 mG. The SNR shell magnetic fields are unlikely to come from compression of ambient ISM fields and must instead be internally generated. We set an upper limit of 7 mG for the ISM magnetic field. The estimated energy in relativistic particles, 2%-20% of the explosion kinetic energy, is consistent with estimates from models that fit the IR-radio correlation in compact starburst galaxies.

We report the results of multi-epoch VLBA observations of the 22.2GHz water maser emission associated with the "water fountain" IRAS 18286-0959. We suggest that this object is the second example of a highly collimated bipolar precessing outflow traced by water maser emission, the other is W43A. The detected water emission peaks are distributed over a velocity range from -50km/s to 150km/s. The spatial distribution of over 70% of the identified maser features is found to be highly collimated along a spiral jet (jet 1) extended southeast to northwest, the remaining features appear to trace another spiral jet (jet 2) with a different orientation. The two jets form a "double-helix" pattern which lies across ~200 milliarcseconds. The maser distribution is reasonably fit by a model consisting of two bipolar precessing jets. The 3D velocities of jet 1 and jet 2 are derived to be 138km/s and 99km/s, respectively. The precession period of jet 1 is about 56 years. For jet 2, three possible models are tested and they give different values for the kinematic parameters. We propose that the appearance of two jets is the result of a single driving source with significant proper motion.

We present multi-epoch, total intensity, high-resolution images of 43GHz, v=1, J=1-0 SiO maser emission toward the Mira variable R Cas. In total we have 23 epochs of data for R Cas at approximate monthly intervals over an optical pulsation phase range from 0.158 to 1.78. These maps show a ring-like distribution of the maser features in a shell, which is assumed to be centred on the star at a radius of 1.6 to 2.3 times the stellar radii. It is clear from these images that the maser emission is significantly extended around the star. At some epochs a faint outer arc can be seen at 2.2 stellar radii. The intensity of the emission waxes and wanes during the stellar phase. Some maser features are seen infalling as well as outflowing. We have made initial comparisons of our data with models by Gray et. al. (2009).

A mechanism for proton acceleration to ~10^21eV is suggested. It may operate in accretion flows onto thin dark matter filaments of cosmic structure formation. The flow compresses the ambient magnetic field to strongly increase and align it with the filament. Particles begin the acceleration by the ExB drift with the accretion flow. The energy gain in the drift regime is limited by the conservation of the adiabatic invariant p_perp^2/B. Upon approaching the filament, the drift turns into the gyro-motion around the filament so that the particle moves parallel to the azimuthal electric field. In this 'betatron' regime the acceleration speeds up to rapidly reach the electrodynamic limit $cp_{max}=eBR$ for an accelerator with magnetic field $B$ and the orbit radius $R$ (Larmor radius). The periodic orbit becomes unstable and the particle slings out of the filament to the region of a weak (uncompressed) magnetic field, which terminates the acceleration. The mechanism requires pre-acceleration that is likely to occur in structure formation shocks upstream or nearby the filament accretion flow. Previous studies identify such shocks as efficient proton accelerators to a firm upper limit ~10^19.5 eV placed by the catastrophic photo-pion losses. The present mechanism combines explosive energy gain in its final (betatron) phase with prompt particle release from the region of strong magnetic field. It is this combination that allows protons to overcome both the photo-pion and the synchrotron-Compton losses and therefore attain energy 10^21 eV. A requirement on accelerator to reach a given E_max placed by the accelerator energy dissipation \propto E_{max}^{2}/Z_0 due to the finite vacuum impedance Z_0 is circumvented by the cyclic operation of the accelerator.

Recent observations of the supernova remnant W44 by the \emph{Fermi} spacecraft observatory strongly support the idea that the bulk of galactic cosmic rays is accelerated in such remnants by a Fermi mechanism, also known as diffusive shock acceleration. However, the W44 expands into weakly ionized dense gas, and so a significant revision of the mechanism is required. In this paper we provide the necessary modifications and demonstrate that strong ion-neutral collisions in the remnant surrounding lead to the steepening of the energy spectrum of accelerated particles by \emph{exactly one power}. The spectral break is caused by Alfven wave evanescence leading to the fractional particle losses. The gamma-ray spectrum generated in collisions of the accelerated protons with the ambient gas is also calculated and successfully fitted to the Fermi Observatory data. The parent proton spectrum is best represented by a classical test particle power law E^-2, steepening to E^-3 at E_br~7GeV due to deteriorated particle confinement.

We present multi-epoch global VLBI observations at 2 cm, 3.6 cm and 6 cm of the compact radio sources in Arp220. We resolve many sources and estimate sizes, expansion velocities and source classes. We find most source properties are consistent with them being either radio supernovae (SNe) or supernova remnants (SNRs). We extend the luminosity-diameter relation for SNRs to very small sources and argue this supports models where shell magnetic fields are internally amplified. We also detect one probable SN/SNR transition object candidate and three highly variable sources with possible superluminal motion (approximately 4c) of jet-like features near rapidly varying almost stationary components. These enigmatic sources, which show similarities to the recently discovered superluminal source in M82, might be associated with an AGN or a new radio source class (e.g. intermediate-mass black holes or micro-blazars).

Both the acceleration of cosmic rays (CR) in supernova remnant shocks and their subsequent propagation through the random magnetic field of the Galaxy deem to result in an almost isotropic CR spectrum. Yet the MILAGRO TeV observatory discovered a sharp ($\sim10^{\circ})$ arrival anisotropy of CR nuclei. We suggest a mechanism for producing a weak and narrow CR beam which operates en route to the observer. The key assumption is that CRs are scattered by a strongly anisotropic Alfven wave spectrum formed by the turbulent cascade across the local field direction. The strongest pitch-angle scattering occurs for particles moving almost precisely along the field line. Partly because this direction is also the direction of minimum of the large scale CR angular distribution, the enhanced scattering results in a weak but narrow particle excess. The width, the fractional excess and the maximum momentum of the beam are calculated from a systematic transport theory depending on a single scale $l$ which can be associated with the longest Alfven wave, efficiently scattering the beam. The best match to all the three characteristics of the beam is achieved at $l\sim1$pc. The distance to a possible source of the beam is estimated to be within a few 100pc. Possible approaches to determination of the scale $l$ from the characteristics of the source are discussed. Alternative scenarios of drawing the beam from the galactic CR background are considered. The beam related large scale anisotropic CR component is found to be energy independent which is also consistent with the observations.

Strong astrophysical shocks, diffusively accelerating cosmic rays (CR) ought to develop CR precursors. The length of such precursor $L_{p}$ is believed to be set by the ratio of the CR mean free path $\lambda$ to the shock speed, i.e., $L_{p}\sim c\lambda/V_{sh}\sim cr_{g}/V_{sh}$, which is formally independent of the CR pressure $P_{c}$. However, the X-ray observations of supernova remnant shocks suggest that the precursor scale may be significantly shorter than $L_{p}$ which would question the above estimate unless the magnetic field is strongly amplified and the gyroradius $r_{g}$ is strongly reduced over a short (unresolved) spatial scale. We argue that while the CR pressure builds up ahead of the shock, the acceleration enters into a strongly nonlinear phase in which an acoustic instability, driven by the CR pressure gradient, dominates other instabilities (at least in the case of low $\beta$ plasma). In this regime the precursor steepens into a strongly nonlinear front whose size scales with \emph{the CR pressure}as $L_{f}\sim L_{p}\cdot(L_{s}/L_{p})^{2}(P_{c}/P_{g})^{2}$, where $L_{s}$ is the scale of the developed acoustic turbulence, and $P_{c}/P_{g}$ is the ratio of CR to gas pressure. Since $L_{s}\ll L_{p}$, the precursor scale reduction may be strong in the case of even a moderate gas heating by the CRs through the acoustic and (possibly also) the other instabilities driven by the CRs.

We present the first view of the magnetic field structure in the OH shell of the extreme OH/IR star OH 26.5+0.6. MERLIN interferometric observations of this object were obtained in December 1993 in full polarisation, at 1612, 1665 and 1667 MHz. The maser spots show a spheroidal distribution both at 1612 and 1667 MHz, while at 1665 MHz emission from the blue-shifted maser peak is concentrated on the stellar position, and the red-shifted peak emission exhibits a filamentary structure oriented on a SE-NW axis. The linear polarisation in both main lines is rather faint, ranging from 9 to 20% at 1665 MHz and from 0 to 30% at 1667 MHz. At 1612 MHz most maser spots exhibit a similar range of linear polarisation although those in the outermost parts of the envelope reach values as high as 66%. This is particularly apparent in the southern part of the shell. The detailed distribution of the polarisation vectors could only be obtained at 1612 MHz. The polarisation vectors show a highly structured distribution indicative of a poloidal magnetic field inclined by 40-60$^\circ$ to the line of sight. The velocity distribution of the maser spots with respect to the radial distance is well explained by an isotropic outflow at constant velocity in the case of a prolate shaped spheroid envelope, also tilted about 45-65$^\circ$ to the line of sight.

We have conducted a Galactic plane survey of methanol masers at 6668 MHz using a 7-beam receiver on the Parkes telescope. Here we present results from the first part, which provides sensitive unbiased coverage of a large region around the Galactic Centre. Details are given for 183 methanol maser sites in the longitude range 345$^{\circ}$ through the Galactic Centre to 6$^{\circ}$. Within 6$^{\circ}$ of the Centre, we found 88 maser sites, of which more than half (48) are new discoveries. The masers are confined to a narrow Galactic latitude range, indicative of many sources at the Galactic Centre distance and beyond, and confined to a thin disk population; there is no high latitude population that might be ascribed to the Galactic Bulge. Within 2$^{\circ}$ of the Galactic Centre the maser velocities all lie between -60 and +77 \kms, a range much smaller than the 540 \kms range observed in CO. Elsewhere, the maser with highest positive velocity (+107 \kms) occurs, surprisingly, near longitude 355$^{\circ}$ and is probably attributable to the Galactic bar. The maser with the most negative velocity (-127 \kms) is near longitude 346$^{\circ}$, within the longitude-velocity locus of the near side of the `3-kpc arm'. It has the most extreme velocity of a clear population of masers associated with the near and far sides of the 3-kpc arm. Closer to the Galactic Centre the maser space density is generally low, except within 0.25 kpc of the Centre itself, the `Galactic Centre Zone', where it is 50 times higher, which is hinted at by the longitude distribution, and confirmed by the unusual velocities.

We have studied the growth of the shell-like radio structure of supernova SN1993J in M81 from September 1993 through October 2003 with very-long-baseline interferometry (VLBI) observations at the wavelengths of 3.6, 6, and 18cm. For this purpose, we have developed a method to accurately determine the outer radius (R) of any circularly symmetric compact radio structure like SN1993J. The source structure of SN1993J remains circularly symmetric (with deviations from circularity under 2%) over almost 4000 days. We characterize the decelerated expansion of SN 1993J through approximately day 1500 after explosion with an expansion parameter $m= 0.845\pm0.005$ ($R \propto t^{m}$). However, from that day onwards the expansion is different when observed at 6 and 18cm. Indeed, at 18cm, the expansion can be well characterized by the same $m$ as before day 1500, while at 6cm the expansion appears more decelerated, and is characterized by another expansion parameter, $m_{6}= 0.788\pm0.015$. Therefore, since about day 1500 on, the radio source size has been progressively smaller at 6cm than at 18cm. These findings are in stark contrast to previous reports by other authors on the details of the expansion. In our interpretation the supernova expands with a single expansion parameter, $m= 0.845\pm0.005$, and the 6cm results beyond day 1500 are due to physical effects, perhaps also coupled to instrumental limitations. Two physical effects may be involved: (a) a changing opacity of the ejecta to the 6cm radiation, and (b) a radial decrease of the magnetic field in the emitting region. (Long abstract cut. Please, read full abstract in manuscript).

We present a time series of synoptic images of the linearly-polarized v=1, J=1-0 SiO maser emission toward the Mira variable, TX Cam. These data comprise 43 individual epochs at an approximate biweekly sampling over an optical pulsation phase range of 0.68 to 1.82. The images have an angular resolution of ~500 microarcsec and were obtained using the Very Long Baseline Array (VLBA), operating in the 43 GHz band in spectral-line, polarization mode. We have previously published the total intensity time series for this pulsation phase range; this paper serves to present the linearly-polarized image sequence and an associated animation representing the evolution of the linear polarization morphology over time. We find a predominantly tangential polarization morphology, a high degree of persistence in linear polarization properties over individual component lifetimes, and stronger linear polarization in the inner projected shell than at larger projected shell radii. We present an initial polarization proper motion analysis examining the possible dynamical influence of magnetic fields in component motions in the extended atmospheres of late-type, evolved stars.

A new 7-beam 6-7 GHz receiver has been built to survey the Galaxy and the Magellanic Clouds for newly forming high-mass stars that are pinpointed by strong methanol maser emission at 6668 MHz. The receiver was jointly constructed by Jodrell Bank Observatory (JBO) and the Australia Telescope National Facility (ATNF) and allows simultaneous coverage at 6668 and 6035 MHz. It was successfully commissioned at Parkes in January 2006 and is now being used to conduct the Parkes-Jodrell multibeam maser survey of the Milky Way. This will be the first systematic survey of the entire Galactic plane for masers of not only 6668-MHz methanol, but also 6035-MHz excited-state hydroxyl. The survey is two orders of magnitude faster than most previous systematic surveys and has an rms noise level of ~0.17 Jy.This paper describes the observational strategy, techniques and reduction procedures of the Galactic and Magellanic Cloud surveys, together with deeper, pointed, follow-up observations and complementary observations with other instruments. It also includes an estimate of the survey detection efficiency. The 111 days of observations with the Parkes telescope have so far yielded >800 methanol sources, of which ~350 are new discoveries. The whole project will provide the first comprehensive Galaxy-wide catalogue of 6668-MHz and 6035-MHz masers.

We present the results of the first complete survey of the Large and Small Magellanic Clouds for 6668-MHz methanol and 6035-MHz excited-state hydroxyl masers. In addition to the survey, higher-sensitivity targeted searches towards known star-formation regions were conducted. The observations yielded the discovery of a fourth 6668-MHz methanol maser in the Large Magellanic Cloud (LMC), found towards the star-forming region N160a, and a second 6035-MHz excited-state hydroxyl maser, found towards N157a. We have also re-observed the three previously known 6668-MHz methanol masers and the single 6035-MHz hydroxyl maser. We failed to detect emission from either transition in the Small Magellanic Cloud. All observations were initially made using the Methanol Multibeam (MMB) survey receiver on the 64-m Parkes telescope as part of the MMB project and accurate positions have been measured with the Australia Telescope Compact Array (ATCA). We compare the maser populations in the Magellanic Clouds with those of our Galaxy and discuss their implications for the relative rates of massive star-formation, heavy metal abundance, and the abundance of complex molecules. The LMC maser populations are demonstrated to be smaller than their Milky Way counterparts. Methanol masers are under-abundant by a factor of ~45, whilst hydroxyl and water masers are a factor of ~10 less abundant than our Galaxy.

NRAO 150 -a compact and bright radio to mm source showing core/jet structure- has been recently identified as a quasar at redshift z=1.52 through a near-IR spectral observation. To study the jet kinematics on the smallest accessible scales and to compute the first estimates of its basic physical properties, we have analysed the ultra-high-resolution images from a new monitoring program at 86 GHz and 43 GHz with the GMVA and the VLBA, respectively. An additional archival and calibration VLBA data set, covering from 1997 to 2007, has been used. Our data shows an extreme projected counter-clock-wise jet position angle swing at an angular rate of up to ~11 deg./yr within the inner ~31 pc of the jet, which is associated with a non-ballistic superluminal motion of the jet within this region. The results suggest that the magnetic field could play an important role in the dynamics of the jet in NRAO 150, which is supported by the large values of the magnetic field strength obtained from our first estimates. The extreme characteristics of the jet swing make NRAO 150 a prime source to study the jet wobbling phenomenon.

MERLIN observations of 6.668-GHz methanol and both 6.031- and 6.035-GHz hydroxyl (OH) emission from the massive star-formation region ON1 are presented. These are the first methanol observations made in full polarization using 5 antennas of MERLIN, giving high resolution and sensitivity to extended emission. Maser features are found to lie at the southern edge of the ultra-compact HII region, following the known distribution of ground-state OH masers. The masers cover a region ~1 arcsec in extent, lying perpendicular to the H13CO+ bipolar outflow. Excited-state OH emission demonstrates consistent polarization angles across the strongest linearly polarized features which are parallel to the overall distribution. The linear polarizations vary between 10.0 and 18.5 per cent, with an average polarization angle of -60 deg +/- 28 deg. The strongest 6.668-GHz methanol features provide an upper limit to linear polarization of ~1 per cent. Zeeman splitting of OH shows magnetic fields between -1.1 to -5.8 mG, and a tentative methanol magnetic field strength of -18 mG is measured.

A gas cloud near a supernova remnant (SNR) provides a target for pp-collisions leading to subsequent gamma-ray emission through neutral pion decay. The assumption of a power-law ambient spectrum of accelerated particles with index near -2 is usually built into models predicting the spectra of very-high energy (VHE) gamma-ray emission from SNRs. However, if the gas cloud is located at some distance from the SNR shock, this assumption is not necessarily correct. In this case, the particles which interact with the cloud are those leaking from the shock and their spectrum is approximately monoenergetic with the injection energy gradually decreasing as the SNR ages. In the GLAST energy range the gamma-ray spectrum resulting from particle interactions with the gas cloud will be flatter than expected, with the cutoff defined by the pion momentum distribution in the laboratory frame. We evaluate the flux of particles escaping from a SNR shock and apply the results to the VHE diffuse emission detected by the HESS at the Galactic centre.

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 report the first detection at multiple radio wavelengths (13, 6 and 3.6 cm) of the compact sources within both nuclei of the Ultra Luminous Infra-Red Galaxy Arp 220. We present the radio spectra of the 18 detected sources. In just over half of the sources we find that these spectra and other properties are consistent with the standard model of powerful Type IIn supernovae interacting with their pre-explosion stellar wind. The rate of appearance of new radio sources identified with these supernova events suggests that an unusually large fraction of core collapse supernovae in Arp 220 are highly luminous; possibly implying a radically different stellar initial mass function or stellar evolution compared to galactic disks. Another possible explanation invokes very short (~3 x 10^5 year) intense (~10^3 M_Sol year^-1) star formation episodes with a duty cycle of ~10%. A second group of our detected sources, consisting of the brightest and longest monitored sources at 18 cm do not easily fit the radio supernova model. These sources show a range of spectral indexes from -0.2 to -1.9. We propose that these are young supernova remnants which have just begun interacting with a surrounding ISM with a density between 10^4 and 10^5 cm^-3. One of these sources is probably resolved at 3.6 cm wavelength with a diameter 0.9 pc. In the western nucleus we estimate that the ionized component of the ISM gives rise to foreground free-free absorption with opacity at 18 cm of <0.6 along the majority of lines of sight. Other sources may be affected by absorption with opacity in the range 1 to 2. These values are consistent with previous models as fitted to the radio recombination lines and the continuum spectrum.

We present observations of the 22 GHz water megamasers in the circumnuclear disk of NGC 3079 obtained with the Green Bank Telescope. The data are analyzed for circular polarization due to the Zeeman-induced splitting of the water maser lines. No circular polarization is detected and we derive a 1 sigma upper limit of 11 mG for the toroidal magnetic field at ~0.64 pc from the central black hole. This is the tightest upper limit for the magnetic field around a black hole to date. We use the magnetic field limit to derive an estimate of the mass accretion onto the central black hole. In addition to the polarimetric results, we present an observation of rapid variability in the maser lines, which we explain as weak interstellar scintillation. From the scintillation parameters, we estimate an intrinsic size of the mostly saturated maser features of ~12 microarcseconds. This is consistent with models assuming a thick, clumpy accretion disk.

Water maser polarization observations in the precessing jet of W43A have revealed that it is magnetically collimated. Here we present a detailed description of the physical properties of the water maser environment in the jet. We discuss the maser saturation level and beaming angle as well as the intrinsic temperatures and densities. Additionally, we show that the polarization angle of the strongest red-shifted maser feature undergoes a fast rotation of 90 degrees across the maser. Along with the variation of linear polarization fraction, this strongly supports the current theoretical description of maser linear polarization.

We present a theory for the generation of mesoscale ($kr_{g}\ll 1$, where $r_{g}$ is the cosmic ray gyroradius) magnetic fields during diffusive shock acceleration. The decay or modulational instability of resonantly excited Alfven waves scattering off ambient density perturbations in the shock environment naturally generates larger scale fields. For a broad spectrum of perturbations, the physical mechanism of energy transfer is random refraction, represented by diffusion of Alfven wave packet in $k-$space. The scattering field can be produced directly by the decay instability or by the Drury instability, a hydrodynamic instability driven by the cosmic ray pressure gradient. This process is of interest to acceleration since it generates waves of longer wavelength, and so enables the confinement and acceleration of higher energy particles. This process also limits the intensity of resonantly generated turbulent magnetic field on $r_{g}$ scales.

We present the first direct measurements of the magnetic field strength and direction in a collimated jet from an evolved star on its way to become a planetary nebula. Very Long Baseline Array (VLBA) observations of the linear and circular polarization of the water masers in the collimated jet of W43A reveal a strong toroidal magnetic field, indicating that the jet is magnetically collimated. The magnetic field strength in the jet extrapolated back to the stellar surface yields a surface field of several Gauss, consistent with the measurements of maser polarization in a large sample of evolved stars. The origin of the magnetic field is yet unknown, although the jet precession might point to the existence of a heavy planet or stellar companion. This is the first direct observational evidence for magnetic collimation in the jets, that likely play an important role in shaping planetary nebulae.

We have used a Very Long Baseline Interferometry (VLBI) array at 18cm wavelength to image the nucleus of the luminous IR galaxy Arp 220 at ~1 pc linear resolution, and with very high sensitivity. The resulting map has an rms of 5.5 microJy/beam, and careful image analysis results in 49 confirmed point sources ranging in flux density from 1.2 mJy down to ~60 microJy. Comparison with high sensitivity data from 12 months earlier reveals at least four new sources. The favored interpretation of these sources is that they are radio supernovae, and if all new supernovae are detectable at this sensitivity, a resulting estimate of the supernova rate in the Arp 220 system is 4 +/- 2 per year. The implied star formation rate is sufficient to power the entire observed far-infrared luminosity of the galaxy. The two nuclei of Arp 220 exhibit striking similarities in their radio properties, though the western nucleus is more compact, and appears to be ~3 times more luminous than the eastern nucleus. There are also some puzzling differences, and differential free-free absorption, synchrotron aging and expansion losses may all be playing a role. Comparison with the nearby starburst galaxy M82 supports the hypothesis that the activity in Arp 220 is essentially a scaled-up version of that in M82.

Planetary nebulae often have asymmetric shapes, which could arise due to collimated jets from evolved stars before evolution to the planetary nebula phase. The source of jet collimation in these stars is unknown. Magnetic fields are thought to collimate outflows that are observed in many other astrophysical sources, such as active galactic nuclei and proto-stars, although hitherto there are no direct observations of both the magnetic field direction and strength in any collimated jet. Theoretical models have shown that magnetic fields could also be the dominant source of collimation of jet in evolved stars. Here we report measurements of the polarization of water vapour masers that trace the precessing jet emanating from the asymptotic giant branch star W43A at 2.6 kpc from the Sun, which is undergoing rapid evolution into a planetary nebula. The masers occur in two clusters at opposing tips of the jets, ~1,000 AU from the star. We find direct evidence that the magnetic field is collimating the jet.

We present linear and circular polarization observations of the water masers in 4 distinct regions spread over 1x2 arcseconds around the HW2 high-mass young stellar object in the Cepheus A star-forming region. We find magnetic fields between 100-500 mG in the central maser region, which has been argued to trace a circumstellar disk. The masers further from HW2 have field strengths between 30-100 mG. In all cases the magnetic field pressure is found to be similar to the dynamic pressure, indicating that the magnetic field is capable of controlling the outflow dynamics around HW2. In addition to several water maser complexes observed before, we also detect a new maser filament, ~1 arcsec (~690 AU) East of HW2, which we interpret as a shocked region between the HW2 outflow and the surrounding medium. We detect a linear polarization gradient along the filament as well as a reversal of the magnetic field direction. This is thought to mark the transition between the magnetic field associated with the outflow and that found in the surrounding molecular cloud. In addition to the magnetic field we determine several other physical properties of the maser region, including density and temperatures as well as the maser beaming angles.

We suggest a physical mechanism whereby the acceleration time of cosmic rays by shock waves can be significantly reduced. This creates the possibility of particle acceleration beyond the knee energy at ~10^15eV. The acceleration results from a nonlinear modification of the flow ahead of the shock supported by particles already accelerated to the knee momentum at p ~ p_*. The particles gain energy by bouncing off converging magnetic irregularities frozen into the flow in the shock precursor and not so much by re-crossing the shock itself. The acceleration rate is thus determined by the gradient of the flow velocity and turns out to be formally independent of the particle mean free path (m.f.p.). The velocity gradient is, in turn, set by the knee-particles at p ~ p_* as having the dominant contribution to the CR pressure. Since it is independent of the m.f.p., the acceleration rate of particles above the knee does not decrease with energy, unlike in the linear acceleration regime. The reason for the knee formation at p ~ p_* is that particles with $p > p_*$ are effectively confined to the shock precursor only while they are within limited domains in the momentum space, while other particles fall into ``loss-islands'', similar to the ``loss-cone'' of magnetic traps. This structure of the momentum space is due to the character of the scattering magnetic irregularities. They are formed by a train of shock waves that naturally emerge from unstably growing and steepening magnetosonic waves or as a result of acoustic instability of the CR precursor. These losses steepen the spectrum above the knee, which also prevents the shock width from increasing with the maximum particle energy.

We present the analysis of the circular polarization, due to Zeeman splitting, of the water masers around a sample of late-type stars to determine the magnetic fields in their circumstellar envelopes. The magnetic field strengths in the water maser regions around the Mira variable stars U Ori and U Her are shown to be several Gauss while those of the supergiants S Per, NML Cyg and VY CMa are several hundred mG. We also show that large scale magnetic fields permeate the CSE of an evolved star; the polarization of the water masers around VX Sgr reveals a dipole field structure. We shortly discuss the coupling of the magnetic field with the stellar outflow, as such fields could possibly be the cause of distinctly aspherical mass-loss and the resulting aspherical planetary nebulae.

The problem of the cosmic ray origin is discussed in connection with their acceleration in supernova remnant shocks. The diffusive shock acceleration mechanism is reviewed and its potential to accelerate particles to the maximum energy of (presumably) galactic cosmic rays ($10^{18}eV $) is considered. It is argued that to reach such energies, a strong magnetic field at scales larger than the particle gyroradius must be created as a result of the acceleration process, itself. One specific mechanism suggested here is based on the generation of Alfven wave at the gyroradius scale with a subsequent transfer to longer scales via interaction with strong acoustic turbulence in the shock precursor. The acoustic turbulence in turn, may be generated by Drury instability or by parametric instability of the Alfven waves.

Cosmic rays (CRs) are thought to be accelerated in SNRs. The most favorable situation for proving that the main, hadronic CR component is accelerated there is when CRs interact with dense gases, such as molecular clouds (MC) which surround the SN shock. Here, a new mechanism of spectrum formation in partially ionized gases near SNRs is proposed. Using an analytic model of nonlinear diffusive shock acceleration, we calculate the spectra of protons and estimate the resulting $\gamma$-ray emission occurring when the SNR shock approaches a MC. We show that the spectrum develops a break in the TeV range and that its GeV component is suppressed. These modifications to the standard theory occur because of the proximity of the partially ionized MC-gas and because of the physics of particle and Alfven wave propagation inside the gas. Possible applications of the new spectra to the recent CANGAROO and HESS observations of the SNR RX J1713.7-3946 are discussed.

We observed v=1, J=1-0 43-GHz SiO maser emission toward the Mira variable IK Tauri (IK Tau) using the Very Long Baseline Array (VLBA). The images resulting from these observations show that SiO masers form a highly elliptical ring of emission approximately 58 x 32 mas with an axial ratio of 1.8:1. The major axis of this elliptical distribution is oriented at position angle of ~59 deg. The line-of-sight velocity structure of the SiO masers has an apparent axis of symmetry consistent with the elongation axis of the maser distribution. Relative to the assumed stellar velocity of 35 km/s, the blue- and red-shifted masers were found to lie to the northwest and southeast of this symmetry axis respectively. This velocity structure suggests a NW-SE rotation of the SiO maser shell with an equatorial velocity, which we determine to be ~3.6 km/s. Such a NW-SE rotation is in agreement with a circumstellar envelope geometry invoked to explain previous water and OH maser observations. In this geometry, water and OH masers are preferentially created in a region of enhanced density along the NE-SW equator orthogonal to the rotation/polar axis suggested by the SiO maser velocities.

Through polarization observations, circumstellar masers are excellent probes of the magnetic field in the envelopes of late-type stars. Whereas observations of the polarization of the SiO masers close to the star and on the OH masers much further out were fairly commonplace, observations of the magnetic field strength in the intermediate density and temperature region where the 22 GHz water masers occur have only recently become possible. Here we present the analysis of the circular polarization, due to Zeeman splitting, of the water masers around the Mira variable stars U Her and U Ori and the supergiant VX Sgr. We present an upper limit of the field around U Her that is lower but consistent with previous measurements, reflecting possible changes in the circumstellar envelope. The field strengths around U Ori and VX Sgr are shown to be of the order of several Gauss. Moreover, we show for the first time that large scale magnetic fields permeate the circumstellar envelopes of an evolved star; the polarization of the water masers around VX Sgr reveals a dipole field structure. We discuss the coupling of the magnetic field with the stellar outflow, as such fields could possibly be the cause of distinctly aspherical mass-loss.

Many observations of late-type M stars show large near-spherical circumstellar envelopes, yet planetary nebulae and supernova remnants are frequently axisymmetric. We present VLBI and MERLIN observations of masers around the red supergiant S Per which show varying degrees of axisymmetry and a dynamically significant magnetic field. There is no evidence for rotation here or in most similar objects. We examine possible origins of the magnetic field.

We present the results of simultaneous high resolution observations of v=1 and v=2, J=1-0 SiO masers toward TX Cam at four epochs covering a stellar cycle. Near maser maximum (Epochs III and IV), the individual components of both masers are distributed in ring-like structures but the ring is severely disrupted near stellar maser minimum (Epochs I and II). In Epochs III and IV there is a large overlap between the radii at which the two maser transitions occur. However in both epochs the average radius of the v=2 maser ring is smaller than for the v=1 maser ring, the difference being larger for Epoch IV. The observed relative ring radii in the two transitions, and the trends on the ring thickness, are close to those predicted by the model of Humphreys et al. (\cite{humphreys02}). In many individual features there is an almost exact overlap in space and velocity of emission from the two transitions, arguing against pure radiative pumping. At both Epochs III and IV in many spectral features only 50% of the flux density is recovered in our images, implying significant smooth maser structure. For both transitions we find that red- and blue-shifted masers occur in all parts of the rings, with relatively few masers at the systemic velocity. Thus there is no evidence for rotation, although the blue-shifted masers are somewhat more prominent to the west. At all four epochs red-shifted components are generally brighter than blue-shifted ones. At Epochs III and IV, we see many filamentary or spoke-like features in both v=1 and v=2 masers, especially in the red-shifted gas. These spokes show systematic velocity gradients consistent with a decelerating outward flow with increasing radius. We outline a possible model to explain why, given the presence of these spokes, there is a deficit of maser features at the systemic velocity.

We present sub-milliarcsecond observations of SiO masers in the late-type stars IRC +10011 and Chi Cyg. We have used the NRAO Very Long Baseline Array (VLBA) to map the 43 GHz (v=1, 2 J=1-0) and the 86 GHz (v=1, 2 J=2-1) SiO masers. All the transitions have been imaged except the v=2 J=2-1 in IRC +10011. We report the first VLBI map of the v=1 J=2-1 28SiO maser in IRC +10011 as well as the first VLBA images of SiO masers in an S-type Mira variable, Chi Cyg. In this paper we have focused on the study of the relative spatial distribution of the different observed lines. We have found that in some cases the observational results are not reproduced by the current theoretical pumping models, either radiative or collisional. In particular, for IRC +10011, the v=1 J=1-0 and J=2-1 28SiO lines have different spatial distributions and emitting region sizes, the J=2-1 emission being located in an outer region of the envelope. For Chi Cyg, the distributions also differ, but the sizes of the masing regions are comparable. We suggest that the line overlaps between ro-vibrational transitions of two abundant molecular species, H2O and 28SiO, is a possible explanation for the discrepancies found between the observations and the theoretical predictions. We have introduced this overlapping process in the calculations of the excitation of the SiO molecule. We conclude that the line overlaps can strongly affect the excitation of SiO and may reproduce the unexpected observational results for the two sources studied.

The high resolution obtained through the use of VLBI gives an unique opportunity to directly observe the interaction of an expanding radio supernova with its surrounding medium. We present here results from our VLBI observations of the young supernovae SN 1979C, SN 1986J, and SN 2001gd.

On this paper, we present simultaneous VLBA observations of the J=1-0 and J=2-1 rotational lines in the v=1 and v=2 vibrationally excited states of SiO (at 7 and 3 mm wavelengths). We have mapped these four maser lines in the circumstellar envelopes of three AGB stars: IRC +10011, Chi Cyg and TX Cam. We study the relative spatial distribution between these maser lines. In particular, for the two v=1 transitions, we found that the J=1-0 and J=2-1 maser spot distributions are unalike, challenging all the current theoretical pumping mechanisms for SiO masers (which predict quite similar distribution for maser lines in the same vibrational excited state).

We describe an observing campaign using the Very Long Baseline Array (VLBA) to monitor the time-evolution of the v=1, J=1-0 SiO maser emission towards the Mira variable TX Cam. The data reported here cover the period 1997 May 24 to 1999 February 19, during which the SiO maser emission was imaged at approximately bi-weekly intervals. The result is an animated movie at an angular resolution of ~500 microarcsec, over a full pulsation period, of the gas motions in the near circumstellar environment of this star, as traced by the SiO maser emission. This paper serves to release the movie and is the first in a series concerning the scientific results from this observing campaign. In this paper, we discuss the global proper motion of the SiO maser emission as a function of pulsation phase. We measure a dominant expansion mode between optical phases $\phi \sim 0.7-1.5$ confirming ballistic deceleration, and compare this to predictions from existing pulsation models for late-type, evolved stars. Local infall and outflow motions are superimposed on the dominant expansion mode, and non-radial local gas motions are also evident for individual SiO maser components. The overall morphology and evolution of the SiO emission deviates significantly from spherical symmetry, with important implications for models of pulsation kinematics in the near-circumstellar environments of Mira variables.

High-mass star formation is not well understood chiefly because examples are deeply embedded, relatively distant, and crowded with sources of emission. Using VLA and VLBA observations of water and SiO maser emission, we have mapped in detail the structure and proper motion of material 20-500 AU from the closest high-mass YSO, radio source-I in the Orion KL region. We observe streams of material driven in a rotating, wide angle, bipolar wind from the surface of an edge-on accretion disk. The example of source-I provides strong evidence that high-mass star formation proceeds via accretion

Water vapor maser emission in the 6_{16}-5_{23} transition towards the narrow-line Seyfert 1 (NLS1) galaxy NGC 4051 has been discovered during an ongoing single-dish extragalactic water maser survey. The Doppler-shifted maser components appear to bracket maser components lying near the systemic velocity of the galaxy symmetrically. The tentative result of a Very Large Array (VLA) snapshot observation is that the masers are confined within 0.1 arcsec (5 pc at a distance D = 9.7 Mpc) of the radio continuum peak seen at 8.4 GHz. The low luminosity of the maser (~2 Lsun) is not typical for masers that coincide with the radio continuum nucleus and appear associated with Active Galactic Nucleus (AGN) activity. A low-luminosity maser in a Type 1 Seyfert nucleus could be explained by a low maser gain resulting from the lower inclination of an obscuring disk around an active nucleus.

With the VLBA, we have observed water masers associated with the OH/IR star, W43A, which trace a birth and growth of a molecular jet. The water masers exhibit the collimated distribution (1700 AU: 20 AU) and fast motions (+/-150 km/s). The maser distribution was well fit by a precessing jet model. The jet length has extended by 800 AU within 8 years, indicating that the extension rate is roughly equivalent to the jet speed (150 km/s). Very likely the jet was born around the year 1960. An elongated planetary nebula will be formed by such a jet during an extremely short period (<<1000 yrs) of the transition to form a ptoto-planetary nebula. Then fading of the water maser jet will be observed during our human life.

The OH main line masers in circumstellar shells were observed with astrometric VLBI in order to accurately determine the position of the star, as well as its motion, including the parallax. Results are presented on 4 Mira stars monitored for 2 - 8 years with the VLBA. The data show that in some stars the VLBI detection is dominated by blue-shifted emission that is associated with the stellar image amplified by the maser shell in front of the star. In other cases the maser is not directly tied to the stellar position, but still persistent enough to measure proper motion and parallax.

We present circular polarization measurements of circumstellar water masers around 2 Mira variables and 4 supergaint stars. Both an LTE and non-LTE method are used to analyze the circular polarization and total intensity spectra and to obtain the magnetic field strength. The non-LTE method is found to be able to reproduce the observations best. Using this, we find fields from a few hundred milliGauss up to a few Gauss, indicating a solar-type r^(-2) dependence of the magnetic field on the distance, which leads to stellar surface magnetic fields of up to several hundred Gauss. No linear polarization is detected to less than 1.5%.

The main-line OH masers around 4 AGB stars have been observed with the NRAO Very Long Baseline Array (VLBA) at 8 epochs over a period of 2.5 years. Using a phase referencing technique, the position of the most compact maser spot of each star was monitored with respect to two extragalactic reference sources. For U Her and W Hya, we observe the most blue-shifted maser spot, while for R Cas and S CrB we only detect a compact red-shifted maser spot. We managed to determine an accurate proper motion and parallax for U Her, R Cas and S CrB, while additional motion of the compact blue-shifted maser of W Hya is shown to possibly be related to the stellar pulsation. The motion and radio position are compared with the stellar trajectory and absolute optical position determined by the Hipparcos satellite. For U Her and W Hya, the most blue-shifted maser is consistent with the amplified stellar image. The new distances are compared with several published P--L relations, and in this respect the VLBI distances seem an improvement upon the Hipparcos distances.

We report 18cm VLBI continuum imaging observations at 5 mas resolution for UGC 5101, NGC 7469, and Mrk 231, all part of a sample of Luminous Infrared Galaxies which have been shown to have strong VLBI radio cores. The radio morphology of these three systems on VLBI scales is AGN-like, with well-defined ridgelines and high-brightness yet spatially resolved components. The structure and flux densities of these VLBI components are not consistent with starburst generated radio supernovae of the type found in Arp 220. On scales of 100pc the radio continuum in all three objects appears to be dominated by an AGN, not a starburst. Radio emission on larger scales may well originate in a less compact circumnuclear star-forming region. Confirming and extending VLBI imaging of Mrk 231 by Ulvestad et al. (1999), our continuum image shows a triple structure, with a core and two lobes, classifying it as a Compact Symmetric Object (CS0). If the southern (primary) lobe/hot-spot in Mrk 231 is confined by ram pressure, we estimate a lobe advance speed, $v_a \sim 10^{-4}c$, and an age for the jet/compact source, $< 10^6 yr$. We have also imaged the 1667 MHz OH maser emission in Mrk 231, which is extended on scales of 50--100 milliarcsec (40--80 pc) and probably coincides with the inner region of the disk which is seen in CO emission and HI absorption. Among OH megamasers studied at high sensitivity with mas resolution, Mrk 231 is unique in the stringent upper limits placed upon the flux density of compact OH structures of the type found in Arp 220 and other LIGs. It is possible that the circumnuclear environment of Mrk 231 has been sufficiently disrupted by the emergent QSO that the cool, dense clouds necessary for such compact masers no longer exist.

We performed VLA observations of the 22 GHz H2O maser emission in the merging galaxy NGC6240, which hosts the well-known double active nuclei. In a previous paper, we reported on the first solid detection of the H2O maser emission in 2001. After two abortive attempts due to the weakness and probable variability of the emission, the maser was detected with the VLA in June 2002. The emission is unresolved at ~ 0.3 arcsecond and coincides with the southern 22 GHz continuum peak to ~ 0.007 arcsecond (~ 3 pc: D = 97 Mpc). The detection of the maser in the southern nucleus indicates that nuclear activity of the galaxy,which is significant in X-ray and far-infrared (FIR) bands, lies mainly in the southern nucleus, and the nucleus without a high brightness peak could be explained by thick dust emitting FIR radiation. We favour the idea that the maser in NGC6240 is associated with the AGN-activity.

We present MERLIN observations of the continuum (both 1.6 and 5 GHz) and OH maser emission towards Arp220. the correct spatial configuration of the various componnents of the galaxy is revealed. In the eastern component the masers are shown to be generally coincident with the larger scale continuum emission; in the west, the masers and continuum do not generally arise from the same location. A velocity gradient (0.32+/-0.03km/s/pc) is found in the eastern nuclear region in MERLIN scales; this gradient is three times smaller than seen in OH and implies that the OH gas lies inside the HI. A re-analysis of previously presented global VLBI data (Lonsdale et al. 1998) reveals a very high velocity gradient (18.67+/-0.12km/s/pc) in one component, possibly the site of a heavily obscured AGN.

We report observations of water masers around the semiregular variable RT Virginis (RT Vir), which have been made with the Very Long Baseline Array (VLBA) of the National Radio Astronomy Observatory (NRAO) at five epochs, each separated by three weeks of time. We detected about 60 maser features at each epoch. Overall, 61 features, detected at least twice, were tracked by their radial velocities and proper motions. The 3-D maser kinematics exhibited a circumstellar envelope that is expanding roughly spherically with a velocity of about 8 km/s. Asymmetries in both the spatial and velocity distributions of the maser features were found in the envelope, but less significant than that found in other semiregular variables. Systematic radial-velocity drifts of individual maser features were found with amplitudes of <= 2 km/s/yr. For one maser feature, we found a quadratic position shift with time along a straight line on the sky. This apparent motion indicates an acceleration with an amplitude of 33 km/s/yr, implying the passage of a shock wave driven by the stellar pulsation of RT Vir. The acceleration motion is likely seen only on the sky plane because of a large velocity gradient formed in the accelerating maser region. We estimated the distance to RT Vir to be about 220 pc on the basis of both the statistical parallax and model-fitting methods for the maser kinematics.

The diffusion of uni-directional magnetic fields by two dimensional turbulent flows in a weakly ionized gas is studied. The fields here are orthogonal to the plane of fluid motion. This simple model arises in the context of the decay of the mean magnetic flux to mass ratio in the interstellar medium. When ions are strongly coupled to neutrals, the transport of a large--scale magnetic field is driven by both turbulent mixing and nonlinear, ambipolar drift. Using a standard homogeneous and Gaussian statistical model for turbulence, we show rigorously that a large-scale magnetic field can decay on at most turbulent mixing time scales when the field and neutral flow are strongly coupled. There is no enhancement of the decay rate by ambipolar diffusion. These results extend the Zeldovich theorem to encompass the regime of two dimensional flows and orthogonal magnetic fields, recently considered by Zweibel (2002). The limitation of the strong coupling approximation and its implications are discussed.

The 22 GHz H_2O masers in the circumstellar envelope of the Mira variable star U Her have been observed with MERLIN using a phase referencing technique to determine accurate astrometric positions. The positions were compared with the optical positions obtained with the Hipparcos satellite to an accuracy of 18 mas. The absolute radio position of the brightest H_2O maser spot is found to match the optical position, indicating that this spot is the stellar image amplified by the maser screen in front of it. The occurrence of an amplified image in the 22 GHz maser can be used to accurately determine the positions of the H_2O with respect to the star as well as with respect to the SiO and OH masers. Our observations seem to indicate that the star is not in the centre of the distribution of maser spots, which has been interpreted as a ring.

We present circular polarization measurements of circumstellar H_2O masers. The circular polarization detected in the (6_{16}-5_{23}) rotational transition of the H_{2}O maser can be attributed to Zeeman splitting in the intermediate temperature and density regime. The magnetic fields are derived using a general, LTE Zeeman analysis as well as a full radiative transfer method (non-LTE), which includes a treatment of all hyperfine components simultaneously as well as the effects of saturation and unequal populations of the magnetic substates. The differences and relevances of these interpretations are discussed extensively. The field strengths are compared with previous detections of the magnetic field on the SiO and OH masers. We show that the magnetic pressure dominates the thermal pressure by a factor of 20 or more.

We present the first circular polarization measurements of circumstellar water masers around a sample of late-type stars. These observations are used to obtain the magnetic field strength in the water maser region with both an LTE and non-LTE analysis. We find fields from a few hundred milliGauss up to a few Gauss, indicating a solar-type r^(-2) dependence of the magnetic field on the distance. No linear polarization is detected to less than 1%.

SN1993J is to date the radio supernova whose evolution has been monitored in greatest detail and the one which holds best promise for a comprehensive theoretical-observational analysis. The shell-like radio structure of SN1993J has expanded in general accord with models of shock excited emission, showing almost circular symmetry for over 8 years, except for a bright feature at the south-eastern region of the shell that has been observed at every epoch. The spectrum of SN1993J has flattened from alpha =-1 to alpha =-0.67 (S_(\nu) propto nu**(alpha)). The decelerated expansion can be modeled well with a single slope but apparently better with two slopes. There are also intriguing hints of structure in the expansion curve. The results by the two VLBI groups carrying out this research show general agreement, but also some differences. A comparison of the optical and VLBI results about the details of the deceleration show some discrepancies.

We present the first VLBI images of SiO masers in the circumstellar envelope of an S-type star obtained with the newly developed capabilities of the Very Long Baseline Array (VLBA) at 86 GHz. We combine these data with those obtained quasy simultaneously at 43 GHz. These observations provide information on the structure and dynamics of the innermost circumstellar shells, where the return of large quantities of stellar material to the interstellar medium starts. Despite of the fact that these are prelimininary results, we report that the maser emission of the $v$=1 J=1--0 and J=2--1 present a very different emission distribution, being not coincident in most case.

A stellar hydrodynamic pulsation model has been combined with a SiO maser model in an attempt to calculate the temporal variability of SiO maser emission in the circumstellar envelope (CE) of a model AGB star. This study investigates whether the variations in local physical conditions brought about by shocks are the predominant contributing factor to SiO maser variability because, in this work, the radiative part of the pump is constant. We find that some aspects of the variability are not consistent with a pump provided by shock-enhanced collisions alone. In these simulations, gas parcels of relatively enhanced SiO abundance are distributed in a model CE by a Monte Carlo method, at a single epoch of the stellar cycle. From this epoch on, Lagrangian motions of individual parcels are calculated according to the velocity fields encountered in the model CE during the stellar pulsation cycle. The potentially masing gas parcels therefore experience different densities and temperatures, and have varying line-of-sight velocity gradients throughout the stellar cycle, which may or may not be suitable to produce maser emission. At each epoch (separated by 16.6 days), emission lines from the parcels are combined to produce synthetic spectra and VLBI-type images. We report here the results for v=1, J=1-0 (43-GHz) and J=2-1 (86-GHz) masers.

We observed SN1979C in M100 on 4 June 1999, about twenty years after explosion, with a very sensitive four-antenna VLBI array at the wavelength of 18cm. The distance to M100 and the expansion velocities are such that the supernova cannot be fully resolved by our Earth-wide array. Model-dependent sizes for the source have been determined and compared with previous results. We conclude that the supernova shock was initially in free expansion for 6 +/- 2 yrs and then experienced a very strong deceleration. The onset of deceleration took place a few years before the abrupt trend change in the integrated radio flux density curves. We estimate the shocked swept-up mass to be about 1.6 solar masses, assuming a standard density profile for the CSM. Such a swept-up mass for SN1979C suggests a mass of the hydrogen-rich envelope ejected at explosion no larger than $M_{\rm env} \sim$ 0.9 solar masses. If SN1979C originated in a binary star the low value of M_env suggests that the companion of the progenitor star stripped off most of the hydrogen-rich envelope mass of the presupernova star prior to the explosion.

We report the result of an on-going survey for 22 GHz H2O maser emission towards infrared luminous galaxies. The observed galaxies were selected primarily from the IRAS bright galaxy sample. The survey has resulted in the detection of one new maser. The new maser was discovered towards the [U]LIRG/merger galaxy NGC6240, which contains a LINER nucleus. This is the first detection of an H2O maser towards this class of galaxy, they are traditionally associated with OH megamaser sources. The detected maser emission is highly redshifted (~ 260-300 km/s) with respect to the adopted systemic velocity of the galaxy, and we identified no other significant emission at velocities < +- 500 km/s relative to the systemic velocity. The presence of high-velocity maser emission implies the possible existence of a rotating maser disk formed in the merging process. The large maser luminosity (~ 40 Lo) suggests that an active galactic nucleus could be the energy source that gives rise to the water emission. Alternatively, the maser emission could be associated with the previously observed double radio source in the centre of the galaxy. Interferometric observations with high angular resolution will be able to clarify the origin of the new maser.

The theory of shock acceleration predicts the maximum particle energy to be limited only by the acceleration time and the size (geometry) of the shock. This led to optimistic estimates for the galactic cosmic ray energy achievable in the SNR shocks. The estimates imply that the accelerated particles, while making no strong impact on the shock structure (test particle approach) are nevertheless scattered by strong self-generated Alfven waves (turbulent boost) needed to accelerate them quickly. We demonstrate that these two assumptions are in conflict when applied to SNRs of the age required for cosmic ray acceleration to the ``knee'' energy. We study the combined effect of acceleration nonlinearity (shock modification by accelerated particles) and wave generation on the acceleration process. We show that the refraction of self-generated waves resulting from the deceleration of the plasma flow by the pressure of energetic particles causes enhanced losses of these particles. This effect slows down the acceleration and changes the shape of particle spectrum near the cut-off. The implications for observations of TeV emission from SNR remnants are also discussed.

The linear theory of shock acceleration predicts the maximum particle energy to be limited only by the acceleration time and the size of the shock. We study the combined effect of acceleration nonlinearity (shock modification by accelerated particles, that must be present in strong astrophysical shocks) and propagation of Alfven waves that are responsible for particle confinement to the shock front. We show that wave refraction to larger wave numbers in the nonlinearly modified flow causes enhanced losses of particles in the momentum range p_max/R < p < p_max, where R>1 is the nonlinear pre-compression of the flow and p_max is a conventional maximum momentum, that could be reached if there was no refraction.

We present the first result of VLBI observations of the blueshifted water maser emission from the type 2 Seyfert galaxy NGC5793, which we combine with new and previous VLBI observations of continuum emission at 1.7, 5.0, 8.4, 15, and 22 GHz. Maser emission was detected earlier in single-dish observations and found to have both red- and blueshifted features relative to the systemic velocity. We could image only the blueshifted emission, which is located 3.6 pc southwest of the 22 GHz continuum peak. The blueshifted emission was found to originate in two clusters that are separated by 0.7 milliarcsecond (0.16 pc). No compact continuum emission was found within 3.6 pc of the maser spot. A compact continuum source showing a marginally inverted spectrum between 1.7 and 5.0 GHz was found 4.2 pc southwest of the maser position. The spectral turnover might be due to synchrotron self-absorption caused by a shock in the jet owing to collision with dense gas, or it might be due to free-free absorption in an ionized screen possibly the inner part of a disk, foreground to the jet. The water maser may be part of a maser disk. If so, it would be rotating in the opposite sense to the highly inclined galactic disk observed in CO emission. We estimate a binding mass within 1 pc of the presumed nucleus to be on the order of 10^7 Msun. Alternatively, the maser emission could result from the amplification of a radio jet by foreground circumnuclear molecular gas. In this case, the high blueshift of the maser emission might mean that the masing region is moving outward away from the molecular gas surrounding an active nucleus.

We present the first circular polarization measurements of circumstellar H2O masers. Previously the magnetic field in circumstellar envelopes has been estimated using polarization observations of SiO and OH masers. SiO masers are probes of the high temperature and density regime close to the central star. OH masers are found at much lower densities and temperatures, generally much further out in the circumstellar envelope. The detection of the circular polarization of the (6_16 - 5_23) rotational transition of the H2O maser could be attributed to Zeeman splitting due to the magnetic field in the intermediate temperature and density regime. The fields inferred here agree well with predicted values for a combination of the r^{-2} dependence of a solar-type magnetic field, and the coupling of the field to the high density masing regions. We also discuss the unexpected narrowing of the circular polarization spectrum.

One of the main features of astrophysical shocks is their ability to accelerate particles to extremely high energies. The leading acceleration mechanism, the diffusive shock acceleration is reviewed. It is demonstrated that its efficiency critically depends on the injection of thermal plasma into acceleration which takes place at the subshock of the collisionless shock structure that, in turn, can be significantly smoothed by energetic particles. Furthermore, their inhomogeneous distribution provides free energy for MHD turbulence regulating the subshock strength and injection rate. Moreover, the MHD turbulence confines particles to the shock front controlling their maximum energy and bootstrapping acceleration. Therefore, the study of the MHD turbulence in a compressive plasma flow near a shock is a key to understanding of the entire process. The calculation of the injection rate became part of the collisionless shock theory. It is argued that the further progress in diffusive shock acceleration theory is impossible without a significant advance in these two areas of plasma physics.

The possibility of fast dynamo action by collisionless kinetic Alfven Wave turbulence is demonstrated. The irreversibility necessary to lock in the generated field is provided by electron Landau damping, so the induced electric field does not vanish with resistivity. Mechanisms for self-regulation of the system and the relation of these results to the theory of alpha quenching are discussed. The dynamo-generated fields have symmetry like to that of zonal flows, and thus are termed zonal fields.

We examine the dynamics of turbulent reconnection in 2D and 3D reduced MHD by calculating the effective dissipation due to coupling between small-scale fluctuations and large-scale magnetic fields. Sweet-Parker type balance relations are then used to calculate the global reconnection rate. Two approaches are employed -- quasi-linear closure and an eddy-damped fluid model. Results indicate that despite the presence of turbulence, the reconnection rate remains inversely proportional to $\sqrt{R_m}$, as in the Sweet-Parker analysis. In 2D, the global reconnection rate is shown to be enhanced over the Sweet-Parker result by a factor of magnetic Mach number. These results are the consequences of the constraint imposed on the global reconnection rate by the requirement of mean square magnetic potential balance. The incompatibility of turbulent fluid-magnetic energy equipartition and stationarity of mean square magnetic potential is demonstrated.

Observations of the starburst galaxy, M82, have been made with a 20-station global VLBI array at $\lambda$18cm. Maps are presented of the brightest young supernova remnants (SNR) in M82 and the wide-field mapping techniques used in making images over a field of view of $\sim$1 arcminute with 3 milliarcsecond resolution are discussed. A limit has been placed on the power law deceleration of the young SNR, 43.31+592 with an index greater than 0.73 $\pm$ 0.11 from observations with the European VLBI Network. Using the global array we have resolved compact knots of radio emission in the source which, with future global observations, will enable better constraints to be placed on the expansion parameters of this SNR. The latest global observations have also provided high resolution images of the most compact radio source in M82, 41.95+575. We determine an upper limit to the radial expansion rate along the major axis of 2000 km s$^{-1}$. However, the new images also show structure resembling that of collimated ejection which brings into question the previous explanation of the source as being due to the confinement of a supernova by a high density circumstellar medium. It is apparent that we are now able to image the brightest supernova remnants in M82 with a linear scale which allows direct comparison with galactic SNR such as Cassiopeia A.

The OH 1667 MHz maser in the circumstellar shell around the Mira variable U Herculis has been observed with the NRAO Very Long Baseline Array (VLBA) at 6 epochs, spread over 4 years. Using phase referencing techniques the position of the most blue-shifted maser spot was monitored with respect to two extra-galactic radio sources. The absolute radio positions of the maser can be compared with the stellar optical position measured by the Hipparcos satellite to 15 mas accuracy. This confirms the model in which one of the maser spots corresponds to the stellar continuum, amplified by the maser. The stellar proper motion and the annual parallax (5.3 +/- 2.1 mas) were measured.

There are two distinct regimes of the first order Fermi acceleration at shocks. The first is a linear (test particle) regime in which most of the shock energy goes into thermal and bulk motion of the plasma. The second is an efficient regime when it goes into accelerated particles. Although the transition region between them is narrow, we identify the factors that drive the system to a {\it self-organized critical state} between those two. Using an analytic solution, we determine this critical state and calculate the spectra and maximum energy of accelerated particles.

The water maser in the mildly active nucleus in the nearby galaxy NGC4258 traces a thin, nearly edge-on, subparsec-scale Keplerian disk. Using the technique of very long baseline interferometry, we have detected the proper motions of these masers as they sweep in front of the central black hole at an orbital velocity of about 1100 km/s. The average maser proper motion of 31.5 microarcseconds per year is used in conjunction with the observed acceleration of the masers to derive a purely geometric distance to the galaxy of 7.2 +- 0.3 Mpc. This is the most precise extragalactic distance measured to date, and, being independent of all other distance indicators, is likely to play an important role in calibrating the extragalactic distance scale.

The spatial structure, fluctuations as well as all state probabilities of self-organized (steady) states of cellular automata can be found (almost) exactly and {\em explicitly} from their Markovian dynamics. The method is shown on an example of a natural sand pile model with a gradient threshold.

The BN/KL region in the Orion molecular cloud is an archetype in the study of the formation of stars much more massive than the Sun. This region contains luminous young stars and protostars, but it is difficult to study because of overlying dust and gas. Our basic expectations are shaped to some extent by the present theoretical picture of star formation, the cornerstone of which is that protostars acrete gas from rotating equatorial disks, and shed angular momentum by ejecting gas in bipolar outflows. The main source of the outflow in the BN/KL region may be an object known as radio source I, which is commonly believed to be surrounded by a rotating disk of molecular material. Here we report high-resolution observations of silicon monoxide (SiO) and water maser emission from the gas surrounding source I; we show that within 60 AU (about the size of the Solar System), the region is dominated by a conical bipolar outflow, rather than the expected disk. A slower outflow, close to the equatorial plane of the protostellar system, extends to radii of 1,000 AU.

We show that in many methanol maser sources the masers are located in lines, with a velocity gradient along them which suggests that the masers are situated in edge-on circumstellar, or protoplanetary, disks. We present VLBI observations of the methanol maser source G309.92+0.48, in the 12.2 GHz transition, which confirm previous observations that the masers in this source lie along a line. We show that such sources are not only linear in space but, in many cases, also have a linear velocity gradient. We then model these and other data in both the 6.7 GHz and the 12.2 GHz transition from a number of star formation regions, and show that the observed spatial and velocity distribution of methanol masers, and the derived Keplerian masses, are consistent with a circumstellar disk rotating around an OB star. We consider this and other hypotheses, and conclude that about half of these methanol masers are probably located in edge-on circumstellar disks around young stars. This is of particular significance for studies of circumstellar disks because of the detailed velocity information available from the masers.

We report VLA 1 sigma upper limits of 1.5% and 3% on the intrinsic circular and linear fractional polarizations, respectively, of the water vapor maser emission 0.2 pc from the central engine of NGC 4258. A corresponding 0.5% upper limit on any Zeeman-splitting-induced circular polarization translates to a 1 sigma upper limit on the parallel, or toroidal, component of the magnetic field of 300 mG. Assuming magnetic and thermal pressure balance in the disk, this magnetic field upper limit corresponds to a model-dependent estimate of the accretion rate through the molecular disk of 10^-1.9 alpha solar masses per year for the case where the magnetic field lies along the line of sight.

A fully nonlinear, time-asymptotic theory of resonant particle trapping in large-amplitude quasi-parallel Alfven waves is presented. The effect of trapped particles on the nonlinear dynamics of quasi-stationary Alfvenic discontinuities and coherent Alfven waves is highly non-trivial and forces to a significant departure of the theory from the conventional DNLS and KNLS equation models. The virial theorem is used to determine the time-asymptotic distribution function.

We report a 3 sigma upper limit of 220 microJanskys on any 22-GHz continuum emission coincident with the central engine in NGC 4258. If NGC 4258 is powered by an advection-dominated accretion flow, this radio upper limit implies that the inner advection-dominated flow cannot extend significantly beyond about 100 Schwarzschild radii.

A rarity among supernova, SN 1993J in M81 can be studied with high spatial resolution. Its radio power and distance permit VLBI observations to monitor the expansion of its angular structure. This radio structure was previously revealed to be shell-like and to be undergoing a self-similar expansion at a constant rate. From VLBI observations at the wavelengths of 3.6 and 6 cm in the period 6 to 42 months after explosion, we have discovered that the expansion is decelerating. Our measurement of this deceleration yields estimates of the density profiles of the supernova ejecta and circumstellar material in standard supernova explosion models.

A simple model collisionless, dissipative, compressible MHD (Alfvenic) turbulence in a magnetized system is investigated. In contrast to more familiar paradigms of turbulence, dissipation arises from Landau damping, enters via nonlinearity, and is distributed over all scales. The theory predicts that two different regimes or phases of turbulence are possible, depending on the ratio of steepening to damping coefficient (m_1/m_2). For strong damping (|m_1/m_2|<1), a regime of smooth, hydrodynamic turbulence is predicted. For |m_1/m_2|>1, steady state turbulence does not exist in the hydrodynamic limit. Rather, spikey, small scale structure is predicted.

The nonlinear dynamics of collisionless Alfven trains, including resonant particle effects is studied using the kinetic nonlinear Schroedinger (KNLS) equation model. Numerical solutions of the KNLS reveal the dynamics of Alfven waves to be sensitive to the sense of polarization as well as the angle of propagation with respect to the ambient magnetic field. The combined effects of both wave nonlinearity and Landau damping result in the evolutionary formation of stationaryOA S- and arc-polarized directional and rotational discontinuities. These waveforms are freqently observed in the interplanetary plasma.

Collisionless regime kinetic models for coherent nonlinear Alfven wave dynamics are studied using fluid moment equations with an approximate closure anzatz. Resonant particle effects are modelled by incorporating an additional term representing dissipation akin to parallel heat conduction. Unlike collisional dissipation, parallel heat conduction is presented by an integral operator. The modified derivative nonlinear Schrodinger equation thus has a spatially nonlocal nonlinear term describing the long-time evolution of the envelope of parallel-propagating Alfven waves, as well. Coefficients in the nonlinear terms are free of the 1/(1-beta) singularity usually encountered in previous analyses, and have very a simple form which clarifies the physical processes governing the large amplitude Alfvenic nonlinear dynamics. The nonlinearity appears via coupling of an Alfvenic mode to a kinetic ion-acoustic mode. Damping of the nonlinear Alfven wave appears via strong Landau damping of the ion-acoustic wave when the electron-to-ion temperature ratio is close to unity. For a (slightly) obliquely propagating wave, there are finite Larmor radius corrections in the dynamical equation. This effect depends on the angle of wave propagation relative to B_0 and vanishes for the limit of strictly parallel propagation. Explicit magnetic perturbation envelope equations amenable to further analysis and numerical solution are obtained. Implications of these models for collisionless shock dynamics are discussed.

The saturation of ion cyclotron Alfven turbulence excited by beam particles is investigated using resonance broadening theory. The stochastic scattering which decorrelates particles, includes both random acceleration by electric fields and a turbulent magnetic mirroring effect. Turbulent mirroring is shown to yield non-Gaussian corrections to the orbits even if the random electric and magnetic fields are Gaussian. The predicted steady-state turbulence level exhibits a peaked anglular distribution, with a maximum near Theta ~ 60 degrees.

The influence of various kinetic effects (e.g. Landau damping, diffusive and collisional dissipation, and finite Larmor radius terms) on the nonlinear evolution of finite amplitude Alfvenic wave trains in a finite-beta environment is systematically investigated using a novel, kinetic nonlinear Schrodinger (KNLS) equation. The dynamics of Alfven waves is sensitive to the sense of polarization as well as the angle of propagation with respect to the ambient magnetic field. Numerical solution for the case with Landau damping reveals the formation of dissipative structures, which are quasi-stationary, S-polarized directional (and rotational) discontinuities which self-organize from parallel propagating, linearly polarized waves. Parallel propagating circularly polarized packets evolve to a few circularly polarized Alfven harmonics on large scales. Stationary arc-polarized rotational discontinuities form from obliquely propagating waves. Collisional dissipation, even if weak, introduces enhanced wave damping when beta is very close to unity. Cyclotron motion effects on resonant particle interactions introduce cyclotron resonance into the nonlinear Alfven wave dynamics.

We present milliarcsecond resolution images of the 6.7 and 12.2 GHz methanol maser emission associated with the well-known star formation region NGC 6334F. The images agree well with previous lower resolution observations, but detect approximately double the number of spots seen in the earlier work. Comparison of the relative positions of the 6.7 and 12.2 GHz maser spots shows that five of them are coincident to within the positional accuracy of these observations (approximately 4 milli-arcsec). Menten et al. (1992) observed similar positional coincidence for W3(OH) and in each case the flux density of the 6.7 GHz maser spot was greater than that of the 12.2 GHz methanol maser spot. However, for NGC 6334F several of the coincident maser spots have a larger flux density at 12.2 GHz than at 6.7 GHz. We also detected several 12.2 GHz methanol maser spots with no coincident 6.7 GHz emission. This implies that, although the 6.7 GHz methanol masers usually have a greater flux density than their 12.2 GHz counterparts, regions within the gas cloud exist where the conditions are more favourable for 12.2 than for 6.7 GHz methanol maser emission. If the conditions which give rise to the various observed ratios can be determined, then VLBI images of class II methanol masers will allow the physical conditions of the star formation region to be probed with unprecedented resolution.

We have made angular broadening measurements in the satellite line of the OH molecule at 1612 MHz of 8 OH/IR stars located near the Galactic Center with the Very Large Array. Early observations with the Very Long Baseline Array had revealed a region of pronounced scattering approximately centered on SgrA*. These new observations show that the compact line emission from these stars is scattered by an amount that is two to three times larger than had been previously seen. Additionally, we measure non-circular scattering disks for these objects, evidence that the electron density variations are distributed anisotropically. The present data suggest that the enhanced turbulence is located near the Galactic Center. However, we cannot completely rule out that the scattering is caused by an unrelated region along the line of sight.

We review mechanisms for the transport of angular momentum in accretion disks that are low mass, in the sense that the gravitational forces produced by the material in these disks has a negligible effect on disk dynamics. There is no established consensus on how this transport takes place. We note that for phenomenological reasons the traditional $\alpha$ model is probably not a good description of real disks. We discuss models in which angular momentum transport is driven by shocks, and by magnetic field instabilities. The latter is more promising, but requires a dynamo. We note that the direction of angular momentum transport due to convection in a conducting disk is not known as competing mechanisms are at work. We briefly discuss a number of possible dynamo mechanisms, and their problems. We then give a detailed exposition of the internal wave driven dynamo model, in which internal waves excited at large radii drive an $\alpha-\Omega$ dynamo.