Abell 2163 at $z \simeq 0.201$ is one of the most massive galaxy clusters known, very likely in a post-merging phase. Data from several observational windows suggest a complex mass structure with interacting subsystems, which makes the reconstruction of a realistic merging scenario very difficult. A missing key element in this sense is unveiling the cluster mass distribution at high resolution. We perform such a reconstruction of the cluster inner total mass through a strong lensing model based on new spectroscopic redshift measurements. We use data from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to confirm 12 multiple images of 4 sources with redshift values from 1.16 to 2.72. We also discover four new multiple images and identify 29 cluster members and 35 foreground and background sources. The resulting galaxy member and image catalogs are used to build five cluster total mass models. The fiducial model consists of 111 small-scale subhalos plus a diffuse component, which is centered $\sim2$ arcseconds away from the BCG belonging to the east Abell 2163 subcluster. We confirm that the latter is well represented by a single, large-scale mass component. Its strong elongation towards a second (west) subcluster confirms the existence of a preferential axis, corresponding to the merging direction. From the fiducial model, we extrapolate the cumulative projected total mass profile and measure a value of $M(<300\,$kpc$) = 1.43^{+0.07}_{-0.06}\times 10^{14}\,$M$_{\odot}$, which has a significantly reduced statistical error compared with previous estimates, thanks to the inclusion of the spectroscopic redshifts. Our strong lensing results are very accurate: the model-predicted positions of the multiple images are, on average, only $0.15$ arcseconds away from the observed ones.

We carry out a detailed strong lensing analysis of a sub-sample of eight galaxy clusters of the CLASH survey, in the redshift range of $ z_{\rm cluster} = [0.23-0.59]$, using extensive spectroscopic information, primarily MUSE archival data complemented with CLASH-VLT redshift measurements. Different models are tested in each cluster depending on the complexity of its mass distribution and on the number of detected multiple images. Four clusters show more than five spectroscopically confirmed multiple image families. In this sample, we do not make use of families that are only photometrically identified, in order to reduce model degeneracies and systematics due to the potential misidentifications of some multiple images. We present spectroscopic confirmation of 27 multiply lensed sources, with no previous spectroscopic measurements, spanning over the redshift range of $z_{\rm src}=[0.7-6.1]$. Moreover, we confirm an average of $48$ galaxy members in the core of each cluster, thanks to the high efficiency and large field of view of MUSE. Despite having different properties (i.e., number of mass components, total mass, redshift, etc), the projected total mass and mass density profiles of all clusters have very similar shapes, when rescaled by independent measurements of $M_{200c}$ and $R_{200c}$. Specifically, we measure the mean value of the projected total mass of our cluster sample within 10 (20)% of $R_{200c}$ to be 0.13 (0.32) of $M_{200c}$, with a remarkably small scatter of 5 (6)%. Furthermore, the large number of high-z sources and the precise magnification maps derived in this work for four clusters add up to the sample of high-quality gravitational telescopes to be used to study the faint and distant Universe. The strong lensing models and the full redshift catalogues from MUSE are publicly available.

We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We use VLT/MUSE integral-field spectroscopy in the cluster cores to measure the stellar velocity dispersion, $\sigma$, of 40-60 member galaxies per cluster, covering 4-5 magnitudes to $m_{F160W}\simeq 21.5$. We verify the robustness and quantify the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determine the normalization and slope of the galaxy $L\mbox{-}\sigma$ Faber-Jackson relation in each cluster and use these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, $\sigma_0$, and truncation radius, $r_{cut}$, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the $\sigma_0\mbox{-}r_{cut}$ scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of $\sigma_0$, or halo masses. We derive galaxy central velocity dispersion functions of the three clusters and found them in agreement with each other. Sub-halo mass functions determined with this method can be compared with those obtained from cosmological simulations.

We present a comparative study of the physical properties and the spatial distribution of column density peaks in two Giant Molecular Clouds (GMC), the Pipe Nebula and Orion A, which exemplify opposite cases of star cluster formation stages. The density peaks were extracted from dust extinction maps constructed from Herschel/SPIRE farinfrared images. We compare the distribution functions for dust temperature, mass, equivalent radius and mean volume density of peaks in both clouds, and made a more fair comparison by isolating the less active Tail region in Orion A and by convolving the Pipe Nebula map to simulate placing it at a distance similar to that of the Orion Complex. The peak mass distributions for Orion A, the Tail, and the convolved Pipe, have similar ranges, sharing a maximum near 5 M$_\odot$, and a similar power law drop above 10 M$_\odot$. Despite the clearly distinct evolutive stage of the clouds, there are very important similarities in the physical and spatial distribution properties of the column density peaks, pointing to a scenario where they form as a result of uniform fragmentation of filamentary structures across the various scales of the cloud, with density being the parameter leading the fragmentation, and with clustering being a direct result of thermal fragmentation at different spatial scales. Our work strongly supports the idea that the formation of clusters in GMC could be the result of the primordial organization of pre-stellar material

In this paper we present XNICER, an optimized multi-band extinction technique based on the extreme deconvolution of the intrinsic colors of objects observed through a molecular cloud. XNICER follows a rigorous statistical approach and provides the full Bayesian inference of the extinction for each observed object. Photometric errors in both the training control field and in the science field are properly taken into account. XNICER improves over the known extinction methods and is computationally fast enough to be used on large datasets of objects. Our tests and simulations show that this method is able to reduce the noise associated with extinction measurements by a factor 2 with respect to the previous NICER algorithm, and it has no evident bias even at high extinctions.

