MUSE Data Cube Research Articles

A MUSE/VLT spatially resolved study of the emission structure of Green Pea galaxies

Green Pea galaxies (GPs) present among the most intense starbursts known in the nearby Universe. These galaxies are regarded as local analogs of high-redshift galaxies, making them a benchmark in the understanding of the star formation processes and the galactic evolution in the early Universe. In this work, we performed an integral field spectroscopic (IFS) study for a set of 24 GPs to investigate the interplay between its ionized interstellar medium (ISM) and the massive star formation that these galaxies present. Observations were taken in the optical spectral range (λ4750 Å–λ9350 Å) with the MUSE spectrograph attached to the 8.2 m telescope VLT. Spatial extension criteria were employed to verify which GPs are spatially resolved in the MUSE data cubes. We created and analyzed maps of spatially distributed emission lines (at different stages of excitation), continuum emission, and properties of the ionized ISM (e.g., ionization structure indicators, physical-chemical conditions, dust extinction). We also took advantage of our IFS data to produce integrated spectra of selected galactic regions in order to study their physical-chemical conditions. Maps of relevant emission lines and emission line ratios show that higher-excitation gas is preferentially located in the center of the galaxy, where the starburst is present. The continuum maps, with an average angular extent of 4″, exhibit more complex structures than the emission line maps. However, the [O III]λ5007 Å emission line maps tend to extend beyond the continuum images (the average angular extent is 5.5″), indicating the presence of low surface brightness ionized gas in the outer parts of the galaxies. Hα/Hβ, [S II]/Hα, and [O I]/Hα maps trace low-extinction, optically thin regions. The line ratios [O III]/Hβ and [N II]/Hα span extensive ranges, with values varying from 0.5 dex to 0.9 dex and from −1.7 dex to −0.8 dex, respectively. Regarding the integrated spectra, the line ratios were fit to derive physical properties including the electron densities ne = 30 − 530 cm−3, and, in six GPs with a measurable [O III]λ4363 Å line, electron temperatures of Te = 11 500 K–15 500 K, so the direct method was applied in these objects to retrieve metallicities 12 + log(O/H)≃8. We found the presence of the high-ionizing nebular He IIλ4686 Å line in three GPs, where two of them present among the highest sSFR values (> 8 × 108 yr−1) in this sample. Non-Wolf-Rayet (WR) features are detected in these galaxy spectra.

The morphology of the HD 163296 jet as a window on its planetary system

Context. HD 163296 is a Herbig Ae star which drives a bipolar knotty jet with a total length of ~6000 au. Strong evidence exists that the disk of HD 163296 harbors planets. Studies have shown that the presence of companions around jet-driving stars could affect the morphology of the jets. This includes a ‘wiggling’ of the jet axis and a periodicity in the positions of the jet knots. Aims. In this study we investigate the morphology (including the jet width and axis position) and proper motions of the HD 163296 jets, and use our results to better understand the whole system. Methods. This study is based on optical integral-field spectroscopy observations obtained with VLT/MUSE in 2017. Using spectro-images and position velocity diagrams extracted from the MUSE data cube, we investigated the number and positions of the jet knots. A comparison was made to X-shooter data collected in 2012 and the knot proper motions were estimated. The jet width and jet axis position with distance from the star were studied from the extracted spectro-images. This was done using Gaussian fitting to spatial profiles of the jet emission extracted perpendicular to the position angle of the jet. The centroid of the fit is taken as the position of the jet axis. Results. We observe the merging of knots and identify two previously undetected knots. We find proper motions that are broadly in agreement with previous studies. The jet width increases with distance from the source and we measure an opening angle of ~5° and 2.5° for the red and blue lobes, respectively. Measurements of the jet axis position, derived from Gaussian centroids of transverse intensity profiles, reveal a similar pattern of deviation in all forbidden emission lines along the first 20 arcsec of the jets. This result is interpreted as being due to asymmetric shocks and not due to a wiggling of the jet axis. Conclusions. The number of new knots detected and their positions challenge the 16-yr knot ejection periodicity proposed in prior studies, arguing for a more complicated jet system than was previously assumed. We use the non-detection of a jet axis wiggling to rule out companions with a mass >0.1 M⊙ and orbits between 1 au and 35 au. Any object inferred at these distances using other methods must be a brown dwarf or planet, otherwise it would have impacted the jet axis position. Both the precession and orbital motion scenarios are considered. Overall it is concluded that it is difficult to detect planets with orbits >1 au through a study of the jet axis.

