We estimate disc-averaged RV variations of the Sun over the last magnetic cycle, from the single Fe I line observed by SDO/HMI, using a physical model for rotationally modulated magnetic activity that was previously validated against HARPS-N solar observations. We estimate the disc-averaged, unsigned magnetic flux and show that a simple linear fit to it reduces the RMS of RV variations by 62%, i.e. a factor of 2.6. We additionally apply the FF' method, which predicts RV variations based on a star's photometric variations. At cycle maximum, we find that additional physical processes must be at play beyond suppression of convective blueshift and velocity imablances resulting from brightness inhomogeneities, in agreement with recent studies of solar RV variations. By modelling RV variations over the magnetic cycle using a linear fit to the unsigned magnetic flux, we recover injected planets at an orbital period of about 300 days with RV semi-amplitudes down to 0.3 m/s. To reach semi-amplitudes of 0.1 m/s, we will need to identify and model additional physical phenomena that are not well traced by the unsigned magnetic flux or FF'. The unsigned magnetic flux is an excellent proxy for rotationally modulated, activity-induced RV variations, and could become a key tool in confirming and characterising Earth analogs orbiting Sun-like stars. The present study motivates ongoing and future efforts to develop observation and analysis techniques to measure the unsigned magnetic flux at high precision in slowly rotating, relatively inactive stars like the Sun.

Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally-driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/non-rocky transition in period-radius space. Here we present the confirmation of TOI-1235 b ($P=3.44$ days, $r_p=1.738^{+0.087}_{-0.076}$ R$_{\oplus}$), a planet whose size and period are intermediate between the competing model predictions thus making the system an important test case for emergence models of the rocky/non-rocky transition around early M dwarfs ($R_s=0.630\pm 0.015$ R$_{\odot}$, $M_s=0.640\pm 0.016$ M$_{\odot}$). We confirm the TESS planet discovery using reconnaissance spectroscopy, ground-based photometry, high-resolution imaging, and a set of 38 precise radial-velocities from HARPS-N and HIRES. We measure a planet mass of $6.91^{+0.75}_{-0.85}$ M$_{\oplus}$, which implies an iron core mass fraction of $20^{+15}_{-12}$% in the absence of a gaseous envelope. The bulk composition of TOI-1235 b is therefore consistent with being Earth-like and we constrain a H/He envelope mass fraction to be $<0.5$% at 90% confidence. Our results are consistent with model predictions from thermally-driven atmospheric mass loss but not with gas-poor formation, suggesting that the former class of processes remain efficient at sculpting close-in planets around early M dwarfs. Our RV analysis also reveals a strong periodicity close to the first harmonic of the photometrically-determined stellar rotation period that we treat as stellar activity, despite other lines of evidence favoring a planetary origin ($P=21.8^{+0.9}_{-0.8}$ days, $m_p\sin{i}=13.0^{+3.8}_{-5.3}$ M$_{\oplus}$) that cannot be firmly ruled out by our data.

We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbital periods of $P_b=0.77$ days, $P_c=12.25$ days and sizes $r_{p,b}=1.33\pm 0.07$ R$_{\oplus}$, $r_{p,c}=2.30\pm 0.16$ R$_{\oplus}$, the two planets span the radius valley in period-radius space around low mass stars thus making the system a laboratory to test competing theories of the emergence of the radius valley in that stellar mass regime. By combining 63 precise radial-velocity measurements from HARPS and HARPS-N, we measure planet masses of $m_{p,b}=2.62^{+0.48}_{-0.46}$ M$_{\oplus}$ and $m_{p,c}=8.6^{+1.6}_{-1.3}$ M$_{\oplus}$, which indicates that LTT 3780b has a bulk composition consistent with being Earth-like, while LTT 3780c likely hosts an extended H/He envelope. We show that the recovered planetary masses are consistent with predictions from both photoevaporation and from core-powered mass loss models. The brightness and small size of LTT 3780, along with the measured planetary parameters, render LTT 3780b and c as accessible targets for atmospheric characterization of planets within the same planetary system and spanning the radius valley.

Exoplanets can evolve significantly between birth and maturity as their atmospheres, orbits, and structures are shaped by their environment. Young planets ($<$1 Gyr) offer the opportunity to probe these sculpting processes. However, most of the known young planets orbit prohibitively faint stars. We present the discovery of two planets transiting HD 63433 (TOI 1726, TIC 130181866), a young Sun-like ($M_*=0.99\pm0.03$) star. Through kinematics, lithium abundance, and rotation, we confirm that HD 63433 is a member of the Ursa Major moving group ($\tau=414\pm23$ Myr). Based on the TESS light curve and updated stellar parameters, the planet radii are $2.15\pm0.10R_\oplus$ and $2.67\pm0.12R_\oplus$, the orbital periods are 7.11 and 20.55 days, and the orbital eccentricities are lower than abut 0.2. Using HARPS-N velocities, we measure the Rossiter-McLaughlin signal of the inner planet, demonstrating the orbit is prograde. Since the host star is bright (V=6.9), both planets are amenable to transmission spectroscopy, radial velocity measurements of their masses, and more precise determination of the stellar obliquity. This system is therefore poised to play an important role in our understanding of planetary system evolution in the first billion years after formation.

Stellar activity is the major roadblock on the path to finding true Earth-analogue planets with the Doppler technique. Thus, identifying new indicators that better trace magnetic activity (i.e. faculae and spots) is crucial to aid in disentangling these signals from that of a planet's Doppler wobble. In this work, we investigate activity related features as seen in disk-integrated spectra from the HARPS-N solar telescope. We divide high-activity spectral echelle orders by low-activity master templates (as defined using both log R'HK and images from the Solar Dynamics Observatory, SDO), creating "relative spectra". With resolved images of the surface of the Sun (via SDO), the faculae and spot filling factors can be calculated, giving a measure of activity independent of, and in addition to, log R'HK. We find pseudo-emission (and pseudo-absorption) features in the relative spectra that are similar to those reported in our previous work on alpha Cen B. In alpha Cen B, the features are shown to correlate better to changes in faculae filling factor than spot filling factor. In this work we more confidently identify changes in faculae coverage of the visible hemisphere of the Sun as the source of features produced in the relative spectra. Finally, we produce trailed spectra to observe the RV component of the features, which show that the features move in a redward direction as one would expect when tracking active regions rotating on the surface of a star.

Efforts to detect low-mass exoplanets using stellar radial velocities (RVs) are currently limited by magnetic photospheric activity. Suppression of convective blueshift is the dominant magnetic contribution to RV variability in low-activity Sun-like stars. Due to convective plasma motions, the magnitude of RV contributions from the suppression of convective blueshift is related to the depth of formation of photospheric spectral lines of a given species used to compute the RV time series. Meunier et al. (2017), used this relation to demonstrate a method for spectroscopic extraction of the suppression of convective blueshift in order to isolate RV contributions, including planetary RVs, that contribute equally to the timeseries for each spectral line. Here, we extract disk-integrated solar RVs from observations over a 2.5 year time span made with the solar telescope integrated with the HARPS-N spectrograph at the Telescopio Nazionale Galileo (La Palma, Canary Islands, Spain). We apply the methods outlined by Meunier et al. (2017). We are not, however, able to isolate physically meaningful contributions of the suppression of convective blueshift from this solar dataset, potentially because our dataset is from solar minimum when the suppression of convective blueshift may not sufficiently dominate activity contributions to RVs. This result indicates that, for low-activity Sun-like stars, one must include additional RV contributions from activity sources not considered in the Meunier et al. (2017) model at different timescales as well as instrumental variation in order to reach the sub-meter per second RV sensitivity necessary to detect low-mass planets in orbit around Sun-like stars.

Ultra-short period (USP) planets are a class of exoplanets with periods shorter than one day. The origin of this sub-population of planets is still unclear, with different formation scenarios highly dependent on the composition of the USP planets. A better understanding of this class of exoplanets will, therefore, require an increase in the sample of such planets that have accurate and precise masses and radii, which also includes estimates of the level of irradiation and information about possible companions. Here we report a detailed characterization of a USP planet around the solar-type star HD 80653 $\equiv$ EP 251279430 using the K2 light curve and 108 precise radial velocities obtained with the HARPS-N spectrograph, installed on the Telescopio Nazionale Galileo. From the K2 C16 data, we found one super-Earth planet ($R_{b}=1.613\pm0.071 R_{\oplus}$) transiting the star on a short-period orbit ($P_{\rm b}=0.719573\pm0.000021$ d). From our radial velocity measurements, we constrained the mass of HD 80653 b to $M_{b}=5.60\pm0.43 M_{\oplus}$. We also detected a clear long-term trend in the radial velocity data. We derived the fundamental stellar parameters and determined a radius of $R_{\star}=1.22\pm0.01 R_{\odot}$ and mass of $M_{\star}=1.18\pm0.04 M_{\odot}$, suggesting that HD 80653, has an age of $2.7\pm1.2$ Gyr. The bulk density ($\rho_{b} = 7.4 \pm 1.1$ g cm$^{-3}$) of the planet is consistent with an Earth-like composition of rock and iron with no thick atmosphere. Our analysis of the K2 photometry also suggests hints of a shallow secondary eclipse with a depth of 8.1$\pm$3.7 ppm. Flux variations along the orbital phase are consistent with zero. The most important contribution might come from the day-side thermal emission from the surface of the planet at $T\sim3480$ K.

Although several thousands of exoplanets have now been detected and characterized, observational biases have led to a paucity of long-period, low-mass exoplanets with measured masses and a corresponding lag in our understanding of such planets. In this paper we report the mass estimation and characterization of the long-period exoplanet Kepler-538b. This planet orbits a Sun-like star (V = 11.27) with M_* = 0.892 +/- (0.051, 0.035) M_sun and R_* = 0.8717 +/- (0.0064, 0.0061) R_sun. Kepler-538b is a 2.215 +/- (0.040, 0.034) R_earth sub-Neptune with a period of P = 81.73778 +/- 0.00013 d. It is the only known planet in the system. We collected radial velocity (RV) observations with HIRES on Keck I and HARPS-N on the TNG. We characterized stellar activity by a Gaussian process with a quasi-periodic kernel applied to our RV and cross correlation function full width at half maximum (FWHM) observations. By simultaneously modeling Kepler photometry, RV, and FWHM observations, we found a semi-amplitude of K = 1.68 +/- (0.39, 0.38) m s^-1 and a planet mass of M_p = 10.6 +/- (2.5, 2.4) M_earth. Kepler-538b is the smallest planet beyond P = 50 d with an RV mass measurement. The planet likely consists of a significant fraction of ices (dominated by water ice), in addition to rocks/metals, and a small amount of gas. Sophisticated modeling techniques such as those used in this paper, combined with future spectrographs with ultra high-precision and stability will be vital for yielding more mass measurements in this poorly understood exoplanet regime. This in turn will improve our understanding of the relationship between planet composition and insolation flux and how the rocky to gaseous transition depends on planetary equilibrium temperature.

Twenty-four years after the discoveries of the first exoplanets, the radial-velocity (RV) method is still one of the most productive techniques to detect and confirm exoplanets. But stellar magnetic activity can induce RV variations large enough to make it difficult to disentangle planet signals from the stellar noise. In this context, HD41248 is an interesting planet-host candidate, with RV observations plagued by activity-induced signals. We report on ESPRESSO observations of HD41248 and analyse them together with previous observations from HARPS with the goal of evaluating the presence of orbiting planets. Using different noise models within a general Bayesian framework designed for planet detection in RV data, we test the significance of the various signals present in the HD41248 dataset. We use Gaussian processes as well as a first-order moving average component to try to correct for activity-induced signals. At the same time, we analyse photometry from the TESS mission, searching for transits and rotational modulation in the light curve. The number of significantly detected Keplerian signals depends on the noise model employed, which can range from 0 with the Gaussian process model to 3 with a white noise model. We find that the Gaussian process alone can explain the RV data while allowing for the stellar rotation period and active region evolution timescale to be constrained. The rotation period estimated from the RVs agrees with the value determined from the TESS light curve. Based on the data that is currently available, we conclude that the RV variations of HD41248 can be explained by stellar activity (using the Gaussian process model) in line with the evidence from activity indicators and the TESS photometry.