We have extended and refined the existing young stellar object (YSO) catalogs for the Orion A molecular cloud, the closest massive star-forming region to Earth. This updated catalog is driven by the large spatial coverage (18.3 deg$^2$, $\sim$950 pc$^2$), seeing limited resolution ($\sim$0.7$"$), and sensitivity ($K_s<19$ mag) of the ESO-VISTA near-infrared survey of the Orion A cloud (VISION). Combined with archival mid- to far-infrared data, the VISTA data allow for a refined and more robust source selection. We estimate that among previously known protostars and pre-main-sequence stars with disks, source contamination levels (false positives) are at least $\sim$6.4% and $\sim$2.3%, respectively, mostly due to background galaxies and nebulosities. We identify 274 new YSO candidates using VISTA/Spitzer based selections within previously analyzed regions, and VISTA/WISE based selections to add sources in the surroundings, beyond previously analyzed regions. The WISE selection method recovers about 59% of the known YSOs in Orion A's low-mass star-forming part L1641, which shows what can be achieved by the all-sky WISE survey in combination with deep near-infrared data in regions without the influence of massive stars. The new catalog contains 2980 YSOs, which were classified based on the de-reddened mid-infrared spectral index into 188 protostars, 185 flat-spectrum sources, and 2607 pre-main-sequence stars with circumstellar disks. We find a statistically significant difference in the spatial distribution of the three evolutionary classes with respect to regions of high dust column-density, confirming that flat-spectrum sources are at a younger evolutionary phase compared to Class IIs, and are not a sub-sample seen at particular viewing angles.

We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhanced star-forming (bent) Head, and a lower density and star-formation quieter $\sim$75 pc long Tail. The true extent of Orion A is not the projected $\sim$40 pc but $\sim$90 pc, making it by far the largest molecular cloud in the local neighborhood. Its aspect ratio ($\sim$30:1) and high column-density fraction ($\sim45\%$) make it similar to large-scale Milky Way filaments ("bones"), despite its distance to the galactic mid-plane being an order of magnitude larger than typically found for these structures.

We present a multi-wavelength study of the massive ($M_{200\textrm{c}} \approx 1$-$2 \times 10^{15} M_\odot$) galaxy clusters RXC J2248.7$-$4431, MACS J0416.1$-$2403, and MACS J1206.2$-$0847 at $z \approx 0.4$. Using the X-ray surface brightness of the clusters from deep Chandra data to model their hot gas, we are able to disentangle this mass term from the diffuse dark matter in our new strong-lensing analysis, with approximately $50$-$100$ secure multiple images per cluster, effectively separating the collisional and collisionless mass components of the clusters. At a radial distance of $10\%$ of $R_{200\textrm{c}}$ (approximately $200$ kpc), we measure a projected total mass of $(0.129 \pm 0.001)$, $(0.131 \pm 0.001)$ and $(0.137 \pm 0.001)\times M_{200\textrm{c}}$, for RXC J2248, MACS J0416 and MACS J1206, respectively. These values are surprisingly similar, considering the large differences in the merging configurations, and, as a consequence, in the mass models of the clusters. Interestingly, at the same radii, the hot gas over total mass fractions differ substantially, ranging from $0.082 \pm 0.001$ to $0.133 \pm 0.001$, reflecting the various dynamical states of the clusters. Moreover, we do not find a statistically significant offset between the positions of the peak of the diffuse dark matter component and of the BCG in the more complex clusters of the sample. We extend to this sample of clusters previous findings of a number of massive sub-halos higher than in numerical simulations. These results highlight the importance of a proper separation of the different mass components to study in detail the properties of dark matter in galaxy clusters.

Multi-wavelength observations in the sub-mm regime provide information on the distribution of both the dust column density and the effective dust temperature in molecular clouds. In this study, we created high-resolution and high-dynamic-range maps of the Pipe nebula region and explored the value of dust-temperature measurements in particular towards the dense cores embedded in the cloud. The maps are based on data from the Herschel and Planck satellites, and calibrated with a near-infrared extinction map based on 2MASS observations. We have considered a sample of previously defined cores and found that the majority of core regions contain at least one local temperature minimum. Moreover, we observed an anti-correlation between column density and temperature. The slope of this anti-correlation is dependent on the region boundaries and can be used as a metric to distinguish dense from diffuse areas in the cloud if systematic effects are addressed appropriately. Employing dust-temperature data thus allows us to draw conclusions on the thermodynamically dominant processes in this sample of cores: external heating by the interstellar radiation field and shielding by the surrounding medium. In addition, we have taken a first step towards a physically motivated core definition by recognising that the column-density-temperature anti-correlation is sensitive to the core boundaries. Dust-temperature maps therefore clearly contain valuable information about the physical state of the observed medium.

Realizing Refsdal's original idea from 1964, we present estimates of the Hubble constant that are complementary to and potentially competitive with those of other cosmological probes. We use the observed positions of 89 multiple images, with extensive spectroscopic information, from 28 background sources and the measured time delays between the images S1-S4 and SX of supernova "Refsdal" ($z = 1.489$), which were obtained thanks to Hubble Space Telescope (HST) deep imaging and Multi Unit Spectroscopic Explorer (MUSE) data. We extend the strong lensing modeling of the Hubble Frontier Fields (HFF) galaxy cluster MACS J1149.5$+$2223 ($z = 0.542$), published by Grillo et al. (2016), and explore different $\Lambda$CDM models. Taking advantage of the lensing information associated to the presence of very close pairs of multiple images at various redshifts and to the extended surface brightness distribution of the SN Refsdal host, we can reconstruct the total mass density profile of the cluster very precisely. The combined dependence of the multiple image positions and time delays on the cosmological parameters allows us to infer the values of $H_{0}$ and $\Omega_{\rm m}$ with relative (1$\sigma$) statistical errors of, respectively, 6% (7%) and 31% (26%) in flat (general) cosmological models, assuming a conservative 3% uncertainty on the final time delay of image SX and, remarkably, no priors from other cosmological experiments. Our best estimate of $H_{0}$, based on the model described in this work, will be presented when the final time-delay measurement becomes available. Our results show that it is possible to utilize time delays in lens galaxy clusters as an important alternative tool for measuring the expansion rate and the geometry of the Universe.

We analyse the Kormendy relations (KRs) of the two Frontier Fields clusters, Abell S1063, at z = 0.348, and MACS J1149.5+2223, at z = 0.542, exploiting very deep Hubble Space Telescope photometry and VLT/MUSE integral field spectroscopy. With this novel dataset, we are able to investigate how the KR parameters depend on the cluster galaxy sample selection and how this affects studies of galaxy evolution based on the KR. We define and compare four different galaxy samples according to: (a) S\'ersic indices: early-type ('ETG'), (b) visual inspection: 'ellipticals', (c) colours: 'red', (d) spectral properties: 'passive'. The classification is performed for a complete sample of galaxies with m$_{\textit{F814W}} \le$ 22.5 ABmag (M$_{*}$ $\gtrsim 10^{10.0}$ M$_{\odot}$). To derive robust galaxy structural parameters, we use two methods: (1) an iterative estimate of structural parameters using images of increasing size, in order to deal with closely separated galaxies and (2) different background estimations, to deal with the Intracluster light contamination. The comparison between the KRs obtained from the different samples suggests that the sample selection could affect the estimate of the best-fitting KR parameters. The KR built with ETGs is fully consistent with the one obtained for ellipticals and passive. On the other hand, the KR slope built on the red sample is only marginally consistent with those obtained with the other samples. We also release the photometric catalogue with structural parameters for the galaxies included in the present analysis.

We have investigated the shape of the extinction curve in the infrared up to ~25 {\mu}m for the Orion A star-forming complex. The basis of this work is near-infrared data acquired with VISTA, in combination with Pan-STARRS and mid-infrared Spitzer photometry. We obtain colour excess ratios for eight passbands by fitting a series of colour-colour diagrams. The fits are performed using Markov chain Monte Carlo methods, together with a linear model under a Bayesian formalism. The resulting colour excess ratios are directly interpreted as a measure of the extinction law. We show that the Orion A molecular cloud is characterized by flat mid-infrared extinction, similar to many other recently studied sightlines. Moreover, we find statistically significant evidence that the extinction law from ~1 {\mu}m to at least ~6 {\mu}m varies across the cloud. In particular, we find a gradient along galactic longitude, where regions near the ONC show a different extinction law compared to L1641 and L1647, the low-mass star-forming sites in the cloud complex. These variations are of the order of only 3% and are most likely caused by the influence of the massive stars on their surrounding medium. While the observed general trends in our measurements are in agreement with model predictions, both well-established and new dust grain models are not able to fully reproduce our infrared extinction curve. We also present a new extinction map featuring a resolution of 1 arcmin and revisit the correlation between extinction and dust optical depth. This analysis shows that cloud substructure, which is not sampled by background sources, affects the conversion factor between these two measures. In conclusion, we argue that specific characteristics of the infrared extinction law are still not well understood, but Orion A can serve as an unbiased template for future studies.

We use Multi Unit Spectroscopic Explorer (MUSE) observations of the galaxy cluster MACS J1149.5+2223 to explore the kinematics of the grand-design spiral galaxy Sp1149 hosting the SN Refsdal. Sp1149 lies at $z\simeq1.49$, has a stellar mass $M_*\simeq5\times10^9 \, \mathrm{M_\odot}$, a star-formation rate $\mathrm{SFR} \simeq1-6 \, \mathrm{M_\odot/yr}$ and represents a likely progenitor of a Milky-Way-like galaxy. All the four multiple images of Sp1149 in our data show strong OII-line emissions pointing to a clear rotation pattern. We take advantage of the gravitational lensing magnification effect ($\simeq 4 \times$) on the OII emission of the least distorted image to fit 3D kinematic models to the MUSE data-cube and derive the rotation curve and the velocity dispersion profile of Sp1149. We find that the rotation curve steeply rises, peaks at $R\simeq1$ kpc and then (initially) declines and flattens to an average $V_\mathrm{flat} = 128^{+29}_{-19}$ km/s. The shape of the rotation curve is well determined but the actual value of $V_\mathrm{flat}$ is quite uncertain because of the nearly face-on configuration of the galaxy. The intrinsic velocity dispersion due to gas turbulence is almost constant across the entire disc with an average of $27\pm5$ km/s. This value is consistent with $z=0$ measurements in the ionized gas component and a factor 2-4 lower than other estimates in different galaxies at similar redshifts. The average stellar-to-total mass fraction is of the order of one fifth. Our kinematic analysis returns the picture of a regular star-forming, mildly turbulent, rotation-dominated ($V / \sigma\simeq5$) spiral galaxy in a 4 Gyr old Universe.

We present a high resolution dissection of the two-dimensional total mass distribution in the core of the Hubble Frontier Fields galaxy cluster MACS J0416.1-2403, at z ~ 0.396. We exploit HST/WFC3 near-IR (F160W) imaging, VLT/MUSE spectroscopy, and Chandra data to separate the stellar, hot gas, and dark-matter mass components in the inner 300 kpc of the cluster. We combine the recent results of our refined strong lensing analysis, which includes the contribution of the intracluster gas, with the modeling of the surface brightness and stellar mass distributions of 193 cluster members, of which 144 are spectroscopically confirmed. We find that moving from 10 to 300 kpc from the cluster center the stellar to total mass fraction decreases from 12% to 1% and the hot gas to total mass fraction increases from 3% to 9%, resulting in a baryon fraction of approximately 10% at the outermost radius. We measure that the stellar component represents ~ 30%, near the cluster center, and 15%, at larger clustercentric distances, of the total mass in the cluster substructures. We subtract the baryonic mass component from the total mass distribution and conclude that within 30 kpc (~ 3 times the effective radius of the BCG) from the cluster center the surface mass density profile of the total mass and global (cluster plus substructures) dark-matter are steeper and that of the diffuse (cluster) dark-matter is shallower than a NFW profile. Our current analysis does not point to a significant offset between the cluster stellar and dark-matter components. This detailed and robust reconstruction of the inner dark-matter distribution in a larger sample of galaxy clusters will set a new benchmark for different structure formation scenarios.

We present a new strong lensing analysis of the galaxy cluster MACS J1206.2-0847 (MACS 1206), at z=0.44, using deep spectroscopy from CLASH-VLT and VLT/MUSE archival data in combination with imaging from the Cluster Lensing and Supernova survey with Hubble. MUSE observations enable the spectroscopic identification of 23 new multiply imaged sources, extending the previous compilations by a factor of approximately five. In total, we use the positional measurements of 82 spectroscopic multiple images belonging to 27 families at z=1.0-6.1 to reconstruct the projected total mass distribution of MACS 1206. Remarkably, 11 multiple images are found within 50 kpc of the brightest cluster galaxy, making this an unprecedented set of constraints for the innermost projected mass distribution of a galaxy cluster. We thus find that, although dynamically relaxed, the smooth matter component (dark matter plus hot gas) of MACS 1206 shows a significant asymmetry, which closely follows the asymmetric distribution of the stellar component (galaxy members and intracluster light). We determine the value of the innermost logarithmic slope of the projected total mass density profile and find it to be close to the canonical Navarro-Frenk-White value. We demonstrate that this quantity is very robust against different parametrizations of the diffuse mass component; however, this is not the case when only one central image is used in the mass reconstruction. We also show that the mass density profile from our new strong lensing model is in very good agreement with dynamical and X-ray measurements at larger radii, where they overlap.

We report the spectroscopic confirmation of 22 new multiply lensed sources behind the Hubble Frontier Field (HFF) galaxy cluster MACS~J0416.1$-$2403 (MACS 0416), using archival data from the Multi Unit Spectroscopic Explorer (MUSE) on the VLT. Combining with previous spectroscopic measurements of 15 other multiply imaged sources, we obtain a sample of 102 secure multiple images with measured redshifts, the largest to date in a single strong lensing system. The newly confirmed sources are largely low-luminosity Lyman-$\alpha$ emitters with redshift in the range [3.08-6.15]. With such a large number of secure constraints, and a significantly improved sample of galaxy members in the cluster core, we have improved our previous strong lensing model and obtained a robust determination of the projected total mass distribution of MACS 0416. We find evidence of three cored dark-matter halos, adding to the known complexity of this merging system. The total mass density profile, as well as the sub-halo population, are found in good agreement with previous works. We update and make public the redshift catalog of MACS 0416 from our previous spectroscopic campaign with the new MUSE redshifts. We also release lensing maps (convergence, shear, magnification) in the standard HFF format.

We present new high resolution and dynamic range dust column density and temperature maps of the California Molecular Cloud derived from a combination of Planck and Herschel dust-emission maps, and 2MASS NIR dust-extinction maps. We used these data to determine the ratio of the 2.2 micron extinction coefficient to the 850 micron opacity and found the value to be close to that found in similar studies of the Orion B and Perseus clouds but higher than that characterizing the Orion A cloud, indicating that variations in the fundamental optical properties of dust may exist between local clouds. We show that over a wide range of extinction, the column density probability distribution function (PDF$_N$) of the cloud can be well described by a simple power law with an index that represents a steeper decline with column density than found in similar studies of the Orion and Perseus clouds. Using only the protostellar population of the cloud and our extinction maps we investigate the Schmidt relation within the cloud. We show that the protostellar surface density, $\Sigma_*$, is directly proportional to the ratio of the protostellar and cloud pdfs. We use the cumulative distribution of protostars to infer the functional forms for both $\Sigma_*$ and PDF$_*$. We find that $\Sigma_*$ is best described by two power-law functions with steeper indicies than found in other local GMCs. We find that the protostellar pdf is a declining function of extinction also best described by two power-laws whose behavior mirrors that of $\Sigma_*$. Our observations suggest that variations both in the slope of the Schmidt relation and in the sizes of the protostellar populations between GMCs are largely driven by variations in the slope of the cloud pdf. This confirms earlier studies suggesting that cloud structure plays a major role in setting the global star formation rates in GMCs.

We present a novel approach for a combined analysis of X-ray and gravitational lensing data and apply this technique to the merging galaxy cluster MACS J0416.1$-$2403. The method exploits the information on the intracluster gas distribution that comes from a fit of the X-ray surface brightness, and then includes the hot gas as a fixed mass component in the strong lensing analysis. With our new technique, we can separate the collisional from the collision-less diffuse mass components, thus obtaining a more accurate reconstruction of the dark matter distribution in the core of a cluster. We introduce an analytical description of the X-ray emission coming from a set of dual Pseudo-Isothermal Elliptical (dPIE) mass distributions, which can be directly used in most lensing softwares. By combining \emph{Chandra} observations with Hubble Frontier Fields imaging and MUSE spectroscopy in MACS J0416.1$-$2403, we measure a projected gas over total mass fraction of approximately $10\%$ at $350$ kpc from the cluster center. Compared to the results of a more traditional cluster mass model (diffuse halos plus member galaxies), we find a significant difference in the cumulative projected mass profile of the dark matter component and that the dark matter to total mass fraction is almost constant, out to more than $350$ kpc. In the coming era of large surveys, these results show the need of multi-probe analyses for detailed dark matter studies in galaxy clusters.

In this paper we present a novel method to identify and characterize stellar clusters deeply embedded in a dark molecular cloud. The method is based on measuring stellar surface density in wide-field infrared images using star counting techniques. It takes advantage of the differing $H$-band luminosity functions (HLFs) of field stars and young stellar populations and is able to statistically associate each star in an image as a member of either the background stellar population or a young stellar population projected on or near the cloud. Moreover, the technique corrects for the effects of differential extinction toward each individual star. We have tested this method against simulations as well as observations. In particular, we have applied the method to 2MASS point sources observed in the Orion A and B complexes, and the results obtained compare very well with those obtained from deep Spitzer and Chandra observations where presence of infrared excess or X-ray emission directly determines membership status for every star. Additionally, our method also identifies unobscured clusters and a low resolution version of the Orion stellar surface density map shows clearly the relatively unobscured and diffuse OB 1a and 1b sub-groups and provides useful insights on their spatial distribution.

The probability distribution function of column density (PDF) has become the tool of choice for cloud structure analysis and star formation studies. Its simplicity is attractive, and the PDF could offer access to cloud physical parameters otherwise difficult to measure, but there has been some confusion in the literature on the definition of its completeness limit and shape at the low column density end. In this Letter we use the natural definition of the completeness limit of a column density PDF, the last closed column-density contour inside a surveyed region, and apply it to a set of large-scale maps of nearby molecular clouds. We conclude that there is no observational evidence for log-normal PDFs in these objects. We find that all studied molecular clouds have PDFs well described by power-laws, including the diffuse cloud Polaris. Our results call for a new physical interpretation for the shape of the column density PDFs. We find that the slope of a cloud PDF is invariant to distance but not to the spatial arrangement of cloud material, and as such it is still a useful tool to investigate cloud structure.

We present a multiwavelength investigation of a region of a nearby giant molecular cloud that is distinguished by a minimal level of star formation activity. With our new 12CO(J=2-1) and 13CO(J=2-1) observations of a remote region within the middle of the California molecular cloud, we aim to investigate the relationship between filaments, cores, and a molecular outflow in a relatively pristine environment. An extinction map of the region from Herschel Space Observatory observations reveals the presence of two 2-pc-long filaments radiating from a high-extinction clump. Using the 13CO observations, we show that the filaments have coherent velocity gradients and that their mass-per-unit-lengths may exceed the critical value above which filaments are gravitationally unstable. The region exhibits structure with eight cores, at least one of which is a starless, prestellar core. We identify a low-velocity, low-mass molecular outflow that may be driven by a flat spectrum protostar. The outflow does not appear to be responsible for driving the turbulence in the core with which it is associated, nor does it provide significant support against gravitational collapse.

Dust extinction is the most robust tracer of the gas distribution in the interstellar medium, but measuring extinction is limited by the systematic uncertainties involved in estimating the intrinsic colors to background stars. In this paper we present a new technique, PNICER, that estimates intrinsic colors and extinction for individual stars using unsupervised machine learning algorithms. This new method aims to be free from any priors with respect to the column density and intrinsic color distribution. It is applicable to any combination of parameters and works in arbitrary numbers of dimensions. Furthermore, it is not restricted to color space. Extinction towards single sources is determined by fitting Gaussian Mixture Models along the extinction vector to (extinction-free) control field observations. In this way it becomes possible to describe the extinction for observed sources with probability densities. PNICER effectively eliminates known biases found in similar methods and outperforms them in cases of deep observational data where the number of background galaxies is significant, or when a large number of parameters is used to break degeneracies in the intrinsic color distributions. This new method remains computationally competitive, making it possible to correctly de-redden millions of sources within a matter of seconds. With the ever-increasing number of large-scale high-sensitivity imaging surveys, PNICER offers a fast and reliable way to efficiently calculate extinction for arbitrary parameter combinations without prior information on source characteristics. PNICER also offers access to the well-established NICER technique in a simple unified interface and is capable of building extinction maps including the NICEST correction for cloud substructure. PNICER is offered to the community as an open-source software solution and is entirely written in Python.

This paper continues our study of the foreground population to the Orion molecular clouds. The goal is to characterize the foreground population north of NGC 1981 and to investigate the star formation history in the large Orion star-forming region. We focus on a region covering about 25 square degrees, centered on the $\epsilon$ Orionis supergiant (HD 37128, B0\,Ia) and covering the Orion Belt asterism. We used a combination of optical (SDSS) and near-infrared (2MASS) data, informed by X-ray (\textit{XMM-Newton}) and mid-infrared (WISE) data, to construct a suite of color-color and color-magnitude diagrams for all available sources. We then applied a new statistical multiband technique to isolate a previously unknown stellar population in this region. We identify a rich and well-defined stellar population in the surveyed region that has about 2\,000 objects that are mostly M stars. We infer the age for this new population to be at least 5\, Myr and likely $\sim10$\,Myr and estimate a total of about 2\,500 members, assuming a normal IMF. This new population, which we call the Orion Belt population, is essentially extinction-free, disk-free, and its spatial distribution is roughly centered near $\epsilon$ Ori, although substructure is clearly present. The Orion Belt population is likely the low-mass counterpart to the Ori OB Ib subgroup. Although our results do not rule out Blaauw's sequential star formation scenario for Orion, we argue that the recently proposed blue streams scenario provides a better framework on which one can explain the Orion star formation region as a whole. We speculate that the Orion Belt population could represent the evolved counterpart of an Orion nebula-like cluster.

In this Letter we compare the abundance of member galaxies of a rich, nearby ($z=0.09$) galaxy cluster, Abell 2142, with that of halos of comparable virial mass extracted from sets of state-of-the-art numerical simulations, both collisionless at different resolutions and with the inclusion of baryonic physics in the form of cooling, star formation, and feedback by active galactic nuclei. We also use two semi-analytical models to account for the presence of orphan galaxies. The photometric and spectroscopic information, taken from the Sloan Digital Sky Survey Data Release 12 (SDSS DR12) database, allows us to estimate the stellar velocity dispersion of member galaxies of Abell 2142. This quantity is used as proxy for the total mass of secure cluster members and is properly compared with that of subhalos in simulations. We find that simulated halos have a statistically significant ($\gtrsim 7$ sigma confidence level) smaller amount of massive (circular velocity above $200\,{\rm km\, s^{-1}}$) subhalos, even before accounting for the possible incompleteness of observations. These results corroborate the findings from a recent strong lensing study of the Hubble Frontier Fields galaxy cluster MACS J0416 \citep{grillo2015} and suggest that the observed difference is already present at the level of dark matter (DM) subhalos and is not solved by introducing baryonic physics. A deeper understanding of this discrepancy between observations and simulations will provide valuable insights into the impact of the physical properties of DM particles and the effect of baryons on the formation and evolution of cosmological structures.

We present VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $\Delta V_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $\Delta V_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.

We present MUSE observations in the core of the HFF galaxy cluster MACS J1149.5+2223, where the first magnified and spatially-resolved multiple images of SN 'Refsdal' at redshift 1.489 were detected. Thanks to a DDT program with the VLT and the extraordinary efficiency of MUSE, we measure 117 secure redshifts with just 4.8 hours of total integration time on a single target pointing. We spectroscopically confirm 68 galaxy cluster members, with redshift values ranging from 0.5272 to 0.5660, and 18 multiple images belonging to 7 background, lensed sources distributed in redshifts between 1.240 and 3.703. Starting from the combination of our catalog with those obtained from extensive spectroscopic and photometric campaigns using the HST, we select a sample of 300 (164 spectroscopic and 136 photometric) cluster members, within approximately 500 kpc from the BCG, and a set of 88 reliable multiple images associated to 10 different background source galaxies and 18 distinct knots in the spiral galaxy hosting SN 'Refsdal'. We exploit this valuable information to build 6 detailed strong lensing models, the best of which reproduces the observed positions of the multiple images with a rms offset of only 0.26". We use these models to quantify the statistical and systematic errors on the predicted values of magnification and time delay of the next emerging image of SN 'Refsdal'. We find that its peak luminosity should should occur between March and June 2016, and should be approximately 20% fainter than the dimmest (S4) of the previously detected images but above the detection limit of the planned HST/WFC3 follow-up. We present our two-dimensional reconstruction of the cluster mass density distribution and of the SN 'Refsdal' host galaxy surface brightness distribution. We outline the roadmap towards even better strong lensing models with a synergetic MUSE and HST effort.

We use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at $z=0.44$ to estimate the value of the ratio $\eta=\Psi/\Phi$ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric.[...] Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only $\Phi$), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum $\Phi+\Psi$). We thus express $\eta$ as a function of the mass profiles and perform our analysis over the radial range $0.5\,Mpc\le r\le r_{200}=1.96\,Mpc$. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain $\eta(r_{200})=1.01\,_{-0.28}^{+0.31}$ at the 68\% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied $f(R)$ modified gravity model, the constraint on $\eta$ translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the $f(R)$ interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.

