We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

We investigate the spatial properties of sources from the GAIA catalogue previously identified as being members of the Taurus star forming region and which appear in the Spitzer catalogue. We study an area of sky of 10x15 deg, centred on Right Ascension (2000.0) = 68.5 deg and Declination (2000.0)= 27.0 deg, this being an area surrounding the Taurus molecular cloud. We use data obtained from the GAIA DR2 release. By using an inversion of GAIA parallax measurements to obtain distance values and by defining limits to the proper motions of the Taurus moving group, we are able to show that there are substantial differences in depth within the Taurus complex. Our results suggest that the Taurus cloud has significant depth and that there are two main associations centred at ~130.6+/-0.7 pc and ~160.2+/-0.9 pc at 1{\sigma}. These two associations also have different proper motions, of 24.5+/-2.8 and 20.1+/-2.4 mas/yr respectively. We here label them the 'Two Horns' of Taurus.

Herschel observations of nearby clouds in the Gould Belt support a paradigm for low-mass star formation, starting with the generation of molecular filaments, followed by filament fragmentation, and the concentration of mass into self-gravitating prestellar cores. With the unique far-infrared and submillimeter continuum imaging capabilities of the Herschel Space observatory, the closeby (d = 139 pc) Ophiuchus cloud was mapped at five wavelengths from 70 ${\mu}m$ to 500 ${\mu}m$ with the aim of providing a complete census of dense cores in this region, including unbound starless cores, bound prestellar cores, and protostellar cores. Taking advantage of the high dynamic range and multi-wavelength nature of the Herschel data, we used the multi-scale decomposition algorithms getsources and getfilaments to identify a complete sample of dense cores and filaments in the cloud and study their properties. The densest clouds of the Ophiuchus complex, L1688 and L1689, which thus far are only indirectly described as filamentary regions owing to the spatial distribution of their young stellar objects (YSOs), are confirmed to be dominated by filamentary structures. The tight correlation observed between prestellar cores and filamentary structures in L1688 and L1689 supports the view that solar-type star formation occurs primarily in dense filaments. While the sub clouds of the complex show disparities, L1689 being less efficient than L1688 at forming stars when considering their total mass budgets, both sub clouds share almost the same prestellar core formation efficiency in dense molecular gas. We also find evidence in the Herschel data for a remarkable concentric geometrical configuration in L1688 which is dominated by up to three arc-like compression fronts and presumably created by shockwave events emanating from the Sco OB2 association, including the neighboring massive (O9V) star ${\sigma}$ Sco.

We present the detection of helium in the extended atmosphere of the sub-Saturn WASP-107b using high resolution ($R \approx 25000$) near-infrared spectra from Keck II/NIRSPEC. We find peak excess absorption of $7.26 \pm 0.24\%$ (30$\sigma$) centered on the HeI triplet at 10833A. The amplitude and shape of the helium absorption profile is in excellent agreement with previous observations of escaping helium from this planet made by CARMENES and HST. This suggests there is no significant temporal variation in the signature of escaping helium from the planet over a two year baseline. This result demonstrates Keck II/NIRSPEC's ability to detect atmospheric escape in exoplanets, making it a useful instrument to further our understanding of the evaporation of exoplanetary atmospheres via ground-based observations of HeI.

We present a new ground-based visual transmission spectrum of the hot Jupiter WASP-43b, obtained as part of the ACCESS Survey. The spectrum was derived from four transits observed between 2015 and 2018, with combined wavelength coverage between 5,300 \r{A}-9,000 \r{A} and an average photometric precision of 708 ppm in 230 \r{A} bins. We perform an atmospheric retrieval of our transmission spectrum combined with literature HST/WFC3 observations to search for the presence of clouds/hazes as well as Na, K, H$\alpha$, and $\rm H_2O$ planetary absorption and stellar spot contamination over a combined spectral range of 5,318 \r{A}-16,420 \r{A}. We do not detect a statistically significant presence of Na I or K I alkali lines, or H$\alpha$ in the atmosphere of WASP-43b. We find that the observed transmission spectrum can be best explained by a combination of heterogeneities on the photosphere of the host star and a clear planetary atmosphere with $\rm H_2O$. This model yields a log-evidence of $8.26\pm0.42$ higher than a flat (featureless) spectrum. In particular, the observations marginally favor the presence of large, low-contrast spots over the four ACCESS transit epochs with an average covering fraction $f_\text{het} = 0.27^{+0.42}_{-0.16}$ and temperature contrast $\Delta T = 132\text{ K} \pm 132\text{ K}$. Within the planet's atmosphere, we recover a log $\rm H_2O$ volume mixing ratio of $-2.78^{+1.38}_{-1.47}$, which is consistent with previous $\rm H_2O$ abundance determinations for this planet.

The detection of low-mass transiting exoplanets in multiple systems brings new constraints to planetary formation and evolution processes and challenges the current planet formation theories. Nevertheless, only a mere fraction of the small planets detected by Kepler and K2 have precise mass measurements, which are mandatory to constrain their composition. We aim to characterise the planets that orbit the relatively bright star K2-138. This system is dynamically particular as it presents the longest chain known to date of planets close to the 3:2 resonance. We obtained 215 HARPS spectra from which we derived the radial-velocity variations of K2-138. Via a joint Bayesian analysis of both the K2 photometry and HARPS radial-velocities (RVs), we constrained the parameters of the six planets in orbit. The masses of the four inner planets, from b to e, are 3.1, 6.3, 7.9, and 13.0 $\mathrm{M}_{\oplus}$ with a precision of 34%, 20%, 18%, and 15%, respectively. The bulk densities are 4.9, 2.8, 3.2, and 1.8 g cm$^{-3}$, ranging from Earth to Neptune-like values. For planets f and g, we report upper limits. Finally, we predict transit timing variations of the order two to six minutes from the masses derived. Given its peculiar dynamics, K2-138 is an ideal target for transit timing variation (TTV) measurements from space with the upcoming CHaracterizing ExOPlanet Satellite (CHEOPS) to study this highly-packed system and compare TTV and RV masses.

We present a ground-based transmission spectrum and comprehensive retrieval analysis of the highly inflated Saturn-mass planet WASP-39b. We obtained low-resolution spectra ($R \approx 400$) of a transit of WASP-39b using the ACAM instrument on the 4.2m William Herschel Telescope as part of the LRG-BEASTS survey. Our transmission spectrum is in good agreement with previous ground- and space-based observations of WASP-39b, and covers a wavelength range of 4000-9000A. Previous analyses of this exoplanet have retrieved water abundances that span more than four orders of magnitude, which in turn lead to conclusions of a subsolar or highly supersolar atmospheric metallicity. In order to determine the cause of the large discrepancies in the literature regarding WASP-39b's atmospheric metallicity, we performed retrieval analyses of all literature data sets. Our retrievals, which assume equilibrium chemistry, recovered highly supersolar metallicities for all data sets. When running our retrievals on a combined spectrum, spanning 0.3-5$\mu$m, we recovered an atmospheric metallicity of $282^{+65}_{-58} \times$ solar. We find that stellar activity has a negligible effect on the derived abundances and instead conclude that different assumptions made during retrieval analyses lead to the reported water abundances that differ by orders of magnitude. This in turn has significant consequences for the conclusions we draw. This is the fourth planet to be observed as part of the LRG-BEASTS survey, which is demonstrating that 4m class telescopes can obtain low-resolution transmission spectra with precisions of around one atmospheric scale height.

We study electron and positron acceleration at the termination shock of a striped pulsar wind. Drift motion along the shock surface keeps either electrons or positrons -but not both, close to the equatorial plane of the pulsar, where they are accelerated by the first-order Fermi process. Their energy spectrum is a power law, and both the X-ray flux and photon index of the Crab Nebula, as measured by NuSTAR, can be reproduced for sufficiently large downstream turbulence levels. The implication that one sign of charge is preferentially accelerated in pulsar wind nebulae is potentially important for the interpretation of the positron fraction in cosmic-rays.

Two-fluid (electron-positron) plasma modelling has shown that inductive acceleration can convert Poynting flux directly into bulk kinetic energy in the relativistic flows driven by rotating magnetized neutron stars and black holes. Here, we generalize this approach by adding an ion fluid. Solutions are presented in which all particles are accelerated as the flow expands, with comparable power channeled into each of the plasma components. In an ion-dominated flow, each species reaches the limiting rigidity, according to Hillas' criterion, in a distance significantly shorter than in a lepton-dominated flow. These solutions support the hypothesis that newly born magnetars and pulsars are potential sources of ultra-high energy cosmic rays. The competing process of Poynting flux dissipation by magnetic reconnection is shown to be ineffective in low-density flows in which the conventionally defined electron multiplicity satisfies $\kappa_{\rm e}\lesssim 10^5\left(4\pi L_{38}/\Omega\right)^{1/4} /\textrm{Max}\left(\eta_{\rm ion}^{1/2},1\right)$, where $L_{38}\times 10^{38}\textrm{erg s}^{-1}$ is the power carried by the flow in a solid angle $\Omega$, and $\eta_{\rm ion}$ is the ratio of the ion to lepton power at launch.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

Herschel observations of interstellar clouds support a paradigm for star formation in which filaments play a central role. One of the foundations of this paradigm is the finding, based on studies of the radial column density profiles, that nearby filaments share a common inner width of ~0.1 pc. The existence of a characteristic filament width has been questioned, however, as it seems inconsistent with the scale-free nature of the power spectrum of cloud images. In an effort to clarify the origin of this apparent discrepancy, we examined the power spectra of the Herschel 250 micron images of the Polaris, Aquila, and Taurus clouds and performed a number of simple numerical experiments by injecting synthetic filaments in both the Herschel images and synthetic background images. We constructed several populations of synthetic filaments of 0.1 pc width with realistic area filling factors ($A_{\rm fil}$) and distributions of column density contrasts ($\delta_c$). After adding synthetic filaments to the Herschel images, we re-computed the image power spectra and compared the results with the original, essentially scale-free power spectra. We used the $\chi^2$ variance of the residuals between the best power-law fit and the output power spectrum as a diagnostic of the presence of a significant departure from a scale-free power spectrum. We found that the $\chi^2$ variance depends primarily on the combined parameter $\delta_c^2\, A_{\rm fil}$. Analysis of the real Herschel data shows that the observed filamentary structure contributes only 1/5 of the power in the power spectra at angular frequencies where an effect of the characteristic filament width is expected. We conclude that the scale-free power spectra of Herschel images remain consistent with the existence of a characteristic filament width ~0.1 pc and do not invalidate the conclusions drawn from studies of the filament profiles.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H~(0.4 - 2.2) x10$^{-6}$, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $\zeta_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We calculate the shape of the large-scale anisotropy of TeV-PeV cosmic-rays (CR) in different models of the interstellar turbulence. In general, the large-scale CR anisotropy (CRA) is not a dipole, and its shape can be used as a new probe of the turbulence. The 400 TeV and 2 PeV data sets of IceTop can be fitted with Goldreich-Sridhar turbulence and a broad resonance function, but other possibilities are not excluded. We then present our first numerical calculations of the CRA down to 3 TeV energies in 3D isotropic Kolmogorov turbulence. At these low energies, the large-scale CRA aligns well with the direction of local magnetic field lines around the observer. In this type of turbulence, the CR intensity is flat in a broad region perpendicular to field lines. Even though the CRA is quite gyrotropic, we show that the local configuration of the turbulence around the observer does result in the appearance of weak, "non-gyrotropic" small-scale anisotropies, which contain information on the local turbulence level.

[Abridged] Molecular filaments have received special attention recently, thanks to new observational results on their properties. In particular, our early analysis of filament properties revealed a narrow distribution of median widths centered at a value of about 0.1 pc. Here, we extend and complement our initial study with the analysis of the filamentary structures observed in eight nearby molecular clouds. We use the column density maps derived from Herschel data and the DisPerSE algorithm to trace a network of filaments in each cloud. We build synthetic maps to assess the completeness limit of our extracted sample and validate our measurements. Our analysis yields a selected sample of 599 filaments with aspect ratios larger than 3 and column density contrasts larger than 0.3. We show that our sample is more than 95% complete for column density contrasts larger than 1, with only 5% of spurious detections. On average, more than 15% of the total gas mass in the clouds, and more than 80% of the dense gas mass, is found to be in the form of filaments, respectively. Analysis of the radial column density profiles of the 599 filaments indicates a narrow distribution of crest-averaged inner widths, with a median value of 0.10 pc and an interquartile range of 0.07 pc. In contrast, the filaments span wide ranges in length, central column density, column density contrast, and mass per unit length. The characteristic filament width is well resolved by Herschel observations, and a median value of 0.1 pc is consistently found using three distinct estimates based on (1) a direct measurement of the width at half power after background subtraction, as well as (2) Gaussian and (3) Plummer fits. The existence of a characteristic filament width is further supported by the presence of a tight correlation between mass per unit length and central column density for the observed filaments.

We report the discovery of the super-Earth K2-265 b detected with K2 photometry. The planet orbits a bright (V_mag = 11.1) star of spectral type G8V with a period of 2.37 days. We obtained high-precision follow-up radial velocity measurements from HARPS, and the joint Bayesian analysis showed that K2-265 b has a radius of 1.71 +/- 0.11 R_earth and a mass of 6.54 +/- 0.84 M_earth, corresponding to a bulk density of 7.1 +/- 1.8 g/cm^3 . Composition analysis of the planet reveals an Earth-like, rocky interior, with a rock mass fraction of 80%. The short orbital period and small radius of the planet puts it below the lower limit of the photoevaporation gap, where the envelope of the planet could have eroded due to strong stellar irradiation, leaving behind an exposed core. Knowledge of the planet core composition allows us to infer the possible formation and evolution mechanism responsible for its current physical parameters.

We study particle acceleration at the termination shock of a striped pulsar wind by integrating trajectories in a prescribed model of the magnetic field and flow pattern. Drift motion on the shock surface maintains either electrons or positrons on "Speiser" orbits in a ring-shaped region close to the equatorial plane of the pulsar, enabling them to be accelerated to very high energy by the first-order Fermi mechanism. A power-law spectrum results: $dN_{\rm e}/d\gamma\propto\gamma^{\alpha_{\rm e}}$, where $\alpha_{\rm e}$ lies in the range $-1.8$ to $-2.4$ and depends on the downstream turbulence level. For sufficiently strong turbulence, we find $\alpha_{\rm e} \simeq -2.2$, and both the photon index and the flux of $1-100$ keV X-rays from the Crab Nebula, as measured by NuSTAR, can be reproduced. The particle spectrum hardens to $\alpha_{\rm e} \simeq -1.8$ at lower turbulence levels, which may explain the hard photon index observed by the Chandra X-ray Observatory in the central regions of the Nebula.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

The Earth, Venus, Mars, and some extrasolar terrestrial planets have a mass and radius that is consistent with a mass fraction of about 30% metallic core and 70% silicate mantle. At the inner frontier of the solar system, Mercury has a completely different composition, with a mass fraction of about 70% metallic core and 30% silicate mantle. Several formation or evolution scenarios are proposed to explain this metal-rich composition, such as a giant impact, mantle evaporation, or the depletion of silicate at the inner-edge of the proto-planetary disk. These scenarios are still strongly debated. Here we report the discovery of a multiple transiting planetary system (K2-229), in which the inner planet has a radius of 1.165+/-0.066 Rearth and a mass of 2.59+/-0.43 Mearth. This Earth-sized planet thus has a core-mass fraction that is compatible with that of Mercury, while it was expected to be similar to that of the Earth based on host-star chemistry. This larger Mercury analogue either formed with a very peculiar composition or it has evolved since, e.g. by losing part of its mantle. Further characterisation of Mercury-like exoplanets like K2-229 b will help putting the detailed in-situ observations of Mercury (with Messenger and BepiColombo) into the global context of the formation and evolution of solar and extrasolar terrestrial planets.

Planets that reside close-in to their host star are subject to intense high-energy irradiation. Extreme-ultraviolet (EUV) and X-ray radiation (together, XUV) is thought to drive mass loss from planets with volatile envelopes. We present $\textit{XMM-Newton}$ observations of six nearby stars hosting transiting planets in tight orbits (with orbital period, $P_\text{orb} < 10\,$d), wherein we characterise the XUV emission from the stars and subsequent irradiation levels at the planets. In order to reconstruct the unobservable EUV emission, we derive a new set of relations from Solar $\textit{TIMED/SEE}$ data that are applicable to the standard bands of the current generation of X-ray instruments. From our sample, WASP-80b and HD$\,$149026b experience the highest irradiation level, but HAT-P-11b is probably the best candidate for Ly$\,\alpha$ evaporation investigations because of the system's proximity to the Solar System. The four smallest planets have likely lost a greater percentage of their mass over their lives than their larger counterparts. We also detect the transit of WASP-80b in the near ultraviolet with the Optical Monitor on $\textit{XMM-Newton}$.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

We present a catalogue of prestellar and starless cores within the Corona Australis molecular cloud using photometric data from the Herschel Space Observatory. At a distance of d~130 pc, Corona Australis is one of the closest star-forming regions. Herschel has taken multi-wavelength data of Corona Australis with both the SPIRE and PACS photometric cameras in a parallel mode with wavelengths in the range 70 {\mu}m to 500 {\mu}m. A complete sample of starless and prestellar cores and embedded protostars is identified. Other results from the Herschel Gould Belt Survey have shown spatial correlation between the distribution of dense cores and the filamentary structure within the molecular clouds. We go further and show correlations between the properties of these cores and their spatial distribution within the clouds, with a particular focus on the mass distribution of the dense cores with respect to their filamentary proximity. We find that only lower-mass starless cores form away from filaments, while all of the higher-mass prestellar cores form in close proximity to, or directly on the filamentary structure. This result supports the paradigm that prestellar cores mostly form on filaments. We analyse the mass distribution across the molecular cloud, finding evidence that the region around the Coronet appears to be at a more dynamically advanced evolutionary stage to the rest of the clumps within the cloud.

We present far-infrared observations of Monoceros R2 (a giant molecular cloud at approximately 830 pc distance, containing several sites of active star formation), as observed at 70 {\mu}m, 160 {\mu}m, 250 {\mu}m, 350 {\mu}m, and 500 {\mu}m by the Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver (SPIRE) instruments on the Herschel Space Observatory as part of the Herschel imaging survey of OB young stellar objects (HOBYS) Key programme. The Herschel data are complemented by SCUBA-2 data in the submillimetre range, and WISE and Spitzer data in the mid-infrared. In addition, C18O data from the IRAM 30-m Telescope are presented, and used for kinematic information. Sources were extracted from the maps with getsources, and from the fluxes measured, spectral energy distributions were constructed, allowing measurements of source mass and dust temperature. Of 177 Herschel sources robustly detected in the region (a detection with high signal-to-noise and low axis ratio at multiple wavelengths), including protostars and starless cores, 29 are found in a filamentary hub at the centre of the region (a little over 1% of the observed area). These objects are on average smaller, more massive, and more luminous than those in the surrounding regions (which together suggest that they are at a later stage of evolution), a result that cannot be explained entirely by selection effects. These results suggest a picture in which the hub may have begun star formation at a point significantly earlier than the outer regions, possibly forming as a result of feedback from earlier star formation. Furthermore, the hub may be sustaining its star formation by accreting material from the surrounding filaments.

The ~400 MeV flaring emission from the Crab Nebula is naturally explained as the result of an abrupt reduction in the mass-loading of the pulsar wind. Very few particles are then available to carry the current required to maintain wave activity, causing them to achieve high Lorentz factors. When they penetrate the Nebula, a tightly beamed, high luminosity burst of hard gamma-rays results, with characteristics similar to the observed flares. This mechanism may operate in other powerful pulsars, such as J0537-6910 (PWN N 157B), B0540-69, B1957+20 and J0205+6449 (3C 58).

We have performed ground-based transmission spectroscopy of the hot Jupiter orbiting the cool dwarf WASP-80 using the ACAM instrument on the William Herschel Telescope (WHT) as part of the LRG-BEASTS programme. This is the third paper of a ground-based transmission spectroscopy survey of hot Jupiters using low-resolution grism spectrographs. We observed two transits of the planet and have constructed transmission spectra spanning a wavelength range of 4640-8840A. Our transmission spectrum is inconsistent with a previously claimed detection of potassium in WASP-80b's atmosphere, and is instead most consistent with a haze. We also do not see evidence for sodium absorption at a resolution of 100A.

We show that the mysterious, rapidly variable emission at ~400 MeV observed from the Crab Nebula by the AGILE and Fermi experiments could be the result of a sudden drop in the mass-loading of the pulsar wind. The current required to maintain wave activity in the wind is then carried by very few particles of high Lorentz factor. On impacting the Nebula, these particles produce a tightly beamed, high luminosity burst of hard gamma-rays, similar to those observed. This implies (i) the emission is synchrotron radiation in the toroidal field of the Nebula, and, therefore, linearly polarized and (ii) this mechanism potentially contributes to the gamma-ray emission from other powerful pulsars, such as the Magellanic Cloud objects J0537-6910 and B0540-69.

We determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-Forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of $B_{\rm pos}=6.6\pm4.7$ mG, where $\delta B_{\rm pos}=4.7$ mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of $\sim1.7\times 10^{-7}$ Jm$^{-3}$ in OMC 1, comparable both to the gravitational potential energy density of OMC 1 ($\sim 10^{-7}$ Jm$^{-3}$), and to the energy density in the Orion BN/KL outflow ($\sim 10^{-7}$ Jm$^{-3}$). We find that neither the Alfv\'{e}n velocity in OMC 1 nor the velocity of the super-Alfv\'{e}nic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the $\sim$500-year lifetime of the outflow. Hence, we propose that the hour-glass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically-symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically-supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.

We report the discovery of the exoplanet K2-110 b (previously EPIC212521166b) from K2 photometry orbiting in a 13.8637d period around an old, metal-poor K3 dwarf star. With a V-band magnitude of 11.9, K2-110 is particularly amenable to RV follow-up. A joint analysis of K2 photometry and high-precision RVs from 28 HARPS and HARPS-N spectra reveal it to have a radius of 2.6$\pm 0.1 R_{\oplus}$ and a mass of 16.7$\pm 3.2$~M$_{\oplus}$, hence a density of $5.2\pm1.2$ g.cm$^{-3}$, making it one of the most massive planets yet to be found with a sub-Neptune radius. When accounting for compression, the resulting Earth-like density is best fitted by a 0.2 M$_{\oplus}$ hydrogen atmosphere over an 16.5 M$_{\oplus}$ Earth-like interior, although the planet could also have significant water content. At 0.1~AU, even taking into account the old stellar age of $8 \pm 3$ Gyr, the planet is unlikely to have been significantly affected by EUV evaporation. However the planet likely disc-migrated to its current position making the lack of a thick H$_2$ atmosphere puzzling. This analysis has made K2-110 b one of the best-characterised mini-Neptunes with density constrained to less than 30%.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We have performed ground-based transmission spectroscopy of the hot Jupiter HAT-P-18b using the ACAM instrument on the William Herschel Telescope (WHT). Differential spectroscopy over an entire night was carried out at a resolution of $R \approx 400$ using a nearby comparison star. We detect a bluewards slope extending across our optical transmission spectrum which runs from 4750 to 9250\AA. The slope is consistent with Rayleigh scattering at the equilibrium temperature of the planet (852K). We do not detect enhanced sodium absorption, which indicates that a high-altitude haze is masking the feature and giving rise to the Rayleigh slope. This is only the second discovery of a Rayleigh scattering slope in a hot Jupiter atmosphere from the ground, and our study illustrates how ground-based observations can provide transmission spectra with precision comparable to the Hubble Space Telescope.

To constrain models of high-mass star formation, the Herschel/HOBYS KP aims at discovering massive dense cores (MDCs) able to host the high-mass analogs of low-mass prestellar cores, which have been searched for over the past decade. We here focus on NGC6334, one of the best-studied HOBYS molecular cloud complexes. We used Herschel PACS and SPIRE 70-500mu images of the NGC6334 complex complemented with (sub)millimeter and mid-infrared data. We built a complete procedure to extract ~0.1 pc dense cores with the getsources software, which simultaneously measures their far-infrared to millimeter fluxes. We carefully estimated the temperatures and masses of these dense cores from their SEDs. A cross-correlation with high-mass star formation signposts suggests a mass threshold of 75Msun for MDCs in NGC6334. MDCs have temperatures of 9.5-40K, masses of 75-1000Msun, and densities of 10^5-10^8cm-3. Their mid-IR emission is used to separate 6 IR-bright and 10 IR-quiet protostellar MDCs while their 70mu emission strength, with respect to fitted SEDs, helps identify 16 starless MDC candidates. The ability of the latter to host high-mass prestellar cores is investigated here and remains questionable. An increase in mass and density from the starless to the IR-quiet and IR-bright phases suggests that the protostars and MDCs simultaneously grow in mass. The statistical lifetimes of the high-mass prestellar and protostellar core phases, estimated to be 1-7x10^4yr and at most 3x10^5yr respectively, suggest a dynamical scenario of high-mass star formation. The present study provides good mass estimates for a statistically significant sample, covering the earliest phases of high-mass star formation. High-mass prestellar cores may not exist in NGC6334, favoring a scenario presented here, which simultaneously forms clouds and high-mass protostars.

We have measured a precise optical transmission spectrum for WASP-52b, a highly inflated hot Jupiter with an equilibrium temperature of 1300 K. Two transits of the planet were observed spectroscopically at low resolution with the auxiliary-port camera (ACAM) on the William Herschel Telescope (WHT), covering a wide range of 4000-8750 \AA. We use a Gaussian process approach to model the correlated noise in the multi-wavelength light curves, resulting in a high precision relative transmission spectrum with errors on the order of a pressure scale height. We attempted to fit a variety of different representative model atmospheres to the transmission spectrum, but did not find a satisfactory match to the entire spectral range. For the majority of the covered wavelength range (4000-7750 \AA) the spectrum is flat, and can be explained by an optically thick and grey cloud layer at 0.1 mbar, but this is inconsistent with a slightly deeper transit at wavelengths $> 7750$ \AA. We were not able to find an obvious systematic source for this feature, so this opacity may be the result of an additional unknown absorber.

We calculate the large-scale cosmic-ray (CR) anisotropies predicted for a range of Goldreich-Sridhar (GS) and isotropic models of interstellar turbulence, and compare them with IceTop data. In general, the predicted CR anisotropy is not a pure dipole; the cold spots reported at 400 TeV and 2 PeV are consistent with a GS model that contains a smooth deficit of parallel-propagating waves and a broad resonance function, though some other possibilities cannot, as yet, be ruled out. In particular, isotropic fast magnetosonic wave turbulence can match the observations at high energy, but cannot accommodate an energy dependence in the shape of the CR anisotropy. Our findings suggest that improved data on the large-scale CR anisotropy could provide a valuable probe of the properties - notably the power-spectrum - of the interstellar turbulence within a few tens of parsecs from Earth.

We study the role of local interstellar turbulence in shaping the large-scale anisotropy in the arrival directions of TeV-PeV cosmic-rays (CRs) on the sky. Assuming pitch-angle diffusion of CRs in a magnetic flux tube containing the Earth, we compute the CR anisotropy for Goldreich-Sridhar turbulence, and for isotropic fast modes. The narrow deficits in the 400 TeV and 2 PeV data sets of IceTop can be fitted for some parameters of the turbulence. The data also rule out a part of the parameter space. The shape of the CR anisotropy may be used as a local probe of the still poorly known properties of the interstellar turbulence and of CR transport.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

We study the acceleration of electrons and positrons at an electromagnetically modified, ultra-relativistic shock in the context of pulsar wind nebulae (PWNe). We simulate the outflow produced by an obliquely rotating pulsar in proximity of its termination shock with a two-fluid code which uses a magnetic shear wave to mimic the properties of the wind. We integrate electron trajectories in the test-particle limit in the resulting background electromagnetic fields to analyse the injection mechanism. We find that the shock-precursor structure energizes and reflects a sizeable fraction of particles, which becomes available for further acceleration. We investigate the subsequent first-order Fermi process sustained by small-scale magnetic fluctuations with a Monte Carlo code. We find that the acceleration proceeds in two distinct regimes: when the gyro-radius $r_{\textrm{g}}$ exceeds the wavelength of the shear $\lambda$, the process is remarkably similar to first-order Fermi acceleration at relativistic, parallel shocks. This regime corresponds to a low density wind which allows the propagation of superluminal waves. When $r_{\textrm{g}}<\lambda$, which corresponds to the scenario of driven reconnection, the spectrum is softer.

We report three newly discovered exoplanets from the SuperWASP survey. WASP-127b is a heavily inflated super-Neptune of mass 0.18 +/- 0.02 M_J and radius 1.37 +/- 0.04 R_J. This is one of the least massive planets discovered by the WASP project. It orbits a bright host star (Vmag = 10.16) of spectral type G5 with a period of 4.17 days. WASP-127b is a low-density planet that has an extended atmosphere with a scale height of 2500 +/- 400 km, making it an ideal candidate for transmission spectroscopy. WASP-136b and WASP-138b are both hot Jupiters with mass and radii of 1.51 +/- 0.08 M_J and 1.38 +/- 0.16 R_J, and 1.22 +/- 0.08 M_J and 1.09 +/- 0.05 R_J, respectively. WASP-136b is in a 5.22-day orbit around an F9 subgiant star with a mass of 1.41 +/- 0.07 M_sun and a radius of 2.21 +/- 0.22 R_sun. The discovery of WASP-136b could help constrain the characteristics of the giant planet population around evolved stars. WASP-138b orbits an F7 star with a period of 3.63 days. Its radius agrees with theoretical values from standard models, suggesting the presence of a heavy element core with a mass of ~10 M_earth. The discovery of these new planets helps in exploring the diverse compositional range of short-period planets, and will aid our understanding of the physical characteristics of both gas giants and low-density planets.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

We have measured the transmission spectrum of the extremely inflated hot Jupiter WASP-52b using simultaneous photometric observations in SDSS u', g' and a filter centred on the sodium doublet (NaI) with the ULTRACAM instrument mounted on the 4.2m William Herschel Telescope. We find that Rayleigh scattering is not the dominant source of opacity within the planetary atmosphere and find a transmission spectrum more consistent with wavelength-independent opacity such as from clouds. We detect an in-transit anomaly that we attribute to the presence of stellar activity and find that this feature can be more simply modelled as a bright region on the stellar surface akin to Solar faculae rather than spots. A spot model requires a significantly larger planet/star radius ratio than that found in previous studies. Our results highlight the precision that can be achieved by ground-based photometry with errors in the scaled planetary radii of less than one atmospheric scale height, comparable to HST observations.

We present a comparison of SCUBA-2 850-$\mu$m and Herschel 70--500-$\mu$m observations of the L1495 filament in the Taurus Molecular Cloud with the goal of characterising the SCUBA-2 Gould Belt Survey (GBS) data set. We identify and characterise starless cores in three data sets: SCUBA-2 850-$\mu$m, Herschel 250-$\mu$m, and Herschel 250-$\mu$m spatially filtered to mimic the SCUBA-2 data. SCUBA-2 detects only the highest-surface-brightness sources, principally detecting protostellar sources and starless cores embedded in filaments, while Herschel is sensitive to most of the cloud structure, including extended low-surface-brightness emission. Herschel detects considerably more sources than SCUBA-2 even after spatial filtering. We investigate which properties of a starless core detected by Herschel determine its detectability by SCUBA-2, and find that they are the core's temperature and column density (for given dust properties). For similar-temperature cores, such as those seen in L1495, the surface brightnesses of the cores are determined by their column densities, with the highest-column-density cores being detected by SCUBA-2. For roughly spherical geometries, column density corresponds to volume density, and so SCUBA-2 selects the densest cores from a population at a given temperature. This selection effect, which we quantify as a function of distance, makes SCUBA-2 ideal for identifying those cores in Herschel catalogues that are closest to forming stars. Our results can now be used by anyone wishing to use the SCUBA-2 GBS data set.

Multi-colour broad-band transit observations offer the opportunity to characterise the atmosphere of an extrasolar planet with small- to medium-sized telescopes. One of the most favourable targets is the hot Jupiter HAT-P-32 b. We combined 21 new transit observations of this planet with 36 previously published light curves for a homogeneous analysis of the broad-band transmission spectrum from the Sloan u' band to the Sloan z' band. Our results rule out cloud-free planetary atmosphere models of solar metallicity. Furthermore, a discrepancy at reddest wavelengths to previously published results makes a recent tentative detection of a scattering feature less likely. Instead, the available spectral measurements of HAT-P-32 b favour a completely flat spectrum from the near-UV to the near-IR. A plausible interpretation is a thick cloud cover at high altitudes.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

We report the discovery of the two hot-Jupiters K2-30b and K2-34b. The two planets were detected transiting their main-sequence star with periods $\sim$ 4.099 and $\sim$ 2.996 days, in campaigns 4 and 5 of the extension of the Kepler mission, K2. Subsequent ground-based radial velocity follow-up with SOPHIE, HARPS-N and CAFE, established the planetary nature of the transiting objects. We analyzed the transit signal, radial velocity and spectral energy distributions of the two systems to characterize their properties. Both planets (K2-30b and K2-34b) are bloated hot-Jupiters (1.25 $R_{\rm Jup}$ and 1.33 $R_{\rm Jup}$) around relatively bright (V =13.5 and V=11.5), slow rotating main-sequence (G8 and F9) stars. Thus, these systems are good candidates for detecting the Rossiter-MacLaughlin effect to measure their obliquity and for atmospheric studies.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

We observed high S/N, high velocity resolution NH$_3$(1,1) and (2,2) emission on an extended map in TMC-1, a filamentary cloud in a nearby quiescent star forming area. By fitting multiple hyperfine-split line profiles to the NH$_3$(1,1) spectra we derived the velocity distribution of the line components and calculated gas parameters on several positions. Herschel SPIRE continuum observations were reduced and used to calculate the physical parameters of the Planck Galactic Cold Clumps in the region. The Herschel-based column density map of TMC-1 shows a main ridge with two local maxima and a separated peak to the south-west. H$_2$-column densities and dust temperatures are in the range of 0.5-3.3 $\times$ 10$^{22}$ cm$^{-2}$ and 10.5-12 K, respectively. NH$_3$-column densities are 2.8-14.2 $\times$ 10$^{14}$ cm$^{-2}$ and and H$_2$-volume densities are 0.4-2.8 $\times$ 10$^4$ cm$^{-3}$. Kinetic temperatures are typically very low with a minimum of 9 K, and a maximum of 13.7 K was found at the Class I protostar IRAS 04381+2540. The kinetic temperatures vary similarly as the dust temperatures in spite of the fact that densities are lower than the critical density for coupling between the gas and dust phase. The k-means clustering method separated four sub-filaments in TMC-1 in the position-velocity-column density parameter space. They have masses of 32.5, 19.6, 28.9 and 45.9 M$_{\odot}$, low turbulent velocity dispersion (0.13-0.2 kms$^{-1}$) and they are close to gravitational equilibrium. We label them TMC-1F1 through F4. TMC-1F1, TMC-1F2 and TMC-1F4 are very elongated, dense and cold. TMC-1F3 is a little less elongated and somewhat warmer, probably heated by IRAS 04381+2540 that is embedded in it. TMC-1F3 is $\approx$ 0.1 pc behind TMC1-F1. Because of its structure, TMC-1 is a good target to test filament evolution scenarios.

We present a catalogue of dense cores in a $\sim 4^\circ\times2^\circ$ field of the Taurus star-forming region, inclusive of the L1495 cloud, derived from Herschel SPIRE and PACS observations in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, and 500 $\mu$m continuum bands. Estimates of mean dust temperature and total mass are derived using modified blackbody fits to the spectral energy distributions. We detect 525 starless cores of which $\sim10$-20% are gravitationally bound and therefore presumably prestellar. Our census of unbound objects is $\sim85$% complete for $M>0.015\,M_\odot$ in low density regions ($A_V\stackrel{<}{_\sim}5$ mag), while the bound (prestellar) subset is $\sim85$% complete for $M>0.1\,M_\odot$ overall. The prestellar core mass function (CMF) is consistent with lognormal form, resembling the stellar system initial mass function, as has been reported previously. All of the inferred prestellar cores lie on filamentary structures whose column densities exceed the expected threshold for filamentary collapse, in agreement with previous reports. Unlike the prestellar CMF, the unbound starless CMF is not lognormal, but instead is consistent with a power-law form below $0.3\,M_\odot$ and shows no evidence for a low-mass turnover. It resembles previously reported mass distributions for CO clumps at low masses ($M\stackrel{<}{_\sim}0.3\,M_\odot$). The volume density PDF, however, is accurately lognormal except at high densities. It is consistent with the effects of self-gravity on magnetized supersonic turbulence. The only significant deviation from lognormality is a high-density tail which can be attributed unambiguously to prestellar cores.

Dust grains play a key role in the physics of star-forming regions, even though they constitute only $\sim$1 % of the mass of the interstellar medium. The derivation of accurate dust parameters such as temperature ($T_{dust}$), emissivity spectral index ($\beta$) and column density requires broadband continuum observations at far-infrared wavelengths. We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) measurements of three starless cores: L1521E, L1521F and L1689B, covering wavelengths between 194 and 671 $\mu$m. This paper is the first to use our recently updated SPIRE-FTS intensity calibration, yielding a direct match with SPIRE photometer measurements of extended sources. In addition, we carefully assess the validity of calibration schemes depending on source extent and on the strength of background emission. The broadband far-infrared spectra for all three sources peak near 250 $\mu$m. Our observations therefore provide much tighter constraints on the spectral energy distribution (SED) shape than measurements that do not probe the SED peak. The spectra are fitted using modified blackbody functions, allowing both $T_{dust}$ and $\beta$ to vary as free parameters. This yields $T_{dust}$ of 9.8$\pm$0.2 K, 15.6$\pm$0.5 K and 10.9$\pm$0.2 K and corresponding $\beta$ of 2.6$\mp$0.9, 0.8$\mp$0.1 and 2.4$\mp$0.8 for L1521E, L1521F and L1689B respectively. The derived core masses are 1.0$\pm$0.1, 0.10$\pm$0.01 and 0.49$\pm$0.05 $M_{\odot}$, respectively. The core mass/Jeans mass ratios for L1521E and L1689B exceed unity indicating that they are unstable to gravitational collapse, and thus pre-stellar cores. By comparison, the elevated temperature and gravitational stability of L1521F support previous arguments that this source is more evolved and likely a protostar.

In the present paper we report the discovery of a new hot Jupiter, EPIC211089792 b, first detected by the Super-WASP observatory and then by the K2 space mission during its campaign 4. The planet has a period of 3.25d, a mass of 0.73 +/- 0.04 Mjup, and a radius of 1.19 +/- 0.02 Rjup. The host star is a relatively bright (V=12.5) G7 dwarf with a nearby K5V companion. Based on stellar rotation and the abundance of Lithium, we find that the system might be as young as about 450 Myr. The observation of the Rossiter-McLaughlin effect shows the planet is aligned with respect to the stellar spin. Given the deep transit (20mmag), the magnitude of the star and the presence of a nearby stellar companion, the planet is a good target for both space- and ground-based transmission spectroscopy, in particular in the near-infrared where the both stars are relatively bright.

We are entering an era of unprecedented quantities of data from current and planned survey telescopes. To maximise the potential of such surveys, automated data analysis techniques are required. Here we implement a new methodology for variable star classification, through the combination of Kohonen Self Organising Maps (SOM, an unsupervised machine learning algorithm) and the more common Random Forest (RF) supervised machine learning technique. We apply this method to data from the K2 mission fields 0-4, finding 154 ab-type RR Lyraes (10 newly discovered), 377 Delta Scuti pulsators, 133 Gamma Doradus pulsators, 183 detached eclipsing binaries, 290 semi-detached or contact eclipsing binaries and 9399 other periodic (mostly spot-modulated) sources, once class significance cuts are taken into account. We present lightcurve features for all K2 stellar targets, including their three strongest detected frequencies, which can be used to study stellar rotation periods where the observed variability arises from spot modulation. The resulting catalogue of variable stars, classes, and associated data features are made available online. We publish our SOM code in Python as part of the open source PyMVPA package, which in combination with already available RF modules can be easily used to recreate the method.

Photometric surveys such as Kepler have the precision to identify exoplanet and eclipsing binary candidates from only a single transit. K2, with its 75d campaign duration, is ideally suited to detect significant numbers of single-eclipsing objects. Here we develop a Bayesian transit-fitting tool ("Namaste: An Mcmc Analysis of Single Transit Exoplanets") to extract orbital information from single transit events. We achieve favourable results testing this technique on known Kepler planets, and apply the technique to 7 candidates identified from a targeted search of K2 campaigns 1, 2 and 3. We find EPIC203311200 to host an excellent exoplanet candidate with a period, assuming zero eccentricity, of $540 ^{+410}_{-230}$ days and a radius of $0.51 \pm 0.05 R_{Jup}$. We also find six further transit candidates for which more follow-up is required to determine a planetary origin. Such a technique could be used in the future with TESS, PLATO and ground-based photometric surveys such as NGTS, potentially allowing the detection of planets in reach of confirmation by Gaia.

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

