Source author record

L. V. E. Koopmans

L. V. E. Koopmans appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.IM astro-ph astro-ph.HE astro-ph.SR astro-ph.EP Databases gr-qc physics.ao-ph

Catalog footprint

What is connected

84works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

The impact of diffuse Galactic emission on direction-independent gain calibration in high-redshift 21 cm observations

This study examines the impact of diffuse Galactic emission (DGE) on sky-based direction-independent (DI) gain calibration using realistic forward simulations of Low-Frequency Array (LOFAR) observations of the high-redshift 21 cm signal of neutral hydrogen during the Epoch of Reionization (EoR). We simulated LOFAR observations between 147 and 159 MHz using a sky model that includes a point source catalog and DGE. The simulated observations were DI-gain calibrated with the point source catalog alone, utilizing the LOFAR-EoR data analysis pipeline. A full power spectrum (PS) analysis was conducted to measure the systematic bias, relative to thermal noise, caused by DI-gain calibration using a point-source-only (PSO) sky model, when applied to simulated data that include both point sources and DGE. The results are compared to a ground truth scenario where both the simulated sky and the calibration model include only point sources. Additionally, the cross-coherence between observation pairs was computed to determine whether DI-gain calibration errors are coherent or incoherent in specific regions of PS space as a function of integration time. We find that DI-gain calibration with a PSO sky model that omits DGE introduces a systematic bias in the PS for $k_{\parallel}$ bins < 0.2 $h\,\mathrm{Mpc}^{-1}$. The PS errors in these bins are coherent in time and frequency; therefore, the resulting bias could be mitigated during the foreground removal step using Gaussian Process Regression, as demonstrated in previous studies. In contrast, errors for $k_{\parallel}$ > 0.2 $h\,\mathrm{Mpc}^{-1}$ are largely incoherent and average down as noise. We conclude that, based on our analysis prior to foreground removal, missing DGE in the sky model during DI-gain calibration is unlikely to be a dominant contributor to the excess noise observed in the current LOFAR-EoR upper limits on the 21 cm signal PS.

preprint2022arXiv

Degree-Scale Galactic Radio Emission at 122 MHz around the North Celestial Pole with LOFAR-AARTFAAC

Aims: Contamination from bright diffuse Galactic thermal and non-thermal radio emission poses crucial challenges in experiments aiming to measure the 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization. If not included in calibration, this diffuse emission can severely impact the analysis and signal extraction in 21-cm experiments. We examine large-scale diffuse Galactic emission at 122~MHz, around the North Celestial Pole, using the AARTFAAC-HBA system. Methods: In this pilot project, we present the first-ever wide-field image produced with a single sub-band of the data recorded with the AARTFAAC-HBA. We demonstrate two methods: multiscale CLEAN and shapelet decomposition, to model the diffuse emission revealed in the image. We use angular power spectrum metrics to quantify different components of the emission and compare the performance of the two diffuse structure modelling approaches. Results: We observe that the point sources dominate the angular power spectrum ($\ell(\ell+1)C_{\ell}/2π\equiv Δ^2(\ell)$) of the emission in the field on scales $\ell\gtrsim 60$ ($\lesssim 3$~degree). The angular power spectrum after subtraction of compact sources is flat within $20\lesssim \ell \lesssim200$ range, suggesting that the residual power is dominated by the diffuse emission on scales $\ell\lesssim200$. The residual diffuse emission has a brightness temperature variance of $Δ^2_{\ell=180} = (145.64 \pm 13.61)~{\rm K}^2$ at 122~MHz on angular scales of 1~degree, and is consistent with a power-law following $C_{\ell}\propto \ell^{-2.0}$ in $20\lesssim \ell \lesssim200$ range. We also find that, in the current setup, the multiscale CLEAN is suitable to model the compact and diffuse structures on a wide range of angular scales, whereas the shapelet decomposition method better models the large scales, which are of the order of a few degrees and wider.

preprint2021arXiv

High-quality strong lens candidates in the final Kilo Degree survey footprint

We present 97 new high-quality strong lensing candidates found in the final $\sim 350\,\rm deg^2$, that completed the full $\sim 1350\,\rm deg^2$ area of the Kilo-Degree Survey (KiDS). Together with our previous findings, the final list of high-quality candidates from KiDS sums up to 268 systems. The new sample is assembled using a new Convolutional Neural Network (CNN) classifier applied to $r$-band (best seeing) and $g,~r,~i$ color-composited images separately. This optimizes the complementarity of the morphology and color information on the identification of strong lensing candidates. We apply the new classifiers to a sample of luminous red galaxies (LRGs) and a sample of bright galaxies (BGs) and select candidates that received a high probability to be a lens from the CNN ($P_{\rm CNN}$). In particular, setting $P_{\rm CNN}>0.8$ for the LRGs, the $1$-band CNN predicts 1213 candidates, while the $3$-band classifier yields 1299 candidates, with only $\sim$30\% overlap. For the BGs, in order to minimize the false positives, we adopt a more conservative threshold, $P_{\rm CNN} >0.9$, for both CNN classifiers. This results in 3740 newly selected objects. The candidates from the two samples are visually inspected by 7 co-authors to finally select 97 "high-quality" lens candidates which received mean scores larger than 6 (on a scale from 0 to 10). We finally discuss the effect of the seeing on the accuracy of CNN classification and possible avenues to increase the efficiency of multi-band classifiers, in preparation of next-generation surveys from ground and space.

preprint2021arXiv

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $μ$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $μ$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit.

preprint2021arXiv

Sub-arcsecond imaging with the International LOFAR Telescope: II. Completion of the LOFAR Long-Baseline Calibrator Survey

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS.

preprint2020arXiv

Cassiopeia A, Cygnus A, Taurus A, and Virgo A at ultra-low radio frequencies

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

preprint2020arXiv

Constraining the intergalactic medium at $z\approx$ 9.1 using LOFAR Epoch of Reionization observations

We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space.

preprint2020arXiv

LOFAR 144-MHz follow-up observations of GW170817

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$σ$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $α^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

preprint2020arXiv

New high-quality strong lens candidates with deep learning in the Kilo Degree Survey

We report new high-quality galaxy scale strong lens candidates found in the Kilo Degree Survey data release 4 using Machine Learning. We have developed a new Convolutional Neural Network (CNN) classifier to search for gravitational arcs, following the prescription by \cite{2019MNRAS.484.3879P} and using only $r-$band images. We have applied the CNN to two "predictive samples": a Luminous red galaxy (LRG) and a "bright galaxy" (BG) sample ($r<21$). We have found 286 new high probability candidates, 133 from the LRG sample and 153 from the BG sample. We have then ranked these candidates based on a value that combines the CNN likelihood to be a lens and the human score resulting from visual inspection (P-value) and we present here the highest 82 ranked candidates with P-values $\ge 0.5$. All these high-quality candidates have obvious arc or point-like features around the central red defector. Moreover, we define the best 26 objects, all with scores P-values $\ge 0.7$ as a "golden sample" of candidates. This sample is expected to contain very few false positives and thus it is suitable for follow-up observations. The new lens candidates come partially from the the more extended footprint adopted here with respect to the previous analyses, partially from a larger predictive sample (also including the BG sample). These results show that machine learning tools are very promising to find strong lenses in large surveys and more candidates that can be found by enlarging the predictive samples beyond the standard assumption of LRGs. In the future, we plan to apply our CNN to the data from next-generation surveys such as the Large Synoptic Survey Telescope, Euclid, and the Chinese Space Station Optical Survey.

preprint2020arXiv

SHARP -- VI. Evidence for CO (1-0) molecular gas extended on kpc-scales in AGN star forming galaxies at high redshift

We present a study of the stellar host galaxy, CO (1$-$0) molecular gas distribution and AGN emission on 50 to 500 pc-scales of the gravitationally lensed dust-obscured AGN MG J0751+2716 and JVAS B1938+666 at redshifts 3.200 and 2.059, respectively. By correcting for the lensing distortion using a grid-based lens modelling technique, we spatially locate the different emitting regions in the source plane for the first time. Both AGN host galaxies have 300 to 500 pc-scale size and surface brightness consistent with a bulge/pseudo-bulge, and 2 kpc-scale AGN radio jets that are embedded in extended molecular gas reservoirs that are 5 to 20 kpc in size. The CO (1$-$0) velocity fields show structures possibly associated with discs (elongated velocity gradients) and interacting objects (off-axis velocity components). There is evidence for a decrement in the CO (1$-$0) surface brightness at the location of the host galaxy, which may indicate radiative feedback from the AGN, or offset star formation.We find CO-H$_2$ conversion factors of around $α_{\rm CO} = 1.5\pm0.5$ (K km s$^{-1}$ pc$^2$)$^{-1}$, molecular gas masses of $> 3\times10^{10}$ M$_{\odot}$, dynamical masses of $\sim 10^{11}$ M$_{\odot}$ and gas fractions of around 60 per cent. The intrinsic CO line luminosities are comparable to those of unobscured AGN and dusty star-forming galaxies at similar redshifts, but the infrared luminosities are lower, suggesting that the targets are less efficient at forming stars. Therefore, they may belong to the AGN feedback phase predicted by galaxy formation models, because they are not efficiently forming stars considering their large amount of molecular gas.

preprint2020arXiv

TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

preprint2020arXiv

The e-MERLIN Galaxy Evolution Survey (e-MERGE): Overview and Survey Description

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5μ$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

preprint2016arXiv

A hierarchical Bayesian approach for reconstructing the Initial Mass Function of Single Stellar Populations

Recent studies based on the integrated light of distant galaxies suggest that the initial mass function (IMF) might not be universal. Variations of the IMF with galaxy type and/or formation time may have important consequences for our understanding of galaxy evolution. We have developed a new stellar population synthesis (SPS) code specifically designed to reconstruct the IMF. We implement a novel approach combining regularization with hierarchical Bayesian inference. Within this approach we use a parametrized IMF prior to regulate a direct inference of the IMF. This direct inference gives more freedom to the IMF and allows the model to deviate from parametrized models when demanded by the data. We use Markov Chain Monte Carlo sampling techniques to reconstruct the best parameters for the IMF prior, the age, and the metallicity of a single stellar population. We present our code and apply our model to a number of mock single stellar populations with different ages, metallicities, and IMFs. When systematic uncertainties are not significant, we are able to reconstruct the input parameters that were used to create the mock populations. Our results show that if systematic uncertainties do play a role, this may introduce a bias on the results. Therefore, it is important to objectively compare different ingredients of SPS models. Through its Bayesian framework, our model is well-suited for this.

preprint2016arXiv

Polarization leakage in epoch of reionization windows: II. Primary beam model and direction dependent calibration

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

preprint2016arXiv

Probing Ionospheric Structures using the LOFAR radio telescope