We present a strong lensing system, composed of 4 multiple images of a source at z = 2.387, created by two lens galaxies, G1 and G2, belonging to the galaxy cluster MACS J1115.9+0129 at z = 0.353. We use observations taken as part of the Cluster Lensing and Supernova survey with Hubble, CLASH, and its spectroscopic follow-up programme at the Very Large Telescope, CLASH-VLT, to estimate the total mass distributions of the two galaxies and the cluster through strong gravitational lensing models. We find that the total projected mass values within the half-light radii, R_{e}, of the two lens galaxies are M_{T,G1}(< R_{e,G1}) = (3.6 +/- 0.4) x 10^{10}M_{Sun} and M_{T,G2}(< R_{e,G2}) = (4.2 +/- 1.6) x 10^{10}M_{Sun}. The effective velocity dispersion values of G1 and G2 are (122 +/- 7) km/s and (137 +/- 27) km/s, respectively. We remark that these values are relatively low when compared to those of ~200-300 km/s, typical of lens galaxies found in the field by previous surveys. By fitting the spectral energy distributions of G1 and G2, we measure projected luminous over total mass fractions within R_{e} of 0.11 +/- 0.03, for G1, and 0.73 +/- 0.32, for G2. The fact that the less massive galaxy, G1, is dark-matter dominated in its inner regions raises the question of whether the dark matter fraction in the core of early-type galaxies depends on their mass. Further investigating strong lensing systems will help us understand the influence that dark matter has on the structure and evolution of the inner regions of galaxies.

Orion A hosts the nearest massive star factory, thus offering a unique opportunity to resolve the processes connected with the formation of both low- and high-mass stars. Here we present the most detailed and sensitive near-infrared (NIR) observations of the entire molecular cloud to date. With the unique combination of high image quality, survey coverage, and sensitivity, our NIR survey of Orion A aims at establishing a solid empirical foundation for further studies of this important cloud. In this first paper we present the observations, data reduction, and source catalog generation. To demonstrate the data quality, we present a first application of our catalog to estimate the number of stars currently forming inside Orion A and to verify the existence of a more evolved young foreground population. We used the European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) to survey the entire Orion A molecular cloud in the NIR $J, H$, and $K_S$ bands, covering a total of $\sim$18.3 deg$^2$. We implemented all data reduction recipes independently of the ESO pipeline. Estimates of the young populations toward Orion A are derived via the $K_S$-band luminosity function. Our catalog (799995 sources) increases the source counts compared to the Two Micron All Sky Survey by about an order of magnitude. The 90% completeness limits are 20.4, 19.9, and 19.0 mag in $J, H$, and $K_S$, respectively. The reduced images have 20% better resolution on average compared to pipeline products. We find between 2300 and 3000 embedded objects in Orion A and confirm that there is an extended foreground population above the Galactic field, in agreement with previous work. The Orion A VISTA catalog represents the most detailed NIR view of the nearest massive star-forming region and provides a fundamental basis for future studies of star formation processes toward Orion.

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $\Lambda$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $\Omega_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $\Lambda$CDM model, and $\Omega_m=0.31^{+0.12}_{-0.13}$ and $\Omega_\Lambda=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

We present optical depth and temperature maps of the Perseus molecular cloud, obtained combining dust emission data from the Herschel and Planck satellites and 2MASS/NIR dust extinction maps. The maps have a resolution of 36 arcsec in the Herschel regions, and of 5 arcmin elsewhere. The dynamic range of the optical depth map ranges from $1\times10^{-2}\, \mathrm{mag}$ up to $20 \,\mathrm{mag}$ in the equivalent K band extinction. We also evaluate the ratio between the $2.2 \,\mathrm{\mu m}$ extinction coefficient and the $850 \,\mathrm{\mu m}$ opacity. The value we obtain is close to the one found in the Orion B molecular cloud. We show that the cumulative and the differential area function of the data (which is proportional to the probability distribution function of the cloud column density) follow power laws with index respectively $\simeq -2$, and $\simeq -3$. We use WISE data to improve current YSO catalogues based mostly on \emph{Spitzer} data and we build an up-to-date selection of Class~I/0 objects. Using this selection, we evaluate the local Schmidt law, $\Sigma_{\mathrm{YSO}} \propto \Sigma_{\mathrm{gas}}^{\beta}$, showing that $\beta=2.4 \pm 0.6$. Finally, we show that the area-extinction relation is important for determining the star formation rate in the cloud, which is in agreement with other recent works.

We present the optical luminosity functions (LFs) of galaxies for the CLASH-VLT cluster MACS J1206.2-0847 at z=0.439, based on HST and SUBARU data, including ~600 spectroscopically confirmed member galaxies. The LFs on the wide SUBARU FoV are well described by a single Schechter function down to M~M*+3, whereas this fit is poor for HST data, due to a faint-end upturn visible down M~M*+7, suggesting a bimodal behaviour. We also investigate the effect of local environment by deriving the LFs in four different regions, according to the distance from the centre, finding an increase in the faint-end slope going from the core to the outer rings. Our results confirm and extend our previous findings on the analysis of mass functions, which showed that the galaxies with stellar mass below 10^10.5, M_sun have been significantly affected by tidal interaction effects, thus contributing to the intra cluster light.

The analysis of galaxy properties and the relations among them and the environment, can be used to investigate the physical processes driving galaxy evolution. We study the cluster A209 by using the CLASH-VLT spectroscopic data combined with Subaru photometry, yielding to 1916 cluster members down to a stellar mass of 10^{8.6} Msun. We determine: i) the stellar mass function of star-forming and passive galaxies; ii) the intra-cluster light and its properties; iii) the orbits of low- and high-mass passive galaxies; and iv) the mass-size relation of ETGs. The stellar mass function of the star-forming galaxies does not depend on the environment, while the slope found for passive galaxies becomes flatter in the densest region. The color distribution of the intra-cluster light is consistent with the color of passive members. The analysis of the dynamical orbits shows that low-mass passive galaxies have tangential orbits, avoiding small pericenters around the BCG. The mass-size relation of low-mass passive ETGs is flatter than that of high mass galaxies, and its slope is consistent with that of field star-forming galaxies. Low-mass galaxies are also more compact within the scale radius of 0.65 Mpc. The ratio between stellar and number density profiles shows a mass segregation in the center. The comparative analysis of the stellar and total density profiles indicates that this effect is due to dynamical friction. Our results are consistent with a scenario in which the "environmental quenching" of low-mass galaxies is due to mechanisms such as harassment out to R200, starvation and ram-pressure stripping at smaller radii, as supported by the analysis of the mass function, of the dynamical orbits and of the mass-size relation of passive early-types in different regions. Our analyses support the idea that the intra-cluster light is formed through the tidal disruption of subgiant galaxies.

(abridged) Methods: We derive maps of submillimeter dust optical depth and effective dust temperature from Herschel data that were calibrated against Planck. After calibration, we then fit a modified blackbody to the long-wavelength Herschel data, using the Planck-derived dust opacity spectral index beta, derived on scales of 30' (or ~1 pc). We use this model to make predictions of the submillimeter flux density at 850 micron, and we compare these in turn with APEX-Laboca observations. Results: A comparison of the submillimeter dust optical depth and near-infrared extinction data reveals evidence for an increased submillimeter dust opacity at high column densities, interpreted as an indication of grain growth in the inner parts of the core. Additionally, a comparison of the Herschel dust model and the Laboca data reveals that the frequency dependence of the submillimeter opacity, described by the spectral index beta, does not change. A single beta that is only slightly different from the Planck-derived value is sufficient to describe the data, beta=1.53+/-0.07. We apply a similar analysis to Barnard 68, a core with significantly lower column densities than FeSt 1-457, and we do not find evidence for grain growth but also a single beta. Conclusions: While we find evidence for grain growth from the dust opacity in FeSt 1-457, we find no evidence for significant variations in the dust opacity spectral index beta on scales 0.02<x<1 pc (or 36"<x<30'). The correction to the Planck-derived dust beta that we find in both cases is on the order of the measurement error, not including any systematic errors, and it would thus be reasonable to directly apply the dust beta from the Planck all-sky dust model. As a corollary, reliable effective temperature maps can be derived which would be otherwise affected by beta variations.

The Orion cloud complex presents a variety of star formation mechanisms and properties and it is still one of the most intriguing targets for star formation studies. We present VISTA/VIRCAM near-infrared observations of the L1630N star forming region, including the stellar clusters NGC 2068 and NGC 2071, in the Orion molecular cloud B and discuss them in combination with Spitzer data. We select 186 young stellar object (YSO) candidates in the region on the basis of multi-colour criteria, confirm the YSO nature of the majority of them using published spectroscopy from the literature, and use this sample to investigate the overall star formation properties in L1630N. The K-band luminosity function of L1630N is remarkably similar to that of the Trapezium cluster, i.e., it presents a broad peak in the range 0.3-0.7 M$_\odot$ and a fraction of sub-stellar objects of $\sim$20%. The fraction of YSOs still surrounded by disk/envelopes is very high ($\sim$85%) compared to other star forming regions of similar age (1-2 Myr), but includes some uncertain corrections for diskless YSOs. Yet, a possibly high disk fraction together with the fact that 1/3 of the cloud mass has a gas surface density above the threshold for star formation ($\sim$129 M$_\odot$ pc$^{-2}$), points towards a still on-going star formation activity in L1630N. The star formation efficiency (SFE), star formation rate (SFR) and density of star formation of L1630N are within the ranges estimated for galactic star forming regions by the Spitzer "core to disk" and "Gould's Belt" surveys. However, the SFE and SFR are lower than the average value measured in the Orion A cloud and, in particular, lower than that in the southern regions of L1630. This might suggest different star formation mechanisms within the L1630 cloud complex.

A deep, wide-field, near-infrared imaging survey was used to construct an extinction map of the southeastern part of the California Molecular Cloud (CMC) with $\sim$ 0.5 arc min resolution. The same region was also surveyed in the $^{12}$CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1) emission lines at the same angular resolution. Strong spatial variations in the abundances of $^{13}$CO and C$^{18}$O were found to be correlated with variations in gas temperature, consistent with temperature dependent CO depletion/desorption on dust grains. The $^{13}$CO to C$^{18}$O abundance ratio was found to increase with decreasing extinction, suggesting selective photodissociation of C$^{18}$O by the ambient UV radiation field. The cloud averaged X-factor is found to be $<$X$_{\rm CO}$$>$ $=$ 2.53 $\times$ 10$^{20}$ ${\rm cm}^{-2}~({\rm K~km~s}^{-1})^{-1}$, somewhat higher than the Milky Way average. On sub-parsec scales we find no single empirical value of the X-factor that can characterize the molecular gas in cold (T$_{\rm k}$ $\lesssim$ 15 K) regions, with X$_{\rm CO}$ $\propto$ A$_{\rm V}$$^{0.74}$ for A$_{\rm V}$ $\gtrsim$ 3 magnitudes. However in regions containing relatively hot (T$_{\rm ex}$ $\gtrsim$ 25 K) gas we find a clear correlation between W($^{12}$CO) and A$_{\rm V}$ over a large (3 $\lesssim$ A$_{\rm V}$ $\lesssim$ 25 mag) extinction range. This suggests a constant X$_{\rm CO}$ $=$ 1.5 $\times$ 10$^{20}$ ${\rm cm}^{-2}~({\rm K~km~s}^{-1})^{-1}$ for the hot gas, a lower value than either the average for the CMC or Milky Way. We find a correlation between X$_{\rm CO}$ and T$_{\rm ex}$ with X$_{\rm CO}$ $\propto$ T$_{\rm ex}$$^{-0.7}$ suggesting that the global X-factor of a cloud may depend on the relative amounts of hot gas within it.

In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 at z~0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from the CLASH-VLT survey. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as part of our ESO Large Programme, sampling the cluster out to a radius ~2R200 (4 Mpc). We classify 412 galaxies as passive, with strong Hdelta absorption (red and blue galaxies, and with emission lines from weak to very strong. A number of tests for substructure detection are applied to analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hdelta galaxies trace the cluster center and the WNW-ESE elongated structure. The red strong Hdelta galaxies also mark a secondary, dense peak ~2 Mpc at ESE. The emission line galaxies cluster in several loose structures, mostly outside R200. The observational scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red strong Hdelta galaxies. The red strong Hdelta galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions.

In this Letter we investigate the shape of the probability distribution of column densities (PDF) in molecular clouds. Through the use of low-noise, extinction-calibrated \textit{Herschel}/\textit{Planck} emission data for eight molecular clouds, we demonstrate that, contrary to common belief, the PDFs of molecular clouds are not described well by log-normal functions, but are instead power laws with exponents close to two and with breaks between $A_K \simeq 0.1$ and $0.2\,\mathrm{mag}$, so close to the CO self-shielding limit and not far from the transition between molecular and atomic gas. Additionally, we argue that the intrinsic functional form of the PDF cannot be securely determined below $A_K \simeq 0.1\,\mathrm{mag}$, limiting our ability to investigate more complex models for the shape of the cloud PDF.

Galaxy intrinsic alignment can be a severe source of error in weak-lensing studies. The problem has been widely studied by numerical simulations and with heuristic models, but without a clear theoretical justification of its origin and amplitude. In particular, it is still unclear whether intrinsic alignment of galaxies is dominated by formation and accretion processes or by the effects of the instantaneous tidal field acting upon them. We investigate this question by developing a simple model of intrinsic alignment for elliptical galaxies, based on the instantaneous tidal field. Making use of the galaxy stellar distribution function, we estimate the intrinsic alignment signal and find that although it has the expected dependence on the tidal field, it is too weak to account for the observed signal. This is an indirect validation of the standard view that intrinsic alignment is caused by formation and/or accretion processes.

Spectra of composite systems (e.g., spectroscopic binaries) contain spatial information that can be retrieved by measuring the radial velocities (i.e., Doppler shifts) of the components in four observations with the slit rotated by 90 degrees in the sky. By using basic concepts of slit spectroscopy we show that the geometry of composite systems can be reliably retrieved by measuring only radial velocity differences taken with different slit angles. The spatial resolution is determined by the precision with which differential radial velocities can be measured. We use the UVES spectrograph at the VLT to observe the known spectroscopic binary star HD 188088 (HIP 97944), which has a maximum expected separation of 23 milli-arcseconds. We measure an astrometric signal in radial velocity of 276 \ms, which corresponds to a separation between the two components at the time of the observations of 18 $\pm2$ milli-arcseconds. The stars were aligned east-west. We describe a simple optical device to simultaneously record pairs of spectra rotated by 180 degrees, thus reducing systematic effects. We compute and provide the function expressing the shift of the centroid of a seeing-limited image in the presence of a narrow slit.The proposed technique is simple to use and our test shows that it is amenable for deriving astrometry with milli-arcsecond accuracy or better, beyond the diffraction limit of the telescope. The technique can be further improved by using simple devices to simultaneously record the spectra with 180 degrees angles.With tachoastrometry, radial velocities and astrometric positions can be measured simultaneously for many double line system binaries in an easy way. The method is not limited to binary stars, but can be applied to any astrophysical configuration in which spectral lines are generated by separate (non-rotational symmetric) regions.

We present a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

We present the first observations of the Frontier Fields Cluster Abell S1063 taken with the newly commissioned Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph. Because of the relatively large field of view (1 arcmin^2), MUSE is ideal to simultaneously target multiple galaxies in blank and cluster fields over the full optical spectrum. We analysed the four hours of data obtained in the Science Verification phase on this cluster and measured redshifts for 53 galaxies. We confirm the redshift of five cluster galaxies, and determine the redshift of 29 other cluster members. Behind the cluster, we find 17 galaxies at higher redshift, including three previously unknown Lyman-alpha emitters at z>3, and five multiply-lensed galaxies. We report the detection of a new z=4.113 multiply lensed galaxy, with images that are consistent with lensing model predictions derived for the Frontier Fields. We detect C III], C IV, and He II emission in a multiply lensed galaxy at z=3.116, suggesting the likely presence of an active galactic nucleus. We also created narrow-band images from the MUSE datacube to automatically search for additional line emitters corresponding to high-redshift candidates, but we could not identify any significant detections other than those found by visual inspection. With the new redshifts, it will become possible to obtain an accurate mass reconstruction in the core of Abell S1063 through refined strong lensing modelling. Overall, our results illustrate the breadth of scientific topics that can be addressed with a single MUSE pointing. We conclude that MUSE is a very efficient instrument to observe galaxy clusters, enabling their mass modelling, and to perform a blind search for high-redshift galaxies.