Transiting extrasolar planets are key objects in the study of the formation, migration, and evolution of planetary systems. In particular, the exploration of the atmospheres of giant planets, through transmission spectroscopy or direct imaging, has revealed a large diversity in their chemical composition and physical properties. Studying these giant planets allows one to test the global climate models that are used for the Earth and other solar system planets. However, these studies are mostly limited either to highly-irradiated transiting giant planets or directly-imaged giant planets at large separations. Here we report the physical characterisation of the planets in a bright multi-planetary system (HIP41378) in which the outer planet, HIP41378 f is a Saturn-sized planet (9.2 $\pm$ 0.1 R$_\oplus$) with an anomalously low density of 0.09 $\pm$ 0.02 g cm$^{-3}$ that is not yet understood. Its equilibrium temperature is about 300 K. Therefore, it represents a planet with a mild temperature, in between the hot Jupiters and the colder giant planets of the Solar System. It opens a new window for atmospheric characterisation of giant exoplanets with a moderate irradiation, with the next-generation space telescopes such as JWST and ARIEL as well as the extremely-large ground-based telescopes. HIP41378 f is thus an important laboratory to understand the effect of the irradiation on the physical properties and chemical composition of the atmosphere of planets.

We present confirmation of the planetary nature of PH-2b, as well as the first mass estimates for the two planets in the Kepler-103 system. PH-2b and Kepler-103c are both long-period and transiting, a sparsely-populated category of exoplanet. We use {\it Kepler} light-curve data to estimate a radius, and then use HARPS-N radial velocities to determine the semi-amplitude of the stellar reflex motion and, hence, the planet mass. For PH-2b we recover a 3.5-$\sigma$ mass estimate of $M_p = 109^{+30}_{-32}$ M$_\oplus$ and a radius of $R_p = 9.49\pm0.16$ R$_\oplus$. This means that PH-2b has a Saturn-like bulk density and is the only planet of this type with an orbital period $P > 200$ days that orbits a single star. We find that Kepler-103b has a mass of $M_{\text{p,b}} = 11.7^{+4.31}_{-4.72}$ M$_{\oplus}$ and Kepler-103c has a mass of $M_{\text{p,c}} = 58.5^{+11.2}_{-11.4}$ M$_{\oplus}$. These are 2.5$\sigma$ and 5$\sigma$ results, respectively. With radii of $R_{\text{p,b}} = 3.49^{+0.06}_{-0.05}$ R$_\oplus$, and $R_{\text{p,c}} = 5.45^{+0.18}_{-0.17}$ R$_\oplus$, these results suggest that Kepler-103b has a Neptune-like density, while Kepler-103c is one of the highest density planets with a period $P > 100$ days. By providing high-precision estimates for the masses of the long-period, intermediate-mass planets PH-2b and Kepler-103c, we increase the sample of long-period planets with known masses and radii, which will improve our understanding of the mass-radius relation across the full range of exoplanet masses and radii.

(Abridged) We perform a detailed study of the main optical activity indicators (Ca II H & K, Balmer lines, Na I D$_{\rm 1}$ D$_{\rm 2}$, and He I D$_{\rm 3}$) measured for the Sun using the data provided by the HARPS-N solar-telescope feed at the Telescopio Nazionale Galileo. The value of the solar rotation period is found in all the activity indicators, with the only exception being H$\delta$. The derived values vary from 26.29 days (H$\gamma$ line) to 31.23 days (He I). From an analysis of sliding periodograms we find that in most of the activity indicators the spectral power is split into several "bands" of periods around 26 and 30 days, that might be explained by the migration of active regions between the equator and a latitude of $\sim$ 30$^{\circ}$, spot evolution or a combination of both effects. In agreement with previous works a typical lifetime of active regions of $\sim$ ten rotation periods is inferred from the pooled variance diagrams. We find that H$\alpha$, H$\beta$, H$\gamma$, H$\epsilon$, and He I show a significant correlation with the S index. Significant correlations between the contrast, bisector span, and the heliocentric radial velocity with the activity indexes are also found. We show that the full width at half maximum, the bisector, and the disc-integrated magnetic field correlate with the radial velocity variations. The correlation of the S index and H$\alpha$ changes with time, increasing with larger sun spot numbers and solar irradiance. A similar tendency with the S index - radial velocity correlation is also present in the data. Our results are consistent with a scenario in which higher activity favours the correlation between the S index and the H$\alpha$ activity indicators and between the S index and radial velocity variations.

The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets' masses, obstructing detection of Earth twins. Since 2015 July we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-velocity measurements and cross-correlation function (CCF) parameters obtained in the first 3 years of observation, interpreting them in the context of spatially-resolved solar observations. We describe a Bayesian mixture-model approach to automated data-quality monitoring. We provide dynamical and daily differential-extinction corrections to place the radial velocities in the heliocentric reference frame, and the CCF shape parameters in the sidereal frame. We achieve a photon-noise limited radial-velocity precision better than 0.43 m s$^{-1}$ per 5-minute observation. The day-to-day precision is limited by zero-point calibration uncertainty with an RMS scatter of about 0.4 m s$^{-1}$. We find significant signals from granulation and solar activity. Within a day, granulation noise dominates, with an amplitude of about 0.4 m s$^{-1}$ and an autocorrelation half-life of 15 minutes. On longer timescales, activity dominates. Sunspot groups broaden the CCF as they cross the solar disc. Facular regions temporarily reduce the intrinsic asymmetry of the CCF. The radial-velocity increase that accompanies an active-region passage has a typical amplitude of 5 m s$^{-1}$ and is correlated with the line asymmetry, but leads it by 3 days. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity.

The number of exoplanet detections continues to grow following the development of better instruments and missions. Key steps for the understanding of these worlds comes from their characterization and its statistical studies. We explore the metallicity-period-mass diagram for known exoplanets by using an updated version of The Stellar parameters for stars With ExoplanETs CATalog (SWEET-Cat), a unique compilation of precise stellar parameters for planet-host stars provided for the exoplanet community. Here we focus on the planets with minimum mass below 30 M$_{\oplus}$ which seems to present a possible correlation in the metallicity-period-mass diagram where the mass of the planet increases with both metallicity and period. Our analysis suggests that the general observed correlation may be not fully explained by observational biases. Additional precise data will be fundamental to confirm or deny this possible correlation.

K2-291 (EPIC 247418783) is a solar-type star with a radius of R_star = 0.899 $\pm$ 0.034 R_sun and mass of M_star=0.934 $\pm$ 0.038 M_sun. From K2 C13 data, we found one super-Earth planet (R_p = 1.589+0.095-0.072 R_Earth) transiting this star on a short period orbit (P = 2.225177 +6.6e-5 -6.8e-5 days). We followed this system up with adaptive-optic imaging and spectroscopy to derive stellar parameters, search for stellar companions, and determine a planet mass. From our 75 radial velocity measurements using HIRES on Keck I and HARPS-N on Telescopio Nazionale Galileo, we constrained the mass of EPIC 247418783b to M_p = 6.49 $\pm$ 1.16 M_Earth. We found it necessary to model correlated stellar activity radial velocity signals with a Gaussian process in order to more accurately model the effect of stellar noise on our data; the addition of the Gaussian process also improved the precision of this mass measurement. With a bulk density of 8.84+2.50-2.03 g cm-3, the planet is consistent with an Earth-like rock/iron composition and no substantial gaseous envelope. Such an envelope, if it existed in the past, was likely eroded away by photo-evaporation during the first billion years of the star's lifetime.

State of the art radial-velocity (RV) exoplanet searches are currently limited by RV signals arising from stellar magnetic activity. We analyze solar observations acquired over a 3-year period during the decline of Carrington Cycle 24 to test models of RV variation of Sun-like stars. A purpose-built solar telescope at the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) provides disk-integrated solar spectra, from which we extract RVs and $\log{R'_{\rm HK}}$. The Solar Dynamics Observatory (SDO) provides disk-resolved images of magnetic activity. The Solar Radiation and Climate Experiment (SORCE) provides near-continuous solar photometry, analogous to a Kepler light curve. We verify that the SORCE photometry and HARPS-N $\log{R'_{\rm HK}}$ correlate strongly with the SDO-derived magnetic filling factor, while the HARPS-N RV variations do not. To explain this discrepancy, we test existing models of RV variations. We estimate the contributions of the suppression of convective blueshift and the rotational imbalance due to brightness inhomogeneities to the observed HARPS-N RVs. We investigate the time variation of these contributions over several rotation periods, and how these contributions depend on the area of active regions. We find that magnetic active regions smaller than $60 \ \rm Mm^2$ do not significantly suppress convective blueshift. Our area-dependent model reduces the amplitude of activity-induced RV variations by a factor of two. The present study highlights the need to identify a proxy that correlates specifically with large, bright magnetic regions on the surfaces of exoplanet-hosting stars.

K2-36 is a K dwarf orbited by two small ($R_{\rm b}=1.43\pm0.08$ $R_\oplus$ and $R_{\rm c}=3.2\pm0.3$ $R_\oplus$), close-in ($a_{\rm b}$=0.022 AU and $a_{\rm c}$=0.054 AU) transiting planets discovered by Kepler/K2. They are representatives of two families of small planets ($R_{\rm p}$<4 $R_\oplus$) recently emerged from the analysis of Kepler data, with likely a different structure, composition and evolutionary pathways. We revise the fundamental stellar parameters and the sizes of the planets, and provide the first measurement of their masses and bulk densities, which we use to infer their structure and composition. We observed K2-36 with the HARPS-N spectrograph over $\sim$3.5 years, collecting 81 useful radial velocity measurements. The star is active, with evidence for increasing levels of magnetic activity during the observing time span. The radial velocity scatter is $\sim$17 \ms due to the stellar activity contribution, which is much larger that the semi-amplitudes of the planetary signals. We tested different methods for mitigating the stellar activity contribution to the radial velocity time variations and measuring the planet masses with good precision. We found that K2-36 is likely a $\sim$1 Gyr old system, and by treating the stellar activity through a Gaussian process regression, we measured the planet masses $m_{\rm b}$=3.9$\pm$1.1 $M_\oplus$ and $m_{\rm c}$=7.8$\pm$2.3 $M_\oplus$. The derived planet bulk densities $\rho_{\rm b}$=7.2$^{+2.5}_{-2.1}$ $g/cm^{3}$ and $\rho_{\rm c}$=1.3$^{+0.7}_{-0.5}$ $g/cm^{3}$ point out that K2-36\,b has a rocky, Earth-like composition, and K2-36\,c is a low-density sub-Neptune. Composed of two planets with similar orbital separations but different densities, K2-36 represents an optimal laboratory for testing the role of the atmospheric escape in driving the evolution of close-in, low-mass planets after $\sim$1 Gyr from their formation.

Measures of exoplanet bulk densities indicate that small exoplanets with radius less than 3 Earth radii ($R_\oplus$) range from low-density sub-Neptunes containing volatile elements to higher density rocky planets with Earth-like or iron-rich (Mercury-like) compositions. Such astonishing diversity in observed small exoplanet compositions may be the product of different initial conditions of the planet-formation process and/or different evolutionary paths that altered the planetary properties after formation. Planet evolution may be especially affected by either photoevaporative mass loss induced by high stellar X-ray and extreme ultraviolet (XUV) flux or giant impacts. Although there is some evidence for the former, there are no unambiguous findings so far about the occurrence of giant impacts in an exoplanet system. Here, we characterize the two innermost planets of the compact and near-resonant system Kepler-107. We show that they have nearly identical radii (about $1.5-1.6~R_\oplus$), but the outer planet Kepler-107c is more than twice as dense (about $12.6~\rm g\,cm^{-3}$) as the innermost Kepler-107b (about $5.3~\rm g\,cm^{-3}$). In consequence, Kepler-107c must have a larger iron core fraction than Kepler-107b. This imbalance cannot be explained by the stellar XUV irradiation, which would conversely make the more-irradiated and less-massive planet Kepler-107b denser than Kepler-107c. Instead, the dissimilar densities are consistent with a giant impact event on Kepler-107c that would have stripped off part of its silicate mantle. This hypothesis is supported by theoretical predictions from collisional mantle stripping, which match the mass and radius of Kepler-107c.