We report the discovery of a new transiting planet from the WASP survey. WASP-135b is a hot Jupiter with a radius of 1.30 pm 0.09 Rjup, a mass of 1.90 pm 0.08 Mjup and an orbital period of 1.401 days. Its host is a Sun-like star, with a G5 spectral type and a mass and radius of 0.98 pm 0.06 Msun and 0.96 pm 0.05 Rsun respectively. The proximity of the planet to its host means that WASP-135b receives high levels of insolation, which may be the cause of its inflated radius. Additionally, we find weak evidence of a transfer of angular momentum from the planet to its star.

K2-19 is the second multi-planetary system discovered with K2 observations. The system is composed of two Neptune size planets close to the 3:2 mean-motion resonance. To better characterise the system we obtained two additional transit observations of K2-19b and five additional radial velocity observations. These were combined with K2 data and fitted simultaneously with the system dynamics (photo-dynamical model) which increases the precision of the transit time measurements. The higher transit time precision allows us to detect the chopping signal of the dynamic interaction of the planets that in turn permits to uniquely characterise the system. Although the reflex motion of the star was not detected, dynamic modelling of the system allowed us to derive planetary masses of $M_b= 44 \pm 12\, M_{\oplus}$ and $M_c = 15.9 \pm 7.0\, M_{\oplus}$ for the inner and the outer planets respectively, leading to densities close to Uranus. We also show that our method allows the derivation of mass ratios using only the 80 days of observations during the first campaign of K2.

The K2 mission has recently begun to discover new and diverse planetary systems. In December 2014 Campaign 1 data from the mission was released, providing high-precision photometry for ~22000 objects over an 80 day timespan. We searched these data with the aim of detecting further important new objects. Our search through two separate pipelines led to the independent discovery of K2-19b \& c, a two-planet system of Neptune sized objects (4.2 and 7.2 $R_\oplus$), orbiting a K dwarf extremely close to the 3:2 mean motion resonance. The two planets each show transits, sometimes simultaneously due to their proximity to resonance and alignment of conjunctions. We obtain further ground based photometry of the larger planet with the NITES telescope, demonstrating the presence of large transit timing variations (TTVs), and use the observed TTVs to place mass constraints on the transiting objects under the hypothesis that the objects are near but not in resonance. We then statistically validate the planets through the \texttt{PASTIS} tool, independently of the TTV analysis.

We present the full catalog of Young Stellar Objects (YSOs) identified in the 18 molecular clouds surveyed by the Spitzer Space Telescope "cores to disks" (c2d) and "Gould Belt" (GB) Legacy surveys. Using standard techniques developed by the c2d project, we identify 3239 candidate YSOs in the 18 clouds, 2966 of which survive visual inspection and form our final catalog of YSOs in the Gould Belt. We compile extinction corrected SEDs for all 2966 YSOs and calculate and tabulate the infrared spectral index, bolometric luminosity, and bolometric temperature for each object. We find that 326 (11%), 210 (7%), 1248 (42%), and 1182 (40%) are classified as Class 0+I, Flat-spectrum, Class II, and Class III, respectively, and show that the Class III sample suffers from an overall contamination rate by background AGB stars between 25% and 90%. Adopting standard assumptions, we derive durations of 0.40-0.78 Myr for Class 0+I YSOs and 0.26-0.50 Myr for Flat-spectrum YSOs, where the ranges encompass uncertainties in the adopted assumptions. Including information from (sub)millimeter wavelengths, one-third of the Class 0+I sample is classified as Class 0, leading to durations of 0.13-0.26 Myr (Class 0) and 0.27-0.52 Myr (Class I). We revisit infrared color-color diagrams used in the literature to classify YSOs and propose minor revisions to classification boundaries in these diagrams. Finally, we show that the bolometric temperature is a poor discriminator between Class II and Class III YSOs.

We present and discuss the results of the Herschel Gould Belt survey observations in a ~11 deg^2 area of the Aquila molecular cloud complex at d~260 pc, imaged with the SPIRE/PACS cameras from 70 to 500 micron. We identify a complete sample of starless dense cores and embedded protostars in this region, and analyze their global properties and spatial distributions. We find a total of 651 starless cores, ~60% of which are gravitationally bound prestellar cores, and they will likely form stars in the future. We also detect 58 protostellar cores. The core mass function (CMF) derived for the prestellar cores is very similar in shape to the stellar initial mass function (IMF), supporting the earlier view that there is a close physical link between the IMF and the CMF. The global shift in mass scale observed between the CMF and the IMF is consistent with a typical star formation efficiency of ~40%. By comparing the numbers of starless cores to the number of young stellar objects, we estimate that the lifetime of prestellar cores is ~1 Myr. We find a strong correlation between the spatial distribution of prestellar cores and the densest filaments. About 90% of the Herschel-identified prestellar cores are located above a background column density corresponding to A_V~7, and ~75% of them lie within filamentary structures with supercritical masses per unit length >~16 M_sun/pc. These findings support a picture wherein the cores making up the peak of the CMF (and probably responsible for the base of the IMF) result primarily from the gravitational fragmentation of marginally supercritical filaments. Given that filaments appear to dominate the mass budget of dense gas at A_V>7, our findings also suggest that the physics of prestellar core formation within filaments is responsible for a characteristic "efficiency" SFR/M_dense ~5+-2 x 10^-8 yr^-1 for the star formation process in dense gas.

Perpendicular shocks are shown to be rapid particle accelerators that perform optimally when the ratio $u_{\rm s}$ of the shock speed to the particle speed roughly equals the ratio $1/\eta$ of the scattering rate to the gyro frequency. We use analytical methods and Monte-Carlo simulations to solve the kinetic equation that governs the anisotropy generated at these shocks, and find, for $\eta u_{\rm s}\approx1$, that the spectral index softens by unity and the acceleration time increases by a factor of two compared to the standard result of diffusive shock acceleration theory. These results provide a theoretical basis for the thirty-year-old conjecture that a supernova exploding into the wind of a Wolf-Rayet star may accelerate protons to an energy exceeding $10^{15}\,$eV.

We have created a catalogue of variable stars found from a search of the publicly available K2 mission data from Campaigns 1 and 0. This catalogue provides the identifiers of 8395 variable stars, including 199 candidate eclipsing binaries with periods up to 60d and 3871 periodic or quasi-periodic objects, with periods up to 20d for Campaign 1 and 15d for Campaign 0. Lightcurves are extracted and detrended from the available data. These are searched using a combination of algorithmic and human classification, leading to a classifier for each object as an eclipsing binary, sinusoidal periodic, quasi periodic, or aperiodic variable. The source of the variability is not identified, but could arise in the non-eclipsing binary cases from pulsation or stellar activity. Each object is cross-matched against variable star related guest observer proposals to the K2 mission, which specifies the variable type in some cases. The detrended lightcurves are also compared to lightcurves currently publicly available. The resulting catalogue is made available online via the MAST archive at https://archive.stsci.edu/prepds/k2varcat/, and gives the ID, type, period, semi-amplitude and range of the variation seen. We also make available the detrended lightcurves for each object.

This is a catalogue of approximately 70 X-ray emitting binary systems in the Small Magellanic Cloud (SMC) that contain a Be star as the mass donor in the system and a clear X-ray pulse signature from a neutron star. The systems are generally referred to as Be/Xray binaries. It lists all their known binary characteristics (orbital period, eccentricity), the measured spin period of the compact object, plus the characteristics of the Be star (spectral type, size of the circumstellar disk, evidence for NRP behaviour). For the first time data from the Spitzer Observatory are combined with ground-based data to provide a view of these systems out into the far-IR. Many of the observational parameters are presented as statistical distributions and compared to other similar similar populations (eg isolated Be & B stars) in the SMC, and to other Be/X-ray systems in the Milky Way. In addition previous important results are re-investigated using this excellently homogeneous sample. In particular, the evidence for a bi-modality in the spin period distribution is shown to be even stronger than first proposed, and the correlation between orbital period and circumstellar disk size seen in galactic sources is shown to be clearly present in the SMC systems and quantised for the first time.

We show Akari data, Herschel data and data from the SCUBA2 camera on JCMT, of molecular clouds. We focus on pre-stellar cores within the clouds. We present Akari data of the L1147-1157 ring in Cepheus and show how the data indicate that the cores are being externally heated. We present SCUBA2 and Herschel data of the Ophiuchus region and show how the environment is also affecting core evolution in this region. We discuss the effects of the magnetic field in the Lupus I region, and how this lends support to a model for the formation and evolution of cores in filamentary molecular clouds.

We present deep NH$_3$ observations of the L1495-B218 filaments in the Taurus molecular cloud covering over a 3 degree angular range using the K-band focal plane array on the 100m Green Bank Telescope. The L1495-B218 filaments form an interconnected, nearby, large complex extending over 8 pc. We observed NH$_3$ (1,1) and (2,2) with a spectral resolution of 0.038 km/s and a spatial resolution of 31$"$. Most of the ammonia peaks coincide with intensity peaks in dust continuum maps at 350 $\mu$m and 500 $\mu$m. We deduced physical properties by fitting a model to the observed spectra. We find gas kinetic temperatures of 8 $-$ 15 K, velocity dispersions of 0.05 $-$ 0.25 km/s, and NH$_3$ column densities of 5$\times$10$^{12}$ $-$ 1$\times$10$^{14}$ cm$^{-2}$. The CSAR algorithm, which is a hybrid of seeded-watershed and binary dendrogram algorithms, identifies a total of 55 NH$_3$ structures including 39 leaves and 16 branches. The masses of the NH$_3$ sources range from 0.05 M$_\odot$ to 9.5 M$_\odot$. The masses of NH$_3$ leaves are mostly smaller than their corresponding virial mass estimated from their internal and gravitational energies, which suggests these leaves are gravitationally unbound structures. 9 out of 39 NH$_3$ leaves are gravitationally bound and 7 out of 9 gravitationally bound NH$_3$ leaves are associated with star formation. We also found that 12 out of 30 gravitationally unbound leaves are pressure-confined. Our data suggest that a dense core may form as a pressure-confined structure, evolve to a gravitationally bound core, and undergo collapse to form a protostar.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

We have searched the K2 campaign 0 data for lightcurve variations associated with stellar variability. The results of this search are presented as a catalogue, giving the identifiers of nearly 2500 variable stars in the dataset. We list the detected range of the variation, periodicity if relevant and semi-amplitude. Lightcurves are classified into strictly periodic, quasi-periodic and aperiodic groups. We do not attempt to identify the source of variability, which may arise from pulsation or stellar activity. However, we cross match the objects against variable star related guest observer proposals, specifying the variable type in many cases. At present eclipsing binary stars are not included. Future releases will address each K2 field as it is made available, and may be improved to include more detailed catalogue information. Detrended object lightcurves are now available for download from the catalogue website.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

In this paper we present a catalogue of Giant Molecular Clouds (GMCs) in the Andromeda (M31) galaxy extracted from the Hershel Exploitation of Local Galaxy Andromeda (HELGA) dataset. GMCs are identified from the Herschel maps using a hierarchical source extraction algorithm. We present the results of this new catalogue and characterise the spatial distribution and spectral energy properties of its clouds based on the radial dust/gas properties found by Smith et al (2012). 236 GMCs in the mass range 10^4-10^7 M_sol are identified, their cumulative mass distribution is found to be proportional to M^-1.45 in agreement with earlier studies. The GMCs appear to follow the same cloud mass to L_CO correlation observed in the Milky Way. However, comparison between this catalogue and interferometry studies also shows that the GMCs are substructured below the Herschel resolution limit suggesting that we are observing associations of GMCs. Following Gordon et al. (2006), we study the spatial structure of M31 by splitting the observed structure into a set of spiral arms and offset rings. We fit radii of 10.5 and 15.5 kpc to the two most prominent rings. We then fit a logarithmic spiral with a pitch angle of 8.9 deg to the GMCs not associated with either ring. Lastly, we comment upon the effects of deprojection on our results and investigate the effect different models for M31's inclination will have upon the projection of an unperturbed spiral arm system.

The density and temperature structures of dense cores in the L1495 cloud of the Taurus star-forming region are investigated using Herschel SPIRE and PACS images in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m and 500 $\mu$m continuum bands. A sample consisting of 20 cores, selected using spectral and spatial criteria, is analysed using a new maximum likelihood technique, COREFIT, which takes full account of the instrumental point spread functions. We obtain central dust temperatures, $T_0$, in the range 6-12 K and find that, in the majority of cases, the radial density falloff at large radial distances is consistent with the $r^{-2}$ variation expected for Bonnor-Ebert spheres. Two of our cores exhibit a significantly steeper falloff, however, and since both appear to be gravitationally unstable, such behaviour may have implications for collapse models. We find a strong negative correlation between $T_0$ and peak column density, as expected if the dust is heated predominantly by the interstellar radiation field. At the temperatures we estimate for the core centres, carbon-bearing molecules freeze out as ice mantles on dust grains, and this behaviour is supported here by the lack of correspondence between our estimated core locations and the previously-published positions of H$^{13}$CO$^+$ peaks. On this basis, our observations suggest a sublimation-zone radius typically $\sim 10^4$ AU. Comparison with previously-published N$_2$H$^+$ data at 8400 AU resolution, however, shows no evidence for N$_2$H$^+$ depletion at that resolution.

Utilizing multi-wavelength dust emission maps acquired with $Herschel$, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using inverse-Abel transform based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight which was previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 $\pm$ 0.5 K and 9.8 $\pm$0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M_sol and 11 M_sol, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, $\kappa_{\lambda} =0.1(\lambda/300 \mu m)^{-2}$, agrees to within 50% with the dust extinction constraints

Strong waves can mediate a shock transition between a pulsar wind and its surroundings, playing the role of an extended precursor, in which the energy is effectively transferred from fields to non-thermal particles. The damping of such precursors results in an essentially unmagnetized shock near the equator. In this context, we discuss the stability of strong waves and its implications for the properties of shocks. Those with stable precursors can exist in the winds of most of isolated pulsars, but the precursors may be unstable if the external pressure in the nebula is high, as in Vela-like pulsars. Pulsar wind shocks in eccentric binary systems, such as B1259-63, can acquire precursors only at certain orbital phases, and this process should be accompanied by enhanced synchro-Compton and inverse Compton emission from the precursor. The same scenario may be at work in the binary HESS J0632+057.

The results from a large field Far-Infrared (FIR) and sub-millimeter (sub-mm) survey of our neighbor galaxy M31 are presented. We have obtained Herschel images of a ~5.5x2.5 degree area centered on Andromeda. Using 21 cm atomic hydrogen maps, we are able to disentangle genuine emission from M31 from that for foreground Galactic cirrus, allowing us to recognize dusty structures out to ~31 kpc from the center. We first characterize the FIR and sub-mm morphology and then, by de-projecting Herschel maps and running an ad--hoc source extraction algorithm, we reconstruct the intrinsic morphology and the spatial distribution of the molecular complexes. Finally, we study the spatially resolved properties of the dust (temperature, emissivity, mass, etc.), by means of a pixel-by-pixel SED fitting approach.

We analyse the damping by radiation reaction and by Compton drag of strong, superluminal electromagnetic waves in the context of pulsar winds. The associated radiation signature is found by estimating the efficiency and the characteristic radiation frequencies. Applying these estimates to the gamma-ray binary containing PSR B1259-63, we show that the GeV flare observed by Fermi-LAT can be understood as inverse Compton emission by particles scattering photons from the companion star, if the pulsar wind termination shock acquires a precursor of superluminal waves roughly 30 days after periastron. This constrains the mass-loading factor of the wind $\mu=L/\dot{N}mc^2$ (where $L$ is the luminosity and $\dot{N}$ the rate of loss of electrons and positrons) to be roughly $6\times 10^4$.

We present the first overview of the Herschel observations of the nearby high-mass star-forming region NGC 7538, taken as part of the Herschel imaging study of OB Young Stellar objects (HOBYS) Key Programme. These PACS and SPIRE maps cover an approximate area of one square degree at five submillimeter and far-infrared wavebands. We have identified 780 dense sources and classified 224 of those. With the intention of investigating the existence of cold massive starless or class 0-like clumps that would have the potential to form intermediate- to high-mass stars, we further isolate 13 clumps as the most likely candidates for followup studies. These 13 clumps have masses in excess of 40 M_sun and temperatures below 15 K. They range in size from 0.4 pc to 2.5 pc and have densities between 3x10^3 cm^-3 to 4x10^4 cm^-3. Spectral energy distributions are then used to characterize their energetics and evolutionary state through a luminosity-mass diagram. NGC 7538 has a highly filamentary structure, previously unseen in the dust continuum of existing submillimeter surveys. We report the most complete imaging to date of a large, evacuated ring of material in NGC 7538 which is bordered by many cool sources.

Nonlinear electromagnetic waves with superluminal phase velocity can propagate in the winds around isolated pulsars, and around some pulsars in binary systems. Using a short-wavelength approximation, we find and analyze an integrable system of equations that govern their evolution in spherical geometry. A confined mode is identified that stagnates to finite pressure at large radius and can form a precursor to the termination shock. Using a simplified criterion, we find this mode is stable for most isolated pulsars, but may be unstable if the external pressure is high, such as in the pulsar wind nebulae in starburst galaxies and in W44. Pulsar winds in eccentric binary systems, such as PSR 1259-63, may go through phases with stable and unstable electromagnetic precursors, as well as phases in which the density is too high for these modes to propagate.

(abridged) We correlated near-infrared stellar H-Ks colour excesses of background stars from NTT/SOFI with the far-IR optical depth map, tauFIR, derived from Herschel 160, 250, 350, and 500 um data. The Herschel maps were also used to construct a model for the cloud to examine the effect of temperature gradients on the estimated optical depths and dust absorption cross-sections. A linear correlation is seen between the colour H-Ks and tauFIR up to high extinctions (AV ~ 25). The correlation translates to the average extinction ratio A250um/AJ = 0.0014 +/- 0.0002, assuming a standard near-infrared extinction law and a dust emissivity index beta=2. Using an empirical NH/AJ ratio we obtain an average absorption cross-section per H nucleus of sigmaH(250um) = (1.8 +/- 0.3) * 10^(-25) cm^2 / H-atom, corresponding to a cross-section per unit mass of gas kappaG(250 um) = 0.08 +/- 0.01 cm^2 / g. The cloud model however suggests that owing to the bias caused by temperature changes along the line-of-sight these values underestimate the true cross-sections by up to 40% near the centre of the core. Assuming that the model describes the effect of the temperature variation on tauFIR correctly, we find that the relationship between H-Ks and tauFIR agrees with the recently determined relationship between sigmaH and NH in Orion A. The derived far-IR cross-section agrees with previous determinations in molecular clouds with moderate column densities, and is not particularly large compared with some other cold cores. We suggest that this is connected to the core not beng very dense (the central density is likely to be ~10^5 cm^-3) and judging from previous molecular line data, it appears to be at an early stage of chemical evolution.

The dynamics of a standing shock front in a Poynting-flux dominated relativistic flow is investigated by using a one-dimensional, relativistic, two-fluid simulation. An upstream flow containing a circularly polarized, sinusoidal magnetic shear wave is considered, mimicking a wave driven by an obliquely rotating pulsar. It is demonstrated that this wave is converted into large amplitude electromagnetic waves with superluminal phase speeds by interacting with the shock when the shock-frame frequency of the wave exceeds the proper plasma frequency. The superluminal waves propagate in the upstream, modify the shock structure substantially, and form a well-developed precursor region ahead of a subshock. Dissipation of Poynting flux occurs in the precursor as well as in the downstream region through a parametric instability driven by the superluminal waves. The Poynting flux remaining in the downstream region is carried entirely by the superluminal waves. The downstream plasma is therefore an essentially unmagnetized, relativistically hot plasma with a non-relativistic flow speed, as suggested by observations of pulsar wind nebulae.

The whole of the Taurus region (a total area of 52 sq. deg.) has been observed by the Herschel SPIRE and PACS instruments at wavelengths of 70, 160, 250, 350 and 500 {\mu}m as part of the Herschel Gould Belt Survey. In this paper we present the first results from the part of the Taurus region that includes the Barnard 18 and L1536 clouds. A new source-finding routine, the Cardiff Source-finding AlgoRithm (CSAR), is introduced, which is loosely based on CLUMPFIND, but that also generates a structure tree, or dendrogram, which can be used to interpret hierarchical clump structure in a complex region. Sources were extracted from the data using the hierarchical version of CSAR and plotted on a mass-size diagram. We found a hierarchy of objects with sizes in the range 0.024-2.7 pc. Previous studies showed that gravitationally bound prestellar cores and unbound starless clumps appeared in different places on the mass-size diagram. However, it was unclear whether this was due to a lack of instrumental dynamic range or whether they were actually two distinct populations. The excellent sensitivity of Herschel shows that our sources fill the gap in the mass-size plane between starless and pre-stellar cores, and gives the first clear supporting observational evidence for the theory that unbound clumps and (gravitationally bound) prestellar cores are all part of the same population, and hence presumably part of the same evolutionary sequence (c.f. Simpson et al. 2011).