LOFAR is the LOw Frequency Radio interferometer ARray located at mid-latitude ($52^{\circ} 53'N$). Here, we present results on ionospheric structures derived from 29 LOFAR nighttime observations during the winters of 2012/2013 and 2013/2014. We show that LOFAR is able to determine differential ionospheric TEC values with an accuracy better than 1 mTECU over distances ranging between 1 and 100 km. For all observations the power law behavior of the phase structure function is confirmed over a long range of baseline lengths, between $1$ and $80$ km, with a slope that is in general larger than the $5/3$ expected for pure Kolmogorov turbulence. The measured average slope is $1.89$ with a one standard deviation spread of $0.1$. The diffractive scale, i.e. the length scale where the phase variance is $1\, \mathrm{rad^2}$, is shown to be an easily obtained single number that represents the ionospheric quality of a radio interferometric observation. A small diffractive scale is equivalent to high phase variability over the field of view as well as a short time coherence of the signal, which limits calibration and imaging quality. For the studied observations the diffractive scales at $150$ MHz vary between $3.5$ and $30\,$ km. A diffractive scale above $5$ km, pertinent to about $90 \%$ of the observations, is considered sufficient for the high dynamic range imaging needed for the LOFAR Epoch of Reionization project. For most nights the ionospheric irregularities were anisotropic, with the structures being aligned with the Earth magnetic field in about $60\%$ of the observations.

preprint2016arXiv

SHARP - II. Mass structure in strong lenses is not necessarily dark matter substructure: A flux ratio anomaly from an edge-on disc

Gravitational lens flux-ratio anomalies provide a powerful technique for measuring dark matter substructure in distant galaxies. However, before using these flux-ratio anomalies to test galaxy formation models, it is imperative to ascertain that the given anomalies are indeed due to the presence of dark matter substructure and not due to some other component of the lensing galaxy halo or to propagation effects. Here we present the case of CLASS~B1555+375, which has a strong radio-wavelength flux-ratio anomaly. Our high-resolution near-infrared Keck~II adaptive optics imaging and archival Hubble Space Telescope data reveal the lensing galaxy in this system to have a clear edge-on disc component that crosses directly over the pair of images that exhibit the flux-ratio anomaly. We find that simple models that include the disc can reproduce the cm-wavelength flux-ratio anomaly without requiring additional dark matter substructure. Although further studies are required, our results suggest the assumption that all flux-ratio anomalies are due to a population of dark matter sub-haloes may be incorrect, and analyses that do not account for the full complexity of the lens macro-model may overestimate the substructure mass fraction in massive lensing galaxies.

preprint2016arXiv

The Herschel-ATLAS: a sample of 500μm-selected lensed galaxies over 600 square degrees

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500μm extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500μm and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

preprint2015arXiv

A Bayesian analysis of redshifted 21-cm HI signal and foregrounds: Simulations for LOFAR

Observations of the EoR with the 21-cm hyperfine emission of neutral hydrogen (HI) promise to open an entirely new window onto the formation of the first stars, galaxies and accreting black holes. In order to characterize the weak 21-cm signal, we need to develop imaging techniques which can reconstruct the extended emission very precisely. Here, we present an inversion technique for LOFAR baselines at NCP, based on a Bayesian formalism with optimal spatial regularization, which is used to reconstruct the diffuse foreground map directly from the simulated visibility data. We notice the spatial regularization de-noises the images to a large extent, allowing one to recover the 21-cm power-spectrum over a considerable $k_{\perp}-k_{\para}$ space in the range of $0.03\,{\rm Mpc^{-1}}<k_{\perp}<0.19\,{\rm Mpc^{-1}}$ and $0.14\,{\rm Mpc^{-1}}<k_{\para}<0.35\,{\rm Mpc^{-1}}$ without subtracting the noise power-spectrum. We find that, in combination with using the GMCA, a non-parametric foreground removal technique, we can mostly recover the spherically average power-spectrum within $2σ$ statistical fluctuations for an input Gaussian random rms noise level of $60 \, {\rm mK}$ in the maps after 600 hrs of integration over a $10 \, {\rm MHz}$ bandwidth.

preprint2015arXiv

Are the total mass density and the low-mass end slope of the IMF anti-correlated?

We conduct a detailed lensing, dynamics and stellar population analysis of nine massive lens early-type galaxies (ETGs) from the X-Shooter Lens Survey (XLENS). Combining gravitational lensing constraints from HST imaging with spatially-resolved kinematics and line-indices constraints from VLT X-Shooter (XSH) spectra, we infer the low-mass slope and the low cut-off mass of the stellar Initial Mass Function (IMF): $x_{250}=2.37^{+0.12}_{-0.12}$ and $M_{{\rm low}, 250}= 0.131^{+0.023}_{-0.026}\, M_{\odot}$, respectively, for a reference point with $σ\equiv 250\, {\rm kms}^{-1}$ and R$_{\rm eff} \equiv 10$ kpc. All the XLENS systems are consistent with an IMF slope steeper than Milky Way-like. We find no significant correlations between IMF slope and any other quantity, except for an anti-correlation between total dynamical mass density and low-mass IMF slope at the 87% CL [$dx/d\log(ρ)$ = $ -0.19^{+0.15}_{-0.15}$]. This anti-correlation is consistent with the low redshift lenses found by Smith et al. (2015) that have high velocity dispersions and high stellar mass densities but surprisingly shallow IMF slopes.

preprint2015arXiv

Galactic interstellar filaments as probed by LOFAR and Planck

Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structures along the line of sight, whereas the Planck dust polarization maps probe the orientation of the sky projected magnetic field component. Hence, no apparent correlation between the two is expected. Here we report a surprising, yet clear, correlation between the filamentary structures, detected with LOFAR, and the magnetic field orientation, probed by the Planck satellite. This finding points to a common, yet unclear, physical origin of the two measurements in this specific area in the sky. A number of follow-up multi- frequency studies are proposed to shed light on this unexpected finding.

preprint2015arXiv

Linear polarization structures in LOFAR observations of the interstellar medium in the 3C196 field

This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. We have used the high band antennas (HBA) of LOFAR to image this region and Rotation Measure (RM) synthesis to unravel the distribution of polarized structures in Faraday depth. The brightness temperature of the detected Galactic emission is $5-15~{\rm K}$ in polarized intensity and covers the range from -3 to +8 ${\rm rad~m^{-2}}$ in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of $+0.5~{\rm rad~m^{-2}}$ running from north to south, right through the centre of the field and parallel to the Galactic plane. There is also an interesting system of linear depolarization canals conspicuous in an image showing the peaks of Faraday spectra. We used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. For the first time, we see some common morphology in the RM cubes made at 150 and 350~{\rm MHz}. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies and previous LOFAR observations. Based on our results, we determined physical parameters of the ISM and proposed a simple model that may explain the observed distribution of the intervening magneto-ionic medium. The mean line-of-sight magnetic field component, $B_\parallel$, is determined to be $0.3\pm0.1~{\rm μG}$ and its spatial variation across the 3C196 field is $0.1~{\rm μG}$. The filamentary structure is probably an ionized filament in the ISM, located somewhere within the Local Bubble. This filamentary structure shows an excess in thermal electron density ($n_e B_\parallel>6.2~{\rm cm^{-3}μG}$) compared to its surroundings.

preprint2015arXiv

Polarization leakage in Epoch of Reionization windows: I. LOFAR observations of the 3C196 field

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

preprint2015arXiv

Scintillation noise in widefield radio interferometry

In this paper, we consider random phase fluctuations imposed during wave propagation through a turbulent plasma (e.g. ionosphere) as a source of additional noise in interferometric visibilities. We derive expressions for visibility variance for the wide field of view case (FOV$\sim10$ deg) by computing the statistics of Fresnel diffraction from a stochastic plasma, and provide an intuitive understanding. For typical ionospheric conditions (diffractive scale $\sim 5-20$ km at $150$ MHz), we show that the resulting ionospheric `scintillation noise' can be a dominant source of uncertainty at low frequencies ($ν\lesssim 200$ MHz). Consequently, low frequency widefield radio interferometers must take this source of uncertainty into account in their sensitivity analysis. We also discuss the spatial, temporal, and spectral coherence properties of scintillation noise that determine its magnitude in deep integrations, and influence prospects for its mitigation via calibration or filtering.

preprint2015arXiv

Scintillation noise power spectrum and its impact on high redshift 21-cm observations

Visibility scintillation resulting from wave propagation through the turbulent ionosphere can be an important sources of noise at low radio frequencies ($ν\lesssim 200$ MHz). Many low frequency experiments are underway to detect the power spectrum of brightness temperature fluctuations of the neutral-hydrogen $21$-cm signal from the Epoch of Reionization (EOR: $12\gtrsim z\gtrsim 7$, $100\lesssim ν\lesssim 175$ MHz). In this paper, we derive scintillation noise power-spectra in such experiments while taking into account the effects of typical data processing operations such as self-calibration and Fourier synthesis. We find that for minimally redundant arrays such as LOFAR and MWA, scintillation noise is of the same order of magnitude as thermal noise, has a spectral coherence dictated by stretching of the snapshot $uv$-coverage with frequency, and thus is confined to the well known wedge-like structure in the cylindrical ($2$-dimensional) power spectrum space. Compact, fully redundant ($d_{\rm core}\lesssim r_{\rm F} \approx 300$ m at $150$ MHz) arrays such as HERA and SKA-LOW (core) will be scintillation noise dominated at all baselines, but the spatial and frequency coherence of this noise will allow it to be removed along with spectrally smooth foregrounds.

preprint2015arXiv

Simulating the 21cm forest detectable with LOFAR and SKA in the spectra of high-z GRBs

We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.

preprint2015arXiv

Strong gravitational lensing with the SKA

Strong gravitational lenses provide an important tool to measure masses in the distant Universe, thus testing models for galaxy formation and dark matter; to investigate structure at the Epoch of Reionization; and to measure the Hubble constant and possibly w as a function of redshift. However, the limiting factor in all of these studies has been the currently small samples of known gravitational lenses (~10^2). The era of the SKA will transform our understanding of the Universe with gravitational lensing, particularly at radio wavelengths where the number of known gravitational lenses will increase to ~10^5. Here we discuss the technical requirements, expected outcomes and main scientific goals of a survey for strong gravitational lensing with the SKA. We find that an all-sky (3pi sr) survey carried out with the SKA1-MID array at an angular resolution of 0.25-0.5 arcsec and to a depth of 3 microJy / beam is required for studies of galaxy formation and cosmology with gravitational lensing. In addition, the capability to carryout VLBI with the SKA1 is required for tests of dark matter and studies of supermassive black holes at high redshift to be made using gravitational lensing.

preprint2015arXiv

Strong lens search in the ESO public Survey KiDS

We have started a systematic search of strong lens candidates in the ESO public survey KiDS based on the visual inspection of massive galaxies in the redshift range $0.1<z<0.5$. As a pilot program we have inspected 100 sq. deg., which overlap with SDSS and where there are known lenses to use as a control sample. Taking advantage of the superb image quality of VST/OmegaCAM, the colour information and accurate model subtracted images, we have found 18 new lens candidates, for which spectroscopic confirmation will be needed to confirm their lensing nature and study the mass profile of the lensing galaxies.

preprint2015arXiv

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

preprint2015arXiv

The X-shooter Lens Survey - II. Sample presentation and spatially resolved kinematics

We present the X-shooter Lens Survey (XLENS) data. The main goal of XLENS is to disentangle the stellar and dark matter content of massive early-type galaxies (ETGs), through combined strong gravitational lensing, dynamics and spectroscopic stellar population studies. The sample consists of 11 lens galaxies covering the redshift range from $0.1$ to $0.45$ and having stellar velocity dispersions between $250$ and $380\,\mathrm{km}\,\mathrm{s}^{-1}$. All galaxies have multi-band, high-quality HST imaging. We have obtained long-slit spectra of the lens galaxies with X-shooter on the VLT. We are able to disentangle the dark and luminous mass components by combining lensing and extended kinematics data-sets, and we are also able to precisely constrain stellar mass-to-light ratios and infer the value of the low-mass cut-off of the IMF, by adding spectroscopic stellar population information. Our goal is to correlate these IMF parameters with ETG masses and investigate the relation between baryonic and non-baryonic matter during the mass assembly and structure formation processes. In this paper we provide an overview of the survey, highlighting its scientific motivations, main goals and techniques. We present the current sample, briefly describing the data reduction and analysis process, and we present the first results on spatially resolved kinematics.

preprint2014arXiv

Cosmology from gravitational lens time delays and Planck data

Under the assumption of a flat Lambda-CDM cosmology, recent data from the Planck satellite point toward a Hubble constant that is in tension with that measured by gravitational lens time delays and by the local distance ladder. Prosaically, this difference could arise from unknown systematic uncertainties in some of the measurements. More interestingly -- if systematics were ruled out -- resolving the tension would require a departure from the flat Lambda-CDM cosmology, introducing for example a modest amount of spatial curvature, or a non-trivial dark energy equation of state. To begin to address these issues, we present here an analysis of the gravitational lens RXJ1131-1231 that is improved in one particular regard: we examine the issue of systematic error introduced by an assumed lens model density profile. We use more flexible gravitational lens models with baryonic and dark matter components, and find that the exquisite Hubble Space Telescope image with thousands of intensity pixels in the Einstein ring and the stellar velocity dispersion of the lens contain sufficient information to constrain these more flexible models. The total uncertainty on the time-delay distance is 6.6% for a single system. We proceed to combine our improved time-delay distance measurements with the WMAP9 and Planck posteriors. In an open Lambda-CDM model, the data for RXJ1131-1231 in combination with Planck favor a flat universe with Omega_k = 0.00+0.01/-0.02 (68% CI). In a flat wCDM model, the combination of RXJ1131-1231 and Planck yields w = -1.52+0.19/-0.20 (68% CI).

preprint2014arXiv

Galaxy Masses

Galaxy masses play a fundamental role in our understanding of structure formation models. This review addresses the variety and reliability of mass estimators that pertain to stars, gas, and dark matter. The different sections on masses from stellar populations, dynamical masses of gas-rich and gas-poor galaxies, with some attention paid to our Milky Way, and masses from weak and strong lensing methods, all provide review material on galaxy masses in a self-consistent manner.

preprint2014arXiv

Inference of the Cold Dark Matter substructure mass function at z=0.2 using strong gravitational lenses

We present the results of a search for galaxy substructures in a sample of 11 gravitational lens galaxies from the Sloan Lens ACS Survey. We find no significant detection of mass clumps, except for a luminous satellite in the system SDSS J0956+5110. We use these non-detections, in combination with a previous detection in the system SDSS J0946+1006, to derive constraints on the substructure mass function in massive early-type host galaxies with an average redshift z ~ 0.2 and an average velocity dispersion of 270 km/s. We perform a Bayesian inference on the substructure mass function, within a median region of about 32 kpc squared around the Einstein radius (~4.2 kpc). We infer a mean projected substructure mass fraction $f = 0.0076^{+0.0208}_{-0.0052}$ at the 68 percent confidence level and a substructure mass function slope $α$ < 2.93 at the 95 percent confidence level for a uniform prior probability density on alpha. For a Gaussian prior based on Cold Dark Matter (CDM) simulations, we infer $f = 0 .0064^{+0.0080}_{-0.0042}$ and a slope of $α$ = 1.90$^{+0.098}_{-0.098}$ at the 68 percent confidence level. Since only one substructure was detected in the full sample, we have little information on the mass function slope, which is therefore poorly constrained (i.e. the Bayes factor shows no positive preference for any of the two models).The inferred fraction is consistent with the expectations from CDM simulations and with inference from flux ratio anomalies at the 68 percent confidence level.

preprint2014arXiv

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

preprint2014arXiv

Lunar occultation of the diffuse radio sky: LOFAR measurements between 35 and 80 MHz

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<ν<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

preprint2014arXiv

Measuring gravitational lens time delays using low-resolution radio monitoring observations

Obtaining lensing time delay measurements requires long-term monitoring campaigns with a high enough resolution (< 1 arcsec) to separate the multiple images. In the radio, a limited number of high-resolution interferometer arrays make these observations difficult to schedule. To overcome this problem, we propose a technique for measuring gravitational time delays which relies on monitoring the total flux density with low-resolution but high-sensitivity radio telescopes to follow the variation of the brighter image. This is then used to trigger high-resolution observations in optimal numbers which then reveal the variation in the fainter image. We present simulations to assess the efficiency of this method together with a pilot project observing radio lens systems with the Westerbork Synthesis Radio Telescope (WSRT) to trigger Very Large Array (VLA) observations. This new method is promising for measuring time delays because it uses relatively small amounts of time on high-resolution telescopes. This will be important because instruments that have high sensitivity but limited resolution, together with an optimum usage of followup high-resolution observations from appropriate radio telescopes may in the future be useful for gravitational lensing time delay measurements by means of this new method.

preprint2014arXiv

Radio Monitoring Campaigns of Six Strongly Lensed Quasars

We observed six strongly lensed, radio-loud quasars (MG 0414+0534, CLASS B0712+472, JVAS B1030+074, CLASS B1127+385, CLASS B1152+199, and JVAS B1938+666) in order to identify systems suitable for measuring cosmological parameters using time delays between their multiple images. These systems are in standard two and four image configurations, with B1938 having a faint secondary pair of images. Two separate monitoring campaigns were carried out using the VLA and upgraded JVLA. Lightcurves were extracted for each individual lensed image and analyzed for signs of intrinsic variability. While it was not possible to measure time delays from these data, $χ^2$-based and structure function tests found evidence for variability in a majority of the lightcurves. B0712 and B1030 had particularly strong variations, exhibiting linear flux trends. These results show that most of these systems should be targeted with followup monitoring campaigns, especially B0712 and B1030.

preprint2013arXiv

Calibrating High-Precision Faraday Rotation Measurements for LOFAR and the Next Generation of Low-Frequency Radio Telescopes

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

preprint2013arXiv

Chromatic effects in the 21 cm global signal from the cosmic dawn

The redshifted 21 cm brightness distribution from neutral hydrogen is a promising probe into the cosmic dark ages, cosmic dawn, and re-ionization. LOFAR's Low Band Antennas (LBA) may be used in the frequency range 45 MHz to 85 MHz (30>z>16) to measure the sky averaged redshifted 21 cm brightness temperature as a function of frequency, or equivalently, cosmic redshift. These low frequencies are affected by strong Galactic foreground emission that is observed through frequency dependent ionospheric and antenna beam distortions which lead to chromatic mixing of spatial structure into spectral structure. Using simple models, we show that (i) the additional antenna temperature due to ionospheric refraction and absorption are at a \sim 1% level--- 2 to 3 orders of magnitude higher than the expected 21 cm signal, and have an approximate ν^{-2} dependence, (ii) ionospheric refraction leads to a knee-like modulation on the sky spectrum at ν\approx 4\times plasma frequency. Using more realistic simulations, we show that in the measured sky spectrum, more than 50% of the 21 cm signal variance can be lost to confusion from foregrounds and chromatic effects. We conclude that foregrounds and chromatic mixing may not be subtracted as generic functions of frequency as previously thought, but must rather be carefully modeled using additional priors and interferometric measurements.

preprint2013arXiv

Dark energy with gravitational lens time delays

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

preprint2013arXiv

Initial deep LOFAR observations of Epoch of Reionization windows: I. The North Celestial Pole

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

preprint2013arXiv

LEDDB: LOFAR Epoch of Reionization Diagnostic Database

One of the key science projects of the Low-Frequency Array (LOFAR) is the detection of the cosmological signal coming from the Epoch of Reionization (EoR). Here we present the LOFAR EoR Diagnostic Database (LEDDB) that is used in the storage, management, processing and analysis of the LOFAR EoR observations. It stores referencing information of the observations and diagnostic parameters extracted from their calibration. This stored data is used to ease the pipeline processing, monitor the performance of the telescope and visualize the diagnostic parameters which facilitates the analysis of the several contamination effects on the signals. It is implemented with PostgreSQL and accessed through the psycopg2 python module. We have developed a very flexible query engine, which is used by a web user interface to access the database, and a very extensive set of tools for the visualization of the diagnostic parameters through all their multiple dimensions.

preprint2013arXiv

LOFAR insights into the epoch of reionization from the cross power spectrum of 21cm emission and galaxies

Using a combination of N-body simulations, semi-analytic models and radiative transfer calculations, we have estimated the theoretical cross power spectrum between galaxies and the 21cm emission from neutral hydrogen during the epoch of reionization. In accordance with previous studies, we find that the 21cm emission is initially correlated with halos on large scales (> 30 Mpc), anti-correlated on intermediate (~ 5 Mpc), and uncorrelated on small (< 3 Mpc) scales. This picture quickly changes as reionization proceeds and the two fields become anti-correlated on large scales. The normalization of the cross power spectrum can be used to set constraints on the average neutral fraction in the intergalactic medium and its shape can be a tool to study the topology of reionization. When we apply a drop-out technique to select galaxies and add to the 21cm signal the noise expected from the LOFAR telescope, we find that while the normalization of the cross power spectrum remains a useful tool for probing reionization, its shape becomes too noisy to be informative. On the other hand, for a Lyalpha Emitter (LAE) survey both the normalization and the shape of the cross power spectrum are suitable probes of reionization. A closer look at a specific planned LAE observing program using Subaru Hyper-Suprime Cam reveals concerns about the strength of the 21cm signal at the planned redshifts. If the ionized fraction at z ~ 7 is lower that the one estimated here, then using the cross power spectrum may be a useful exercise given that at higher redshifts and neutral fractions it is able to distinguish between two toy models with different topologies.

preprint2013arXiv

LOFAR: The LOw-Frequency ARray

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

preprint2013arXiv

Probing reionization with LOFAR using 21-cm redshift space distortions

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios.

preprint2013arXiv

Studying Galactic interstellar turbulence through fluctuations in synchrotron emission: First LOFAR Galactic foreground detection

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

preprint2013arXiv

The brightness and spatial distributions of terrestrial radio sources

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

preprint2013arXiv

The Scale of the Problem : Recovering Images of Reionization with GMCA

The accurate and precise removal of 21-cm foregrounds from Epoch of Reionization redshifted 21-cm emission data is essential if we are to gain insight into an unexplored cosmological era. We apply a non-parametric technique, Generalized Morphological Component Analysis or GMCA, to simulated LOFAR-EoR data and show that it has the ability to clean the foregrounds with high accuracy. We recover the 21-cm 1D, 2D and 3D power spectra with high accuracy across an impressive range of frequencies and scales. We show that GMCA preserves the 21-cm phase information, especially when the smallest spatial scale data is discarded. While it has been shown that LOFAR-EoR image recovery is theoretically possible using image smoothing, we add that wavelet decomposition is an efficient way of recovering 21-cm signal maps to the same or greater order of accuracy with more flexibility. By comparing the GMCA output residual maps (equal to the noise, 21-cm signal and any foreground fitting errors) with the 21-cm maps at one frequency and discarding the smaller wavelet scale information, we find a correlation coefficient of 0.689, compared to 0.588 for the equivalently smoothed image. Considering only the central 50% of the maps, these coefficients improve to 0.905 and 0.605 respectively and we conclude that wavelet decomposition is a significantly more powerful method to denoise reconstructed 21-cm maps than smoothing.

preprint2013arXiv

Two accurate time-delay distances from strong lensing: Implications for cosmology

Strong gravitational lenses with measured time delays between the multiple images and models of the lens mass distribution allow a one-step determination of the time-delay distance, and thus a measure of cosmological parameters. We present a blind analysis of the gravitational lens RXJ1131-1231 incorporating (1) the newly measured time delays from COSMOGRAIL, (2) archival HST imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323+/-20km/s based on Keck spectroscopy, and (4) a characterization of the line-of-sight structures via observations of the lens' environment and ray tracing through the Millennium Simulation. Our blind analysis is designed to prevent experimenter bias. The joint analysis of the data sets allows a time-delay distance measurement to 6% precision that takes into account all known systematic uncertainties. In combination with the WMAP7 data set in flat wCDM cosmology, our unblinded cosmological constraints for RXJ1131-1231 are: H_0=80.0+5.8/-5.7km/s/Mpc, OmegaDE=0.79+/-0.03 and w=-1.25+0.17/-0.21. We find the results to be statistically consistent with those from the analysis of the gravitational lens B1608+656. The joint constraints from the two lenses and WMAP7 are H_0=75.2+4.4/-4.2km/s/Mpc, OmegaDE=0.76+0.02/-0.03 and w=-1.14+0.17/-0.20 in flat wCDM, and H_0=73.1+2.4/-3.6km/s/Mpc, OmegaL=0.75+0.01/-0.02 and OmegaK=0.003+0.005/-0.006 in open LCDM. Time-delay lenses constrain especially tightly the Hubble constant (5.7% and 4.0% respectively in wCDM and open LCDM) and curvature of the universe. They complement well other cosmological probes, and provide an independent check of unknown systematics. Our measurement of the Hubble constant is completely independent of those based on the local distance ladder method, providing an important consistency check of the standard cosmological model and of general relativity.

preprint2012arXiv

Disentangling Baryons and Dark Matter in the Spiral Gravitational Lens B1933+503

Measuring the relative mass contributions of luminous and dark matter in spiral galaxies is important for understanding their formation and evolution. The combination of a galaxy rotation curve and strong lensing is a powerful way to break the disk-halo degeneracy that is inherent in each of the methods individually. We present an analysis of the 10-image radio spiral lens B1933+503 at z_l=0.755, incorporating (1) new global VLBI observations, (2) new adaptive-optics assisted K-band imaging, (3) new spectroscopic observations for the lens galaxy rotation curve and the source redshift. We construct a three-dimensionally axisymmetric mass distribution with 3 components: an exponential profile for the disk, a point mass for the bulge, and an NFW profile for the halo. The mass model is simultaneously fitted to the kinematics and the lensing data. The NFW halo needs to be oblate with a flattening of a/c=0.33^{+0.07}_{-0.05} to be consistent with the radio data. This suggests that baryons are effective at making the halos oblate near the center. The lensing and kinematics analysis probe the inner ~10 kpc of the galaxy, and we obtain a lower limit on the halo scale radius of 16 kpc (95% CI). The dark matter mass fraction inside a sphere with a radius of 2.2 disk scale lengths is f_{DM,2.2}=0.43^{+0.10}_{-0.09}. The contribution of the disk to the total circular velocity at 2.2 disk scale lengths is 0.76^{+0.05}_{-0.06}, suggesting that the disk is marginally submaximal. The stellar mass of the disk from our modeling is log_{10}(M_{*}/M_{sun}) = 11.06^{+0.09}_{-0.11} assuming that the cold gas contributes ~20% to the total disk mass. In comparison to the stellar masses estimated from stellar population synthesis models, the stellar initial mass function of Chabrier is preferred to that of Salpeter by a probability factor of 7.2.

preprint2012arXiv

Evidence For A Mild Steepening And Bottom-Heavy IMF In Massive Galaxies From Sodium And Titanium-Oxide Indicators

We measure equivalent widths (EW) - focussing on two unique features (NaI and TiO2) of low-mass stars (<0.3M\odot) - for luminous red galaxy spectra from the the Sloan Digital Sky Survey (SDSS) and X-Shooter Lens Survey (XLENS) in order to study the low-mass end of the initial mass function (IMF). We compare these EWs to those derived from simple stellar population models computed with different IMFs, ages, [α/Fe], and elemental abundances. We find that models are able to simultaneously reproduce the observed NaD λ5895 and Na I λ8190 features for lower-mass (\sim σ\ast) early-type galaxies (ETGs) but deviate increasingly for more massive ETGs, due do strongly mismatching NaD EWs. The TiO2 λ6230 and the Na I λ8190 features together appear to be a powerful IMF diagnostic, with age and metallicity effects orthogonal to the effect of IMF. We find that both features correlate strongly with galaxy velocity dispersion. The XLENS ETG (SDSSJ0912+0029) and an SDSS ETG (SDSSJ0041-0914) appear to require both an extreme dwarf-rich IMF and a high sodium enhancement ([Na/Fe] = +0.4). In addition, lensing constraints on the total mass of the XLENS system within its Einstein radius limit a bottom-heavy IMF with a power-law slope to x \leq 3.0 at the 90% C.L. We conclude that NaI and TiO features, in comparison with state-of-the-art SSP models, suggest a mildly steepening IMF from Salpeter (dn/dm \propto m-x with x = 2.35) to x \approx 3.0 for ETGs in the range σ = 200 - 335 km s-1.

preprint2012arXiv

First LOFAR observations at very low frequencies of cluster-scale non-thermal emission: the case of Abell 2256

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

preprint2012arXiv

Foreground Removal using FastICA: A Showcase of LOFAR-EoR

We introduce a new implementation of the FastICA algorithm on simulated LOFAR EoR data with the aim of accurately removing the foregrounds and extracting the 21-cm reionization signal. We find that the method successfully removes the foregrounds with an average fitting error of 0.5 per cent and that the 2D and 3D power spectra are recovered across the frequency range. We find that for scales above several PSF scales the 21-cm variance is successfully recovered though there is evidence of noise leakage into the reconstructed foreground components. We find that this blind independent component analysis technique provides encouraging results without the danger of prior foreground assumptions.

preprint2012arXiv

Gravitational detection of a low-mass dark satellite at cosmological distance

The mass-function of dwarf satellite galaxies that are observed around Local Group galaxies substantially differs from simulations based on cold dark matter: the simulations predict many more dwarf galaxies than are seen. The Local Group, however, may be anomalous in this regard. A massive dark satellite in an early-type lens galaxy at z = 0.222 was recently found using a new method based on gravitational lensing, suggesting that the mass fraction contained in substructure could be higher than is predicted from simulations. The lack of very low mass detections, however, prohibited any constraint on their mass function. Here we report the presence of a 1.9 +/- 0.1 x 10^8 M_sun dark satellite in the Einstein-ring system JVAS B1938+666 at z = 0.881, where M_sun denotes solar mass. This satellite galaxy has a mass similar to the Sagittarius galaxy, which is a satellite of the Milky Way. We determine the logarithmic slope of the mass function for substructure beyond the local Universe to be alpha = 1.1^+0.6_-0.4, with an average mass-fraction of f = 3.3^+3.6_-1.8 %, by combining data on both of these recently discovered galaxies. Our results are consistent with the predictions from cold dark matter simulations at the 95 per cent confidence level, and therefore agree with the view that galaxies formed hierarchically in a Universe composed of cold dark matter.

preprint2012arXiv

Imaging neutral hydrogen on large-scales during the Epoch of Reionization with LOFAR

The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (\sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it hard to disentangle the two on such large scales. The poor sensitivity on small scales on the one hand, and the foregrounds effect on large scales on the other hand, make direct imaging of the Epoch of Reionization of the Universe very difficult, and detection of the signal therefore is expected to be statistical.Despite these hurdles, in this paper we argue that for many reionization scenarios low resolution images could be obtained from the expected data. This is because at the later stages of the process one still finds very large pockets of neutral regions in the IGM, reflecting the clustering of the large-scale structure, which stays strong up to scales of \sim 120 comoving Mpc/h (\sim 1 degree). The coherence of the emission on those scales allows us to reach sufficient S/N (\sim 3) so as to obtain reionization 21 cm images. Such images will be extremely valuable for answering many cosmological questions but above all they will be a very powerful tool to test our control of the systematics in the data. The existence of this typical scale (\sim 120 comoving Mpc/h) also argues for designing future EoR experiments, e.g., with SKA, with a field of view of at least 4 degree.

preprint2012arXiv

Prospects for detecting the 21cm forest from the diffuse intergalactic medium with LOFAR

We discuss the feasibility of the detection of the 21cm forest in the diffuse IGM with the radio telescope LOFAR. The optical depth to the 21cm line has been derived using simulations of reionization which include detailed radiative transfer of ionizing photons. We find that the spectra from reionization models with similar total comoving hydrogen ionizing emissivity but different frequency distribution look remarkably similar. Thus, unless the reionization histories are very different from each other (e.g. a predominance of UV vs. x-ray heating) we do not expect to distinguish them by means of observations of the 21cm forest. Because the presence of a strong x-ray background would make the detection of 21cm line absorption impossible, the lack of absorption could be used as a probe of the presence/intensity of the x-ray background and the thermal history of the universe. Along a random line of sight LOFAR could detect a global suppression of the spectrum from z>12, when the IGM is still mostly neutral and cold, in contrast with the more well-defined, albeit broad, absorption features visible at lower redshift. Sharp, strong absorption features associated with rare, high density pockets of gas could be detected also at z~7 along preferential lines of sight.

preprint2012arXiv

SHARP - I. A high-resolution multi-band view of the infra-red Einstein ring of JVAS B1938+666

We present new mass models for the gravitational lens system B1938+666, using multi-wavelength data acquired from Keck adaptive optics (AO) and Hubble Space Telescope (HST) observations. These models are the first results from the Strong-lensing at High Angular Resolution Program (SHARP), a project designed to study known quadruple-image and Einstein ring lenses using high-resolution imaging, in order to probe their mass distributions in unprecedented detail. Here, we specifically highlight differences between AO- and HST-derived lens models, finding that -- at least when the lens and source galaxies are both bright and red, and the system has a high degree of circular symmetry -- AO-derived models place significantly tighter constraints on model parameters. Using this improved precision, we infer important physical properties about the B1938+666 system, including the mass density slope of the lensing galaxy (gamma = 2.045), the projected dark matter mass fraction within the Einstein radius (M_dark/M_lens = 0.55), and the total magnification factor of the source galaxy (~ 13). Additionally, we measure an upper-limit constraint on luminous substructure (M_V > 16.2), based on the non-detection of bright satellite galaxies in all data sets. Finally, we utilize the improved image resolution of the AO data to reveal the presence of faint arcs outside of the primary Einstein ring. The positions and orientations of these arcs raise the intriguing possibility that B1938+666 has a second source galaxy, located at a more distant redshift. However, future work is needed to verify this hypothesis.

preprint2012arXiv

The LOFAR radio environment

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

preprint2012arXiv

The SWELLS survey - V. A Salpeter stellar initial mass function in the bulges of massive spiral galaxies

Recent work has suggested that the stellar initial mass function (IMF) is not universal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or morphological type. In this paper, we investigate variations of the IMF within individual galaxies. For this purpose, we use strong lensing and gas kinematics to measure independently the normalisation of the IMF of the bulge and disk components of a sample of 5 massive spiral galaxies with substantial bulge components taken from the SWELLS survey. We find that the stellar mass of the bulges are tightly constrained by the lensing and kinematic data. A comparison with masses based on stellar population synthesis models fitted to optical and near infrared photometry favors a Salpeter-like normalisation of the IMF. Conversely, the disk masses are less well constrained due to degeneracies with the dark matter halo, but are consistent with Milky Way type IMFs in agreement with previous studies. The disks are submaximal at 2.2 disk scale lengths, but due to the contribution of the bulges, the galaxies are baryon dominated at 2.2 disk scale lengths. Globally, our inferred IMF normalisation is consistent with that found for early-type galaxies of comparable stellar mass (> 10^11 M_sun). Our results suggest a non-universal IMF within the different components of spiral galaxies, adding to the well-known differences in stellar populations between disks and bulges.

preprint2012arXiv

Wide-band Simultaneous Observations of Pulsars: Disentangling Dispersion Measure and Profile Variations

Dispersion in the interstellar medium is a well known phenomenon that follows a simple relationship, which has been used to predict the time delay of dispersed radio pulses since the late 1960s. We performed wide-band simultaneous observations of four pulsars with LOFAR (at 40-190 MHz), the 76-m Lovell Telescope (at 1400 MHz) and the Effelsberg 100-m Telescope (at 8000 MHz) to test the accuracy of the dispersion law over a broad frequency range. In this paper we present the results of these observations which show that the dispersion law is accurate to better than 1 part in 100000 across our observing band. We use this fact to constrain some of the properties of the ISM along the line-of-sight and use the lack of any aberration or retardation effects to determine upper limits on emission heights in the pulsar magnetosphere. We also discuss the effect of pulse profile evolution on our observations, and the implications that it could have for precision pulsar timing projects such as the detection of gravitational waves with pulsar timing arrays.

preprint2011arXiv

Galaxy Number Counts and Implications for Strong Lensing

We compare galaxy number counts in HST/ACS fields containing moderate-redshift (0.2<z<1.0) strong gravitational lenses with those in two control samples: (1) the first square degree of the COSMOS survey, comprising 259 ACS fields and (2) 20 "pure parallel" fields randomly located on the sky. Through a Bayesian analysis we determine the expectation values (mu_0) and confidence levels of the underlying number counts for a range of apertures and magnitude bins. Our analysis has produced the following results: (i) We infer that our control samples are not consistent, with the number counts in the COSMOS sample being significantly higher than in the pure parallel sample for 21 <= F814W <= 23. This result matches those found in previous analyses of COSMOS data using different techniques. (ii) We find that small-size apertures, centered on strong lenses, are overdense compared with randomly placed apertures in the control samples, even compared to the COSMOS sample. Correcting for the local clustering of elliptical galaxies, based on the average two-point correlation function reduces this overdensity to the 1-2 sigma level. Thus, the overdensity of galaxies seen along a typical line of sight to a lens can be explained mostly by the natural clustering of galaxies, rather than being due to lenses lying along otherwise biased lines of sight. However, a larger sample of lenses is needed to see if the remaining bias persists when the lens-field uncertainties are smaller. (iii) There is considerable scatter in the lines of sight to _individual_ lens systems, but quantities that are linearly dependent on the external convergence (e.g., H_0) should become unbiased if the extra galaxies that cause the bias can be accounted for in the lens models either through direct modeling or via an informed prior on the external convergence. The number counts can used to set such an informed prior.

preprint2011arXiv

Observing pulsars and fast transients with LOFAR

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

preprint2011arXiv

Radio Interferometric Calibration Using The SAGE Algorithm

The aim of the new generation of radio synthesis arrays such as LOFAR and SKA is to achieve much higher sensitivity, resolution and frequency coverage than what is available now, especially at low frequencies. To accomplish this goal, the accuracy of the calibration techniques used is of considerable importance. Moreover, since these telescopes produce huge amounts of data, speed of convergence of calibration is a major bottleneck. The errors in calibration are due to system noise (sky and instrumental) as well as the estimation errors introduced by the calibration technique itself, which we call solver noise. We define solver noise as the distance between the optimal solution (the true value of the unknowns, uncorrupted by the system noise) and the solution obtained by calibration. We present the Space Alternating Generalized Expectation Maximization (SAGE) calibration technique, which is a modification of the Expectation Maximization algorithm, and compare its performance with the traditional Least Squares calibration based on the level of solver noise introduced by each technique. For this purpose, we develop statistical methods that use the calibrated solutions to estimate the level of solver noise. The SAGE calibration algorithm yields very promising results both in terms of accuracy and speed of convergence. The comparison approaches we adopt introduce a new framework for assessing the performance of different calibration schemes.

preprint2011arXiv

The SWELLS Survey. I. A large spectroscopically selected sample of edge-on late-type lens galaxies

The relative contribution of baryons and dark matter to the inner regions of spiral galaxies provides critical clues to their formation and evolution, but it is generally difficult to determine. For spiral galaxies that are strong gravitational lenses, however, the combination of lensing and kinematic observations can be used to break the disk-halo degeneracy. In turn, such data constrain fundamental parameters such as i) the mass density profile slope and axis ratio of the dark matter halo, and by comparison with dark matter-only numerical simulations the modifications imposed by baryons; ii) the mass in stars and therefore the overall star formation efficiency, and the amount of feedback; iii) by comparison with stellar population synthesis models, the normalization of the stellar initial mass function. In this first paper of a series, we present a sample of 16 secure, 1 probable, and 6 possible strong lensing spiral galaxies, for which multi-band high-resolution images and rotation curves were obtained using the Hubble Space Telescope and Keck-II Telescope as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS). The sample includes 8 newly discovered secure systems. [abridged] We find that the SWELLS sample of secure lenses spans a broad range of morphologies (from lenticular to late-type spiral), spectral types (quantified by Halpha emission), and bulge to total stellar mass ratio (0.22-0.85), while being limited to M_*>10^{10.5} M_sun. The SWELLS sample is thus well-suited for exploring the relationship between dark and luminous matter in a broad range of galaxies. We find that the deflector galaxies obey the same size-mass relation as that of a comparison sample of elongated non-lens galaxies selected from the SDSS survey. We conclude that the SWELLS sample is consistent with being representative of the overall population of high-mass high-inclination disky galaxies.

preprint2011arXiv

The SWELLS survey. II. Breaking the disk-halo degeneracy in the spiral galaxy gravitational lens SDSS J2141-0001

The degeneracy among the disk, bulge and halo contributions to galaxy rotation curves prevents an understanding of the distribution of baryons and dark matter in disk galaxies. In an attempt to break this degeneracy, we present an analysis of the spiral galaxy strong gravitational lens SDSS J2141-0001, discovered as part of the SLACS survey. We present new Hubble Space Telescope multicolor imaging, gas and stellar kinematics data derived from long-slit spectroscopy, and K-band LGS adaptive optics imaging, both from the Keck telescopes. We model the galaxy as a sum of concentric axisymmetric bulge, disk and halo components and infer the contribution of each component, using information from gravitational lensing and gas kinematics. This analysis yields a best-fitting total (disk plus bulge) stellar mass of log_{10}(Mstar/Msun) = 10.99(+0.11,-0.25). The photometric data combined with stellar population synthesis models yield log_{10}(Mstar/Msun) = 10.97\pm0.07, and 11.21\pm0.07 for the Chabrier and Salpeter IMFs, respectively. Accounting for the expected gas fraction of \simeq 20% reduces the lensing plus kinematics stellar mass by 0.10\pm0.05 dex, resulting in a Bayes factor of 11.9 in favor of a Chabrier IMF. The dark matter halo is roughly spherical, with minor to major axis ratio q_{halo}=0.91(+0.15,-0.13). The dark matter halo has a maximum circular velocity of V_{max}=276(+17,-18) km/s, and a central density parameter of log_{10}Δ_{V/2}=5.9(+0.9,-0.5). This is higher than predicted for uncontracted dark matter haloes in LCDM cosmologies, log_{10}Δ_{V/2}=5.2, suggesting that either the halo has contracted in response to galaxy formation, or that the halo has a higher than average concentration. At 2.2 disk scale lengths the dark matter fraction is f_{DM}=0.55(+0.20,-0.15), suggesting that SDSS J2141-0001 is sub-maximal.

preprint2011arXiv

The X-Shooter Lens Survey - I. Dark-Matter Domination and a Salpeter-type IMF in a Massive Early-type Galaxy

We present the first results from the X-shooter Lens Survey (XLENS): an analysis of the massive early-type galaxy SDSS J1148+1930 at redshift z=0.444. We combine its extended kinematic profile -derived from spectra obtained with X-shooter on the ESO VLT- with strong gravitational lensing and multi-color information derived from SDSS images. Our main results are (i) The luminosity-weighted stellar velocity dispersion is σ*(<Reff) = 351\pm10 km/s, extracted from a rectangular aperture of 1.8"x1.6", more accurate and considerably lower than a previously published value of \sim450km/s. (ii) A single-component (stellar plus dark) mass model of the lens galaxy yields a logarithmic total-density slope of γ' = 1.72+0.05-0.06. (iii) The projected stellar mass fraction, derived solely from lensing, is f*(<Reff) = 0.19+0.04-0.09 inside the Einstein radius for a Hernquist profile and no anisotropy. The dark-matter fraction inside the effective radius f_DM(<Reff) = 0.60+0.15\pm0.1. (iv) Based on the SDSS colors, we find f*(<Reff)=0.17 \pm 0.06 for a Salpeter IMF and f*(<Reff)=0.07 \pm 0.02 for a Chabrier IMF. The lensing and dynamics constraints on the stellar mass fraction agree well with those independently derived from the SDSS colors for a Salpeter IMF, preferred over Chabrier IMF at variance with standard results for lower mass galaxies. Dwarf-rich IMFs in the lower mass range of 0.1-0.7 solar mass, with α >= 3 (with dN/dM \propto M^(-α)) -such as that recently suggested for massive ETGs- are excluded at the > 90 % C.L. and in some cases violate the total lensing-derived mass limit. We conclude that this very massive ETG is dark-matter dominated inside one effective radius, consistent with the trend recently found from massive SLACS galaxies, with a IMF normalization consistent with Salpeter.

preprint2010arXiv

A method for efficient measurement of gravitational lens time delays

The Hubble constant value is currently known to 10% accuracy unless assumptions are made for the cosmology (Sandage et al. 2006). Gravitational lens systems provide another probe of the Hubble constant using time delay measurements. However, current investigations of ~20 time delay lenses, albeit of varying levels of sophistication, have resulted in different values of the Hubble constant ranging from 50-80 km/s/Mpc. In order to reduce uncertainties, more time delay measurements are essential together with better determined mass models (Oguri 2007, Saha et al. 2006). We propose a more efficient technique for measuring time delays which does not require regular monitoring with a high-resolution interferometer array. The method uses double image and long-axis quadruple lens systems in which the brighter component varies first and dominates the total flux density. Monitoring the total flux density with low-resolution but high sensitivity radio telescopes provides the variation of the brighter image and is used to trigger high-resolution observations which can then be used to see the variation in the fainter image. We present simulations of this method together with a pilot project using the WSRT (Westerbork Radio Synthesis Telescope) to trigger VLA (Very Large Array) observations. This new method is promising for measuring time delays because it uses relatively small amounts of time on high-resolution telescopes. This will be important because many SKA pathfinder telescopes, such as MeerKAT (Karoo Array Telescope) and ASKAP (Australian Square Kilometre Array Pathfinder), have high sensitivity but limited resolution.

preprint2010arXiv

Dissecting the Gravitational Lens B1608+656. II. Precision Measurements of the Hubble Constant, Spatial Curvature, and the Dark Energy Equation of State

Strong gravitational lens systems with measured time delays between the multiple images provide a method for measuring the "time-delay distance" to the lens, and thus the Hubble constant. We present a Bayesian analysis of the strong gravitational lens system B1608+656, incorporating (i) new, deep Hubble Space Telescope (HST) observations, (ii) a new velocity dispersion measurement of 260+/-15 km/s for the primary lens galaxy, and (iii) an updated study of the lens' environment. When modeling the stellar dynamics of the primary lens galaxy, the lensing effect, and the environment of the lens, we explicitly include the total mass distribution profile logarithmic slope gamma' and the external convergence kappa_ext; we marginalize over these parameters, assigning well-motivated priors for them, and so turn the major systematic errors into statistical ones. The HST images provide one such prior, constraining the lens mass density profile logarithmic slope to be gamma'=2.08+/-0.03; a combination of numerical simulations and photometric observations of the B1608+656 field provides an estimate of the prior for kappa_ext: 0.10 +0.08/-0.05. This latter distribution dominates the final uncertainty on H_0. Compared with previous work on this system, the new data provide an increase in precision of more than a factor of two. In combination with the WMAP 5-year data set, we find that the B1608+656 data set constrains the curvature parameter to be -0.031 < Omega_k < 0.009 (95% CL), a level of precision comparable to that afforded by the current Type Ia SNe sample. Asserting a flat spatial geometry, we find that, in combination with WMAP, H_0 = 69.7 +4.9/-5.0 km/s/Mpc and w=-0.94 +0.17/-0.19 (68% CL), suggesting that the observations of B1608+656 constrain w as tightly as do the current Baryon Acoustic Oscillation data. (abridged)

preprint2010arXiv

Foregrounds for observations of the cosmological 21 cm line: II. Westerbork observations of the fields around 3C196 and the North Celestial Pole

In the coming years a new insight into galaxy formation and the thermal history of the Universe is expected to come from the detection of the highly redshifted cosmological 21 cm line. The cosmological 21 cm line signal is buried under Galactic and extragalactic foregrounds which are likely to be a few orders of magnitude brighter. Strategies and techniques for effective subtraction of these foreground sources require a detailed knowledge of their structure in both intensity and polarization on the relevant angular scales of 1-30 arcmin. We present results from observations conducted with the Westerbork telescope in the 140-160 MHz range with 2 arcmin resolution in two fields located at intermediate Galactic latitude, centred around the bright quasar 3C196 and the North Celestial Pole. They were observed with the purpose of characterizing the foreground properties in sky areas where actual observations of the cosmological 21 cm line could be carried out. The polarization data were analysed through the rotation measure synthesis technique. We have computed total intensity and polarization angular power spectra. Total intensity maps were carefully calibrated, reaching a high dynamic range, 150000:1 in the case of the 3C196 field. [abridged]

preprint2010arXiv

Ionospheric Power-Spectrum Tomography in Radio Interferometry

A tomographic method is described to quantify the three-dimensional power-spectrum of the ionospheric electron-density fluctuations based on radio-interferometric observations by a two-dimensional planar array. The method is valid to first-order Born approximation and might be applicable to correct observed visibilities for phase variations due to the imprint of the full three-dimensional ionosphere. It is shown that not the ionospheric electron density distribution is the primary structure to model in interferometry, but its autocorrelation function or equivalent its power-spectrum. An exact mathematical expression is derived that provides the three dimensional power-spectrum of the ionospheric electron-density fluctuations directly from a rescaled scattered intensity field and an incident intensity field convolved with a complex unit phasor that depends on the w-term and is defined on the full sky pupil plane. In the limit of a small field of view, the method reduces to the single phase screen approximation. Tomographic self-calibration can become important in high-dynamic range observations at low radio frequencies with wide-field antenna interferometers, because a three-dimensional ionosphere causes a spatially varying convolution of the sky, whereas a single phase screen results in a spatially invariant convolution. A thick ionosphere can therefore not be approximated by a single phase screen without introducing errors in the calibration process. By applying a Radon projection and the Fourier projection-slice theorem, it is shown that the phase-screen approach in three dimensions is identical to the tomographic method. Finally we suggest that residual speckle can cause a diffuse intensity halo around sources, due to uncorrectable ionospheric phase fluctuations in the short integrations, which could pose a fundamental limit on the dynamic range in long-integration images.

preprint2010arXiv

The mass distribution of a moderate redshift galaxy group and brightest group galaxy from gravitational lensing and kinematics

The gravitational lens system CLASS B2108+213 has two radio-loud lensed images separated by 4.56 arcsec. The relatively large image separation implies that the lensing is caused by a group of galaxies. In this paper, new optical imaging and spectroscopic data for the lensing galaxies of B2108+213 and the surrounding field galaxies are presented. These data are used to investigate the mass and composition of the lensing structure. The redshift and stellar velocity dispersion of the main lensing galaxy (G1) are found to be z = 0.3648 +/- 0.0002 and sigma_v = 325 +/- 25 km/s, respectively. The optical spectrum of the lensed quasar shows no obvious emission or absorption features and is consistent with a BL Lac type radio source. However, the tentative detection of the G-band and Mg-b absorption lines, and a break in the spectrum of the host galaxy of the lensed quasar gives a likely source redshift of z = 0.67. Spectroscopy of the field around B2108+213 finds 51 galaxies at a similar redshift to G1, thus confirming that there is a much larger structure at z ~ 0.365 associated with this system. The width of the group velocity distribution is 694 +/- 93 km/s, but is non-Gaussian, implying that the structure is not yet viralized. The main lensing galaxy is also the brightest group member and has a surface brightness profile consistent with a typical cD galaxy. A lensing and dynamics analysis of the mass distribution, which also includes the newly found group members, finds that the logarithmic slope of the mass density profile is on average isothermal inside the Einstein radius, but steeper at the location of the Einstein radius. This apparent change in slope can be accounted for if an external convergence gradient, representing the underlying parent halo of the galaxy group, is included in the mass model.

preprint2009arXiv

Detection of a Dark Substructure through Gravitational Imaging

We report the detection of a dark substructure through direct gravitational imaging - undetected in the HST-ACS F814W image - in the gravitational lens galaxy of SLACS SDSSJ0946+1006 (the "Double Einstein Ring"). The detection is based on a Bayesian grid reconstruction of the two-dimensional surface density of the galaxy inside an annulus around its Einstein radius (few kpc). [...] We confirm this detection by modeling the system including a parametric mass model with a tidally truncated pseudo-Jaffe density profile; in that case the substructure mass is M_sub=(3.51+-0.15)x10^9 Msun, located at (-0.651+-0.038,1.040+-0.034)'', precisely where also the surface density map shows a strong convergence peak. [...] We set a lower limit of (M/L)_V}>=120 (Msun/L}_V,sun (3-sigma) inside a sphere of 0.3 kpc centred on the substructure (r_tidal=1.1kpc). The result is robust under substantial changes in the model and the data-set (e.g. PSF, pixel number and scale, source and potential regularization, rotations and galaxy subtraction). Despite being at the limits of detectability, it can therefore not be attributed to obvious systematic effects. Our detection implies a dark matter mass fraction at the radius of the inner Einstein ring of f_CDM=2.15^{+2.05}_{-1.25} percent (68 percent C.L) in the mass range 4x10^6 Msun to 4x10^9 Msun assuming alpha=1.9+-0.1 (with dN/dm ~ m^-alpha). Assuming a flat prior on alpha, between 1.0 and 3.0, increases this to f_CDM=2.56^{+3.26}_{-1.50} percent (68 percent C.L). The likelihood ratio is 0.51 between our best value (f_CDM=0.0215) and that from simulations (f_sim=0.003). Hence the inferred mass fraction, admittedly based on a single lens system, is large but still consistent with predictions. [...]

preprint2009arXiv

Dissecting the Gravitational Lens B1608+656. I. Lens Potential Reconstruction

Strong gravitational lensing is a powerful technique for probing galaxy mass distributions and for measuring cosmological parameters. We present a pixelated approach to modeling simultaneously the lens potential and source intensity of strong gravitational lens systems with extended source-intensity distributions. For systems with sources of sufficient extent such that the separate lensed images are connected by intensity measurements, the accuracy in the reconstructed potential is solely limited by the quality of the data. We apply this potential reconstruction technique to deep HST observations of B1608+656, a four-image gravitational lens system formed by a pair of interacting lens galaxies. We present a comprehensive Bayesian analysis of the system that takes into account the extended source-intensity distribution, dust extinction, and the interacting lens galaxies. Our approach allows us to compare various models of the components of the lens system, which include the point-spread function (PSF), dust, lens galaxy light, source-intensity distribution, and lens potential. Using optimal combinations of the PSF, dust, and lens galaxy light models, we successfully reconstruct both the lens potential and the extended source-intensity distribution of B1608+656. The resulting reconstruction can be used as the basis of a measurement of the Hubble constant. We use our reconstruction of the gravitational potential to study the relative distribution of mass and light in the lensing galaxies. We find that the mass-to-light ratio for the primary lens galaxy is (2.0+/-0.2)h M_{\sun} L_{B,\sun}^{-1} within the Einstein radius 3.9 h^{-1} kpc, in agreement with what is found for noninteracting lens galaxies at the same scales. (Abridged)

preprint2009arXiv

Foregrounds for observations of the cosmological 21 cm line: I. First Westerbork measurements of Galactic emission at 150 MHz in a low latitude field

We present the first results from a series of observations conducted with the Westerbork telescope in the 140--160 MHz range with a 2 arcmin resolution aimed at characterizing the properties of the foregrounds for epoch of reionization experiments. For the first time we have detected fluctuations in the Galactic diffuse emission on scales greater than 13 arcmin at 150 MHz, in the low Galactic latitude area known as Fan region. Those fluctuations have an $rms$ of 14 K. The total intensity power spectrum shows a power--law behaviour down to $\ell \sim 900$ with slope $β^I_\ell = -2.2 \pm 0.3$. The detection of diffuse emission at smaller angular scales is limited by residual point sources. We measured an $rms$ confusion noise of $\sim$3 mJy beam$^{-1}$. Diffuse polarized emission was also detected for the first time at this frequency. The polarized signal shows complex structure both spatially and along the line of sight. The polarization power spectrum shows a power--law behaviour down to $\ell \sim 2700$ with slope $β^P_\ell = -1.65 \pm 0.15$. The $rms$ of polarization fluctuations is 7.2 K on 4 arcmin scales. By extrapolating the measured spectrum of total intensity emission, we find a contamination on the cosmological signal of $δT= \sqrt{\ell (\ell+1) C^I_\ell / 2π} \sim 5.7$ K on 5 arcmin scales and a corresponding $rms$ value of $\sim$18.3 K at the same angular scale. The level of the polarization power spectrum is $δT \sim 3.3$ K on 5 arcmin scales. Given its exceptionally bright polarized signal, the Fan region is likely to represent an upper limit on the sky brightness at moderate and high Galactic latitude.

preprint2009arXiv

Radio counterpart of the lensed sub-mm emission in the cluster MS0451.6-0305: new evidence for the merger scenario

SMM J04542-0301 is an extended (~1 arcmin) sub-mm source located near the core of the cluster MS0451.6-0305. It has been suggested that part of its emission arises from the interaction between a LBG and two EROs at z~2.9 that are multiply-imaged. However, the dramatic resolution difference between the sub-mm map and the optical/NIR images make it difficult to confirm this hypothesis. In this paper, we present a deep (~ 10 microJy/beam), high resolution (~2 arcsec) 1.4 GHz radio map of the cluster core, in which we have identified 6 sources located within SMM J04542-0301. The strong lensing effect in the radio data has been quantified by constructing a new lens model of the cluster. The brightest and most extended of these sources (RJ) is located in the middle of the sub-mm emission, and has no obvious counterpart in the optical/NIR. Three other detections (E1, E2 and E3) seem to be associated with the images of one of the EROs. The last two detections (CR1 and CR2), for which no optical/NIR counterpart have been found, seem to constitute two relatively compact emitting regions embedded in a ~5 arcsec extended radio source located at the position of the sub-mm peak. The presence of this extended component can only be explained if it is being produced by a lensed region of dust obscured star formation in the center of the merger. A comparison between the radio and sub-mm data at the same resolution suggests that E1, E2, E3, CR1 and CR2 are associated with the sub-mm emission. The radio observations provide strong observational evidence in favor of the merger hypothesis. However, the question if RJ is also contributing to the observed sub-mm emission remains open. These results illustrate the promising prospects for radio interferometry and strong gravitational lensing to study the internal structure of SMGs.

preprint2009arXiv

Stars and dark matter in the spiral gravitational lens 2237+0305

We construct a mass model for the spiral lens galaxy 2237+0305, at redshift z_l=0.04, based on gravitational-lensing constraints, HI rotation, and new stellar-kinematic information, based on data taken with the ESI spectrograph on the 10m Keck-II Telescope. High resolution rotation curves and velocity dispersion profiles along two perpendicular directions, close to the major and minor axes of the lens galaxy, were obtained by fitting the Mgb-Fe absorption line region. The stellar rotation curve rises slowly and flattens at r~1.5" (~1.1 kpc). The velocity dispersion profile is approximately flat. A combination of photometric, kinematic and lensing information is used to construct a mass model for the four major mass components of the system -- the dark matter halo, disc, bulge, and bar. The best-fitting solution has a dark matter halo with a logarithmic inner density slope of gamma=0.9+/-0.3 for rho_DM propto r^-gamma, a bulge with M/L_B=6.6+/-0.3 Upsilon_odot, and a disc with M/L_B =1.2+/-0.3 Upsilon_odot, in agreement with measurements of late-type spirals. The bulge dominates support in the inner regions where the multiple images are located and is therefore tightly constrained by the observations. The disc is sub-maximal and contributes 45+/-11 per cent of the rotational support of the galaxy at 2.2r_d. The halo mass is (2.0+/-0.6) x 10^12 M_odot, and the stellar to virial mass ratio is 7.0+/-2.3 per cent, consistent with typical galaxies of the same mass.

preprint2009arXiv

Statistics of mass substructure from strong gravitational lensing: quantifying the mass fraction and mass function

A Bayesian statistical formalism is developed to quantify the level at which the mass function slope (alpha) and the projected cumulative mass fraction (f) of (CDM) substructure in strong gravitational-lens galaxies, with arcs or Einstein rings, can be recovered as function of the lens-survey parameters and the detection threshold of the substructure mass. The method is applied to different sets of mock data to explore a range of observational limits: (i) the number of lens galaxies in the survey, (ii) the mass threshold, Mlow, for the detection of substructures and (iii) the uncertainty of the measured substructure masses. We explore two different priors on the mass function slope: a uniform prior and a Gaussian prior with alpha = 1.90+-0.1. With a substructure detection threshold Mlow=3x10^8 Msun, the number of lenses available now (n_l=30), a true dark-matter mass fraction in (CDM) substructure <=1.0% and a prior of alpha = 1.90+-0.1, we find that the upper limit of f can be constrained down to a level <=1.0% (95% CL). In the case of a Gaussian prior on alpha, it is always possible to set stringent constraints on both parameters. We also find that lowering the detection threshold has the largest impact on the ability to recover alpha, because of the (expected) steep mass-function slope. In the future, thanks to new surveys with telescopes, such as SKA, LSST and JDEM and follow-up telescopes with high-fidelity data, a significant increase in the number of known lenses will allow us to recover the satellite population in its completeness. For example, a sample of 200 lenses, equivalent in data-quality to the Sloan Lens ACS Survey and a detection threshold of 10^8 Msun, allows one to determine f=0.5+-0.1% (68% CL) and alpha=1.90+-0.2 (68% CL).

preprint2009arXiv

The initial mass function of early-type galaxies

We determine an absolute calibration of the initial mass function (IMF) of early-type galaxies, by studying a sample of 56 gravitational lenses identified by the SLACS Survey. Under the assumption of standard Navarro, Frenk & White dark matter halos, a combination of lensing, dynamical, and stellar population synthesis models is used to disentangle the stellar and dark matter contribution for each lens. We define an "IMF mismatch" parameter α=M*(L+D)/M*(SPS) as the ratio of stellar mass inferred by a joint lensing and dynamical models (M*(L+D)) to the current stellar mass inferred from stellar populations synthesis models (M*(SPS)). We find that a Salpeter IMF provides stellar masses in agreement with those inferred by lensing and dynamical models (<\log α>=0.00+-0.03+-0.02), while a Chabrier IMF underestimates them (<\log α>=0.25+-0.03+-0.02). A tentative trend is found, in the sense that αappears to increase with galaxy velocity dispersion. Taken at face value, this result would imply a non universal IMF, perhaps dependent on metallicity, age, or abundance ratios of the stellar populations. Alternatively, the observed trend may imply non-universal dark matter halos with inner density slope increasing with velocity dispersion. While the degeneracy between the two interpretations cannot be broken without additional information, the data imply that massive early-type galaxies cannot have both a universal IMF and universal dark matter halos.

preprint2009arXiv

The Sloan Lens ACS Survey. IX. Colors, Lensing and Stellar Masses of Early-type Galaxies

We present the current photometric dataset for the Sloan Lens ACS (SLACS) Survey, including HST photometry from ACS, WFPC2, and NICMOS. These data have enabled the confirmation of an additional 15 grade `A' (certain) lens systems, bringing the number of SLACS grade `A' lenses to 85; including 13 grade `B' (likely) systems, SLACS has identified nearly 100 lenses and lens candidates. Approximately 80% of the grade `A' systems have elliptical morphologies while ~10% show spiral structure; the remaining lenses have lenticular morphologies. Spectroscopic redshifts for the lens and source are available for every system, making SLACS the largest homogeneous dataset of galaxy-scale lenses to date. We have developed a novel Bayesian stellar population analysis code to determine robust stellar masses with accurate error estimates. We apply this code to deep, high-resolution HST imaging and determine stellar masses with typical statistical errors of 0.1 dex; we find that these stellar masses are unbiased compared to estimates obtained using SDSS photometry, provided that informative priors are used. The stellar masses range from 10^10.5 to 10^11.8 M$_\odot$ and the typical stellar mass fraction within the Einstein radius is 0.4, assuming a Chabrier IMF. The ensemble properties of the SLACS lens galaxies, e.g. stellar masses and projected ellipticities, appear to be indistinguishable from other SDSS galaxies with similar stellar velocity dispersions. This further supports that SLACS lenses are representative of the overall population of massive early-type galaxies with M* >~ 10^11 M$_\odot$, and are therefore an ideal dataset to investigate the kpc-scale distribution of luminous and dark matter in galaxies out to z ~ 0.5.

preprint2008arXiv

The Sloan Lens ACS Survey. VIII. The relation between environment and internal structure of early-type galaxies

We study the relation between the internal structure of early-type galaxies and their environment using 70 strong gravitational lenses from the Sloan ACS Lens Survey. The Sloan database is used to determine two measures of overdensity of galaxies around each lens: the projected number density of galaxies inside the tenth nearest neighbor (Σ_{10}) and within a cone of radius one h^{-1} Mpc (D_1). Our main results are: 1) The average overdensity is somewhat larger than unity, consistent with lenses preferring overdense environments as expected for massive early-type galaxies (12/70 lenses are in known groups/clusters). 2) The distribution of overdensities is indistinguishable from that of "twin" non-lens galaxies selected from SDSS to have the same redshift and stellar velocity dispersion σ_*. Thus, within our errors, lens galaxies are an unbiased population, and the SLACS results can be generalized to the overall population of early-type galaxies. 3) Typical contributions from external mass distribution are no more than a few per cent, reaching 10-20% (~0.05-0.10 external convergence) only in the most extreme overdensities. 4) No significant correlation between overdensity and slope of the mass density profile of the lens is found. 5) Satellite galaxies (those with a more luminous companion) have marginally steeper mass density profiles than central galaxies (those without). This result suggests that tidal stripping may affect the mass structure of early-type galaxies down to kpc scales probed by strong lensing, when they fall into larger structures [ABRIDGED].