We investigated the discrepancy between planetary mass determination using the transit timing variations (TTVs) and radial velocities (RVs), by analysing the multi-planet system Kepler-9. Despite being the first system characterised with TTVs, there are several discrepant solutions in the literature, with those reporting lower planetary densities being apparently in disagreement with high-precision RV observations. To resolve this, we gathered HARPS-N RVs at epochs that maximised the difference between the predicted RV curves from discrepant solutions in the literature. We also re-analysed the full Kepler data-set and performed a dynamical fit, within a Bayesian framework, using the newly derived central and duration times of the transits. We compared these results with the RV data and found that our solution better describes the RV observations, despite the masses of the planets being nearly half that presented in the discovery paper. We therefore confirm that the TTV method can provide mass determinations that agree with those determined using high-precision RVs. The low densities of the planets place them in the scarcely populated region of the super-Neptunes / inflated sub-Saturns in the mass-radius diagram.

Super-Earths belong to a class of planet not found in the Solar System, but which appear common in the Galaxy. Given that some super-Earths are rocky, while others retain substantial atmospheres, their study can provide clues as to the formation of both rocky planets and gaseous planets, and - in particular - they can help to constrain the role of photo-evaporation in sculpting the exoplanet population. GJ 9827 is a system already known to host 3 super-Earths with orbital periods of 1.2, 3.6 and 6.2 days. Here we use new HARPS-N radial velocity measurements, together with previously published radial velocities, to better constrain the properties of the GJ 9827 planets. Our analysis can't place a strong constraint on the mass of GJ 9827 c, but does indicate that GJ 9827 b is rocky with a composition that is probably similar to that of the Earth, while GJ 9827 d almost certainly retains a volatile envelope. Therefore, GJ 9827 hosts planets on either side of the radius gap that appears to divide super-Earths into pre-dominantly rocky ones that have radii below $\sim 1.5 R_\oplus$, and ones that still retain a substantial atmosphere and/or volatile components, and have radii above $\sim 2 R_\oplus$. That the less heavily irradiated of the 3 planets still retains an atmosphere, may indicate that photoevaporation has played a key role in the evolution of the planets in this system.

Context: Exoplanets have now been proven to be very common. The number of its detections continues to grow following the development of better instruments and missions. One key step for the understanding of these worlds is their characterization, which mostly depend on their host stars. Aims:We perform a significant update of the Stars With ExoplanETs CATalog (SWEET-Cat), a unique compilation of precise stellar parameters for planet-host stars provided for the exoplanet community. Methods: We made use of high-resolution spectra for planet-host stars, either observed by our team or found in several public archives. The new spectroscopic parameters were derived for the spectra following the same homogeneous process (ARES+MOOG). The host star parameters were then merged together with the planet properties listed in exoplanet.eu to perform simple data analysis. Results: We present new spectroscopic homogeneous parameters for 106 planet-host stars. Sixty-three planet hosts are also reviewed with new parameters. We also show that there is a good agreement between stellar parameters derived for the same star but using spectra obtained from different spectrographs. The planet-metallicity correlation is reviewed showing that the metallicity distribution of stars hosting low-mass planets (below 30 M$_{\oplus}$) is indistinguishable from that from the solar neighborhood sample in terms of metallicity distribution.

This paper reports on the validation and mass measurement of K2-263b, a sub-Neptune orbiting a quiet G9V star. Using K2 data from campaigns C5 and C16, we find this planet to have a period of $50.818947\pm 0.000094$ days and a radius of $2.41\pm0.12$ R$_{\oplus}$. We followed this system with HARPS-N to obtain 67 precise radial velocities. A combined fit of the transit and radial velocity data reveals that K2-263b has a mass of $14.8\pm3.1$ M$_{\oplus}$. Its bulk density ($5.7_{-1.4}^{+1.6}$ g cm$^{-3}$) implies that this planet has a significant envelope of water or other volatiles around a rocky core. EPIC211682544b likely formed in a similar way as the cores of the four giant planets in our own Solar System, but for some reason, did not accrete much gas. The planetary mass was confirmed by an independent Gaussian process-based fit to both the radial velocities and the spectroscopic activity indicators. K2-263b belongs to only a handful of confirmed K2 exoplanets with periods longer than 40 days. It is among the longest periods for a small planet with a precisely determined mass using radial velocities.

(shorter version)The aim of this work is to search for planets around intermediate-mass stars in open clusters by using RV data obtained with HARPS from an extensive survey with more than 15 years of observations for a sample of 142 giant stars in 17 open clusters. We present the discovery of a periodic RV signal compatible with the presence of a planet candidate in the 1.15 Gyr open cluster IC4651 orbiting the 2.06 M$_\odot$ star No. 9122. If confirmed, the planet candidate would have a minimum mass of 7.2 M$_{J}$ and a period of 747 days. However, we also find that the FWHM of the CCF varies with a period close to the RV, casting doubts on the planetary nature of the signal. We also provide refined parameters for the previously discovered planet around NGC2423 No. 3 but show evidence that the BIS of the CCF is correlated with the RV during some of the observing periods. This fact advises us that this might not be a real planet and that the RV variations could be caused by stellar activity and/or pulsations. Finally, we show that the previously reported signal by a brown dwarf around NGC4349 No. 127 is presumably produced by stellar activity modulation. The long-term monitoring of several red giants in open clusters has allowed us to find periodic RV variations in several stars. However, we also show that the follow-up of this kind of stars should last more than one orbital period to detect long-term signals of stellar origin. This work warns that although it is possible to detect planets around red giants, large-amplitude, long-period RV modulations do exist in such stars that can mimic the presence of an orbiting planetary body. Therefore, we need to better understand how such RV modulations behave as stars evolve along the RGB and perform a detailed study of all the possible stellar-induced signals (e.g. spots, pulsations, granulation) to comprehend the origin of RV variations.

We present the confirmation of a small, moderately-irradiated (F = 155 +/- 7 Fearth) Neptune with a substantial gas envelope in a P=11.8728787+/-0.0000085-day orbit about a quiet, Sun-like G0V star Kepler-1655. Based on our analysis of the Kepler light curve, we determined Kepler-1655b's radius to be 2.213+/-0.082 Rearth. We acquired 95 high-resolution spectra with TNG/HARPS-N, enabling us to characterize the host star and determine an accurate mass for Kepler-1655b of 5.0+3.1/-2.8 Mearth via Gaussian-process regression. Our mass determination excludes an Earth-like composition with 98\% confidence. Kepler-1655b falls on the upper edge of the evaporation valley, in the relatively sparsely occupied transition region between rocky and gas-rich planets. It is therefore part of a population of planets that we should actively seek to characterize further.

M-dwarf stars are promising targets for identifying and characterizing potentially habitable planets. K2-3 is a nearby (45 pc), early-type M dwarf hosting three small transiting planets, the outermost of which orbits close to the inner edge of the stellar (optimistic) habitable zone. The K2-3 system is well suited for follow-up characterization studies aimed at determining accurate masses and bulk densities of the three planets. Using a total of 329 radial velocity measurements collected over 2.5 years with the HARPS-N and HARPS spectrographs and a proper treatment of the stellar activity signal, we aim to improve measurements of the masses and bulk densities of the K2-3 planets. We use our results to investigate the physical structure of the planets. We analyse radial velocity time series extracted with two independent pipelines by using Gaussian process regression. We adopt a quasi-periodic kernel to model the stellar magnetic activity jointly with the planetary signals. We use Monte Carlo simulations to investigate the robustness of our mass measurements of K2-3\,c and K2-3\,d, and to explore how additional high-cadence radial velocity observations might improve them. Despite the stellar activity component being the strongest signal present in the radial velocity time series, we are able to derive masses for both planet b ($M_{\rm b}=6.6\pm1.1$ $M_{\rm \oplus}$) and planet c ($M_{\rm c}=3.1^{+1.3}_{-1.2}$ $M_{\rm \oplus}$). The Doppler signal due to K2-3\,d remains undetected, likely because of its low amplitude compared to the radial velocity signal induced by the stellar activity. The closeness of the orbital period of K2-3\,d to the stellar rotation period could also make the detection of the planetary signal complicated. [...]

Ultra-short period (USP) planets are a class of low mass planets with periods shorter than one day. Their origin is still unknown, with photo-evaporation of mini-Neptunes and in-situ formation being the most credited hypotheses. Formation scenarios differ radically in the predicted composition of USP planets, it is therefore extremely important to increase the still limited sample of USP planets with precise and accurate mass and density measurements. We report here the characterization of an USP planet with a period of 0.28 days around K2-141 (EPIC 246393474), and the validation of an outer planet with a period of 7.7 days in a grazing transit configuration. We derived the radii of the planets from the K2 light curve and used high-precision radial velocities gathered with the HARPS-N spectrograph for mass measurements. For K2-141b we thus inferred a radius of $1.51\pm0.05~R_\oplus$ and a mass of $5.08\pm0.41~M_\oplus$, consistent with a rocky composition and lack of a thick atmosphere. K2-141c is likely a Neptune-like planet, although due to the grazing transits and the non-detection in the RV dataset, we were not able to put a strong constraint on its density. We also report the detection of secondary eclipses and phase curve variations for K2-141b. The phase variation can be modeled either by a planet with a geometric albedo of $0.30 \pm 0.06$ in the Kepler bandpass, or by thermal emission from the surface of the planet at $\sim$3000K. Only follow-up observations at longer wavelengths will allow us to distinguish between these two scenarios.

We present precise radial velocity observations of WASP-47, a star known to host a hot Jupiter, a distant Jovian companion, and, uniquely, two additional transiting planets in short-period orbits: a super-Earth in a ~19 hour orbit, and a Neptune in a ~9 day orbit. We analyze our observations from the HARPS-N spectrograph along with previously published data to measure the most precise planet masses yet for this system. When combined with new stellar parameters and reanalyzed transit photometry, our mass measurements place strong constraints on the compositions of the two small planets. We find unlike most other ultra-short-period planets, the inner planet, WASP-47 e, has a mass (6.83 +/- 0.66 Me) and radius (1.810 +/- 0.027 Re) inconsistent with an Earth-like composition. Instead, WASP-47 e likely has a volatile-rich envelope surrounding an Earth-like core and mantle. We also perform a dynamical analysis to constrain the orbital inclination of WASP-47 c, the outer Jovian planet. This planet likely orbits close to the plane of the inner three planets, suggesting a quiet dynamical history for the system. Our dynamical constraints also imply that WASP-47 c is much more likely to transit than a geometric calculation would suggest. We calculate a transit probability for WASP-47 c of about 10%, more than an order of magnitude larger than the geometric transit probability of 0.6%.

HD 3167 is a bright (V = 8.9), nearby K0 star observed by the NASA K2 mission (EPIC 220383386), hosting two small, short-period transiting planets. Here we present the results of a multi-site, multi-instrument radial velocity campaign to characterize the HD 3167 system. The masses of the transiting planets are 5.02+/-0.38 MEarth for HD 3167 b, a hot super-Earth with a likely rocky composition (rho_b = 5.60+2.15-1.43 g/cm^3), and 9.80+1.30-1.24 MEarth for HD 3167 c, a warm sub-Neptune with a likely substantial volatile complement (rho_c = 1.97+0.94-0.59 g/cm^3). We explore the possibility of atmospheric composition analysis and determine that planet c is amenable to transmission spectroscopy measurements, and planet b is a potential thermal emission target. We detect a third, non-transiting planet, HD 3167 d, with a period of 8.509+/-0.045 d (between planets b and c) and a minimum mass of 6.90+/-0.71 MEarth. We are able to constrain the mutual inclination of planet d with planets b and c: we rule out mutual inclinations below 1.3 degrees as we do not observe transits of planet d. From 1.3-40 degrees, there are viewing geometries invoking special nodal configurations which result in planet d not transiting some fraction of the time. From 40-60 degrees, Kozai-Lidov oscillations increase the system's instability, but it can remain stable for up to 100Myr. Above 60 degrees, the system is unstable. HD 3167 promises to be a fruitful system for further study and a preview of the many exciting systems expected from the upcoming NASA TESS mission.

