The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to major uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary pathways provide a new route to understanding planetary interiors. We present the discovery of TOI-849b, the remnant core of a giant planet, with a radius smaller than Neptune but an anomalously high mass $M_p=40.8^{+2.4}_{-2.5}M_{\oplus}$ and density of $5.5\pm0.8$gcm$^{-3}$, similar to the Earth. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$% of the total mass of the planet. TOI-849b transits a late G type star (T$_{\rm mag}=11.5$) with an orbital period of 18.4 hours, leading to an equilibrium temperature of 1800K. The planet's mass is larger than the theoretical threshold mass for runaway gas accretion. As such, the planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it avoided substantial gas accretion, perhaps through gap opening or late formation. Photoevaporation rates cannot provide the mass loss required to reduce a Jupiter-like gas giant, but can remove a few $M_\oplus$ hydrogen and helium envelope on timescales of several Gyr, implying that any remaining atmosphere is likely to be enriched by water or other volatiles from the planetary interior. TOI-849b represents a unique case where material from the primordial core is left over from formation and available to study.

Interferometric observations of the mm dust distribution in protoplanetary discs are now showing a ubiquity of annular gap and ring substructures. Their identification and accurate characterization is critical to probing the physical processes responsible. We present Frankenstein (frank), an open source code that recovers axisymmetric disc structures at sub-beam resolution. By fitting the visibilities directly, the model reconstructs a disc's 1D radial brightness profile nonparametrically using a fast (<~1 min) Gaussian process. The code avoids limitations of current methods that obtain the radial brightness profile by either extracting it from the disc image via nonlinear deconvolution at the cost of reduced fit resolution, or by assumptions placed on the functional forms of disc structures to fit the visibilities parametrically. We use mock ALMA observations to quantify the method's intrinsic capability and its performance as a function of baseline-dependent signal-to-noise. Comparing the technique to profile extraction from a CLEAN image, we motivate how our fits accurately recover disc structures at a sub-beam resolution. Demonstrating the model's utility in fitting real high and moderate resolution observations, we conclude by proposing applications to address open questions on protoplanetary disc structure and processes.

When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $\alpha/St$ and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission lines data-cubes. As a proof-of-concept, we then apply the method to two discs with prominent sub-structure, HD163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good ($\alpha/St \sim 0.1$). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the $\alpha$ turbulent parameter ($\alpha \sim 10^{-2}$). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.

Exoplanets orbiting M dwarfs present a valuable opportunity for their detection and atmospheric characterisation. This is evident from recent inferences of H$_2$O in such atmospheres, including that of the habitable-zone exoplanet K2-18b. With a bulk density between Earth and Neptune, K2-18b may be expected to possess a H/He envelope. However, the extent of such an envelope and the thermodynamic conditions of the interior remain unexplored. In the present work, we investigate the atmospheric and interior properties of K2-18b based on its bulk properties and its atmospheric transmission spectrum. We constrain the atmosphere to be H$_2$-rich with a H$_2$O volume mixing ratio of $0.02-14.8$%, consistent with previous studies, and find a depletion of CH$_4$ and NH$_3$, indicating chemical disequilibrium. We do not conclusively detect clouds/hazes in the observable atmosphere. We use the bulk parameters and retrieved atmospheric properties to constrain the internal structure and thermodynamic conditions in the planet. The constraints on the interior allow multiple scenarios between rocky worlds with massive H/He envelopes and water worlds with thin envelopes. We constrain the mass fraction of the H/He envelope to be $\lesssim 6$%; spanning $\lesssim 10^{-5}$ for a predominantly water world to $\sim6$% for a pure iron interior. The thermodynamic conditions at the surface of the H$_2$O layer range from the super-critical to liquid phases, with a range of solutions allowing for habitable conditions on K2-18b. Our results demonstrate that the potential for habitable conditions is not necessarily restricted to Earth-like rocky exoplanets.

The Sun shows a $\sim 10$% depletion in refractory elements relative to nearby solar twins. It has been suggested that this depletion is a signpost of planet formation. The exoplanet statistics are now good enough to show that the origin of this depletion does not arise from the sequestration of refractory material inside the planets themselves. This conclusion arises because most sun-like stars host close-in planetary systems that are on average more massive than the Sun's. Using evolutionary models for the protoplanetary discs that surrounded the young Sun and solar twins we demonstrate that the origin of the depletion likely arises due to the trapping of dust exterior to the orbit of a forming giant planet. In this scenario a forming giant planet opens a gap in the gas disc, creating a pressure trap. If the planet forms early enough, while the disc is still massive, the planet can trap $\gtrsim 100 $M$_\oplus$ of dust exterior to its orbit, preventing the dust from accreting onto the star in contrast to the gas. Forming giant planets can create refractory depletions of $\sim 5-15$%, with the larger values occurring for initial conditions that favour giant planet formation (e.g. more massive discs, that live longer). The incidence of solar-twins that show refractory depletion matches both the occurrence of giant planets discovered in exoplanet surveys and "transition" discs that show similar depletion patterns in the material that is accreting onto the star.

Protoplanetary discs form and evolve in a wide variety of stellar environments and are accordingly exposed to a wide range of ambient far ultraviolet (FUV) field strengths. Strong FUV fields are known to drive vigorous gaseous flows from the outer disc. In this paper we conduct the first systematic exploration of the evolution of the solid component of discs subject to external photoevaporation. We find that the main effect of photoevaporation is to reduce the reservoir of dust at large radii and this leads to more efficient subsequent depletion of the disc dust due to radial drift. Efficient radial drift means that photoevaporation causes no significant increase of the dust to gas ratio in the disc. We show that the disc lifetime in both dust and gas is strongly dependent on the level of the FUV background and that the relationship between these two lifetimes just depends on the Shakura-Sunyaev $\alpha$ parameter, with the similar lifetimes observed for gas and dust in discs pointing to higher $\alpha$ values ($\sim 10^{-2}$). On the other hand the distribution of observed discs in the plane of disc size versus flux at $850~\mu$m is better reproduced by lower $\alpha$ ($\sim 10^{-3}$). We find that photoevaporation does not assist rocky planet formation but need not inhibit mechanisms (such as pebble accretion at the water snow line) which can be effective sufficiently early in the disc's lifetime (i.e. well within a Myr).

Collisional fragmentation is a ubiquitous phenomenon arising in a variety of astrophysical systems, from asteroid belts to debris and protoplanetary disks. Numerical studies of fragmentation typically rely on discretizing the size distribution of colliding objects into a large number N of bins in mass space, usually logarithmically spaced. A standard approach for redistributing the debris produced in collisions into the corresponding mass bins results in O(N^3) calculation, which leads to significant computational overhead when N is large. Here we formulate a more efficient explicit O(N^2) fragmentation algorithm, which works when the size spectrum of fragments produced in an individual collision has a self-similar shape with only a single characteristic mass scale (which can have arbitrary dependence on the energy and masses of colliding objects). Fragment size spectra used in existing fragmentation codes typically possess this property. We also show that our O(N^2) approach can be easily extended to work with non-self-similar fragment size distributions, for which we provide a worked example. This algorithm offers a substantial speedup of fragmentation calculations for large N> 100, even over the implicit methods, making it an attractive tool for studying collisionally evolving systems.

Binary stars often move through an ambient medium from which they accrete material and angular momentum, as in triple-star systems, star-forming clouds, young globular clusters and in the centres of galaxies. A binary form of Bondi-Hoyle-Lyttleton accretion results whereby the accretion rate depends on the binary properties: the stellar masses and separation, and the relative wind speed. We present the results of simulations performed with the hydrodynamic code GANDALF, to determine the mass accretion rates over a range of binary separations, inclinations and mass ratios. When the binary separation is short, the binary system accretes like a single star, while accretion onto stars in wide binaries is barely affected by their companion. We investigate intermediate-separation systems in some detail, finding that as the binary separation is increased, accretion rates smoothly decrease from the rate equal to that of a single star to the rate expected from two isolated stars. The form of this decrease depends on the relative centre-of-mass velocity of the binary and the gas, with faster-moving binaries showing a shallower decrease. Accretion rates vary little with orbital inclination, except when the orbit is side-on and the stars pass through each others' wakes. The specific angular momentum accretion rate also depends on the inclination but is never sufficient to prevent the binary orbit from contracting. Our results may be applied to accretion onto protostars, pollution of stars in globular and nuclear clusters, and wind mass-transfer in multiple stellar systems.

Stellar magnetic activity decays over the main-sequence life of cool stars due to the stellar spin-down driven by magnetic braking. The evolution of chromospheric emission is well-studied for younger stars, but difficulties in determining the ages of older cool stars on the main sequence have complicated such studies for older stars in the past. Here we report on chromospheric Ca II H and K line measurements for 26 main-sequence cool stars with asteroseismic ages older than a gigayear and spectral types F and G. We find that for the G stars and the cooler F-type stars which still have convective envelopes the magnetic activity continues to decrease at stellar ages above one gigayear. Our magnetic activity measurements do not show evidence for a stalling of the magnetic braking mechanism, which has been reported for stellar rotation versus age for G and F type stars. We also find that the measured R'HK indicator value for the cool F stars in our sample is lower than predicted by common age-activity relations that are mainly calibrated on data from young stellar clusters. We conclude that, within individual spectral type bins, chromospheric magnetic activity correlates well with stellar age even for old stars.

The Galactic Centre contains a supermassive black hole with a mass of 4 million suns within an environment that differs markedly from that of the Galactic disk. While the black hole is essentially quiescent in the broader context of active galactic nuclei, X-ray observations have provided evidence for energetic outbursts from its surroundings. Also, while the levels of star formation in the Galactic Centre have been approximately constant over the last few hundred Myr, there is evidence of elevated short-duration bursts, strongly influenced by interaction of the black hole with the enhanced gas density present within the ring-like Central Molecular Zone at Galactic longitude |l| < 0.7 degrees and latitude |b| < 0.2 degrees. The inner 200 pc region is characterized by large amounts of warm molecular gas, a high cosmic ray ionization rate, unusual gas chemistry, enhanced synchrotron emission, and a multitude of radio-emitting magnetised filaments, the origin of which has not been established. Here we report radio imaging that reveals bipolar bubbles spanning 1 degree x 3 degrees (140 parsecs x 430 parsecs), extending above and below the Galactic plane and apparently associated with the Galactic Centre. The structure is edge-brightened and bounded, with symmetry implying creation by an energetic event in the Galactic Centre. We estimate the age of the bubbles to be a few million years, with a total energy of 7 x 10^52 ergs. We postulate that the progenitor event was a major contributor to the increased cosmic-ray density in the Galactic Centre, and is in turn the principal source of the relativistic particles required to power the synchrotron emission of the radio filaments within and in the vicinity of the bubble cavities.