Context. The study of the galaxy stellar mass function (SMF) in relation to the galaxy environment and the stellar mass density profile, rho(r), is a powerful tool to constrain models of galaxy evolution. Aims. We determine the SMF of the z=0.44 cluster of galaxies MACS J1206.2-0847 separately for passive and star-forming (SF) galaxies, in different regions of the cluster, from the center out to approximately 2 virial radii. We also determine rho(r) to compare it to the number density and total mass density profiles. Methods. We use the dataset from the CLASH-VLT survey. Stellar masses are obtained by SED fitting on 5-band photometric data obtained at the Subaru telescope. We identify 1363 cluster members down to a stellar mass of 10^9.5 Msolar. Results. The whole cluster SMF is well fitted by a double Schechter function. The SMFs of cluster SF and passive galaxies are statistically different. The SMF of the SF cluster galaxies does not depend on the environment. The SMF of the passive population has a significantly smaller slope (in absolute value) in the innermost (<0.50 Mpc), highest density cluster region, than in more external, lower density regions. The number ratio of giant/subgiant galaxies is maximum in this innermost region and minimum in the adjacent region, but then gently increases again toward the cluster outskirts. This is also reflected in a decreasing radial trend of the average stellar mass per cluster galaxy. On the other hand, the stellar mass fraction, i.e., the ratio of stellar to total cluster mass, does not show any significant radial trend. Conclusions. Our results appear consistent with a scenario in which SF galaxies evolve into passive galaxies due to density-dependent environmental processes, and eventually get destroyed very near the cluster center to become part of a diffuse intracluster medium.

(abridged) [...] Methods: In a continued study of the molecular core population of the Pipe Nebula, we present a molecular-line survey of 52 cores. Previous research has shown a variety of different chemical evolutionary stages among the cores. Using the Mopra radio telescope, we observed the ground rotational transitions of HCO+, H13CO+, HCN, H13CN, HNC, and N2H+. These data are complemented with near-infrared extinction maps to constrain the column densities, effective dust temperatures derived from Herschel data, and NH3-based gas kinetic temperatures. Results: The target cores are located across the nebula, span visual extinctions between 5 and 67 mag, and effective dust temperatures (averaged along the lines of sight) between 13 and 19 K. The extinction-normalized integrated line intensities, a proxy for the abundance in constant excitation conditions of optically thin lines, vary within an order of magnitude for a given molecule. The effective dust temperatures and gas kinetic temperatures are correlated, but the effective dust temperatures are consistently higher than the gas kinetic temperatures. Combining the molecular line and temperature data, we find that N2H+ is only detected toward the coldest and densest cores while other lines show no correlation with these core properties. Conclusions: Within this large sample, N2H+ is the only species to exclusively trace the coldest and densest cores, in agreement with chemical considerations. In contrast, the common high-density tracers HCN and HNC are present in a majority of cores, demonstrating the utility of these molecules to characterize cores over a large range of extinctions. The correlation between the effective dust temperatures and the gas kinetic temperatures suggests that the former are dominated by dust that is both dense and thermodynamically coupled to the dense gas traced by NH3. [...]

We present high-resolution, high dynamic range column-density and color-temperature maps of the Orion complex using a combination of Planck dust-emission maps, Herschel dust-emission maps, and 2MASS NIR dust-extinction maps. The column-density maps combine the robustness of the 2MASS NIR extinction maps with the resolution and coverage of the Herschel and Planck dust-emission maps and constitute the highest dynamic range column-density maps ever constructed for the entire Orion complex, covering $0.01 \, \mathrm{mag} < A_K < 30 \,\mathrm{mag}$, or $2 \times 10^{20} \, \mathrm{cm}^{-2} < N < 5 \times 10^{23} \,\mathrm{cm}^{-2}$. We determined the ratio of the 2.2 microns extinction coefficient to the 850 microns opacity and found that the values obtained for both Orion A and B are significantly lower than the predictions of standard dust models, but agree with newer models that incorporate icy silicate-graphite conglomerates for the grain population. We show that the cloud projected pdf, over a large range of column densities, can be well fitted by a simple power law. Moreover, we considered the local Schmidt-law for star formation, and confirm earlier results, showing that the protostar surface density $\Sigma_*$ follows a simple law $\Sigma_* \propto \Sigma_{gas}^\beta$, with $\beta \sim 2$.

We aim at constraining the assembly history of clusters by studying the intra cluster light (ICL) properties, estimating its contribution to the fraction of baryons in stars, f*, and understanding possible systematics/bias using different ICL detection techniques. We developed an automated method, GALtoICL, based on the software GALAPAGOS to obtain a refined version of typical BCG+ICL maps. We applied this method to our test case MACS J1206.2-0847, a massive cluster located at z=0.44, that is part of the CLASH sample. Using deep multi-band SUBARU images, we extracted the surface brightness (SB) profile of the BCG+ICL and we studied the ICL morphology, color, and contribution to f* out to R500. We repeated the same analysis using a different definition of the ICL, SBlimit method, i.e. a SB cut-off level, to compare the results. The most peculiar feature of the ICL in MACS1206 is its asymmetric radial distribution, with an excess in the SE direction and extending towards the 2nd brightest cluster galaxy which is a Post Starburst galaxy. This suggests an interaction between the BCG and this galaxy that dates back to t <= 1.5 Gyr. The BCG+ICL stellar content is 8% of M_(*,500) and the (de-) projected baryon fraction in stars is f*=0.0177 (0.0116), in excellent agreement with recent results. The SBlimit method provides systematically higher ICL fractions and this effect is larger at lower SB limits. This is due to the light from the outer envelopes of member galaxies that contaminate the ICL. Though more time consuming, the GALtoICL method provides safer ICL detections that are almost free of this contamination. This is one of the few ICL study at redshift z > 0.3. At completion, the CLASH/VLT program will allow us to extend this analysis to a statistically significant cluster sample spanning a wide redshift range: 0.2<z<0.6.

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2\rho)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $\rho$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

We present a near-infrared extinction map of a large region ($\sim$870 deg$^2$) covering the isolated Corona Australis complex of molecular clouds. We reach a 1-$\sigma$ error of 0.02 mag in the K-band extinction with a resolution of 3 arcmin over the entire map. We find that the Corona Australis cloud is about three times as large as revealed by previous CO and dust emission surveys. The cloud consists of a 45 pc long complex of filamentary structure from the well known star forming Western-end (the head, $N \geq10^{23}$ cm$^{-2}$) to the diffuse Eastern-end the tail, ($N \leq10^{21}$ cm$^{-2}$). Remarkably, about two thirds of the complex both in size and mass lie beneath A$_V\sim1$ mag. We find that the PDF of the cloud cannot be described by a single log-normal function. Similar to prior studies, we found a significant excess at high column densities, but a log-normal + power-law tail fit does not work well at low column densities. We show that at low column densities near the peak of the observed PDF, both the amplitude and shape of the PDF are dominated by noise in the extinction measurements making it impractical to derive the intrinsic cloud PDF below A$_K <$ 0.15 mag. Above A$_K \sim 0.15$ mag, essentially the molecular component of the cloud, the PDF appears to be best described by a power-law with index $-3$, but could also described as the tail of a broad and relatively low amplitude, log-normal PDF that peaks at very low column densities.

We present VIsible Multi-Object Spectrograph (VIMOS) observations of a z 6 galaxy quintuply imaged by the Frontier Fields galaxy cluster RXC J2248.7-4431 (z=0.348). This sub-L^*, high-z galaxy has been recently discovered by Monna et al. (2013) using dropout techniques with the 16-band HST photometry acquired as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). Obtained as part of the CLASH-VLT survey, the VIMOS medium-resolution spectra of this source show a very faint continuum between ~8700A and ~9300A and a prominent emission line at 8643A, which can be readily identified with Lyman-alpha at z=6.110. The emission line exhibits an asymmetric profile, with a more pronounced red wing. The rest-frame equivalent width of the line is EW=79+-10A. After correcting for magnification, the star-formation rate (SFR) estimated from the Lya line is SFR(Lya)=11 M_{sol}/yr and that estimated from the UV data is SFR(UV)=3 M_{sol}/yr. We estimate that the effective radius of the source is R_e<~0.4 kpc, which implies a star formation surface mass density Sigma_{SFR}>6 M_{sol}/yr/kpc^2 and, using the Kennicutt-Schmidt relation, a gas surface mass density Sigma_{gas}>10^3 M_{sol}/pc^2. Our results support the idea that this magnified, distant galaxy is a young and compact object with 0.4 L^* at z=6, with comparable amount of mass in gas and stars. Future follow-up observations with ALMA will provide valuable insight into the SFR and molecular gas content of this source. In the spirit of the Frontier Fields initiative, we also publish the redshifts of several multiply imaged sources and other background objects which will help improving the strong lensing model of this galaxy cluster.

We consider the general problem of fitting a parametric density model to discrete observations, taken to follow a non-homogeneous Poisson point process. This class of models is very common, and can be used to describe many astrophysical processes, including the distribution of protostars in molecular clouds. We give the expression for the likelihood of a given spatial density distribution of protostars and apply it to infer the most probable dependence of the protostellar surface density on the gas surface density. Finally, we apply this general technique to model the distribution of protostars in the Orion molecular cloud and robustly derive the local star formation scaling (Schmidt) law for a molecular cloud. We find that in this cloud the protostellar surface density, $\Sigma_\mathrm{YSO}$, is directly proportional to the square gas column density, here expressed as infrared extinction in the $K$-band, $A_K$: more precisely, $\Sigma_\mathrm{YSO} = (1.65 \pm 0.19) A_K^{(2.03 \pm 0.15)}$ stars pc$^{-2}$.

We investigate Schmidt's conjecture (i.e., that the star formation rate scales in a power-law fashion with the gas density) for four well-studied local molecular clouds (GMCs). Using the Bayesian methodology we show that a local Schmidt scaling relation of the form Sigma*(A_K) = kappa x (A_K)^{beta} (protostars pc^{-2}) exists within (but not between) GMCs. Further we find that the Schmidt scaling law, by itself, does not provide an adequate description of star formation activity in GMCs. Because the total number of protostars produced by a cloud is given by the product of Sigma*(A_K) and S'(> A_K), the differential surface area distribution function, integrated over the entire cloud, the cloud's structure plays a fundamental role in setting the level of its star formation activity. For clouds with similar functional forms of Sigma*(A_K), observed differences in their total SFRs are primarily due to the differences in S'(> A_K) between the clouds. The coupling of Sigma*(A_K) with the measured S'(> A_K) in these clouds also produces a steep jump in the SFR and protostellar production above A_K ~ 0.8 magnitudes. Finally, we show that there is no global Schmidt law that relates the star formation rate and gas mass surface densities between GMCs. Consequently, the observed Kennicutt-Schmidt scaling relation for disk galaxies is likely an artifact of unresolved measurements of GMCs and not a result of any underlying physical law of star formation characterizing the molecular gas.

We use an unprecedented data-set of about 600 redshifts for cluster members, obtained as part of a VLT/VIMOS large programme, to constrain the mass profile of the z=0.44 cluster MACS J1206.2-0847 over the radial range 0-5 Mpc (0-2.5 virial radii) using the MAMPOSSt and Caustic methods. We then add external constraints from our previous gravitational lensing analysis. We invert the Jeans equation to obtain the velocity-anisotropy profiles of cluster members. With the mass-density and velocity-anisotropy profiles we then obtain the first determination of a cluster pseudo-phase-space density profile. The kinematics and lensing determinations of the cluster mass profile are in excellent agreement. This is very well fitted by a NFW model with mass M200=(1.4 +- 0.2) 10^15 Msun and concentration c200=6 +- 1, only slightly higher than theoretical expectations. Other mass profile models also provide acceptable fits to our data, of (slightly) lower (Burkert, Hernquist, and Softened Isothermal Sphere) or comparable (Einasto) quality than NFW. The velocity anisotropy profiles of the passive and star-forming cluster members are similar, close to isotropic near the center and increasingly radial outside. Passive cluster members follow extremely well the theoretical expectations for the pseudo-phase-space density profile and the relation between the slope of the mass-density profile and the velocity anisotropy. Star-forming cluster members show marginal deviations from theoretical expectations. This is the most accurate determination of a cluster mass profile out to a radius of 5 Mpc, and the only determination of the velocity-anisotropy and pseudo-phase-space density profiles of both passive and star-forming galaxies for an individual cluster [abridged]

Context. The properties of dust grains, in particular their size distribution, are expected to differ from the interstellar medium to the high-density regions within molecular clouds. Aims. We measure the mid-infrared extinction law produced by dense material in molecular cloud cores. Since the extinction at these wavelengths is caused by dust, the extinction law in cores should depart from that found in low-density environments if the dust grains have different properties. Methods. We use the unbiased LINES method to measure the slope of the reddening vectors in color-color diagrams. We derive the mid-infrared extinction law toward the dense cores B59 and FeSt 1-457 in the Pipe Nebula over a range of visual extinction between 10 and 50 magnitudes, using a combination of Spitzer/IRAC, and ESO NTT/VLT data. Results. The mid-infrared extinction law in both cores departs significantly from a power-law between 3.6 and 8 micron, suggesting that these cores contain dust with a considerable fraction of large dust grains. We find no evidence for a dependence of the extinction law with column density up to 50 magnitudes of visual extinction in these cores, and no evidence for a variation between our result and those for other clouds at lower column densities reported elsewhere in the literature. This suggests that either large grains are present even in low column density regions, or that the existing dust models need to be revised at mid-infrared wavelengths. We find a small but significant difference in the extinction law of the two cores, that we tentatively associate with the onset of star formation in B59.

The spatial variation of the colour of a galaxy may introduce a bias in the measurement of its shape if the PSF profile depends on wavelength. We study how this bias depends on the properties of the PSF and the galaxies themselves. The bias depends on the scales used to estimate the shape, which may be used to optimise methods to reduce the bias. Here we develop a general approach to quantify the bias. Although applicable to any weak lensing survey, we focus on the implications for the ESA Euclid mission. Based on our study of synthetic galaxies we find that the bias is a few times 10^-3 for a typical galaxy observed by Euclid. Consequently, it cannot be neglected and needs to be accounted for. We demonstrate how one can do so using spatially resolved observations of galaxies in two filters. We show that HST observations in the F606W and F814W filters allow us to model and reduce the bias by an order of magnitude, sufficient to meet Euclid's scientific requirements. The precision of the correction is ultimately determined by the number of galaxies for which spatially-resolved observations in at least two filters are available. We use results from the Millennium Simulation to demonstrate that archival HST data will be sufficient for the tomographic cosmic shear analysis with the Euclid dataset.

