Different studies have reported a power-law mass-size relation $M \propto R^q$ for ensembles of molecular clouds. In the case of nearby clouds, the index of the power-law $q$ is close to 2. However, for clouds spread all over the Galaxy, indexes larger than 2 are reported. We show that indexes larger than 2 could be the result of line-of-sight superposition of emission that does not belong to the cloud itself. We found that a random factor of gas contamination, between 0.001\%\ and 10\%\ of the line-of-sight, allows to reproduce the mass-size relation with $q \sim 2.2-2.3$ observed in Galactic CO surveys. Furthermore, for dense cores within a single cloud, or molecular clouds within a single galaxy, we argue that, even in these cases, there is observational and theoretical evidence that some degree of superposition may be occurring. However, additional effects may be present in each case, and are briefly discussed. We also argue that defining the fractal dimension of clouds via the mass-size relation is not adequate, since the mass is not {necessarily} a proxy to the area, and the size reported in $M-R$ relations is typically obtained from the square root of the area, rather than from an estimation of the size independent from the area. Finally, we argue that the statistical analysis of finding clouds satisfying the Larson's relations does not mean that each individual cloud is in virial equilibrium.

Aims:We take advantage of the second data release of the Gaia space mission and the state-of-the-art astrometry delivered from very long baseline interferometry observations to revisit the structure and kinematics of the nearby Taurus star-forming region. Methods: We apply a hierarchical clustering algorithm for partitioning the stars in our sample into groups (i.e., clusters) that are associated with the various molecular clouds of the complex, and derive the distance and spatial velocity of individual stars and their corresponding molecular clouds. Results: We show that the molecular clouds are located at different distances and confirm the existence of important depth effects in this region reported in previous studies. For example, we find that the L 1495 molecular cloud is located at $d=129.9^{+0.4}_{-0.3}$ pc, while the filamentary structure connected to it (in the plane of the sky) is at $d=160.0^{+1.2}_{-1.2}$ pc. We report B 215 and L 1558 as the closest ($d=128.5^{+1.6}_{-1.6}$ pc) and most remote ($d=198.1^{+2.5}_{-2.5}$ pc) substructures of the complex, respectively. The median inter-cloud distance is 25 pc and the relative motion of the subgroups is on the order of a few km/s. We find no clear evidence for expansion (or contraction) of the Taurus complex, but signs of the potential effects of a global rotation. Finally, we compare the radial velocity of the stars with the velocity of the underlying $^{13}$CO molecular gas and report a mean difference of $0.04\pm0.12$ km/s (with r.m.s. of 0.63 km/s) confirming that the stars and the gas are tightly coupled.

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

Recent analyses of Gaia data have provided direct evidence that most young stellar clusters are in a state of expansion, with velocities of the order of ~0.5 km/s. Traditionally, expanding young clusters have been pictured as entities that became unbound due to the lack of gravitational binding once the gas from the parental cloud that formed the cluster has been expelled by the stellar radiation of the massive stars in the cluster. In the present contribution, we used radiation-magnetohydrodynamic numerical simulations of molecular cloud formation and evolution to understand how stellar clusters form and disperse. We found that the ionising feedback from the newborn massive stars expels the gas from the collapse centre, flipping-up the gravitational potential as a consequence of the mass removal from the inside-out. Since neither the parental clouds nor the formed shells are distributed symmetrically around the HII region, net forces pulling out the stars are present, accelerating them towards the edges of the cavity. We call this mechanism ``gravitational feedback", in which the gravity from the expelled gas appears to be the crucial mechanism producing unbound clusters that expand away from their formation centre in an accelerated way in young stellar clusters. This mechanism naturally explains the "Hubble flow-like" expansion observed in several young clusters.

Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from APOGEE-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of $\sim$0.05--1.5 \msun, we find 399 binaries, including both RV variables and SB2s. The mass ratio distribution of SB2s is consistent with a uniform for $q<0.95$ with an excess of twins with $q>0.95$. The period distribution is consistent with what has been observed in close binaries ($<10$ AU) in the evolved populations. Three systems are found to have $q\sim$0.1, with a companion located within the brown dwarf desert. There are not any strong trends in the multiplicity fraction (MF) as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at $\Sigma_*\sim30$ stars/pc$^{-2}$, and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of $\sim$2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal mass systems that is less effective in binaries with lower $q$.

The stellar initial mass function (IMF) is an essential input for many astrophysical studies but only in a few cases it has been determined over the whole cluster mass range, limiting the conclusions about its nature. The 25 Orionis group (25 Ori) is an excellent laboratory to investigate the IMF across the entire mass range of the population, from planetary-mass objects to intermediate/high-mass stars. We combine new deep optical photometry with optical and near-infrared data from the literature to select 1687 member candidates covering a 1.1$^\circ$ radius area in 25 Ori. With this sample we derived the 25 Ori system IMF from 0.012 to 13.1 $M_\odot$. This system IMF is well described by a two-segment power-law with $\Gamma=-0.74\pm0.04$ for $m<0.4\ M_\odot$ and $\Gamma=1.50\pm0.11$ for $m\ge0.4\ M_\odot$. It is also well described over the whole mass range by a tapered power-law function with $\Gamma=1.10\pm0.09$, $m_p=0.31\pm0.03$ and $\beta=2.11\pm0.09$. The best lognormal representation of the system IMF has $m_c=0.31\pm0.04$ and $\sigma=0.46\pm0.05$ for $m<1\ M_\odot$. This system IMF does not present significant variations with the radii. We compared the resultant system IMF as well as the BD/star ratio of $0.16\pm0.03$ we estimated for 25 Ori with that of other stellar regions with diverse conditions and found no significant discrepancies. These results support the idea that general star formation mechanisms are probably not strongly dependent to environmental conditions. We found that the substellar and stellar objects in 25 Ori have similar spatial distributions and confirmed that 25 Ori is a gravitationally unbound stellar association.