The presence of planets in binary systems poses interesting problems for planet formation theories, both in cases where planets must have formed in very compact discs around the individual stars and where they are located near the edge of the stable circumbinary region, where in situ formation is challenging. Dust dynamics is expected to play an important role in such systems, since dust trapping at the inner edge of circumbinary discs could aid in situ formation, but would simultaneously starve the circumstellar discs of the solid material needed to form planets. Here we investigate the dynamics of dust in binary systems using Smooth Particle Hydrodynamics. We find that all our simulations tend towards dust trapping in the circumbinary disc, but the timescale on which trapping begins depends on binary mass ratio ($q$) and eccentricity as well as the angular momentum of the infalling material. For $q \gtrsim 0.1$, we find that dust can initially accrete onto the circumstellar discs, but as the circumbinary cavity grows in radius, dust eventually becomes trapped in the circumbinary disc. For $q = 0.01$, we find that increasing the binary eccentricity increases the time required for dust trapping to begin. However, even this longer timescale is likely to be shorter than the planet formation timescale in the inner disc and is insufficient to explain the observed pre-transitional discs. This indicates that increase in companion eccentricity alone is not enough to allow significant transfer of solids from the outer to the inner disc.

The composition of gas and solids in protoplanetary discs sets the composition of planets that form out of them. Recent chemical models have shown that the composition of gas and dust in discs evolves on Myr time-scales, with volatile species disappearing from the gas phase. However, discs evolve due to gas accretion and radial drift of dust on time-scales similar to these chemical time-scales. Here we present the first model coupling the chemical evolution in the disc mid-planes with the evolution of discs due to accretion and radial drift of dust. Our models show that transport will always overcome the depletion of CO$_2$ from the gas phase, and can also overcome the depletion of CO and CH$_4$ unless both transport is slow (viscous $\alpha \lesssim 10^{-3}$) and the ionization rate is high ($\zeta \approx 10^{-17}$). Including radial drift further enhances the abundances of volatile species because they are carried in on the surface of grains before evaporating left at their ice lines. Due to large differences in the abundances within 10 au for models with and without efficient radial drift, we argue that composition can be used to constrain models of planet formation via pebble accretion.

A correlation between proto-planetary disc radii and sub-mm fluxes has been recently reported. In this Letter we show that the correlation is a sensitive probe of grain growth processes. Using models of grain growth and drift, we have shown in a companion paper that the observed disc radii trace where the dust grains are large enough to have a significant sub-mm opacity. We show that the observed correlation emerges naturally if the maximum grain size is set by radial drift, implying relatively low values of the viscous $\alpha$ parameter $ \lesssim 0.001$. In this case the relation has an almost universal normalisation, while if the grain size is set by fragmentation the flux at a given radius depends on the dust-to-gas ratio. We highlight two observational consequences of the fact that radial drift limits the grain size. The first is that the dust masses measured from the sub-mm could be overestimated by a factor of a few. The second is that the correlation should be present also at longer wavelengths (e.g. 3mm), with a normalisation factor that scales as the square of the observing frequency as in the optically thick case.

Proto-planetary disc surveys conducted with ALMA are measuring disc radii in multiple star forming regions. The disc radius is a fundamental quantity to diagnose whether discs undergo viscous spreading, discriminating between viscosity or angular momentum removal by winds as drivers of disc evolution. Observationally, however, the sub-mm continuum emission is dominated by the dust, which also drifts inwards, complicating the picture. In this paper we investigate, using theoretical models of dust grain growth and radial drift, how the radii of dusty viscous proto-planetary discs evolve with time. Despite the existence of a sharp outer edge in the dust distribution, we find that the radius enclosing most of the dust $\textit{mass}$ increases with time, closely following the evolution of the gas radius. This behaviour arises because, although dust initially grows and drifts rapidly onto the star, the residual dust retained on Myr timescales is relatively well coupled to the gas. Observing the expansion of the dust disc requires using definitions based on high fractions of the disc $\textit{flux}$ (e.g. 95 per cent) and very long integrations with ALMA, because the dust grains in the outer part of the disc are small and have a low sub-mm opacity. We show that existing surveys lack the sensitivity to detect viscous spreading. The disc radii they measure do not trace the mass radius or the sharp outer edge in the dust distribution, but the outer limit of where the grains have significant sub-mm opacity. We predict that these observed radii should shrink with time.

The Galactic magnetic field is an integral constituent of the interstellar medium (ISM), and knowledge of its structure is crucial to understanding Galactic dynamics. The Rotation Measures (RM) of extragalactic (EG) sources have been the basis of comprehensive Galactic magnetic field models. Polarised extended emission (XE) is also seen along lines of sight through the Galactic disk, and also displays the effects of Faraday rotation. Our aim is to investigate and understand the relationship between EG and XE RMs near the Galactic plane, and to determine how the XE RMs, a hitherto unused resource, can be used as a probe of the large-scale Galactic magnetic field. We used polarisation data from the Canadian Galactic Plane Survey (CGPS), observed near 1420 MHz with the Dominion Radio Astrophysical Observatory (DRAO) Synthesis Telescope. We calculated RMs from a linear fit to the polarisation angles as a function of wavelength squared in four frequency channels, for both the EG sources and the XE. Across the CGPS area, $55^{\circ} < {\ell} <193^{\circ}, -3^{\circ} < b < 5^{\circ}$, the RMs of the XE closely track the RMs of the EG sources, with XE RMs about half the value of EG-source RMs. The exceptions are places where large local HII complexes heavily depolarise more distant emission. We conclude that there is valuable information in the XE RM dataset. The factor of 2 between the two types of RM values is close to that expected from a Burn slab model of the ISM. This result indicates that, at least in the outer Galaxy, the EG and XE sources are likely probing similar depths, and that the Faraday rotating medium and the synchrotron emitting medium have similar variation with galactocentric distance.

Recent observations show that rings and gaps are ubiquitous in protoplanetary discs. These features are often interpreted as being due to the presence of planets; however, the effect of planetary migration on the observed morphology has not been investigated hitherto. In this work we investigate whether multiwavelength mm/submm observations can detect signatures of planet migration, using 2D dusty hydrodynamic simulations to model the structures generated by migrating planets and synthesising ALMA continuum observations at 0.85 and 3 mm. We identify three possible morphologies for a migrating planet: a slowly migrating planet is associated with a single ring outside the planet's orbit, a rapidly migrating planet is associated with a single ring inside the planet's orbit while a planet migrating at intermediate speed generates one ring on each side of the planet's orbit. We argue that multiwavelength data can distinguish multiple rings produced by a migrating planet from other scenarios for creating multiple rings, such as multiple planets or discs with low viscosity. The signature of migration is that the outer ring has a lower spectral index, due to larger dust grains being trapped there. Of the recent ALMA observations revealing protoplanetary discs with multiple rings and gaps, we suggest that Elias 24 is the best candidate for a planet migrating in the intermediate speed regime.

We have investigated the properties of gravito-turbulent discs in 3D using high-resolution shearing-box simulations. For large enough domain sizes, $L_y \gtrsim 60H$, the disc settles down into a quasi-steady state, showing no long term trends in properties or variation with box size. For smaller boxes, we find that the azimuthal wavelength of the dominant spiral modes are limited to the domain size. This is associated with a bursty behaviour that differs from the quasi-steady dynamics at larger sizes. We point out that a similar transition may be expected in global simulations at the point where the range of azimuthal wavelengths is limited by the finite disc size. This condition (i.e. when $60 H \sim 2 \pi R$, i.e. $H/R \sim 0.1$) correctly predicts the transition to bursty behaviour previously found in global simulations for disc-to-star mass ratios in excess of 0.25. We recover a transition in the dynamics from two- to three-dimensional behaviour, characterized by a turbulence that becomes more isotropic on small scales. This turbulence likely plays an important role in the evolution of dust in self-gravitating discs, potentially dominating the collision velocity for particles with Stokes number $< 1$. Finally, we consider the stability of gravito-turbulence against fragmentation, finding that discs which cool faster than a few dynamical times fragment immediately, supporting previous results. We also find hints of stochastic fragmentation at marginally longer cooling times, in which a fragment forms from a quasi-equilibrium state. However, this makes little practical difference to region where planet formation via gravitational instability may occur.

The formation of planets occurs within protoplanetary disks surrounding young stars, resulting in perturbation of the gas and dust surface densities. Here, we report the first evidence of spatially resolved gas surface density ($\Sigma_{g}$) perturbation towards the AS~209 protoplanetary disk from the optically thin C$^{18}$O ($J=2-1$) emission. The observations were carried out at 1.3~mm with ALMA at a spatial resolution of about 0.3$\arcsec$ $\times$ 0.2$\arcsec$ (corresponding to $\sim$ 38 $\times$ 25 au). The C$^{18}$O emission shows a compact ($\le$60~au), centrally peaked emission and an outer ring peaking at 140~au, consistent with that observed in the continuum emission and, its azimuthally averaged radial intensity profile presents a deficit that is spatially coincident with the previously reported dust map. This deficit can only be reproduced with our physico-thermochemical disk model by lowering $\Sigma_{gas}$ by nearly an order of magnitude in the dust gaps. Another salient result is that contrary to C$^{18}$O, the DCO$^{+}$ ($J=3-2$) emission peaks between the two dust gaps. We infer that the best scenario to explain our observations (C$^{18}$O deficit and DCO$^{+}$ enhancement) is a gas perturbation due to forming-planet(s), that is commensurate with previous continuum observations of the source along with hydrodynamical simulations. Our findings confirm that the previously observed dust gaps are very likely due to perturbation of the gas surface density that is induced by a planet of at least 0.2~M$\rm_{Jupiter}$ in formation. Finally, our observations also show the potential of using CO isotopologues to probe the presence of saturn mass planet(s).