We report a detailed characterization of the Kepler-19 system. This star was previously known to host a transiting planet with a period of 9.29 days, a radius of 2.2 R$_\oplus$ and an upper limit on the mass of 20 M$_\oplus$. The presence of a second, non-transiting planet was inferred from the transit time variations (TTVs) of Kepler-19b, over 8 quarters of Kepler photometry, although neither mass nor period could be determined. By combining new TTVs measurements from all the Kepler quarters and 91 high-precision radial velocities obtained with the HARPS-N spectrograph, we measured through dynamical simulations a mass of $8.4 \pm 1.6$ M$_\oplus$ for Kepler-19b. From the same data, assuming system coplanarity, we determined an orbital period of 28.7 days and a mass of $13.1 \pm 2.7$ M$_\oplus$ for Kepler-19c and discovered a Neptune-like planet with a mass of $20.3 \pm 3.4$ M$_\oplus$ on a 63 days orbit. By comparing dynamical simulations with non-interacting Keplerian orbits, we concluded that neglecting interactions between planets may lead to systematic errors that could hamper the precision in the orbital parameters when the dataset spans several years. With a density of $4.32 \pm 0.87$ g cm$^{-3}$ ($0.78 \pm 0.16$ $\rho_\oplus$) Kepler-19b belongs to the group of planets with a rocky core and a significant fraction of volatiles, in opposition to low-density planets characterized by transit-time variations only and the increasing number of rocky planets with Earth-like density. Kepler-19 joins the small number of systems that reconcile transit timing variation and radial velocity measurements.

Context: Thanks to the importance that the star-planet relation has to our understanding of the planet formation process, the precise determination of stellar parameters for the ever increasing number of discovered extrasolar planets is of great relevance. Furthermore, precise stellar parameters are needed to fully characterize the planet properties. It is thus important to continue the efforts to determine, in the most uniform way possible, the parameters for stars with planets as new discoveries are announced. Aims: In this paper we present new precise atmospheric parameters for a sample of 50 stars with planets. The results are presented in the catalogue: SWEET-Cat. Methods: Stellar atmospheric parameters and masses for the 50 stars were derived assuming local thermodynamic equilibrium (LTE) and using high-resolution and high signal-to-noise spectra. The methodology used is based on the measurement of equivalent widths with ARES2 for a list of iron lines. The line abundances were derived using MOOG. We then used the curve of growth analysis to determine the parameters. We implemented a new minimization procedure which significantly improves the computational time. Results: The stellar parameters for the 50 stars are presented and compared with previously determined literature values. For SWEET-Cat, we compile values for the effective temperature, surface gravity, metallicity, and stellar mass for almost all the planet host stars listed in the Extrasolar Planets Encyclopedia. This data will be updated on a continuous basis. The data can be used for statistical studies of the star-planet correlation, and for the derivation of consistent properties for known planets.

Context. Distinguishing between a signal induced by either stellar activity or a planet is currently the main challenge in radial velocity searches for low-mass exoplanets. Even when the presence of a transiting planet and hence its period are known, stellar activity can be the main barrier to measuring the correct amplitude of the radial velocity signal. Several tools are being used to help understand which signals come from stellar activity in the data. Aims. We aim to present a new tool that can be used for the purpose of identifying periodicities caused by stellar activity, and show how it can be used to track the signal-to-noise ratio (SNR) of the detection over time. The tool is based on the principle that stellar activity signals are variable and incoherent. Methods. We calculate the Bayesian general Lomb-Scargle periodogram for subsets of data and by adding one extra data point we track what happens to the presence and significance of periodicities in the data. Publicly available datasets from HARPS and HARPS-N were used for this purpose. Additionally, we analysed a synthetic dataset that we created with SOAP2.0 to simulate pure stellar activity and a mixture of stellar activity and a planet. Results. We find that this tool can easily be used to identify unstable and incoherent signals, such as those introduced by stellar activity. The SNR of the detection grows approximately as the square root of the number of data points, in the case of a stable signal. This can then be used to make decisions on whether it is useful to keep observing a specific object. The tool is relatively fast and easy to use, and thus lends itself perfectly to a quick analysis of the data.

HD 179070, aka Kepler-21, is a V = 8.25 F6IV star and the brightest exoplanet host discovered by Kepler. An early detailed analysis by Howell et al. (2012) of the first thirteen months (Q0 - Q5) of Kepler light curves revealed transits of a planetary companion, Kepler-21b, with a radius of about 1.60 +/- 0.04 R_earth and an orbital period of about 2.7857 days. However, they could not determine the mass of the planet from the initial radial velocity observations with Keck-HIRES, and were only able to impose a 2-sigma upper limit of 10 M_earth. Here we present results from the analysis of 82 new radial velocity observations of this system obtained with HARPS-N, together with the existing 14 HIRES data points. We detect the Doppler signal of Kepler-21b with a radial velocity semi-amplitude K = 2.00 +/- 0.65 m/s, which corresponds to a planetary mass of 5.1 +/- 1.7 M_earth. We also measure an improved radius for the planet of 1.639 (+0.019, -0.015) R_earth, in agreement with the radius reported by Howell et al. (2012). We conclude that Kepler-21b, with a density of 6.4 +/- 2.1 g/cm^3, belongs to the population of terrestrial planets with iron, magnesium silicate interiors, which have lost the majority of their envelope volatiles via stellar winds or gravitational escape. The radial velocity analysis presented in this paper serves as example of the type of analysis that will be necessary to confirm the masses of TESS small planet candidates.

Kepler-20 is a solar-type star (V = 12.5) hosting a compact system of five transiting planets, all packed within the orbital distance of Mercury in our own Solar System. A transition from rocky to gaseous planets with a planetary transition radius of ~1.6 REarth has recently been proposed by several publications in the literature (Rogers 2015;Weiss & Marcy 2014). Kepler-20b (Rp ~ 1.9 REarth) has a size beyond this transition radius, however previous mass measurements were not sufficiently precise to allow definite conclusions to be drawn regarding its composition. We present new mass measurements of three of the planets in the Kepler-20 system facilitated by 104 radial velocity measurements from the HARPS-N spectrograph and 30 archival Keck/HIRES observations, as well as an updated photometric analysis of the Kepler data and an asteroseismic analysis of the host star (MStar = 0.948+-0.051 Msun and Rstar = 0.964+-0.018 Rsun). Kepler-20b is a 1.868+0.066-0.034 REarth planet in a 3.7 day period with a mass of 9.70+1.41-1.44 MEarth resulting in a mean density of 8.2+1.5-1.3 g/cc indicating a rocky composition with an iron to silicate ratio consistent with that of the Earth. This makes Kepler-20b the most massive planet with a rocky composition found to date. Furthermore, we report the discovery of an additional non-transiting planet with a minimum mass of 19.96+3.08-3.61 MEarth and an orbital period of ~34 days in the gap between Kepler-20f (P ~ 11 days) and Kepler-20d (P ~ 78 days).

We present the discovery of three new transiting giant planets, first detected with the WASP telescopes, and establish their planetary nature with follow up spectroscopy and ground-based photometric lightcurves. WASP-92 is an F7 star, with a moderately inflated planet orbiting with a period of 2.17 days, which has $R_p = 1.461 \pm 0.077 R_{\rm J}$ and $M_p = 0.805 \pm 0.068 M_{\rm J}$. WASP-93b orbits its F4 host star every 2.73 days and has $R_p = 1.597 \pm 0.077 R_{\rm J}$ and $M_p = 1.47 \pm 0.029 M_{\rm J}$. WASP-118b also has a hot host star (F6) and is moderately inflated, where $R_p = 1.440 \pm 0.036 R_{\rm J}$ and $M_p = 0.513 \pm 0.041 M_{\rm J}$ and the planet has an orbital period of 4.05 days. They are bright targets (V = 13.18, 10.97 and 11.07 respectively) ideal for further characterisation work, particularly WASP-118b, which is being observed by K2 as part of campaign 8. WASP-93b is expected to be tidally migrating outwards, which is divergent from the tidal behaviour of the majority of hot Jupiters discovered.

With about 2000 extrasolar planets confirmed, the results show that planetary systems have a whole range of unexpected properties. We present a full investigation of the HD219828 system, a bright metal-rich star for which a hot neptune has previously been detected. We used a set of HARPS, SOPHIE, and ELODIE radial velocities to search for the existence of orbiting companions to HD219828. A dynamical analysis is also performed to study the stability of the system and to constrain the orbital parameters and planet masses. We announce the discovery of a long period (P=13.1years) massive (msini=15.1MJup) companion (HD219828c) in a very eccentric orbit (e=0.81). The same data confirms the existence of a hot-neptune, HD219828b, with a minimum mass of 21 MEarth and a period of 3.83days. The dynamical analysis shows that the system is stable. The HD219828 system is extreme and unique in several aspects. First, among all known exoplanet systems it presents an unusually high mass ratio. We also show that systems like HD219828, with a hot neptune and a long-period massive companion are more frequent than similar systems with a hot jupiter instead. This suggests that the formation of hot neptunes follows a different path than the formation of their hot jovian counterparts. The high mass, long period, and eccentricity of HD219828c also make it a good target for Gaia astrometry as well as a potential target for atmospheric characterisation, using direct imaging or high-resolution spectroscopy. Astrometric observations will allow us to derive its real mass and orbital configuration. If a transit of HD219828b is detected, we will be able to fully characterise the system, including the relative orbital inclinations. With a clearly known mass, HD219828c may become a benchmark object for the range in between giant planets and brown dwarfs.

Stellar metallicity -- as a probe of the metallicity of proto-planetary disks -- is an important ingredient for giant planet formation, likely through its effect on the timescales in which rocky/icy planet cores can form. Giant planets have been found to be more frequent around metal-rich stars, in agreement with predictions based on the core-accretion theory. In the metal-poor regime, however, the frequency of planets, especially low-mass planets, and how it depends on metallicity are still largely unknown. As part of a planet search programme focused on metal-poor stars, we study the targets from this survey that were observed with HARPS on more than 75 nights. The main goals are to assess the presence of low-mass planets and provide a first estimate of the frequency of Neptunes and super-Earths around metal-poor stars. We perform a systematic search for planetary companions, both by analysing the periodograms of the radial-velocities and by comparing, in a statistically-meaningful way, models with an increasing number of Keplerians. A first constraint on the frequency of planets in our metal-poor sample is calculated considering the previous detection (in our sample) of a Neptune-sized planet around HD175607 and one candidate planet (with an orbital period of 68.42d and minimum mass $M_p \sin i = 11.14 \pm 2.47 M_{\oplus}$) for HD87838, announced in the present study. This frequency is determined to be close to 13% and is compared with results for solar-metallicity stars.