preprint2006arXiv

Three Gravitational Lenses for the Price of One: Enhanced Strong Lensing through Galaxy Clustering

We report the serendipitous discovery of two strong gravitational lens candidates (ACS J160919+6532 and ACS J160910+6532) in deep images obtained with the Advanced Camera for Surveys on the Hubble Space Telescope, each less than 40 arcsec from the previously known gravitational lens system CLASS B1608+656. The redshifts of both lens galaxies have been measured with Keck and Gemini: one is a member of a small galaxy group at z~0.63, which also includes the lensing galaxy in the B1608+656 system, and the second is a member of a foreground group at z~0.43. By measuring the effective radii and surface brightnesses of the two lens galaxies, we infer their velocity dispersions based on the passively evolving Fundamental Plane (FP) relation. Elliptical isothermal lens mass models are able to explain their image configurations within the lens hypothesis, with a velocity dispersion compatible with that estimated from the FP for a reasonable source-redshift range. Based on the large number of massive early-type galaxies in the field and the number-density of faint blue galaxies, the presence of two additional lens systems around CLASS B1608+656 is not unlikely in hindsight. Gravitational lens galaxies are predominantly early-type galaxies, which are clustered, and the lensed quasar host galaxies are also clustered. Therefore, obtaining deep high-resolution images of the fields around known strong lens systems is an excellent method of enhancing the probability of finding additional strong gravitational lens systems.