Planet migration in protoplanetary discs plays an important role in the longer term evolution of planetary systems, yet we currently have no direct observational test to determine if a planet is migrating in its gaseous disc. We explore the formation and evolution of dust rings - now commonly observed in protoplanetary discs by ALMA - in the presence of relatively low mass (12-60 Earth masses) migrating planets. Through two dimensional hydrodynamical simulations using gas and dust we find that the importance of perturbations in the pressure profile interior and exterior to the planet varies for different particle sizes. For small sizes a dust enhancement occurs interior to the planet, whereas it is exterior to it for large particles. The transition between these two behaviours happens when the dust drift velocity is comparable to the planet migration velocity. We predict that an observational signature of a migrating planet consists of a significant outwards shift of an observed midplane dust ring as the wavelength is increased.

We present high resolution millimeter continuum imaging of the disc surrounding the young star CI Tau, a system hosting the first hot Jupiter candidate in a protoplanetary disc system. The system has extended mm emission on which are superposed three prominent annular gaps at radii ~ 13, 39 and 100 au. We argue that these gaps are most likely to be generated by massive planets so that, including the hot Jupiter, the system contains four gas giant planets at an age of only 2 Myr. Two of the new planets are similarly located to those inferred in the famous HL Tau protoplanetary disc; in CI Tau, additional observational data enables a more complete analysis of the system properties than was possible for HL Tau. Our dust and gas dynamical modeling satisfies every available observational constraint and points to the most massive ensemble of exo-planets ever detected at this age, with its four planets spanning a factor 1000 in orbital radius.Our results show that the association between hot Jupiters and gas giants on wider orbits, observed in older stars, is apparently in place at an early evolutionary stage.

Theory and observations suggest that star formation occurs hierarchically due to the fragmentation of giant molecular clouds. In this case we would expect substructure and enhanced stellar multiplicity in the primordial cluster. This substructure is expected to decay quickly in most environments, however historic stellar encounters might leave imprints in a protoplanetary disc (PPD) population. In a low density environment such as Taurus, tidal tails from violent star-disc or disc-disc encounters might be preserved over time-scales sufficient to be observed. In this work, we investigate the possibility that just such an event occured between HV Tau C (itself a component of a triple system) and DO Tau $\sim 0.1$ Myr ago, as evidenced by an apparent `bridge' structure evident in the $160$ $\mu$m emission. By modelling the encounter using smoothed particle hydrodynamics (SPH) we reproduce the main features of the observed extended structure (`V'-shaped emission pointing west of HV Tau and a tail-like structure extending east of DO Tau). We suggest that HV Tau and DO Tau formed together in a quadruple system on a scale of $\sim 5000$ au ($0.025$ pc).

We present high-resolution ALMA Band 6 and 7 observations of the tidally disrupted protoplanetary disks of the RW Aurigae binary. Our observations reveal the presence of additional tidal streams to the previously observed tidal arm around RW Aur A. The observed configuration of tidal streams surrounding RW Aur A and B is incompatible with a single star--disk tidal encounter, suggesting that the RW Aurigae system has undergone multiple fly-by interactions. We also resolve the circumstellar disks around RW Aur A and B, with CO radii of 58 au and 38 au consistent with tidal truncation, and 2.5 times smaller dust emission radii. The disks appear misaligned by 12$^{\circ}$ or 57$^{\circ}$. Using new photometric observations from the American Association of Variable Star Observers (AAVSO) and All Sky Automated Survey for SuperNovae (ASAS-SN) archives, we have also identified an additional dimming event of the primary that began in late 2017 and is currently ongoing. With over a century of photometric observations, we are beginning to explore the same spatial scales as ALMA.

New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation.

Protoplanetary discs around brown dwarfs and very low mass stars offer some of the best prospects for forming Earth-sized planets in their habitable zones. To this end, we study the nature of the disc around the very low mass star V410 X-ray 1, whose SED is indicative of an optically thick and very truncated dust disc, with our modelling suggesting an outer radius of only 0.6 au. We investigate two scenarios that could lead to such a truncation, and find that the observed SED is compatible with both. The first scenario involves the truncation of both the dust and gas in the disc, perhaps due to a previous dynamical interaction or the presence of an undetected companion. The second scenario involves the fact that a radial location of 0.6 au is close to the expected location of the H$_2$O snowline in the disc. As such, a combination of efficient dust growth, radial migration, and subsequent fragmentation within the snowline leads to an optically thick inner dust disc and larger, optically thin outer dust disc. We find that a firm measurement of the CO $J=2$--1 line flux would enable us to distinguish between these two scenarios, by enabling a measurement of the radial extent of gas in the disc. Many models we consider contain at least several Earth-masses of dust interior to 0.6 au, suggesting that V410 X-ray 1 could be a precursor to a system with tightly-packed inner planets, such as TRAPPIST-1.

The majority of stars form in a clustered environment. This has an impact on the evolution of surrounding protoplanetary discs (PPDs) due to either photoevaporation or tidal truncation. Consequently, the development of planets depends on formation environment. Here we present the first thorough investigation of tidally induced angular momentum loss in PPDs in the distant regime, partly motivated by claims in the literature for the importance of distant encounters in disc evolution. We employ both theoretical predictions and dynamical/hydrodynamical simulations in 2D and 3D. Our theoretical analysis is based on that of Ostriker (1994) and leads us to conclude that in the limit that the closest approach distance $x_{min} \gg r$, the radius of a particle ring, the fractional change in angular momentum scales as $(x_{min}/r)^{-5}$. This asymptotic limit ensures that the cumulative effect of distant encounters is minor in terms of its influence on disc evolution. The angular momentum transfer is dominated by the $m=2$ Lindblad resonance for closer encounters and by the $m=1$, $\omega = 0$ Lindblad resonance at large $x_{min}/r$. We contextualise these results by comparing expected angular momentum loss for the outer edge of a PPD due to distant and close encounters. Contrary to the suggestions of previous works we do not find that distant encounters contribute significantly to angular momentum loss in PPDs. We define an upper limit for closest approach distance where interactions are significant as a function of arbitrary host to perturber mass ratio $M_2/M_1$.

For grain growth to proceed effectively and lead to planet formation a number of barriers to growth must be overcome. One such barrier, relevant for compact grains in the inner regions of the disc, is the `bouncing barrier' in which large grains ($\sim$ mm size) tend to bounce off each other rather than sticking. However, by maintaining a population of small grains it has been suggested that cm-size particles may grow rapidly by sweeping up these small grains. We present the first numerically resolved investigation into the conditions under which grains may be lucky enough to grow beyond the bouncing barrier by a series of rare collisions leading to growth (so-called `breakthrough'). Our models support previous results, and show that in simple models breakthrough requires the mass ratio at which high velocity collisions transition to growth instead of causing fragmentation to be low, $\phi \lesssim 50$. However, in models that take into account the dependence of the fragmentation threshold on mass-ratio, we find breakthrough occurs more readily, even if mass transfer is relatively inefficient. This suggests that bouncing may only slow down growth, rather than preventing growth beyond a threshold barrier. However, even when growth beyond the bouncing barrier is possible, radial drift will usually prevent growth to arbitrarily large sizes.

During the process of planet formation, the planet-discs interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long ($t\sim 3\times 10^5$ orbits) numerical simulations: (a) one (with a relatively light disc, $M_{\rm d}/M_{\rm p}=0.2$) where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, $M_{\rm d}/M_{\rm p}=0.65$) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.

The paper presents new high angular resolution ALMA 1.3 mm dust continuum observations of the protoplanetary system AS 209 in the Ophiuchus star forming region. The dust continuum emission is characterized by a main central core and two prominent rings at $r = 75\,$au and $r = 130\,$au intervaled by two gaps at at $r = 62\,$au and $r = 103\,$au. The two gaps have different widths and depths, with the inner one being narrower and shallower. We determined the surface density of the millimeter dust grains using the 3D radiative transfer disk code \textsc{dali}. According to our fiducial model the inner gap is partially filled with millimeter grains while the outer gap is largely devoid of dust. The inferred surface density is compared to 3D hydrodynamical simulations (FARGO-3D) of planet-disk interaction. The outer dust gap is consistent with the presence of a giant planet ($M_{\rm planet} \sim 0.8\,M_{\rm Staturn}$); the planet is responsible for the gap opening and for the pile-up of dust at the outer edge of the planet orbit. The simulations also show that the same planet can give origin to the inner gap at $r = 62\,$au. The relative position of the two dust gaps is close to the 2:1 resonance and we have investigated the possibility of a second planet inside the inner gap. The resulting surface density (including location, width and depth of the two dust gaps) are in agreement with the observations. The properties of the inner gap pose a strong constraint to the mass of the inner planet ($M_{\rm planet} < 0.1\,M_{\rm J}$). In both scenarios (single or pair of planets), the hydrodynamical simulations suggest a very low disk viscosity ($\alpha < 10^{-4}$). Given the young age of the system (0.5 - 1 Myr), this result implies that the formation of giant planets occurs on a timescale of $\lesssim$ 1\,Myr.

We present the discovery of NGTS-1b, a hot-Jupiter transiting an early M-dwarf host ($T_{eff}=3916^{+71}_{-63}~K$) in a P=2.674d orbit discovered as part of the Next Generation Transit Survey (NGTS). The planet has a mass of $0.812^{+0.066}_{-0.075}~M_{J}$, making it the most massive planet ever discovered transiting an M-dwarf. The radius of the planet is $1.33^{+0.61}_{-0.33}~R_{J}$. Since the transit is grazing, we determine this radius by modelling the data and placing a prior on the density from the population of known gas giant planets. NGTS-1b is the third transiting giant planet found around an M-dwarf, reinforcing the notion that close-in gas giants can form and migrate similar to the known population of hot Jupiters around solar type stars. The host star shows no signs of activity, and the kinematics hint at the star being from the thick disk population. With a deep (2.5%) transit around a $K=11.9$ host, NGTS-1b will be a strong candidate to probe giant planet composition around M-dwarfs via JWST transmission spectroscopy.

MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities.