Context. The presence of a small-mass planet (M$_p<$0.1\,M$_{Jup}$) seems, to date, not to depend on metallicity, however, theoretical simulations have shown that stars with subsolar metallicities may be favoured for harbouring smaller planets. A large, dedicated survey of metal-poor stars with the HARPS spectrograph has thus been carried out to search for Neptunes and super-Earths. Aims. In this paper, we present the analysis of \object{HD175607}, an old G6 star with metallicity [Fe/H] = -0.62. We gathered 119 radial velocity measurements in 110 nights over a time span of more than nine years. Methods. The radial velocities were analysed using Lomb-Scargle periodograms, a genetic algorithm, a Markov chain Monte Carlo analysis, and a Gaussian processes analysis. The spectra were also used to derive stellar properties. Several activity indicators were analysed to study the effect of stellar activity on the radial velocities. Results. We find evidence for the presence of a small Neptune-mass planet (M$_{p}\sin i = 8.98\pm1.10$\,M$_{\oplus}$) orbiting this star with an orbital period $P = 29.01\pm0.02$\, days in a slightly eccentric orbit ($e=0.11\pm0.08$). The period of this Neptune is close to the estimated rotational period of the star. However, from a detailed analysis of the radial velocities together with the stellar activity, we conclude that the best explanation of the signal is indeed the presence of a planetary companion rather than stellar related. An additional longer period signal ($P\sim 1400$\,d) is present in the data, for which more measurements are needed to constrain its nature and its properties. Conclusions. HD\,175607 is the most metal-poor FGK dwarf with a detected low-mass planet amongst the currently known planet hosts. This discovery may thus have important consequences for planet formation and evolution theories.

We present here the detection of a system of four low-mass planets around the bright (V=5.5) and close-by (6.5 pc) star HD219134. This is the first result of the Rocky Planet Search program with HARPS-N on the TNG in La Palma. The inner planet orbits the star in 3.0937 +/-0.0004 days, on a quasi-circular orbit with a semi-major axis of 0.0382 +/- 0.0003 AU. Spitzer observations allowed us to detect the transit of the planet in front of the star making HD219134b the nearest known transiting planet to date. From the amplitude of the radial-velocity variation (2.33 +/- 0.24 m/s) and observed depth of the transit (359 +/- 38 ppm), the planet mass and radius are estimated to be 4.46 +/- 0.47 M_{\oplus} and 1.606 +/- 0.086 R_{\oplus} leading to a mean density of 5.89 +/- 1.17 g/cc, suggesting a rocky composition. One additional planet with minimum mass of 2.67 +/- 0.59 M_{\oplus} moves on a close-in, quasi-circular orbit with a period of 6.765 +/- 0.005 days. The third planet in the system has a period of 46.78 +/- 0.16 days and a minimum mass of 8.7 +/- 1.1 M{\oplus}, at 0.234 +/- 0.002 AU from the star. Its eccentricity is 0.32 +/- 0.14. The period of this planet is close to the rotational period of the star estimated from variations of activity indicators (42.3 +/- 0.1 days). The planetary origin of the signal is, however, the preferred solution as no indication of variation at the corresponding frequency is observed for activity-sensitive parameters. Finally, a fourth additional longer-period planet of mass of 62 +/- 6 M_{\oplus} orbits the star in 1190 days, on an eccentric orbit (e=0.27 +/- 0.11) at a distance of 2.14 +/- 0.27 AU.

We explore the possibility that the stellar relative abundances of different species can be used to constrain the bulk abundances of known transiting rocky planets. We use high resolution spectra to derive stellar parameters and chemical abundances for Fe, Si, Mg, O, and C in three stars hosting low mass, rocky planets: CoRoT-7, Kepler-10, and Kepler-93. These planets follow the same line along the mass-radius diagram, pointing toward a similar composition. The derived abundance ratios are compared with the solar values. With a simple stoichiometric model, we estimate the iron mass fraction in each planet, assuming stellar composition. We show that in all cases, the iron mass fraction inferred from the mass-radius relationship seems to be in good agreement with the iron abundance derived from the host star's photospheric composition. The results suggest that stellar abundances can be used to add constraints on the composition of orbiting rocky planets.

We performed a uniform and detailed abundance analysis of 12 refractory elements (Na, Mg, Al, Si, Ca, Ti, Cr, Ni, Co, Sc, Mn, and V) for a sample of 257 G- and K-type evolved stars from the CORALIE planet search program. To date, only one of these stars is known to harbor a planetary companion. We aimed to characterize this large sample of evolved stars in terms of chemical abundances and kinematics, thus setting a solid base for further analysis of planetary properties around giant stars. This sample, being homogeneously analyzed, can be used as a comparison sample for other planet-related studies, as well as for different type of studies related to stellar and Galaxy astrophysics. The abundances of the chemical elements were determined using an LTE abundance analysis relative to the Sun, with the spectral synthesis code MOOG and a grid of Kurucz ATLAS9 atmospheres. To separate the Galactic stellar populations both a purely kinematical approach and a chemical method were applied. We confirm the overabundance of Na in giant stars compared to the field FGK dwarfs. This enhancement might have a stellar evolutionary character, but departures from LTE may also produce a similar enhancement. Our chemical separation of stellar populations also suggests a "gap" in metallicity between the thick-disk and high-alpha metal-rich stars, as previously observed in dwarfs sample from HARPS. The present sample, as most of the giant star samples, also suffers from the B - V colour cut-off, which excludes low-log g stars with high metallicities, and high-logg star with low-[Fe/H]. For future studies of planet occurrence dependence on stellar metallicity around these evolved stars we suggest to use a sub-sample of stars in a "cut-rectangle" in the logg - [Fe/H] diagram to overcome the aforementioned issue.

Aims. In this work we derive new precise and homogeneous parameters for 37 stars with planets. For this purpose, we analyze high resolution spectra obtained by the NARVAL spectrograph for a sample composed of bright planet host stars in the northern hemisphere. The new parameters are included in the SWEET-Cat online catalogue. Methods. To ensure that the catalogue is homogeneous, we use our standard spectroscopic analysis procedure, ARES+MOOG, to derive effective temperatures, surface gravities, and metallicities. These spectroscopic stellar parameters are then used as input to compute the stellar mass and radius, which are fundamental for the derivation of the planetary mass and radius. Results. We show that the spectroscopic parameters, masses, and radii are generally in good agreement with the values available in online databases of exoplanets. There are some exceptions, especially for the evolved stars. These are analyzed in detail focusing on the effect of the stellar mass on the derived planetary mass. Conclusions. We conclude that the stellar mass estimations for giant stars should be managed with extreme caution when using them to compute the planetary masses. We report examples within this sample where the differences in planetary mass can be as high as 100% in the most extreme cases.

We present the evolution of dust and molecular gas properties in a sample of 11 $z\sim0.03$ starburst to post-starburst (PSB) galaxies selected to span an age sequence from ongoing starburst to 1 Gyr after the starburst ended. All PSBs harbour significant molecular gas and dust reservoirs and residual star formation, indicating that complete quenching of the starburst due to exhaustion or expulsion of gas has not occurred during this timespan. As the starburst ages, we observe a clear decrease in the star-formation efficiency, molecular gas and SFR surface density, and effective dust temperature, from levels coincident with starburst galaxies to those of normal star-forming galaxies. These trends are consistent with a natural decrease in the SFR following consumption of molecular gas by the starburst, and corresponding decrease in the interstellar radiation field strength as the starburst ages. The gas and dust contents of the PSBs are coincident with those of star-forming galaxies and molecular gas-rich early-type galaxies, and are not consistent with galaxies on the red-sequence. We find no evidence that the global gas reservoir is expelled by stellar winds or AGN feedback. Our results show that although a strong starburst in a low-redshift galaxy may cause the galaxy to ultimately have a lower specific SFR and be of an earlier morphological type, the galaxy will remain in the "green valley" for an extended time. Multiple such episodes may be needed to complete migration of the galaxy from the blue- to red-sequence.

We find that hot jupiter host stars within the T$_\mathrm{eff}$ range 5900-6300K show lower Li abundances, by 0.14 dex, than stars without detected planets. This offset has a significance at the level 7$\sigma$, pointing to a stronger effect of planet formation on Li abundances when the planets are more massive and migrate close to the star. However, we also find that the average v \textit{sin}i of (a fraction of) stars with hot jupiters is higher on average than for single stars in the same T$_\mathrm{eff}$ region, suggesting that rotationally-induced mixing (and not the presence of planets) might be the cause for a greater depletion of Li. We confirm that the mass-metallicity dependence of the Li dip is extended towards [Fe/H] $\sim$ 0.4 dex (beginning at [Fe/H] $\sim$ -0.4 dex for our stars) and that probably reflects the mass-metallicity correlation of stars of the same T$_\mathrm{eff}$ on the Main Sequence. We find that for the youngest stars ($<$ 1.5 Gyr) around the Li dip, the depletion of Li increases with v \textit{sin}i values, as proposed by rotationally-induced depletion models. This suggests that the Li dip consists of fast rotators at young ages whereas the most Li-depleted old stars show lower rotation rates (probably caused by the spin-down during their long lifes). We have also explored the Li evolution with [Fe/H] taking advantage of the metal-rich stars included in our sample. We find that Li abundance reaches its maximum around solar metallicity but decreases in the most metal-rich stars, as predicted by some models of Li Galactic production.

Context. Frequency analyses are very important in astronomy today, not least in the ever-growing field of exoplanets, where short-period signals in stellar radial velocity data are investigated. Periodograms are the main (and powerful) tools for this purpose. However, recovering the correct frequencies and assessing the probability of each frequency is not straightforward. Aims. We provide a formalism that is easy to implement in a code, to describe a Bayesian periodogram that includes weights and a constant offset in the data. The relative probability between peaks can be easily calculated with this formalism. We discuss the differences and agreements between the various periodogram formalisms with simulated examples. Methods. We used the Bayesian probability theory to describe the probability that a full sine function (including weights derived from the errors on the data values and a constant offset) with a specific frequency is present in the data. Results. From the expression for our Baysian generalised Lomb-Scargle periodogram (BGLS), we can easily recover the expression for the non-Bayesian version. In the simulated examples we show that this new formalism recovers the underlying periods better than previous versions. A Python-based code is available for the community.

Precise stellar parameters are crucial for several reasons, amongst which are the precise characterization of orbiting exoplanets and the correct determination of galactic chemical evolution. The atmospheric parameters are extremely important because all the other stellar parameters depend on them. Using our standard equivalent-width method on high-resolution spectroscopy, good precision can be obtained for the derived effective temperature and metallicity. The surface gravity, however, is usually not well constrained with spectroscopy. We use two different samples of FGK dwarfs to study the effect of the stellar surface gravity on the precise spectroscopic determination of the other atmospheric parameters. Furthermore, we present a straightforward formula for correcting the spectroscopic surface gravities derived by our method and with our linelists. Our spectroscopic analysis is based on Kurucz models in LTE, performed with the MOOG code to derive the atmospheric parameters. The surface gravity was either left free or fixed to a predetermined value. The latter is either obtained through a photometric transit light curve or derived using asteroseismology. We find first that, despite some minor trends, the effective temperatures and metallicities for FGK dwarfs derived with the described method and linelists are, in most cases, only affected within the errorbars by using different values for the surface gravity, even for very large differences in surface gravity, so they can be trusted. The temperatures derived with a fixed surface gravity continue to be compatible within 1 sigma with the accurate results of the InfraRed Flux Method (IRFM), as is the case for the unconstrained temperatures. Secondly, we find that the spectroscopic surface gravity can easily be corrected to a more accurate value using a linear function with the effective temperature.

Planetary studies demand precise and accurate stellar parameters as input to infer the planetary properties. Different methods often provide different results that could lead to biases in the planetary parameters. In this work, we present a refinement of the spectral synthesis technique designed to treat better more rapidly rotating FGK stars. This method is used to derive precise stellar parameters, namely effective temperature, surface gravity, metallicitity and rotational velocity. This procedure is tested for samples of low and moderate/fast rotating FGK stars. The spectroscopic analysis is based on the spectral synthesis package Spectroscopy Made Easy (SME), assuming Kurucz model atmospheres in LTE. The line list where the synthesis is conducted, is comprised of iron lines and the atomic data are derived after solar calibration. The comparison of our stellar parameters shows good agreement with literature values, both for low and for higher rotating stars. In addition, our results are on the same scale with the parameters derived from the iron ionization and excitation method presented in our previous works. We present new atmospheric parameters for 10 transiting planet-hosts as an update to the SWEET-Cat catalogue. We also re-analyse their transit light curves to derive new updated planetary properties.