The nearest massive star-forming complex, Cygnus-X, is widely used as a laboratory for star cluster formation and feedback processes, under the implicit assumption that all its components are located roughly at the same distance. We present a multi-wavelength study of a 15' x 15' field in southern Cygnus-X, where different components involving clustered star formation are overlapped. Preliminary results indicate that the Berkeley 87 and [DB2001] CL05 clusters are actually located at very different distances, invalidating previous claims of physical interaction between them. This shows the importance of a careful treatment of extinction and distance calculations for cluster formation studies, particularly in Cygnus-X.

We investigate the prompt emission and the afterglow properties of short duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multi-wavelength facilities including the GTC 10.4m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL/SPI/ACS, Swift/BAT and Fermi/GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB130603B, our new afterglow photometric data constraints the pre jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new as well as previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for 8 other sGRBs we are able to; place tight constraints on the non-detection of the afterglow, host galaxy or any underlying kilonova emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs.

We present optical photometry of the afterglow of the long GRB 180205A with the COATLI telescope from 217 seconds to about 5 days after the {\itshape Swift}/BAT trigger. We analyse this photometry in the conjunction with the X-ray light curve from {\itshape Swift}/XRT. The late-time light curves and spectra are consistent with the standard forward-shock scenario. However, the early-time optical and X-ray light curves show non-typical behavior; the optical light curve exhibits a flat plateau while the X-ray light curve shows a flare. We explore several scenarios and conclude that the most likely explanation for the early behavior is late activity of the central engine.

In this letter we present evidence of coherent outward motion of a sample of young stars ($t<$30 Myr) in the Perseus Arm, whose apparent origin is located in the vicinity of the W3/W4/W5 complex. Using astrometric and photometric data from the Gaia DR2 catalog of an 8$^\circ$ radius field centered near W4, we selected a sample of young, intermediate to high-mass star candidates. The sample is limited to sources with parallax uncertainties below 20% and Bayesian distance estimates within 1800 and 3100 pc. The selection includes embedded stellar populations as well as young open clusters. Projected velocities derived from perspective-corrected proper motions clearly suggest that the young star population emerge from the Perseus arm, with a possible convergence zone near W3/W4/W5 region, tracing a front that expands away at a rate of about $15~{\rm km~s}^{-1}~{\rm kpc}^{-1}$.

We report new, $\sim$1000 AU spatial resolution observations of 225 GHz dust continuum emission towards the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive ($\gtrsim$100 $M_{\odot}$) cores. From the new, higher resolution observations, we identified 28 localized, spatially compact dust continuum emission sources, which may be candidates of young stellar objects. Only one of them is not embedded within known arm-like (or elongated) dense gas structures. The spatial separations of these compact sources can be very well explained by Jeans lengths. We found that G33.92+0.11 may be consistently described by a marginally centrifugally supported, Toomre unstable accretion flow which is approximately in a face-on projection. The arm-like overdensities are natural consequence of the Toomre instability, which can fragment to form young stellar objects in shorter time scales than the timescale of the global clump contraction. On our resolved spatial scales, there is not yet evidence that the fragmentation is halted by turbulence, magnetic field, or stellar feedback.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.

We present an analysis of spectrosopic and astrometric data from APOGEE-2 and Gaia DR2 to identify structures towards the Orion Complex. By applying a hierarchical clustering algorithm to the 6-dimensional stellar data, we identify spatially and/or kinematically distinct groups of young stellar objects with ages ranging from 1 to 12 Myr. We also investigate the star forming history within the Orion Complex, and identify peculiar sub-clusters. With this method we reconstruct the older populations in the regions that are presently largely devoid of molecular gas, such as Orion C (which includes the $\sigma$ Ori cluster), and Orion D (the population that traces Ori OB1a, OB1b, and Orion X). We report on the distances, kinematics, and ages of the groups within the Complex. The Orion D groups is in the process of expanding. On the other hand, Orion B is still in the process of contraction. In $\lambda$ Ori the proper motions are consistent with a radial expansion due to an explosion from a supernova; the traceback age from the expansion exceeds the age of the youngest stars formed near the outer edges of the region, and their formation would have been triggered when they were half-way from the cluster center to their current positions. We also present a comparison between the parallax and proper motion solutions obtained by Gaia DR2, and those obtained towards star-forming regions by Very Long Baseline Array.