GANDALF is a new hydrodynamics and N-body dynamics code designed for investigating planet formation, star formation and star cluster problems. GANDALF is written in C++, parallelised with both OpenMP and MPI and contains a python library for analysis and visualisation. The code has been written with a fully object-oriented approach to easily allow user-defined implementations of physics modules or other algorithms. The code currently contains implementations of Smoothed Particle Hydrodynamics, Meshless Finite-Volume and collisional N-body schemes, but can easily be adapted to include additional particle schemes. We present in this paper the details of its implementation, results from the test suite, serial and parallel performance results and discuss the planned future development. The code is freely available as an open source project on the code-hosting website github at https://github.com/gandalfcode/gandalf and is available under the GPLv2 license.

We study the sub-Keplerian rotation and dust content of the circumstellar material around the asymptotic giant branch (AGB) star L$_2$ Puppis. We find that the thermal pressure gradient alone cannot explain the observed rotation profile. We find that there is a family of possible dust populations for which radiation pressure can drive the observed sub-Keplerian rotation. This set of solutions is further constrained by the spectral energy distribution (SED) of the system, and we find that a dust-to-gas mass ratio of $\sim10^{-3}$ and a maximum grain size that decreases radially outwards can satisfy both the rotation curve and SED. These dust populations are dynamically tightly coupled to the gas azimuthally. However grains larger than $\sim0.5\mu$m are driven outward radially by radiation pressure at velocities $\sim5$kms$^{-1}$, which implies a dust replenishment rate of $\sim3\times10^{-9}$M$_\odot$yr$^{-1}$. This replenishment rate is consistent with observational estimates to within uncertainties. Coupling between the radial motion of the dust and gas is weak and hence the gas does not share in this rapid outward motion. Overall we conclude that radiation pressure is a capable and necessary mechanism to explain the observed rotation profile of L$_2$ Pup, and offers other additional constraints on the dust properties.

The chemical conditions in the planet forming regions of protoplanetary discs remain difficult to observe directly. Gas accreting from the disc on to the star provides a way to measure the elemental abundances because even refractory species are in an atomic gaseous form. Here we compare the abundance ratios derived from UV lines probing T Tauri accretion streams to simple models of disc evolution. Although the interpretation of line ratios in terms of abundances is highly uncertain, discs with large cavities in mm images tend to have lower Si emission. Since this can naturally be explained by the suppressed accretion of dust, this suggests that abundance variations are at least partially responsible for the variations seen in the line ratios. Our models of disc evolution due to grain growth, radial drift and the flux of volatile species carried as ices on grain surfaces, give rise to a partial sorting of the atomic species based on the volatility of their dominant molecular carriers. This arises because volatiles are left behind at their snow lines while the grains continue to drift. We show that this reproduces the main features seen in the accretion line ratio data, such as C/N ratios that are a few times solar and the correlation between the Si to volatile ratio with mm-flux. We highlight the fact that developing a more robust linkage between line ratios and abundance ratios and acquiring data for larger samples has the potential to cast considerable light on the chemical history of protoplanetary discs.

In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.

Stars with convective envelopes display magnetic activity, which decreases over time due to the magnetic braking of the star. This age-dependence of magnetic activity is well-studied for younger stars, but the nature of this dependence for older stars is not well understood. This is mainly because absolute stellar ages for older stars are hard to measure. However, relatively accurate stellar ages have recently come into reach through asteroseismology. In this work we present X-ray luminosities, which are a measure for magnetic activity displayed by the stellar coronae, for 24 stars with well-determined ages older than a gigayear. We find 14 stars with detectable X-ray luminosities and use these to calibrate the age-activity relationship. We find a relationship between stellar X-ray luminosity, normalized by stellar surface area, and age that is steeper than the relationships found for younger stars, with an exponent of $-2.80 \pm 0.72$. Previous studies have found values for the exponent of the age-activity relationship ranging between -1.09 to -1.40, dependent on spectral type, for younger stars. Given that there are recent reports of a flattening relationship between age and rotational period for old cool stars, one possible explanation is that we witness a strong steepening of the relationship between activity and rotation.

Chemical compositions of giant planets provide a means to constrain how and where they form. Traditionally, super-stellar elemental abundances in giant planets were thought to be possible due to accretion of metal-rich solids. Such enrichments are accompanied by oxygen-rich compositions (i.e. C/O below the disc's value, assumed to be solar, C/O=0.54). Without solid accretion the planets are expected to have sub-solar metallicity, but high C/O ratios. This arises because the solids are dominated by oxygen-rich species, e.g. H$_2$O and CO$_2$, which freeze out in the disk earlier than CO, leaving the gas metal poor but carbon-rich. Here we demonstrate that super-solar metallicities can be achieved by gas accretion alone when growth and radial drift of pebbles are considered in protoplanetary discs. Through this mechanism planets may simultaneously acquire super-solar metallicities and super-solar C/O ratios. This happens because the pebbles transport volatile species inward as they migrate through the disc, enriching the gas at snow lines where the volatiles sublimate. Furthermore, the planet's composition can be used to constrain where it formed. Since high C/H and C/O ratios cannot be created by accreting solids, it may be possible to distinguish between formation via pebble accretion and planetesimal accretion by the level of solid enrichment. Finally, we expect that Jupiter's C/O ratio should be near or above solar if its enhanced carbon abundance came through accreting metal rich gas. Thus Juno's measurement of Jupiter's C/O ratio should determine whether Jupiter accreted its metals from carbon rich gas or oxygen rich solids.

We examine the gender balance of the 18th and 19th meetings of the Cambridge Workshop on Cool Stellar Systems and the Sun (CS18 and CS19). The percent of female attendees at both meetings (31% at CS18 and 37% at CS19) was higher than the percent of women in the American Astronomical Society (25%) and the International Astronomical Union (18%). The representation of women in Cool Stars as SOC members, invited speakers, and contributed speakers was similar to or exceeded the percent of women attending the meetings. We requested that conference attendees assist in a project to collect data on the gender of astronomers asking questions after talks. Using this data, we found that men were over-represented (and women were under-represented) in the question sessions after each talk. Men asked 79% of the questions at CS18 and 75% of the questions at CS19, but were 69% and 63% of the attendees respectively. Contrary to findings from previous conferences, we did not find that the gender balance of questions was strongly affected by the session chair gender, the speaker gender, or the length of the question period. We also found that female and male speakers were asked a comparable number of questions after each talk. The contrast of these results from previous incarnations of the gender questions survey indicate that more data would be useful in understanding the factors that contribute to the gender balance of question askers. We include a preliminary set of recommendations based on this and other work on related topics, but also advocate for additional research on the demographics of conference participants. Additional data on the intersection of gender with race, seniority, sexual orientation, ability and other marginalized identities is necessary to fully address the role of gender in asking questions at conferences.

The young star Elias 2-27 has recently been observed to posses a massive circumstellar disc with two prominent large-scale spiral arms. In this Letter we perform three-dimensional Smoothed Particle Hydrodynamics simulations, radiative transfer modelling, synthetic ALMA imaging and an unsharped masking technique to explore three possibilities for the origin of the observed structures -- an undetected companion either internal or external to the spirals, and a self-gravitating disc. We find that a gravitationally unstable disc and a disc with an external companion can produce morphology that is consistent with the observations. In addition, for the latter, we find that the companion could be a relatively massive planetary mass companion (less than approximately 10 - 13 MJup) and located at large radial distances (between approximately 300 - 700 au). We therefore suggest that Elias 2-27 may be one of the first detections of a disc undergoing gravitational instabilities, or a disc that has recently undergone fragmentation to produce a massive companion.

RW Aur is a young binary system showing strong signatures of a recent tidal encounter between the circumprimary disk and the secondary star. The primary star has recently undergone two major dimming events ($\Delta$mag $\approx$ 2 in V-band), whose origin is still under debate. To shed light on the mechanism leading to the dimming events, we study the extinction properties, accretion variability, and gas kinematics using absorption lines from the material obscuring star RW Aur A. We compare our moderate resolution X-Shooter spectra of the dim state of RW Aur A with other spectral observations. In particular, we analyse archival high resolution UVES spectra obtained during the bright state of the system, in order to track the evolution of the spectral properties across the second dimming event. The spectrum obtained during the dim state shows narrow absorption lines in the Na and K optical doublets, where the former is saturated. With a velocity of -60 km/s these lines indicate that during the dim state the disk wind is either enhanced, or significantly displaced into the line of sight. The photometric evolution across the dimming event shows a gray extinction, and is correlated with a significant reduction of the EW of all photospheric lines. Emission lines tracing accretion do not vary significantly across the dimming. We conclude that the dimming event is related to a major perturbation on the inner disk. We suggest that the inner disk is occulting (most of) the star, and thus its photosphere, but is not occulting the accretion regions within a few stellar radii. Since observations of the outer disk indicate that the disk is modestly inclined (45 - 60 deg), we propose that the inner disk might be warped by a yet unseen (sub-)stellar companion, which may also explain the 2.77 day periodic variability of the spectral lines.

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

Following the recent discovery of the first radial velocity planet in a star still possessing a protoplanetary disc (CI Tau), we examine the origin of the planet's eccentricity (e $\sim 0.3$). We show through long timescale ($10^5$ orbits) simulations that the planetary eccentricity can be pumped by the disc, even when its local surface density is well below the threshold previously derived from short timescale integrations. We show that the disc may be able to excite the planet's orbital eccentricity in $<$ a Myr for the system parameters of CI Tau. We also perform two planet scattering experiments and show that alternatively the observed planet may plausibly have acquired its eccentricity through dynamical scattering of a migrating lower mass planet, which has either been ejected from the system or swallowed by the central star. In the latter case the present location and eccentricity of the observed planet can be recovered if it was previously stalled within the disc's magnetospheric cavity.

The construction of the KAT-7 array in the Karoo region of the Northern Cape in South Africa was intended primarily as an engineering prototype for technologies and techniques applicable to the MeerKAT telescope. This paper looks at the main engineering and scien- tific highlights from this effort, and discusses their applicability to both MeerKAT and other next-generation radio telescopes. In particular we found that the composite dish surface works well, but it becomes complicated to fabricate for a dish lacking circular symmetry; the Stir- ling cycle cryogenic system with ion pump to achieve vacuum works but demands much higher maintenance than an equivalent Gifford-McMahon cycle system; the ROACH (Recon- figurable Open Architecture Computing Hardware)-based correlator with SPEAD (Stream- ing Protocol for Exchanging Astronomical Data) protocol data transfer works very well and KATCP (Karoo Array Telescope Control Protocol) control protocol has proven very flexible and convenient. KAT-7 has also been used for scientific observations where it has a niche in mapping low surface-brightness continuum sources, some extended HI halos and OH masers in star-forming regions. It can also be used to monitor continuum source variability, observe pulsars, and make VLBI observations

We investigate the minimum planet mass that produces observable signatures in infrared scattered light and submm continuum images and demonstrate how these images can be used to measure planet masses to within a factor of about two. To this end we perform multi-fluid gas and dust simulations of discs containing low mass planets, generating simulated observations at $1.65 \mu$m, $10 \mu$m and $850 \mu$m. We show that the minimum planet mass that produces a detectable signature is $\sim 15 M_\oplus$: this value is strongly dependent on disc temperature and changes slightly with wavelength (favouring the submm). We also confirm previous results that there is a minimum planet mass of $\sim 20 M_\oplus$ that produces a pressure maximum in the disc: only planets above this threshold mass generate a dust trap that can eventually create a hole in the submm dust. Below this mass, planets produce annular enhancements in dust outward of the planet and a reduction in the vicinity of the planet. These features are in steady state and can be understood in terms of variations in the dust radial velocity, imposed by the perturbed gas pressure radial profile, analogous to a traffic jam. We also show how planet masses can be derived from structure in scattered light and sub-mm images. We emphasise that simulations with dust need to be run over thousands of planetary orbits so as to allow the gas profile to achieve a steady state and caution against the estimation of planet masses using gas only simulations.

We have conducted the first comprehensive numerical investigation of the relative velocity distribution of dust particles in self-gravitating protoplanetary discs with a view to assessing the viability of planetesimal formation via direct collapse in such environments. The viability depends crucially on the large sizes that are preferentially collected in pressure maxima produced by transient spiral features (Stokes numbers, $St \sim 1$); growth to these size scales requires that collision velocities remain low enough that grain growth is not reversed by fragmentation. We show that, for a single sized dust population, velocity driving by the disc's gravitational perturbations is only effective for $St > 3$, while coupling to the gas velocity dominates otherwise. We develop a criterion for understanding this result in terms of the stopping distance being of order the disc scale height. Nevertheless, the relative velocities induced by differential radial drift in multi-sized dust populations are too high to allow the growth of silicate dust particles beyond $St \sim 10^{-2}$ or $10^{-1}$ ($10\,\mathrm{cm}$ to $\mathrm{m}$ sizes at $30\,\mathrm{au}$), such Stokes numbers being insufficient to allow concentration of solids in spiral features. However, for icy solids (which may survive collisions up to several $10\,\mathrm{m\,s}^{-1}$), growth to $St \sim 1$ ($10\,\mathrm{m}$ size) may be possible beyond $30\,\mathrm{au}$ from the star. Such objects would be concentrated in spiral features and could potentially produce larger icy planetesimals/comets by gravitational collapse. These planetesimals would acquire moderate eccentricities and remain unmodified over the remaining lifetime of the disc.

Recent interpretations of narrow, variable absorption lines detected in some Type Ia supernovae suggest that their progenitors are surrounded by dense, circumstellar material. Similar variations detected in the symbiotic recurrent nova system RS Oph, which undergoes thermonuclear outbursts every ~20 years, making it an ideal candidate to investigate the origin of these lines. To this end, we present simulations of multiple mass transfer-nova cycles in RS Oph. We find that the quiescent mass transfer produces a dense, equatorial outflow, i.e., concentrated towards the binary orbital plane, and an accretion disc forms around the white dwarf. The interaction of a spherical nova outburst with these aspherical circumstellar structures produces a bipolar outflow, similar to that seen in HST imaging of the 2006 outburst. In order to produce an ionization structure that is consistent with observations, a mass-loss rate of $5 \times 10^{-7}\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1}$ from the red giant is required. The simulations also produce a polar accretion flow, which may explain the broad wings of the quiescent H {\alpha} line and hard X-rays. By comparing simulated absorption line profiles to observations of the 2006 outburst, we are able to determine which components arise in the wind and which are due to the novae. We explore the possible behaviour of absorption line profiles as they may appear should a supernova occur in a system like RS Oph. Our models show similarities to supernovae like SN 2006X, but require a high mass-loss rate, $\dot{M} \sim 10^{-6}$ to $10^{-5}\,\mathrm{M_\odot}\,\mathrm{yr}^{-1}$, to explain the variability in SN 2006X.

