We investigate a time-delayed epidemic model for multi-strain diseases with temporary immunity. In the absence of cross-immunity between strains, dynamics of each individual strain exhibits emergence and anni- hilation of limit cycles due to a Hopf bifurcation of the endemic equilibrium, and a saddle-node bifurcation of limit cycles depending on the time delay associated with duration of temporary immunity. Effects of all-to-all and non-local coupling topologies are systematically investigated by means of numerical simulations, and they suggest that cross-immunity is able to induce a diverse range of complex dynamical behaviors and synchro- nization patterns, including discrete traveling waves, solitary states, and amplitude chimeras. Interestingly, chimera states are observed for narrower cross-immunity kernels, which can have profound implications for understanding the dynamics of multi-strain diseases.

We study the environments of low- and high- excitation radio galaxies (LERGs and HERGs respectively) in the redshift range $0.01 < z < 0.4$, using a sample of 399 radio galaxies and environmental measurements from the Galaxy And Mass Assembly (GAMA) survey. In our analysis we use the fifth nearest neighbour density ($\Sigma_{5}$) and the GAMA galaxy groups catalogue (G3Cv6) and construct control samples of galaxies matched in {\update stellar mass and colour} to the radio-detected sample. We find that LERGs and HERGs exist in different environments and that this difference is dependent on radio luminosity. High-luminosity LERGs ($L_{\rm NVSS} \gtrsim 10^{24}$ W Hz$^{-1}$) lie in much denser environments than a matched radio-quiet control sample (about three times as dense, as measured by $\Sigma_{5}$), and are more likely to be members of galaxy groups ($82^{+5}_{-7}$ percent of LERGs are in GAMA groups, compared to $58^{+3}_{-3}$ percent of the control sample). In contrast, the environments of the HERGs and lower luminosity LERGs are indistinguishable from that of a matched control sample. Our results imply that high-luminosity LERGs lie in more massive haloes than non-radio galaxies of similar stellar mass and colour, in agreement with earlier studies (Wake et al. 2008; Donoso et al. 2010). When we control for the preference of LERGs to be found in groups, both high- and low- luminosity LERGs are found in higher-mass haloes ($\sim 0.2$ dex; at least 97 percent significant) than the non-radio control sample.

Granger causality is well established within the neurosciences for inference of directed functional connectivity from neurophysiological data. These data usually consist of time series which subsample a continuous-time biophysiological process. While it is well-known that subsampling can lead to imputation of spurious causal connections where none exist, here we address the equally important issue of the effects of subsampling on the ability to reliably detect causal connections which do exist. Neurophysiological processes typically feature signal propagation delays on multiple time scales; accordingly, we base our analysis on a distributed-lag, continuous-time stochastic model, and consider Granger causality in continuous time at finite prediction horizons. Via exact analytical solutions, we identify relationships among sampling frequency, underlying causal time scales and detectability of causalities. Our analysis reveals complex interactions between the time scale(s) of neural signal propagation and sampling frequency: we demonstrate that Granger causality decays exponentially as the sample time interval increases beyond causal delay times, identify detectability "black spots" and "sweet spots", and show that subsampling may sometimes improve detectability. We also demonstrate that the invariance of Granger causality under causal, invertible filtering fails at finite prediction horizons. We discuss the implications of our results for inference of Granger causality at the neural level from various neurophysiological recording modes, and emphasise that sampling rates for causal analysis of neurophysiological time series should be informed by domain-specific time scales.

Leakage of diffuse polarized emission into Stokes I caused by the polarized primary beam of the instrument might mimic the spectral structure of the 21-cm signal coming from the epoch of reionization (EoR) making their separation difficult. Therefore, understanding polarimetric performance of the antenna is crucial for a successful detection of the EoR signal. Here, we have calculated the accuracy of the nominal model beam of LOFAR in predicting the leakage from Stokes I to Q, U by comparing them with the corresponding leakage of compact sources actually observed in the 3C295 field. We have found that the model beam has errors of less than or equal to 10% on the predicted levels of leakage of ~1% within the field of view, i. e. if the leakage is taken out perfectly using this model the leakage will reduce to $10^{-3}$ of the Stokes I flux. If similar levels of accuracy can be obtained in removing leakage from Stokes Q, U to I, we can say, based on the results of our previous paper, that the removal of this leakage using this beam model would ensure that the leakage is well below the expected EoR signal in almost the whole instrumental k-space of the cylindrical power spectrum. We have also shown here that direction dependent calibration can remove instrumentally polarized compact sources, given an unpolarized sky model, very close to the local noise level.