SATELLITE: Application to Planetary Nebulae IFU Data

The integral field unit (IFU) spectroscopic view of extended ionized nebulae, such as planetary nebulae (PNe), H II regions, and galaxies, has changed the approach of studying these objects, providing a simultaneous characterization in both spatial directions. However, the spatial spaxel-by-spaxel analysis of such nebulae through IFUs is not directly comparable with the results obtained from the traditional slit-aperture spectroscopy or the predictions from 1D modelling. The new Python software called “satellite: Spectroscopic Analysis Tool for intEgraL fieLd unIt daTacubEs” is used in the analysis of the VIMOS and MUSE datacubes of four Galactic PNe. The 2D analysis of line ratio maps has shown important variations from one to another nebular component in NGC 7009 and NGC 6778. In particular, the knots in both PNe are characterized by strong emission from neutral gas that is weak or even absent from the main nebula, indicating significant variation in the ionization state and density structure among the nebular components. The far-UV radiation from the central star results in the photo-evaporation of the dense molecular knots resembling the spectrum of photodissociation regions.

FSR 1776: A new globular cluster in the Galactic bulge?

Context. Recent near-IR surveys have uncovered a plethora of new globular cluster (GC) candidates towards the Milky Way bulge. These new candidates need to be confirmed as real GCs and properly characterised. Aims. We investigate the physical nature of FSR 1776, a very interesting star cluster projected towards the Galactic bulge. This object was originally classified as an intermediate-age open cluster, and has recently been re-discovered independently and classified as a GC candidate (Minni 23). Firstly, our aim is to confirm its GC nature; secondly, we determine its physical parameters. Methods. The confirmation of the cluster existence is checked using the radial velocity (RV) distribution of a MUSE data cube centred at FSR 1776. The cluster parameters are derived from isochrone fitting to the RV-cleaned colour-magnitude diagrams (CMDs) from visible and near-IR photometry taken from VVV, 2MASS, DECAPS, and Gaia all together. Results. The predicted RV distribution for the FSR 1776 coordinates, considering only contributions from the bulge and disc field stars, is not enough to explain the observed MUSE RV distribution. The extra population (12% of the sample) is FSR 1776 with an average RV of −103.7 ± 0.4 km s−1. The CMDs reveal that it is 10 ± 1 Gyr metal-rich population with [Fe/H]phot ≈ +0.2 ± 0.2, [Fe/H]spec = +0.02 ± 0.01(σ = 0.14 dex), located at the bulge distance of 7.24 ± 0.5 kpc with AV ≈ 1.1 mag. The mean cluster proper motions are (⟨μα⟩,⟨μδ⟩) = (−2.3 ± 1.1, −2.6 ± 0.8) mas yr−1. Conclusions. FSR 1776 is an old GC located in the Galactic bulge with a super-solar metallicity, among the highest for a Galactic GC. This is consistent with predictions for the age–metallicity relation of the bulge, being FSR 1776 the probable missing link between typical GCs and the metal-rich bulge field. High-resolution spectroscopy of a larger field of view and deeper CMDs are now required for a full characterisation.

Recovery and analysis of rest-frame UV emission lines in 2052 galaxies observed with MUSE at 1.5 <z< 6.4