We analyse the geometry of the Pipe Nebula, drawn by the distribution ($Q$-spatial parameter) and hierarchy ($\Lambda$ spatial segregation) of column density peaks previously detected and catalogued. By analysing the mass and volume density of the cores, we determine that both variables shown to be spatially segregated and with a high degree of substructure. Given the early evolutionary state of the Pipe Nebula, our results suggest that both, mass and volume-density segregations, may be primordial, in the sense of appearing in the early phases of the chain of physical mechanisms which conform the cluster-formation process. We also propose that volume density, and not mass, is the parameter that most clearly determines the initial spatial distribution of pre-stellar cores.

The Orion Star Forming Complex (OSFC) is a central target for the APOGEE-2 Young Cluster Survey. Existing membership catalogs span limited portions of the OSFC, reflecting the difficulty of selecting targets homogeneously across this extended, highly structured region. We have used data from wide field photometric surveys to produce a less biased parent sample of young stellar objects (YSOs) with infrared (IR) excesses indicative of warm circumstellar material or photometric variability at optical wavelengths across the full 420 square degrees extent of the OSFC. When restricted to YSO candidates with H < 12.4, to ensure S/N ~100 for a six visit source, this uniformly selected sample includes 1307 IR excess sources selected using criteria vetted by Koenig & Liesawitz and 990 optical variables identified in the Pan-STARRS1 3$\pi$ survey: 319 sources exhibit both optical variability and evidence of circumstellar disks through IR excess. Objects from this uniformly selected sample received the highest priority for targeting, but required fewer than half of the fibers on each APOGEE-2 plate. We fill the remaining fibers with previously confirmed and new color-magnitude selected candidate OSFC members. Radial velocity measurements from APOGEE-1 and new APOGEE-2 observations taken in the survey's first year indicate that ~90% of the uniformly selected targets have radial velocities consistent with Orion membership.The APOGEE-2 Orion survey will include >1100 bona fide YSOs whose uniform selection function will provide a robust sample for comparative analyses of the stellar populations and properties across all sub-regions of Orion.

Recently we have witnessed the first multi-messenger detection of colliding neutron stars through Gravitational Waves (GWs) and Electromagnetic (EM) waves (GW170817), thanks to the joint efforts of LIGO/Virgo and Space/Ground-based telescopes. In this paper, we report on the RATIR followup observation strategies and show the results for the trigger G194575. This trigger is not of astrophysical interest; however, is of great interests to the robust design of a followup engine to explore large sky error regions. We discuss the development of an image-subtraction pipeline for the 6-color, optical/NIR imaging camera RATIR. Considering a two band ($i$ and $r$) campaign in the Fall of 2015, we find that the requirement of simultaneous detection in both bands leads to a factor $\sim$10 reduction in false alarm rate, which can be further reduced using additional bands. We also show that the performance of our proposed algorithm is robust to fluctuating observing conditions, maintaining a low false alarm rate with a modest decrease in system efficiency that can be overcome utilizing repeat visits. Expanding our pipeline to search for either optical or NIR detections (3 or more bands), considering separately the optical $riZ$ and NIR $YJH$ bands, should result in a false alarm rate $\approx 1\%$ and an efficiency $\approx 90\%$. RATIR's simultaneous optical/NIR observations are expected to yield about one candidate transient in the vast 100 $\mathrm{deg^2}$ LIGO error region for prioritized followup with larger aperture telescopes.

In this work we make use of DR14 APOGEE spectroscopic data to study a sample of 92 known OB stars. We developed a near-infrared semi-empirical spectral classification method that was successfully used in case of four new exemplars, previously classified as later B-type stars. Our results agree well with those determined independently from ECHELLE optical spectra, being in line with the spectral types derived from the "canonical" MK blue optical system. This confirms that the APOGEE spectrograph can also be used as a powerful tool in surveys aiming to unveil and study large number of moderately and highly obscured OB stars still hidden in the Galaxy.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

The Orion OB1a sub-association is a rich low mass star (LMS) region. Previous spectroscopic studies have confirmed 160 LMSs in the 25 Orionis stellar group (25 Ori), which is the most prominent overdensity of Orion OB1a. Nonetheless, the current census of the 25 Ori members is estimated to be less than 50% complete, leaving a large number of members to be still confirmed. We retrieved 172 low-resolution stellar spectra in Orion OB1a observed as ancillary science in the SDSS-III/BOSS survey, for which we classified their spectral types and determined physical parameters. To determine memberships, we analyzed the H$_\alpha$ emission, LiI$\lambda$6708 absorption, and NaI$\lambda\lambda$8183, 8195 absorption as youth indicators in stars classified as M-type. We report 50 new LMSs spread across the 25 Orionis, ASCC 18, and ASCC 20 stellar groups with spectral types from M0 to M6, corresponding to a mass range of 0.10$\le m/\textrm{M}_\odot \le$0.58. This represents an increase of 50% in the number of known LMSs in the area and a net increase of 20% in the number of 25 Ori members in this mass range. Using parallax values from the Gaia DR1 catalog, we estimated the distances to these three stellar groups and found that they are all co-distant, at 338$\pm$66 pc. We analyzed the spectral energy distributions of these LMSs and classified their disks by evolutionary classes. Using H-R diagrams, we found a suggestion that 25 Ori could be slightly older that the other two observed groups in Orion OB1a.