We present deep neutral hydrogen (HI) observations of the starburst/Seyfert galaxy NGC 3079 and its environment, obtained with the Westerbork Synthesis Radio Telescope. Our observations reveal previously unknown components, both in HI emission and in absorption, that show that NGC 3079 is going through a hectic phase in its evolution. The HI disk appears much more extended than previously observed and is morphologically and kinematically lopsided on all scales with evidence for strong non-circular motions in the central regions. Our data reveal prominent gas streams encircling the entire galaxy suggesting strong interaction with its neighbours. A 33-kpc long HI bridge is detected between NGC 3079 and MCG 9-17-9, likely caused by ram-pressure stripping of MGC 9-17-9 by the halo of hot gas of NGC 3079. The cometary HI tail of the companion NGC 3073, earlier discovered by Irwin et al., extends about twice as long in our data, while a shorter, second tail is also found. This tail is likely caused by ram-pressure stripping by the strong, starburst driven wind coming from NGC 3079. We also detect, in absorption, a nuclear HI outflow extending to velocities well outside what expected for gravitational motion. This is likely an atomic counterpart of the well-studied outflow of ionised gas present in this galaxy. This may indicate that also large amounts of cold gas are blown out of NGC 3079 by the starburst/AGN. Our estimates of the jet energy and kinetic power suggest that both the AGN and the starburst in NGC 3079 are powerful enough to drive the atomic outflow.

We present a 'two-fluid' implementation of dust in smoothed particle hydrodynamics (SPH) in the test particle limit. The scheme is able to handle both short and long stopping times and reproduces the short friction time limit, which is not properly handled in other implementations. We apply novel tests to verify its accuracy and limitations, including multi-dimensional tests that have not been previously applied to the drag-coupled dust problem and which are particularly relevant to self-gravitating protoplanetary discs. Our tests demonstrate several key requirements for accurate simulations of gas-dust mixtures. Firstly, in standard SPH particle jitter can degrade the dust solution, even when the gas density is well reproduced. The use of integral gradients, a Wendland kernel and a large number of neighbours can control this, albeit at a greater computational cost. Secondly, when it is necessary to limit the artificial viscosity we recommend using the Cullen & Dehnen (2010) switch, since the alternative, using ${\alpha} \sim 0.1$, can generate a large velocity noise up to ${{\sigma}_v} \lesssim 0.3 c_s$ in the dust particles. Thirdly, we find that an accurate dust density estimate requires $>400$ neighbours, since, unlike the gas, the dust particles do not feel regularization forces. This density noise applies to all particle-based two-fluid implementations of dust, irrespective of the hydro solver and could lead to numerically induced fragmentation. Although our tests show accurate dusty gas simulations are possible, care must be taken to minimize the contribution from numerical noise.

Eight new class II methanol masers selected from the 6.7 GHz Methanol Multibeam survey catalogues I and II were monitored at 6.7 GHz with the 26m Hartebeesthoek Radio Astronomy Observatory (HartRAO) radio telescope for three years and seven months, from February 2011 to September 2014. The sources were also observed at 12.2 GHz and two were sufficiently bright to permit monitoring. One of the eight sources, namely G358.460-0.391, was found to show periodic variations at 6.7 GHz. The period was determined and tested for significance using the Lomb-Scargle, epoch-folding and Jurkevich methods, and by fitting a simple analytic function. The best estimate for the period of the 6.7 GHz class II methanol maser line associated with G358.460-0.391 is 220.0 $\pm$ 0.2 day.

The advent of international wideband communication by optical fibre has produced a revolution in communications and the use of the internet. Many African countries are now connected to undersea fibre linking them to other African countries and to other continents. Previously international communication was by microwave links through geostationary satellites. These are becoming redundant in some countries as optical fibre takes over, as this provides 1000 times the bandwidth of the satellite links. In the 1970's and 1980's some two dozen large (30 m diameter class) antennas were built in various African countries to provide the satellite links. Twenty six are currently known in 19 countries. As these antennas become redundant, the possibility exists to convert them for radio astronomy at a cost of roughly one tenth that of a new antenna of similar size. HartRAO, SKA Africa and the South African Department of Science and Technology (DST) have started exploring this possibility with some of the African countries.

Herbig-Haro objects are regions of shocked gas and dust which are produced when collimated outflows from a protostar interact with the surrounding dense gas. They have many similarities to supernova remnants which are interacting with molecular clouds. 1720-MHz OH masers have been identified towards a number of interacting supernova remnants. Observations and models indicate that these masers are shock excited and are produced behind C-type shocks. If conditions behind the shock fronts of Herbig-Haro objects are similarly able to support 1720-MHz OH masers they could be a useful diagnostic tool for star formation. We therefore searched for 1720-MHz OH maser emission towards a sample of 97 Herbig-Haro objects using the Green Bank radio telescope. We detected 1720-MHz OH lines in emission in 17 of them, but neither their spectral signature nor follow-up observations with the Very Large Array showed any conclusive evidence of maser emission. We conclude that the emission detected from our single-dish observations must be extended and most likely originates from thermal or quasi-thermal excitation processes. We also investigated the properties of Herbig-Haro shocks more closely and conclude that despite the overall similarities to supernova remnants, the conditions required for maser emission, in particular, a sufficient velocity-coherent column density, are not likely to occur in Herbig-Haro objects.