Despite over 30 years of study, the mass-area relationship within and among clouds is still poorly understood both observationally and theoretically. Modern extinction datasets should have sufficient resolution and dynamic range to characterize this relationship for nearby molecular clouds, although recent papers using extinction data seem to yield different interpretations regarding the nature and universality of this aspect of cloud structure. In this paper we try to unify these various results and interpretations by accounting for the different ways cloud properties are measured and analyzed. We interpret the mass-area relationship in terms of the column density distribution function and its possible variation within and among clouds. We quantitatively characterize regional variations in the column density PDF. We show that structures both within and among clouds possess the same degree of "universality", in that their PDF means do not systematically scale with structure size. Because of this, mass scales linearly with area.

Context. The properties of dust grains, in particular their size distribution, are expected to differ from the interstellar medium to the high-density regions within molecular clouds. Since the extinction at near-infrared wavelengths is caused by dust, the extinction law in cores should depart from that found in low-density environments if the dust grains have different properties. Aims. We explore methods to measure the near-infrared extinction law produced by dense material in molecular cloud cores from photometric data. Methods. Using controlled sets of synthetic and semi-synthetic data, we test several methods for linear regression applied to the specific problem of deriving the extinction law from photometric data. We cover the parameter space appropriate to this type of observations. Results. We find that many of the common linear-regression methods produce biased results when applied to the extinction law from photometric colors. We propose and validate a new method, LinES, as the most reliable for this effect. We explore the use of this method to detect whether or not the extinction law of a given reddened population has a break at some value of extinction.

In this article I present IEAD, a new interface for astronomical science databases. It is based on a powerful, yet simple, syntax designed to completely abstract the user from the structure of the underlying database. The programming language chosen for its implementation, JavaScript, makes it possible to interact directly with the user and to provide real-time information on the parsing process, error messages, and the name resolution of targets; additionally, the same parsing engine is used for context-sensitive autocompletion. Ultimately, this product should significantly simplify the use of astronomical archives, inspire more advanced uses of them, and allow the user to focus on what scientific research to perform, instead of on how to instruct the computer to do it.

In this paper we investigate scaling relations between star formation rates and molecular gas masses for both local Galactic clouds and a sample of external galaxies. We specifically consider relations between the star formation rates and measurements of dense, as well as total, molecular gas masses. We argue that there is a fundamental empirical scaling relation that directly connects the local star formation process with that operating globally within galaxies. Specifically, the total star formation rate in a molecular cloud or galaxy is linearly proportional to the mass of dense gas within the cloud or galaxy. This simple relation, first documented in previous studies, holds over a span of mass covering nearly nine orders of magnitude and indicates that the rate of star formation is directly controlled by the amount of dense molecular gas that can be assembled within a star formation complex. We further show that the star formation rates and total molecular masses, characterizing both local clouds and galaxies, are correlated over similarly large scales of mass and can be described by a family of linear star formation scaling laws, parameterized by $f_{DG}$, the fraction of dense gas contained within the clouds or galaxies. That is, the underlying star formation scaling law is always linear for clouds and galaxies with the same dense gas fraction. These considerations provide a single unified framework for understanding the relation between the standard (non-linear) extragalactic Schmidt-Kennicutt scaling law, that is typically derived from CO observations of the gas, and the linear star formation scaling law derived from HCN observations of the dense gas.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

Centroid positions of peaks identified in weak lensing mass maps often show offsets with respect to other means of identifying halo centres, like position of the brightest cluster galaxy or X-ray emission centroid. Here we study the effect of projected large-scale structure (LSS), smoothing of mass maps, and shape noise on the weak lensing peak positions. Additionally we compare the offsets in mass maps to those found in parametric model fits. Using ray-tracing simulations through the Millennium Run $N$-body simulation, we find that projected LSS does not alter the weak-lensing peak position within the limits of our simulations' spatial resolution, which exceeds the typical resolution of weak lensing maps. We conclude that projected LSS, although a major contaminant for weak-lensing mass estimates, is not a source of confusion for identifying halo centres. The typically reported offsets in the literature are caused by a combination of shape noise and smoothing alone. This is true for centroid positions derived both from mass maps and model fits.

We present a near-infrared extinction map of a large region (approximately 2200 deg^2) covering the Orion, the Monoceros R2, the Rosette, and the Canis Major molecular clouds. We used robust and optimal methods to map the dust column density in the near-infrared (NICER and NICEST) towards ~19 million stars of the Two Micron All Sky Survey (2MASS) point source catalog. Over the relevant regions of the field, we reached a 1-sigma error of 0.03 mag in the K-band extinction with a resolution of 3 arcmin. We measured the cloud distances by comparing the observed density of foreground stars with the prediction of galactic models, thus obtaining d_{Orion A} = (371 +/- 10) pc, d_{Orion B} = (398 +/- 12) pc, $d_{Mon R2} = (905 +/- 37) pc, $d_{Rosette} = (1330 +/- 48) pc, and $d_{CMa} = (1150 +/- 64) pc, values that compare very well with independent estimates.

Weak lensing is one of the best available diagnostic tools to measure the total density profiles of distant clusters of galaxies. Unfortunately, it suffers from the well-known mass-sheet degeneracy, so that weak lensing analyses cannot lead to fully reliable determinations of the total mass of the clusters. One possible way to set the relevant scale of the density profile would be to make a direct measurement of the magnification produced by the clusters as gravitational lenses. In this paper we revisit a suggestion made a few years ago for this general purpose, based on the use of the Fundamental Plane as a standard rod for early-type galaxies. Here we move one step further, beyond the simple outline of the idea given earlier, and quantify some statistical properties of this innovative diagnostic tool, with the final goal of identifying clear guidelines for a future observational test on concrete cases, which turns out to be well within the current instrument capabilities. The study is carried out by discussing the statistical properties of Fundamental Plane measurements for a sample of early-type source galaxies behind a massive cluster, for which a weak lensing analysis is assumed to be available. We proceed to study the best strategy to use Fundamental Plane measurements to determine the mass scale of a given cluster and find that the optimal choice is that of a sample of early-type galaxies behind the cluster distributed approximately uniformly in the sky. We also show that substructures do not contribute much to the magnification signal that we are looking for, but only add a modest amount of scatter. We find that for a massive cluster (M > 10^15 M_Sun) located at redshift 0.3 +/- 0.1, a set of about 20 Fundamental Plane measurements, combined with a good weak lensing analysis, should be able to lead to a mass determination with a precision of 20% or better.

A public release of slitless spectra, obtained with ACS/WFC and the G800L grism, is presented. Spectra were automatically extracted in a uniform way from 153 archival fields (or "associations") distributed across the two Galactic caps, covering all observations to 2008. The ACS G800L grism provides a wavelength range of 0.55-1.00 \mu$m, with a dispersion of $40 \ \AA / pixel$ and a resolution of $\sim 80\ \AA$ for point-like sources. The ACS G800L images and matched direct images were reduced with an automatic pipeline that handles all steps from archive retrieval, alignment and astrometric calibration, direct image combination, catalogue generation, spectral extraction and collection of metadata. The large number of extracted spectra (73,581) demanded automatic methods for quality control and an automated classification algorithm was trained on the visual inspection of several thousand spectra. The final sample of quality controlled spectra includes 47,919 datasets (65% of the total number of extracted spectra) for $32,149$ unique objects, with a median $i_{\rm AB}$-band magnitude of 23.7, reaching 26.5 AB for the faintest objects. Each released dataset contains science-ready 1D and 2D spectra, as well as multi-band image cutouts of corresponding sources and a useful preview page summarising the direct and slitless data, astrometric and photometric parameters. In order to characterize the slitless spectra, emission-line flux and equivalent width sensitivity of the ACS data were compared with public ground-based spectra in the GOODS-South field. An example list of emission line galaxies with two or more identified lines is also included, covering the redshift range $0.2-4.6$.

Lyman-break galaxies are now regularly found in the high redshift Universe by searching for the break in the galaxy spectrum caused by the Lyman-limit redshifted into the optical or even near-IR. At lower redshift, this break is covered by the GALEX UV channels and small samples of z ~ 1 LBGs have been presented in the literature. Here we give results from fitting the spectral energy distributions of a small sub-set of low redshift LBGs and demonstrate the advantage of including photometric points derived from HST ACS slitless grism observations. The results show these galaxies to have very young, star forming populations, while still being massive and dusty. LBGs at low and high redshift show remarkable similarities in their properties, indicating that the LBG selection method picks similar galaxies throughout the Universe.

We present a new set of high resolution dust extinction maps of the nearby and essentially starless Pipe Nebula molecular cloud. The maps were constructed from a concerted deep near-infrared imaging survey with the ESO-VLT, ESO-NTT, CAHA 3.5m telescopes, and 2MASS data. The new maps have a resolution three times higher than the previous extinction map of this cloud by Lombardi et al. (2006) and are able to resolve structure down to 2600 AU. We detect 244 significant extinction peaks across the cloud. These peaks have masses between 0.1 and 18.4 M_sun, diameters between 1.2 and 5.7e4 AU (0.06 and 0.28 pc), and mean densities of about 1e4 cm${^-3}$, all in good agreement with previous results. From the analysis of the Mean Surface Density of Companions we find a well defined scale near 1.4e4 AU below which we detect a significant decrease in structure of the cloud. This scale is smaller than the Jeans Length calculated from the mean density of the peaks. The surface density of peaks is not uniform but instead it displays clustering. Extinction peaks in the Pipe Nebula appear to have a spatial distribution similar to the stars in Taurus, suggesting that the spatial distribution of stars evolves directly from the primordial spatial distribution of high density material.

In this paper we investigate the level of star formation activity within nearby molecular clouds. We employ a uniform set of infrared extinction maps to provide accurate assessments of cloud mass and structure and compare these with inventories of young stellar objects within the clouds. We present evidence indicating that both the yield and rate of star formation can vary considerably in local clouds, independent of their mass and size. We find that the surface density structure of such clouds appears to be important in controlling both these factors. In particular, we find that the star formation rate (SFR) in molecular clouds is linearly proportional to the cloud mass (M_{0.8}) above an extinction threshold of A_K approximately equal to 0.8 magnitudes, corresponding to a gas surface density threshold of approximaely 116 solar masses per square pc. We argue that this surface density threshold corresponds to a gas volume density threshold which we estimate to be n(H_2) approximately equal to 10^4\cc. Specifically we find SFR (solar masses per yr) = 4.6 +/- 2.6 x 10^{-8} M_{0.8} (solar masses) for the clouds in our sample. This relation between the rate of star formation and the amount of dense gas in molecular clouds appears to be in excellent agreement with previous observations of both galactic and extragalactic star forming activity. It is likely the underlying physical relationship or empirical law that most directly connects star formation activity with interstellar gas over many spatial scales within and between individual galaxies. These results suggest that the key to obtaining a predictive understanding of the star formation rates in molecular clouds and galaxies is to understand those physical factors which give rise to the dense components of these clouds.

Larson (1981) first noted a scaling relation between masses and sizes in molecular clouds that implies that these objects have approximately constant column densities. This original claim, based upon millimeter observations of carbon monoxide lines, has been challenged by many theorists, arguing that the apparent constant column density observed is merely the result of the limited dynamic range of observations, and that in reality clouds have column density variations over two orders of magnitudes. In this letter we investigate a set of nearby molecular clouds with near-infrared excess methods, which guarantee very large dynamic ranges and robust column density measurements, to test the validity of Larson's third law. We verify that different clouds have almost identical average column densities above a given extinction threshold; this holds regardless of the extinction threshold, but the actual average surface mass density is a function of the specific threshold used. We show that a second version of Larson's third law, involving the mass-radius relation for single clouds and cores, does not hold in our sample, indicating that individual clouds are not objects that can be described by constant column density. Our results instead indicate that molecular clouds are characterized by a universal structure. Finally we point out that this universal structure can be linked to the log-normal nature of cloud column density distributions.

We present a lensing and photometric study of the exceptional system SDSS J1538+5817, identified by the SLACS survey. The lens is a luminous elliptical at redshift z=0.143. Using HST public images in two different filters, the presence of two background sources lensed into an Einstein ring and a double system is ascertained. Our new spectroscopic observations, performed at the NOT, reveal that the two sources are located at the same redshift z=0.531. We investigate the total mass distribution of the lens between 1 and 4 kpc from the galaxy center by means of parametric and non-parametric lensing codes that describe the multiple images as point-like objects. Several disparate lensing models agree on: (1) reproducing accurately the observed image positions; (2) predicting a nearly axisymmetric total mass distribution, centered and oriented as the light distribution; (3) measuring a value of 8.11 x 10^{10} M_{Sun} for the total mass projected within the Einstein radius of 2.5 kpc; (4) estimating a total mass density profile slightly steeper than an isothermal one. A fit of the SDSS multicolor photometry with CSP models provides a value of 20 x 10^{10} M_{Sun} for the total stellar mass of the galaxy and of 0.9 for the fraction of projected luminous over total mass enclosed inside the Einstein radius. By combining lensing and photometric mass measurements, we differentiate the lens mass content in terms of luminous and dark matter components. This two-component modeling, which is viable only in extraordinary systems like SDSS J1538+5817, leads to a description of the global properties of the galaxy dark matter halo. Extending these results to a larger number of lenses would improve considerably our understanding of galaxy formation and evolution processes in the LCDM scenario.

The Pipe Nebula, a large nearby molecular cloud lacks obvious signposts of star formation in all but one of more than 130 dust extinction cores that have been identified within it. In order to quantitatively determine the current level of star formation activity in the Pipe Nebula, we analyzed 13 square degrees of sensitive mid-infrared maps of the entire cloud, obtained with the Multiband Imaging Photometer for Spitzer (MIPS) at wavelengths of 24 micron and 70 micron to search for candidate Young Stellar Objects (YSOs) in the high-extinction regions. We argue that our search is complete for class I and typical class II YSOs with luminosities of Lbol~0.2 Lo and greater. We find only 18 candidate YSOs in the high-extinction regions of the entire Pipe cloud. Twelve of these sources are previously known members of a small cluster associated with Barnard 59, the largest and most massive dense core in the cloud. With only six candidate class I and class II YSOs detected towards extinction cores outside of this cluster, our findings emphatically confirm the notion of an extremely low level of star formation activity in the Pipe Nebula. The resulting star formation efficiency for the entire cloud mass is only ~0.06 %.

We present an analysis of wide-field infrared extinction maps of a region in Perseus just north of the Taurus-Auriga dark cloud complex. From this analysis we have identified a massive, nearby, but previously unrecognized, giant molecular cloud (GMC). From comparison of foreground star counts with Galactic models we derive a distance of 450 +/- 23 parsecs to the cloud. At this distance the cloud extends over roughly 80 pc and has a mass of approximately 10^5 solar masses, rivaling the Orion (A) Molecular Cloud as the largest and most massive GMC in the solar neighborhood. Although surprisingly similar in mass and size to the more famous Orion Molecular Cloud (OMC) the newly recognized cloud displays significantly less star formation activity with more than an order of magnitude fewer young stellar objects than found in the OMC, suggesting that both the level of star formation and perhaps the star formation rate in this cloud are an order of magnitude or more lower than in the OMC. Analysis of extinction maps of both clouds shows that the new cloud contains only 10% the amount of high extinction (A_K > 1.0 mag) material as is found in the OMC. This, in turn, suggests that the level of star formation activity and perhaps the star formation rate in these two clouds may be directly proportional to the total amount of high extinction material and presumably high density gas within them and that there might be a density threshold for star formation on the order of a few times 10^4 hydrogen molecules per cc.

