We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$\alpha$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

Sgr A* represents a unique laboratory for the detailed study of accretion processes around a low-luminosity supermassive black hole (SMBH). Recent X-ray observations have allowed for spatially resolved modeling of the emission from the accretion flow around the SMBH, placing tight constraints on the flux and spectral shape of the accretion from the inner region with $r<10^3~r_g$, where $r_g$ is the gravitational radius of the black hole. We present here the first modeling of the multi-band spectral energy distribution (SED) of this inner region to better constrain the physical condition of the innermost accretion flow. Our modeling uses the Markov chain Monte Carlo (MCMC) method to fit the SED, accounting for the limitations on the accretion rate at the outer radius of $10^3~r_g$ from the earlier works and the domination of the accretion flow within $30~r_g$ to the sub-mm bump. It is found that the fitting results of the outflow index could be very different. If only the most luminous part of the SED, the sub-mm bump, is considered, the outflow index is about 0, while if low-frequency radio data and X-ray data are also included, the outflow index could be 0.37 or even higher. The great difference of the fitting results indicates that the outflow index should be variable along radius, with a strong outflow in the outer region and a weak outflow in the innermost region. Such weak outflow agrees with numerical simulations and is possible to explain the multi-band SED even better.

The next generation of large scale imaging surveys (such as those conducted with the Large Synoptic Survey Telescope and Euclid) will require accurate photometric redshifts in order to optimally extract cosmological information. Gaussian Processes for photometric redshift estimation (GPz) is a promising new method that has been proven to provide efficient, accurate photometric redshift estimations with reliable variance predictions. In this paper, we investigate a number of methods for improving the photometric redshift estimations obtained using GPz (but which are also applicable to others). We use spectroscopy from the Galaxy and Mass Assembly Data Release 2 with a limiting magnitude of r<19.4 along with corresponding Sloan Digital Sky Survey visible (ugriz) photometry and the UKIRT Infrared Deep Sky Survey Large Area Survey near-IR (YJHK) photometry. We evaluate the effects of adding near-IR magnitudes and angular size as features for the training, validation and testing of GPz and find that these improve the accuracy of the results by ~15-20 per cent. In addition, we explore a post-processing method of shifting the probability distributions of the estimated redshifts based on their Quantile-Quantile plots and find that it improves the bias by ~40 per cent. Finally, we investigate the effects of using more precise photometry obtained from the Hyper Suprime-Cam Subaru Strategic Program Data Release 1 and find that it produces significant improvements in accuracy, similar to the effect of including additional features.