Unpulsed, high-energy emission from pulsar binaries can be attributed to the interaction of a pulsar wind with that of a companion star. At the shock between the outflows, particles carried away from the pulsar magnetosphere are accelerated and radiate both in synchrotron and inverse Compton processes. This emission constitutes a significant fraction of the pulsar spin-down luminosity. It is not clear however, how the highly magnetized pulsar wind could convert its mainly electromagnetic energy into the particles with such high efficiency. Here we investigate a scenario in which a pulsar striped wind converts into a strong electromagnetic wave before reaching the shock. This mode can be thought of as a shock precursor that is able to accelerate particles to ultrarelativistic energies at the expense of the electromagnetic energy it carries. Radiation of the particles leads to damping of the wave. The efficiency of this process depends on the physical conditions imposed by the external medium. Two regimes can be distinguished: a high density one, where the EM wave cannot be launched at all and the current sheets in the striped wind are first compressed by an MHD shock and subsequently dissipate by reconnection, and a low density one, where the wind can first convert into an electromagnetic wave in the shock precursor, which then damps and merges into the surroundings. Shocks in binary systems can transit from one regime to another according to binary phase. We discuss possible observational implications for these objects.

The termination shock of a pulsar wind is located roughly where the ram pressure matches that of the surrounding medium. Downstream of the shock, MHD models of the diffuse nebular emission suggest the plasma is weakly magnetized. However, the transition from a Poynting-dominated MHD wind to a particle-dominated flow is not well understood. We discuss a solution of this "sigma problem" in which a striped wind converts into a strong, superluminal electromagnetic wave. This mode slows down as it propagates radially, and its ram pressure tends to a constant value at large radius, a property we use to match the solution to the surrounding nebula. The wave thus forms a pre-cursor to the termination shock, which occurs at the point where the wave dissipates. Possible damping and dissipation mechanisms are discussed qualitatively.

We present a detailed study of how the Star Formation Rate (SFR) relates to the interstellar medium (ISM) of M31 at ~140pc scales. The SFR is calculated using the far-ultraviolet and 24um emission, corrected for the old stellar population in M31. We find a global value for the SFR of 0.25+/-0.05Msun/yr and compare this with the SFR found using the total far-infrared (FIR) luminosity. There is general agreement in regions where young stars dominate the dust heating. Atomic hydrogen (HI) and molecular gas (traced by carbon monoxide, CO) or the dust mass is used to trace the total gas in the ISM. We show that the global surface densities of SFR and gas mass place M31 amongst a set of low-SFR galaxies in the plot of Kennicutt (1998b). The relationship between SFR and gas surface density is tested in six radial annuli across M31, assuming a power law relationship with index, N. The star formation law using total gas traced by HI and CO gives a global index of N=2.03+/-0.04, with a significant variation with radius; the highest values are observed in the 10kpc ring. We suggest that this slope is due to HI turning molecular at ~10Msun/pc2. When looking at H2 regions, we measure a higher mean SFR suggesting a better spatial correlation between H2 and SF. We find N~0.6 with consistent results throughout the disk - this is at the low end of values found in previous work and argues against a superlinear SF law on small scales.

Context. G29.96-0.02 is a high-mass star-forming cloud observed at 70, 160, 250, 350, and 500 microns as part of the Herschel survey of the Galactic Plane during the Science Demonstration Phase. Aims. We wish to conduct a far-infrared study of the sources associated with this star-forming region by estimating their physical properties and evolutionary stage, and investigating the clump mass function, the star formation efficiency and rate in the cloud. Methods. We have identified the Hi-GAL sources associated with the cloud, searched for possible counterparts at centimeter and infrared wavelengths, fitted their spectral energy distribution and estimated their physical parameters. Results. A total of 198 sources have been detected in all 5 Hi-GAL bands, 117 of which are associated with 24 microns emission and 87 of which are not associated with 24 microns emission. We called the former sources 24 microns-bright and the latter ones 24 microns-dark. The [70-160] color of the 24 microns-dark sources is smaller than that of the 24 microns-bright ones. The 24 microns-dark sources have lower L_bol and L_bol/M_env than the 24 microns-bright ones for similar M_env, which suggests that they are in an earlier evolutionary phase. The G29-SFR cloud is associated with 10 NVSS sources and with extended centimeter continuum emission well correlated with the 70 microns emission. Most of the NVSS sources appear to be early B or late O-type stars. The most massive and luminous Hi-GAL sources in the cloud are located close to the G29-UC region, which suggests that there is a privileged area for massive star formation towards the center of the G29-SFR cloud. Almost all the Hi-GAL sources have masses well above the Jeans mass but only 5% have masses above the virial mass, which indicates that most of the sources are stable against gravitational collapse. The sources with M_env > M_virial and that ...

Motivated by the long-standing "luminosity problem" in low-mass star formation whereby protostars are underluminous compared to theoretical expectations, we identify 230 protostars in 18 molecular clouds observed by two Spitzer Space Telescope Legacy surveys of nearby star-forming regions. We compile complete spectral energy distributions, calculate Lbol for each source, and study the protostellar luminosity distribution. This distribution extends over three orders of magnitude, from 0.01 Lsun - 69 Lsun, and has a mean and median of 4.3 Lsun and 1.3 Lsun, respectively. The distributions are very similar for Class 0 and Class I sources except for an excess of low luminosity (Lbol < 0.5 Lsun) Class I sources compared to Class 0. 100 out of the 230 protostars (43%) lack any available data in the far-infrared and submillimeter (70 um < wavelength < 850 um) and have Lbol underestimated by factors of 2.5 on average, and up to factors of 8-10 in extreme cases. Correcting these underestimates for each source individually once additional data becomes available will likely increase both the mean and median of the sample by 35% - 40%. We discuss and compare our results to several recent theoretical studies of protostellar luminosities and show that our new results do not invalidate the conclusions of any of these studies. As these studies demonstrate that there is more than one plausible accretion scenario that can match observations, future attention is clearly needed. The better statistics provided by our increased dataset should aid such future work.

We present first results from the Herschel Gould Belt survey for the B211/L1495 region in the Taurus molecular cloud. Thanks to their high sensitivity and dynamic range, the Herschel images reveal the structure of the dense, star-forming filament B211 with unprecedented detail, along with the presence of striations perpendicular to the filament and generally oriented along the magnetic field direction as traced by optical polarization vectors. Based on the column density and dust temperature maps derived from the Herschel data, we find that the radial density profile of the B211 filament approaches a power-law behavior {\rho} {\propto} r^(-2.0{\pm}0.4) at large radii and that the temperature profile exhibits a marked drop at small radii. The observed density and temperature profiles of the B211 filament are in good agreement with a theoretical model of a cylindrical filament undergoing gravitational contraction with a polytropic equation of state: P {\propto} {\rho}^{\gamma} and T {\propto} {\rho}^({\gamma}-1), with {\gamma}=0.97{\pm}0.01<1 (i.e. not strictly isothermal). The morphology of the column density map, where some of the perpendicular striations are apparently connected to the B211 filament, further suggests that the material may be accreting along the striations onto the main filament. The typical velocities expected for the infalling material in this picture are ~0.5-1 km/s, which are consistent with the existing kinematical constraints from previous CO observations.

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with spatial resolution of ~140pc and find that a variable dust-emissivity index (beta) is required to fit the data. We find no significant long-wavelength excess above this model suggesting there is no cold dust component. We show that the gas-to-dust ratio varies radially, increasing from ~20 in the center to ~70 in the star-forming ring at 10kpc, consistent with the metallicity gradient. In the 10kpc ring the average beta is ~1.9, in good agreement with values determined for the Milky Way (MW). However, in contrast to the MW, we find significant radial variations in beta, which increases from 1.9 at 10kpc to ~2.5 at a radius of 3.1kpc and then decreases to 1.7 in the center. The dust temperature is fairly constant in the 10kpc ring (ranging from 17-20K), but increases strongly in the bulge to ~30K. Within 3.1kpc we find the dust temperature is highly correlated with the 3.6 micron flux, suggesting the general stellar population in the bulge is the dominant source of dust heating there. At larger radii, there is a weak correlation between the star formation rate and dust temperature. We find no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value of (1.9+/-0.4)x10^20 cm^-2 [K kms^-1]^-1.

With recent Herschel observations, the northern filament of the Corona Australis cloud has now been mapped in a number of bands from 1.2um to 870um. The data set provides a good starting point for the study of the cloud over several orders of magnitude in density. We wish to examine the differences of the column density distributions derived from dust extinction, scattering, and emission, and to determine to what extent the observations are consistent with the standard dust models. From Herschel data, we calculate the column density distribution that is compared to the corresponding data derived in the near-infrared regime from the reddening of the background stars, and from the surface brightness attributed to light scattering. We construct three-dimensional radiative transfer models to describe the emission and the scattering. The scattered light traces low column densities of A_V~1mag better than the dust emission, remaining useful to A_V ~ 10-15 mag. Based on the models, the extinction and the level of dust emission are surprisingly consistent with a sub-millimetre dust emissivity typical of diffuse medium. However, the intensity of the scattered light is very low at the centre of the densest clump and this cannot be explained without a very low grain albedo. Both the scattered light and dust emission indicate an anisotropic radiation field. The modelling of the dust emission suggests that the radiation field intensity is at least three times the value of the normal interstellar radiation field. The inter-comparison between the extinction, light scattering, and dust emission provides very stringent constraints on the cloud structure, the illuminating radiation field, and the grain properties.

In a previous paper we identified cores within infrared dark clouds (IRDCs). We regarded those without embedded sources as the least evolved, and labelled them starless. Here we identify the most isolated starless cores and model them using a three-dimensional, multi-wavelength, Monte Carlo, radiative transfer code. We derive the cores' physical parameters and discuss the relation between the mass, temperature, density, size and the surrounding interstellar radiation field (ISRF) for the cores. The masses of the cores were found not to correlate with their radial size or central density. The temperature at the surface of a core was seen to depend almost entirely on the level of the ISRF surrounding the core. No correlation was found between the temperature at the centre of a core and its local ISRF. This was seen to depend, instead, on the density and mass of the core.

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

A growing body of evidence indicates that the formation of filaments in interstellar clouds is a key component of the star formation process. In this paper, we present new Herschel PACS and SPIRE observations of the B59 and Stem regions in the Pipe Nebula complex, revealing a rich, organized network of filaments. The asymmetric column density profiles observed for several filaments, along with the bow-like edge of B59, indicates that the Pipe Nebula is being compressed from its western side, most likely by the winds from the nearby Sco OB2 association. We suggest that this compressive flow has contributed to the formation of some of the observed filamentary structures. In B59, the only region of the entire Pipe complex showing star formation activity, the same compressive flow has likely enhanced the initial column density of the clump, allowing it to become globally gravitationally unstable. Although more speculative, we propose that gravity has also been responsible for shaping the converging filamentary pattern observed in B59. While the question of the relative impact of large-scale compression and gravity remains open in B59, large-scale compression appears to be a plausible mechanism for the initial formation of filamentary structures in the rest of the complex

Pulsar winds are the ideal environment for the study of non-linear electromagnetic waves. It is generally thought that a pulsar launches a striped wind, a magnetohydrodynamic entropy wave, where plasma sheets carried along with the flow separate regions of alternating magnetic field. But when the density drops below a critical value, or equivalently for distances from the pulsar greater than a critical radius, a strong superluminal wave can also propagate. In this contribution we discuss the conversion of the equatorial striped wind into a linearly polarized superluminal wave, and we argue that this mode is important for the conversion of Poynting flux to kinetic energy flux before the outflow reaches the termination shock.

As part of the Herschel guaranteed time key program 'HOBYS', we present the photometric survey of the star forming region Vela-C, one of the nearest sites of low-to-high-mass star formation in the Galactic plane. Vela-C has been observed with PACS and SPIRE in parallel mode between 70 um and 500 um over an area of about 3 square degrees. A photometric catalogue has been extracted from the detections in each band, using a threshold of 5 sigma over the local background. Out of this catalogue we have selected a robust sub-sample of 268 sources, of which 75% are cloud clumps and 25% are cores. Their Spectral Energy Distributions (SEDs) have been fitted with a modified black body function. We classify 48 sources as protostellar and 218 as starless. For two further sources, we do not provide a secure classification, but suggest they are Class 0 protostars. From SED fitting we have derived key physical parameters. Protostellar sources are in general warmer and more compact than starless sources. Both these evidences can be ascribed to the presence of an internal source(s) of moderate heating, which also causes a temperature gradient and hence a more peaked intensity distribution. Moreover, the reduced dimensions of protostellar sources may indicate that they will not fragment further. A virial analysis of the starless sources gives an upper limit of 90% for the sources gravitationally bound and therefore prestellar. We fit a power law N(logM) prop M^-1.1 to the linear portion of the mass distribution of prestellar sources. This is in between that typical of CO clumps and those of cores in nearby star-forming regions. We interpret this as a result of the inhomogeneity of our sample, which is composed of comparable fractions of clumps and cores.

We have used data taken as part of the Herschel infrared Galactic Plane survey (Hi-GAL) to study 3171 infrared-dark cloud (IRDC) candidates that were identified in the mid-infrared (8 {\mu}m) by Spitzer (we refer to these as 'Spitzer-dark' regions). They all lie in the range l=300 - 330 \circ and |b| 6 1 \circ. Of these, only 1205 were seen in emission in the far-infrared (250-500 {\mu}m) by Herschel (we call these 'Herschel-bright' clouds). It is predicted that a dense cloud will not only be seen in absorption in the mid-infrared, but will also be seen in emission in the far-infrared at the longest Herschel wavebands (250-500 {\mu}m). If a region is dark at all wavelengths throughout the mid-infrared and far-infrared, then it is most likely to be simply a region of lower background infrared emission (a 'hole in the sky'). Hence, it appears that previous surveys, based on Spitzer and other mid-infrared data alone, may have over-estimated the total IRDC population by a factor of 2. This has implications for estimates of the star formation rate in IRDCs in the Galaxy.We studied the 1205 Herschel-bright IRDCs at 250 {\mu}m, and found that 972 of them had at least one clearly defined 250-{\mu}m peak, indicating that they contained one or more dense cores. Of these, 653 (67 per cent) contained an 8-{\mu}m point source somewhere within the cloud, 149 (15 per cent) contained a 24-{\mu}m point source but no 8-{\mu}m source, and 170 (18 per cent) contained no 24-{\mu}m or 8-{\mu}m point sources. We use these statistics to make inferences about the lifetimes of the various evolutionary stages of IRDCs.

Recent work on the properties of superluminal waves in pulsar winds is summarized. It is speculated that these waves play an important role in the termination shock that divides the wind from the surrounding nebula.

The filament IRDC G035.39--00.33 in the W48 molecular complex is one of the darkest infrared clouds observed by \textit{Spitzer}. It has been observed by the PACS (70 and 160\,$\micron$) and SPIRE (250, 350, and 500\,$\micron$) cameras of the \textit{Herschel} Space Observatory as part of the W48 molecular cloud complex in the framework of the HOBYS key programme. The observations reveal a sample of 28 compact sources (deconvolved FWHM sizes $<$0.3 pc) complete down to $\sim$$5 \msun$ in G035.39--00.33 and its surroundings. Among them, 13 compact sources are massive dense cores with masses $>$$20 \msun$. The cloud characteristics we derive from the analysis of their spectral energy distributions are masses of $20-50 \msun$, sizes of 0.1--0.2 pc, and average densities of $2-20 \times 10^{5} \cmc$, which make these massive dense cores excellent candidates to form intermediate- to high-mass stars. Most of the massive dense cores are located inside the G035.39--00.33 ridge and host IR-quiet high-mass protostars. The large number of protostars found in this filament suggests that we are witnessing a mini-burst of star formation with an efficiency of $\sim$15% and a rate density of $\sim$$40 \msun\,$yr$^{-1}\,$kpc$^{-2}$ within $\sim$8 pc$^2$, a large area covering the full ridge. Part of the extended SiO emission observed towards G035.39--00.33 is not associated with obvious protostars and may originate from low-velocity shocks within converging flows, as advocated by previous studies.

The energy lost by a rotation-powered pulsar is carried by a relativistic flow containing a mixture of electromagnetic fields and particles. In the inner regions, this is thought to be a magnetically dominated, cold, electron-positron wind that is well described by the MHD equations. However, beyond a critical radius r_{cr}, the same particle, energy and momentum fluxes can be transported by a strong, transverse electromagnetic wave with superluminal phase speed. We analyze the nonlinear dispersion relation of these waves for linear and circular polarization, and find the dependence of r_{cr} on the mass-loading, magnetization and luminosity of the flow, as well as on the net magnetic flux. We show that, for most isolated pulsars, the wind lies well outside r_{cr}, and speculate that superluminal modes play an important role in the dissipation of electromagnetic energy into nonthermal particles at the termination shock.

We present the first Herschel PACS and SPIRE results of the Vela C molecular complex in the far-infrared and submillimetre regimes at 70, 160, 250, 350, and 500 um, spanning the peak of emission of cold prestellar or protostellar cores. Column density and multi-resolution analysis (MRA) differentiates the Vela C complex into five distinct sub-regions. Each sub-region displays differences in their column density and temperature probability distribution functions (PDFs), in particular, the PDFs of the `Centre-Ridge' and `South-Nest' sub-regions appear in stark contrast to each other. The Centre-Ridge displays a bimodal temperature PDF representative of hot gas surrounding the HII region RCW 36 and the cold neighbouring filaments, whilst the South-Nest is dominated by cold filamentary structure. The column density PDF of the Centre-Ridge is flatter than the South-Nest, with a high column density tail, consistent with formation through large-scale flows, and regulation by self-gravity. At small to intermediate scales MRA indicates the Centre-Ridge to be twice as concentrated as the South-Nest, whilst on larger scales, a greater portion of the gas in the South-Nest is dominated by turbulence than in the Centre-Ridge. In Vela C, high-mass stars appear to be preferentially forming in ridges, i.e., dominant high column density filaments.

High energy emission from blazars is thought to arise in a relativistic jet launched by a supermassive black hole. The emission site must be far from the hole and the jet relativistic, in order to avoid absorption of the photons. In extreme cases, rapid variability of the emission suggests that structures of length-scale smaller than the gravitational radius of the central black hole are imprinted on the jet as it is launched, and modulate the radiation released after it has been accelerated to high Lorentz factor. We propose a mechanism which can account for the acceleration of the jet, and for the rapid variability of the radiation, based on the propagation characteristics of large-amplitude waves in charge-starved, polar jets. Using a two-fluid (electron-positron) description, we find the outflows exhibit a delayed acceleration phase, that starts at roughly 1pc, where the inertia associated with the wave currents becomes important. The time-structure imprinted on the jet at launch modulates photons produced by the accelerating jet provided the pair multiplicity in the black-hole magnetosphere is sufficiently small, suggesting that very rapid variability is confined to sources in which the electromagnetic cascade in the black-hole magnetosphere is not prolific.

We present high angular resolution observations, taken with the Very Large Array (VLA) and Multiple Element Radio Linked Interferometer Network (MERLIN) radio telescopes, at 7mm and 4.4cm respectively, of the prototype Class 0 protostar VLA1623. At 7mm we detect two sources (VLA1623A & B) coincident with the two previously detected components at the centre of this system. The separation between the two is 1.2arcsec, or ~170AU at an assumed distance of 139pc. The upper limit to the size of the source coincident with each component of VLA1623 is ~0.7arcsec, in agreement with previous findings. This corresponds to a diameter of ~100AU at an assumed distance of 139pc. Both components show the same general trend in their broadband continuum spectra, of a steeper dust continuum spectrum shortward of 7mm and a flatter spectrum longward of this. We estimate an upper limit to the VLA1623A disc mass of <0.13Msol and an upper limit to its radius of ~50AU. The longer wavelength data have a spectral index of \alpha~0.6+/-0.3. This is too steep to be explained by optically thin free-free emission. It is most likely due to optically thick free-free emission. Alternatively, we speculate that it might be due to the formation of larger grains or planetesimals in the circumstellar disc. We estimate the mass of VLA1623B to be <0.15M$sol. We can place a lower limit to its size of ~30x7 AU, and an upper limit to its diameter of ~100AU. The longer wavelength data of VLA1623B also have a spectral index of \alpha~0.6+/-0.3. The nature of VLA1623B remains a matter of debate. It could be a binary companion to the protostar, or a knot in the radio jet from VLA1623A.

We provide a first look at the results of the Herschel Gould Belt survey toward the IC5146 molecular cloud and present a preliminary analysis of the filamentary structure in this region. The column density map, derived from our 70-500 micron Herschel data, reveals a complex network of filaments, and confirms that these filaments are the main birth sites of prestellar cores. We analyze the column density profiles of 27 filaments and show that the underlying radial density profiles fall off as r^{-1.5} to r^{-2.5} at large radii. Our main result is that the filaments seem to be characterized by a narrow distribution of widths having a median value of 0.10 +- 0.03 pc, which is in stark contrast to a much broader distribution of central Jeans lengths. This characteristic width of ~0.1 pc corresponds to within a factor of ~2 to the sonic scale below which interstellar turbulence becomes subsonic in diffuse gas, supporting the argument that the filaments may form as a result of the dissipation of large-scale turbulence.