We seek to understand the impact on shape estimators obtained from circular and elliptical shapelet models under two realistic conditions: (a) only a limited number of shapelet modes is available for the model, and (b) the intrinsic galactic shapes are not restricted to shapelet models. We create a set of simplistic simulations, in which the galactic shapes follow a Sersic profile. By varying the Sersic index and applied shear, we quantify the amount of bias on shear estimates which arises from insufficient modeling. Additional complications due to PSF convolution, pixelation and pixel noise are also discussed. Steep and highly elliptical galaxy shapes cannot be accurately modeled within the circular shapelet basis system and are biased towards shallower and less elongated shapes. This problem can be cured partially by allowing elliptical basis functions, but for steep profiles elliptical shapelet models still depend critically on accurate ellipticity priors. As a result, shear estimates are typically biased low. Independently of the particular form of the estimator, the bias depends on the true intrinsic galaxy morphology, but also on the size and shape of the PSF. As long as the issues discussed here are not solved, the shapelet method cannot provide weak-lensing measurements with an accuracy demanded by upcoming missions and surveys, unless one can provide an accurate and reliable calibration, specific for the dataset under investigation.

By using the SDSS multicolor photometry and lens modeling, we study stellar mass properties and the luminous and dark matter composition of the 57 early-type lens galaxies analyzed by the SLACS Survey. We fit the lens SEDs composed of ugriz magnitudes with a three-parameter grid of different composite stellar population models, computed by adopting solar metallicity and various IMFs, and we employ the best-fit values of the total projected mass enclosed within the Einstein ring of each system. We measure size-stellar mass and surface stellar mass density-stellar mass relations consistent with those determined for non-lens galaxies. We find statistically significant evidence that more massive lens galaxies tend to form in regions of higher galaxy density, as for all early-type galaxies. We measure for the corresponding stellar quantities the same scaling law between effective mass-to-light ratio and mass, that is used to explain the "tilt" of the FP, and the same evolution of the effective mass-to-light ratio with redshift, that can be derived from the FP. We conclude that the lens total mass is linearly proportional to the luminous mass, at a more than 99% CL. In addition, by assuming that the lens galaxies are homologous systems, we study their distribution of dark matter and estimate a value on the order of 30% for the dark over total projected mass fraction contained within the average Einstein radius of approximately 4 kpc. These results suggest that early-type lens galaxies are an unbiased subsample representative of the whole sample of early-type galaxies. This allows us to generalize our findings on dark matter in lens galaxies to the population of massive early-type galaxies. We note that, for the assumed metallicity, a Salpeter IMF describes well the lens sample.

In this paper we derive an improved core mass function (CMF) for the Pipe Nebula from a detailed comparison between measurements of visual extinction and molecular-line emission. We have compiled a refined sample of 201 dense cores toward the Pipe Nebula using a 2-dimensional threshold identification algorithm informed by recent simulations of dense core populations. Measurements of radial velocities using complimentary C18O (1-0) observations enable us to cull out from this sample those 43 extinction peaks that are either not associated with dense gas or are not physically associated with the Pipe Nebula. Moreover, we use the derived C18O, central velocities to differentiate between single cores with internal structure and blends of two or more physically distinct cores, superposed along the same line-of-sight. We then are able to produce a more robust dense core sample for future follow-up studies and a more reliable CMF than was possible previously. We confirm earlier indications that the CMF for the Pipe Nebula departs from a single power-law like form with a break or knee at M ~ 2.7 +/- 1.3 Msun. Moreover, we also confirm that the CMF exhibits a similar shape to the stellar IMF, but is scaled to higher masses by a factor of ~4.5. We interpret this difference in scaling to be a measure of the star formation efficiency (22 +/- 8%). This supports earlier suggestions that the stellar IMF may originate more or less directly from the CMF.

Optically selected clusters of galaxies display a relation between their optical mass estimates and their X-ray luminosities Lx that has a large scatter. A substantial fraction of optically selected clusters have Lx estimates or upper limits significantly below the values expected from the Lx-mass relation established for X-ray selected clusters, i.e., these clusters are X-ray underluminous for their mass. We attempt to confirm or falsify the X-ray underluminous nature of two clusters, Abell 315 and Abell 1456, by using weak gravitational lensing as a third and independent measure of the clusters' masses. We obtained optical wide-field imaging data and selected background galaxies using their colors and measured the shear exerted by the tidal field of the foreground galaxy clusters. We then fitted parametrized models to our shear catalogs. After accounting for projections of large-scale structure and halo triaxiality, we find that A 315 is significantly X-ray underluminous for its mass, while no significant lensing signal was detected for A 1456. We re-evaluate earlier kinematic and X-ray analyses of these two clusters and discuss the nature of the X-ray underluminous cluster A 315 and why A 1456 was probably erroneously identified as being X-ray underluminous.

We present an extinction map of a ~1,700 deg sq region that encloses the Ophiuchus, the Lupus, and the Pipe dark complexes using 42 million stars from the Two Micron All Sky Survey (2MASS) point source catalog. The use of a robust and optimal near-infrared method to map dust column density (Nicer, described in Lombardi & Alves 2001) allow us to detect extinction as low as A_K = 0.05 mag with a 2-sigma significance, and still to have a resolution of 3 arcmin on our map. We also present a novel, statistically sound method to characterize the small-scale inhomogeneities in molecular clouds. Finally, we investigate the cloud structure function, and show that significant deviations from the results predicted by turbulent models are observed.

Observational data and theoretical calculations show that significant small-scale substructures are present in dark molecular clouds. These inhomogeneities can provide precious hints on the physical conditions inside the clouds, but can also severely bias extinction measurements. We present NICEST, a novel method to account and correct for inhomogeneities in molecular cloud extinction studies. The method, tested against numerical simulations, removes almost completely the biases introduced by sub-pixel structures and by the contamination of foreground stars. We applied NICEST to 2MASS data of the Pipe molecular complex. The map thereby obtained shows significantly higher (up to 0.41 mag in A_K) extinction peaks than the standard NICER (Lombardi et al. 2001) map. This first application confirms that substructures in nearby molecular clouds, if not accounted for, can significantly bias extinction measurements in regions with A_K > 1 mag; the effect, moreover, is expected to increase in more distant molecular cloud, because of the poorer physical resolution achievable.

The Hubble Legacy Archive (HLA) aims to create calibrated science data from the Hubble Space Telescope archive and make them accessible via user-friendly and Virtual Observatory (VO) compatible interfaces. It is a collaboration between the Space Telescope Science Institute (STScI), the Canadian Astronomy Data Centre (CADC) and the Space Telescope - European Coordinating Facility (ST-ECF). Data produced by the Hubble Space Telescope (HST) instruments with slitless spectroscopy modes are among the most difficult to extract and exploit. As part of the HLA project, the ST-ECF aims to provide calibrated spectra for objects observed with these HST slitless modes. In this paper, we present the HLA NICMOS G141 grism spectra. We describe in detail the calibration, data reduction and spectrum extraction methods used to produce the extracted spectra. The quality of the extracted spectra and associated direct images is demonstrated through comparison with near-IR imaging catalogues and existing near-IR spectroscopy. The output data products and their associated metadata are publicly available through a web form at http://hla.stecf.org and via VO interfaces. In total, 2470 spectra of 1923 unique targets are included in the current release.

Gravitational lens modeling is presented for the first discovered example of a three-component source for which each component is quadruply imaged. The lens is a massive galaxy member of the cluster Cl J0152.7-1357 at z ~ 0.84. Taking advantage of this exceptional configuration and of the excellent angular resolution of the HST/ACS, we measure the properties of the lens. Several parametric macroscopic models were developed for the lens galaxy, starting from pointlike to extended image models. For a lens model in terms of a singular isothermal sphere with external shear, the Einstein radius is found to be R_{E} = (9.54 +/- 0.15) kpc. The external shear points to the cluster's northern mass peak. The unknown redshift of the source is determined to be higher than 1.9 and lower than 2.9. Our estimate of the lensing projected total mass inside the Einstein radius, M_{len}(R < 9.54 kpc), depends on the source distance and lies between 4.6 and 6.2 x 10^{11} M_{Sun}. This result turns out to be compatible with the dynamical estimate based on an isothermal model. By considering the constraint on the stellar mass-to-light ratio that comes from the evolution of the Fundamental Plane, we can exclude the possibility that at more than 4 sigma level the total mass enclosed inside the Einstein ring is only luminous matter. Moreover, the photometric-stellar mass measurement within the Einstein radius gives a minimum value of 50% (1 sigma) for the dark-to-total matter fraction. The lensing analysis has allowed us to investigate the distribution of mass of the deflector, also providing some interesting indications on scales that are larger and smaller than the Einstein radius of the lens galaxy. The combination of different diagnostics has proved to be essential in determining quantities that otherwise would have not been directly measurable with the current data.

The galaxy cluster RX J1347-1145 is one of the most X-ray luminous and most massive clusters known. Its extreme mass makes it a prime target for studying issues addressing cluster formation and cosmology. In this paper we present new high-resolution HST/ACS and Chandra X-ray data. The high resolution and sensitivity of ACS enabled us to detect and quantify several new multiply imaged sources, we now use a total of eight for the strong lensing analysis. Combining this information with shape measurements of weak lensing sources in the central regions of the cluster, we derive a high-resolution, absolutely-calibrated mass map. This map provides the best available quantification of the total mass of the central part of the cluster to date. We compare the reconstructed mass with that inferred from the new Chandra X-ray data, and conclude that both mass estimates agree extremely well in the observed region, namely within 400 / h_70 kpc of the cluster center. In addition we study the major baryonic components (gas and stars) and hence derive the dark matter distribution in the center of the cluster. We find that the dark matter and baryons are both centered on the BCG within the uncertainties (alignment is better than <10 kpc). We measure the corresponding 1-D profiles and find that dark matter distribution is consistent with both NFW and cored profiles, indicating that a more extended radial analysis is needed to pinpoint the concentration parameter, and hence the inner slope of the dark matter profile.

We combine extinction maps from the Two Micron All Sky Survey (2MASS) with Hipparcos and Tycho parallaxes to obtain reliable and high-precision estimates of the distance to the Ophiuchus and Lupus dark complexes. Our analysis, based on a rigorous maximum-likelihood approach, shows that the rho-Ophiuchi cloud is located at (119 +/- 6) pc and the Lupus complex is located at (155 +/- 8) pc; in addition, we are able to put constraints on the thickness of the clouds and on their orientation on the sky (both these effects are not included in the error estimate quoted above). For Ophiuchus, we find some evidence that the streamers are closer to us than the core. The method applied in this paper is currently limited to nearby molecular clouds, but it will find many natural applications in the GAIA-era, when it will be possible to pin down the distance and three-dimensional structure of virtually every molecular cloud in the Galaxy.

High resolution HST/ACS images of the galaxy cluster RX J1347-1145 have enabled us to identify several new multiple image candidates in the cluster, including a 5 image system with a central image. The multiple images allow us to construct an accurate 2-dimensional mass map of the central part of the cluster. The modelling of the cluster mass includes the most prominent cluster galaxies modelled as truncated isothermal spheres and a smooth halo component that is described with 2 parametric profiles. The mass reconstruction is done using a Markov chain Monte Carlo method that provides us with a total projected mass density as well as estimates for the parameters of interest and their respective errors. The mass profile is in reasonable agreement with previous mass estimates based on the X-ray emission from the hot intra-cluster gas, however the X-ray mass estimates are systematically lower than what we obtain with gravitational lensing.

In this paper we compare two different diagnostics for estimating stellar masses in early-type galaxies and we establish their level of reliability. In particular, we consider the well-studied sample of 15 field elliptical galaxies selected from the Sloan Lens ACS (SLACS) Survey (z = 0.06-0.33). We examine here the correlation between the stellar mass values, enclosed inside the Einstein radius of each lens, based on analyses of lensing and stellar dynamics combined and based on multiwavelength photometry spectral template fitting. The lensing+dynamics stellar mass is obtained from the published SLACS Survey results, assuming a two-component density distribution model and a prior from the fundamental plane on the mass-to-light ratio for the lens galaxies. The photometric stellar mass is measured by fitting the observed spectral energy distribution of the galaxies (from the SDSS multi-band photometry over 354-913 nm) with composite stellar population templates, under the assumption that light traces stellar mass. The two methods are completely independent. They rely on several different assumptions, and so, in principle, both can have significant biases. Based on our sample of massive galaxies, we find consistency between the lensing+dynamics and the photometric mass estimates. We obtain a Pearson linear correlation coefficient of 0.94 and a median value of the ratio between the former and the latter mass measurements of 1.1+/-0.1. This suggests that both methods can separately yield reliable stellar masses of early-type galaxies, and confirms that photometric mass estimates are accurate, as long as optical/near-IR rest frame photometry is available.

We show how the combination of observations related to strong gravitational lensing and stellar dynamics in ellipticals offers a new way to measure the cosmological matter and dark-energy density parameters. A gravitational lensing estimate of the mass enclosed inside the Einstein circle can be obtained by measuring the Einstein angle, once the critical density of the system is known. A model-dependent dynamical estimate of this mass can also be obtained by measuring the central velocity dispersion of the stellar component. By assuming the well-tested homologous 1/r^{2} profile for the total density distribution in the lens elliptical galaxies, these two mass measurements can be properly compared. Thus, a relation between the Einstein angle and the central stellar velocity dispersion is derived, and the cosmological matter and the dark-energy density parameters can be estimated from this. We determined the accuracy of the cosmological parameter estimates by means of simulations that include realistic measurement uncertainties on the relevant quantities. Interestingly, the expected constraints on the cosmological parameter plane are complementary to those coming from other observational techniques. Then, we applied the method to the data sets of the Sloan Lens ACS and the Lenses Structure and Dynamics Surveys, and showed that the concordance value between 0.7 and 0.8 for the dark-energy density parameter is included in our 99% confidence regions. The small number of lenses available to date prevents us from precisely determining the cosmological parameters, but it still proves the feasibility of the method. When applied to samples made of hundreds of lenses that are expected to become available from forthcoming surveys, this technique will be an important tool for measuring the geometry of the Universe.

In this paper we present the results of a systematic investigation of an entire population of starless dust cores within a single molecular cloud. Analysis of extinction data shows the cores to be dense objects characterized by a narrow range of density. Analysis of C18O and NH3 molecular-line observations reveals very narrow lines. The non-thermal velocity dispersions measured in both these tracers are found to be subsonic for the large majority of the cores and show no correlation with core mass (or size). Thermal pressure is thus the dominate source of internal gas pressure and support for most of the core population. The total internal gas pressures of the cores are found to be roughly independent of core mass over the entire range of the core mass function (CMF) indicating that the cores are in pressure equilibrium with an external source of pressure. This external pressure is most likely provided by the weight of the surrounding Pipe cloud within which the cores are embedded. Most of the cores appear to be pressure confined, gravitationally unbound entities whose nature, structure and future evolution are determined by only a few physical factors which include self-gravity, the fundamental processes of thermal physics and the simple requirement of pressure equilibrium with the surrounding environment. The observed core properties likely constitute the initial conditions for star formation in dense gas. The entire core population is found to be characterized by a single critical Bonnor-Ebert mass. This mass coincides with the characteristic mass of the Pipe CMF indicating that most cores formed in the cloud are near critical stability. This suggests that the mass function of cores (and the IMF) has its origin in the physical process of thermal fragmentation in a pressurized medium.