preprint2005arXiv

Gravitational Imaging of CDM Substructure

We propose the novel method of ``gravitational imaging'' to detect and quantify luminous and dark-matter substructure in gravitational-lens galaxies. The method utilizes highly-magnified Einstein rings and arcs as sensitive probes of small perturbations in the lens potential (due to the presence of mass substructure), reconstructing the gravitational lens potential non-parametrically. Numerical simulations show that the implemented algorithm can reconstruct the smooth mass distribution of a typical lens galaxy - exhibiting reasonable signal-to-noise Einstein rings - as well as compact substructure with masses as low as M_sub~10^-3 M_lens, if present. ``Gravitational imaging'' of pure dark-matter substructure around massive galaxies can provide a new window on the standard cold-dark-matter paradigm, using very different physics than ground-based direct-detection experiments, and probe the hierarchical structure-formation model which predicts this substructure to exist in great abundance.

preprint2002arXiv

The Structure and Dynamics of Luminous and Dark Matter in the Early-Type Lens Galaxy of 0047-281 at z=0.485

We have measured the kinematic profile of the early-type (E/S0) lens galaxy in the system 0047-281 (z=0.485) with the {\sl Echelle Spectrograph and Imager} (ESI) on the W.M. Keck--II Telescope, as part of the {\sl Lenses Structure and Dynamics (LSD) Survey}. The central velocity dispersion is σ=229\pm 15 \kms, and the dispersion profile is nearly flat to beyond one effective radius (R_e). Surface photometry of the lens galaxy is measured from {\sl Hubble Space Telescope} images. From the offset from the local Fundamental Plane (FP), we measure an evolution of the effective mass-to-light ratio of Δ\log M/L_B=-0.37\pm0.06 between z=0 and z=0.485, consistent with the observed evolution of field E/S0 galaxies. (We assume h_{65}=1,Ω_{\rm m}=0.3 and Ω_Λ=0.7 throughout.) Gravitational lens models provide a mass of M_{\rmE}=(4.06\pm0.20)\times 10^{11} h_{65}^{-1}M_\odot inside the Einstein radius of R_{\rm E}=(8.70\pm0.07) h_{65}^{-1} kpc. This allows us to break the degeneracy between velocity anisotropy and density profile, typical of dynamical models for E/S0 galaxies. We find that constant M/L model, even with strongly tangential anisotropy of the stellar velocity ellipsoid, are excluded at >99.9%CL. The total mass distribution inside R_{\rm E} can be described by a single power-law density profile, ρ_t\propto r^{-γ'}, with an effective slope γ'=1.90^{+0.05}_{-0.23} (68%CL; \pm0.1 systematic error). Two-component models yield an upper limit (68% CL) of γ\le 1.55(1.12) on the power-law slope of the dark-matter density profile and a projected dark-matter mass fraction of 0.41(0.54)^{+0.15}_{-0.05}(^{+0.09}_{-0.06}) (68% CL) inside R_{\rm E}, for Osipkov--Merritt models with anisotropy radius r_i=\infty(R_e).

preprint1999arXiv

B0218+357: Time Delays and New MERLIN/VLA 5GHz Maps of the Einstein Ring

This poster presents a new 5 GHz combined MERLIN/VLA map of B0218+357 which shows for the first time believable substructure in the Einstein ring. This will now be exploited for further constraints on the model which presently dominates the error on the estimate of H_0 derived from the time delay (10.5+/-0.4 days) measured for this system.

preprint1999arXiv

Measuring Cosmological Parameters with the JVAS and CLASS Gravitational Lens Surveys

The JVAS (Jodrell Bank-VLA Astrometric Survey) and CLASS (Cosmic Lens All-Sky Survey) are well-defined surveys containing about ten thousand flat-spectrum radio sources. For many reasons, flat-spectrum radio sources are particularly well-suited as a population from which one can obtain unbiased samples of gravitational lenses. These are by far the largest gravitational (macro)lens surveys, and particular attention was paid to constructing a cleanly-defined sample for the survey itself and for the underlying luminosity function. Here we present the constraints on cosmological parameters, particularly the cosmological constant, derived from JVAS and combine them with constraints from optical gravitational lens surveys, `direct' measurements of $Ω_{0}$, $H_{0}$ and the age of the universe, and constraints derived from CMB anisotropies, before putting this final result into the context of the latest results from other, independent cosmological tests.

L. V. E. Koopmans

What is connected

Connect this record

See the researcher in context

Building this map preview

84 published item(s)

The impact of diffuse Galactic emission on direction-independent gain calibration in high-redshift 21 cm observations

Degree-Scale Galactic Radio Emission at 122 MHz around the North Celestial Pole with LOFAR-AARTFAAC

High-quality strong lens candidates in the final Kilo Degree survey footprint

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

Sub-arcsecond imaging with the International LOFAR Telescope: II. Completion of the LOFAR Long-Baseline Calibrator Survey

Cassiopeia A, Cygnus A, Taurus A, and Virgo A at ultra-low radio frequencies

Constraining the intergalactic medium at $z\approx$ 9.1 using LOFAR Epoch of Reionization observations

LOFAR 144-MHz follow-up observations of GW170817

New high-quality strong lens candidates with deep learning in the Kilo Degree Survey

SHARP -- VI. Evidence for CO (1-0) molecular gas extended on kpc-scales in AGN star forming galaxies at high redshift

TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography

The e-MERLIN Galaxy Evolution Survey (e-MERGE): Overview and Survey Description

A hierarchical Bayesian approach for reconstructing the Initial Mass Function of Single Stellar Populations

Polarization leakage in epoch of reionization windows: II. Primary beam model and direction dependent calibration

Probing Ionospheric Structures using the LOFAR radio telescope

SHARP - II. Mass structure in strong lenses is not necessarily dark matter substructure: A flux ratio anomaly from an edge-on disc

The Herschel-ATLAS: a sample of 500μm-selected lensed galaxies over 600 square degrees

A Bayesian analysis of redshifted 21-cm HI signal and foregrounds: Simulations for LOFAR

Are the total mass density and the low-mass end slope of the IMF anti-correlated?

Galactic interstellar filaments as probed by LOFAR and Planck

Linear polarization structures in LOFAR observations of the interstellar medium in the 3C196 field

Polarization leakage in Epoch of Reionization windows: I. LOFAR observations of the 3C196 field

Scintillation noise in widefield radio interferometry

Scintillation noise power spectrum and its impact on high redshift 21-cm observations

Simulating the 21cm forest detectable with LOFAR and SKA in the spectra of high-z GRBs

Strong gravitational lensing with the SKA

Strong lens search in the ESO public Survey KiDS

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

The X-shooter Lens Survey - II. Sample presentation and spatially resolved kinematics

Cosmology from gravitational lens time delays and Planck data

Galaxy Masses

Inference of the Cold Dark Matter substructure mass function at z=0.2 using strong gravitational lenses

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

Lunar occultation of the diffuse radio sky: LOFAR measurements between 35 and 80 MHz

Measuring gravitational lens time delays using low-resolution radio monitoring observations

Radio Monitoring Campaigns of Six Strongly Lensed Quasars

Calibrating High-Precision Faraday Rotation Measurements for LOFAR and the Next Generation of Low-Frequency Radio Telescopes

Chromatic effects in the 21 cm global signal from the cosmic dawn

Dark energy with gravitational lens time delays

Initial deep LOFAR observations of Epoch of Reionization windows: I. The North Celestial Pole

LEDDB: LOFAR Epoch of Reionization Diagnostic Database

LOFAR insights into the epoch of reionization from the cross power spectrum of 21cm emission and galaxies

LOFAR: The LOw-Frequency ARray

Probing reionization with LOFAR using 21-cm redshift space distortions

Studying Galactic interstellar turbulence through fluctuations in synchrotron emission: First LOFAR Galactic foreground detection

The brightness and spatial distributions of terrestrial radio sources

The Scale of the Problem : Recovering Images of Reionization with GMCA

Two accurate time-delay distances from strong lensing: Implications for cosmology

Disentangling Baryons and Dark Matter in the Spiral Gravitational Lens B1933+503

Evidence For A Mild Steepening And Bottom-Heavy IMF In Massive Galaxies From Sodium And Titanium-Oxide Indicators

First LOFAR observations at very low frequencies of cluster-scale non-thermal emission: the case of Abell 2256

Foreground Removal using FastICA: A Showcase of LOFAR-EoR

Gravitational detection of a low-mass dark satellite at cosmological distance

Imaging neutral hydrogen on large-scales during the Epoch of Reionization with LOFAR

Prospects for detecting the 21cm forest from the diffuse intergalactic medium with LOFAR

SHARP - I. A high-resolution multi-band view of the infra-red Einstein ring of JVAS B1938+666

The LOFAR radio environment

The SWELLS survey - V. A Salpeter stellar initial mass function in the bulges of massive spiral galaxies

Wide-band Simultaneous Observations of Pulsars: Disentangling Dispersion Measure and Profile Variations

Galaxy Number Counts and Implications for Strong Lensing

Observing pulsars and fast transients with LOFAR

Radio Interferometric Calibration Using The SAGE Algorithm

The SWELLS Survey. I. A large spectroscopically selected sample of edge-on late-type lens galaxies

The SWELLS survey. II. Breaking the disk-halo degeneracy in the spiral galaxy gravitational lens SDSS J2141-0001

The X-Shooter Lens Survey - I. Dark-Matter Domination and a Salpeter-type IMF in a Massive Early-type Galaxy

A method for efficient measurement of gravitational lens time delays

Dissecting the Gravitational Lens B1608+656. II. Precision Measurements of the Hubble Constant, Spatial Curvature, and the Dark Energy Equation of State

Foregrounds for observations of the cosmological 21 cm line: II. Westerbork observations of the fields around 3C196 and the North Celestial Pole

Ionospheric Power-Spectrum Tomography in Radio Interferometry

The mass distribution of a moderate redshift galaxy group and brightest group galaxy from gravitational lensing and kinematics

Detection of a Dark Substructure through Gravitational Imaging

Dissecting the Gravitational Lens B1608+656. I. Lens Potential Reconstruction

Foregrounds for observations of the cosmological 21 cm line: I. First Westerbork measurements of Galactic emission at 150 MHz in a low latitude field

Radio counterpart of the lensed sub-mm emission in the cluster MS0451.6-0305: new evidence for the merger scenario