The densest infrared dark clouds (IRDCs) may represent the earliest observable stage of high-mass star formation. These clouds are very cold, hence they emit mainly at far-infrared and sub-mm wavelengths. For the first time, Herschel has provided multi-wavelength, spatially resolved observations of cores within IRDCs, which, when combined with radiative transfer modelling, can constrain their properties, such as mass, density profile and dust temperature. We use a 3D, multi-wavelength Monte Carlo radiative transfer code to model in detail the emission from six cores in three typical IRDCs seen in the Hi-GAL survey (G030.50+00.95, G031.03+00.26 and G031.03+00.76), and thereby to determine the properties of these cores and compare them with their low-mass equivalents. We found masses ranging from 90 to 290 solar masses with temperatures from 8 to 11K at the centre of each core and 18 to 28K at the surface. The maximum luminosity of an embedded star within each core was calculated, and we rule out the possibility of significant high mass star formation having yet occurred in three of our cores.

Particle-in-cell (PIC) simulations of relativistic shocks are in principle capable of predicting the spectra of photons that are radiated incoherently by the accelerated particles. The most direct method evaluates the spectrum using the fields given by the Lienard-Wiechart potentials. However, for relativistic particles this procedure is computationally expensive. Here we present an alternative method, that uses the concept of the photon formation length. The algorithm is suitable for evaluating spectra both from particles moving in a specific realization of a turbulent electromagnetic field, or from trajectories given as a finite, discrete time series by a PIC simulation. The main advantage of the method is that it identifies the intrinsic spectral features, and filters out those that are artifacts of the limited time resolution and finite duration of input trajectories.

The radio, optical, X-ray and gamma-ray nebulae that surround many pulsars are thought to arise from synchrotron and inverse Compton emission. The energy powering this emission, as well as the magnetic fields and relativistic particles, are supplied by a "wind" driven by the central object. The inner parts of the wind can be described using the equations of MHD, but these break down in the outer parts, when the density of charge carriers drops below a critical value. This paper reviews the wave properties of the inner part (striped wind), and uses a relativistic two-fluid model (cold electrons and positrons) to re-examine the nonlinear electromagnetic modes that propagate in the outer parts. It is shown that in a radial wind, two solutions exist for circularly polarised electromagnetic modes. At large distances one of them turns into a freely expanding flow containing a vacuum wave, whereas the other decelerates, corresponding to a confined flow.

We demonstrate the use of the 3D Monte Carlo radiative transfer code PHAETHON to model infrared-dark clouds (IRDCs) that are externally illuminated by the interstellar radiation field (ISRF). These clouds are believed to be the earliest observed phase of high-mass star formation, and may be the high-mass equivalent of lower-mass prestellar cores. We model three different cases as examples of the use of the code, in which we vary the mass, density, radius, morphology and internal velocity field of the IRDC. We show the predicted output of the models at different wavelengths chosen to match the observing wavebands of Herschel and Spitzer. For the wavebands of the long- wavelength SPIRE photometer on Herschel, we also pass the model output through the SPIRE simulator to generate output images that are as close as possible to the ones that would be seen using SPIRE. We then analyse the images as if they were real observations, and compare the results of this analysis with the results of the radiative transfer models. We find that detailed radiative transfer modelling is necessary to accurately determine the physical parameters of IRDCs (e.g. dust temperature, density profile). This method is applied to study G29.55+00.18, an IRDC observed by the Herschel Infrared Galactic Plane survey (Hi-GAL), and in the future it will be used to model a larger sample of IRDCs from the same survey.

We summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class 0 protostars can be identified in the Aquila field, while ~ unbound starless cores and no protostars are observed in the Polaris field. The prestellar core mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability.

The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result.

We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

We present a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30" to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.

Our PACS and SPIRE images of the Aquila Rift and part of the Polaris Flare regions, taken during the science demonstration phase of Herschel discovered fascinating, omnipresent filamentary structures that appear to be physically related to compact cores. We briefly describe a new multi-scale, multi-wavelength source extraction method used to detect objects and measure their parameters in our Herschel images. All of the extracted starless cores (541 in Aquila and 302 in Polaris) appear to form in the long and very narrow filaments. With its combination of the far-IR resolution and sensitivity, Herschel directly reveals the filaments in which the dense cores are embedded; the filaments are resolved and have deconvolved widths of 35 arcsec in Aquila and 59 arcsec in Polaris (9000 AU in both regions). Our first results of observations with Herschel enable us to suggest that in general dense cores may originate in a process of fragmentation of complex networks of long, thin filaments, likely formed as a result of an interplay between gravity, interstellar turbulence, and magnetic fields. To unravel the roles of the processes, one has to obtain additional kinematic and polarization information; these follow-up observations are planned.

We present far-infrared spectra and maps of the DR21 molecular cloud core between 196 and 671 microns, using the Herschel-SPIRE spectrometer. Nineteen molecular lines originating from CO, 13CO, HCO+ and H2O, plus lines of [N II] and [CI] were recorded, including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nH2 ~ 7 x 10^4 cm-3, CO column density N(CO) ~ 3.5 x 10^18 cm^-2 and a filling factor of ~ 12%, and appear to trace gas associated with an outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5 x 10^21 cm^-2. Astronomy & Astrophysics HERSCHEL special Issue, in press.

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index {\beta} and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 {\mu}m we determine T_dust and {\beta} throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of {\beta} ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between {\beta} and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

The origin and possible universality of the stellar initial mass function (IMF) is a major issue in astrophysics. One of the main objectives of the Herschel Gould Belt Survey is to clarify the link between the prestellar core mass function (CMF) and the IMF. We present and discuss the core mass function derived from Herschel data for the large population of prestellar cores discovered with SPIRE and PACS in the Aquila Rift cloud complex at d ~ 260 pc. We detect a total of 541 starless cores in the entire ~11 deg^2 area of the field imaged at 70-500 micron with SPIRE/PACS. Most of these cores appear to be gravitationally bound, and thus prestellar in nature. Our Herschel results confirm that the shape of the prestellar CMF resembles the stellar IMF, with much higher quality statistics than earlier submillimeter continuum ground-based surveys.

The Polaris Flare cloud region contains a great deal of extended emission. It is at high declination and high Galactic latitude. It was previously seen strongly in IRAS Cirrus emission at 100 microns. We have detected it with both PACS and SPIRE on Herschel. We see filamentary and low-level structure. We identify the five densest cores within this structure. We present the results of a temperature, mass and density analysis of these cores. We compare their observed masses to their virial masses, and see that in all cases the observed masses lie close to the lower end of the range of estimated virial masses. Therefore, we cannot say whether they are gravitationally bound prestellar cores. Nevertheless, these are the best candidates to be potentialprestellar cores in the Polaris cloud region.

As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has been observed. The complete ~ 3.3deg x 3.3deg imaging with SPIRE 250/350/500 micron and PACS 70/160 micron allows a deep investigation of embedded protostellar phases, probing of the dust emission from warm inner regions at 70 and 160 micron to the bulk of the cold envelopes between 250 and 500 micron. We used a systematic detection technique operating simultaneously on all Herschel bands to build a sample of protostars. Spectral energy distributions are derived to measure luminosities and envelope masses, and to place the protostars in an M_env - L_bol evolutionary diagram. The spatial distribution of protostars indicates three star-forming sites in Aquila, with W40/Sh2-64 HII region by far the richest. Most of the detected protostars are newly discovered. For a reduced area around the Serpens South cluster, we could compare the Herschel census of protostars with Spitzer results. The Herschel protostars are younger than in Spitzer with 7 Class 0 YSOs newly revealed by Herschel. For the entire Aquila field, we find a total of ~ 45-60 Class 0 YSOs discovered by Herschel. This confirms the global statistics of several hundred Class~0 YSOs that should be found in the whole Gould Belt survey.

Aims. To follow the species chemistry arising in diverse sources of the Galaxy with Herschel. Methods. SPIRE FTS sparse sampled maps of the Orion bar & compact HII regions G29.96-0.02 and G32.80+0.19 have been analyzed. Results. Beyond the wealth of atomic and molecular lines detected in the high-resolution spectra obtained with the FTS of SPIRE in the Orion Bar, one emission line is found to lie at the position of the fundamental rotational transition of CH+ as measured precisely in the laboratory (Pearson & Drouion 2006). This coincidence suggests that it is the first detection of the fundamental rotational transition of CH+. This claim is strengthened by the observation of the lambda doublet transitions arising from its relative, CH, which are also observed in the same spectrum. The broad spectral coverage of the SPIRE FTS allows for the simultaneous measurement of these closely related chemically species, under the same observing conditions. The importance of these lines are discussed and a comparison with results obtained from models of the Photon Dominated Region (PDR) of Orion are presented. The CH+ line also appears in absorption in the spectra of the two galactic compact HII regions G29.96-0.02 and G32.80+0.19, which is likely due to the presence of CH+ in the the Cold Neutral Medium of the galactic plane. These detections will shed light on the formation processes and on the existence of CH+, which are still outstanding questions in astrophysics.

We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum of the Ultra-compact HII region G29.96-0.02. Spectral lines from species including 13CO, CO, [CI], and [NII] are detected. A sparse map of the [NII] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. The quality of the spectrum longwards of 100 {\mu}m allows us to fit a single temperature greybody with temperature 80.3\pm0.6K and dust emissivity index 1.73\pm0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 Msol for the clump containing the HII region. The clump's bolometeric luminosity of 4 x 10^6 Lsol is comparable to, or slightly greater than, the known O-star powering the UCHII region.

Context. RCW 120 is a well-studied, nearby Galactic HII region with ongoing star formation in its surroundings. Previous work has shown that it displays a bubble morphology at mid-infrared wavelengths and has a massive layer of collected neutral material seen at sub-mm wavelengths. Given the well-defined photo-dissociation region (PDR) boundary and collected layer, it is an excellent laboratory to study the "collect and collapse" process of triggered star formation. Using Herschel Space Observatory data at 100, 160, 250, 350, and 500 micron, in combination with Spitzer and APEX-LABOCA data, we can for the first time map the entire spectral energy distribution of an HII region at high angular resolution. Aims. We seek a better understanding of RCW120 and its local environment by analysing its dust temperature distribution. Additionally, we wish to understand how the dust emissivity index, beta, is related to the dust temperature. Methods. We determine dust temperatures in selected regions of the RCW 120 field by fitting their spectral energy distribution (SED), derived using aperture photometry. Additionally, we fit the SED extracted from a grid of positions to create a temperature map. Results. We find a gradient in dust temperature, ranging from >30 K in the interior of RCW 120, to ~20K for the material collected in the PDR, to ~10K toward local infrared dark clouds and cold filaments. Our results suggest that RCW 120 is in the process of destroying the PDR delineating its bubble morphology. The leaked radiation from its interior may influence the creation of the next generation of stars. We find support for an anti-correlation between the fitted temperature and beta, in rough agreement with what has been found previously. The extended wavelength coverage of the Herschel data greatly increases the reliability of this result.

Infrared dark clouds (IRDCs) are cold and dense reservoirs of gas potentially available to form stars. Many of these clouds are likely to be pristine structures representing the initial conditions for star formation. The study presented here aims to construct and analyze accurate column density and dust temperature maps of IRDCs by using the first Herschel data from the Hi-GAL galactic plane survey. These fundamental quantities, are essential for understanding processes such as fragmentation in the early stages of the formation of stars in molecular clouds. We have developed a simple pixel-by-pixel SED fitting method, which accounts for the background emission. By fitting a grey-body function at each position, we recover the spatial variations in both the dust column density and temperature within the IRDCs. This method is applied to a sample of 22 IRDCs exhibiting a range of angular sizes and peak column densities. Our analysis shows that the dust temperature decreases significantly within IRDCs, from background temperatures of 20-30 K to minimum temperatures of 8-15 K within the clouds, showing that dense molecular clouds are not isothermal. Temperature gradients have most likely an important impact on the fragmentation of IRDCs. Local temperature minima are strongly correlated with column density peaks, which in a few cases reach NH2 = 1 x 10^{23} cm^{-2}, identifying these clouds as candidate massive prestellar cores. Applying this technique to the full Hi-GAL data set will provide important constraints on the fragmentation and thermal properties of IRDCs, and help identify hundreds of massive prestellar core candidates.

It has been shown that by means of different physical mechanisms the expansion of HII regions can trigger the formation of new stars of all masses. This process may be important to the formation of massive stars but has never been quantified in the Galaxy. We use Herschel-PACS and -SPIRE images from the Herschel Infrared survey of the Galactic plane, Hi-GAL, to perform this study. We combine the Spitzer-GLIMPSE and -MIPSGAL, radio-continuum and sub-millimeter surveys such as ATLASGAL with Hi-GAL to study Young Stellar Objects (YSOs) observed towards Galactic HII regions. We select a representative HII region, N49, located in the field centered on l=30 degr observed as part of the Hi-GAL Science Demonstration Phase, to demonstrate the importance Hi-GAL will have to this field of research. Hi-GAL PACS and SPIRE images reveal a new population of embedded young stars, coincident with bright ATLASGAL condensations. The Hi-GAL images also allow us, for the first time, to constrain the physical properties of the newly formed stars by means of fits to their spectral energy distribution. Massive young stellar objects are observed at the borders of the N49 region and represent second generation massive stars whose formation has been triggered by the expansion of the ionized region. Hi-GAL enables us to detect a population of young stars at different evolutionary stages, cold condensations only being detected in the SPIRE wavelength range. The far IR coverage of Hi-GAL strongly constrains the physical properties of the YSOs. The large and unbiased spatial coverage of this survey offers us a unique opportunity to lead, for the first time, a global study of star formation triggered by HII regions in our Galaxy.

(abbrev.) Linear acceleration emission occurs when a charged particle is accelerated parallel to its velocity. We evaluate the spectral and angular distribution of this radiation for several special cases, including constant acceleration (hyperbolic motion) of finite duration. Based on these results, we find the following general properties of the emission from an electron in a linear accelerator that can be characterized by an electric field E acting over a distance L: (i) the spectrum extends to a cut-off photon energy ~ LE^2 MeV, where E is in units of the Schwinger critical field and L in units of the Compton wavelength of the electron. (ii) the total energy emitted by a particle traversing the accelerator is in agreement with the standard Larmor formula (iii) the low frequency spectrum is flat for hyperbolic trajectories, but in general depends on the details of the accelerator. We also show that linear acceleration emission complements curvature radiation in the strongly magnetized pair formation regions in pulsar magnetospheres. It dominates when the length L of the accelerator is less than the formation length of curvature photons, which is given by the ratio of the radius of curvature of the magnetic field lines to the Lorentz factor of the particle. In standard static models of pair creating regions linear acceleration emission is negligible, but it is important in more realistic dynamical models in which the accelerating field fluctuates on a short length-scale.

(Abbreviated) Particle-in-cell simulations of relativistic, weakly magnetized collisionless shocks show that particles can gain energy by repeatedly crossing the shock front. This requires scattering off self-generated small length-scale magnetic fluctuations. The radiative signature of this first-order Fermi acceleration mechanism is important for models of both the prompt and afterglow emission in gamma-ray bursts and depends on the strength parameter "a" of the fluctuations. For electrons (and positrons), acceleration saturates when the radiative losses produced by the scattering cannot be compensated by the energy gained on crossing the shock. We show that this sets an upper limit on both the electron Lorentz factor and on the energy of the photons radiated during the scattering process. This rules out "jitter" radiation on self-excited fluctuations with a < 1 as a source of gamma-rays, although high-energy photons might still be produced when the jitter photons are upscattered in an analog of the synchrotron self-Compton process. In fluctuations with a > 1, radiation is generated by the standard synchrotron mechanism, and the maximum photon energy rises linearly with a, until saturating at approximately 70 MeV.

We have used the Berkeley-Illinois-Maryland Array (BIMA) to make deep N2H+ 1-0 maps of the pre-stellar core L183, in order to study the spatial and kinematic substructure within the densest part of the core. Three spatially and kinematically distinct clumps are detected, which we label L183-N1, L183-N2 and L183-N3. L183-N2 is approximately coincident with the submillimetre dust peak and lies at the systemic velocity of L183. Thus we conclude that L183-N2 is the central dense core of L183. L183-N1 and 3 are newly-discovered fragments of L183, which are marked by velocity gradients that are parallel to, but far stronger than, the velocity gradient of L183 as a whole, as detected in previous single-dish data. Furthermore, the ratio of the large-scale and small-scale velocity gradients, and the ratio of their respective size-scales, are consistent with the conservation of angular momentum for a rotating, collapsing core undergoing spin-up. The inferred axis of rotation is parallel to the magnetic field direction, which is offset from its long axis, as we have seen in other pre-stellar cores. Therefore, we propose that we have detected a fragmenting, collapsing, filamentary, pre-stellar core, rotating about its B-field, which is spinning up as it collapses. It will presumably go on to form a multiple protostellar system.

We present optical and submillimetre polarimetry data of the Bok globule CB3 and optical polarimetry data of the Bok globule CB246. We use each set of polarimetry data to infer the B-field orientation in each of the clouds. The optical data can only be used in the low density, low extinction edge regions of clouds. The submillimetre data can only be used in the high column-density, central regions of the clouds. It has previously been found that near-infrared polarisation mapping of background stars does not accurately trace the magnetic field in dense cloud regions. This may be due to a lack of aligned grains in dense regions. We test this by comparing the field orientations measured by our two independent polarimetry methods. We find that the field orientation deduced from the optical data matches up well with the orientation estimated from the submillimetre data. We therefore claim that both methods are accurately tracing the same magnetic field in CB3. Hence, in this case, there must be significant numbers of aligned dust grains in the high density region, and they do indeed trace the magnetic field in the submillimetre. We find an offset of 40$\pm$14 degrees between the magnetic field orientation and the short axis of the globule. This is consistent with the mean value of 31$\pm$3 degrees found in our previous work on prestellar cores, even though CB3 is a protostellar core. Taken together, the six prestellar cores that we have now studied in this way show a mean offset between magnetic field orientation and core short axis of $\sim30\pm$3 degrees, in apparent contradiction with some models of magnetically dominated star formation.

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the "knee" at 10^{15.5}eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized non-propagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves, can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times we observe saturation of the instability in the simulation, but the mechanism responsible is an artefact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

We present Spitzer IRAC (~2 deg^2) and MIPS (~8 deg^2) observations of the Cepheus Flare which is associated with the Gould Belt, at an approximate distance of ~300 pc. Around 6500 sources are detected in all four IRAC bands, of which ~900 have MIPS 24 micron detections. We identify 133 YSO candidates using color-magnitude diagram techniques, a large number of the YSO candidates are associated with the NGC 7023 reflection nebula. Cross identifications were made with the Guide Star Catalog II and the IRAS Faint Source Catalog, and spectral energy distributions (SED) were constructed. SED modeling was conducted to estimate the degree of infrared excess. It was found that a large majority of disks were optically thick accreting disks, suggesting that there has been little disk evolution in these sources. Nearest-neighbor clustering analysis identified four small protostellar groups (L1228, L1228N, L1251A, and L1251B) with 5-15 members each and the larger NGC 7023 association with 32 YSO members. The star formation efficiency for cores with clusters of protostars and for those without clusters was found to be ~8% and ~1% respectively. The cores L1155, L1241, and L1247 are confirmed to be starless down to our luminosity limit of L_bol=0.06 L_sol.

Stationary solutions to the equations of non-linear diffusive shock acceleration play a fundamental role in the theory of cosmic-ray acceleration. Their existence usually requires that a fraction of the accelerated particles be allowed to escape from the system. Because the scattering mean-free-path is thought to be an increasing function of energy, this condition is conventionally implemented as an upper cut-off in energy space -- particles are then permitted to escape from any part of the system, once their energy exceeds this limit. However, because accelerated particles are responsible for substantial amplification of the ambient magnetic field in a region upstream of the shock front, we examine an alternative approach in which particles escape over a spatial boundary. We use a simple iterative scheme that constructs stationary numerical solutions to the coupled kinetic and hydrodynamic equations. For parameters appropriate for supernova remnants, we find stationary solutions with efficient acceleration when the escape boundary is placed at the point where growth and advection of strongly driven non-resonant waves are in balance. We also present the energy dependence of the distribution function close to the energy where it cuts off - a diagnostic that is in principle accessible to observation.