We present molecular-line observations of 94 dark cloud cores identified in the Pipe nebula through near-IR extinction mapping. Using the Arizona Radio Observatory 12m telescope, we obtained spectra of these cores in the J=1-0 transition of C18O. We use the measured core parameters, i.e., antenna temperature, linewidth, radial velocity, radius and mass, to explore the internal kinematics of these cores as well as their radial motions through the larger molecular cloud. We find that the vast majority of the dark extinction cores are true cloud cores rather than the superposition of unrelated filaments. While we identify no significant correlations between the core's internal gas motions and the cores' other physical parameters, we identify spatially correlated radial velocity variations that outline two main kinematic components of the cloud. The largest is a 15pc long filament that is surprisingly narrow both in spatial dimensions and in radial velocity. Beginning in the Stem of the Pipe, this filament displays uniformly small C18O linewidths (dv~0.4kms-1) as well as core to core motions only slightly in excess of the gas sound speed. The second component outlines what appears to be part of a large (2pc; 1000 solar mass) ring-like structure. Cores associated with this component display both larger linewidths and core to core motions than in the main cloud. The Pipe Molecular Ring may represent a primordial structure related to the formation of this cloud.

Recent extinction studies of the Pipe Nebula (d=130 pc) reveal many cores spanning a range in mass from 0.2 to 20.4 Msun. These dense cores were identified via their high extinction and comprise a starless population in a very early stage of development. Here we present a survey of NH3 (1,1), NH3 (2,2), CCS (2_1,1_0), and HC5N (9,8) emission toward 46 of these cores. An atlas of the 2MASS extinction maps is also presented. In total, we detect 63% of the cores in NH3 (1,1) 22% in NH3 (2,2), 28% in CCS, and 9% in HC5N emission. We find the cores are associated with dense gas (~10^4 cm-3) with 9.5 < T_k < 17 K. Compared to C18O, we find the NH3 linewidths are systematically narrower, implying that the NH3 is tracing the dense component of the gas and that these cores are relatively quiescent. We find no correlation between core linewidth and size. The derived properties of the Pipe cores are similar to cores within other low-mass star-forming regions: the only differences are that the Pipe cores have weaker NH3 emision and most show no current star formation as evidenced by the lack of embedded infrared sources. Such weak NH3 emission could arise due to low column densities and abundances or reduced excitation due to relatively low core volume densities. Either alternative implies that the cores are relatively young. Thus, the Pipe cores represent an excellent sample of dense cores in which to study the initial conditions for star formation and the earliest stages of core formation and evolution.

We present a detailed mass reconstruction of the cluster RCS0224-0002 at $z=0.773$ from the strong lensing features observed with HST/WFPC2. The mass profile is reconstructed using a parametric approach. We introduce a novel method to fit extended multiple images based on the Modified Hausdorff Distance between observed arcs and the arcs reproduced by the model. We perform the detailed error analysis of the model parameter using the MCMC method. Our model reproduces all the observed strong lensing features of the RCS0224-0002 and predicts the redshift of one of the arcs systems to be $z\approx 2.65$ (the other system has an spectroscopic redshift of $z=4.87$). The reconstructed inner mass profile is well fitted by a non-singular isothermal sphere, rather than with an NFW model. Dark matter substructure, derived from the light distribution of the most luminous cluster members, is crucial for reproducing the complexity of the quadrupole image system, which could not be achieved otherwise. The reconstructed mass distribution closely follows the light, however it is significantly shifted from the X-ray emission of the gas. The mass of RCS0224-0002 derived from the lensing model, $\approx 2\times10^{14} M_\odot$ is in a very good agreement with the one obtained from the X-ray temperature measured with deep Chandra observations.

Context: Stars form in the cold dense cores of interstellar molecular clouds and the detailed knowledge of the spectrum of masses of such cores is clearly a key for the understanding of the origin of the IMF. To date, observations have presented somewhat contradictory evidence relating to this issue. Aims: In this paper we propose to derive the mass function of a complete sample of dense molecular cores in a single cloud employing a robust method that uses uses extinction of background starlight to measure core masses and enables the reliable extension of such measurements to lower masses than previously possible. Methods: We use a map of near-infrared extinction in the nearby Pipe dark cloud to identify the population of dense cores in the cloud and measure their masses. Results: We identify 159 dense cores and construct the mass function for this population. We present the first robust evidence for a departure from a single power-law form in the mass function of a population of cores and find that this mass function is surprisingly similar in shape to the stellar IMF but scaled to a higher mass by a factor of about 3. This suggests that the distribution of stellar birth masses (IMF) is the direct product of the dense core mass function and a uniform star formation efficiency of 30%+/-10%, and that the stellar IMF may already be fixed during or before the earliest stages of core evolution. These results are consistent with previous dust continuum studies which suggested that the IMF directly originates from the core mass function. The typical density of ~10^4/cm^3 measured for the dense cores in this cloud suggests that the mass scale that characterizes the dense core mass function may be the result of a simple process of thermal (Jeans) fragmentation.

We consider the Fundamental Plane of elliptical galaxies lensed by the gravitational field of a massive deflector (typically, a cluster of galaxies). We show that the Fundamental Plane relation provides a straightforward measurement of the projected mass distribution of the lens with a typical accuracy of ~0.15 in the dimensionless column density kappa. The proposed technique breaks the mass-sheet degeneracy completely and is thus expected to serve as an important complement to other lensing-based analyses. Moreover, its ability to measure directly the mass distribution on the small pencil beams that characterize the size of background galaxies may lead to crucial tests for current scenarios of structure formation.

We present a 8 deg x 6 deg, high resolution extinction map of the Pipe nebula using 4.5 million stars from the Two Micron All Sky Survey (2MASS) point source catalog. The use of NICER, a robust and optimal technique to map the dust column density, allows us to detect a Av = 0.5 mag extinction at a 3-sigma level with a 1 arcmin resolution. We find for the Pipe nebula a normal reddening law, E(J-H) = (1.85 +/- 0.15) E(H-K). We measure the cloud distance using Hipparchos and Tycho parallaxes, and obtain ~130 pc. This, together with the total estimated mass, 10^4 Msun, makes the Pipe the closest massive cloud complex to Earth. We compare the NICER extinction map to the NANTEN 12CO observations and derive with unprecedented accuracy the relationship between the near-infrared extinction and the 12CO column density and hence (indirectly) the 12CO X-factor, that we estimate to be 2.91 10^20 cm^-2 K^-1 km^-1 s in the range Av <- [0.9, 5.4] mag. We identify ~1500 OH/IR stars located within the Galactic bulge in the direction of the Pipe field. This represents a significant increase of the known numbers of such stars in the Galaxy. Our analysis confirms the power and simplicity of the color excess technique to study molecular clouds. The comparison with the NANTEN 12CO data corroborates the insensitivity of CO observations to low column densities (up to approximately 2 mag in Av), and shows also an irreducible uncertainty in the dust-CO correlation of about 1 mag of visual extinction.

We present an overview of data available for the Ophiuchus and Perseus molecular clouds from ``Phase I'' of the COMPLETE Survey of Star-Forming Regions. This survey provides a range of data complementary to the Spitzer Legacy Program ``From Molecular Cores to Planet Forming Disks.'' Phase I includes: Extinction maps derived from 2MASS near-infrared data using the NICER algorithm; extinction and temperature maps derived from IRAS 60 and 100um emission; HI maps of atomic gas; 12CO and 13CO maps of molecular gas; and submillimetre continuum images of emission from dust in dense cores. Not unexpectedly, the morphology of the regions appears quite different depending on the column-density tracer which is used, with IRAS tracing mainly warmer dust and CO being biased by chemical, excitation and optical depth effects. Histograms of column-density distribution are presented, showing that extinction as derived from 2MASS/NICER gives the closest match to a log-normal distribution as is predicted by numerical simulations. All the data presented in this paper, and links to more detailed publications on their implications are publically available at the COMPLETE website.

We consider from a general point of view the problem of determining the extinction in dense molecular clouds. We use a rigorous statistical approach to characterize the properties of the most widely used optical and infrared techniques, namely the star count and the color excess methods. We propose a new maximum-likelihood method that takes advantage of both star counts and star colors to provide an optimal estimate of the extinction. Detailed numerical simulations show that our method performs optimally under a wide range of conditions and, in particular, is significantly superior to the standard techniques for clouds with high column-densities and affected by contamination by foreground stars.

We have shown that the cluster-mass reconstruction method which combines strong and weak gravitational lensing data, developed in the first paper in the series, successfully reconstructs the mass distribution of a simulated cluster. In this paper we apply the method to the ground-based high-quality multi-colour data of RX J1347.5-1145, the most X-ray luminous cluster to date. A new analysis of the cluster core on very deep, multi-colour data analysis of VLT/FORS data reveals many more arc candidates than previously known for this cluster. The combined strong and weak lensing reconstruction confirms that the cluster is indeed very massive. If the redshift and identification of the multiple-image system as well as the redshift estimates of the source galaxies used for weak lensing are correct, we determine the enclosed cluster mass in a cylinder to M(<360 h^-1 kpc)= (1.2 +/- 0.3) 10^15 Msun. In addition the reconstructed mass distribution follows the distribution found with independent methods (X-ray measurements, SZ). With higher resolution (e.g. HST imaging data) more reliable multiple imaging information can be obtained and the reconstruction can be improved to accuracies greater than what is currently possible with weak and strong lensing techniques.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, the mass-sheet degeneracy transformation has limited its success. We present a novel method for a cluster mass reconstruction which combines weak and strong lensing information on common scales and can, as a consequence, break the mass-sheet degeneracy. We extend the weak lensing formalism to the inner parts of the cluster and combine it with the constraints from multiple image systems. We demonstrate the feasibility of the method with simulations, finding an excellent agreement between the input and reconstructed mass also on scales within and beyond the Einstein radius. Using a single multiple image system and photometric redshift information of the background sources used for weak and strong lensing analysis, we find that we are effectively able to break the mass-sheet degeneracy, therefore removing one of the main limitations on cluster mass estimates. We conclude that with high resolution (e.g. HST) imaging data the method can more accurately reconstruct cluster masses and their profiles than currently existing lensing techniques.

We study the relation of the three-point cosmic shear statistics to the third-order statistical properties of the underlying convergence, expressed in terms of its bispectrum. Explicit relations for the natural components of the shear three-point correlation function in terms of the bispectrum are derived. The behavior of the correlation function under parity transformation is obtained and found to agree with previous results. We find that in contrast to the two-point shear correlation function, the three-point function at a given angular scale \theta is not affected by power in the bispectrum on much larger scales. These relations are then inverted to obtain the bispectrum in terms of the three-point shear correlator; two different expressions, corresponding to different natural components of the shear correlator, are obtained and can be used to separate E and B-mode shear contributions. These relations allow us to explicitly show that correlations containing an odd power of B-mode shear vanish for parity-symmetric fields. Generalizing a recent result by Jarvis et al., we derive expressions for the third-order aperture measures, employing multiple angular scales, in terms of the (natural components of the) three-point shear correlator and show that they contain essentially all the information about the underlying bispectrum. We discuss the many useful features these (generalized) aperture measures have that makes them convenient for future analyses of the skewness of the cosmic shear field (and any other polar field, such as the polarization of the Cosmic Microwave Background). (Abridged)