Rest-frame ultraviolet (UV) emission lines probe electron densities, gas-phase abundances, metallicities, and ionization parameters of the emitting star-forming galaxies and their environments. The strongest main UV emission line, Lyα, has been instrumental in advancing the general knowledge of galaxy formation in the early universe. However, observing Lyαemission becomes increasingly challenging atz ≳ 6 when the neutral hydrogen fraction of the circumgalactic and intergalactic media increases. Secondary weaker UV emission lines provide important alternative methods for studying galaxy properties at high redshift. We present a large sample of rest-frame UV emission line sources at intermediate redshift for calibrating and exploring the connection between secondary UV lines and the emitting galaxies’ physical properties and their Lyαemission. The sample of 2052 emission line sources with 1.5 < z < 6.4 was collected from integral field data from the MUSE-Wide and MUSE-Deep surveys taken as part of Guaranteed Time Observations. The objects were selected through untargeted source detection (i.e., no preselection of sources as in dedicated spectroscopic campaigns) in the three-dimensional MUSE data cubes. We searched optimally extracted one-dimensional spectra of the full sample for UV emission features via emission line template matching, resulting in a sample of more than 100 rest-frame UV emission line detections. We show that the detection efficiency of (non-Lyα) UV emission lines increases with survey depth, and that the emission line strength of He IIλ1640 Å, [O III]λ1661 + O III]λ1666, and [Si III]λ1883 + Si III]λ1892 correlate with the strength of [C III]λ1907 + C III]λ1909. The rest-frame equivalent width (EW0) of [C III]λ1907 + C III]λ1909 is found to be roughly 0.22 ± 0.18 of EW0(Lyα). We measured the velocity offsets of resonant emission lines with respect to systemic tracers. For C IVλ1548 + C IVλ1551 we find that ΔvC IV≲ 250 km s−1, whereas ΔvLyαfalls in the range of 250−500 km s−1which is in agreement with previous results from the literature. The electron densitynemeasured from [Si III]λ1883 + Si III]λ1892 and [C III]λ1907 + C III]λ1909 line flux ratios is generally < 105cm−3and the gas-phase abundance is below solar at 12 + log10(O/H)≈8. Lastly, we used “PhotoIonization Model Probability Density Functions” to infer physical parameters of the full sample and individual systems based on photoionization model parameter grids and observational constraints from our UV emission line searches. This reveals that the UV line emitters generally have ionization parameter log10(U) ≈ −2.5 and metal mass fractions that scatter aroundZ ≈ 10−2, that isZ ≈ 0.66 Z⊙. Value-added catalogs of the full sample of MUSE objects studied in this work and a collection of UV line emitters from the literature are provided with this paper.

Physics of ULIRGs with MUSE and ALMA: The PUMA project

Context. Ultraluminous infrared galaxies (ULIRGs) are characterised by extreme starburst (SB) and active galactic nucleus (AGN) activity, and are therefore ideal laboratories for studying the outflow phenomena and their feedback effects. We have recently started a project called Physics of ULIRGs with MUSE and ALMA (PUMA), which is a survey of 25 nearby (z < 0.165) ULIRGs observed with the integral field spectrograph MUSE and the interferometer ALMA. This sample includes systems with both AGN and SB nuclear activity in the pre- and post-coalescence phases of major mergers. Aims. The main goals of the project are (i) to study the prevalence of (ionised, neutral, and molecular) outflows as a function of the galaxy properties, (ii) to constrain the driving mechanisms of the outflows (e.g. distinguish between SB and AGN winds), and (iii) to identify and characterise feedback effects on the host galaxy. In this first paper, we present details on the sample selection, MUSE observations, and data reduction, and derive first high-level data products. Methods. MUSE data cubes were analysed to study the dynamical status of each of the 21 ULIRGs observed so far, taking the stellar kinematics and the morphological properties inferred from MUSE narrow-band images into account. We also located the ULIRG nuclei, taking advantage of near-infrared (HST) and millimeter (ALMA) data, and studied their optical spectra to infer (i) the ionisation state through standard optical line ratio diagnostics, and (ii) outflows in both atomic ionised ([O III], Hα) and neutral (Na ID) gas. Results. We show that the morphological and stellar kinematic classifications are consistent: post-coalescence systems are more likely associated with ordered motions, while interacting (binary) systems are dominated by non-ordered and streaming motions. We also find broad and asymmetric [O III] and Na ID profiles in almost all nuclear spectra, with line widths in the range [300 − 2000] km s−1, possibly associated with AGN- and SB-driven winds. This result reinforces previous findings that indicated that outflows are ubiquitous during the pre- and post-coalescence phases of major mergers.

The ultra-diffuse galaxy NGC 1052-DF2 with MUSE