Circinus X-1 is a bright and highly variable X-ray binary which displays strong and rapid evolution in all wavebands. Radio flaring, associated with the production of a relativistic jet, occurs periodically on a ~17-day timescale. A longer-term envelope modulates the peak radio fluxes in flares, ranging from peaks in excess of a Jansky in the 1970s to an historic low of milliJanskys during the years 1994 to 2007. Here we report first observations of this source with the MeerKAT test array, KAT-7, part of the pathfinder development for the African dish component of the Square Kilometre Array (SKA), demonstrating successful scientific operation for variable and transient sources with the test array. The KAT-7 observations at 1.9 GHz during the period 13 December 2011 to 16 January 2012 reveal in temporal detail the return to the Jansky-level events observed in the 1970s. We compare these data to contemporaneous single-dish measurements at 4.8 and 8.5 GHz with the HartRAO 26-m telescope and X-ray monitoring from MAXI. We discuss whether the overall modulation and recent dramatic brightening is likely to be due to an increase in the power of the jet due to changes in accretion rate or changing Doppler boosting associated with a varying angle to the line of sight.

We have mapped faint 1667 OH line emission (TA \approx 20 - 40 mK in our \approx 30' beam) along many lines of sight in the Galaxy covering an area of \approx 4\circ \times 4\circ in the general direction of l \approx 108\circ, b \approx 5\circ. The OH emission is widespread, similar in extent to the local HI (r </= 2 kpc) both in space and in velocity. The OH profile amplitudes show a good general correlation with those of HI in spectral channels of \approx 1 km/s; this relation is described by TA(OH) \approx 1.50 \times 10^{-4} TB(HI) for values of TB(HI) </\approx 60 - 70 K. Beyond this the HI line appears to "saturate", and few values are recorded above \approx 90 K. However, the OH brightness continues to rise, by a further factor \approx 3. The OH velocity profiles show multiple features with widths typically 2 - 3 km/s, but less than 10% of these features are associated with CO(1-0) emission in existing surveys of the area smoothed to comparable resolution.

We present a detailed multi-epoch analysis of 31 potential southern hemisphere radio calibrators that were originally observed as part of a program to maintain the International Celestial Reference Frame (ICRF). At radio wavelengths, the primary calibrators are Active Galactic Nuclei (AGN), powerful radio emitters which exist at the centre of most galaxies. These are known to vary at all wavelengths at which they have been observed. By determining the amount of radio source structure and variability of these AGN, we determine their suitability as phase calibrators for long baseline radio interferometry at 2.3 GHz. For this purpose, we have used a set of complementary metrics to classify these 31 southern sources into five categories pertaining to their suitability as VLBI calibrators. We find that all of the sources in our sample would be good interferometric calibrators and almost ninety per cent would be very good calibrators.

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

One of the defining characteristics of BL Lacertae objects is their strong variability across the electromagnetic spectrum. PKS 0537-441 is one such object and is one of the most luminous blazars from radio to gamma-ray wavelengths. It was detected as a strong and highly variable source by EGRET and has been reported several times to be in an active state by Fermi . It is one of the brightest gamma-ray blazars detected in the southern sky so far. The TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) program is monitoring PKS 0537-441 at VLBI resolutions. We present 8.4 GHz and 22 GHz images of the milliarcsecond scale structure. We also present our ongoing analysis of the spectral and temporal changes in this object.

We introduce the TANAMI program (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) which is monitoring an initial sample of 43 extragalactic jets located south of -30 degrees declination at 8.4 GHz and 22 GHz since 2007. All aspects of the program are discussed. First epoch results at 8.4 GHz are presented along with physical parameters derived therefrom. We present first epoch images for 43 sources, some observed for the first time at milliarcsecond resolution. Parameters of these images as well as physical parameters derived from them are also presented and discussed. These and subsequent images from the TANAMI survey are available at http://pulsar.sternwarte.uni-erlangen.de/tanami/ We obtain reliable, high dynamic range images of the southern hemisphere AGN. All the quasars and BL Lac objects in the sample have a single-sided radio morphology. Galaxies are either double-sided, single-sided or irregular. About 28% of the TANAMI sample has been detected by LAT during its first three months of operations. Initial analysis suggests that when galaxies are excluded, sources detected by LAT have larger opening angles than those not detected by LAT. Brightness temperatures of LAT detections and non-detections seem to have similar distributions. The redshift distributions of the TANAMI sample and sub-samples are similar to those seen for the bright gamma-ray AGN seen by LAT and EGRET but none of the sources with a redshift above 1.8 have been detected by LAT.

We present a summary of the observation strategy of TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry), a monitoring program to study the parsec-scale structure and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Australian Long Baseline Array (LBA) and the trans-oceanic antennas Hartebeesthoek, TIGO, and O'Higgins. TANAMI is focusing on extragalactic sources south of -30 degrees declination with observations at 8.4 GHz and 22 GHz every ~2 months at milliarcsecond resolution. The initial TANAMI sample of 43 sources has been defined before the launch of the Fermi Gamma Ray Space Telescope to include the most promising candidates for bright gamma-ray emission to be detected with its Large Area Telescope (LAT). Since November 2008, we have been adding new sources to the sample, which now includes all known radio- and gamma-ray bright AGN of the Southern Hemisphere. The combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN and the nature of relativistic jets.

The TANAMI (Tracking AGN with Austral Milliarcsecond Interferometry) and associated programs provide comprehensive radio monitoring of extragalactic gamma-ray sources south of declination -30 degrees. Joint quasi-simultaneous observations between the Fermi Gamma-ray Space Telescope and ground based observatories allow us to discriminate between competing theoretical blazar emission models. High resolution VLBI observations are the only way to spatially resolve the sub-parsec level emission regions where the high-energy radiation originates. The gap from radio to gamma-ray energies is spanned with near simultaneous data from the Swift satellite and ground based optical observatories. We present early results from the TANAMI program in the context of this panchromatic suite of observations.

TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) is a monitoring program to study the parsec-scale structures and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Long Baseline Array and associated telescopes. Extragalactic jets south of -30 degrees declination are observed at 8.4 GHz and 22 GHz every two months at milliarcsecond resolution. The initial TANAMI sample is a hybrid radio and gamma-ray selected sample since the combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN.

We present the specifications of the MeerKAT Karoo Array Telescope, the South African Square Kilometre Array Precursor. Some of the key science for MeerKAT is described in this document. We invite the community to submit proposals for Large Key Projects.

Very deep neutral hydrogen (HI) observations of the edge-on spiral galaxy NGC 3079 with the Westerbork Synthesis Radio Telescope (WSRT) are presented. The galaxy has been studied extensively in different wavelengths and known for its several unique and complex features. However, the new data still revealed several new features and show that the galaxy is even more complicated and interesting than previously thought. In the new data a large stream of gas, encircling the entire galaxy, was discovered, while the data also show, for the first time, that not only hot gas is blown into space by the starburst/AGN, but also large amounts of cold gas, despite the high energies involved in the outflow.

A fundamental assumption inherent in the standard Lambda+CDM Hot Big Bang (HBB) model is that photons lose energy as they are redshifted due to the expansion of the universe. We show that for the Quasi-Static Universe (QSU) model, in which photon energy is an invariant in the cosmological reference frame, the photon number density in the universe today is a factor of approximately 1600 less than in the standard model. We examine some of the consequences for a number of processes that occur during the thermal history of the early universe, including primordial nucleosynthesis, the formation of neutral hydrogen (recombination), and the evolution of the Cosmic Microwave Background (CMB) radiation. We show that the QSU model predicts that the measured CMB photon noise level will be a factor of 40 higher than the level that would be observed assuming the standard HBB model. The CMB data that will be collected by the recently launched Planck satellite mission provides an ideal opportunity to test the validity of this prediction.

In a study of the radio emission mechanism of the FR-I AGN NGC 5128 (Centaurus A)}, we have determined the centimeter and millimeter continuum spectrum of the whole Centaurus A radio source and measured the continuum emission from the active galxy nucleus at various times between 1989 and 2005 at frequencies between 86 GHz (3.5 mm) and 345 GHz (0.85 mm). The integral Cen A spectrum becomes steeper at frequencies above 5 GHz, where the spectral index changes from -0.70 to -0.82. Millimeter emission from the core of Centaurus A is variable, and correlates appreciably better with the 20-200 keV than the 2 - 10 keV X-ray variability. In its quiescent state, the core spectral index is -0.3, which steepens when the core brightens. The variability appears to be mostly associated with the inner nuclear jet components that have been detected in VLBI measurements. The densest nuclear components are optically thick below 45-80 GHz.

We have carried out sub-mm 12CO(J=3-2) observations of 6 giant molecular clouds (GMCs) in the Large Magellanic Cloud (LMC) with the ASTE 10m sub-mm telescope at a spatial resolution of 5 pc and very high sensitivity. We have identified 32 molecular clumps in the GMCs and revealed significant details of the warm and dense molecular gas with n(H2) $\sim$ 10$^{3-5}$ cm$^{-3}$ and Tkin $\sim$ 60 K. These data are combined with 12CO(J=1-0) and 13CO(J=1-0) results and compared with LVG calculations. We found that the ratio of 12CO(J=3-2) to 12CO(J=1-0) emission is sensitive to and is well correlated with the local Halpha flux. We interpret that differences of clump propeties represent an evolutionary sequence of GMCs in terms of density increase leading to star formation.Type I and II GMCs (starless GMCs and GMCs with HII regions only, respectively) are at the young phase of star formation where density does not yet become high enough to show active star formation and Type III GMCs (GMCs with HII regions and young star clusters) represents the later phase where the average density is increased and the GMCs are forming massive stars. The high kinetic temperature correlated with \Halpha flux suggests that FUV heating is dominant in the molecular gas of the LMC.

We have searched for emission from the 557 GHz ortho-water line in the interstellar medium of six nearby starburst galaxies. We used the Odin satellite to observe the 1_10-1_01 transition of o-H2O in the galaxies NGC253, IC342, M82, NGC4258, CenA, and M51. None of the galaxies in our sample was detected. We derive three sigma upper limits to the H2O abundance relative to H2 ranging from 2e-9 to 1e-8. The best of these upper limits are comparable to the measured abundance of H$_2$O in the Galactic star forming region W3. However, if only 10% of the molecular gas is in very dense cores, then the water abundance limits in the cores themselves would be larger by a factor of 10 i.e. 2e-8 to 1e-7. These observations suggest that detections of H2O emission in galaxies with the upcoming Herschel Space Observatory are likely to require on-source integration times of an hour or more except in the very brightest extragalactic targets such as M82 and NGC253.