We present a weak lensing analysis of one of the most distant massive galaxy cluster known, RDCS 1252.9-2927 at z=1.24, using deep images from the Advanced Camera for Survey (ACS) on board the Hubble Space Telescope (HST). By taking advantage of the depth and of the angular resolution of the ACS images, we detect for the first time at z>1 a clear weak lensing signal in both the i (F775W) and z (F850LP) filters. We measure a 5-\sigma signal in the i band and a 3-\sigma signal in the shallower z band image. The two radial mass profiles are found to be in very good agreement with each other, and provide a measurement of the total mass of the cluster inside a 1Mpc radius of M(<1Mpc) = (8.0 +/- 1.3) x 10^14 M_\odot in the current cosmological concordance model h =0.70, \Omega_m=0.3, \Omega_\Lambda=0.7, assuming a redshift distribution of background galaxies as inferred from the Hubble Deep Fields surveys. A weak lensing signal is detected out to the boundary of our field (3' radius, corresponding to 1.5Mpc at the cluster redshift). We detect a small offset between the centroid of the weak lensing mass map and the brightest cluster galaxy, and we discuss the possible origin of this discrepancy. The cumulative weak lensing radial mass profile is found to be in good agreement with the X-ray mass estimate based on Chandr and XMM-Newton observations, at least out to R_500=0.5Mpc.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, weak lensing mass reconstructions are plagued by the so-called mass-sheet degeneracy--the surface mass density \kappa of the cluster can be determined only up to a degeneracy transformation \kappa \to \kappa' = \lambda \kappa + (1 -\lambda), where \lambda is an arbitrary constant. This transformation fundamentally limits the accuracy of cluster mass determinations if no further assumptions are made. We discuss here a possibility to break the mass-sheet degeneracy in weak lensing mass maps using distortion and redshift information of background galaxies. Compared to other techniques proposed in the past, it does not rely on any assumptions on cluster potential and does not make use of weakly constrained information (such as the source number counts, used in the magnification effect).Our simulations show that we are effectively able to break the mass-sheet degeneracy for supercritical lenses and that for undercritical lenses the mass-sheet degeneracy is very difficult to be broken, even under idealised conditions.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, weak lensing mass reconstructions are plagued by the so-called mass-sheet degeneracy--the surface mass density \kappa of the cluster can be determined only up to a degeneracy transformation \kappa \to \kappa' = \lambda \kappa + (1 -\lambda), where \lambda is an arbitrary constant. This transformation fundamentally limits the accuracy of cluster mass determinations if no further assumptions are made. We describe here a method to break the mass-sheet degeneracy in weak lensing mass maps using distortion and redshift information of background galaxies and illustrate this by two simple toy models. Compared to other techniques proposed in the past, it does not rely on any assumptions on cluster potential; it can be easily applied to non-parametric mass-reconstructions and no assumptions on boundary conditions have to be made. In addition it does not make use of weakly constrained information (such as the source number counts, used in the magnification effect). Our simulations show that we are effectively able to break the mass-sheet degeneracy for supercritical lenses, but that for undercritical lenses the mass-sheet degeneracy is very difficult to be broken, even under idealised conditions.

We present a study of the lens properties of quadruply imaged systems, lensed by numerically simulated galaxies. We investigate a simulated elliptical and disc galaxy drawn from high resolution simulations of galaxy formation in a concordance $\Lambda$CDM universe. The simulations include the effects of gas dynamics, star formation and feedback processes. Flux-ratio anomalies observed in strong gravitational lensing potentially provide an indicator for the presence of mass substructure in lens galaxies as predicted from CDM simulations. We particularly concentrate on the prediction that, for an ideal cusp caustic, the sum of the signed magnifications of the three highly magnified images should vanish when the source approaches the cusp. Strong violation of this cusp relation indicates the presence of substructure, regardless of the global, smooth mass model of the lens galaxy. We draw the following conclusions: (1) the level of substructure present in simulations produces violations of the cusp relation comparable to those observed, (2) higher-order catastrophes (e.g. swallowtails) can also cause changes of the order of 0.6 in the cusp relation as predicted by a smooth model, (3) the flux anomaly distribution depends on the image parity and flux and both the brightest minimum and saddle-point images are more affected by substructure than the fainter images. In addition, the brightest saddle point is demagnified w.r.t. the brightest minimum. Our results are fully numerical and properly include all mass scales, without making semi-analytic assumptions. They are ultimately limited by the mass resolution of single particles in the simulation determined by current computational limits, however show that our results are not affected by shot-noise due to the finite number of particles.

The Lynx arc, with a redshift of 3.357, was discovered during spectroscopic follow-up of the z=0.570 cluster RX J0848+4456 from the ROSAT Deep Cluster Survey. The arc is characterized by a very red R-K color and strong, narrow emission lines. Analysis of HST WFPC2 imaging and Keck optical and infrared spectroscopy shows that the arc is an \hii galaxy magnified by a factor of ~10 by a complex cluster environment. The high intrinsic luminosity, the emission line spectrum, the absorption components seen in Ly\alpha and C IV, and the restframe ultraviolet continuum are all consistent with a simple HII region model containing ~ 10^6 hot O stars. The best fit parameters for this model imply a very hot ionizing continuum (T_BB ~ 80,000 K), high ionization parameter (log U ~ -1), and low nebular metallicity (Z / Z_\odot ~ 0.05). The narrowness of the emission lines requires a low mass-to-light ratio for the ionizing stars, suggestive of an extremely low metallicity stellar cluster. The apparent overabundance of silicon in the nebula could indicate enrichment by past pair instability supernov\ae, requiring stars more massive than ~140 M_\odot.

In a series of papers (Lombardi & Schneider 2001, 2002) we studied in detail the statistical properties of an interpolation technique widely used in astronomy. In particular, we considered the average interpolated map and its covariance under the hypotheses that the map is obtained by smoothing unbiased measurements of an unknown field, and that the measurements are uniformly distributed on the sky. In this paper we generalize the results obtained to the case of observations carried out only on a finite field and distributed on the field with a non-uniform density. These generalizations, which are required in many astronomically relevant cases, still allow an exact, analytical solution of the problem. We also consider a number of properties of the interpolated map, and provide asymptotic expressions for the average map and the two-point correlation function which are valid at high densities.

We analyse the Chandra dataset of the galaxy cluster MS1008.1-1224 to recover an estimate of the gravitating mass as function of the radius and compare these results with the weak lensing reconstruction of the mass distribution obtained from deep FORS1-VLT multicolor imaging. Even though the X-ray morphology is disturbed with a significant excess in the northern direction suggesting that the cluster is not in a relaxed state, we are able to match the two mass profiles both in absolute value and in shape within 1 sigma uncertainty and up to 1100 h50^-1 kpc. The recovered X-ray mass estimate does not change by using either the azimuthally averaged gas density and temperature profiles or the results obtained in the northern sector alone where the signal-to-noise ratio is higher.

Previous work on three-point statistics of cosmic shear has mainly concentrated on the convergence, or on aperture measures of the shear. However, as has become clear recently for the two-point statistics of cosmic shear, the basic quantity that should be used is the correlation function: first, it is much easier to measure from observational data, since it is immune against complicated geometries of data fields; second, all other (linear) 2-pt statistics can be expressed as integrals over the correlation function. The situation is the same for the 3-pt statistics. However, in contrast to the 2-pt correlation function, the invariants of the shear 3-pt correlation function have not been employed yet. Here we consider the transformation properties of the shear 3-pt correlation function under rotations. We show that there are four complex linear combinations of components of the 3-pt correlation function, which we shall call `natural components', since they are multiplied just by a phase factor for arbitrary rotations, but do not mix. In particular, their moduli are invariant under rotations and thus (non-linear) invariants of the 3-pt correlation function. In terms of these natural components, the invariance of the statistical properties of the shear field under parity transformations are easily obtained. Our results do not apply only to cosmic shear, but also to other quantities with the same mathematical properties -- that of a polar. For example, practically every relation derived here applies also to the polarization of the cosmic microwave background radiation.

Smoothing is omnipresent in astronomy, because almost always measurements performed at discrete positions in the sky need to be interpolated into a smooth map for subsequent analysis. Still, the statistical properties of different interpolation techniques are very poorly known. In this paper, we consider the general problem of interpolating discrete data whose location measurements are distributed on the sky according to a known density distribution (with or without clustering). We derive expressions for the expectation value and for the covariance of the smoothed map for many interpolation techniques, and obtain a general method that can be used to obtain these quantities for any linear smoothing. Moreover, we show that few basic properties of smoothing procedures have important consequences on the statistical properties of the smoothed map. Our analysis allows one to obtain the statistical properties of an arbitrary interpolation procedure, and thus to optimally choose the technique that is most suitable for one's needs.

Gravitationally lensed systems with multiply imaged quasars are an excellent tool for studying the properties of distant galaxies. In particular, they provide the most accurate mass measures for the lensing galaxy. The system B1422+231 is a well studied example of a quadruply imaged quasar, with high-quality VLBI data available. Very accurate data on image positions, fluxes and deconvolved image sizes provide good constraints for lensing models. We discuss here the failure of smooth models in fitting the data. Since it is intuitively clear that the mass of a lens galaxy is not a smooth entity, we have investigated how deviation from a smooth model can influence lensing phenomena, especially the image flux ratios. To explore expectations about the level of substructure in galaxies and its influence on strong lensing, N-body simulations of a model galaxy are employed. By using the mass distribution of this model galaxy as a lens, synthetic data sets of different four image system configurations are generated. Their analysis can possibly provide evidence for the presence and strong influence of substructure in the primary lens galaxy.

Interpolation techniques play a central role in Astronomy, where one often needs to smooth irregularly sampled data into a smooth map. In a previous article (Lombardi & Schneider 2001), we have considered a widely used smoothing technique and we have evaluated the expectation value of the smoothed map under a number of natural hypotheses. Here we proceed further on this analysis and consider the variance of the smoothed map, represented by a two-point correlation function. We show that two main sources of noise contribute to the total error budget and we show several interesting properties of these two noise terms. The expressions obtained are also specialized to the limiting cases of low and high densities of measurements. A number of examples are used to show in practice some of the results obtained.

In this work we investigate the gravitationally lensed system B1422+231. High--quality VLBI image positions, fluxes and shapes as well as an optical HST lens galaxy position are used. First, two simple and smooth models for the lens galaxy are applied to fit observed image positions and fluxes; no even remotely acceptable model was found. Such models also do not accurately reproduce the image shapes. In order to fit the data successfully, mass substructure has to be added to the lens, and its level is estimated. To explore expectations about the level of substructure in galaxies and its influence on strong lensing, N-body simulation results of a model galaxy are employed. By using the mass distribution of this model galaxy as a lens, synthetic data sets of different four image system configurations are generated and simple lens models are again applied to fit them. The difficulties in fitting these lens systems turn out to be similar to the case of some real gravitationally lensed systems, thus possibly providing evidence for the presence and strong influence of substructure in the primary lens galaxy.

The parameter-free reconstruction of the surface-mass density of clusters of galaxies is one of the principal applications of weak gravitational lensing. From the observable ellipticities of images of background galaxies, the tidal gravitational field (shear) of the mass distribution is estimated, and the corresponding surface mass density is constructed. The noise of the resulting mass map is investigated here, generalizing previous work which included mainly the noise due to the intrinsic galaxy ellipticities. Whereas this dominates the noise budget if the lens is very weak, other sources of noise become important, or even dominant, for the medium-strong lensing regime close to the center of clusters. In particular, shot noise due to a Poisson distribution of galaxy images, and increased shot noise owing to the correlation of galaxies in angular position and redshift, can yield significantly larger levels of noise than that from the intrinsic ellipticities only. We estimate the contributions from these various effects for two widely used smoothing operations, showing that one of them effectively removes the Poisson and the correlation noises related to angular positions of galaxies. Noise sources due to the spread in redshift of galaxies are still present in the optimized estimator and are shown to be relevant in many cases. We show how (even approximate) redshift information can be profitably used to reduce the noise in the mass map. The dependence of the various noise terms on the relevant parameters (lens redshift, strength, smoothing length, redshift distribution of background galaxies) are explicitly calculated and simple estimates are provided.

We generalize the technique of Lada et al. (1994) to map dust column density through a molecular cloud (NICE) to an optimized multi-band technique (NICER) that can be applied to any multi-band survey of molecular clouds. We present a first application to a ~625 deg^2 subset of the Two Micron All Sky Survey (2MASS) data and show that when compared to NICE, the optimized NICER technique (i) achieves the same extinction peak values, (ii) improves the noise variance of the map by a factor of 2 and (iii) is able to reach 3 \sigma dust extinction measurements as low as A_V = 0.5 magnitudes, better than or equivalent to classical optical star count techniques and below the threshold column density for the formation of CO, the brightest H2 tracer in radio-spectroscopy techniques. The application of the NICER techniques to near-infrared data obtained with a 8 meter-class telescope with a state-of-the-art NIR camera, such as the VLT-ISAAC combination, will be able to achieve dynamic ranges from below 10^21 protons cm^-2 to over 10^23 protons cm^-2 (A_V in the range [0.3, 60]) and spatial resolutions < 10'', making use of a single and straightforward column density tracer, extinction by interstellar dust.

We have re-examined the classical problem of the macroscopic equation of state for a hydrostatic isothermal self-gravitating gas cloud bounded by an external medium at constant pressure. We have obtained analytical conditions for its equilibrium and stability without imposing any specific shape and symmetry to the cloud density distribution. The equilibrium condition can be stated in the form of an upper limit to the cloud mass; this is found to be inversely proportional to the power 3/2 of a form factor \mu characterizing the shape of the cloud. In this respect, the spherical solution, associated with the maximum value of the form factor, \mu = 1, turns out to correspond to the shape that is most difficult to realize. Surprisingly, the condition that defines the onset of the Bonnor instability (or gravothermal catastrophe) can be cast in the form of an upper limit to the density contrast within the cloud that is independent of the cloud shape. We have then carried out a similar analysis in the two-dimensional case of infinite cylinders, without assuming axisymmetry. The results obtained in this paper generalize well-known results available for spherical or axisymmetric cylindrical isothermal clouds that have had wide astrophysical applications, especially in the study of the interstellar medium.

We study an estimator for smoothing irregularly sampled data into a smooth map. The estimator has been widely used in astronomy, owing to its low level of noise; it involves a weight function -- or smoothing kernel -- w(\theta). We show that this estimator is not unbiased, in the sense that the expectation value of the smoothed map is not the underlying process convolved with $w$, but a convolution with a modified kernel w_eff(\theta). We show how to calculate w_eff for a given kernel w and investigate its properties. In particular, it is found that (1) w_eff is normalized, (2) has a shape `similar' to the original kernel w, (3) converges to w in the limit of high number density of data points, and (4) reduces to a top-hat filter in the limit of very small number density of data points. Hence, although the estimator is biased, the bias is well understood analytically, and since w_eff has all the desired properties of a smoothing kernel, the estimator is in fact very useful. We present explicit examples for several filter functions which are commonly used, and provide a series expression valid in the limit of large density of data points.

We present an in-depth weak lensing analysis of the cluster MS1008 based on deep multicolor imaging obtained during the Science Verification of FORS1 at the VLT. The image quality (half arcsec seeing) and depth of the VLT images allow the shear signal to be mapped with high signal-to-noise and to be traced out to 1.2 h_50^-1 Mpc, near the edge of the 6'.8x6'.8 field of view. Using BVRI color information, as well as 81 redshifts in the field from the CNOC survey, background galaxies can be effectively separated from cluster and foreground objects. PSF distorsions are found to be moderate across the FORS images and thus easily removed. Due to the small statistical errors in the mass reconstruction, this dataset provides a testing ground where several systematic effects (e.g. mass-sheet degeneracy, redshift distribution of the background sources, cluster galaxy contamination) can be quantified. Several methods are used to remove the mass-sheet degeneracy which is found to dominate the systematic error budget. We measure a lower limit to the mass of 2.6 10^14 h_50^-1 M_sun within 1 h_50^-1 Mpc and a ``total'' mass of 5.3 10^14 h_50^-1 M_sun by fitting a softened isothermal sphere. The availability of the CNOC redshift data and X-ray observations on this cluster allow a comparison of different determinations of the mass radial profile. We find the lensing and X-ray measurements in excellent agreement, while the mass derived from the virial analysis is marginally (1-2 sigma) in agreement at radii where both methods are reliable.

The analysis of the shear induced by a single cluster on the images of a large number of background galaxies is all centered around the curl-free character of a well-known vector field that can be derived from the data. Such basic property breaks down when the source galaxies happen to be observed through two clusters at different redshifts, partially aligned along the line of sight. In this paper we address the study of double lenses and obtain five main results. (i) First we generalize the procedure to extract the available information, contained in the observed shear field, from the case of a single lens to that of a double lens. (ii) Then we evaluate the possibility of detecting the signature of double lensing given the known properties of the distribution of clusters of galaxies. (iii) As a different astrophysical application, we demonstrate how the method can be used to detect the presence of a dark cluster that might happen to be partially aligned with a bright cluster studied in terms of statistical lensing. (iv) In addition, we show that the redshift distribution of the source galaxies, which in principle might also contribute to break the curl-free character of the shear field, actually produces systematic effects typically two orders of magnitude smaller than the double lensing effects we are focusing on. (v) Remarkably, a discussion of relevant contributions to the noise of the shear measurement has brought up an intrinsic limitation of weak lensing analyses, since one specific contribution, associated with the presence of a non-vanishing two-galaxy correlation function, turns out not to decrease with the density of source galaxies (and thus with the depth of the observations).

Statistical analyses of observed galaxy distortions are often used to reconstruct the mass distribution of an intervening cluster responsible for gravitational lensing. In current projects, distortions of thousands of source galaxies have to be handled efficiently; much larger data bases and more massive investigations are envisaged for new major observational initiatives. In this article we present an efficient mass reconstruction procedure, a direct method that solves a variational principle noted in an earlier paper, which, for rectangular fields, turns out to reduce the relevant execution time by a factor from 100 to 1000 with respect to the fastest methods currently used, so that for grid numbers N = 400 the required CPU time on a good workstation can be kept within the order of 1 second. The acquired speed also opens the way to some long-term projects based on simulated observations (addressing statistical or cosmological questions) that would be, at present, practically not viable for intrinsically slow reconstruction methods.

Recently, it has been shown that it is possible to reconstruct the projected mass distribution of a cluster from weak lensing provided that both the geometry of the universe and the probability distribution of galaxy redshifts are known; actually, when additional photometric data are taken to be available, the galaxy redshift distribution could be determined jointly with the cluster mass from the weak lensing analysis. In this paper we develop, in the spirit of a ``thought experiment,'' a method to constrain the geometry of the universe from weak lensing, provided that the redshifts of the source galaxies are measured. The quantitative limits and merits of the method are discussed analytically and with a set of simulations, in relation to point estimation, interval estimation, and test of hypotheses for homogeneous Friedmann-Lemaitre models. The constraints turn out to be significant when a few thousand source galaxies are used.

In this paper we provide a statistical analysis of the parameter-free method often used in weak lensing mass reconstructions. It is found that a proper assessment of the errors involved in such a non-local analysis requires the study of the relevant two-point correlation functions. After calculating the two-point correlation function for the reduced shear, we determine the expected error on the inferred mass distribution and on other related quantities, such as the total mass, and derive the error power spectrum. This allows us to optimize the reconstruction method, with respect to the kernel used in the inversion procedure. In particular, we find that curl-free kernels are bound to lead to more accurate mass reconstructions. Our analytical results clarify the arguments and the numerical simulations by Seitz & Schneider (1996).

Different options can be used in order to measure the shear from observations in the context of weak lensing. Here we introduce new methods where the isotropy assumption for the distribution of the source galaxies is implemented directly on the observed quadrupole moments. A quantitative analysis of the error associated with the finite number of source galaxies and with their ellipticity distribution is provided, applicable even when the shear is not weak. Monte Carlo simulations based on a realistic sample of source galaxies show that our procedure generally leads to errors ~30% smaller than those associated with the standard method of Kaiser and Squires (1993).