SNO+ is a large liquid scintillator-based experiment located 2km underground at SNOLAB, Sudbury, Canada. It reuses the Sudbury Neutrino Observatory detector, consisting of a 12m diameter acrylic vessel which will be filled with about 780 tonnes of ultra-pure liquid scintillator. Designed as a multipurpose neutrino experiment, the primary goal of SNO+ is a search for the neutrinoless double-beta decay (0$\nu\beta\beta$) of 130Te. In Phase I, the detector will be loaded with 0.3% natural tellurium, corresponding to nearly 800 kg of 130Te, with an expected effective Majorana neutrino mass sensitivity in the region of 55-133 meV, just above the inverted mass hierarchy. Recently, the possibility of deploying up to ten times more natural tellurium has been investigated, which would enable SNO+ to achieve sensitivity deep into the parameter space for the inverted neutrino mass hierarchy in the future. Additionally, SNO+ aims to measure reactor antineutrino oscillations, low-energy solar neutrinos, and geoneutrinos, to be sensitive to supernova neutrinos, and to search for exotic physics. A first phase with the detector filled with water will begin soon, with the scintillator phase expected to start after a few months of water data taking. The 0$\nu\beta\beta$ Phase I is foreseen for 2017.

Detection of the 21-cm signal coming from the epoch of reionization (EoR) is challenging especially because, even after removing the foregrounds, the residual Stokes $I$ maps contain leakage from polarized emission that can mimic the signal. Here, we discuss the instrumental polarization of LOFAR and present realistic simulations of the leakages between Stokes parameters. From the LOFAR observations of polarized emission in the 3C196 field, we have quantified the level of polarization leakage caused by the nominal model beam of LOFAR, and compared it with the EoR signal using power spectrum analysis. We found that at 134--166 MHz, within the central 4$^\circ$ of the field the $(Q,U)\rightarrow I$ leakage power is lower than the EoR signal at $k<0.3$ Mpc$^{-1}$. The leakage was found to be localized around a Faraday depth of 0, and the rms of the leakage as a fraction of the rms of the polarized emission was shown to vary between 0.2-0.3\%, both of which could be utilized in the removal of leakage. Moreover, we could define an `EoR window' in terms of the polarization leakage in the cylindrical power spectrum above the PSF-induced wedge and below $k_\parallel\sim 0.5$ Mpc$^{-1}$, and the window extended up to $k_\parallel\sim 1$ Mpc$^{-1}$ at all $k_\perp$ when 70\% of the leakage had been removed. These LOFAR results show that even a modest polarimetric calibration over a field of view of $\lesssim 4^\circ$ in the future arrays like SKA will ensure that the polarization leakage remains well below the expected EoR signal at the scales of 0.02-1 Mpc$^{-1}$.

When the cosmic star formation history peaks (z ~ 2), galaxies vigorously fed by cosmic reservoirs are gas dominated and contain massive star-forming clumps, thought to form by violent gravitational instabilities in highly turbulent gas-rich disks. However, a clump formation event has not been witnessed yet, and it is debated whether clumps survive energetic feedback from young stars, thus migrating inwards to form galaxy bulges. Here we report spatially resolved spectroscopy of a bright off-nuclear emission line region in a galaxy at z = 1.987. Although this region dominates the star formation in the galaxy disk, its stellar continuum remains undetected in deep imaging, revealing an extremely young (age < 10 Myr) massive clump, forming through the gravitational collapse of > 10$^9$ M$_{\odot}$ of gas. Gas consumption in this young clump is > 10 times faster than in the host galaxy, displaying high star formation efficiency during this phase, in agreement with our hydrodynamic simulations. The frequency of older clumps with similar masses coupled with our initial estimate of their formation rate (~ 2.5 Gyr$^{-1}$) supports long lifetimes (~ 500 Myr), favouring scenarios where clumps survive feedback and grow the bulges of present-day galaxies.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

This is the first in a series of papers in which we study the application of spectroastrometry in the context of gas kinematical studies aimed at measuring the mass of supermassive black holes. The spectroastrometrical method consists in measuring the photocenter of light emission in different wavelength or velocity channels. In particular we explore the potential of spectroastrometry of gas emission lines in galaxy nuclei to constrain the kinematics of rotating gas disks and to measure the mass of putative supermassive black holes. By means of detailed simulations and test cases, we show that the fundamental advantage of spectroastrometry is that it can provide information on the gravitational potential of a galaxy on scales significantly smaller (~ 1/10) than the limit imposed by the spatial resolution of the observations. We then describe a simple method to infer detailed kinematical informations from spectroastrometry in longslit spectra and to measure the mass of nuclear mass concentrations. Such method can be applied straightforwardly to integral field spectra, which do not have the complexities due to a partial spatial covering of the source in the case of longslit spectra.

We use the emission line measurements of 3CR radio sources with redshift < 0.3, to explore their spectroscopic properties. The 3CR sources show a bimodal distribution of Excitation Index, a new spectroscopic indicator that measures the relative intensity of low and high excitation lines. This unveils the presence of two main sub-populations of radio-loud AGN, High and Low Excitation Galaxies (HEG and LEG, respectively). All broad-line objects are HEG from the point of view of their narrow emission line ratios and all HEG are FRII radio-galaxies with log L(178) [erg/s] > 32.8. Conversely LEG cover the whole range of radio power encompassed by this 3CR subsample (30.7 < log L(178) < 35.4) and they are of both FRI and FRII type. The brightest LEG are all FRII. HEG and LEG obey to two (quasi) linear correlations between the optical line and extended radio luminosities, with HEG being brighter than LEG in the [OIII] line by a factor of ~10. HEG and LEG are offset also in a plane that compares the black hole mass and the ionizing nuclear luminosity. However, although HEG are associated with higher nuclear luminosities, we find LEG among the brightest radio sources of the sample and with a clear FRII morphology, indistinguishable from those seen in HEG. This suggests that LEG are not simply objects with a lower level of accretion. We speculate that the differences between LEG and HEG are related to a different mode of accretion: LEG are powered by hot gas, while HEG require the presence of cold accreting material. The high temperature of the accreting gas in LEG accounts for the lack of "cold" structures (i.e. molecular torus and Broad Line Region), for the reduced radiative output of the accretion disk, and for the lower gas excitation. [ABRIDGED]

An array of sixteen 1 cm^3 CdZnTe semiconductor detectors was operated at the Gran Sasso Underground Laboratory (LNGS) to further investigate the feasibility of double-beta decay searches with such devices. As one of the double-beta decay experiments with the highest granularity the 4 x 4 array accumulated an overall exposure of 18 kg days. The setup and performance of the array is described. Half-life limits for various double-beta decay modes of Cd, Zn and Te isotopes are obtained. No signal has been found, but several limits beyond 10^20 years have been performed. They are an order of magnitude better than those obtained with this technology before and comparable to most other experimental approaches for the isotopes under investigation. An improved limit for the beta^+/EC decay of Te 120 is given.

We present luminosity and surface brightness distributions of 40,111 galaxies with K-band photometry from the United Kingdom Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS), Data Release 3 and optical photometry from Data Release 5 of the Sloan Digital Sky Survey (SDSS). Various features and limitations of the new UKIDSS data are examined, such as a problem affecting Petrosian magnitudes of extended sources. Selection limits in K- and r-band magnitude, K-band surface brightness and K-band radius are included explicitly in the 1/Vmax estimate of the space density and luminosity function. The bivariate brightness distribution in K-band absolute magnitude and surface brightness is presented and found to display a clear luminosity--surface brightness correlation that flattens at high luminosity and broadens at low luminosity, consistent with similar analyses at optical wavelengths. Best fitting Schechter function parameters for the K-band luminosity function are found to be M*-5 log h=-23.19 +/- 0.04, alpha=-0.81 +/- 0.04 and phi*=(0.0166 +/- 0.0008)h^3 Mpc^{-3}, although the Schechter function provides a poor fit to the data at high and low luminosity, while the luminosity density in the K band is found to be j = (6.305 +/- 0.067) x 10^8 L_sun h Mpc^{-3}. However, we caution that there are various known sources of incompleteness and uncertainty in our results. Using mass-to-light ratios determined from the optical colours we estimate the stellar mass function, finding good agreement with previous results. Possible improvements are discussed that could be implemented when extending this analysis to the full LAS.

Bivariate luminosity functions (LFs) are computed for galaxies in the New York Value-Added Galaxy Catalogue, based on the Sloan Digital Sky Survey Data Release 4. The galaxy properties investigated are the morphological type, inverse concentration index, Sersic index, absolute effective surface brightness, reference frame colours, absolute radius, eClass spectral type, stellar mass and galaxy environment. The morphological sample is flux-limited to galaxies with r < 15.9 and consists of 37,047 classifications to an RMS accuracy of +/- half a class in the sequence E, S0, Sa, Sb, Sc, Sd, Im. These were assigned by an artificial neural network, based on a training set of 645 eyeball classifications. The other samples use r < 17.77 with a median redshift of z ~ 0.08, and a limiting redshift of z < 0.15 to minimize the effects of evolution. Other cuts, for example in axis ratio, are made to minimize biases. A wealth of detail is seen, with clear variations between the LFs according to absolute magnitude and the second parameter. They are consistent with an early type, bright, concentrated, red population and a late type, faint, less concentrated, blue, star forming population. This bimodality suggests two major underlying physical processes, which in agreement with previous authors we hypothesize to be merger and accretion, associated with the properties of bulges and discs respectively. The bivariate luminosity-surface brightness distribution is fit with the Choloniewski function (a Schechter function in absolute magnitude and Gaussian in surface brightness). The fit is found to be poor, as might be expected if there are two underlying processes.

A theoretical investigation for implementing a scheme of forced evaporative cooling in radio-frequency (rf) adiabatic potentials is presented. Supposing the atoms to be trapped by a rf field RF1, the cooling procedure is facilitated using a second rf source RF2. This second rf field produces a controlled coupling between the spin states dressed by RF1. The evaporation is then possible in a pulsed or continuous mode. In the pulsed case, atoms with a given energy are transferred into untrapped dressed states by abruptly switching off the interaction. In the continuous case, it is possible for energetic atoms to adiabatically follow the doubly-dressed states and escape out of the trap. Our results also show that when the frequencies of the fields RF1 and RF2 are separated by at least the Rabi frequency associated with RF1, additional evaporation zones appear which can make this process more efficient.

The Geneva group has reported two Saturn-mass planets orbiting HD 83443 (K0V) with periods of 2.98 and 29.8 d. The two planets have raised interest in their dynamics because of the possible 10:1 orbital resonance and the strong gravitational interactions. We report precise Doppler measurements of HD 83443 obtained with the Keck/HIRES and the AAT/UCLES spectrometers. These measurements strongly confirm the inner planet with period of 2.985 d, with orbital parameters in very good agreement with those of the Geneva group. However these Doppler measurements show no evidence of the outer planet, at thresholds of 1/4 (3 m/s) of the reported velocity amplitude of 13.8 m/s. Thus, the existence of the outer planet is in question. Indeed, the current Doppler measurements reveal no evidence of any second planet with periods less than a year.

It has been shown that Linear Indexed Grammars can be processed in polynomial time by exploiting constraints which make possible the extensive use of structure-sharing. This paper describes a formalism that is more powerful than Linear Indexed Grammar, but which can also be processed in polynomial time using similar techniques. The formalism, which we refer to as Partially Linear PATR manipulates feature structures rather than stacks.