We examine the problem of particle acceleration at a relativistic shocks assuming pitch-angle scattering and using a Hartree-Fock method to approximate the associated eigenfunctions. This leads to a simple transcendental equation determining the power-law index, $s$, given the up and downstream velocities. We compare our results with accurate numerical solutions obtained using the eigenfunction method. In addition to the power-law index this method yields the angular and spatial distributions upstream of the shock.

Using recent polarimetric observations of the Crab Nebula in the hard X-ray band by INTEGRAL, we show that the absence of vacuum birefringence effects constrains O(E/M) Lorentz violation in QED to the level |\xi| < 9x10^{-10} at three sigma CL, tightening by more than three orders of magnitude previous constraints. We show that planned X-ray polarimeters have the potential the potential to probe |\xi|~ 10^{-16} by detecting polarization in active galaxies at red-shift ~1.

We investigate the so-called non-resonant cosmic-ray streaming instability, first discussed by Bell (2004). The extent to which thermal damping and ion-neutral collisions reduce the growth of this instability is calculated. Limits on the growth of the non-resonant mode in SN1006 and RX J1713.7-3946 are presented.

The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic-ray driven instability, in which the unstable waves are excited beyond the linear regime. In a snapshot of the resultant field, particle transport simulations are carried out. The use of a static snapshot of the field is reasonable given that the Larmor period for particles is typically very short relative to the instability growth time. The diffusion rate is found to be close to, or below, the Bohm limit for a range of energies. This provides the first explicit demonstration that self-excited turbulence reduces the diffusion coefficient and has important implications for cosmic ray transport and acceleration in supernova remnants.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

We present continuum data from the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), and the Mid-Infrared Photometer for Spitzer (MIPS) on the Spitzer Space Telescope, at submillimetre and infrared wavelengths respectively. We study the Taurus molecular cloud 1 (TMC1), and in particular the region of the Taurus Molecular Ring (TMR). In the continuum data we see no real evidence for a ring, but rather we see one side of it only, appearing as a filament. We name the filament `the bull's tail'. The filament is seen in emission at 850, 450 and 160um, and in absorption at 70um. We compare the data with archive data from the Infra-Red Astronomical Satellite (IRAS) at 12, 25, 60, 100um, in which the filament is also seen in absorption. We find that the emission from the filament consists of two components: a narrow, cold (~8K), central core; and a broader, slightly warmer (~12K), shoulder of emission. We use a radiative transfer code to model the filament's appearance, either in emission or absorption, simultaneously at each of the different wavelengths. Our best fit model uses a Plummer-like density profile and a homogeneous interstellar dust grain population. Unlike previous work on a similar, but different filament in Taurus, we require no grain coagulation to explain our data.

Rapid inverse Compton cooling sets in when the brightness temperature (T_B) of a self-absorbed synchrotron source with power-law electrons reaches ~10^{12} K. However, T_B inferred from observations of intra-day variable sources (IDV) are well above the "Compton catastrophe" limit. This can be understood if the underlying electron distribution cuts off at low energy. We approximate a low-energy cut-off with monoenergetic electrons. We compute the synchrotron self-Compton (SSC) spectrum of such distribution, and using the IDV source S5~0716+714 as an example, we compare it to the observed SED of S5~0716+714. The hard radio spectrum is well-fitted by this model, and the optical data can be accommodated by a power-law extension to the electron spectrum. We therefore examine the scenario of an injection of electrons that is a double power law in energy with a hard low-energy component that does not contribute to the synchrotron opacity. We show that the double power-law injection model is in good agreement with the observed SED of S5~0716+714. For intrinsic variability, we find that a Doppler factor of D\geq30 can explain the observed SED provided that low-frequency (<32 GHz) emission originates from a larger region than the higher-frequency emission. To fit the entire spectrum, D\geq65 is needed. We find the constraint imposed by induced Compton scattering at high T_B is insignificant in our model. We confirm that electron distribution with a low-energy cut-off can explain the high T_B in compact radio sources. We show that synchrotron spectrum from such distributions naturally accounts for the observed hard radio continuum with a softer optical component, without the need for an inhomogeneous source.

We present an algorithm for solving the linear dispersion relation in an inhomogeneous, magnetised, relativistic plasma. The method is a generalisation of a previously reported algorithm that was limited to the homogeneous case. The extension involves projecting the spatial dependence of the perturbations onto a set of basis functions that satisfy the boundary conditions (spectral Galerkin method). To test this algorithm in the homogeneous case, we derive an analytical expression for the growth rate of the Weibel instability for a relativistic Maxwellian distribution and compare it with the numerical results. In the inhomogeneous case, we present solutions of the dispersion relation for the relativistic tearing mode, making no assumption about the thickness of the current sheet, and check the numerical method against the analytical expression.

We set constraints on O(E/M) Lorentz Violation in QED in an effective field theory framework. A major consequence of such assumptions is the modification of the dispersion relations for electrons/positrons and photons, which in turn can affect the electromagnetic output of astrophysical objects. We compare the information provided by multiwavelength observations with a full and self-consistent computation of the broad-band spectrum of the Crab Nebula. We cast constraints of order 10^{-5} at 95% confidence level on the lepton Lorentz Violation parameters.

We investigate the growth of hydromagnetic waves driven by streaming cosmic rays in the precursor environment of a supernova remnant shock. It is known that transverse waves propagating parallel to the mean magnetic field are unstable to anisotropies in the cosmic ray distribution, and may provide a mechanism to substantially amplify the ambient magnetic field. We quantify the extent to which temperature and ionisation fractions modify this picture. Using a kinetic description of the plasma we derive the dispersion relation for a collisionless thermal plasma with a streaming cosmic ray current. Fluid equations are then used to discuss the effects of neutral-ion collisions. We calculate the extent to which the environment into which the cosmic rays propagate influences the growth of the magnetic field, and determines the range of possible growth rates. If the cosmic ray acceleration is efficient, we find that very large neutral fractions are required to stabilise the growth of the non-resonant mode. For typical supernova parameters in our galaxy, thermal effects do not significantly alter the growth rates. For weakly driven modes, ion-neutral damping can dominate over the instability at more modest ionisation fractions. In the case of a supernova shock interacting with a molecular clouds, such as in RX J1713.7-3946, with high density and low ionisation, the modes can be rapidly damped.

(abridged) Models of many astrophysical gamma-ray sources assume they contain a homogeneous distribution of electrons that are injected as a power-law in energy and evolve by interacting with radiation fields, magnetic fields and particles in the source and by escaping. This problem is particularly complicated if the radiation fields have higher energy density than the magnetic field and are sufficiently energetic that inverse Compton scattering is not limited to the Thomson regime. We present a simple, time-dependent, semi-analytical solution of the electron kinetic equation that treats both continuous and impulsive injection, cooling via synchrotron and inverse Compton radiation, (taking into account Klein-Nishina effects) and energy dependent particle escape. The kinetic equation for an arbitrary, time-dependent source function is solved by the method of Laplace transformations. Using an approximate expression for the energy loss rate that takes into account synchrotron and inverse Compton losses including Klein-Nishina effects for scattering off an isotropic photon field with either a power-law or black-body distribution, we find explicit expressions for the cooling time and escape probability of individual electrons. This enables the full, time-dependent solution to be reduced to a single quadrature. From the electron distribution, we then construct the time-dependent, multi-wavelength emission spectrum. We compare our solutions with several limiting cases and discuss the general appearance and temporal behaviour of spectral features (i.e., cooling breaks, bumps etc.). As a specific example, we model the broad-band energy spectrum of the open stellar association Westerlund-2 at different times of its evolution, and compare it with observations.

We review current theoretical ideas on pulsar winds and their surrounding nebulae. Relativistic MHD models of the wind of the aligned rotator, and of the striped wind, together with models of magnetic dissipation are discussed. It is shown that the observational signature of this dissipation is likely to be point-like, rather than extended, and that pulsed emission may be produced. The possible pulse shapes and polarisation properties are described. Particle acceleration at the termination shock of the wind is discussed, and it is argued that two distinct mechanisms must be operating, with the first-order Fermi mechanism producing the high-energy electrons (above 1 TeV) and either magnetic annihilation or resonant absorption of ion cyclotron waves responsible for the 100 MeV to 1 TeV electrons. Finally, MHD models of the morphology of the nebula are discussed and compared with observation.

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

We investigate a magnetized plasma in which injected high energy gamma-rays annihilate on a soft photon field, that is provided by the synchrotron radiation of the created pairs. For a very wide range of magnetic fields, this process involves gamma-rays between 0.3 GeV and 30 TeV. We derive a simple dynamical system for this process, analyze its stability to runaway production of soft photons and paris, and find conditions for it to automatically quench by reaching a steady state with an optical depth to photon-photon annihilation larger than unity. We discuss applications to broad-band gamma-ray emitters, in particular supermassive black holes. Automatic quenching limits the gamma-ray luminosity of these objects and predicts substantial pair loading of the jets of less active sources.

We describe an algorithm that computes the linear dispersion relation of waves and instabilities in relativistic plasmas within a Vlasov-Maxwell description. The method used is fully relativistic and involves explicit integration of particle orbits along the unperturbed equilibrium trajectories. We check the algorithm against the dispersion curves for a single component magnetised plasma and for an unmagnetised plasma with counter-streaming components in the non-relativistic case. New results on the growth rate of the Weibel or two-stream instability in a hot unmagnetised pair plasma consisting of two counter-streaming relativistic Maxwellians are presented. These are relevant to the physics of the relativistic plasmas found in gamma-ray bursts, relativistic jets and pulsar winds.

Recent observations of the polarisation of the optical pulses from the Crab pulsar motivated detailed comparative studies of the emission predicted by the polar cap, the outer gap and the two-pole caustics models. In this work, we study the polarisation properties of the synchrotron emission emanating from the striped wind model. We use an explicit asymptotic solution for the large-scale field structure related to the oblique split monopole and valid for the case of an ultra-relativistic plasma. This is combined with a crude model for the emissivity of the striped wind and of the magnetic field within the dissipating stripes themselves. We calculate the polarisation properties of the high-energy pulsed emission and compare our results with optical observations of the Crab pulsar. The resulting radiation is linearly polarised. In the off-pulse region, the electric vector lies in the direction of the projection on the sky of the rotation axis of the pulsar, in good agreement with the data. Other properties such as a reduced degree of polarisation and a characteristic sweep of the polarisation angle within the pulses are also reproduced.

The occurrence of the inverse Compton catastrophe when the synchrotron brightness temperature exceeds a threshold value, usually estimated to be 10^{12} K, appears to be in contradiction with observation because: (i) the threshold is substantially exceeded by several intra-day variable radio sources, but the inverse Compton emission is not observed, (ii) powerful, extra-galactic radio sources of known angular size do not appear to congregate close to the predicted maximum brightness temperature. We re-examine the parameter space available to synchrotron sources using a population of monoenergetic electrons, in order to see whether the revised threshold temperature is consistent with the data. The electron distribution and the population of each generation of scattered photons are computed using spatially averaged equations. We confirm our previous finding that intrinsic brightness temperatures T_{\rm B}~10^{14} K can occur without catastrophic cooling. We show that substantially higher temperatures cannot be achieved either in transitory solutions or in solutions that balance losses with a powerful acceleration mechanism. Depending on the observing frequency, we find strong cooling can set in at a range of threshold temperatures and the imposition of the additional constraint of equipartition between particle and magnetic field energy is not warranted by the data. Postulating a monoenergetic electron distribution, which approximates one that is truncated below a certain Lorentz factor, \gamma_{min}, alleviates several theoretical difficulties associated with the inverse Compton catastrophe, including those mentioned above.

We have retrieved Spitzer archive data of pre-stellar cores taken with the Multiband Imaging Photometer for Spitzer (MIPS) at a wavelength of 160 um. Seventeen images, containing eighteen cores, were constructed. Flux densities were measured for each core, and background estimates were made. Mean off-source backgrounds were found to be 48\pm10 MJy/sr in Taurus and 140\pm55 MJy/sr in Ophiuchus. Consistency was found between the MIPS 170um and ISOPHOT 160um calibration. Fourteen cores were detected both by MIPS and our previous submillimetre surveys. Spectral Energy Distributions (SEDs) were made for each core, using additional 24- and 70-um data from the Spitzer data archive, as well as previous infra-red and submillimetre data. Previous temperature estimates were refined, and new temperature estimates were made where no Infrared Space Observatory (ISO) data exist. A temperature range of 8-18K was found for the cores, with most lying in the range 10-13K. We discount recent claims that a large number of pre-stellar cores may have been misclassified and in fact contain low luminosity protostars detectable only by Spitzer. We find no new protostars in our sample other than that previously reported in L1521F. It is shown that this has a negligible effect on pre-stellar lifetime estimates.

Recently, Bell has reanalysed the problem of wave excitation by cosmic rays propagating in the pre-cursor region of a supernova remnant shock front. He pointed out a strong, non-resonant, current-driven instability that had been overlooked in the kinetic treatments, and suggested that it is responsible for substantial amplification of the ambient magnetic field. Magnetic field amplification is also an important issue in the problem of the formation and structure of relativistic shock fronts, particularly in relation to models of gamma-ray bursts. We have therefore generalised the linear analysis to apply to this case, assuming a relativistic background plasma and a monoenergetic, unidirectional incoming proton beam. We find essentially the same non-resonant instability noticed by Bell, and show that also under GRB conditions, it grows much faster than the resonant waves. We quantify the extent to which thermal effects in the background plasma limit the maximum growth rate.

We performed spectral line observations of CO J=2-1, 13CO J=1-0, and C18O J=1-0 and polarimetric observations in the 1.3 mm continuum and CO J=2-1 toward a multiple protostar system, L1448 IRS 3, in the Perseus molecular complex at a distance of ~250 pc, using the BIMA array. In the 1.3 mm continuum, two sources (IRS 3A and 3B) were clearly detected with estimated envelope masses of 0.21 and 1.15 solar masses, and one source (IRS 3C) was marginally detected with an upper mass limit of 0.03 solar masses. In CO J=2-1, we revealed two outflows originating from IRS 3A and 3B. The masses, mean number densities, momentums, and kinetic energies of outflow lobes were estimated. Based on those estimates and outflow features, we concluded that the two outflows are interacting and that the IRS 3A outflow is nearly perpendicular to the line of sight. In addition, we estimated the velocity, inclination, and opening of the IRS 3B outflow using Bayesian statistics. When the opening angle is ~20 arcdeg, we constrain the velocity to ~45 km/s and the inclination angle to ~57 arcdeg. Linear polarization was detected in both the 1.3 mm continuum and CO J=2-1. The linear polarization in the continuum shows a magnetic field at the central source (IRS 3B) perpendicular to the outflow direction, and the linear polarization in the CO J=2-1 was detected in the outflow regions, parallel or perpendicular to the outflow direction. Moreover, we comprehensively discuss whether the binary system of IRS 3A and 3B is gravitationally bound, based on the velocity differences detected in 13CO J=1-0 and C18O J=1-0 observations and on the outflow features. The specific angular momentum of the system was estimated as ~3e20 cm^2/s, comparable to the values obtained from previous studies on binaries and molecular clouds in Taurus.

We have mapped linearly polarized dust emission from the prestellar cores L1498 and L1517B with the James Clerk Maxwell Telescope (JCMT) using the Submillimetre Common User Bolometer Array (SCUBA) and its polarimeter SCUBAPOL at a wavelength of 850um. We use these measurements to determine the plane-of-sky magnetic field orientation in the cores. In L1498 we see a magnetic field across the peak of the core that lies at an offset of 19 degrees to the short axis of the core. This is similar to the offsets seen in previous observations of prestellar cores. To the southeast of the peak, in the filamentary tail of the core, we see that the magnetic field has rotated to lie almost parallel to the long axis of the filament. We hypothesise that the field in the core may have decoupled from the field in the filament that connects the core to the rest of the cloud. We use the Chandrasekhar-Fermi (CF) method to measure the plane-of-sky field strength in the core of L1498 to be 10 +/- 7 uG. In L1517B we see a more gradual turn in the field direction from the northern part of the core to the south. This appears to follow a twist in the filament in which the core is buried, with the field staying at a roughly constant 25 degree offset to the short axis of the filament, also consistent with previous observations of prestellar cores. We again use the CF method and calculate the magnetic field strength in L1517B also to be 30 +/- 10 uG. Both cores appear to be roughly virialised. Comparison with our previous work on somewhat denser cores shows that, for the denser cores, thermal and non-thermal (including magnetic) support are approximately equal, while for the lower density cores studied here, thermal support dominates.

Some rapidly variable extra-galactic radio sources show very high brightness temperatures T_B>10^{12}K and high degrees of circular polarisation (1%). Standard synchrotron models that assume a power-law electron distribution cannot produce such high temperatures and have much lower degrees of intrinsic circular polarisation. We examine the synchrotron and inverse Compton radiation from a monoenergetic electron distribution using standard synchrotron theory. Constraints on the source parameters are found by formulating the results as functions of the source size, Doppler boosting factor, optical depth to synchrotron self-absorption, maximum frequency of synchrotron emission, and the strength of the inverse Compton radiation. The model gives brightness temperatures T_B=10^{13}K to 10^{14}K for moderate (<10) Doppler boosting factors and intrinsic degrees of circular polarisation at the percent level. It predicts a spectrum I_\nu\propto\nu^{1/3} between the radio and the infra-red as well as emission in the MeV to GeV range. We find the conditions under which electrons do not cool within the source, enabling the GHz emission to emerge without absorption and the potentially catastrophic energy losses by inverse Compton scattering to be avoided. We suggest that sources such as PKS 1519 -273, PKS 0405 -385 and J 1819 +3845 can be understood within this scenario without invoking high Doppler boosting factors, coherent emission mechanisms, or the dominance of proton synchrotron radiation.

The high-energy gamma-ray emission discovered using the H.E.S.S. telescopes from the binary system PSR B1259 -63, is modelled using an extension of the approach that successfully predicted it. We find that the simultaneous INTEGRAL and H.E.S.S. data permit both a model with dominant radiative losses, high pulsar wind Lorentz factor and modest efficiency as well as one with dominant adiabatic losses, a slower wind and higher efficiency. Additional, simultaneous, X-ray and TeV data sets are needed to lift this degeneracy.

The high-energy gamma-ray light curve of the binary system PSR B1259 -63, is computed using the approach that successfully predicted the spectrum at periastron. The simultaneous INTEGRAL and H.E.S.S. spectra taken 16 days after periastron currently permit both a model with dominant radiative losses, high pulsar wind Lorentz factor and modest efficiency as well as one with dominant adiabatic losses, a slower wind and higher efficiency. In this paper we shown how the long-term light curve may help to lift this degeneracy.

Two-dimensional, relativistic, MHD simulations of pulsar-wind powered nebulae provide strong constraints on the properties of the winds themselves. In particular, they confirm that Poynting flux must be converted into particle energy close to or inside the termination shock front, emphasising the puzzle known as the $\sigma$ paradox. To distinguish between the different possible resolutions of this paradox, additional observational constraints are required. In this paper, I briefly discuss two recents developments in this respect: the modelling of high time-resolution optical polarimetry of the Crab pulsar, and the detection of the pulsar/Be star binary PSR 1259-63 in TeV energy gamma-rays.

The Stokes parameters of the pulsed synchrotron radiation produced in the striped pulsar wind model are computed and compared with optical observations of the Crab pulsar. We assume the main contribution to the wind emissivity comes from a thin transition layer where the dominant toroidal magnetic field reverses its polarity. The radial component of the field is neglected, but a small meridional component is added. The resulting radiation is linearly polarized (Stokes V=0). In the off-pulse region, the electric vector lies in the direction of the projection on the sky of the rotation axis of the pulsar. This property is unique to the wind model and in good agreement with the data. Other properties such as a reduced degree of polarization and a characteristic sweep of the polarization angle within the pulses are also reproduced. These properties are qualitatively unaffected by variations of the wind Lorentz factor, the electron injection power law index and the inclination of the line of sight.

Observations have been carried out with SCUBA at the JCMT of 52 molecular cloud cores that do not contain any sign of protostellar activity. These are all therefore candidate prestellar cores, which are believed to represent the stage of star formation that precedes the formation of a protostar. 29 of the 52 cores were detected at 850 microns at varying levels of signal-to-noise ratio greater than 3 sigma at peak. The detected cores were split into 'bright' cores and `intermediate' cores, depending on their peak flux density at 850 microns. Cores with peak 850 microns flux densities greater than 170 mJy/beam were designated 'bright' (13 cores), while those flux densities below this value were designated 'intermediate' (16 cores). This dividing line corresponds to A_v~50 under typical assumptions. The data are combined with our previously published ISO data, and the physical parameters of the cores, such as density and temperature, are calculated. Detailed fitting of the bright core radial profiles shows that they are not critical Bonnor-Ebert spheres, in agreement with previous findings. However, we find that intermediate cores, such as B68 (which has previously been claimed to be a Bonnor-Ebert sphere), may in fact be consistent with the Bonnor-Ebert criterion, suggesting perhaps that cores pass through such a phase during their evolution. We make rough estimates of core lifetimes based on the statistics of detections and find that the lifetime of a prestellar core is roughly ~3x10^5 years, while that of a bright core is \~1.5x10^5 years. Comparisons with some magnetic and turbulence regulated collapse models show that no model can match all of the data. Models that are tuned to fit the total prestellar core lifetime, do not predict the relative numbers of cores seen at each stage.