The so-called ultra-diffuse galaxy NGC 1052-DF2 was announced to be a galaxy lacking dark matter based on a spectroscopic study of its constituent globular clusters. Here we present the first spectroscopic analysis of the stellar body of this galaxy using the MUSE integral-field spectrograph at the (ESO) Very Large Telescope. The MUSE datacube simultaneously provides DF2’s stellar velocity field and systemic velocities for seven globular clusters (GCs). We further discovered three planetary nebulae (PNe) that are likely part of this galaxy. While five of the clusters had velocities measured in the literature, we were able to confirm the membership of two more candidates through precise radial velocity measurements, which increases the measured specific frequency of GCs in DF2. The mean velocity of the diffuse stellar body, 1792.9+1.4−1.8 km s−1, is consistent with the mean globular cluster velocity. We detect a weak but significant velocity gradient within the stellar body, with a kinematic axis close to the photometric major axis, making it a prolate-like rotator. We estimate a velocity dispersion from the clusters and PNe of σint = 10.6−2.3+3.9 km s−1. The velocity dispersion σDF2⋆(Re) for the stellar body within one effective radius is 10.8+3.2−4.0 km s−1. Considering various sources of systemic uncertainties, this central value varies between 5 and 13 km s−1, and we conservatively report a 95% confidence upper limit to the dispersion within one Re of 21 km s−1. We provide updated mass estimates based on these dispersions corresponding to the different distances to NGC 1052-DF2 that have been reported in the recent literature.

The environment of the SN-less GRB 111005A at z = 0.0133

The collapsar model has proved highly successful in explaining the properties of long γ-ray bursts (GRBs), with the most direct confirmation being the detection of a supernova (SN) coincident with the majority of nearby long GRBs. Within this model, a long GRB is produced by the core-collapse of a metal-poor, rapidly rotating, massive star. The detection of some long GRBs in metal-rich environments, and more fundamentally the three examples of long GRBs (GRB 060505, GRB 060614 and GRB 111005A) with no coincident SN detection down to very deep limits is in strong contention with theoretical expectations. In this paper we present MUSE observations of the host galaxy of GRB 111005A, which is the most recent and compelling example yet of a SN-less, long GRB. At z = 0.01326, GRB 111005A is the third closest GRB ever detected, and second closest long duration GRB, enabling the nearby environment to be studied at a resolution of 270 pc. From the analysis of the MUSE data cube, we find GRB 111005A to have occurred within a metal-rich environment with little signs of ongoing star formation. Spectral analysis at the position of the GRB indicates the presence of an old stellar population (τ ≥ 10 Myr), which limits the mass of the GRB progenitor to MZAMS < 15 M⊙, in direct conflict with the collapsar model. Our deep limits on the presence of any SN emission combined with the environmental conditions at the position of GRB 111005A necessitate the exploration of a novel long GRB formation mechanism that is unrelated to massive stars.

Mapping the core of the Tarantula Nebula with VLT-MUSE

We introduce VLT-MUSE observations of the central 2′× 2′ (30 × 30 pc) of the Tarantula Nebula in the Large Magellanic Cloud. The observations provide an unprecedented spectroscopic census of the massive stars and ionised gas in the vicinity of R136, the young, dense star cluster located in NGC 2070, at the heart of the richest star-forming region in the Local Group. Spectrophotometry and radial-velocity estimates of the nebular gas (superimposed on the stellar spectra) are provided for 2255 point sources extracted from the MUSE datacubes, and we present estimates of stellar radial velocities for 270 early-type stars (finding an average systemic velocity of 271 ± 41 km s−1). We present an extinction map constructed from the nebular Balmer lines, with electron densities and temperatures estimated from intensity ratios of the [S II], [N II], and [S III] lines. The interstellar medium, as traced by Hα and [N II] λ6583, provides new insights in regions where stars are probably forming. The gas kinematics are complex, but with a clear bi-modal, blue- and red-shifted distribution compared to the systemic velocity of the gas centred on R136. Interesting point-like sources are also seen in the eastern cavity, western shell, and around R136; these might be related to phenomena such as runaway stars, jets, formation of new stars, or the interaction of the gas with the population of Wolf–Rayet stars. Closer inspection of the core reveals red-shifted material surrounding the strongest X-ray sources, although we are unable to investigate the kinematics in detail as the stars are spatially unresolved in the MUSE data. Further papers in this series will discuss the detailed stellar content of NGC 2070 and its integrated stellar and nebular properties.

Deep MUSE observations in the HDFS