We observed the afterglow of GRB 130427A with the RATIR instrument on the 1.5-m Harold L. Johnson telescope of the Observatorio Astron\'omico Nacional on Sierra San Pedro M\'artir. Our homogenous $griZYJH$ photometry extends from the night of burst to three years later. We fit a model for the afterglow. There is a significant positive residual which matches the behavior of SN 1998bw in the $griZ$ filters; we suggest that this is a photometric signature of the supernova SN 2013cq associated with the GRB. The peak absolute magnitude of the supernova is $M_r=-18.43\pm0.11$.

We report on the first open-use based Atacama Large Millimeter/submm Array (ALMA) 345-GHz observation for the late afterglow phase of GRB131030A. The ALMA observation constrained a deep limit at 17.1 d for the afterglow and host galaxy. We also identified a faint submillimeter source (ALMAJ2300-0522) near the GRB131030A position. The deep limit at 345 GHz and multifrequency observations obtained using {\it Swift} and RATIR yielded forward shock modeling with a two-dimensional relativistic hydrodynamic jet simulation and described X-ray excess in the afterglow. The excess was inconsistent with the synchrotron self-inverse Compton radiation from the forward shock. The host galaxy of GRB131030A and optical counterpart of ALMAJ2300-0522 were also identified in the SUBARU image. Based on the deep ALMA limit for the host galaxy, the 3-$\sigma$ upper limits of IR luminosity and the star formation rate (SFR) is estimated as $L_{IR}<1.11\times10^{11} L_{\odot}$ and SFR$<18.7$ ($M_{\odot}$~yr$^{-1}$), respectively. Although the separation angle from the burst location (3.5 arcsec) was rather large, ALMAJ2300-0522 may be one component of the GRB131030A host galaxy, according to previous host galaxy cases.

COATLI will provide 0.3 arcsec FWHM images from 550 to 900 nm over a large fraction of the sky. It consists of a robotic 50-cm telescope with a diffraction-limited fast-guiding imager. Since the telescope is small, fast guiding will provide diffraction-limited image quality over a field of at least 1 arcmin and with coverage of a large fraction of the sky, even in relatively poor seeing. The COATLI telescope will be installed at the at the Observatorio Astron\'omico Nacional in Sierra San Pedro M\'artir, M\'exico, during 2016 and the diffraction-limited imager will follow in 2017.