I review the constraints placed on relativistic pulsar winds by comparing optical and X-ray images of the inner Crab Nebula on the one hand with two-dimensional MHD simulations on the other. The various proposals in the literature for achieving the low magnetisation required at the inner edge of the Nebula, are then discussed, emphasising that of dissipation in the striped-wind picture. The possibility of direct observation of the wind is examined. Based on the predicted orientation of the polarisation vector, I outline a new argument suggesting that the off-pulse component of the optical emission of the Crab pulsar originates in the wind.

A property common to several different astrophysical sources of high-energy gamma-rays is the presence of bulk motion at relativistic speed. The intrinsic spectra of the nonthermal radiating particles also show interesting similarities, with a pronounced hardening towards lower energies. This suggests two distinct acceleration mechanisms could be at work in these sources. At high energies, the stochastic first-order Fermi process at shocks seems to provide a reasonable explanation. I will briefly review the status of this mechanism before discussing the possibility that, at lower energies, non-stochastic acceleration in the induced electric field of a relativistic current sheet plays a role.

We present calculations of the spectral and temporal radiative signatures expected from ultrarelativistic protons in compact sources. The coupling between the protons and the leptonic component is assumed to occur via Bethe-Heitler pair production. This process is treated by modeling the results of Monte-Carlo simulations and incorporating them in a time-dependent kinetic equation, that we subsequently solve numerically. Thus, the present work is, in many respects, an extension of the leptonic `one-zone' models to include hadrons. Several examples of astrophysical importance are presented, such as the signature resulting from the cooling of relativistic protons on an external black-body field and that of their cooling in the presence of radiation from injected electrons. We also investigate and refine the threshold conditions for the 'Pair Production/Synchrotron' feedback loop which operates when relativistic protons cool efficiently on the synchrotron radiation of the internally produced Bethe-Heitler pairs. We demonstrate that an additional component of injected electrons lowers the threshold for this instability.

Relativistic current sheets have been proposed as the sites of dissipation in pulsar winds, jets in active galaxies and other Poynting-flux dominated flows. It is shown that the steady versions of these structures differ from their nonrelativistic counterparts because they do not permit transformation to a de Hoffmann/Teller reference frame, in which the electric field vanishes. Instead, their generic form is that of a true neutral sheet: one in which the linking magnetic field component normal to the sheet is absent. Taken together with Alfven's limit on the total cross-field potential, this suggests plasma is ejected from the sheet in the cross-field direction rather than along it. The maximum energy to which such structures can accelerate particles is derived, and used to compute the maximum frequency of the subsequent synchrotron radiation. This can be substantially in excess of standard estimates. In the magnetically driven gamma-ray burst scenario, acceleration of electrons is possible to energies sufficient to enable photon-photon pair production after an inverse Compton scattering event.

We have mapped linearly polarized dust emission from L183 with the JCMT SCUBA polarimeter and have analyzed these and our previously published data for the prestellar cores L183, L1544, and L43 in order to estimate magnetic field strengths in the plane of the sky, $B_{pos}$. The analysis used the Chandrasekhar-Fermi technique, which relates the dispersion in polarization position angles to $B_{pos}$. We have used these estimates of the field strengths (neglecting the unmeasured line-of-sight component) to find the mass-to-magnetic flux ratios $\lambda$ (in units of the critical ratio for magnetic support). Results are $B_{pos} \approx 80$ $\mu$G and $\lambda \approx 2.6$ for L183, $B_{pos} \approx 140$ $\mu$G and $\lambda \approx 2.3$ for L1544, and $B_{pos} \approx 160$ $\mu$G and $\lambda \approx 1.9$ for L43. Hence, without correction for geometrical biases, for all three cores the mass-to-flux ratios are supercritical by a factor of $\sim 2$, and magnetic support cannot prevent collapse. However, a statistical mean correction for geometrical bias may be up to a factor of three; this correction would reduce the individual $\lambda$'s to $\lambda_{cor} \approx 0.9$, 0.8, and 0.6, respectively; these values are approximately critical or slightly subcritical. These data are consistent with models of star formation driven by ambipolar diffusion in a weakly turbulent medium, but cannot rule out models of star formation driven by turbulence.

A brief summary is presented of our current knowledge of the structure of cold molecular cloud cores that do not contain protostars, sometimes known as starless cores. The most centrally condensed starless cores are known as pre-stellar cores. These cores probably represent observationally the initial conditions for protostellar collapse that must be input into all models of star formation. The current debate over the nature of core density profiles is summarised. A cautionary note is sounded over the use of such profiles to ascertain the equilibrium status of cores. The magnetic field structure of pre-stellar cores is also discussed.

The conversion of the Crab pulsar wind from one dominated by Poynting flux close to the star to one dominated by particle-born energy at the termination shock is considered. The idea put forward by Coroniti (1990) and criticised by Lyubarsky & Kirk (2001) that reconnection in a striped wind is responsible, is generalised to include faster prescriptions for the a priori unknown dissipation rate. Strong acceleration of the wind is confirmed, and the higher dissipation rates imply complete conversion of Poynting flux into particle-born flux within the unshocked wind.

Flows in which energy is transported predominantly as Poynting flux are thought to occur in pulsars, gamma-ray bursts and relativistic jets from compact objects. The fluctuating component of the magnetic field in such a flow can in principle be dissipated by magnetic reconnection, and used to accelerate the flow. We investigate how rapidly this transition can take place, by implementing into a global MHD model, that uses a thermodynamic description of the plasma, explicit, physically motivated prescriptions for the dissipation rate: a lower limit on this rate is given by limiting the maximum drift speed of the current carriers to that of light, an upper limit follows from demanding that the dissipation zone expand only subsonically in the comoving frame and a further prescription is obtained by assuming that the expansion speed is limited by the growth rate of the relativistic tearing mode. In each case, solutions are presented which give the Lorentz factor of a spherical wind containing a transverse, oscillating magnetic field component as a function of radius. In the case of the Crab pulsar, we find that the Poynting flux can be dissipated before the wind reaches the inner edge of the Nebula if the pulsar emits electron positron pairs at a rate >1.E40 per second, thus providing a possible solution to the sigma-problem.

We discuss the results of modelling the TeV gamma-ray spectra of two AGNs, Mkn 501 and Mkn 421 that have almost the same redshifts: z=0.031 and z=0.034, respectively. The effect of intergalactic gamma-ray absorption is treated as an uncertainty in the measurement of the intrinsic spectrum. Although the objects differ, we obtain satisfactory fits for both of them in a synchrotron self-Compton scenario. Compared to previous models, our fits are characterised by higher values of the Doppler factor (>= 50) and an electron injection spectrum extending to higher energies (Gmax = 1.5x10^5). In the case of Mkn 421, the observed difference in spectral slope in X-rays and TeV gamma-rays between the high and low states can be explained as a variation of a single parameter - the maximum energy Gmax mc^2 at which electrons are injected.

A solution to the sigma problem - that of finding a mechanism capable of converting Poynting energy flux to particle-borne energy flux in a pulsar wind - was proposed several years ago by Coroniti and Michel who considered a particular prescription for magnetic reconnection in a striped wind. This prescription was later shown to be ineffective. In this paper, we discuss the basic microphysics of the reconnection process and conclude that a more rapid prescription is permissible. Assuming dissipation to set in at some distance outside the light-cylinder, we compute the resulting radiation signature and find that the synchrotron emission of heated particles appears periodic, in general showing both a pulse and an interpulse. The predicted spacing of these agrees well with observation in the case of the Crab and Vela pulsars. Using parameters appropriate for the Crab pulsar - magnetization parameter at the light cylinder sigma_L = 6 x 10^4, Lorentz factor Gamma=250 - reasonable agreement is found with the observed total pulsed luminosity. This suggest that the high-energy pulses from young pulsars originate not in the co-rotating magnetosphere within the light cylinder (as in all other models) but from the radially directed wind well outside it.

The spectral evolution of powerful double radio galaxies (FR II's) is thought to be determined by the acceleration of electrons at the termination shock of the jet, their transport through the bright head region into the lobes and the production of the radio emission by synchrotron radiation in the lobes. Models presented to date incorporate some of these processes in prescribing the electron distribution which enters the lobes. We have extended these models to include a description of electron acceleration at the relativistic termination shock and a selection of transport models for the head region. These are coupled to the evolution of the electron spectrum in the lobes under the influence of losses due to adiabatic expansion, by inverse Compton scattering on the cosmic background radiation and by synchrotron radiation. The evolutionary tracks predicted by this model are compared to observation using the power/source-size (P-D) diagram. We find that the simplest scenario, in which accelerated particles suffer adiabatic losses in the head region which become more severe as the source expands produces P-D-tracks which conflict with observation, because the power is predicted to decline too steeply with increasing size. Agreement with observation can be found by assuming that adiabatic losses are compensated during transport between the termination shock and the lobe by a re-acceleration process distributed throughout the head region.

The radiation of a pulsar wind is computed assuming that at roughly 10 to 100 light cylinder radii from the star, magnetic energy is dissipated into particle energy. The synchrotron emission of heated particles appears periodic, with, in general, both a pulse and an interpulse. The predicted spacing agrees well with the Crab and Vela pulse profiles.Using parameters appropriate for the Crab pulsar (magnetisation parameter at the light cylinder $\sigma=6\times10^4$, Lorentz factor $\Gamma=250$) agreement is found with the observed total pulsed luminosity. This suggests that the high-energy pulses from young pulsars originate not in the corotating magnetosphere within the light cylinder (as in all other models) but from the radially directed wind well outside it.

The spin-down power of a pulsar is thought to be carried away in an MHD wind in which, at least close to the star, the energy transport is dominated by Poynting flux. The pulsar drives a low-frequency wave in this wind, consisting of stripes of toroidal magnetic field of alternating polarity, propagating in a region around the equatorial plane. The current implied by this configuration falls off more slowly with radius than the number of charged particles available to carry it, so that the MHD picture must, at some point, fail. Recently, magnetic reconnection in such a structure has been shown to accelerate the wind significantly. This reduces the magnetic field in the comoving frame and, consequently, the required current, enabling the solution to extend to much larger radius. This scenario is discussed and, for the Crab Nebula, the range of validity of the MHD solution is compared with the radius at which the flow appears to terminate. For sufficiently high particle densities, it is shown that a low frequency entropy wave can propagate out to the termination point. In this case, the "termination shock" itself must be responsible for dissipating the wave.

We review recent theoretical developments on pulsar winds, their nebulae and relativistic shock acceleration, and show how they illuminate unsolved problems in plerion spectra, in particular the multiple spectral breaks in the Crab and the low-frequency breaks of plerions such as G 21.5-0.9 and 3C 58. Recent work on Fermi acceleration theory at relativistic shocks shows that a particle spectral index of 2.2-2.3, compatible with the X-ray spectra of plerions, results under a wide variety of assumptions. If pulsar winds contain ions as well as electrons and positrons, the mechanism of Hoshino et al. (1992), which yields harder spectra, would operate at lower energies and may explain the flat radio spectral indices of plerions. This scenario implies wind parameters in the Crab compatible with the pulsar wind acceleration model of Lyubarsky & Kirk (2001). Recent hydrodynamical simulations of plerion evolution inside SNR blast waves demonstrate that the passage of the reverse shock rapidly compresses the plerion. Using a simple isobaric model, we investigate the influence of the resulting magnetic field compression and decrease in shock radius on the evolution of the plerion spectrum. We suggest that the passage of the reverse shock may explain the low-frequency breaks in 3C 58 and G 21.5-0.9, as well as the increase in 3C 58's radio flux.

Galactic binary systems that contain a black hole candidate emit hard X-rays in their low luminosity mode. We show that this emission can be understood as due to the Compton scattering of photons from the companion star and/or the accretion disk by relativistic electrons in a jet. The same electrons are also responsible for the radio emission. Two sources -- XTE J1118+480 and Cygnus X-1 -- are modelled as representatives of black holes with low and high luminosity companion stars respectively. We further show that the ultraluminous compact X-ray sources observed in nearby galaxies have the properties expected of stellar mass black holes with high luminosity companions in which the jet is oriented close to our line of sight.

ISO data taken with the long-wavelength imaging photo-polarimeter ISOPHOT are presented of 18 pre-stellar cores at three far-infrared wavelengths - 90, 170 and 200 microns. Most of the cores are detected clearly at 170 and 200 but only one is detected strongly at 90 microns, indicating that mostly they are very cold, with typical temperatures of only 10-20K. Colour temperature images are constructed for each of the cores. Most of the cores are seen to be either isothermal, or to have associated temperature gradients from the core centres to their edges, with all except one being cooler at the centre. We compare the data with previous ISOCAM absorption data and calculate the energy balance for those cores in common between the two samples. We find that the energy radiated by each core in the far-IR is similar to that absorbed at shorter wavelengths. Hence there is no evidence for a central heating source in any of the cores - even those for which previous evidence for core contraction exists. This is all consistent with external heating of the cores by the local interstellar radiation field, confirming their pre-stellar nature.

We consider the acceleration of charged particles near ultra-relativistic shocks, with Lorentz factor Gamma_s >> 1. We present simulations of the acceleration process and compare these with results from semi-analytical calculations. We show that the spectrum that results from acceleration near ultra-relativistic shocks is a power law, N(E) \propto E^{-s}, with a nearly universal value s \approx 2.2 - 2.3 for the slope of this power law. We confirm that the ultra-relativistic equivalent of Fermi acceleration at a shock differs from its non-relativistic counterpart by the occurence of large anisotropies in the distribution of the accelerated particles near the shock. In the rest frame of the upstream fluid, particles can only outrun the shock when their direction of motion lies within a small loss cone of opening angle theta_c \approx 1/Gamma_s around the shock normal. We also show that all physically plausible deflection or scattering mechanisms can change the upstream flight direction of relativistic particles originating from downstream by only a small amount: Delta theta ~ 1/Gamma_s. This limits the energy change per shock crossing cycle to Delta E ~ E, except for the first cycle where particles originate upstream. In that case the upstream energy is boosted by a factor ~ Gamma_s^2 for those particles that are scattered back across the shock into the upstream region.

By including Klein-Nishina effects, we generalize previous calculations of the beaming pattern of photons produced by inverse Compton scattering. For an isotropic distribution of soft photons upscattered by nonthermal electrons with a power-law density distribution embedded in a plasma moving with relativistic bulk speed we show that the observed radiation intensity is proportional to D^(3+p), where D is the Doppler boosting factor, and p the electron distribution index. This agrees with previous computations performed in the Thomson limit, where the observed spectral index is a=(p-1)/2 and the beaming pattern D^(4+2a). Independent of D, Klein-Nishina effects limit the location of the peak energy of the observed spectral energy distribution such that e_{peak} less than 1/e_0, where e_0 is the energy of the seed photons in units of mc^2. Assuming that the seed photons originate in the broad line region, we demonstrate that the GeV emission of blazars is significantly modified by Klein-Nishina effects, the spectrum being softer than that calculated in the Thomson limit. We further show that the change in spectral index of the inverse Compton emission across the peak energy can exceed the value of 0.5 predicted by computations performed in the Thomson limit. The model spectra agree with OSSE and COMPTEL limits on this break without invoking the effects of differential absorption at the edge of a gamma-ray photosphere.

Motivated by recent unsuccessful efforts to detect the predicted flux of TeV gamma-rays from supernova remnants, we present a critical examination of the theory on which these predictions are based. Three crucial problems are identified: injection, maximum achievable particle energy and spectral index. In each case significant new advances in understanding have been achieved, which cast doubt on prevailing paradigms such as Bohm diffusion and single-fluid MHD. This indicates that more realistic analytical models, backed by more sophisticated numerical techniques should be employed to obtain reliable predictions. Preliminary work on incorporating the effects of anomalous transport suggest that the resulting spectrum should be significantly softer than that predicted by conventional theory.

We extend the eigenfunction method of computing the power-law spectrum of particles accelerated at a relativistic shock fronts to apply to shocks of arbitrarily high Lorentz factor. In agreement with the findings of Monte-Carlo simulations, we find the index of the power-law distribution of accelerated particles which undergo isotropic diffusion in angle at an ultrarelativistic, unmagnetized shock is s=4.23 (where s=-d(ln f)/dp with f the Lorentz invariant phase-space density and p the momentum). This corresponds to a synchrotron index for uncooled electrons of a=0.62 (taking cooling into account a=1.12), where a=-d(ln F)/dn, F is the radiation flux and n the frequency. We also present an approximate analytic expression for the angular distribution of accelerated particles, which displays the effect of particle trapping by the shock: compared with the non-relativistic case the angular distribution is weighted more towards the plane of the shock and away from its normal. We investigate the sensitivity of our results to the transport properties of the particles and the presence of a magnetic field. Shocks in which the ratio of Poynting to kinetic energy flux upstream is not small are less compressive and lead to larger values of $s$.

We present preliminary results of our work on blazar unification. We assume that all blazars have a broad line region (BLR) and that the size of the BLR scales with the power of the source in a manner similar to that derived through reverberation mapping in radio quiet active galactic nuclei (AGNs). Using a self-consistent emission model that includes particle acceleration we show that according to this scaling, in weak sources like MKN 421, the inverse Compton (IC) scattering losses are dominated by synchrotron-self Compton scattering (SSC), while in powerful sources, like 3C 279, they are dominated by external Compton (EC) scattering of BLR photons. In agreement with other workers, we show that even in the powerful sources that are dominated by EC scattering, the hard X-ray emission is due to SSC. Finally, we show that this scaling reproduces well the observed sequence of blazar properties with luminosity.

It is generally thought that most of the spin-down power of a pulsar is carried away in an MHD wind dominated by Poynting flux. In the case of an oblique rotator, a significant part of this energy can be considered to be in a low-frequency wave, consisting of stripes of toroidal magnetic field of alternating polarity, propagating in a region around the equatorial plane. Magnetic reconnection in such a structure has been proposed as a mechanism for transforming the Poynting flux into particle energy in the pulsar wind. We have re-examined this process and conclude that the wind accelerates significantly in the course of reconnection. This dilates the timescale over which the reconnection process operates, so that the wind requires a much larger distance than was previously thought in order to convert the Poynting flux to particle flux. In the case of the Crab, the wind is still Poynting-dominated at the radius at which a standing shock is inferred from observation. An estimate of the radius of the termination shock for other pulsars implies that all except the milli-second pulsars have Poynting-flux dominated winds all the way out to the shock front.

We present the first published maps of magnetic fields in pre-stellar cores, to test theoretical ideas about the way in which the magnetic field geometry affects the star formation process. The observations are JCMT-SCUBA maps of 850 micron thermal emission from dust. Linear polarizations at typically ten or more independent positions in each of three objects, L1544, L183 and L43 were measured, and the geometries of the magnetic fields in the plane of the sky were mapped from the polarization directions. The observed polarizations in all three objects appear smooth and fairly uniform. In L1544 and L183 the mean magnetic fields are at an angle of around 30 degrees to the minor axes of the cores. The L43 B-field appears to have been influenced in its southern half, such that it is parallel to the wall of a cavity produced by a CO outflow from a nearby T Tauri star, whilst in the northern half the field appears less disturbed and has an angle of 44 degrees to the core minor axis. We briefly compare our results with published models of magnetized cloud cores and conclude that no current model can explain these observations simultaneously with previous ISOCAM data.

We consider particle acceleration to high energy via diffusive shock acceleration in a simple, self-consistent shock in jet model for blazars. Electrons are assumed to be accelerated at a shock front in relativistic jets and radiate synchrotron emission in a post-shock region. The full time, space and momentum dependence of the electron distribution function is used for a calculation of the nonthermal synchrotron spectra. We discuss the evolution of the spectral index by varying the rate at which particles enter the acceleration process. The results indicate that the synchrotron spectral index displays a characteristic looplike behaviour with intensity (as has been observed in several blazars), where the orientation of the loop depends on whether the acceleration time scale is comparable to the synchrotron cooling time scale or not. We show that our model provides a good fit to the observed evolution of the spectral index of Mkn 421 during a flare in 1994.

We adapt and modify the eigenfunction method of computing the power-law spectrum of particles accelerated at a relativistic shock front via the first-order Fermi process (Kirk, J.G., Schneider, P., Astrophysical Journal 315, 425 (1987)) to apply to shocks of arbitrarily high Lorentz factor. The power-law index of accelerated particles undergoing isotropic small-angle scattering at an ultrarelativistic, unmagnetized shock is found to be s=4.23 +/- 0.2 (where s=d\ln f/ d\ln p, with f the Lorentz-invariant phase-space density and p the momentum), in agreement with the results of Monte-Carlo simulations. We present results for shocks in plasmas with different equations of state and for Lorentz factors ranging from 5 to infinity.

We examine Fermi-type acceleration at relativistic shocks, and distinguish between the initial boost of the first shock crossing cycle, where the energy gain per particle can be very large, and the Fermi process proper with repeated shock crossings, in which the typical energy gain is of order unity. We calculate by means of numerical simulations the spectrum and angular distribution of particles accelerated by this Fermi process, in particular in the case where particle dynamics can be approximated as small-angle scattering. We show that synchrotron emission from electrons or positrons accelerated by this process can account remarkably well for the observed power-law spectra of GRB afterglows and Crab-like supernova remnants. In the context of a decelerating relativistic fireball, we calculate the maximum particle energy attainable by acceleration at the external blast wave, and discuss the minimum energy for this acceleration process and its consequences for the observed spectrum.

We investigate the effects of inverse Compton scattering by electrons and positrons in the unshocked winds of rotationally-powered binary pulsars. This process can scatter low energy target photons to produce gamma rays with energies from MeV to TeV. The binary radio pulsars PSR B1259-63 and PSR J0045-73 are both in close eccentric orbits around bright main sequence stars which provide a huge density of low energy target photons. The inverse Compton scattering process transfers momentum from the pulsar wind to the scattered photons, and therefore provides a drag which tends to decelerate the pulsar wind. We present detailed calculations of the dynamics of a pulsar wind which is undergoing inverse Compton scattering, showing that the deceleration of the wind of PSR B1259-63 due to `inverse Compton drag' is small, but that this process may confine the wind of PSR J0045-73 before it attains pressure balance with the outflow of its companion star. We calculate the spectra and light curves of the resulting inverse Compton emission from PSR B1259-63 and show that if the size of the pulsar wind nebula is comparable to the binary separation, then the gamma-ray emission from the unshocked wind may be detectable by atmospheric Cerenkov detectors or by the new generation of satellite-borne gamma-ray detectors such as INTEGRAL and GLAST. This mechanism may therefore provide a direct probe of the freely-expanding regions of pulsar winds, previously thought to be invisible.