Molecular oxygen, O2 has been expected historically to be an abundant component of the chemical species in molecular clouds and, as such, an important coolant of the dense interstellar medium. However, a number of attempts from both ground and from space have failed to detect O2 emission. The work described here uses heterodyne spectroscopy from space to search for molecular oxygen in the interstellar medium. The Odin satellite carries a 1.1 m sub-millimeter dish and a dedicated 119 GHz receiver for the ground state line of O2. Starting in 2002, the star forming molecular cloud core rho Oph A was observed with Odin for 34 days during several observing runs. We detect a spectral line at v(LSR) = 3.5 km/s with dv(FWHM) = 1.5 km/s, parameters which are also common to other species associated with rho Ohp A. This feature is identified as the O2 (N_J = 1_1 - 1_0) transition at 118 750.343 MHz. The abundance of molecular oxygen, relative to H2,, is 5E-8 averaged over the Odin beam. This abundance is consistently lower than previously reported upper limits.

Context: Methanol masers at 6.7GHz are recognised markers of high-mass star formation regions. The study of their distribution in the Galaxy gives important insights into the star formation activity of the Milky Way. We present a statistical analysis on the General Catalogue of 6.7GHz methanol masers in the Galaxy with the aim of extracting global properties of the masers. Aims: We provide constraints on the luminosity function of 6.7GHz methanol masers and on their total number in the Galaxy. Methods: We model the spatial distribution of the masers in the Milky Way by using their distribution in galactocentric distance which is unambiguous once a rotation curve for the Galaxy is assumed. This is the starting point for determining the luminosity function of the masers. Results: The luminosity function of 6.7GHz methanol masers is modelled as a power-law with sharp cutoffs and having an index lying between -1.5 and -2. We also predict the number of detections of methanol masers assuming different sensitivity limits in the observations.

We present a 39 h integration with the Odin satellite on the ground-state 118.75 GHz line of O2 towards the region of strongest molecular emission in the Small Magellanic Cloud. Our 3sigma upper limit to the O2 integrated intensity of <0.049 K km/s in a 9'(160 pc) diameter beam corresponds to an upper limit on the O2/H2 abundance ratio of <1.3E-6. Although a factor of 20 above the best limit on the O2 abundance obtained for a Galactic source, our result has interesting implications for understanding oxygen chemistry at sub-solar metal abundances. We compare our abundance limit to a variety of astrochemical models and find that, at low metallicities, the low O2 abundance is most likely produced by the effects of photo-dissociation on molecular cloud structure. Freeze-out of molecules onto dust grains may also be consistent with the observed abundance limit, although such models have not yet been run at sub-solar initial metallicities.

Several new main line OH masers have been detected in the nearby starburst galaxy M82. Eight masers have been detected to 5 sigma, six of which are new detections. Observations covering both the 1665 and 1667 MHz lines with both the Very Large Array (VLA) and the Multi-Element Radio Linked Interferometer Network (MERLIN) have been used to accurately measure the positions and velocities of these features. Following analysis of the data, another six objects below 5 sigma, but with velocities consistent with the distribution inferred from the more certain detections, have been detected. These are classified as possible detections.

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.

Methanol masers are often detected in regions of intense star formation. Several studies in the last decade indicate that they may even be the earliest signpost of a high-mass star-forming region. Their powerful emission make them very good candidates for observations using both single--dish telescopes and interferometers, the latter allows detailed structural and dynamical studies of these objects. We have prepared a catalogue of all known 6.7-GHz methanol masers, discovered both by surveys that targeted possible associated objects and unbiased surveys covering a large fraction of galactic longitudes across the Galactic plane (-0.5deg <= b <= 0.5deg for most of the regions). The catalogue contains 519 sources which are listed with their kinematic (galactocentric and heliocentric) distances as well as possibly associated IR objects. We find that 6.7-GHz methanol masers clearly trace the molecular ring of our Galaxy, where most of the OB associations are located. The present list of masers also reports detections of other masing transitions of methanol as further information for the study of the maser phenomenon. In a further publication we will address some statistical considerations.