The search for planets orbiting metal-poor stars is of uttermost importance for our understanding of the planet formation models. However, no dedicated searches have been conducted so far for very low mass planets orbiting such objects. Only a few cases of low mass planets orbiting metal-poor stars are thus known. Amongst these, HD41248 is a metal-poor, solar-type star on which a resonant pair of super-Earth like planets has In the present paper we present a new planet search program that is using the HARPS spectrograph to search for Neptunes and Super-Earths orbiting a sample of metal-poor FGK dwarfs. We then present a detailed analysis of an additional 162 radial velocity measurements of HD41248, obtained within this program, with the goal of confirming the existence of the proposed planetary system. We analyzed the precise radial velocities, obtained with the HARPS spectrograph, together with several stellar activity diagnostics and line profile indicators. A careful analysis shows no evidence for the planetary system previously announced. One of the signals, with a period of about 25 days, is shown to be related to the rotational period of the star, and is clearly seen in some of the activity proxies. The remaining signal (P~18 days) could not be convincingly retrieved in the new data set. We discuss possible causes for the complex (evolving) signals observed in the data of HD41248, proposing that they may be explained by the appearance and disappearance of active regions on the surface of a star with strong differential rotation, or by a combination of the sparse data sampling and active region evolution.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

We want to study the effects of the formation of planets and planetary systems on the atmospheric Li abundance of planet host stars. In this work we present new determinations of lithium abundances for 326 Main Sequence stars with and without planets in the T$_\mathrm{eff}$ range 5600-5900 K. 277 stars come from the HARPS sample, the remaining targets have been observed with a variety of high resolution spectrographs. We confirm significant differences in the Li distribution of solar twins (T$_\mathrm{eff}$ = T$_{\odot} \pm$ 80 K, log g = log g$_{\odot}$ $\pm$ 0.2 and [Fe/H] = [Fe/H]$_{\odot} \pm$ 0.2): the full sample of planet host stars (22) shows Li average values lower than "single" stars with no detected planets (60). If we focus in subsamples with narrower ranges in metallicity and age, we observe indications of a similar result though it is not so clear for some of the studied subsamples. Furthermore, we compare the observed spectra of several couples of stars with very similar parameters which show different Li abundances up to 1.6 dex. Therefore we show that neither age, nor mass nor metallicity of a parent star is the only responsible for enhanced Li depletion in solar analogues. We conclude that another variable must account for that difference and suggest that this could be the presence of planets which causes additional rotationally induced mixing in the external layers of planet host stars. Moreover, we find indications that the amount of depletion of Li in planet host solar-type stars is higher when the planets are more massive than Jupiter.

We explore the relations between physical and orbital properties of planets and properties of their host stars to identify the main observable signatures of the formation and evolution processes of planetary systems. We use a large sample of FGK dwarf planet hosts with stellar parameters derived in a homogeneous way from the SWEET-Cat database to study the relation between stellar metallicity and position of planets in the period-mass diagram. In the second part we use all the RV-detected planets orbiting FGK stars to explore the role of planet-disk and planet-planet interaction on the evolution of orbital properties of planets with masses above 1MJup. We show that planets orbiting metal-poor stars have longer periods than those in metal-rich systems. This trend is valid for masses at least from 10MEarth to 4MJup. Earth-like planets orbiting metal-rich stars always show shorter periods (fewer than 20 days) than those orbiting metal-poor stars. We also found statistically significant evidence that very high mass giants have on average more eccentric orbits than giant planets with lower mass.Finally, we show that the eccentricity of planets with masses higher than 4MJup tends to be lower for planets with shorter periods. Our results suggest that the planets in the P-MP diagram are evolving differently because of a mechanism that operates over a wide range of planetary masses. This mechanism is stronger or weaker depending on the metallicity of the respective system. One possibility is that planets in metal-poor disks form farther out from their central star and/or they form later and do not have time to migrate as far as the planets in metal-rich systems. The trends and dependencies obtained for very high mass planetary systems suggest that planet-disk interaction is a very important and orbit-shaping mechanism for planets in the high-mass domain. Shortened.

Context. Precise stellar parameters are crucial in exoplanet research for correctly determining of the planetary parameters. For stars hosting a transiting planet, determining of the planetary mass and radius depends on the stellar mass and radius, which in turn depend on the atmospheric stellar parameters. Different methods can provide different results, which leads to different planet characteristics.}%Spectroscopic surface gravities have shown to be poorly constrained, but the photometry of the transiting planet can provide an independent measurement of the surface gravity. Aims. In this paper, we use a uniform method to spectroscopically derive stellar atmospheric parameters, chemical abundances, stellar masses, and stellar radii for a sample of 90 transit hosts. Surface gravities are also derived photometrically using the stellar density as derived from the light curve. We study the effect of using these different surface gravities on the determination of the chemical abundances and the stellar mass and radius. Methods. A spectroscopic analysis based on Kurucz models in LTE was performed through the MOOG code to derive the atmospheric parameters and the chemical abundances. The photometric surface gravity was determined through isochrone fitting and the use of the stellar density, directly determined from the light curve. Stellar masses and radii are determined through calibration formulae. Results. Spectroscopic and photometric surface gravities differ, but this has very little effect on the precise determination of the stellar mass in our spectroscopic analysis. The stellar radius, and hence the planetary radius, is most affected by the surface gravity discrepancies. For the chemical abundances, the difference is, as expected, only noticable for the abundances derived from analyzing of lines of ionized species.

We present $^{12}$CO \textit{J=3$\rightarrow$2} maps of NGC 2976 and NGC 3351 obtained with the James Clerk Maxwell Telescope (JCMT), both early targets of the JCMT Nearby Galaxy Legacy Survey (NGLS). We combine the present observations with $^{12}$CO \textit{J=1$\rightarrow$0} data and find that the computed $^{12}$CO \textit{J=3$\rightarrow$2} to $^{12}$CO \textit{J=1$\rightarrow$0} line ratio ($R_{31}$) agrees with values measured in other NGLS field galaxies. We compute the M$_{\rm{H_2}}$ value and find that it is robust against the value of $R_{31}$ used. Using \textsc{Hi} data from the The \textsc{Hi} Nearby Galaxy Survey (THINGS) survey, we find a tight correlation between surface density of H$_2$ and star formation rate density for NGC 3351 when $^{12}$CO \textit{J=3$\rightarrow$2} data are used. Finally, we compare the $^{12}$CO \textit{J=3$\rightarrow$2} intensity with the PAH 8 $\mu$m surface brightness and find a good correlation in the high surface brightness regions. We extend this study to include all 25 \textit{Spitzer} Infrared Nearby Galaxies Survey (SINGS) galaxies within the NGLS sample and find a tight correlation at large spatial scales. We suggest that both PAH 8 $\mu$m and $^{12}$CO \textit{J=3$\rightarrow$2} are likely to originate in regions of active star formation.

It is still being debated whether the well-known metallicity - giant planet correlation for dwarf stars is also valid for giant stars. For this reason, having precise metallicities is very important. Different methods can provide different results that lead to discrepancies in the analysis of planet hosts. To study the impact of different analyses on the metallicity scale for evolved stars, we compare different iron line lists to use in the atmospheric parameter derivation of evolved stars. Therefore, we use a sample of 71 evolved stars with planets. With these new homogeneous parameters, we revisit the metallicity - giant planet connection for evolved stars. A spectroscopic analysis based on Kurucz models in local thermodynamic equilibrium (LTE) was performed through the MOOG code to derive the atmospheric parameters. Two different iron line list sets were used, one built for cool FGK stars in general, and the other for giant FGK stars. Masses were calculated through isochrone fitting, using the Padova models. Kolmogorov-Smirnov tests (K-S tests) were then performed on the metallicity distributions of various different samples of evolved stars and red giants. All parameters compare well using a line list set, designed specifically for cool and solar-like stars to provide more accurate temperatures. All parameters derived with this line list set are preferred and are thus adopted for future analysis. We find that evolved planet hosts are more metal-poor than dwarf stars with giant planets. However, a bias in giant stellar samples that are searched for planets is present. Because of a colour cut-off, metal-rich low-gravity stars are left out of the samples, making it hard to compare dwarf stars with giant stars. Furthermore, no metallicity enhancement is found for red giants with planets ($\log g < 3.0$\,dex) with respect to red giants without planets.

Due to the importance that the star-planet relation has to our understanding of the planet formation process, the precise determination of stellar parameters for the ever increasing number of discovered extra-solar planets is of great relevance. Furthermore, precise stellar parameters are needed to fully characterize the planet properties. It is thus important to continue the efforts to determine, in the most uniform way possible, the parameters for stars with planets as new discoveries are announced. In this paper we present new precise atmospheric parameters for a sample of 48 stars with planets. We then take the opportunity to present a new catalogue of stellar parameters for FGK and M stars with planets detected by radial velocity, transit, and astrometry programs. Stellar atmospheric parameters and masses for the 48 stars were derived assuming LTE and using high resolution and high signal-to-noise spectra. The methodology used is based on the measurement of equivalent widths for a list of iron lines and making use of iron ionization and excitation equilibrium principles. For the catalog, and whenever possible, we used parameters derived in previous works published by our team, using well defined methodologies for the derivation of stellar atmospheric parameters. This set of parameters amounts to over 65% of all planet host stars known, including more than 90% of all stars with planets discovered through radial velocity surveys. For the remaining targets, stellar parameters were collected from the literature.

The asteroseismic and planetary studies, like all research related to stars, need precise and accurate stellar atmospheric parameters as input. We aim at deriving the effective temperature (Teff), the surface gravity (log g), the metallicity ([Fe/H]), the projected rotational velocity (v sin i) and the MK type for 169 F, G, K, and M-type Kepler targets which were observed spectroscopically from the ground with five different instruments. We use two different spectroscopic methods to analyse 189 high-resolution, high-signal-to-noise spectra acquired for the 169 stars. For 67 stars, the spectroscopic atmospheric parameters are derived for the first time. KIC 9693187 and 11179629 are discovered to be double-lined spectroscopic binary systems. The results obtained for those stars for which independent determinations of the atmospheric parameters are available in the literature are used for a comparative analysis. As a result, we show that for solar-type stars the accuracy of present determinations of atmospheric parameters is +/- 150 K in Teff, +/- 0.15 dex in [Fe/H], and +/-? 0.3 dex in log g. Finally, we confirm that the curve-of-growth analysis and the method of spectral synthesis yield systematically different atmospheric parameters when they are applied to stars hotter than 6,000 K.

Temperature, surface gravity, and metallicitity are basic stellar atmospheric parameters necessary to characterize a star. We aim to improve the description of the spectroscopic temperatures especially for the cooler stars where the differences with the Infrared Flux Method are higher, as presented in previous work. Our spectroscopic analysis is based on the iron excitation and ionization balance, assuming Kurucz model atmospheres in LTE. The abundance analysis is determined using the code MOOG. We optimize the line list using a cool star with high resolution and high signal-to-noise spectrum, as a reference in order to check for weak, isolated lines. We test the quality of the new line list by re-deriving stellar parameters for 451 stars with high resolution and signal-to-noise HARPS spectra, that were analyzed in a previous work with a larger line list. The comparison in temperatures between this work and the latest IRFM shows that the differences for the cooler stars are significantly smaller and more homogeneously distributed than in previous studies for stars with temperatures below 5000 K. We use the new line list to re-derive parameters for some of the cooler stars that host planets. Finally, we present the impact of the new temperatures on the [Cr I/Cr II] and [Ti I/Ti II] abundance ratios that previously showed systematic trends with temperature.

It is generally accepted that the presence of a giant planet is strongly dependent on the stellar metallicity. A stellar mass dependence has also been investigated, but this dependence does not seem as strong as the metallicity dependence. Even for metallicity, however, the exact form of the correlation has not been established. In this paper, we test several scenarios for describing the frequency of giant planets as a function of its host parameters. We perform this test on two volume-limited samples (from CORALIE and HARPS). By using a Bayesian analysis, we quantitatively compared the different scenarios. We confirm that giant planet frequency is indeed a function of metallicity. However, there is no statistical difference between a constant or an exponential function for stars with subsolar metallicities contrary to what has been previously stated in the literature. The dependence on stellar mass could neither be confirmed nor be discarded.