DDOTI will be a wide-field robotic imager consisting of six 28-cm telescopes with prime focus CCDs mounted on a common equatorial mount. Each telescope will have a field of view of 12 square degrees, will have 2 arcsec pixels, and will reach a 10-sigma limiting magnitude in 60 seconds of r = 18.7 in dark time and r = 18.0 in bright time. The set of six will provide an instantaneous field of view of about 72 square degrees. DDOTI uses commercial components almost entirely. The first DDOTI will be installed at the Observatorio Astron\'omico Nacional in Sierra San Pedro Mart\'ir, Baja California, M\'exico in early 2017. The main science goals of DDOTI are the localization of the optical transients associated with GRBs detected by the GBM instrument on the Fermi satellite and with gravitational-wave transients. DDOTI will also be used for studies of AGN and YSO variability and to determine the occurrence of hot Jupiters. The principal advantage of DDOTI compared to other similar projects is cost: a single DDOTI installation costs only about US$500,000. This makes it possible to contemplate a global network of DDOTI installations. Such geographic diversity would give earlier access and a higher localization rate. We are actively exploring this option.

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 present optical and near-infrared light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be $E_{\gamma,iso}=6.4_{-1.0}^{+1.3}\times10^{50}$erg (1keV-10MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed $g'r'i'z'$ light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves ~20% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined Type Ic supernovae, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of ~21,000km s$^{-1}$. We construct a quasi-bolometric ($g'r'i'z'yJ$) light curve for SN 2013dx, only the fifth GRB-associated SN with extensive NIR coverage and the third with a bolometric light curve extending beyond $\Delta t>40$d. Together with the measured photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a $^{56}$Ni mass of $M_{\mathrm{Ni}}=0.37\pm0.01$M$_{\odot}$, an ejecta mass of $M_{\mathrm{ej}}=3.1\pm0.1$M$_{\odot}$, and a kinetic energy of $E_{\mathrm{K}}=(8.2\pm0.43)\times10^{51}$ erg (statistical uncertainties only), consistent with previous GRB-associated SNe. When considering the ensemble population of GRB-associated SNe, we find no correlation between the mass of synthesized $^{56}$Ni and high-energy properties, despite clear predictions from numerical simulations that $M_{\mathrm{Ni}}$ should correlate with the degree of asymmetry. On the other hand, $M_{\mathrm{Ni}}$ clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.

We present a multi-wavelength study of the young stellar population in the Cygnus-X DR15 region. We studied young stars forming or recently formed at and around the tip of a prominent molecular pillar and an infrared dark cloud. Using a combination of ground based near-infrared, space based infrared and X-ray data, we constructed a point source catalog from which we identified 226 young stellar sources, which we classified into evolutionary classes. We studied their spatial distribution across the molecular gas structures and identified several groups possibly belonging to distinct young star clusters. We obtained samples of these groups and constructed K-band luminosity functions that we compared with those of artificial clusters, allowing us to make first order estimates of the mean ages and age spreads of the groups. We used a $^{13}$CO(1-0) map to investigate the gas kinematics at the prominent gaseous envelope of the central cluster in DR15, and we infer that the removal of this envelope is relatively slow compared to other cluster regions, in which gas dispersal timescale could be similar or shorter than the circumstellar disk dissipation timescale. The presence of other groups with slightly older ages, associated with much less prominent gaseous structures may imply that the evolution of young clusters in this part of the complex proceeds in periods that last 3 to 5 Myr, perhaps after a slow dissipation of their dense molecular cloud birthplaces.

Gamma-ray bursts (GRBs) are the most luminous explosions in the universe, yet the nature and physical properties of their energy sources are far from understood. Very important clues, however, can be inferred by studying the afterglows of these events. We present optical and X-ray observations of GRB 130831A obtained by Swift, Chandra, Skynet, RATIR, Maidanak, ISON, NOT, LT and GTC. This burst shows a steep drop in the X-ray light-curve at $\simeq 10^5$ s after the trigger, with a power-law decay index of $\alpha \sim 6$. Such a rare behaviour cannot be explained by the standard forward shock (FS) model and indicates that the emission, up to the fast decay at $10^5$ s, must be of "internal origin", produced by a dissipation process within an ultrarelativistic outflow. We propose that the source of such an outflow, which must produce the X-ray flux for $\simeq 1$ day in the cosmological rest frame, is a newly born magnetar or black hole. After the drop, the faint X-ray afterglow continues with a much shallower decay. The optical emission, on the other hand, shows no break across the X-ray steep decrease, and the late-time decays of both the X-ray and optical are consistent. Using both the X-ray and optical data, we show that the emission after $\simeq 10^5$ s can be explained well by the FS model. We model our data to derive the kinetic energy of the ejecta and thus measure the efficiency of the central engine of a GRB with emission of internal origin visible for a long time. Furthermore, we break down the energy budget of this GRB into the prompt emission, the late internal dissipation, the kinetic energy of the relativistic ejecta, and compare it with the energy of the associated supernova, SN 2013fu.

We present a multi-wavelength analysis of the history of star formation in the W3 complex. Using deep, near-infrared ground-based images, combined with images obtained with Spitzer and Chandra observatories, we identified and classified young embedded sources. We identified the principal clusters in the complex, and determined their structure and extension. We constructed extinction-limited samples for five principal clusters, and constructed K-band luminosity functions (KLF) that we compare with those of artificial clusters with varying ages. This analysis provided mean ages and possible age spreads for the clusters. We found that IC 1795, the centermost cluster of the complex, still hosts a large fraction of young sources with circumstellar disks. This indicates that star formation was active in IC 1795 as recently as 2 Myr ago, simultaneous to the star forming activity in the flanking embedded clusters, W3-Main and W3(OH). A comparison with carbon monoxide emission maps indicates strong velocity gradients in the gas clumps hosting W3-Main and W3(OH) and show small receding clumps of gas at IC 1795, suggestive of rapid gas removal (faster than the T Tauri timescale) in the cluster forming regions. We discuss one possible scenario for the progression of cluster formation in the W3 complex. We propose that early processes of gas collapse in the main structure of the complex could have defined the progression of cluster formation across the complex with relatively small age differences from one group to another. However, triggering effects could act as catalysts for enhanced efficiency of formation at a local level, in agreement with previous studies.

How rapidly collapsing parsec-scale massive molecular clumps feed high-mass stars, and how they fragment to form OB clusters, have been outstanding questions in the field of star-formation. In this work, we report the resolved structures and kinematics of the approximately face-on, rotating massive molecular clump, G33.92+0.11. Our high resolution Atacama Large Millimeter/submillimeter Array (ALMA) images show that the spiral arm-like gas overdensities form in the eccentric gas accretion streams. First, we resolved that the dominant part of the $\sim$0.6 pc scale massive molecular clump (3.0$^{+2.8}_{-1.4}$$\cdot$10$^{3}$ $M_{\odot}$) G33.92+0.11 A is tangled with several 0.5-1 pc size molecular arms spiraling around it, which may be connected further to exterior gas accretion streams. Within G33.92+0.11 A, we resolved the $\sim$0.1 pc width gas mini-arms connecting with the two central massive (100-300 $M_{\odot}$) molecular cores. The kinematics of arms and cores elucidate a coherent accretion flow continuing from large to small scales. We demonstrate that the large molecular arms are indeed the cradles of dense cores, which are likely current or future sites of high-mass star formation. Since these deeply embedded massive molecular clumps preferentially form the highest mass stars in the clusters, we argue that dense cores fed by or formed within molecular arms play a key role in making the upper end of the stellar and core mass functions.

We present a study of 15 new brown dwarfs belonging to the $\sim7$ Myr old 25 Orionis group and Orion OB1a sub-association with spectral types between M6 and M9 and estimated masses between $\sim0.07$M$_\odot$ and $\sim0.01$ M$_\odot$. By comparing them through a Bayesian method with low mass stars ($0.8\lesssim$ M/M$_\odot\lesssim0.1$) from previous works in the 25 Orionis group, we found statistically significant differences in the number fraction of classical T Tauri stars, weak T Tauri stars, class II, evolved discs and purely photospheric emitters at both sides of the sub-stellar mass limit. Particularly we found a fraction of $3.9^{+2.4}_{-1.6}~\%$ low mass stars classified as CTTS and class II or evolved discs, against a fraction of $33.3^{+10.8}_{-9.8}~\%$ in the sub-stellar mass domain. Our results support the suggested scenario in which the dissipation of discs is less efficient for decreasing mass of the central object.

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.

We present optical and near-infrared (NIR) photometry of 28 gamma-ray bursts (GRBs) detected by the \textit{Swift} satellite and rapidly observed by the Reionization and Transients Infrared/Optical (RATIR) camera. We compare the optical flux at fiducial times of 5.5 and 11 hours after the high-energy trigger to that in the X-ray regime to quantify optical darkness. 46$\pm$9 per cent (13/28) of all bursts in our sample and 55$\pm$10 per cent (13/26) of long GRBs are optically dark, which is statistically consistently with previous studies. Fitting RATIR optical and NIR spectral energy distributions (SEDs) of 19 GRBs, most (6/7) optically dark GRBs either occur at high-redshift ($z>4.5$) or have a high dust content in their host galaxies ($A_{\rm V} > 0.3$). Performing K-S tests, we compare the RATIR sample to those previously presented in the literature, finding our distributions of redshift, optical darkness, host dust extinction and X-ray derived column density to be consistent. The one reported discrepancy is with host galaxy dust content in the BAT6 sample, which appears inconsistent with our sample and other previous literature. Comparing X-ray derived host galaxy hydrogen column densities to host galaxy dust extinction, we find that GRBs tend to occur in host galaxies with a higher metal-to-dust ratio than our own Galaxy, more akin to the Large and Small Magellanic Clouds. Finally, to mitigate time evolution of optical darkness, we measure $\beta_{\rm OX,rest}$ at a fixed rest frame time, $t_{\rm rest}=1.5$ hours and fixed rest frame energies in the X-ray and optical regimes. Choosing to evaluate optical flux at $\lambda_{\rm rest}=0.25~\mu$m, we remove high-redshift as a source of optical darkness, demonstrating that optical darkness must result from either high-redshift, dust content in the host galaxy along the GRB sight line, or a combination of the two.

Context. The Pipe nebula is a molecular cloud that lacks star formation feedback and has a relatively simple morphology and velocity structure. This makes it an ideal target to test cloud evolution through collisions. Aims. We aim at drawing a comprehensive picture of this relatively simple cloud to better understand the formation and evolution of molecular clouds on large scales. Methods. We use archival data to compare the optical polarization properties, the visual extinction, and the 13CO velocities and linewidths of the entire cloud in order to identify trends among the observables. Results. The Pipe nebula can be roughly divided in two filaments with different orientations and gas velocity ranges: E-W at 2-4 km s-1 and N-S at 6-7 km s-1. The two filaments overlap at the bowl, where the gas shows a velocity gradient spanning from 2 to 7 km s-1. Compared to the rest of the Pipe nebula, the bowl gas appears to be denser and exhibits larger linewidths. In addition, the polarization data at the bowl shows lower angular dispersion and higher polarization degree. Cores in the bowl tend to cluster in space and tend to follow the 13CO velocity gradient. In the stem, cores tend to cluster in regions with properties similar to those of the bowl. Conclusions. The velocity pattern points to a collision between the filaments in the bowl region. The magnetic field seems to be compressed and strengthened in the shocked region. The proportional increase of density and magnetic field strength by a factor similar to the Alfv\'enic Mach number suggests a continuous shock at low Alfv\'enic Mach number under flux-freezing. Shocked regions seem to enhance the formation and clustering of dense cores.

The infrared data from the Spitzer Space Telescope has provided an invaluable tool for identifying physical processes in star formation. In this study we calculate the IRAC color space of UV fluorescent molecular hydrogen (H$_2$) and Polycyclic Aromatic Hydrocarbon (PAH) emission in photodissociation regions (PDRs) using the Cloudy code with PAH opacities from Draine & Li 2007. We create a set of color diagnostics that can be applied to study the structure of PDRs and to distinguish between FUV excited and shock excited H$_2$ emission. To test this method we apply these diagnostics to Spitzer IRAC data of NGC 2316. Our analysis of the structure of the PDR is consistent with previous studies of the region. In addition to UV excited emission, we identify shocked gas that may be part of an outflow originating from the cluster.

We present a template fitting algorithm for determining photometric redshifts, $z_{\rm phot}$, of candidate high-redshift gamma-ray bursts (GRBs). Using afterglow photometry, obtained by the Reionization And Transients InfraRed (RATIR) camera, this algorithm accounts for the intrinsic GRB afterglow spectral energy distribution (SED), host dust extinction and the effect of neutral hydrogen (local and cosmological) along the line of sight. We present the results obtained by this algorithm and RATIR photometry of GRB 130606A, finding a range of best fit solutions $5.6 < z_{\rm phot} < 6.0$ for models of several host dust extinction laws (none, MW, LMC and SMC), consistent with spectroscopic measurements of the redshift of this GRB. Using simulated RATIR photometry, we find our algorithm provides precise measures of $z_{\rm phot}$ in the ranges $4 < z_{\rm phot} \lesssim 8$ and $9 < z_{\rm phot} < 10$ and can robustly determine when $z_{\rm phot}>4$. Further testing highlights the required caution in cases of highly dust extincted host galaxies. These tests also show that our algorithm does not erroneously find $z_{\rm phot} < 4$ when $z_{\rm sim}>4$, thereby minimizing false negatives and allowing us to rapidly identify all potential high-redshift events.

Sh 2-29 is a conspicuous star-forming region marked by the presence of massive embedded stars as well as several notable interstellar structures. In this research, our goals were to determine the role of magnetic fields and to study the size distribution of interstellar dust particles within this turbulent environment. We have used a set of optical and near-infrared polarimetric data obtained at OPD/LNA (Brazil) and CTIO (Chile), correlated with extinction maps, 2MASS data and images from DSS and Spitzer. The region's most striking feature is a swept out interstellar cavity whose polarimetric maps indicate that magnetic field lines were dragged outwards, pilling up along its borders. This led to a higher magnetic strength value ($\approx400\,\mu$G) and an abrupt increase in polarization degree, probably due to an enhancement in alignment efficiency. Furthermore, dense cloud fragmentations with peak $A_{V}$ between 20 and 37 mag were probably triggered by its expansion. The presence of $24\,\mu$m point-like sources indicates possible newborn stars inside this dense environment. A statistical analysis of the angular dispersion function revealed areas where field lines are aligned in a well-ordered pattern, seemingly due to compression effects from the HII region expansion. Finally, Serkowski function fits were used to study the ratio of the total-to-selective extinction, reveling a dual population of anomalous grain particles sizes. This trend suggests that both effects of coagulation and fragmentation of interstellar grains are present in the region.

Using Spitzer Space Telescope and Chandra X-ray Observatory data, we identify YSOs in the Rosette Molecular Cloud (RMC). By being able to select cluster members and classify them into YSO types, we are able to track the progression of star formation locally within the cluster environments and globally within the cloud. We employ nearest neighbor method (NNM) analysis to explore the density structure of the clusters and YSO ratio mapping to study age progressions in the cloud. We find a relationship between the YSO ratios and extinction which suggests star formation occurs preferentially in the densest parts of the cloud and that the column density of gas rapidly decreases as the region evolves. This suggests rapid removal of gas may account for the low star formation efficiencies observed in molecular clouds. We find that the overall age spread across the RMC is small. Our analysis suggests that star formation started throughout the complex around the same time. Age gradients in the cloud appear to be localized and any effect the HII region has on the star formation history is secondary to that of the primordial collapse of the cloud.

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 dense molecular clump at the center of the Barnard 59 (B59) complex is the only region in the Pipe Nebula that has formed a small,stellar cluster. The previous analysis of a high resolution near-IR dust extinction map revealed that the nuclear region in B59 is a massive, mostly quiescent clump of 18.9 solar masses. The clump shows a monolithic profile, possibly indicating that the clump is on the way to collapse, with no evident fragmentation that could lead to another group of star systems. In this paper we present new analysis that compares the dust extinction map with a new dust emission radio-continuum map of higher spatial resolution. We confirm that the clump does not show any significant evidence for prestellar fragmentation at scales smaller than those probed previously.

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.

We present a study of dense structures in the L1495 filament in the Taurus Molecular Cloud and examine its star-forming properties. In particular we construct a dust extinction map of the filament using deep near-infrared observations, exposing its small-scale structure in unprecedented detail. The filament shows highly fragmented substructures and a high mass-per-length value of M_line=17 Msun/pc, reflecting star-forming potential in all parts of it. However, a part of the filament, namely B211, is remarkably devoid of young stellar objects. We argue that in this region the initial filament collapse and fragmentation is still taking place and star formation is yet to occur. In the star-forming part of the filament, we identify 39 cores with masses from 0.4...10 Msun and preferred separations in agreement with the local Jeans length. Most of these cores exceed the Bonnor-Ebert critical mass, and are therefore likely to collapse and form stars. The Dense Core Mass Function follows a power law with exponent Gamma=1.2, a form commonly observed in star-forming regions.

We present Chandra X-ray images of the NGC 2237 young star cluster on the periphery of the Rosette Nebula. We detect 168 X-ray sources, 80% of which have counterparts in USNO, 2MASS, and deep FLAMINGOS images. These constitute the first census of the cluster members with 0.2<~M<~2 Msun. Star locations in near-infrared color-magnitude diagrams indicate a cluster age around 2 Myr with a visual extinction of 1<Av<3 at 1.4 kpc, the distance of the Rosette Nebula's main cluster NGC 2244. We derive the K-band luminosity function and the X-ray luminosity function of the cluster, which indicate a population ~400-600 stars. The X-ray-selected sample shows a K-excess disk frequency of 13%. The young Class II counterparts are aligned in an arc ~3 pc long suggestive of a triggered formation process induced by the O stars in NGC 2244. The diskless Class III sources are more dispersed. Several X-ray emitting stars are located inside the molecular cloud and around gaseous pillars projecting from the cloud. These stars, together with a previously unreported optical outflow originating inside the cloud, indicate that star formation is continuing at a low level and the cluster is still growing. This X-ray view of young stars on the western side of the Rosette Nebula complements our earlier studies of the central cluster NGC 2244 and the embedded clusters on the eastern side of the Nebula. The large scale distribution of the clusters and molecular material is consistent with a scenario in which the rich central NGC 2244 cluster formed first, and its expanding HII region triggered the formation of the now-unobscured clusters RMC XA and NGC 2237. A large swept-up shell material around the HII region is now in a second phase of collect-and-collapse fragmentation, leading to the recent formation of subclusters. Other clusters deeper in the molecular cloud appear unaffected by the Nebula expansion.

We present our analysis of a fully sampled, high resolution dust extinction map of the Barnard 59 complex in the Pipe Nebula. The map was constructed with the infrared color excess technique applied to a photometric catalog that combines data from both ground and space based observations. The map resolves for the first time the high density center of the main core in the complex, that is associated with the formation of a small cluster of stars. We found that the central core in Barnard 59 shows an unexpected lack of significant substructure consisting of only two significant fragments. Overall, the material appears to be consistent with being a single, large core with a density profile that can be well fit by a King model. A series of NH$_3$ pointed observations towards the high column density center of the core appear to show that the core is still thermally dominated, with sub-sonic non-thermal motions. The stars in the cluster could be providing feedback to support the core against collapse, but the relatively narrow radio lines suggest that an additional source of support, for example a magnetic field, may be required to stabilize the core. Outside the central core our observations reveal the structure of peripheral cores and resolve an extended filament into a handful of significant substructures whose spacing and masses appear to be consistent with Jeans fragmentation.

The Rosette Complex in the constellation of Monoceros is a magnificent laboratory for the study of star formation. The region presents an interesting scenario, in which an expanding HII region generated by the large OB association NGC 2244 is interacting with a giant molecular cloud. Inside the cloud a number of stellar clusters have formed recently. In this chapter we present a review of past and present research on the region, and discuss investigations relevant to the physics of the nebula and the molecular cloud. We also review recent work on the younger embedded clusters and individual nebulous objects located across this important star forming region.

In studies of star-forming regions, near-infrared excess (NIRX) sources--objects with intrinsic colors redder than normal stars--constitute both signal (young stars) and noise (e.g. background galaxies). We hunt down (identify) galaxies using near-infrared observations in the Perseus star-forming region by combining structural information, colors, and number density estimates. Galaxies at moderate redshifts (z = 0.1 - 0.5) have colors similar to young stellar objects (YSOs) at both near- and mid-infrared (e.g. Spitzer) wavelengths, which limits our ability to identify YSOs from colors alone. Structural information from high-quality near-infrared observations allows us to better separate YSOs from galaxies, rejecting 2/5 of the YSO candidates identified from Spitzer observations of our regions and potentially extending the YSO luminosity function below K of 15 magnitudes where galaxy contamination dominates. Once they are identified we use galaxies as valuable extra signal for making extinction maps of molecular clouds. Our new iterative procedure: the Galaxies Near Infrared Color Excess method Revisited (GNICER), uses the mean colors of galaxies as a function of magnitude to include them in extinction maps in an unbiased way. GNICER increases the number of background sources used to probe the structure of a cloud, decreasing the noise and increasing the resolution of extinction maps made far from the galactic plane.

We combined sensitive near-infrared data obtained with ground-based imagers on the ESO NTT and VLT telescopes with space mid-infrared data acquired with the IRAC imager on the Spitzer Space Telescope to calculate the extinction law A_\lambda/A_K as a function of \lambda between 1.25 and 7.76 micron to an unprecedented depth in Barnard 59, a star forming, dense core located in the Pipe Nebula. The ratios A_\lambda/A_K were calculated from the slopes of the distributions of sources in color-color diagrams \lambda-K vs. H-K. The distributions in the color-color diagrams are fit well with single slopes to extinction levels of A_K ~ 7 (A_V ~ 59 mag). Consequently, there appears to be no significant variation of the extinction law with depth through the B59 line of sight. However, when slopes are translated into the relative extinction coefficients A_\lambda/A_K, we find an extinction law which departs from the simple extrapolation of the near-infrared power law extinction curve, and agrees more closely with a dust extinction model for a cloud with a total to selective absorption R_V=5.5 and a grain size distribution favoring larger grains than those in the diffuse ISM. Thus, the difference we observe could be possibly due to the effect of grain growth in denser regions. Finally, the slopes in our diagrams are somewhat less steep than those from the study of Indebetouw et al. (2005) for clouds with lower column densities, and this indicates that the extinction law between 3 and 8 micron might vary slightly as a function of environment.