We summarize the present status of VLBI experiments at 3 mm (86 GHz), 2 mm (129-150 GHz) and 1.3 mm (215-230 GHz). We present and discuss a new 3 mm VLBI map of M87 (Virgo A), which has a spatial resolution of only approx. 20 Schwarzschild radii. We discuss recent VLBI results for SgrA* and argue in favor of new observations within an extended European mm-VLBI network, in order to search for variability. We discuss the possibilities to image the `event horizon' of a super-massive black hole at wavelengths < 2mm, and conclude that the addition of large and sensitive millimetre telescopes such as CARMA, the SMA, the LMT and ALMA will be crucial for this.

Global millimetre VLBI allows detailed studies of the most central jet regions of AGN with unprecedent spatial resolution of a few 100-1000 Schwartzschild radii to be made. Study of these regions will help to answer the question how the highly relativistic AGN jets are launched and collimated. Since the early 1990s, bright mm-sources have been observed with global 3 mm VLBI. Here we present new images from an ongoing systematic analysis of the available observations. In particular, we focus on the structure and structural evolution of the best observed AGN jets, taking 3C 454.3 as a characteristic example. This core-dominated and highly variable quasar shows a complex morphology with individual jet components accelerating superluminally towards the outer structure. We briefly discuss the X-ray properties of 3C 454.3 and present its radio- to X-ray large-scale brightness distribution.

We present an edge-on Keplerian disk model to explain the main component of the 12.2 and 6.7 GHz methanol maser emission detected toward NGC7538-IRS1 N. The brightness distribution and spectrum of the line of bright masers are successfully modeled with high amplification of background radio continuum emission along velocity coherent paths through a maser disk. The bend seen in the position-velocity diagram is a characteristic signature of differentially rotating disks. For a central mass of 30 solar masses, suggested by other observations, our model fixes the masing disk to have inner and outer radii of about 270 AU and 750 AU.

We present an edge-on Keplerian disc model to explain the main component of the 12.2 and 6.7 GHz methanol maser emission detected toward NGC7538-IRS1 N. The brightness distribution and spectrum of the line of bright masers are successfully modeled with high amplification of background radio continuum emission along velocity coherent paths through a maser disc. The bend seen in the position-velocity diagram is a characteristic signature of differentially rotating discs. For a central mass of 30 solar masses, suggested by other observations, our model fixes the masing disc to have inner and outer radii of about 350 AU and 1000 AU.

A prediction of standard Big Bang cosmology is that the observed UHECR (ultra-high-energy cosmic rays) spectrum will exhibit a cutoff at the GKZ limit, resulting from interaction with the photons that constitute the cosmic microwave background. We show that for the Quasi-Static Universe (QSU) model, in which photon energy is an invariant in the cosmological reference frame, the photon number density in the universe today is a factor of 10^9 less than in the standard model. As a consequence, the mean free path of UHECRs will exceed the horizon distance of the universe, rendering it essentially transparent to UHECRs. The QSU model therefore predicts that no cutoff will be observed in the UHECR spectrum.

We present J=1-0 and J=2-1 12CO maps of several star-forming regions in both the Large and the Small Magellanic Cloud, and briefly discuss their structure. Many of the detected molecular clouds are relatively isolated and quite small with dimensions of typically 20 pc. Some larger complexes have been detected, but in all cases the extent of the molecular clouds sampled by CO emission is significantly less than the extent of the ionized gas of the star-formation region. Very little diffuse extended CO emission was seen; diffuse CO in between or surrounding the detected discrete clouds is either very weak or absent. The majority of all LMC lines of sight detected in 13CO has an isotopic emission ratio I(12CO)/I(13CO) of about 10, i.e. twice higher than found in Galactic star-forming complexes. At the lowest 12CO intensities, the spread of isotopic emission ratios rapidly increases, low ratios representing relatively dense and cold molecular gas and high ratios marking CO photo-dissociation at cloud edges.

We have mapped the submillimeter emission from the 1(10)-1(01) transition of ortho-water in the W3 star-forming region. A 5'x5' map of the W3 IRS4 and W3 IRS5 region reveals strong water lines at half the positions in the map. The relative strength of the Odin lines compared to previous observations by SWAS suggests that we are seeing water emission from an extended region. Across much of the map the lines are double-peaked, with an absorption feature at -39 km/s; however, some positions in the map show a single strong line at -43 km/s. We interpret the double-peaked lines as arising from optically thick, self-absorbed water emission near the W3 IRS5, while the narrower blue-shifted lines originate in emission near W3 IRS4. In this model, the unusual appearance of the spectral lines across the map results from a coincidental agreement in velocity between the emission near W3 IRS4 and the blue peak of the more complex lines near W3 IRS5. The strength of the water lines near W3 IRS4 suggests we may be seeing water emission enhanced in a photon-dominated region.

Odin has successfully observed the molecular core rho Oph A in the 572.5 GHz rotational ground state line of ammonia, NH3 (J,K = 1,0 - 0,0). The interpretation of this result makes use of complementary molecular line data obtained from the ground (C17O and CH3OH) as part of the Odin preparatory work. Comparison of these observations with theoretical model calculations of line excitation and transfer yields a quite ordinary abundance of methanol, X(CH3OH) = 3e-9. Unless NH3 is not entirely segregated from C17O and CH3OH, ammonia is found to be significantly underabundant with respect to typical dense core values, viz. X(NH3) = 8e-10.

Odin has successfully observed three regions in the Orion A cloud, i.e. Ori KL, Ori S and the Orion Bar, in the 572.5 GHz rotational ground state line of ammonia, ortho-NH3 (J,K) = (1,0) -> (0,0), and the result for the Orion Bar represents the first detection in an ammonia line. Several velocity components are present in the data. Specifically, the observed line profile from the Orion Bar can be decomposed into two components, which are in agreement with observations in high-J CO lines by Wilson et al. 2001. Using the source model for the Orion Bar by these authors, our Odin observation implies a total ammonia abundance of NH3/H2 = 5E-9.

We present the first VLBI maps of H2O maser emission (lambda 1.3cm) in the nucleus of the Circinus Galaxy, constructed from data obtained with the Australia Telescope Long Baseline Array. The maser emission traces a warped, edge-on accretion disk between radii of 0.11+/-0.02 and ~0.40 pc, as well as a wide-angle outflow that extends up to ~1 pc from the estimated disk center. The disk rotation is close to Keplerian (v varies as 1/sqrt(r)), the maximum detected rotation speed is 260 km/s, and the inferred central mass is 1.7+/-0.3 x 10^6 solar masses. The outflowing masers are irregularly distributed above and below the disk, with relative outflow velocities up to ~+/-160 km/s, projected along the line of sight. The flow probably originates closer than 0.1 pc to the central engine, possibly in an inward extension of the accretion disk, though there is only weak evidence of rotation in the outward moving material. We observe that the warp of the disk appears to collimate the outflow and to fix the extent of the ionization cone observed on larger angular scales. This study provides the first direct evidence (i.e., through imaging) of dusty, high-density, molecular material in a nuclear outflow <1 pc from the central engine of a Seyfert galaxy, as well as the first graphic evidence that warped accretion disks can channel outflows and illumination patterns in AGN. We speculate that the same arrangement, which in some ways obviates the need for a geometrically thick, dusty torus, may apply to other type-2 AGN.

We have used the Australia Telescope Compact Array to image the neutral hydrogen in the merging system NGC 3256, to test the idea that globular clusters (GC) form during the interactions and mergers of disk galaxies. We compare our observations with hydrodynamical numerical simulations, from the literature, to examine the hypothesis that the H I fragments with masses greater than 10$^{7\pm1} \mo$ are sites of GC formation. We appear to have detected detached H I fragments in the vicinity of NGC 3256. These fragments, with masses ~10^7 solar masses, may have little dark matter content which is also a characteristic of globular clusters, and so our observations support the hypothesis that globular clusters form in the type of interaction that resulted in NGC 3256.

We present maps and a catalogue containing the J=1-0 12CO parameters of 29 individual molecular clouds in the second-brightest LMC star formation complex, N11. In the southwestern part of N11, molecular clouds occur in a ring or shell surrounding the major OB star association LH9. In the northeastern part, a chain of molecular clouds delineates the rim of one of the so-called supergiant shells in the LMC. The well-defined clouds have dimensions close to those of the survey beam (diameters of 25 pc or less). Some of the clouds were also observed in J=2-1 12CO, and in the lower two transitions of 13CO. Clouds mapped with a twice higher angular resolution in J=2-1 12CO show substructure with dimensions once again comparable to those of the mapping beam. The few clouds for which we could model physical parameters have fairly warm (T(kin) = 60 - 150 K) and moderately dense (n(H2) = 3000 cm-3) gas. The northeastern chain of CO clouds, although lacking in diffuse intercloud emission, is characteristic of the more quiescent regions of the LMC and appears to have been subject to relatively little photo-processing. The clouds forming part of the southwestern shell or ring, however, are almost devoid of diffuse intercloud emission and also exhibit other characteristics of an extreme photon-dominated region (PDR).

At 147 GHz (2mm wavelength), we detected three prominent AGN (NRAO150, 3C279, 1633+382) with Very Long Baseline Interferometry (VLBI) with an angular resolution of only 18 micro-arcseconds on the baseline between two antennas in Arizona (10m HHT and 12m KittPeak) and the IRAM 30m antenna on Pico Veleta in Spain. This is a new world record in radio interferometry and astronomical imaging and opens fascinating future possibilities to directly image and study the innermost regions in Quasars and other Active Galactic Nuclei.

We report SIMBA 1.2 mm dust continuum observations of the environments of eight methanol maser sources, all discovered during spatially fully-sampled, untargeted surveys of the galactic plane. We summarise our search for possible associations of the masers with IR sources (IRAS and MSX) and find that it is not always possible to make definite associations. A preliminary characterisation of the IR sources found in the maser neighbourhood is given according their position in the [60-25] -- [25-12] colour-colour diagram.

Observations of an apparent acceleration in the expansion rate of the universe, derived from measurements of high-redshift supernovae, have been used to support the hypothesis that the universe is permeated by some form of dark energy. We show that an alternative cosmological model, based on a linearly expanding universe with Omega=1, can fully account for these observations. This model is also able to resolve the other problems associated with the standard Big Bang model. A scale-invariant form of the field equations of General Relativity is postulated, based on the replacement of the Newtonian gravitational constant by an formulation based explicitly on Mach's principle. We derive the resulting dynamical equations, and show that their solutions describe a linearly expanding universe. Potential problems with non-zero divergencies in the modified field equations are shown to be resolved by adopting a radically different definition of time, in which the unit of time is a function of the scale-factor. We observe that the effects of the modified field equations are broadly equivalent to a Newtonian gravitational constant that varies both in time and in space, and show that this is also equivalent to Varying Speed of Light (VSL) theories using a standard definition of time. Some of the implications of this observation are discussed in relation to Black Holes, Planck scale phenomena, and the ultimate fate of the Universe.

We review the state of the Onsala blind survey of the galactic plane, searching for new 6.7 GHz methanol masers. We also describe preliminary results of millimeter follow-up observations of the new detections and high resolution observations using the EVN. We conclude that blind surveys are important to complement targeted searches done until now and give the possibility to detect new classes of objects.

The first VLBI images of water maser emission in the Circinus Galaxy AGN show both a warped, edge-on accretion disk and an outflow 0.1 to 1 pc from the central engine. The inferred central mass is 1.3 million suns, while the disk mass may be on the order of 0.1 million suns, based on a nearly Keplerian rotation curve. The bipolar, wide-angle outflow appears to contain ``bullets'' ejected from within <0.1 pc of the central mass. The positions of filaments and bullets observed in the AGN ionization cone on kpc-scales suggest that the disk channels the flow to a radius of about 0.4 pc, at which the flow appears to disrupt the disk.

We present the analysis of a survey of atomic and molecular gas in interacting and merging galaxies (Horellou & Booth 1997). The sample is optically selected and contains all interacting galaxies (~ 60 systems) in a well-defined region of the Southern sky (Bergvall 1981). In order to take into account the upper limits due to non-detections, we have carried out a survival analysis. The results are compared with the ones obtained when ignoring the non-detections. We have found evidence for a lower atomic gas content per unit area in interacting galaxies compared to isolated ones. Except for a few systems with high far-infrared luminosity and correspondingly high CO fluxes, the interacting galaxies are not unusually CO-bright. The observed HI-deficiency is therefore not due to a conversion from HI into H2. Atomic gas from the outer part of the galaxies may have been swept by tidal interaction, or the optical extent of the galaxies may have increased due to the tidal interaction. The interacting galaxies seem to be more efficient at transforming their molecular gas into stars, as indicated by their higher LFIR/M(H2) ratio. We found no correlation between molecular gas fraction and degree of morphological distortion.

We report simultaneous observations at 1.2 and 2 mm, with a double channel photometer on the SEST Telescope, of the X-ray cluster RXJ0658-5557 in search for the Sunyaev-Zel'dovich (S-Z). The S-Z data were analyzed using the relativistically correct expression for the Comptonization parameter and we find from the detected decrement (2.60 +/- 0.79) ~ 10^{-4}, which is consistent with that computed using the X-ray (ROSAT and ASCA) observations. The uncertainty includes contributions due to statistical uncertainty in the detection and systematics and calibration. The 1.2 {mm} channel data alone gives rise to a larger Comptonization parameter and this result is discussed in terms of contamination from foreground sources and/or dust in the cluster or from a possible systematic effect. We then make use of the combined analysis of the ROSAT and ASCA X-ray satellite observations to determine an isothermal model for the S-Z surface brightness. Since the cluster is asymmetrical and probably in a merging process, models are only approximate. The associated uncertainty can, however, be estimated by exploring a set of alternative models. We then find as the global uncertainty on the Comptonization parameter a factor of 1.3. Combining the S-Z and the X-ray measurements, we determine a value for the Hubble constant. The 2 mm data are consistent with H_0(q_0 = 1/2)= 53^{+38}_{-28} km/s Mpc^{-1}, where the uncertainty is dominated by the uncertainty in models of the X-ray plasma halo.

The chemical structure of neutral clouds in low metallicity environments is examined with particular emphasis on the H to H_2 and C+ to CO transitions. We observed near-IR H_2 lines and the CO J=1-0 line from 30 Doradus and N159/N160 in the Large Magellanic Cloud and from DEM S 16, DEM S 37, and LI-SMC 36 in the Small Magellanic Cloud. We find that the H_2 emission is UV-excited and that (weak) CO emission always exists (in our surveyed regions) toward positions where H_2 and [CII] emission have been detected. Using a PDR code and a radiative transfer code, we simulate the emission of line radiation from spherical clouds and from large planar clouds. Because the [CII] emission and H_2 emission arise on the surface of the cloud and the lines are optically thin, these lines are not affected by changes in the relative sizes of the neutral cloud and the CO bearing core, while the optically thick CO emission can be strongly affected. The sizes of clouds are estimated by measuring the deviation of CO emission strength from that predicted by a planar cloud model of a given size. The average cloud column density and therefore size increases as the metallicity decreases. Our result agrees with the photoionization regulated star formation theory by Mc Kee (1989).

An observing campaign was devoted to the search for the Sunyaev-Zeldovich (S-Z) effect towards X-ray ROSAT Clusters in the millimetric spectral domain. A double channel (1.2 and 2 {\it mm}) photometer was installed at the focus of the 15m Swedish ESO Submillimeter Telescope (SEST) in Chile in september 1994 and 1995 and observations of the targets S1077, A2744, S295 and RXJ0658-5557 were gathered. Detections were found for A2744 at 1 {\it mm} and in both channels (at 1.2 and 2 {\it mm}) towards RXJ0658-5557. For the first time there is evidence for the S-Z enhancement and both the latter and the decrement were detected on the same source. We discuss astrophysical and systematic effects which could give origin to these signals.

We report on observations of the Sunyaev-Zeldovich effect towards X-ray ROSAT clusters taken with a double channel (1.2 and 2 mm) photometer installed at the focus of the 15m SEST antenna in Chile. This paper describes the first results obtained for the high-z clusters S1077, A2744 and S295. Marginal detections were found for A2744 and at 1 mm for S1077. We discuss these data in terms of contamination of sources along the line of sight and give a constraint on the amplitude of the kinematic effect.