(Abridged) We make use of the deepest VLT/MUSE observations performed so far on the Hubble Deep Field South (HDFS) to characterize the low-mass (< $10^{10}$M$_\odot$) galaxy population at intermediate redshift. We identify a sample of 28 spatially-resolved emission-line galaxies in the deep (27h integration time) MUSE data cube, spread over a redshift interval of 0.2 < z < 1.4. The public HST images and multi-band photometry over the HDFS are used to constrain the stellar mass and star formation rate (SFR) of the galaxies and to perform a morphological analysis. We derive the resolved ionized gas properties of these galaxies from the MUSE data and model the disk (both in 2D and with GalPaK$^{\rm 3D}$) to retrieve their intrinsic gas kinematics. We build a sample of resolved emission-line galaxies of much lower stellar mass and SFR (by $\sim$1-2 orders of magnitude) than previous 3D spectroscopic surveys. Most of the spatially-resolved MUSE-HDFS galaxies have gas kinematics consistent with disk-like rotation, but about 20% have velocity dispersions larger than the rotation velocities, and 30% are part of a close pair and/or show clear signs of recent gravitational interactions. In the high-mass regime, the MUSE-HDFS galaxies follow the Tully-Fisher relation defined from previous surveys in a similar redshift range. This scaling relation extends also to lower masses/velocities but with a higher dispersion. The MUSE-HDFS galaxies follow the scaling relations defined in the local universe between the specific angular momentum and the stellar mass. However, we find that intermediate-redshift star-forming galaxies fill a continuum transition from the spiral to elliptical local scaling relations, according to the dynamical state of the gas. This indicates that some galaxies may lose their angular momentum and become dispersion-dominated prior to becoming passive.

The mysterious morphology of MRC0943-242 as revealed by ALMA and MUSE

We present a pilot study of the z=2.923 radio galaxy MRC0943-242, where we for the first time combine information from ALMA and MUSE data cubes. Even with modest integration times, we disentangle an AGN and a starburst dominated set of components. These data reveal a highly complex morphology, as the AGN, starburst, and molecular gas components show up as widely separated sources in dust continuum, optical continuum and CO line emission observations. CO(1-0) and CO(8-7) line emission suggest that there is a molecular gas reservoir offset from both the dust and the optical continuum that is located ~90kpc from the AGN. The UV line emission has a complex structure in emission and absorption. The line emission is mostly due to i) a large scale ionisation cone energised by the AGN, ii) a Ly-alpha emitting bridge of gas between the radio galaxy and a heavily star-forming set of components. Strangely, the ionisation cone has no Ly-alpha emission. We find this is due to an optically thick layer of neutral gas with unity covering fraction spread out over a region of at least ~100kpc from the AGN. Other, less thick absorption components are associated with Ly-alpha emitting gas within a few tens of kpc from the radio galaxy and are connected by a bridge of emission. We speculate that this linear structure of dust, Ly-alpha and CO emission, and the redshifted absorption seen in the circum-nuclear region may represent an accretion flow feeding gas into this massive AGN host galaxy.

The MUSE 3D view of theHubbleDeep Field South

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

Processing MUSE hyperspectral data: Denoising, deconvolution and detection of astrophysical sources

We address the problem of processing astronomical hyperspectral data cubes in the context of the forthcoming MUSE instrument. MUSE, which is under construction, will provide massive hyperspectral data with about 300 × 300 pixels at approximately 4000 wavelengths. One of its main astrophysical objectives concerns the observation of extragalactic deep fields, where MUSE should be able to detect and characterize galaxies much fainter than the ones currently observed by other ground-based instruments. The data will suffer, however, from very powerful and spectrally variable perturbations. In this paper, MUSE data cubes are first considered as a collection of spectra, which are processed independently. A restoration method is proposed, based on the hypothesis that data can be approximated by appropriate sparse representations. Sparsity can be naturally expressed in the spectral domain, where a galaxy spectrum is mainly the superposition of an emission and absorption line spectrum, which is naturally sparse, on a continuum, which is supposed to have a sparse discrete cosine transform. The problem is addressed within the ℓ 1 -norm penalization setting. The original features of the model consist, first, in taking into account observational specificities such as the spectrally variable instrumental response and non-identically distributed noise, and, second, in tuning regularization parameters in this setting, which are fixed in order to obtain uniform false alarm rates for decomposition coefficients. In a second step, such sparse decompositions are used as an input to an object detection and characterization method. The decomposed spectra are first used for spatial segmentation. Then, once a group of pixels has been identified as belonging to the same object, the corresponding spectrum and amplitude map are jointly estimated under the former sparsity assumption. Applications to object identification, the amplitude map and spectrum estimation are presented for realistic deep field simulated data cubes provided by the MUSE consortium.