In this work we investigate the problem concerning the presence of additional bodies gravitationally bounded with the WASP-3 system. We present eight new transits of this planet and analyse all the photometric and radial velocity data published so far. We did not observe significant periodicities in the Fourier spectrum of the observed minus calculated (O-C) transit timing and radial velocity diagrams (the highest peak having false-alarm probabilities of 56 per cent and 31 per cent, respectively) or long-term trends. Combining all the available information, we conclude that the radial velocity and transit timing techniques exclude, at 99 per cent confidence limit, any perturber more massive than M \gtrsim 100 M_Earth with periods up to 10 times the period of the inner planet. We also investigate the possible presence of an exomoon on this system and determined that considering the scatter of the O-C transit timing residuals a coplanar exomoon would likely produce detectable transits. This hypothesis is however apparently ruled out by observations conducted by other researchers. In case the orbit of the moon is not coplanar the accuracy of our transit timing and transit duration measurements prevents any significant statement. Interestingly, on the basis of our reanalysis of SOPHIE data we noted that WASP-3 passed from a less active (log R'_hk=-4.95) to a more active (log R'_hk=-4.8) state during the 3 yr monitoring period spanned by the observations. Despite no clear spot crossing has been reported for this system, this analysis claims for a more intensive monitoring of the activity level of this star in order to understand its impact on photometric and radial velocity measurements.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

The AzTEC 1.1 mm survey of the SCUBA HAlf Degree Extragalactic Survey (SHADES) fields is the largest (0.7 deg2) blank-field millimetre-wavelength survey undertaken to date at a resolution of ~18" and a depth of ~1 mJy. We have used the deep optical-to-radio multi-wavelength data in the SHADES Lockman Hole East and SXDF/UDS fields to obtain galaxy identifications for ~64% (~80% with tentative identifications) of the 148 AzTEC-SHADES 1.1 mm sources, exploiting deep radio and 24 um data complemented by methods based on 8 um flux-density and red optical-infrared (i-K) colour. This unusually high identification rate can be attributed to the relatively bright millimetre-wavelength flux-density threshold, combined with the relatively deep supporting multi-frequency data. We have further exploited the optical-mid-infrared-radio data to derive a ~60% (~75% with tentative identifications) complete redshift distribution for the AzTEC-SHADES sources, yielding a median redshift of z~2.2, with a high-redshift tail extending to at least z~4. Despite the larger area probed by the AzTEC survey relative to the original SCUBA SHADES imaging, the redshift distribution of the AzTEC sources is consistent with that displayed by the SCUBA sources, and reinforces tentative evidence that the redshift distribution of mm/sub-mm sources in the Lockman Hole field is significantly different from that found in the SXDF/UDS field. Comparison with simulated surveys of similar scale extracted from semi-analytic models based on the Millennium simulation indicates that this is as expected if the mm/sub-mm sources are massive (M>10^11 Mo) star-forming galaxies tracing large-scale structures over scales of 10-20 Mpc. This confirms the importance of surveys covering several square degrees (as now underway with SCUBA2) to obtain representative samples of bright (sub)mm-selected galaxies.

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

Context. The discovery of about 700 extrasolar planets, so far, has lead to the first statistics concerning extrasolar planets. The presence of giant planets seems to depend on stellar metallicity and mass. For example, they are more frequent around metal-rich stars,with an exponential increase in planet occurrence rates with metallicity. Aims. We analyzed two samples of metal-poor stars (-2.0 \leq [Fe/H] \leq 0.0) to see if giant planets are indeed rare around these objects. Radial velocity datasets were obtained with two different spectrographs (HARPS and HIRES). Detection limits for these data,expressed in minimum planetary mass and period, are calculated. These produce trustworthy numbers for the planet frequency. Methods. A general Lomb Scargle (GLS) periodogram analysis was used together with a bootstrapping method to produce the detection limits. Planet frequencies were calculated based on a binomial distribution function within metallicity bins. Results. Almost all hot Jupiters and most giant planets should have been found in these data. Hot Jupiters around metal-poor stars have a frequency lower than 1.0% at one sigma. Giant planets with periods up to 1800 days, however, have a higher frequency of $f_p = 2.63^{+2.5}_{-0.8}%$. Taking into account the different metallicities of the stars, we show that giant planets appear to be very frequent $(f_p = 4.48^{+4.04}_{-1.38}%)$ around stars with [Fe/H] > -0.7, while they are rare around stars with [Fe/H] \leq -0.7 (\leq 2.36% at one sigma). Conclusions. Giant planet frequency is indeed a strong function of metallicity, even in the low-metallicity tail. However, the frequencies are most likely higher than previously thought.

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.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

The results of an optical spectroscopic survey of a sample of young stellar objects (YSOs) and pre-main sequence (PMS) stars in the Lupus Clouds are presented. 92 objects were observed with VLT/FLAMES. All of those objects show IR excess as discovered by the Spitzer Legacy Program "From Molecular Cores to Planet-Forming Disks" (c2d). After reduction, 54 spectra with good signal-to-noise ratio are spectrally classified. Effective temperatures and luminosities are derived for these objects, and used to construct H-R diagrams for the population. The sample consists mostly of M-type stars, with 10% K-type stars. Individual ages and masses are inferred for the objects according to theoretical evolutionary models. The mean population age is found to be between 3.6 and 4.4 Myr, depending on the model, while the mean mass is found to be ~0.3 M for either model. Together with literature data, the distribution of spectral types is found to be similar to that in Chamaeleon I and IC348. The H{\alpha} line in emission, found in 49% of the sample, is used to distinguish between classical and weak-line T Tauri stars. 56% of the objects show H{\alpha} in emission and are accreting T Tauri stars. Mass accretion rates between 10-8 and 10-11 M yr-1 are determined from the full width at 10% of the H{\alpha} peak intensity. These mass accretion rates are, within a large scatter, consistent with the Mac \propto M 2 relation found in the literature.

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

An analysis of large-area CO J=3-2 maps from the James Clerk Maxwell Telescope for 12 nearby spiral galaxies reveals low velocity dispersions in the molecular component of the interstellar medium. The three lowest luminosity galaxies show a relatively flat velocity dispersion as a function of radius while the remaining nine galaxies show a central peak with a radial fall-off within 0.2-0.4 r(25). Correcting for the average contribution due to the internal velocitydispersions of a population of giant molecular clouds, the average cloud-cloud velocity dispersion across the galactic disks is 6.1 +/- 1.0 km/s (standard deviation 2.9 km/s), in reasonable agreement with previous measurements for the Galaxy andM33. The cloud-cloud velocity dispersion derived from the CO data is on average two times smaller than the HI velocity dispersion measured in the same galaxies. The low cloud-cloudvelocity dispersion implies that the molecular gas is the critical component determining the stability of the galactic disk against gravitational collapse, especially in those regions of the disk which are H2 dominated. The cloud-cloud velocity dispersion shows a significant positivecorrelation with both the far-infrared luminosity, which traces the star formation activity, and the K-band absolute magnitude, which traces the total stellar mass. For three galaxies in the Virgo cluster, smoothing the data to a resolution of 4.5 kpc (to match the typical resolution of high redshift CO observations) increases the measured velocity dispersion by roughly a factor of two, comparable to the dispersion measured recently in a normal galaxy at z=1. This comparison suggests that the mass and star formation rate surface densities may be similar in galaxies from z=0-1 and that the high star formation rates seen at z=1 may be partly due to the presence of physically larger molecular gas disks.

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