Supernova remnants are expected to contain braided (or stochastic) magnetic fields, which are in some regions directed mainly perpendicular to the shock normal. For particle acceleration due to repeated shock crossings, the transport in the direction of the shock normal is crucial. The mean squared deviation along the shock normal is then proportional to the square root of the time. This kind of anomalous transport is called sub-diffusion. We use a Monte-Carlo method to examine this non-Markovian transport and the acceleration. As a result of this simulation we are able to examine the propagator, density and pitch-angle distribution of accelerated particles, and especially the spectral properties. These are in broad agreement with analytic predictions for both the sub-diffusive and the diffusive regimes, but the steepening of the spectrum predicted when changing from diffusive to sub-diffusive transport is found to be even more pronounced than predicted.

The present work considers diffusive shock acceleration at non-relativistic shocks using a system of stochastic differential equations (SDE) equivalent to the Fokker-Planck equation. We compute approximate solutions of the transport of cosmic particles at shock fronts with a SDE numerical scheme. The momentum gain is given by implicit calculations of the fluid velocity gradients using a linear interpolation between two consecutive time steps. We validate our procedure in the case of single shock acceleration with different shock thickness, with or without synchrotron losses. A comparative discussion of implicit and explicit schemes for different shock thickness shows that implicit calculations extend the range of applicability of SDE schemes to infinitely thin 1D shocks. The method is then applied to multiple shock acceleration for a system of identical shocks which free parameters are the distance between two consecutive shocks, the synchrotron losses time and the escape time of the particles. The stationary distribution reproduces quite well the flat differential logarithm energy distribution produced by multiple shock effect, and also the piling-up effect due synchrotron losses at a momentum where they equilibrate the acceleration rate. At higher momenta particle losses dominate and the spectrum drops. The competition between acceleration and loss effects leads to a pile-up shaped distribution which appears to be effective only in a restrict range of inter-shock distances of 10-100 diffusion lengths. We finally compute the optically thin synchrotron spectrum produced such periodic pattern which can explain flat and/or inverted spectra observed in Flat Radio spectrum Quasars and in the galactic centre.

Observations of both gamma-ray burst sources and certain classes of active galaxy indicate the presence of relativistic shock waves and require the production of high energy particles to explain their emission. In this paper we review the basic theory of shock waves in relativistic hydrodynamics and magneto-hydrodynamics, emphasising the astrophysically interesting cases. This is followed by an overview of the theory of particle acceleration at such shocks. We summarise the applications to the astrophysics of relativistic jets and fireball models of gamma-ray-bursts.

The recent HEGRA observations of the blazar Mkn 501 show strong curvature in the very high energy gamma-ray spectrum. Applying the gamma-ray opacity derived from an empirically based model of the intergalactic infrared background radiation field (IIRF), to these observations, we find that the intrinsic spectrum of this source is consistent with a power-law: dN/dE~ E^-alpha with alpha=2.00 +/- 0.03 over the range 500 GeV - 20 TeV. Within current synchrotron self-Compton scenarios, the fact that the TeV spectral energy distribution of Mkn 501 does not vary with luminosity, combined with the correlated, spectrally variable emission in X-rays, as observed by the BeppoSAX and RXTE instruments, also independently implies that the intrinsic spectrum must be close to alpha=2. Thus, the observed curvature in the spectrum is most easily understood as resulting from intergalactic absorption.

We present the expected variability features in the context of the Synchro- Self-Compton (SSC) model of emission from Active Galactic Nuclei (AGNs) and show that the homogeneous SSC model describes well the observed multi- wavelength spectrum and variability of objects such as Mkn 421 and Mkn 501; however in order to explain the very fast flaring timescales occasionally observed from these objects one should use a laminar, rather than spherical, source geometry. This picture leads naturally to a shock-in- jet model which we approach as a site of diffusive particle acceleration. We show that this acceleration scheme can explain certain characteristic features of AGN X-ray spectra and provides us with further observational tests.

In the theory of diffusive acceleration at oblique shock fronts the question of the existence of a discontinuity of energetic particle density is contentious. The resolution of this problem is interesting from a theoretical point of view, and potentially for the interpretation of observations of particle densities at heliospheric shocks and of high-resolution radio observations of the rims of supernova remnants. It can be shown analytically that an isotropic particle distribution at a shock front implies continuity of the particle density -- whether or not the shock is oblique. However, if the obliquity of the shock induces an anisotropy, a jump is permitted. Both semi-analytic computations and Monte-Carlo simulations are used to show that, for interesting parameter ranges, a jump is indeed produced, with accelerated particles concentrated in a precursor ahead of the shock front.

We derive light curves for the hard gamma-ray emission, at energies up to several TeV, expected from the unique pulsar/Be-star binary system PSR B1259-63. This is the only known system in our galaxy in which a radio pulsar is orbiting a main sequence star. We show that inverse Compton emission from the electrons and positrons in the shocked pulsar wind, scattering target photons from the Be star, produces a flux of hard gamma-rays that should be above the sensitivity threshold of present day atmospheric Cerenkov detectors. Furthermore, we predict that the flux of hard gamma-rays produced via this mechanism has a characteristic variation with orbital phase that should be observable, and which is not expected from any other mechanism.

We model the acceleration of electrons at a shock front in a relativistic blazar jet and compute the radiation they emit in a post-shock region which contains a homogeneous magnetic field. The full space, time and momentum dependence of the electron distribution is used in this calculation. It is shown that the "homogeneous" synchrotron model is recovered, provided the downstream speed of the plasma away from the shock front is nonrelativistic, and provided that the light travel times across the face of the shock front is unimportant. By varying the rate at which particles are picked up by the acceleration process, we calculate the time-dependence of the spectra. Since the magnetic field strength is assumed constant within the emission region, each frequency band can be identified with electrons of a particular energy. We find that for a band in which the electrons are accelerated rapidly compared to the rate at which they cool, the spectra typically harden during phases of rising flux, and soften during phases of falling flux, as has been observed in the objects PKS 2155-304 and Mkn 421. However, in a frequency band in which the timescales are comparable, the reverse behaviour is to be expected. We discuss the extent to which observations of both the stationary spectrum and the spectral variability of the synchrotron component of blazar emission can be used to constrain the model.

The theory of particle acceleration at shock fronts is briefly reviewed, with special emphasis on the production of the particles responsible for the nonthermal emission from blazars. The flat radio/IR spectra of these sources cannot be produced by diffusive acceleration at a simple nonrelativistic shock front propagating in a homogeneous medium. It can, however, be produced by a single unmodified mildly relativistic shock, if the pressure in the shocked gas is provided by the leptonic component, or, independently of the equation of state, by a relativistic shock which is oblique to the magnetic field. The analytic theory of these shocks makes several simplifications, but Monte-Carlo simulations exist which extend the range of validity. Of particular interest is acceleration in a tangled magnetic field. Here, however, the Monte-Carlo simulations have not yet yielded unambiguous results. The "homogeneous" models of blazar emission are discussed, and is it shown that they imply a geometry of the emitting region which is laminar in form, with an aspect ratio of d/R< 3% in the case of Mkn~421. Identifying these with relativistic shock fronts, a model of acceleration is described, which displays characteristic variations in the synchrotron spectral index with intensity.

Anomalous transport processes in which the variance of the distance travelled does not necessarily increase linearly with time are modelled using the formalism of continuous time random walks. We compute particle propagators which have the required dependence on space and time and use these to find the spatial dependence of the synchrotron radiation emitted by a population of continuously injected electrons. As the electrons are transported away from the source they cool, and the synchrotron spectrum softens. Sub-diffusive transport -- corresponding to stochastic trapping, or restriction of the transport across the average direction of a stochastic magnetic field -- produces a much slower rate of change of spectral index than does supra-diffusion -- which occurs when particles move almost without scattering, in a field containing large ordered regions. Application to the diffuse emission of the outer parts of the Coma cluster favours an interpretation in terms of supra-diffusion.

We present a semi-analytic approach to solving the Boltzmann equation describing the comptonisation of low frequency input photons by a thermal distribution of electrons in the Thomson limit. Our work is based on the formulation of the problem by Titarchuk & Lyubarskij (ApJ 450 (1995) 876), but extends their treatment by accommodating an arbitrary anisotropy of the source function. To achieve this, we expand the eigenfunctions of the integro/differential eigenvalue problem defining the spectral index of comptonised radiation in terms of Legendre polynomials and Chebyshev polynomials. The resulting algebraic eigenvalue problem is then solved by numerical means, yielding the spectral index and the full angular and spatial dependence of the specific intensity of radiation. For a thin (\tau_0 < 1) plasma disk, the radiation is strongly collimated along the disk surface - for an optical thickness of \tau_0 = 0.05, the radiation intensity along the surface is roughly ten times that along the direction of the normal, and varies only slightly with the electron temperature. Our results for the spectral index confirm those of Titarchuk & Lyubarskij over a wide range of electron temperature and optical depth; the largest difference we find is roughly 10% and occurs at low optical depth.

We present a model of the spectra of gamma-ray emitting blazars in which a single homogeneous emission region both emits synchrotron photons directly and scatters them to high (gamma-ray) energy before emission (a ``synchrotron self-Compton'' or SSC model). In contrast to previous work, we follow the full time dependent evolution of the electron and photon spectra, assuming a power-law form of the electron injection and examine the predictions of the model with regard to variability of the source. We apply these computations to the object Mkn 421, which displayed rapid variability in its X-ray and TeV emission during a multiwavelength campaign in 1994. This observation strongly implies that the same population of electrons produces the radiation in both energy bands. By fitting first the observed quiescent spectrum over all 18 orders of magnitude in frequency, we show that the time dependence of the keV/TeV flare could have been the result of a sudden increase in the maximum energy of the injected electrons. We show also that different types of flare may occur in this object and others, and that the energy band most sensitive to the properties of the acceleration mechanism is the X-ray band.

The theory of diffusive acceleration of energetic particles at shock fronts assumes charged particles undergo spatial diffusion in a uniform magnetic field. If, however, the magnetic field is not uniform, but has a stochastic or braided structure, the transport of charged particles across the average direction of the field is more complicated. Assuming quasi-linear behaviour of the field lines, the particles undergo sub-diffusion on short time scales. We derive the propagator for such motion, which differs from the Gaussian form relevant for diffusion, and apply it to a configuration with a plane shock front whose normal is perpendicular to the average field direction. Expressions are given for the acceleration time as a function of the diffusion coefficient of the wandering magnetic field lines and the spatial diffusion coefficient of the charged particles parallel to the local field. In addition we calculate the spatial dependence of the particle density in both the upstream and downstream plasmas. In contrast to the diffusive case, the density of particles at the shock front is lower than it is far downstream. This is a consequence of the partial trapping of particles by structures in the magnetic field. As a result, the spectrum of accelerated particles is a power-law in momentum which is steeper than in the diffusive case. For a phase- space density $f\propto p^{-s}$, we find $s=\sdiff[1+1/(2\rcomp)]$, where $\rcomp$ is the compression ratio of the shock front and $\sdiff$ is the standard result of diffusive acceleration: $\sdiff=3\rcomp/(\rcomp-1)$. A strong shock in a monatomic ideal gas yields a spectrum of $s=4.5$. In the case of electrons, this corresponds to a radio synchrotron spectral index of $\alpha=0.75$.

The limits imposed on diffusive shock acceleration by upstream ion-neutral Alfven wave damping, and by ionisation and Coulomb losses of low energy particles, are calculated. Analytic solutions are given for the steady upstream wave excitation problem with ion-neutral damping and the resulting escaping upstream flux calculated. The time dependent problem is discussed and numerical solutions presented. Finally the significance of these results for possible observational tests of shock acceleration in supernova remnants is discussed.

The variable flux of TeV gamma-rays detected from Mkn 421 and Mkn 501 requires the presence of high energy electrons, which could in principle produce large numbers of electron/positron pairs, leading to an electromagnetic cascade. We point out that this scenario can be avoided if electrons are accelerated to high energy rectilinearly, rather than being injected isotropically into a blob, as in most of the models of the GeV gamma-ray emission. By balancing linear acceleration by an electric field against inverse Compton losses in the radiation field of the accretion disk we calculate the emitted spectra and find the conditions which must be fulfilled in order to exclude the development of electromagnetic cascades during acceleration. Assuming these to be fulfilled, we show that the maximum possible photon energy is approximately $10M_8^{2/5}\,$TeV, where $M_8$ is the mass of the central black hole in units of $10^8\,\msolar$. In addition we compute the optical depth to absorption of TeV photons on a possible isotropic scattered component and on the observed nonthermal radiation (in the case of Mkn 421) and find that TeV photons can escape provided the nonthermal X-rays originate in a jet moving with a Lorentz factor of around 8.

We propose a model for radio supernovae (RSN) based on the synchrotron emission from relativistic electrons which are diffusively accelerated at the expanding supernova shock. This model was originally developed for application to the optically thin emission observed from SN1987A. Here we generalise it by including the effects of free-free absorption from both an external screen and from material internal to the source, and by relaxing the restriction to an azimuthal B-field. We find a good fit to the entire set of radio data for the best observed highly-luminous RSN -- SN1986J -- with a reduced \chi^2 of 3.85. Applying the new model to SN1988Z, another intrinsically bright RSN, also yields a good fit (\chi^2_{\rm red}\approx 2) but this is less significant, because of the limited data on this distant (z=0.02) source. These fits suggest that the shock expands at constant speed, that the magnetic field within the source decreases with time according to t^{-2}, and that the compression ratio of the shock front is close to the value expected of a strong shock in an ideal gas of adiabatic index 5/3 -- indicating a relatively low value of the cosmic ray pressure compared with SN1987A. In the case of SN1986J we derive an explosion date in August/September 1982, a magnetic field at the position of the shock 1000 days after explosion of 4nT and a spatial diffusion coefficient of the electrons of 4 10^{19} cm^2/s, four orders of magnitude greater than the Bohm value. In addition, we obtain the optical depths to external and internal absorption, and derive an estimate of the mass-loss rate.

The theory of first order Fermi acceleration at shocks assumes that particles diffuse due to scattering off slow-moving magnetic irregularities. However, cosmic rays are closely tied to magnetic field lines, and the transport process, particularly across the direction of the field, is likely to be more complicated. To describe cross field transport we employ recent extensions of the Rechester-Rosenbluth theory in localised stochastic regions of magnetic field. During acceleration at a shock and when the motion along the field is diffusive, there is a transition at a critical energy from \lq\lq sub-diffusive\rq\rq\ motion, where the mean square displacement of a particle increases with time as $t^{1/2}$, to compound diffusion, a combined process involving diffusion along a magnetic field which is itself wandering. Requiring this critical energy to be less than the system cut-off in a SNR of radius $R$ places an upper limit on the coherence length, $\lambda$, of the magnetic field for which diffusive shock acceleration can occur at quasi-perpendicular shocks when field line wandering dominates the effective transport of particles: $\lambda<R(U/v)^{1/2}$ where $U$ is the shock speed and $v$ the particle speed. The acceleration rate for particles in the sub-diffusive regime is derived and its implication for the origin of high energy cosmic rays is discussed.

The observations of radio emission from SNR~1987A can be accounted for on the basis of diffusive shock acceleration of electrons by the supernova blast wave. However, with this interpretation the observed spectral index implies that the compression ratio of the gas subshock is roughly $2.7$ rather than the value of $4$ expected of a strong shock front. We propose that in SNR~1987A, ions also undergo diffusive acceleration at the shock, a process that is likely to be rapid. Unlike the electron population, the accelerated ions can have an important effect on the gas dynamics. We calculate this coupled gas and energetic particle dynamics on the basis of the two-fluid model, in which the accelerated ions provide an additional component to the total pressure acting on the fluid. By accelerating and possibly heating the upstream plasma, the initially strong shock is modified and a weaker subshock with an upstream precursor results. The electrons behave as test particles. They are accelerated at the evolving subshock, escape downstream, and emit synchrotron radiation in the swept up magnetic field. Two models are considered for the surroundings of the progenitor: that of a freely expanding wind of number density $n\propto r^{-2}$, and that of a wind confined by a shell of denser material, creating a stagnation zone of roughly constant density beyond the standing shock which terminates the free wind. We model the observed radio light curves and the relatively steep spectrum of SNR~1987A using similar values for the ion acceleration parameters to those used in models of cosmic ray acceleration in older SNRs which can also contain high Mach number shocks, and find a good fit for the case in which the termination shock is located at about $2\times 10^{15}\,$m from the progenitor.

When applied to the blast wave formed by the explosion of a massive star as a supernova (SN), the theory of diffusive particle acceleration at shock fronts predicts a very high energy density in cosmic rays. Almost immediately after particles begin to be injected into the process, the cosmic ray pressure rises until comparable to the ram-pressure encountered by the shock front. Those supernovae which are observed in the radio band i.e., radio supernovae (RSNe), provide direct evidence of particle acceleration in the form of synchrotron emitting electrons. Furthermore, these objects are particularly interesting, since they are usually surrounded by a relatively dense confining medium. The acceleration of cosmic rays can then lead to the production of very high energy (VHE) gamma-rays which arise from collisions between energetic particles and target nuclei. We estimate the cosmic ray energy density assuming a fraction $\phi\lesim1$ of the energy available at the shock front is converted into cosmic rays. Combining this with the parameters describing the environment of the SN progenitor, as deduced from observations, and from more detailed modelling, we compute the flux at Earth $\tevflux$ of photons of energy above $1\,$TeV. For the relatively weak but nearby supernova SN1987A we predict $\tevflux=2\times10^{-13}\,\fluxunits$ before the shock front encounters the ring of dense matter seen by the Hubble Space Telescope. Subsequently, the flux is expected to rise further. The medium around SN1993J in M81 is thought to have a density

We analyse the propagation of particles in the narrow funnel of a thick accretion disk. It is assumed that: (1) the funnel walls emit black body radiation with temperature decreasing outwards; (2) the magnetic and electric fields are longitudinal in the funnel. Such a scenario is consistent with models in which a large potential drop is induced by a rotating massive black hole threaded by a magnetic field (Blandford & Znajek 1977, Macdonald & Thorne 1982). The interaction of relativistic protons with thermal photons from the funnel results in direct production of pairs and/or pions. We discuss the relative importance of these processes for different conditions in the funnel (temperature and electric field profiles). Injected pairs interact with thermal photons in the Thomson or Klein-Nishina regime. Under some conditions, a pair avalanche results, which we assume saturates in the production of stable bunches, containing almost equal numbers of electrons and positrons. As a consequence, highly collimated gamma-ray photons are produced. We obtain gamma-ray spectra from our model in order to test its applicability to the gamma-ray emitting AGNs recently reported by EGRET.