We present ARC2 (Astrophysically Robust Correction 2), an open-source Python-based systematics-correction pipeline to correct for the Kepler prime mission long cadence light curves. The ARC2 pipeline identifies and corrects any isolated discontinuities in the light curves, then removes trends common to many light curves. These trends are modelled using the publicly available co-trending basis vectors, within an (approximate) Bayesian framework with `shrinkage' priors to minimise the risk of over-fitting and the injection of any additional noise into the corrected light curves, while keeping any astrophysical signals intact. We show that the ARC2 pipeline's performance matches that of the standard Kepler PDC-MAP data products using standard noise metrics, and demonstrate its ability to preserve astrophysical signals using injection tests with simulated stellar rotation and planetary transit signals. Although it is not identical, the ARC2 pipeline can thus be used as an open source alternative to PDC-MAP, whenever the ability to model the impact of the systematics removal process on other kinds of signal is important.

The interplay between supermassive black holes (SMBHs) and their environments is believed to command an essential role in galaxy evolution. The majority of these SMBHs are in the radiative inefficient accretion phase where this interplay remains elusive, but suggestively important, due to few observational constraints. To remedy this, we directly fit 2-D hydrodynamic simulations to Chandra observations of Sgr A* with Markov Chain Monte Carlo sampling, self-consistently modelling the 2-D inflow-outflow solution for the first time. We find the temperature and density at flow onset are consistent with the origin of the gas in the stellar winds of massive stars in the vicinity of Sgr A*. We place the first observational constraints on the angular momentum of the gas and estimate the centrifugal radius, r$_c$ $\approx$ 0.056 r$_b$ $\approx8\times10^{-3}$ pc, where r$_b$ is the Bondi radius. Less than 1\% of the inflowing gas accretes onto the SMBH, the remainder being ejected in a polar outflow. For the first time, we decouple the quiescent point-like emission from the spatially extended flow. We find this point-like emission, accounting for $\sim4$\% of the quiescent flux, is spectrally too steep to be explained by unresolved flares, nor Bremsstrahlung, but is likely a combination of a relatively steep synchrotron powerlaw and the high energy tail of inverse-Compton emission. With this self-consistent model of the accretion flow structure, we make predictions for the flow dynamics and discuss how future X-ray spectroscopic observations can further our understanding of the Sgr A* accretion flow.

In this paper, we present a framework for assessing the effect of non-stationary Gaussian noise and radio frequency interference (RFI) on the signal to noise ratio, the number of false positives detected per true positive and the sensitivity of standard pulsar search pipelines. The results highlight the necessity to develop algorithms that are able to identify and remove non-stationary variations from the data before RFI excision and searching is performed in order to limit false positive detections. The results also show that the spectrum whitening algorithms currently employed, severely affect the efficiency of pulsar search pipelines by reducing their sensitivity to long period pulsars.

The Robo-AO \textit{Kepler} Planetary Candidate Survey is observing every \textit{Kepler} planet candidate host star with laser adaptive optics imaging to search for blended nearby stars, which may be physically associated companions and/or responsible for transit false positives. We present in this paper the results of our search for stars nearby 1629 \textit{Kepler} planet candidate hosts. With survey sensitivity to objects as close as $\sim$0.15" and magnitude differences $\Delta$m$\le$6, we find 223 stars in the vicinity of 206 target KOIs; 209 of these nearby stars have not previously been imaged in high resolution. We measure an overall nearby-star probability for \textit{Kepler} planet candidates of 12.6\%$\pm$0.9\% at separations between 0.15" and 4.0". Particularly interesting KOI systems are discussed, including 23 stars with detected companions which host rocky, habitable zone candidates, and five new candidate planet-hosting quadruple star systems. We explore the broad correlations between planetary systems and stellar binarity using the combined dataset of Baranec et al. (2016) and this paper. Our previous 2$\sigma$ result of a low binary fraction of KOIs hosting close-in giant planets is less apparent in this larger dataset. We also find a significant correlation between binary fraction and KOI number, suggesting possible variation between early and late \textit{Kepler} data releases.

The next generation of cosmology experiments will be required to use photometric redshifts rather than spectroscopic redshifts. Obtaining accurate and well-characterized photometric redshift distributions is therefore critical for Euclid, the Large Synoptic Survey Telescope and the Square Kilometre Array. However, determining accurate variance predictions alongside single point estimates is crucial, as they can be used to optimize the sample of galaxies for the specific experiment (e.g. weak lensing, baryon acoustic oscillations, supernovae), trading off between completeness and reliability in the galaxy sample. The various sources of uncertainty in measurements of the photometry and redshifts put a lower bound on the accuracy that any model can hope to achieve. The intrinsic uncertainty associated with estimates is often non-uniform and input-dependent, commonly known in statistics as heteroscedastic noise. However, existing approaches are susceptible to outliers and do not take into account variance induced by non-uniform data density and in most cases require manual tuning of many parameters. In this paper, we present a Bayesian machine learning approach that jointly optimizes the model with respect to both the predictive mean and variance we refer to as Gaussian processes for photometric redshifts (GPz). The predictive variance of the model takes into account both the variance due to data density and photometric noise. Using the SDSS DR12 data, we show that our approach substantially outperforms other machine learning methods for photo-z estimation and their associated variance, such as TPZ and ANNz2. We provide a Matlab and Python implementations that are available to download at https://github.com/OxfordML/GPz .

Pulsars are known to display short-term variability. Recently, examples of longer-term emission variability have emerged that are often correlated with changes in the rotational properties of the pulsar. To further illuminate this relationship, we have developed techniques to identify emission and rotation variability in pulsar data, and determine correlations between the two. Individual observations may be too noisy to identify subtle changes in the pulse profile. We use Gaussian process (GP) regression to model noisy observations and produce a continuous map of pulse profile variability. Generally, multiple observing epochs are required to obtain the pulsar spin frequency derivative ($\dot{\nu}$). GP regression is, therefore, also used to obtain $\dot{\nu}$, under the hypothesis that pulsar timing noise is primarily caused by unmodelled changes in $\dot{\nu}$. Our techniques distinguish between two types of variability: changes in the total flux density versus changes in the pulse shape. We have applied these techniques to 168 pulsars observed by the Parkes radio telescope, and see that although variations in flux density are ubiquitous, substantial changes in the shape of the pulse profile are rare. We reproduce previously published results and present examples of profile shape changing in seven pulsars; in particular, a clear new example of correlated changes in profile shape and rotation is found in PSR~J1602$-$5100. In the shape changing pulsars, a more complex picture than the previously proposed two state model emerges. We conclude that our simple assumption that all timing noise can be interpreted as $\dot{\nu}$ variability is insufficient to explain our dataset.

Accurate photometric redshifts are a lynchpin for many future experiments to pin down the cosmological model and for studies of galaxy evolution. In this study, a novel sparse regression framework for photometric redshift estimation is presented. Simulated and real data from SDSS DR12 were used to train and test the proposed models. We show that approaches which include careful data preparation and model design offer a significant improvement in comparison with several competing machine learning algorithms. Standard implementations of most regression algorithms have as the objective the minimization of the sum of squared errors. For redshift inference, however, this induces a bias in the posterior mean of the output distribution, which can be problematic. In this paper we directly target minimizing $\Delta z = (z_\textrm{s} - z_\textrm{p})/(1+z_\textrm{s})$ and address the bias problem via a distribution-based weighting scheme, incorporated as part of the optimization objective. The results are compared with other machine learning algorithms in the field such as Artificial Neural Networks (ANN), Gaussian Processes (GPs) and sparse GPs. The proposed framework reaches a mean absolute $\Delta z = 0.0026(1+z_\textrm{s})$, over the redshift range of $0 \le z_\textrm{s} \le 2$ on the simulated data, and $\Delta z = 0.0178(1+z_\textrm{s})$ over the entire redshift range on the SDSS DR12 survey, outperforming the standard ANNz used in the literature. We also investigate how the relative size of the training set affects the photometric redshift accuracy. We find that a training set of \textgreater 30 per cent of total sample size, provides little additional constraint on the photometric redshifts, and note that our GP formalism strongly outperforms ANNz in the sparse data regime for the simulated data set.

We re-analyse the publicly available radial velocity (RV) measurements for Alpha Cen B, a star hosting an Earth-mass planet candidate, Alpha Cen Bb, with 3.24 day orbital period. We demonstrate that the 3.24 d signal observed in the Alpha Cen B data almost certainly arises from the window function (time sampling) of the original data. We show that when stellar activity signals are removed from the RV variations, other significant peaks in the power spectrum of the window function are coincidentally suppressed, leaving behind a spurious yet apparently-significant 'ghost' of a signal that was present in the window function's power spectrum to begin with. Even when fitting synthetic data with time sampling identical to the original data, but devoid of any genuine periodicities close to that of the planet candidate, the original model used to infer the presence of Alpha Cen Bb leads to identical conclusions: viz., the 3$\sigma$ detection of a half-a-metre-per-second signal with 3.236 day period. Our analysis underscores the difficulty of detecting weak planetary signals in RV data, and the importance of understanding in detail how every component of an RV data set, including its time sampling, influences final statistical inference.

We study the unresolved X-ray emission in three Local Group dwarf elliptical galaxies (NGC 147, NGC 185 and NGC 205) using XMM-Newton observations, which most likely originates from a collection of weak X-ray sources, mainly cataclysmic variables and coronally active binaries. Precise knowledge of this stellar X-ray emission is crucial not only for understanding the relevant stellar astrophysics but also for disentangling and quantifying the thermal emission from diffuse hot gas in nearby galaxies.We find that the integrated X-ray emissivities of the individual dwarf ellipticals agree well with that of the Solar vicinity, supporting an often assumed but untested view that the X-ray emissivity of old stellar populations is quasi-universal in normal galactic environments, in which dynamical effects on the formation and destruction of binary systems are not important. The average X-ray emissivity of the dwarf ellipticals, including M32 studied in the literature, is measured to be $L_{0.5-2\ \rm {keV}}/M_{\ast} = (6.0 \pm 0.5 \pm 1.8) \times 10^{27} \ \rm{erg \ s^{-1} \ M_\odot^{-1}}$. We also compare this value to the integrated X-ray emissivities of Galactic globular clusters and old open clusters and discuss the role of dynamical effects in these dense stellar systems.

To date, the radial velocity (RV) method has been one of the most productive techniques for detecting and confirming extrasolar planetary candidates. Unfortunately, stellar activity can induce RV variations which can drown out or even mimic planetary signals - and it is notoriously difficult to model and thus mitigate the effects of these activity-induced nuisance signals. This is expected to be a major obstacle to using next-generation spectrographs to detect lower mass planets, planets with longer periods, and planets around more active stars. Enter Gaussian processes (GPs) which, we note, have a number of attractive features that make them very well suited to disentangling stellar activity signals from planetary signals. We present here a GP framework we developed to model RV time series jointly with ancillary activity indicators (e.g. bisector velocity spans, line widths, chromospheric activity indices), allowing the activity component of RV time series to be constrained and disentangled from e.g. planetary components. We discuss the mathematical details of our GP framework, and present results illustrating its encouraging performance on both synthetic and real RV datasets, including the publicly-available Alpha Centauri B dataset.

Fast Radio Bursts (FRBs), are millisecond radio signals that exhibit dispersion larger than what the Galactic electron density can account for. We have conducted a 1446 hour survey for Fast Radio Bursts (FRBs) at 145~MHz, covering a total of 4193 sq. deg on the sky. We used the UK station of the LOFAR radio telescope -- the Rawlings Array -- , accompanied for a majority of the time by the LOFAR station at Nan\c{c}ay, observing the same fields at the same frequency. Our real-time search backend, ARTEMIS, utilizes graphics processing units to search for pulses with dispersion measures up to 320 cm$^{-3}$ pc. Previous derived FRB rates from surveys around 1.4~GHz, and favoured FRB interpretations, motivated this survey, despite all previous detections occurring at higher dispersion measures. We detected no new FRBs above a signal-to-noise threshold of 10, leading to the most stringent upper limit yet on the FRB event rate at these frequencies: 29 sky$^{-1}$ day$^{-1}$ for 5~ms-duration pulses above 62~Jy. The non-detection could be due to scatter-broadening, limitations on the volume and time searched, or the shape of FRB flux density spectra. Assuming the latter and that FRBs are standard candles, the non-detection is compatible with the published FRB sky rate, if their spectra follow a power law with frequency ($\propto \nu^{\alpha}$), with $\alpha\gtrsim+0.1$, demonstrating a marked difference from pulsar spectra. Our results suggest that surveys at higher frequencies, including the low frequency component of the Square Kilometre Array, will have better chances to detect, estimate rates and understand the origin and properties of FRBs.

The Cygnus Loop has been the focus of substantial debate concerning the contribution of charge exchange (CX) to supernova remnant (SNR) X-ray emission. We take advantage of a distinct feature of CX, enhanced K{\alpha} forbidden line emission, and employ the energy centroid of the OVII K{\alpha} triplet as a diagnostic. Based on X-ray spectra extracted from an extensive set of Suzaku observations, we measure the energy centroid shifts of the triplet on and off the shock rim of the remnant. We find that enhanced forbidden to resonance line emission exists throughout much of the rim and this enhancement azimuthally correlates with non-radiative H{\alpha} filaments, a tracer of strong neutral-plasma interaction in the optical. We also show that alternative mechanisms cannot explain the enhancement observed. These results demonstrate the need to model the CX contribution to the X-ray emission of SNRs, particularly for shocks propagating in a partially neutral medium. Such modeling may be critically important to the correct measurements of the ionization, thermal, and chemical properties of SNRs.

The recently approved NASA K2 mission has the potential to multiply by an order of magnitude the number of short-period transiting planets found by Kepler around bright and low-mass stars, and to revolutionise our understanding of stellar variability in open clusters. However, the data processing is made more challenging by the reduced pointing accuracy of the satellite, which has only two functioning reaction wheels. We present a new method to extract precise light curves from K2 data, combining list-driven, soft-edged aperture photometry with a star-by-star correction of systematic effects associated with the drift in the roll-angle of the satellite about its boresight. The systematics are modelled simultaneously with the stars' intrinsic variability using a semi-parametric Gaussian process model. We test this method on a week of data collected during an engineering test in January 2014, perform checks to verify that our method does not alter intrinsic variability signals, and compute the precision as a function of magnitude on long-cadence (30-min) and planetary transit (2.5-hour) timescales. In both cases, we reach photometric precisions close to the precision reached during the nominal Kepler mission for stars fainter than 12th magnitude, and between 40 and 80 parts per million for brighter stars. These results confirm the bright prospects for planet detection and characterisation, asteroseismology and stellar variability studies with K2. Finally, we perform a basic transit search on the light curves, detecting 2 bona fide transit-like events, 7 detached eclipsing binaries and 13 classical variables.

Debris disks and asteroid belts are expected to form around young pulsars due to fallback material from their original supernova explosions. Disk material may migrate inwards and interact with a pulsar's magnetosphere, causing changes in torque and emission. Long term monitoring of PSR J0738-4042 reveals both effects. The pulse shape changes multiple times between 1988 and 2012. The torque, inferred via the derivative of the rotational period, changes abruptly from September 2005. This change is accompanied by an emergent radio component that drifts with respect to the rest of the pulse. No known intrinsic pulsar processes can explain these timing and radio emission signatures. The data lead us to postulate that we are witnessing an encounter with an asteroid or in-falling debris from a disk.

We outline a proposal to use the Kepler spacecraft in two-wheel mode to monitor a handful of young associations and open clusters, for a few weeks each. Judging from the experience of similar projects using ground-based telescopes and the CoRoT spacecraft, this program would transform our understanding of early stellar evolution through the study of pulsations, rotation, activity, the detection and characterisation of eclipsing binaries, and the possible detection of transiting exoplanets. Importantly, Kepler's wide field-of-view would enable key spatially extended, nearby regions to be monitored in their entirety for the first time, and the proposed observations would exploit unique synergies with the GAIA ESO spectroscopic survey and, in the longer term, the GAIA mission itself. We also outline possible strategies for optimising the photometric performance of Kepler in two-wheel mode by modelling pixel sensitivity variations and other systematics.

Space-based transit search missions such as Kepler are collecting large numbers of stellar light curves of unprecedented photometric precision and time coverage. However, before this scientific goldmine can be exploited fully, the data must be cleaned of instrumental artefacts. We present a new method to correct common-mode systematics in large ensembles of very high precision light curves. It is based on a Bayesian linear basis model and uses shrinkage priors for robustness, variational inference for speed, and a de-noising step based on empirical mode decomposition to prevent the introduction of spurious noise into the corrected light curves. After demonstrating the performance of our method on a synthetic dataset, we apply it to the first month of Kepler data. We compare the results, which are publicly available, to the output of the Kepler pipeline's pre-search data conditioning, and show that the two generally give similar results, but the light curves corrected using our approach have lower scatter, on average, on both long and short timescales. We finish by discussing some limitations of our method and outlining some avenues for further development. The trend-corrected data produced by our approach are publicly available.

Classifier combination methods need to make best use of the outputs of multiple, imperfect classifiers to enable higher accuracy classifications. In many situations, such as when human decisions need to be combined, the base decisions can vary enormously in reliability. A Bayesian approach to such uncertain combination allows us to infer the differences in performance between individuals and to incorporate any available prior knowledge about their abilities when training data is sparse. In this paper we explore Bayesian classifier combination, using the computationally efficient framework of variational Bayesian inference. We apply the approach to real data from a large citizen science project, Galaxy Zoo Supernovae, and show that our method far outperforms other established approaches to imperfect decision combination. We go on to analyse the putative community structure of the decision makers, based on their inferred decision making strategies, and show that natural groupings are formed. Finally we present a dynamic Bayesian classifier combination approach and investigate the changes in base classifier performance over time.

We investigate the variability properties of main sequence stars in the first month of Kepler data, using a new astrophysically robust systematics correction, and find that 60% of stars are more variable then the active Sun. We define low and high variability samples, with a cut corresponding to twice the variability level of the active Sun, and compare the properties of the stars belonging to each sample. We show tentative evidence that the more active stars have lower proper motions and may be located closer to the galactic plane. We also investigate the frequency content of the variability, finding clear evidence for periodic or quasi-periodic behaviour in 16% of stars, and showing that there exist significant differences in the nature of variability between spectral types. Of the periodic objects, most A and F stars have short periods (< 2 days) and highly sinusoidal variability, suggestive of pulsations, whilst G, K and M stars tend to have longer periods (> 5 days, with a trend towards longer periods at later spectral types) and show a mixture of periodic and stochastic variability, indicative of activity. Finally, we use auto-regressive models to characterise the stochastic component of the variability, and show that its typical amplitude and time-scale both increase towards later spectral types, which we interpret as a corresponding increase in the characteristic size and life-time of active regions.

Within the efforts to bring frontline interactive astrophysics and astronomy to the classroom, the Hands on Universe (HOU) developed a set of exercises and platform using real data obtained by some of the most advanced ground and space observatories. The backbone of this endeavour is a new free software Web tool - Such a Lovely Software for Astronomy based on Image J (Salsa J). It is student-friendly and developed specifically for the HOU project and targets middle and high schools. It allows students to display, analyze, and explore professionally obtained astronomical images, while learning concepts on gravitational dynamics, kinematics, nuclear fusion, electromagnetism. The continuous evolving set of exercises and tutorials is being completed with real (professionally obtained) data to download and detailed tutorials. The flexibility of the Salsa J platform tool enables students and teachers to extend the exercises with their own observations. The software developed for the HOU program has been designed to be a multi-platform, multi-lingual experience for image manipulation and analysis in the classroom. Its design enables easy implementation of new facilities (extensions and plugins), minimal in-situ maintenance and flexibility for exercise plugin. Here, we describe some of the most advanced exercises about astrophysics in the classroom, addressing particular examples on gravitational dynamics, concepts currently introduced in most sciences curricula in middle and high schools.

Astrophysical radio transients are excellent probes of extreme physical processes originating from compact sources within our Galaxy and beyond. Radio frequency signals emitted from these objects provide a means to study the intervening medium through which they travel. Next generation radio telescopes are designed to explore the vast unexplored parameter space of high time resolution astronomy, but require High Performance Computing (HPC) solutions to process the enormous volumes of data that are produced by these telescopes. We have developed a combined software /hardware solution (code named ARTEMIS) for real-time searches for millisecond radio transients, which uses GPU technology to remove interstellar dispersion and detect millisecond radio bursts from astronomical sources in real-time. Here we present an introduction to ARTEMIS. We give a brief overview of the software pipeline, then focus specifically on the intricacies of performing incoherent de-dispersion. We present results from two brute-force algorithms. The first is a GPU based algorithm, designed to exploit the L1 cache of the NVIDIA Fermi GPU. Our second algorithm is CPU based and exploits the new AVX units in Intel Sandy Bridge CPUs.

Transmission spectroscopy, which consists of measuring the wavelength-dependent absorption of starlight by a planet's atmosphere during a transit, is a powerful probe of atmospheric composition. However, the expected signal is typically orders of magnitude smaller than instrumental systematics, and the results are crucially dependent on the treatment of the latter. In this paper, we propose a new method to infer transit parameters in the presence of systematic noise using Gaussian processes, a technique widely used in the machine learning community for Bayesian regression and classification problems. Our method makes use of auxiliary information about the state of the instrument, but does so in a non-parametric manner, without imposing a specific dependence of the systematics on the instrumental parameters, and naturally allows for the correlated nature of the noise. We give an example application of the method to archival NICMOS transmission spectroscopy of the hot Jupiter HD 189733, which goes some way towards reconciling the controversy surrounding this dataset in the literature. Finally, we provide an appendix giving a general introduction to Gaussian processes for regression, in order to encourage their application to a wider range of problems.

One of the tenets of the radio pulsar observational picture is that the integrated pulse profiles are constant with time. This assumption underpins much of the fantastic science made possible via pulsar timing. Over the past few years, however, this assumption has come under question with a number of pulsars showing pulse shape changes on a range of timescales. Here, we show the dramatic appearance of a bright component in the pulse profile of PSR J0738-4042 (B0736-40). The component arises on the leading edge of the profile. It was not present in 2004 but strongly present in 2006 and all observations thereafter. A subsequent search through the literature shows the additional component varies in flux density over timescales of decades. We show that the polarization properties of the transient component are consistent with the picture of competing orthogonal polarization modes. Faced with the general problem of identifying and characterising average profile changes, we outline and apply a statistical technique based on a Hidden Markov Model. The value of this technique is established through simulations, and is shown to work successfully in the case of low signal-to-noise profiles.

The nature of the dwarf galaxy population as a function of location in the cluster and within different environments is investigated. We have previously described the results of a search for low surface brightness objects in data drawn from an East-West strip of the Virgo cluster (Sabatini et al., 2003) and have compared this to a large area strip outside of the cluster (Roberts et al., 2004). In this talk I compare the East-West data (sampling sub-cluster A and outward) to new data along a North-South cluster strip that samples a different region (part of sub-cluster A, and the N,M clouds) and with data obtained for the Ursa Major cluster and fields around the spiral galaxy M101. The sample of dwarf galaxies in different environments is obtained from uniform datasets that reach central surface brightness values of ~26 B mag/arcsec^2 and an apparent B magnitude of 21 (M_B=-10 for a Virgo Cluster distance of 16 Mpc). We discuss and interpret our results on the properties and distribution of dwarf low surface brightness galaxies in the context of variuos physical processes that are thought to act on galaxies as they form and evolve.

Galaxy harassment has been proposed as a physical process that morphologically transforms low surface density disc galaxies into dwarf elliptical galaxies in clusters. It has been used to link the observed very different morphology of distant cluster galaxies (relatively more blue galaxies with 'disturbed' morphologies) with the relatively large numbers of dwarf elliptical galaxies found in nearby clusters. One prediction of the harassment model is that the remnant galaxies should lie on low surface brightness tidal streams or arcs. We demonstrate in this paper that we have an analysis method that is sensitive to the detection of arcs down to a surface brightness of 29 B mag/arcsec^2 and then use this method to search for arcs around 46 Virgo cluster dwarf elliptical galaxies. We find no evidence for tidal streams or arcs and consequently no evidence for galaxy harassment as a viable explanation for the relatively large numbers of dwarf galaxies found in the Virgo cluster.

Many observations indicate that dark matter dominates the extra-galactic Universe, yet no totally dark structure of galactic proportions has ever been convincingly identified. Previously we have suggested that VIRGOHI 21, a 21-cm source we found in the Virgo Cluster using Jodrell Bank, was a possible dark galaxy because of its broad line-width (~200 km/s) unaccompanied by any visible gravitational source to account for it. We have now imaged VIRGOHI 21 in the neutral-hydrogen line and find what could be a dark, edge-on, spinning disk with the mass and diameter of a typical spiral galaxy. Moreover, VIRGOHI 21 has unquestionably been involved in an interaction with NGC 4254, a luminous spiral with an odd one-armed morphology, but lacking the massive interactor normally linked with such a feature. Numerical models of NGC 4254 call for a close interaction ~10^8 years ago with a perturber of ~10^11 solar masses. This we take as additional evidence for the massive nature of VIRGOHI 21 as there does not appear to be any other viable candidate. We have also used the Hubble Space Telescope to search for stars associated with the HI and find none down to an I band surface brightness limit of 31.1 +/- 0.2 mag/sq. arcsec.

Dark Matter supposedly dominates the extragalactic, yet no totally dark structure of galactic proportions has ever been convincingly identified. Earlier (Minchin et al. 2005) we suggested that VIRGOHI 21, a 21-cm source we found in the Virgo Cluster at Jodrell Bank using single-dish observations (Davies et al. 2004), was probably such a dark galaxy because of its broad line-width (~ 200 km/s) unaccompanied by any visible gravitational source to account for it. Now we have managed to image VIRGOHI 21 in the neutral-hydrogen line, and indeed we find what appears to be a dark, edge-on, spinning disc with the mass and diameter of a typical spiral galaxy. Moreover the disc has unquestionably interacted with NGC 4254, a luminous spiral with an odd one-armed morphology, but lacking the massive interactor invariably responsible for such a feature. Published numerical models (Vollmer, Huchtmeier & van Driel 2005) of NGC 4254 call for a close interaction ~ 10^8 years ago with a perturber of 10^11 solar masses. This we take as completely independent evidence for the massive nature of VIRGOHI 21.

The Monte Carlo method is a powerful tool for performing radiative equilibrium calculations, even in complex geometries. The main drawback of the standard Monte Carlo radiative equilibrium methods is that they require iteration, which makes them numerically very demanding. Bjorkman & Wood recently proposed a frequency distribution adjustment scheme, which allows radiative equilibrium Monte Carlo calculations to be performed without iteration, by choosing the frequency of each re-emitted photon such that it corrects for the incorrect spectrum of the previously re-emitted photons. Although the method appears to yield correct results, we argue that its theoretical basis is not completely transparent, and that it is not completely clear whether this technique is an exact rigorous method, or whether it is just a good and convenient approximation. We critically study the general problem of how an already sampled distribution can be adjusted to a new distribution by adding data points sampled from an adjustment distribution. We show that this adjustment is not always possible, and that it depends on the shape of the original and desired distributions, as well as on the relative number of data points that can be added. Applying this theorem to radiative equilibrium Monte Carlo calculations, we provide a firm theoretical basis for the frequency distribution adjustment method of Bjorkman & Wood, and we demonstrate that this method provides the correct frequency distribution through the additional requirement of radiative equilibrium. We discuss the advantages and limitations of this approach, and show that it can easily be combined with the presence of additional heating sources and the concept of photon weighting. However, the method may fail if small dust grains are included... (abridged)

VIRGOHI21 is an HI source detected in the Virgo Cluster survey of Davies et al. (2004) which has a neutral hydrogen mass of 10^8 M_solar and a velocity width of Delta V_20 = 220 km/s. From the Tully-Fisher relation, a galaxy with this velocity width would be expected to be 12th magnitude or brighter; however deep CCD imaging has failed to turn up a counterpart down to a surface-brightness level of 27.5 B mag/sq. arcsec. The HI observations show that it is extended over at least 16 kpc which, if the system is bound, gives it a minimum dynamical mass of ~10^11 M_solar and a mass to light ratio of M_dyn/L_B > 500 M_solar/L_solar. If it is tidal debris then the putative parents have vanished; the remaining viable explanation is that VIRGOHI21 is a dark halo that does not contain the expected bright galaxy. This object was found because of the low column density limit of our survey, a limit much lower than that achieved by all-sky surveys such as HIPASS. Further such sensitive surveys might turn up a significant number of the dark matter halos predicted by Dark Matter models.

In order to investigate the nature of dwarf Low Surface Brightness (LSB) galaxies we have undertaken a deep B and I band CCD survey of a 14 sq degree strip in the Virgo Cluster and applied a Fourier convolution technique to explore its dwarf galaxy population down to a central surface brightness of 26 B mag/sq arcsec and a total absolute B mag of ~ -10. In this paper we carry out an analysis of their morphology, (B-I) colours and atomic hydrogen content. We compare these properties with those of dwarf galaxies in other environments to try and assess how the cluster environment has influenced their evolution. Field dwarfs are generally of a more irregular morphology, are bluer and contain relatively more gas. We assess the importance that various physical processes have on the evolution of cluster dwarf galaxies (ram pressure stripping, tidal interactions, supernova driven gas loss). We suggest that enhanced star formation triggered by tidal interactions is the major reason for the very different general properties of cluster dwarfs: they have undergone accelerated evolution.

We present R-band and H-alpha+[N II] narrow-band imaging of FCC032, FCC206 and FCCB729, three dwarf elliptical galaxies (dEs) in the Fornax Cluster. These dEs contain significant amounts of ionized gas. FCC032 harbours a large ionized gas complex, consisting of several individual clouds, a superbubble and a filament that extends away from the galaxy centre. The ionized gas structures observed in FCC032 bear a strong resemblance to those observed in more gas-rich and more fiercely star-forming dwarf galaxies. FCC206, a very low surface brightness dE, contains one faint extended emission region, and two compact clouds. In FCCB729, the only nucleated galaxy in this sample, one of the ionized gas clouds coincides with the stellar nucleus. We derive ionized gas masses of a few 100 to 1000 solar masses for these galaxies. This brings our sample of dEs with ionized gas with H-alpha+[N II] narrow-band imaging to five objects. The ionized gas morphologies in these galaxies range from pure nuclear emission peaks to extended emission complexes. This morphological diversity could also indicate a diversity in ionizing processes in dEs with a warm interstellar medium (active galactic nuclei, starbursts, post-AGB stars,...). Inside FCC206, four star clusters appear to be on the verge of merging to form a nucleus in this presently non-nucleated dE. Understanding the formation of nuclei in dEs could give us important clues to the formation of super-massive black holes.

The HIDEEP survey (Minchin et al. 2003) was done in an attempt to find objects having low inferred neutral hydrogen column densities, yet they found a distribution which was strongly peaked at 10^20.65 cm^-2. In an attempt to understand this distribution and similar survey results, we model HI profiles of gas discs and use simple simulations of objects having a wide range of HI properties in the presence of an ionizing background. We find that inferred column density (N_HI^o) values, which are found by averaging total HI masses over some disc area, do not vary strongly with central column density (N_max) for detectable objects, so that even a population having a wide range of N_max values will give rise to a strongly peaked distribution of N_HI^o values. We find that populations of objects, having a wide range of model parameters, give rise to inferred column density distributions around 10^20.6+/-0.3 cm^-2. However, populations of fairly massive objects having a wide range of central column densities work best in reproducing the HIDEEP data, and these populations are also consistent with observed Lyman limit absorber counts. It may be necessary to look two orders of magnitude fainter than HIDEEP limits to detect ionized objects having central column densities <10^20 cm^-2, but the inferred column densities of already detected objects might be lower if their radii could be estimated more accurately.

According to the Cold Dark Matter (CDM) hierarchical clustering theory of galaxy and large scale structure formation, there should be numerous low mass dark matter haloes present in the Universe today. If these haloes contain sufficient stars they should be detectable as low luminosity stellar systems or dwarf galaxies. We have previously described a new detection method for faint low surface brightness objects and shown that there are relatively large numbers of very faint dwarf galaxies in the nearby Virgo cluster. In this paper we present results from a similar survey carried out on the Millennium Galaxy strip which runs along the celestial equator and samples a very different galaxy environment. We show that the dwarf-to-giant galaxy number ratio along this strip ranges from 0.7:1 to, at most, 6:1, corresponding to a flat luminosity function (alpha approx -0.8 to -1.0). This is very different to our value of 20:1 for the Virgo cluster. There is no population of low surface brightness dwarf galaxies in the field that have gone undetected by the redshift surveys. This result is exactly opposite to what CDM models predict for the environmental dependence of the dark matter mass function which is that there are proportionally more small dark matter haloes in lower density environments.

We have carried out a fully sampled large area ($4^{\circ} \times 8^{\circ}$) 21cm \HI line survey of part of the Virgo cluster using the Jodrell Bank multi-beam instrument. The survey has a sensitivity some 3 times better than the standard HIJASS and HIPASS surveys. We detect 31 galaxies, 27 of which are well known cluster members. The four new detections have been confirmed in the HIPASS data and by follow up Jodrell Bank pointed observations. One object lies behind M86, but the other 3 have no obvious optical counter parts upon inspection of the digital sky survey fields. These 3 objects were mapped at Arecibo with a smaller \am{3}{6} HPBW and a 4 times better sensitivity than the Jodrell Bank data, which allow an improved determination of the dimensions and location of two of the objects, but surprisingly failed to detect the third. The two objects are resolved by the Arecibo beam giving them a size far larger than any optical images in the nearby field. To our mass limit of $5 \times 10^{7}$ $\frac{\Delta v}{50 km s^{-1}}$ $M_{\odot}$ and column density limit of $3 \times 10^{18}$ $\frac{\Delta v}{50 km s^{-1}}$ atoms cm$^{-2}$ these new detections represent only about 2% of the cluster atomic hydrogen mass. Our observations indicate that the \HI mass function of the cluster turns down at the low mass end making it very different to the field galaxy \HI mass function. This is quite different to the Virgo cluster optical luminosity function which is much steeper than that in the general field. Many of the sample galaxies are relatively gas poor compared to \HI selected samples of field galaxies, confirming the 'anaemic spirals' view of Virgo cluster late type galaxies.

The varying dwarf galaxy populations in different environments poses a problem for Cold Dark Matter (CDM) hierarchical clustering models. In this paper we present results from a survey conducted in different environments to search for low surface brightness (LSB) dwarf galaxies.

We present SKIRT (Stellar Kinematics Including Radiative Transfer), a new Monte Carlo radiative transfer code that allows the calculation of the observed stellar kinematics of a dusty galaxy. The code incorporates the effects of both absorption and scattering by interstellar dust grains, and calculates the Doppler shift of the emerging radiation exactly by taking into account the velocities of the emitting stars and the individual scattering dust grains. The code supports arbitrary distributions of dust through a cellular approach, whereby the integration through the dust is optimized by means of a novel efficient trilinear interpolation technique. We apply our modelling technique to calculate the observed kinematics of realistic models for dusty disc galaxies. We find that the effects of dust on the mean projected velocity and projected velocity dispersion are severe for edge-on galaxies. For galaxies which deviate more than a few degrees from exactly edge-on, the effects are already strongly reduced. As a consequence, dust attenuation cannot serve as a possible way to reconcile the discrepancy between the observed shallow slopes of the inner rotation curves of LSB galaxies and the predictions of CDM cosmological models. For face-on galaxies, the velocity dispersion increases with increasing dust mass due to scattering, but the effects are limited, even for extended dust distributions. Finally, we show that serious errors can be made when the individual velocities of the dust grains are neglected in the calculations.

Some Lyman continuum photons are likely to escape from most galaxies, and these can play an important role in ionizing gas around and between galaxies, including gas that gives rise to Lyman alpha absorption. Thus the gas surrounding galaxies and in the intergalactic medium will be exposed to varying amounts of ionizing radiation depending upon the distances, orientations, and luminosities of any nearby galaxies. The ionizing background can be recalculated at any point within a simulation by adding the flux from the galaxies to a uniform quasar contribution. Normal galaxies are found to almost always make some contribution to the ionizing background radiation at redshift zero, as seen by absorbers and at random points in space. Assuming that about 2 percent of ionizing photons escape from a galaxy like the Milky Way, we find that normal galaxies make a contribution of at least 30 to 40 percent of the assumed quasar background. Lyman alpha absorbers with a wide range of neutral column densities are found to be exposed to a wide range of ionization rates, although the distribution of photoionization rates for absorbers is found to be strongly peaked. On average, less highly ionized absorbers are found to arise farther from luminous galaxies, while local fluctuations in the ionization rate are seen around galaxies having a wide range of properties.

The widely varying dwarf galaxy counts in different environments provide a strong challenge to standard hierarchical clustering models. The luminosity function is not universal, but seems to be strongly dependent upon environment. In this paper we describe an automated procedure for detecting and measuring very low surface brightness (LSB) features in deep CCD data. We apply this procedure to large area CCD survey fields of the Virgo cluster. We show that there are many more faint ($-10 \geq M_{B} \geq -14$) low surface brightness galaxies than what would be predicted from extrapolation of the Virgo cluster catalogue luminosity function. Over our limited range of measurement the faint end slope of the luminosity function becomes $\alpha=- 1.6$. The luminosity function is flatter in the inner regions of the cluster than it is in the outer regions. Although these galaxies contribute a small fraction of the total stellar light of the cluster, they may contribute significantly to the mass in galaxies if they have large mass-to-light ratios similar to those recently measured for Local Group dwarf galaxies.