Spitzer spectroscopy has revealed that ~80% of submm galaxies (SMGs) are starburst (SB) dominated in the mid-infrared. Here we focus on the remaining ~20% that show signs of harboring powerful active galactic nuclei (AGN). We have obtained Spitzer-IRS spectroscopy of a sample of eight SMGs which are candidates for harboring powerful AGN on the basis of IRAC color-selection (S8/S4.5>2; i.e. likely power-law mid-infrared SEDs). SMGs with an AGN dominating (>50%) their mid-infrared emission could represent `missing link' sources in an evolutionary sequence involving a major merger. First of all, we detect PAH features in all of the SMGs, indicating redshifts from 2.5-3.4, demonstrating the power of the mid-infrared to determine redshifts for these optically faint dusty galaxies. Secondly, we see signs of both star-formation (from the PAH features) and AGN activity (from continuum emission) in our sample: 62% of the sample are AGN-dominated in the mid-infrared with a median AGN content of 56%, compared with <30% on average for typical SMGs, revealing that our IRAC color selection has successfully singled out sources with proportionately more AGN emission than typical SB-dominated SMGs. However, we find that only about 10% of these AGN dominate the bolometric emission of the SMG when the results are extrapolated to longer infrared wavelengths, implying that AGN are not a significant power source to the SMG population overall, even when there is evidence in the mid-infrared for substantial AGN activity. When existing samples of mid-infrared AGN-dominated SMGs are considered, we find that S8/S4.5>1.65 works well at selecting mid-infrared energetically dominant AGN in SMGs, implying a duty cycle of ~15% if all SMGs go through a subsequent mid-infrared AGN-dominated phase in the proposed evolutionary sequence.

We present a comparison of the SCUBA Half Degree Extragalactic Survey (SHADES) at 450, 850 and 1100 microns with deep guaranteed time 15 microns AKARI FU-HYU survey data and Spitzer guaranteed time data at 3.6-24 microns in the Lockman Hole East. The AKARI data was analysed using bespoke software based in part on the drizzling and minimum-variance matched filtering developed for SHADES, and was cross calibrated against Infrared Space Observatory (ISO) fluxes. Our stacking analyses find AKARI 15um galaxies with >~200 microJy contribute >10% of the 450 micron background, but only <4% of the 1100 micron background, suggesting that different populations contribute at mm-wavelengths. We confirm our earlier result that the ultra-deep 450 micron SCUBA-2 Cosmology Survey will be dominated by populations already detected by AKARI and Spitzer mid-infrared surveys. The superb mid-infrared wavelength coverage afforded by combining Spitzer and AKARI photometry is an excellent diagnostic of AGN contributions, and we find that (23-52)% of submm-selected galaxies have AGN bolometric fractions f_AGN>0.3.

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.

We have studied the uJy radio properties of K-selected galaxies detected in the Ultra-Deep Survey portion of UKIDSS using 610- and 1,400-MHz images from the VLA and GMRT. These deep radio mosaics, combined with the largest and deepest K image currently available, allow high-S/N detections of many K-selected sub-populations, including sBzK and pBzK star-forming and passive galaxies. We find a strong correlation between the radio and K fluxes and a linear relationship between SFR and K luminosity. We find no evidence, from either radio spectral indices or a comparison with submm-derived SFRs, that the full sample is strongly contaminated by AGN. The sBzK and pBzK galaxies have similar levels of radio flux, SFR and specific SFR (SSFR) at z < 1.4, suggesting there is strong contamination of the pBzK sample by star-forming galaxies. At z > 1.4, pBzK galaxies become difficult to detect in the radio stack, though the implied SFRs are still much higher than expected for passively evolving galaxies. Their radio emission may come from low-luminosity AGN. EROs straddle the passive and star-forming regions of the BzK diagram and also straddle the two groups in terms of their radio properties. K-bright ERO samples are dominated by passive galaxies and faint ERO samples contain more star-forming galaxies. The star-formation history (SFH) from stacking all K sources in the UDS agrees well with that derived for other wavebands and other radio surveys, at least out to z ~ 2. The radio-derived SFH then appears to fall more steeply than that measured at other wavelengths. The SSFR for K-selected sources rises strongly with redshift at all stellar masses, and shows a weak dependence on stellar mass. High- and low-mass galaxies show a similar decline in SSFR since z ~ 2 (abridged).

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near the reionization epoch at z=6.595. The giant LAE is found from the extensive 1 deg^2 Subaru narrow-band survey for z=6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 LAE candidates, this LAE is not only the brightest narrow-band object with L(Lya) = 3.9+/-0.2 x 10^43 erg/s in our survey volume of 10^6 Mpc^3, but also a spatially extended Lya nebula with the largest isophotal area whose major axis is at least ~3''. This object is more likely to be a large Lya nebula with a size of >~17 kpc than to be a strongly-lensed galaxy by a foreground object. Our Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is v(FWHM)=251+/-21 km/s, and that the line-center velocity changes by ~60 km/s in a 10-kpc range. The stellar mass and star-formation rate are estimated to be 0.9-5.0 x 10^10 Mo and >34 Mo/yr, respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this could be an important object for studying cooling clouds accreting onto a massive halo, or forming-massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment of inter-galactic medium.

We present large-area maps of the CO J=3-2 emission obtained at the James Clerk Maxwell Telescope for four spiral galaxies in the Virgo Cluster. We combine these data with published CO J=1-0, 24 micron, and Halpha images to measure the CO line ratios, molecular gas masses, and instantaneous gas depletion times. For three galaxies in our sample (NGC 4254, NGC4321, and NGC 4569), we obtain molecular gas masses of 7E8-3E9 Msun and disk-averaged instantaneous gas depletion times of 1.1-1.7 Gyr. We argue that the CO J=3-2 line is a better tracer of the dense star forming molecular gas than the CO J=1-0 line, as it shows a better correlation with the star formation rate surface density both within and between galaxies. NGC 4254 appears to have a larger star formation efficiency(smaller gas depletion time), perhaps because it is on its first passage through the Virgo Cluster. NGC 4569 shows a large-scale gradient in the gas properties traced by the CO J=3-2/J=1-0 line ratio, which suggests that its interaction with the intracluster medium is affecting the dense star-forming portion of the interstellar medium directly. The fourth galaxy in our sample, NGC 4579, has weak CO J=3-2 emission despite having bright 24 micron emission; however, much of the central luminosity in this galaxy may be due to the presence of a central AGN.

We present a comparison between the SCUBA Half Degree Extragalactic Survey (SHADES) at 450 and 850 microns in the Lockman Hole East with a deep Spitzer Space Telescope survey at 3.6-24 microns conducted in Guaranteed Time. Using stacking analyses we demonstrate a striking correspondence between the galaxies contributing the submm extragalactic background light, with those likely to dominate the backgrounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photometric redshifts, we show that at least a third of the Spitzer-identified submm galaxies at 1<z<1.5 appear to reside in overdensities when the density field is smoothed at 0.5-2 Mpc comoving diameters, supporting the high-redshift reversal of the local star formation - galaxy density relation. We derive the dust-shrouded cosmic star formation history of galaxies as a function of assembled stellar masses. For model stellar masses <10^11 Msun, this peaks at lower redshifts than the ostensible z~2.2 maximum for submm point sources, adding to the growing consensus for ``downsizing'' in star formation. Our surveys are also consistent with ``downsizing'' in mass assembly. Both the mean star formation rates <dM/dt> and specific star formation rates <(1/M)dM/dt> are in striking disagreement with some semi-analytic predictions from the Millenium simulation. The discrepancy could either be resolved with a top-heavy initial mass function, or a significant component of the submm flux heated by the interstellar radiation field.

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 microns has been carried out to map the regions around several 850 micron-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared luminosities and hence star-formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350 micron flux densities for 24 850 micron sources, seven of which are detected at >2.5-sigma within a 10 arcsec search radius of the 850 micron positions. When results from the literature are included the total number of 350 micron flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-infrared (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich (Mdust~9x10^8 Msun), luminous (Lfir~2x10^12 Lsun), star-forming galaxies with intrinsic dust temperatures of ~35 K and star-formation rates of ~400 Msun/yr. We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350 micron source counts, N350(>25mJy)~200-500 deg^-2.

We have studied the submillimetre (submm) properties of the following classes of near-infrared (NIR)-selected massive galaxies at high redshifts: BzK-selected star-forming galaxies (BzKs); distant red galaxies (DRGs); and extremely red objects (EROs). We used the SCUBA HAlf Degree Extragalactic Survey (SHADES), the largest uniform submm survey to date. Partial overlap of SIRIUS/NIR images and SHADES in SXDF has allowed us to identify 4 submm-bright NIR-selected galaxies, which are detected in the mid-infrared, 24 micron, and the radio, 1.4 GHz. We find that all of our submm-bright NIR-selected galaxies satisfy the BzK selection criteria, except for one galaxy whose B-z and z-K colours are however close to the BzK colour boundary. Two of the submm-bright NIR-selected galaxies satisfy all of the selection criteria we considered, i.e. they belong to the BzK-DRG-ERO overlapping population, or `extremely red' BzKs. Although these extremely red BzKs are rare (0.25 arcmin^{-2}), up to 20 % of this population could be submm galaxies. This fraction is significantly higher than that found for other galaxy populations studied here. Via a stacking analysis, we have detected the 850 micron flux of submm-faint BzKs and EROs in our SCUBA maps. While the contribution of z~2 BzKs to the submm background is about 10--15 % and similar to that from EROs typically at z~1, BzKs have a higher fraction (~30 %) of submm flux in resolved sources compared with EROs and submm sources as a whole. From the SED fitting analysis for both submm-bright and submm-faint BzKs, we found no clear signature that submm-bright BzKs are experiencing a specifically luminous evolutionary phase, compared with submm-faint BzKs. An alternative explanation might be that submm-bright BzKs are more massive than submm-faint ones.

We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.

Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

The modest significance of most sources detected in current submm surveys can potentially compromise some analyses due to the inclusion of spurious sources in catalogues typically selected at >3.0-3.5sigma. Here, we develop and apply a dual-survey extraction technique to SCUBA and MAMBO images of the Lockman Hole. Cut above 5sigma, our catalogue of SMGs is more robust than previous samples, with a reduced likelihood of real, but faint SMGs (beneath and around the confusion limit) entering via superposition with noise. The effective flux limit of the survey is well matched to our deep 1.4-GHz image. The former is sensitive to luminous, dusty galaxies at extreme redshifts whilst the latter probes the z<3 regime. A high fraction of our robust SMGs (~80%) have radio counterparts which, given the ~10% contamination by spurious sources, suggests that very distant SMGs (z>>3) are unlikely to make up more than ~10% of the bright SMG population. This implies that almost all of the bright SMG population is amenable to study via the deepest current radio imaging. We use these radio counterparts to provide an empirical calibration of the positional uncertainty in SMG catalogues. We then go on to outline the acquisition of redshifts for radio-identified SMGs, from sample selection in the submm, to counterpart selection in the radio and optical/IR, to slit placement on spectrograph masks. We determine a median of z=2.05 +/- 0.41 from a sample of six secure redshifts for unambigious radio-identified submm sources and z=2.14 +/- 0.27 when we include submm sources with multiple radio counterparts and/or less reliable redshifts. These figures are consistent with previous estimates, suggesting that our knowledge of the median redshift of bright SMGs population has not been biased by the low significance of the source catalogues employed.

The SCUBA HAlf Degree Extragalactic Survey (SHADES) is a major blank-field extragalactic sub-mm survey underway at the James Clerk Maxwell telescope.SHADES aims to cover half a square degree at 450+850microns to a 4sigma depth of ~8mJy at 850microns. Two fields are being observed, the Subaru/XMM--Newton Deep Field (SXDF) and the Lockman Hole East. The survey has 3 main aims to investigate: i) the cosmic history of massive dust enshrouded star-formation activity, ii) the clustering properties of sub-mm-selected galaxies in order to determine whether these objects could be progenitors of present-day massive ellipticals, and iii) the fraction of sub-mm-selected sources that harbour AGN. To achieve these aims requires the sub-mm data be combined with co-spatial information spanning the radio-to-X-ray frequency range, the resulting extensive multi-wavelength dataset providing complete photometric redshift information accurate to delta z <~0.5 and detailed SEDs for the vast majority of the sub-mm-selected sources. In this paper, the first of a series on SHADES, we present an overview of the motivation for the survey, describe the survey strategy and provide a detailed description of the primary data analysis pipeline. As of Feb. 2004, 720arcmin^2 had been mapped with SCUBA (about 40% of the anticipated final total) to a median 1sigma depth of 2.2mJy per beam at 850microns (25mJy per beam at 450microns), and a source density of 650+/-50 sources/deg^2 >3sigma at 850microns, uncorrected for Eddington bias. A refined re-analysis of the original 8-mJy survey Lockman hole data was carried out in order to evaluate the new data reduction pipeline. Of the 17 most secure sources in the original sample, 12 have been re-confirmed, including 10 of the 11 for which radio IDs were previously secured.(Abridged)

We present Spitzer IRAC (3.6, 4.5, 5.8, 8.0 um) and MIPS (24 um) observations of the SCUBA submillimeter sources and uJy VLA radio sources in a 5'x5' area in the Lockman Hole East region. Out of the ~40 SCUBA/VLA sources in the field, Spitzer counterparts were detected for nearly all except for the few low-weight SCUBA detections. We show that the majority (80-90%) of the detected sources are cold (i.e., starburst-like) infrared-luminous galaxies (L_{IR} > 10^{11} Lsun) at redshift 0.5<z<3.5, whose star-formation rate density (SFRD) is comparable to that of the optically-selected star-forming galaxies.

We present Spitzer observations in five wavebands between 3.6 and 24um of an unbiased sample of 9 luminous, dusty galaxies selected at 1200um by the MAMBO camera on the IRAM 30-m telescope, a population akin to the well-known submm or `SCUBA' galaxies (hereafter SMGs). Owing to the coarse resolution of submm/mm instrumentation, SMGs have traditionally been difficult to identify at other wavelengths. We compare our multi-wavelength catalogs to show that the overlap between 24 and 1200um must be close to complete at these flux levels. We find that all (4/4) of the most secure >=4sigma SMGs have robust >=4sigma counterparts at 1.4GHz, while the fraction drops to 7/9 using all >=3sigma SMGs. We show that combining mid-IR and marginal (>=3sigma) radio detections provides plausible identifications in the remaining cases, enabling us to identify the complete sample. Accretion onto an obscured central engine is betrayed by the shape of the mid-IR continuum emission for several sources, confirming Spitzer's potential to weed out active galaxies. We demonstrate the power of a S(24um)/S(8um) vs S(8um)/S(4.5um) color-color plot as a diagnostic for this purpose. However, we conclude that the majority (~75%) of SMGs have rest-frame mid-/far-IR SEDs commensurate with obscured starbursts. Sensitive 24-um observations are clearly a useful route to identify and characterize reliable counterparts to high-redshift far-IR-bright galaxies, complementing what is possible via deep radio imaging.

We present sub-millimeter statistical detections of galaxies discovered in the 5'x5' Spitzer Early Release Observations (to 4-15 microJy 5 sigma at 3.6-8 microns, 170 microJy at 24 microns) through a stacking analysis of our reanalysed SCUBA 8mJy survey maps, and a Spitzer identification of a new sub-millimeter point source in the 8mJy survey region. For sources detected at 5.8 or 8 microns (154 and 111 sources respectively), we detect positive skews in the sub-millimeter flux distributions at 99.2-99.8% confidence using Kolmogorov-Smirnov tests, at both 850 microns and 450 microns. We also marginally detect the Spitzer 24 micron galaxies at 850 microns at 97% confidence, and place limits on the mean sub-millimeter fluxes of the 3.6 and 4.5 micron sources. Integrating the sub-millimeter fluxes of the Spitzer populations, we find the 5.8 micron galaxies contribute 0.12 +/- 0.05 nW/m^2/sr to the 850 micron background, and 2.4 +/- 0.7 nW/m^2/sr to the 450 micron background; similar contributions are made by the 8 micron-selected sample. We infer that the populations dominating the 5.8 and 8 micron extragalactic background light also contribute around a quarter of the 850 micron background and the majority of the 450 micron background.