Issue 
A&A
Volume 654, October 2021



Article Number  A49  
Number of page(s)  21  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/202040225  
Published online  08 October 2021 
The MUSE Hubble Ultra Deep Field Survey
XVI. The angular momentum of lowmass starforming galaxies: A cautionary tale and insights from TNG50
^{1}
Univ Lyon, Univ Lyon1, ENS de Lyon, CNRS, Centre de Recherche Astrophysique de Lyon (CRAL) UMR5574, 69230 SaintGenisLaval, France
email: nicolas.bouche@univlyon1.fr
^{2}
Center for Computational Astrophysics, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010, USA
^{3}
Columbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, New York, NY 10027, USA
^{4}
Institut de Recherche en Astrophysique et Planétologie (IRAP), Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
^{5}
Laboratoire Lagrange, Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Blvd de l’Observatoire, 06304 Nice cedex 4, France
^{6}
Aix Marseille Univ, CNRS, CNES, LAM, Marseille, France
^{7}
MaxPlanckInstitut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany
^{8}
LeibnizInstitut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Postdam, Germany
^{9}
Universität Heidelberg, Zentrum für Astronomie, Institut für theoretische Astrophysik, AlbertUeberleStr. 2, 69120 Heidelberg, Germany
^{10}
Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
^{11}
Department of Physics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
^{12}
Department of Astronomy, University of WisconsinMadison, 475 N. Charter Street, Madison, WI 53706, USA
Received:
23
December
2020
Accepted:
19
June
2021
We investigate the specific angular momentum (sAM) j(< r) profiles of intermediate redshift (0.4 < z < 1.4) starforming galaxies (SFGs) in the relatively unexplored regime of low masses (down to M_{⋆} ∼ 10^{8} M_{⊙}) and small sizes (down to R_{e} ∼ 1.5 kpc), and we characterize the sAM scaling relation (i.e., Fall relation) and its redshift evolution. We have developed a 3D methodology to constrain sAM profiles of the starforming gas using a forward modeling approach with GAlPAK^{3D} that incorporates the effects of beam smearing, yielding the intrinsic morphokinematic properties even with limited spatial resolution data. Using mock observations from the TNG50 simulation, we find that our 3D methodology robustly recovers the star formation rate (SFR)weighted j̃_{⋆}(<_{r}) profiles down to a low effective signaltonoise ratio of ⪆3. We applied our methodology blindly to a sample of 494 [O II]selected SFGs in the MUSE Ultra Deep Field (UDF) 9 arcmin^{2} mosaic data, covering the unexplored 8 < log M_{⋆}/M_{⊙} < 9 mass range. We find that the (SFRweighted) sAM relation follows j̃_{⋆} ∝ M_{⋆}^{α} with an index α varying from α = 0.3 to α = 0.5, from log M_{⋆}/M_{⊙} = 8 to log M_{⋆}/M_{⊙} = 10.5. The UDF sample supports a redshift evolution j̃_{⋆} ∝(1+z)^{a}, with a = −0.27_{−0.56}^{+0.42} which is consistent with the (1 + z)^{−0.5} expectation from a universe in expansion. The scatter of the sAM sequence is a strong function of the dynamical state with logj_{M⋆} ∝ 0.65_{−0.08}^{+0.06} × log(V_{max}/σ), where σ is the velocity dispersion at 2R_{e}. In TNG50, SFGs also form a j̃_{⋆}−M_{⋆}−(V/σ) plane, but it correlates more with galaxy size than with morphological parameters. Our results suggest that SFGs might experience a dynamical transformation, and lose their sAM, before their morphological transformation to becoming passive via either merging or secular evolution.
Key words: galaxies: highredshift / galaxies: evolution / galaxies: kinematics and dynamics / galaxies: structure
© N. F. Bouché et al. 2021
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
In a Λ cold dark matter (ΛCDM) universe, baryons cool, fall inwards, and form centrifugally supported disks in the centers of halos. The specific angular momentum (sAM) sets the pressure, instabilities, gas fractions (Obreschkow et al. 2016; Romeo & Mogotsi 2018; Romeo 2020; Li et al. 2020), and most importantly determines the disk properties, such as size (e.g., White & Rees 1978; Fall & Efstathiou 1980; Fall 1983; Mo et al. 1998; Dalcanton et al. 1997; van den Bosch et al. 2003; Dutton & van den Bosch 2009; Somerville et al. 2018). As disks evolve from z = 2 to the present, they must grow from a vast reservoir of corotating cold accreting material in the circumgalactic medium (CGM) as argued in Renzini (2020). This is also strongly supported by hydrodynamical simulations which predict roughly coplanar gaseous structures (Stewart et al. 2013, 2017; Danovich et al. 2015; Ho et al. 2019; Kretschmer et al. 2020; DeFelippis et al. 2021) embedded in a rotating CGM (DeFelippis et al. 2020). These coplanar structures were initially found by Barcons et al. (1995) and Steidel et al. (2002) in a handful of background quasar sightlines, but they are now routinely observed on scales of 20–80 kpc (as in Bouché et al. 2013, 2016; Ho et al. 2017; Lopez et al. 2018; Martin et al. 2019; Zabl et al. 2019). These structures are important as they not only bring fuel for starformation but also angular momentum (Renzini 2020; DeFelippis et al. 2021).
In this context, one of the most fundamental properties of disk galaxies is their sAM (as argued in Fall & Efstathiou 1980). The sAM of halos (j_{h} = J_{h}/M_{h}) is tightly correlated to the halo mass M_{h} () as a result of tidal torques from the largescale structure (e.g., Peebles 1969; Efstathiou & Jones 1979). It is often assumed that the sAM of baryons and of dark matter are equal at accretion, and that as disks form in the centers of halos the baryonic sAM is conserved (e.g., Fall & Efstathiou 1980; Mo et al. 1998; Burkert et al. 2010). While this assumption of j conservation during collapse has been appealing since the late 1970s (e.g., White & Rees 1978; Fall & Efstathiou 1980), it is not trivial given the infalling baryons can both lose and gain angular momentum from the virial radius to the inner disk (e.g., Navarro & Steinmetz 2000; Sharma et al. 2012; Danovich et al. 2015; Genel et al. 2015; DeFelippis et al. 2017; Jiang et al. 2019). In particular, Genel et al. (2015) and DeFelippis et al. (2017) found that galactic winds are essential for producing latetype galaxies with sufficient angular momentum.
Observationally, the stellar sAM of massive disks (with log M_{⋆}/M_{⊙}>9) follows a similar scaling relation as the dark matter j_{h} − M_{h} scaling relation with with α ≃ 0.6, as reported by M. Fall and collaborators (Fall 1983; Romanowsky & Fall 2012; Fall & Romanowsky 2013, hereafter FR13) and recently revisited by Fall & Romanowsky (2018, hereafter FR18). This similarity between the halo and baryonic sAM scaling relations, both in normalization (for galaxy disks) and slope, highlights a potential intimate link between galaxies and their host DM haloes, as discussed extensively in Posti et al. (2018b). However, hydrodynamical simulations tend to show a more complex relation between the DM and baryonic sAM (e.g., Genel et al. 2015; Jiang et al. 2019).
Since the pioneering study of Fall (1983), the j_{⋆} − M_{⋆} relation has become a standard tool in galaxy evolution research thanks to recent Integral Field Spectrograph (IFS) surveys and to more accurate simulations of disk formation (e.g., Vogelsberger et al. 2012, 2013; Schaye et al. 2015; Ceverino et al. 2017; Grand et al. 2019; Nelson et al. 2019; Pillepich et al. 2019), see review in Vogelsberger et al. (2020). Hence, understanding the slope and normalization of this important scaling relation is of paramount importance as it is intimately related to galaxy morphological diversity and holds clues to accretion of baryons onto galaxies (Genel et al. 2015; Teklu et al. 2015; Zavala et al. 2016; Lagos et al. 2017; ElBadry et al. 2018).
In the local universe, several groups have expanded on the work of Fall and collaborators to confirm the j_{⋆} − M_{⋆} relation with increasingly large samples of galaxies as in Cortese et al. (2016) and Lapi et al. (2018) which each had ∼500 galaxies, albeit mostly with M_{⋆} > 10^{9.5} M_{⊙}. Some of these local surveys have shown that the scatter of the j_{⋆} − M_{⋆} relation strongly correlates with morphology (Obreschkow & Glazebrook 2014; Cortese et al. 2016; Rizzo et al. 2018; Sweet et al. 2018). At lower masses M_{⋆} < 10^{9} M_{⊙}, the j_{⋆} − M_{⋆} relation is not understood as well either theoretically (e.g., Stevens et al. 2016; Mitchell et al. 2018; ElBadry et al. 2018) or observationally (but see Butler et al. 2017; Posti et al. 2018a; Mancera Piña et al. 2021). Both the sample of 14 dIrr galaxies in the 10^{6}–10^{9} M_{⊙} mass range from Butler et al. (2017) and the local disks with M_{⋆} > 10^{7} M_{⊙} from Posti et al. (2018a) and Mancera Piña et al. (2021) show that the j_{⋆} − M_{⋆} relation is well described by a single powerlaw at all scales with a slope α = 0.6 ± 0.1 (see Fall & Romanowsky 2018).
At high redshifts, studying the j − M relation has become particularly important in order to constrain the redshift evolution of this fundamental scaling relation. Thanks to recent advancements with second generation IFS instruments, such as the MultiUnit Spectroscopic Explorer (MUSE; Bacon et al. 2010) and the Kband MultiObject Spectrograph (KMOS; Sharples et al. 2013), several groups have constrained the j − M relation using increasingly larger samples of z > 1 starforming galaxies (SFGs). Among these, there is the sample of 360 SFGs from the SINS/KMOS3D survey (Burkert et al. 2010, 2016; Wisnioski et al. 2019) at redshifts 0.8 < z < 2.3, the sample of ≈400 SFGs at z ∼ 0.8 of Swinbank et al. (2017) collected in MUSE and KMOS surveys, and the sample of 486 z = 0.6−1.0 SFGs (Harrison et al. 2017) from the KROSS Hα survey (Stott et al. 2016). These surveys (and others, e.g. Alcorn et al. 2018; Gillman et al. 2020) are often limited to relatively massive SFGs with a typical limit of M_{⋆} > 10^{9.5} M_{⊙}.
The j − M relation has also not been explored in the low mass regime at M_{⋆} of 10^{8} M_{⊙}–10^{9.5} M_{⊙} beyond the local universe (Posti et al. 2018a). Thanks to the MUSE IFS (Bacon et al. 2010) and its exquisite sensitivity, it is now possible to study the kinematics of SFGs at low masses down to M_{⋆} = 10^{8} M_{⊙}, as first demonstrated by Contini et al. (2016). The lowmass regime (M_{⋆} = 10^{8−9} M_{⊙}) is challenging as galaxies become marginally resolved (when smaller than 3 kpc), and it becomes difficult to recover reliable measurements of both sizes and rotation velocity.
This paper aims to show that the j − M relation can be characterized in small and lowmass galaxies at redshifts from z = 0.4–1.4, opening up new parameter space thanks to the MUSE instrument and to our 3D modeling approach. Contrary to galaxies in the local universe where the j − M relation can be studied independently for the stellar j_{⋆} or gas components j_{gas} (e.g., Obreschkow & Glazebrook 2014; Obreschkow et al. 2016; Li et al. 2020; Mancera Piña et al. 2021), this is currently extremely challenging in distant galaxies at z ≃ 1.0 or beyond. One should note that the gas sAM j_{gas} in this context is derived solely from H I data at z = 0, and it is often assumed that the sAM for the starforming gas (molecular or ionized) is similar to j_{⋆} because of the similar extent (Leroy et al. 2008; Nelson et al. 2016; Wilman et al. 2020) and kinematics (Martinsson et al. 2013; Lang et al. 2020)^{1}. As discussed in the next section, we focus our analysis on the sAM j derived from the starforming gas, denoted , both for a large sample of 494 SFGs from the MUSE Hubble Ultra Deep Field (UDF) observations presented in Bacon et al. (2017) and for SFGs from the TNG50 simulations (Nelson et al. 2019; Pillepich et al. 2019).
This paper is organized as follows. In Sect. 2, we present our methodology that is designed to be more accurate in the lowmass regime. In Sect. 3, we validate our methodology on 158 simulated galaxies taken from TNG50. In Sect. 4, we apply our methodology on 494 galaxies in the MUSE UDF. In Sect. 5, we compare our results to ∼5000 SFGs in TNG50. We discuss the implication of our results in Sect. 6. Finally, we present our conclusions in Sect. 7. Throughout this paper, we use a “Planck 2015” cosmology (Planck Collaboration XIII 2016) with Ω_{M} = 0.307, Λ = 0.693, H_{0} = 67 km s^{−1} Mpc^{−1}, yielding 8.23 kpc arcsec^{−1}, and we consistently use “log” for the base10 logarithm.
2. Measuring angular momentum
2.1. Definitions
The sAM j ≡ J/M is, in general, given by
In practice, the twodimensional version of the specific angular momentum, appropriate for an axisymmetric disk, is used
where R is the enclosed radius, Σ(r) the surface brightness, and v(r) the rotation velocity profile.
Regarding Eq. (2), the total angular momentum j_{tot} ≡ j(< ∞) is not observable directly as inherent surface brightness limitations preclude measurements beyond R_{e} or at best 2R_{e} (e.g., Marasco et al. 2019). For massive galaxies, probing j up to R_{e} or 2R_{e}, where the j profile reaches its plateau, might be sufficient, but in low mass galaxies, that regime might be difficult to reach even in the local universe, as discussed in Posti et al. (2018a), Marasco et al. (2019). Thus, previous work often resorted to estimators for j_{tot} such as Eq. (3), which is discussed in Sect. 2.2. Here, we use a 3D methodology (Sect. 2.3) that bypasses common approximations and similar in essence to the analysis of Mancera Piña et al. (2021) who used ^{3D}Barolo (Di Teodoro & Fraternali 2015). In Sect. 3, we show with mock cubes from TNG50 that our 3D method combined with Eq. (2) does reproduce the true j values, computed directly from Eq. (1).
It is worth noting that in order to determine j_{⋆}(r) from Eq. (2), one needs the stellar mass profile Σ_{⋆}(r) and the stellar rotation curve v_{⋆}(r). However, determining stellar kinematics beyond z = 0 is challenging in SFGs (but see Guérou et al. 2017, for a first attempt) or even in passive galaxies with strong continuum (e.g., van de Sande et al. 2013; Bezanson et al. 2018; Cole et al. 2020). As a result, most studies of the j_{⋆} − M_{⋆} relation use the starforming gas v_{Hα} kinematics as a proxy for v_{⋆} (e.g., Burkert et al. 2016; Harrison et al. 2017; Swinbank et al. 2017). Two studies attempted to convert v_{Hα} to v_{⋆}, Marasco et al. (2019) at z = 1 and Posti et al. (2018a) at z = 0 using the asymmetric drift correction (Meurer et al. 1996; Dalcanton & Stilp 2010; Burkert et al. 2010).
Regarding the stellar mass profile Σ_{⋆}(r), it acts as the normalization factor in Eq. (2), and thus one only need to know the profile shape, namely the Sersic index and the halflight radius R_{e}. For SFGs, one can use the starforming gas profile Σ_{Hα} (or Σ_{[O II]}) as a proxy for Σ_{⋆} provided that the SF gas and stars have similar scale length^{2}, which is expected from the almost linear KennicutSchmidt relation (Kennicutt 1998; Genzel et al. 2010; Daddi et al. 2010; Bacchini et al. 2019). In fact, the Hα size R_{e, Hα} to continuum size ratio R_{e, ⋆} has been found to be ≈1^{3} from a very large sample of 3200 galaxies at 0.7 < z < 1.5 with WFC3 imaging in Nelson et al. (2016), confirmed in our data as discussed in Sect. 4.2.
In this paper, we measure the star formation rate (SFR)weighted sAM j_{sfr} determined from the starforming gas properties, namely v_{[O II]}(r) and Σ_{[O II}(r), as a proxy for j_{⋆} following other studies (e.g., Burkert et al. 2016; Contini et al. 2016; Swinbank et al. 2017). This has the advantage that the measurement is selfcontained (from Hα or [O II]) and that a comparison to similar quantities from simulations can be performed. We use the notation to designate j_{sfr} and in Sect. 3.4, we quantify the potential bias between and j_{⋆} in the TNG50 galaxies.
2.2. 2D Methodology: An Approximate sAM
Given the difficulties in measuring j(< r) in the outskirts of galaxies, a method that is commonly used in the literature (e.g., Burkert et al. 2016; Swinbank et al. 2017; Harrison et al. 2017; Gillman et al. 2020) for estimating j_{tot} is the following formulae (Romanowsky & Fall 2012, hereafter RF12):
where R_{e} is the halflight radius, v_{s} the rotation velocity at 2R_{e}, v(2 R_{e}), and k_{n} = 1.15 + 0.029 ⋅ n + 0.062 ⋅ n^{2} is a constant that depends on the Sérsic (1963) index n. Equation (3) reduces to the wellknown j_{tot, RF} = 2 v_{s} R_{d} for disks with R_{d} being the exponential scale length.
However, we caution the reader against using Eq. (3) for low mass galaxies with log M_{⋆}/M_{⊙}) < 9.5–10, for the following reasons. Equation (3) assumes that the velocity profile v(r) is a constant from r = 0, which might be appropriate for galaxies with log M/M_{⋆}) > 10, but is less suitable for intermediate and low mass galaxies, which are known to have slowly rising rotation curves (e.g., Persic et al. 1996; Catinella et al. 2006). Consequently, Eq. (3) can lead to an overestimation of the total angular momentum as discussed in Obreschkow & Glazebrook (2014).
To quantify the overestimation in the lowmass regime, we show in Fig. 1(left) several sAM profiles (solid lines), corresponding to “tanh” rotation curves, v(r) = V_{max}tanh(r/R_{t}) with varying inner slopes (parameterized by the turnover radius R_{t}), normalized to the RF12 approximation (represented by the blue horizontal line). This figure shows that, for galaxies with slowly rising rotation curves, that is with R_{t}/R_{e} of order unity, Eq. (3) results in an overestimation of the angular momentum by up to 15%. Obreschkow & Glazebrook (2014) made a similar point using rotation curves v(r)∝(1 − exp(−r/R_{t})).
Fig. 1.
Systematics at play in measuring j(< r). Left: sAM profile j(< r) normalized by the RF12 approximation (Eq. (3)) as a function of radius for disks with Sérsic n = 1 and a tanh rotation curve. The solid blue horizontal line represents the Eq. (3) approximation. The solid curves represent the sAM for various rotation curves with various turnover radii R_{t} where the opacity of the line is inversely proportional to R_{t}/R_{e}. The vertical dotted line represents 2R_{e}. For galaxies with slowly rising rotation curves, with larger R_{t}, the RF12 approximation leads to a large overestimation of the total angular momentum (represented by the horizontal dashed lines), by 10–20%. Right: relative difference for the k_{n} coefficient between the k_{n} formula from RF12 and Eq. (4) as a function of Sérsic index n. One sees that Eq. (3) under or overestimate k_{n} by up to ∼10% for 0.5 < n < 4. 
Moreover, there can be significant differences between the often used k_{n} approximation, namely k_{n} = 1.15 + 0.029 ⋅ n + 0.062 ⋅ n^{2} and the exact expression (derived in Appendix A) given by:
that can be as large as 15% for Sérsic indices n ≲ 1 as indicated in Fig. 1(right).
In conclusion, combining these two effects, Eqs. (3) and (4), can lead to an overestimation of the sAM (up to 20–25%) for galaxies with n ≲ 1 and/or when the rotation curve is slowly rising, that is in the lowmass regime at M_{⋆} < 10^{9.5} M_{⊙}. In addition, such galaxies with a large turnover radius R_{t} tend to have large sizes and hence a low surface brightness, and thus their true V_{max} is more difficult to measure beyond R_{t}.
2.3. 3D Methodology
In order to estimate j directly from Eq. (2), one needs the intrinsic (deprojected and free from instrumental and resolution effects) velocity and surface brightness flux profile and these can be measured directly from 3D (positionpositionvelocity) data using the GAlPAK^{3D} algorithm^{4} described in Bouché et al. (2015).
Briefly, GAlPAK^{3D} performs a parametric fit on emission line data for the morphology and kinematics simultaneously directly onto the 3D data using a 3D (x, y, λ) disk model which specifies the morphology and kinematic profiles. For the morphology, the model assumes a Sérsic (1963) surface brightness profile Σ(r), with Sérsic index n. For the kinematics, the model assumes a parametric form for the rotation curve v(r) and for the dispersion σ(r) profile. The model parameters are the centroid position x, y, z, the total line flux, the inclination i, the majoraxis position angle P.A, the halflight radius (R_{e}), the turnover radius R_{t} and maximum velocity V_{max} for the rotation curve, and the velocity dispersion σ_{0}, described below.
The rotation curve v(r) can be “arctan”, ‘tanh’ or other analytical forms as described in the documentation^{5}. As described in Bouché et al. (2015), the velocity dispersion profile σ(r) is made of three components added in quadrature. The first component σ_{1} is that of a thick disk (Genzel et al. 2008), namely σ_{1}(r)/v(r) = h_{z}/r where h_{z} is the disk thickness, taken to be 0.2 × R_{e}. The second component σ_{2} is comes from the natural broadning of the 3D disk model of finite thickness. The third component σ_{3} is an adjustable constant term σ_{0}, as in Genzel et al. (2008), Förster Schreiber et al. (2006), Cresci et al. (2009), Wisnioski et al. (2015), Übler et al. (2019) and describes any additional turbulent component to the kinematic dispersions.
GAlPAK^{3D} convolves the model with the Point Spread Function (PSF) and the instrumental Line Spread Function (LSF), which implies that the parameters are “intrinsic”, that is corrected for beam smearing and instrumental effects. Here, we use a Moffat PSF as characterized in Bacon et al. (2017), and we use Eq. (7) of Bacon et al. (2017) for the LSF.
GAlPAK^{3D} uses a simple MetropolisHasting (MH) algorithm to optimize the parameters whose posterior distributions are given by the chain posteriors. As described in Bouché et al. (2015), the GAlPAK^{3D} implementation of the MH uses a nonstandard Cauchy proposal distribution which shortens the burnin phase considerably, but suffers, like most MH algorithms, from the need to manually tune the width of the proposal distribution.
Since its initial release, the following improvements and new features have been implemented in GAlPAK^{3D}. First, in order to quantify the chain convergence for each parameter, GAlPAK^{3D} performs the Geweke (1992)’s convergence diagnostic test for each of the parameters, which is a Ztest of equality of means (of two parts of the chain) where autocorrelation in the samples is taken into account. Second, GAlPAK^{3D} can now adapt the widths of the Cauchy proposal distribution automatically. Third, and most importantly, GAlPAK^{3D} is now (since v1.20) compatible with the EMCEE code (ForemanMackey et al. 2013). EMCEE has the advantage that its multiwalkers remove the need to tune the proposal distribution, but it has the disadvantage to retain the Normal proposal distribution which is very slow and inefficient in reaching convergence because the Gaussian function does not have broad wings. For our purposes, we combine the selftuning of our fast MH algorithm implemented in GAlPAK^{3D} with the EMCEE sampler once convergence is reached as follows: we first run the selftuning of our fast MH algorithm with Cauchy sampling, and then use the EMCEE sampler once convergence is reached.
Furthermore, as illustrated in Fig. 2, we use the following twostep process. We first determine the morphological parameters by fitting a 2D Sérsic model to the collapsed [O II] cube (i.e. a twodimensional flux map), keeping the inclination i free. This yields estimates of the morphological parameters, namely inclination, size and Sérsic index (i_{2d}, R_{e,2d}, n_{2d}). Then, we fit a 3D model with the Sérsic n fixed to 0.5, 1, 2 or 4 based on n_{2d}, and tophat priors on the halflight radius R_{e} and inclination i from the results of the 2D analysis (Fig. 2). For the kinematics, we choose a tanh(r/r_{t}) rotation curve v(r) after finding that it gives more robust measurements of V_{max} than a arctan profile. Finally, from the morphological and kinematics information, namely Σ(r) and v(r) we can compute j using Eq. (2). The total sAM is then calculated from the j(< r) profile at 10R_{e}^{6}.
Fig. 2.
Schematic description of our 3D methodology using a 2step process. First, we determine the morphological parameters by fitting a 2D Sérsic model for the flux profile Σ(r) on a [O II] image from MUSE. Then, we fit the MUSE [O II] data with 3D model with the Sérsic n fixed to 0.5, 1, 2 or 4 based on n_{2d} and priors on the halflight radius R_{e} and inclination i from the results of the 2D analysis. We run the MCMC algorithm in two steps, first with Cauchy sampling which reaches convergence more efficiently, then with EMCEE which better sample the posteriors. 
For the velocity dispersion, we evaluate the intrinsic (deconvolved) velocity dispersion profile σ(r) at 2R_{e}, namely σ_{t} ≡ σ(2R_{e}). This is a more reliable estimator of the dispersion in the outer regions from our parametrization discussed earlier. Indeed, the fitted parameter, σ_{0}, can approach zero in some galaxies when the kinematic dispersions (σ_{1, 2}) is sufficient to describe the data. In the remainder of the paper, we use σ(2R_{e}) as the total velocity dispersion σ_{t} of SFGs.
In summary, our 3D method has several advantages and avoids the shortcomings mentioned in the previous section. In particular, it includes all of the information available by combining the many spaxels with very low signaltonoise ratio (S/N), down to ∼0.1 in the outer regions, which are inaccessible with 2D velocity maps. Second, it avoids the assumptions of the Eq. (3) formula, namely of flat rotation curves, which can lead to biases in the lowmass regime. Third, it takes into account the effect of the PSF convolution. And fourth, it is selfconsistent and does not rely on external data to break the i − V_{max} degeneracy. It does, however, assume that SFGs have regular velocity fields. Our modeling is similar to the DYsmal code used in Bouché et al. (2007), Cresci et al. (2009), Davies et al. (2011), Genzel et al. (2017, 2020), but GAlPAK^{3D} adjusts the model parameters directly to the raw 3D cube while DYsmal fits the 1D flux, velocity and dispersion maps generated from a 3D model.
3. Methodology validation with the TNG50 simulation
As our source of simulated galaxies we employ the TNG50 simulation (Nelson et al. 2019; Pillepich et al. 2019) from the IllustrisTNG project (Springel et al. 2018; Pillepich et al. 2018a; Nelson et al. 2018; Marinacci et al. 2018; Naiman et al. 2018). TNG50 is a cosmological hydrodynamical simulation following a volume of (51.7 Mpc)^{3} at a baryonic (DM) mass resolution of 8.5 × 10^{4} M_{⊙} (4.5 × 10^{5} M_{⊙}) and a spatial resolution of ∼ 300 pc. The simulation is evolved with the AREPO code (Springel 2010), and includes physical models for radiative cooling, starformation, stellar evolution, black hole growth, as well as feedback channels in the form of galactic winds and several modes of AGN feedback. We refer the reader to Vogelsberger et al. (2013), Weinberger et al. (2017) and Pillepich et al. (2018b) for full details on these models and only describe here the one most directly relevant to this work, namely the starformation model. Starformation in TNG50 is based on Springel & Hernquist (2003) and prescribes the starformation rate as a function of volumetric gas density with a threshold of n_{H,th} = 0.13 cm^{−3} and a gas consumption timescale of . The TNG simulations reproduce rather well several basic properties of galaxies that are relevant to the topic of this paper, such as galaxy stellar mass functions (Pillepich et al. 2018a) and sizemass relations (Genel et al. 2018) at various redshifts.
The angular momentum of galaxies was studied by Genel et al. (2015) in the precursor Illustris simulation (Vogelsberger et al. 2014a,b; Genel et al. 2014) and was found to match z = 0 observations within the uncertainties.
3.1. Sample selection
We select 100 SFGs from TNG50 at each of the z = 0.5 and z = 1.0 snapshots with stellar masses in the range 10^{8.5} M_{⊙} < M_{⋆} < 10^{10.5} M_{⊙}. We selected SFGs with specific SFR within ±2σ from the median sSFR, excluding galaxies with sSFR< 10^{−11} yr^{−1} in order to mimic the observational selection on SFGs with [O II] emission. From this sample, we randomly select the final sample to have a uniform distribution in log M_{⋆} and limited the sample to the range of sizes matching the UDF sample, namely from 1.5 to 6.5 kpc (see Sect. 4). This resulted in a sample of 158 SFGS.
Figure 3 shows the location of the sample of 158 SFGs in the size–mass and sSFR–mass scaling relations. In both panels, the gray points and contours represent the distributions for the z = 1 TNG50 SFGs and the red symbols represent the 158 selected SFGs at z = 0.5 and z = 1.0. On the left panel, the gray band represents the size–mass relation from van der Wel et al. (2014) and on the right panel, the gray band represents the mainsequence for UDF SFGs (Boogaard et al. 2018). In the left (right) panel, the z = 0.5 galaxies are rescaled to z = 1 according to the size (SFR) evolution found by van der Wel et al. (2014) (Whitaker et al. 2014), respectively. This figure shows that we selected SFGs with properties representative to the full underlying population around the size and SFR scaling relations for SFGs.
Fig. 3.
Sizes and SFRs of the 158 galaxies selected from the overall TNG50 galaxy population. Left: sizemass relation for the selected TNG50 for the mocks where the sizes are halfSFR sizes (calculated in faceon projection and without any dust attenuation). The open (solid) squares represent the z = 0.5 (z = 1) galaxies where the sizes for the z = 0.5 galaxies are adjusted to z = 1 according to the (1 + z)^{−0.5} evolution found by van der Wel et al. (2014). The background gray points and contours represent the sizemass distribution for the z = 1 TNG50 galaxy population. The solid line (gray band) represents the sizemass relation for late types from Mowla et al. (2019) (van der Wel et al. 2014), respectively. Right: SFRmass or “mainsequence” relation for the selected TNG50 galaxies (red dots, where SFRs for the z = 0.5 galaxies are adjusted to z = 1 according to the (1 + z)^{2.2} evolution found by Whitaker et al. 2014). The background gray points and contours represent the SFRmass distribution for the z = 1 TNG50 galaxy population. For comparison, we show with the gray band the z = 1 mainsequence from the MUSE UDF analysis of Boogaard et al. (2018). 
Finally, Fig. 4 shows highresolution images of three representative starforming galaxies where the two left columns show the stellar light (dustfree color composite of the (SDSS)rg(Johnson)B bands) and the two right columns show the SFR in the gas phase, faceon and edgeon.
Fig. 4.
Visual representations of three example TNG50 starforming galaxies from the set analyzed in this work. The two columns on the left represent stellar light (dustfree color composite of the (SDSS)rg(Johnson)B bands) and the two on the right the SFR in the gas phase. Each panel is 30 kpc on a side. 
3.2. Simulated MUSE cubes
Each of the 158 simulated SFGs is processed in the following manner to generate several mock MUSE minicubes, for a total number of mock cubes of 790. First, a cube of two spatial dimensions with bounds [ − N_{r}/2Δ_{r}, N_{r}/2Δ_{r}] and one velocity dimension with bounds [ − N_{v}/2Δ_{v}, N_{v}/2Δ_{v}] is generated with N_{r} = 30 pixels in each spatial dimension and N_{v} = 40 pixels in the velocity dimension. For z = 1 (z = 0.5), we use voxels of size Δ_{r} = 1.65 kpc (1.25 kpc) and Δ_{v} = 50 km s^{−1} (65 km s^{−1}), respectively, which corresponds to the MUSE instrument specifications of 0.2″ and 1.25 Å at a restframe wavelength of λ_{rest} = 3728 Å (Bacon et al. 2017). Second, the center of this cube is placed at the position of the most bound resolution element of the simulated galaxy, which generally is very close to the center of rotation, and at the velocity of the galaxy centerofmass. Third, the cube is rotated such that the plane defined by the two spatial dimensions is perpendicular to the eigenvector of the moment of inertia tensor of the galaxy that has the smallest eigenvalue, such that the galaxy appears “faceon”.
For each galaxy, we apply further rotations to generate mock cubes for 5 inclinations (i = 15, 30, 45, 60 and 75°), each at a random azimuthal orientation ϕ^{7}. The SFR in each voxel of the cube is then calculated by summing up the contributions from all simulation gas cells using a Gaussian kernel representing the PointSpreadFunction (PSF) and LineSpreadFunction (LSF) appropriate for nonAO MUSE observations. In particular, for z = 1 (z = 0.5) the PSF has a Moffat profile with β = 2.5 and a FWHM_{r} of 0.625 arcsec (0.675 arcsec) (Bacon et al. 2017), and the LSF has a Gaussian profile with FWHM_{v} = 2.55 Å (2.75 Å), respectively. The PSF and LSF FWHMs correspond to 5.15 kpc (4.25 kpc) and to 100 km s^{−1} (143 km s^{−1}), for z = 1 (z = 0.5), respectively.
Finally, we add noise to the mock MUSE cubes in order to match the range of S/N (spaxel^{−1}) of our MUSE observations, which have S/N_{max} of a few to ∼40 (see Sect. 4). For galaxies with log M_{⋆}/M_{⊙}>9, we set the noise to 0.001 (0.0002) SFR yr^{−1}, respectively.
3.3. Reference measurements
In TNG50, we measure the following quantities: the sAM, the maximum rotation velocity V_{max} and velocity dispersion σ_{t}, as follows.
Specific angularmomentum. The SFRweighted is calculated using Eq. (1) with ρ(x) as the starforming gas density, which corresponds to a [O II] weighting, using the 3D positions, velocities and densities of all SF gas at a given radius. We then use the cumulative profile to then define the total sAM from max.
Velocity dispersion. The gas velocity dispersion σ_{t} is defined as the 1D lineofsight dispersion assuming isotropy, namely it is equal to . Here, σ_{3D} is the average (SFRweighted) of the internal velocity dispersion of starforming gas inside kpcsized patches at galactocentric radii of 1.5–2.5R_{e}, a method similar to Pillepich et al. (2019). In order to account for thermal broadening for T = 10^{4} K gas, we incorporate a σ_{th} ≈ 15.7 km s^{−1} component into σ_{3D} by adding it in quadrature^{8}. This definition of lineofsight SFRweighted velocity dispersion at ∼2R_{e} has the advantage that it closely mimics our observations σ_{t}(2R_{e}) (see Sect. 2.3).
Rotational velocity. The rotation velocity V is calculated as the SFRweighted mean of the azimuthal velocities of the same patches used for the calculation of the velocity dispersion.
3.4. Robustness of (3D) from TNG50
We now apply the methodology described in Sect. 2.3 to the sample of 790 mock cubes of 158 simulated galaxies in order to test the performances of GAlPAK^{3D} in recovering the morphological (sizes, inclination) and kinematic parameters of realistic galaxies. We find that the inclination is recovered with a relative error of ±12° provided that i ≥ 30°, as found in Contini et al. (2016) (see also Bouché et al. 2015). Consequently, we restrict the analysis to mock cubes with i ≥ 30°.
Figure 5 shows the sAM profile found when applying Eq. (2) and the 3D methodology described in Sect. 2.3 on the massive galaxy shown at the top of Fig. 4. The 5 panels show the Σ_{sfr} map, the SB profile, the intrinsic velocity field, the intrinsic velocity profile, and the intrinsic profile, respectively. The velocity and SB profile are computed on the pseudoslit, represented by the two horizontal gray lines. The light gray solid lines represent the profiles recovered from our algorithm for each inclination (i ≥ 30°), each performed on 5 orientations in this case. The red solid lines in columns 2, 3 and 4 represent the Σ_{sfr}, v(r) and profiles determined directly from the TNG50 data. In the right most panel, the blue horizontal tick marks show our estimates , while the red horizontal tick marks show from RF12 (Eq. (3)). Depending on the steepness of the v(r) profile, the RF12 approximation can lead to an overestimation of the total angular momentum by 15–25% for the reasons described in Sect. 2. Figure 5 suggests that our 3D method is able to recover the intrinsic profile with less bias, even when the rotation curve deviates from our assumed parametrization v(r) = tanh. Next, we investigate the difference between the estimated and true as a function of S/N.
Fig. 5.
Example of GAlPAK^{3D} fits for three SFGs from TNG50 shown in Fig. 4. The columns show the flux map, the flux profile Σ(r), the intrinsic velocity field, the deprojected velocity profile and the deprojected angular momentum profile, from left to right respectively. The 1D profiles are computed along the pseudoslit shown in the first and third columns by the gray lines. In the second column, the SB profile for a single orientation is shown (black points). In the second, fourth and fifth columns, the solid gray curves show the model fits for the various inclinations and orientations. The red solid lines represent the true profiles determined directly from the TNG50 data. The vertical dotted lines represent twice the halflight radius (R_{1/2,SFR}). The horizontal blue lines show the total modeled angular momentum j_{3D}, while the horizontal red lines show the total j from RF12 (Eq. (3)), which can lead to overestimation of the total angular momentum by 10–20%. 
Figure 6(left) shows the measured total angular momentum, (3D), as a function of the true angular momentum (SFRweighted) as measured directly from the simulation output (bottom subpanel), and their relative difference, (top subpanel). The points are colored according to M_{⋆}. The red circles represent the median of dj/j whose errorbars represent the 1σ standard deviation.
Fig. 6.
Reliability of our sAM estimator. Left: bottom panel: total SFRweighted sAM, (in km s^{−1} kpc), as a function of the true angular momentum (SFRweighted) for the 790 mock galaxies. Top panel: relative error ( as a function of the effective S/N (see text) along with the histogram of . The points are colored according to the stellar mass in both panels. The red circles represent the median whose errorbar represent 1σ. Right: relative error ( for each mass bin of 1 dex centered on log M_{⋆} = 8.5, 9.5, 10.5, from top to bottom. The histograms are weighted by the errors. Left (right) panels: distributions of , in three mass bins, resulting from the 3D method (the Eq. (3) approximation), respectively. The red solid line represents the continuous density estimation of the distributions using a Gaussian kernel with bandwidth of ∼0.2. This shows that, for massive galaxies (M_{⋆} > 10^{10} M_{⊙}), the Eq. (3) approximation is less biased, while for galaxies with M_{⋆} < 10^{10} M_{⊙}, the 3D method is less biased. 
The histogram in the top panel of Fig. 6(left) shows the distribution of the relative errors dj/j. This histogram shows that globally our method does recover the intrinsic angular momentum with no apparent bias. We can quantify this bias further as a function of galaxy mass in Fig. 6(right). Figure 6(right) shows, for 1dex bins in log M_{⋆} (at 8.5, 9.5 and 10.5), from top to bottom), the distributions of dj/j resulting from our 3D method (left column) or the RF12 method (right column). The red solid line represents the continuous density estimation of the distributions using a Gaussian kernel with bandwidth of 0.2. Quantitatively, the mean relative error dj/j for our 3D methodology is < 5%, while for the intermediate and lower mass bins, the RF12 approximation (Eq. (3)) can lead to a significant bias of ≈10–15% for the reasons discussed in Sect. 2. For the most massive galaxies with M_{⋆} > 10^{10} M_{⊙}, and with , our method underestimates compared to RF12, but it is important to note that our sample in the UDF (Sect. 4) has only a handful of SFGs in this regime.
The width of the distributions shown in Fig. 6 is consistent with the errors obtained with our methodology. Indeed, the width of the ratios of the errors to dj is consistent with a (Cauchy) distribution^{9} of width unity. This implies that the errors obtained from our methodology are broadly representative of the true error dj on the measurement.
We return now to the top panel of Fig. 6 where the relative difference dj/j is plotted as a function of the effective S/N, S/N_{eff}, following the discussion in Bouché et al. (2015). This effective S/N accounts for the fact that surface brightness (SB) alone is not sufficient to determine the accuracy in the fitted parameters, given that the compactness of the galaxy with respect to the PSF plays also an important role. As discussed in Refregier et al. (2012), in the presence of a PSF convolution, the relative errors on the majoraxis a scales as
where A_{o} is the observed central surface brightness in the central voxel, R_{PSF} is the radius of the PSF (R_{PSF} ≡ FWHM/2) and R_{e} the intrinsic halflight radius. After performing a Taylor expansion around R_{PSF}/R_{e} ∼ (1 − x) with x ≡ (R_{e} − R_{PSF})/R_{e} and x≪1, one finds that Eq. (5) becomes in the regime where R_{PSF}/R_{e} is ≃1.0 (see Bouché et al. 2015):
The last step follows from the simple argument that a given exposure time sets the surface brightness limit per voxel, and hence the maximum S/N (in the central voxel or spaxel), S/N_{max}, is directly proportional to A_{o}.
Considering the sAM is essentially determined by the product of R_{e} and V_{max}, one can thus expect that the accuracy of j follows Eq. (6), which is represented by the gray curve in the top panel of Fig. 6.
Finally, we end this section by quantifying the potential difference between the SFRweighted and the stellar massweighted j_{⋆} because in observations of highredshift galaxies, as in Sect. 4, one is often limited to the SFRweighted . Using TNG50 z = 1 SFGs in the parent sample, Fig. 7 shows the relation between and j_{⋆}. While this figure shows that the relation becomes non linear in the low mass regime, with log j_{⋆} < 1.5, (as found by ElBadry et al. 2018, in the FIRE simulations), it also shows that, on average, is an estimate of log j_{⋆} + 0.25 dex over a wide range of sAM from 50 to 2000 km s^{−1} kpc, which is the relevant range for our study.
Fig. 7.
Relation between the stellar sAM j_{⋆} and for z = 1 TNG50 SFGs. The solid points with errorbars represent the median difference between the two as a function of the stellar angular momentum, where the error bars represent the standard deviation. Thus, the SFRweighted sAM, , is higher by ≈ 0.25 dex (dotted line) compared to j_{⋆} over a wide range of sAM from 50 to 2000 km s^{−1} kpc. A second order polynomial fit is shown with the dotdashed black line. 
4. Application to UDF galaxies
4.1. Sample
We exploit the 3′ × 3′ observations taken with the MUSE instrument over the Hubble Ultra Deep Field (HUDF; Bacon et al. 2017) to investigate the angular momentum of intermediate redshift galaxies using the methodology presented in Sect. 2.3. The observations were performed on a mosaic of nine 1′ × 1′ 10 h pointings and contain a single 1′ × 1′ 27 h pointing, but we concentrate on the mosaic sample to ensure a homogenous selection. In the MUSE UDF mosaic, there are 494 [O II] emitters (Inami et al. 2017). These [O II] emitters have stellar masses ranging from 10^{7.5} to 10^{11} M_{⊙}, where the stellar masses are estimated (as in Boogaard et al. 2018) using the stellar population synthesis (SPS) code FAST (Kriek et al. 2009) assuming a Chabrier (2003) initial mass function with an exponentially declining star formation history. As discussed in the next section, most [O II] emitters are rather faint, with only about ∼50% of these having a S/N above ∼3.
Figure 8 (left) shows the S/N of the total [O II] flux as a function of the “resolved size” (R_{e} divided by the PSF FWHM) for the sample of SFGs, where the blue (red) points are SFGs with a S/N_{eff} > 3.0 (< 3.0), respectively. For comparison, the selection used in AbrilMelgarejo et al. (2021) is shown with the gray bands. The right panel compares the effective S/N and the total S/N, where the points are colorcoded as a function of size R_{e}, displaying the common surfacebrightness biases. This figure shows that a S/N_{eff} 3, as determined in the previous section, opens the possibility to study the kinematics of faint galaxies, enlarging the sample by 30–50%.
Fig. 8.
Left: following AbrilMelgarejo et al. (2021), we show the S/N of the total [O II] flux as a function of the resolved ratio (R_{e, hst}/PSF FWHM) for the sample of SFGs. The gray bands represent the selection of resolved galaxies in AbrilMelgarejo et al. (2021), consisting of S/N > 30 and R_{e, hst}/FWHM > 0.5. Right: comparison between the total [O II] S/N and the effective S/N (Eq. (6)), colorcoded with R_{e}. In both panels, the red points are galaxies with S/N_{eff} < 3.0, while the blue points correspond to galasxes with S/N_{eff} > 3.0. 
Figure 9 shows the properties of the UDF sample in terms of sizes and SFRs as in Fig. 3. The sizes are the continuumbased R_{e} obtained from HST/F160W from van der Wel et al. (2014) and the SFRs are derived from SED fitting. The left panel shows the sample in the sizemass plane color coded as a function of the effective S/N, where the sizes were adjusted to z = 1. For comparison, we show the sizemass scaling relation of van der Wel et al. (2014) and Mowla et al. (2019) down to M_{⋆} ∼ 10^{8} M_{⊙}. The TNG50 mocks used in Sect. 3 are represented with the red small circles. The right panel shows the mainsequence or the SFRmass plane where the points are colorcode as a function of size. The points are scaled to z = 1 where we used the mainsequence redshift evolution of Boogaard et al. (2018).
Fig. 9.
Comparison of selection techniques. Left: sizemass relation for the SFGs in the UDF mosaic sample (blue circles). The sizes are adjusted to z = 1 according to the (1 + z)^{−0.5} evolution found by van der Wel et al. (2014). The solid line (gray band) represent the sizemass relation of Mowla et al. (2019) (van der Wel et al. 2014). Right: mainsequence relation for the SFGs in the UDF mosaic sample. The SFRs are adjusted to z = 1 according to the (1 + z)^{1.74} evolution of Boogaard et al. (2018) whose z = 1 main sequence is represented by the gray band. In both panels, the points are colorcoded according to effective S/N. For comparison we show the properties of the mock 158 galaxies (Sect. 3.2) as red circles. 
4.2. Size estimates
The first part of the methodology outlined in Sect. 2.3 is to determine the morphological parameters (size, inclination) by performing two dimensional Sérsic fits to the [O II] flux maps. Given that the UDF data set has extensive Hubble Space Telescope (HST) coverage, we first ask whether our Sérsic fits on [O II] (from groundbased MUSE observations at a resolution of 0.62″) are able to recover the HST morphological parameters (at a resolution of 0.15″) from their stellar continuum.
Given that the galaxy halflight radius plays a crucial role in determining j, we show in Fig. 10(left) the intrinsic [O II]based R_{e} found from our MUSE data compared to the continuumbased R_{e} obtained from HST/F160W from van der Wel et al. (2014). This figure shows that there is a good agreement between the two quantities, provided that the sizecorrected “effective S/N” (Eq. (6)) in the brightest voxel, namely S/N_{max} × (R_{e}/R_{psf}) where S/N_{max} is the S/N in the central voxel and R_{psf} is the PSF radius (halfFWHM), is ⪆3. The S/N level of ∼3 is our fiducial threshold for the remainder of our analysis. From the 494 SFGs in the 9 arcmin^{2} UDF mosaic, there are 182 galaxies meeting this criterion.
Fig. 10.
Left: comparison between [O II]based sizes R_{e} from MUSE (at FWHM ≈ 0.7″ resolution) versus Hband continuum R_{e, ⋆} from HST/F160W at ≈0.15″ resolution, for the full sample of 494 SFGs, where the points are colorcoded according to Eq. (6). The top panel shows the relative error δR_{e}/R_{e} as a function of the effective S/N (see text) and the side histogram shows its distribution. The gray errorbars represent the 3σ_{t} statistical uncertainties. This shows that our MUSE[O II] sizes are of sufficient accuracy down to ≈1 kpc, provided that the effective S/N is ⪆3. Right: sizes R_{e} (from HST/F160W) as a function of redshift for 182 galaxies in the MUSE UDFmosaic field whose effective S/N is sufficient. The gray region represents the region probed by Swinbank et al. (2017). The dashed line represents ∝(1 + z)^{1}. The vertical error bars are 3σ. 
This demonstrates that our algorithm is able to determine R_{e} with sufficient accuracy down to ≈1 kpc. This is beyond the 3 kpc limit of other IFS surveys as the parameter space between 1 < R_{e} < 3 kpc (and 10^{8} < M_{⋆} < 10^{9.5} M_{⊙}) has been relatively unexplored so far with any of the large IFS surveys (e.g., Burkert et al. 2016; Swinbank et al. 2017; Harrison et al. 2017). Figure 10 shows the location of halflight radii R_{e} of our subsample of 182 galaxies as a function of redshift where the parameter space probed by Swinbank et al. (2017) is shown as the gray band. This figure shows that, compared to Swinbank et al. (2017), our sample includes galaxies with sizes in the range 1–3 kpc, previously inaccessible from groundbased observations.
4.3. The angular momentum of SFGs at 0.4 < z < 1.4
We apply the twostep methodology outlined in Sect. 2.3 on the MUSE UDF sample of 494 [O II] emitters. We find that 182 meet our S/N threshold discussed above, namely with effective S/N ≳ 3. The sample is thus made of all [O II] emitters in the UDF MUSE mosaic with sufficient S/N and is composed of SFGs with Sérsic profiles that are mostly exponential. Two thirds of the sample has n = 1.0, while the remainder has n = 0.5. The distribution of inclinations i is consistent with an unbiased sample, namely P(cos(i)) ∼ b/a is flat for i > 30 (b/a < 0.86). As discussed in Sect. 2.3, we compute the angular momentum of the starforming gas by integrating Eq. (2) from the flux Σ(r) and kinematic v(r) GAlPAK^{3D} fits.
Figure 11(left) shows the resulting sequence for the 182 galaxies in our sample. This figure shows the sAM (SFR weighted) as a function of stellar mass M_{⋆} where the points are colorcoded according to log(V_{max}/σ_{t}) where σ_{t} is the outer velocity dispersion (see Sect. 2.3) measured with GAlPAK^{3D} (defined as in Cresci et al. 2009). The right axis shows the corresponding stellar massweighted j_{⋆} from the typical offset found in Fig. 7.
Fig. 11.
Specific angular momentum for the 182 galaxies in the MUSE UDF mosaic fields with a max S/N ≳ 3. Left: total specific angular momentum (SFR weighted) as a function of stellar mass M_{⋆}. The points are colorcoded as a function of V/σ_{t}. The rightaxis shows the corresponding j_{⋆} according to the 0.25dex offset found in Fig. 7. The blue (red) dashed line show the sequence for late (early) type galaxies at z = 0 from Fall & Romanowsky (2018, FR18). The dotted line (solid squares) show the z = 0 sAM j_{⋆} − M_{⋆} sequence from Posti et al. (2018a) (Lapi et al. 2018), respectively. The dotted line show the sequence for z = 0 disk galaxies of Posti et al. (2018a). The solid black line shows the expected j_{⋆} − M_{⋆} relation from Eq. (8). The solid blue and red lines show the fitted fiducial model from Eq. (9) and the gray band represents the 95% confidence interval. Right: normalized AM sequence as a function of (V/σ_{t})_{toz1} (normalized to z = 1 see text), i.e. after taking into account for the redshift evolution of σ_{t}(z), where the points are colorcoded with redshift. In the inset, we show the redshift evolution of the sAM relation at M_{⋆} = 10^{10} M_{⊙} (solid line) at fixed (V/σ_{t})_{toz1} (=5) along with the 95% predictive interval (gray band), with the z = 0 ZP from Fall & Romanowsky (2018) (red square), the Swinbank et al. (2017) binned data (see text), and the expected (1 + z)^{−0.5} evolution (dotted line) from Eq. (8). 
In Fig. 11(left), the dotteddashed blue (red) line show the j_{⋆} − M_{⋆} sequence for late (early) type galaxies at z = 0 from Fall & Romanowsky (2018), respectively. The dotted line shows the z = 0 sequence for the sample of disk galaxies of Posti et al. (2018a) which has a pure powerlaw sAM (see also Mancera Piña et al. 2021). The thin black dotted curve represents the selection limit of our observations of described in Appendix B.
In a cosmological framework, the j_{⋆} of a galaxy with stellar mass M_{⋆} that resides in a virialized halo of mass M_{vir} at a given redshift z can be written as (following Mo et al. 1998; Navarro & Steinmetz 2000; Romanowsky & Fall 2012; Obreschkow et al. 2015; Burkert et al. 2010, 2016; Stevens et al. 2016; Swinbank et al. 2017):
where λ is the halo spin parameter, f_{m} ≡ M_{⋆}/M_{h} is the stellartohalo mass ratio, f_{j} ≡ j_{⋆}/j_{h} is the j retention factor (or the stellartohalo sAM ratio), and the expected (1 + z)^{−1/2} redshift evolution follows from H(z), the Hubble parameter and Δ_{c}(z), the halo overdensity relative to the mean density of the universe (Bryan & Norman 1998). For constant f_{m} in massive galaxies, the j_{⋆} − M_{⋆} relation is expected to be , while in lowmass galaxies, j_{⋆} should be given that abundance matching shows that (Moster et al. 2010; Dutton et al. 2010; Behroozi et al. 2013) . In other words, the j_{⋆} − M_{⋆} relation should have a slope of 2/3 (1/3) in the high(low) mass regime, respectively, provided that the retention fraction f_{j} is constant (see Posti et al. 2019, for detailed discussions).
Following Romanowsky & Fall (2012), Eq. (7) becomes
where M_{⋆,10} ≡ M_{⋆}/10^{10} M_{⊙}, λ the halo spin parameter and (1 + z)_{1} = (1 + z)/2. Eq. (8) is shown as the solid black line in Fig. 11(left), assuming f_{m} from abundance matching (Dutton et al. 2010), λ = 0.034 from DM simulations (Macciò et al. 2007; Bett et al. 2007, 2010), and f_{j} = 1.
4.3.1. The j − M − (V/σ) relation
In the local universe, a transition between the late to early type j − M sequence has also been observed (Romanowsky & Fall 2012; Fall & Romanowsky 2013, 2018; Obreschkow & Glazebrook 2014; Murugeshan et al. 2020), where the j − M relation is a function of the bulgetototal ratio B/T. Similarly, Harrison et al. (2017) found, on a sample of 586 SFGs at z ≃ 0.9, that the j − M relation is a function of Sérsic index n, where sAM decreases with increasing n. Cortese et al. (2016), using SAMI, showed that the j − M relation is a function of the stellar spin parameter, λ_{R} (as defined in Emsellem et al. 2007), thus defining a continuous sequence in the sAMstellar mass plane (j − M − λ_{R}).
At intermediate redshifts, Contini et al. (2016), with only 28 0.4 < z < 1.4 galaxies in the MUSE HDFS, found that z ≃ 1 SFGs fill a continuous transition between rotation dominated (with V/σ_{t} > 1) to dispersion dominated galaxies (with V/σ_{t} < 1) where the two populations followed the late/early type sAM sequence of Fall & Romanowsky (2013). This V/σ trend was also present in the high mass sample of Burkert et al. (2016). Here, with our larger sample of intermediate redshift SFGs from the MUSE UDF mosaic 9 arcmin^{2} fields, we revisit this question.
To quantify this j − M − (V/σ) relation, we show in Fig. 11(right) the normalized AM sequence as a function of V/σ_{t} where the points are colorcoded with the redshift z. The normalization f(M) removes the mass dependence obtained from the fit (see below). Figure 11 clearly shows that the z ∼ 1 AM sequence is actually a 3dimensional j − M − (V/σ_{t}) sequence, meaning that the sAM is a strong function of the galaxy dynamical state V/σ_{t}. This confirms the early results of Contini et al. (2016) and extends the results of Burkert et al. (2016) (their Fig. 2) for highmass SFGs with M_{⋆} > 10^{10} M_{⊙} to low mass galaxies with M_{⋆} ∼ 10^{8} M_{⊙}.
Following Obreschkow & Glazebrook (2014) and Fall & Romanowsky (2018) who fit j_{⋆}(M_{⋆}, B/T) as a 2dimensional plane in a 3D space, we fit the AM sequence in Fig. 11 with a 3dimensional function in M_{⋆}, V/σ_{t} and redshift z. For this purpose, we use a quadratic polynomial f(log M_{⋆}) in log M_{⋆}, a linear redshift dependence in log(1 + z) and a linear relation with log(V/σ_{t}) such that:
where log M_{⋆,10} = log(M_{⋆}/10^{10} M_{⊙}), (V/σ_{t})_{5} is V/σ_{t}/5, and (1 + z)_{1} = (1 + z)/2. In this parametrization, a captures the redshift evolution (1 + z)^{a} and b is the zeropoint of the sAM relation at z = 1 for rotation dominated SFGs with V/σ_{t} = 5 and M_{⋆} = 10^{10} M_{⊙}.
In fitting for the redshift evolution in Eq. (9), a key assumption is often made (e.g., Swinbank et al. 2017), namely that the distribution of V/σ_{t} is redshiftindependent. However, V/σ_{t} is known to be strongly evolving with redshift (e.g., Förster Schreiber et al. 2006; Genzel et al. 2008; Epinat et al. 2012; Kassin et al. 2012; Wisnioski et al. 2015; Übler et al. 2019). This has two (but related) important consequences. First, any particular definition of disk galaxies (e.g., with a constant cut at V/σ_{t} > 1 or > 2) becomes a redshiftdependent selection of the underlying popution. Indeed, a population with a V/σ_{t}(z) distribution evolving with redshift means that a fixed V/σ_{t} cut (commonly used in the litterature) will select a changing fraction of the underlying disk population. Second, this evolution should be taken into account if the redshift evolution (1 + z)^{a} captured in Eq. (9) is to mean the redshift evolution at a fixed dynamical state or for the median of the popultation.
Since we aim to compare the sAM sequence across our wide redshift range 0.4–1.4, the redshift evolution of V/σ_{t} can bias the j − M redshift evolution, and cannot be neglected. Thus, the k term in Eq. (9) ought to be k = γ log[V/σ_{t})_{toz1}/5], where (V/σ_{t})_{toz1} is (V/σ_{t}) corrected to z = 1. This ensures that we are describing the evolution of the sequence without being impacted from the redshift evolution of V/σ_{t}.
To estimate (V/σ_{t})_{toz1}, we use the redshift evolution of V/σ_{t}, discussed extensively in Übler et al. (2019) and Pillepich et al. (2019), which comes primarily from the evolution of σ_{t}(z) given the negligible evolution of the TullyFisher relation (e.g., Di Teodoro et al. 2016; Übler et al. 2017; Tiley et al. 2019). It can be parametrized as
Ignoring this effect will increase bias the redshift evolution of the j − M relation, specifically the a coefficient in Eq. (9) will be biased by approximately q × γ.
We find that the parameter q in Eq. (10) can be found by minimizing the V/σ_{t} distribution over the full redshift range (as discussed in Appendix D). Figure D.1(right) shows that q ≃ −0.5 is most appropriate for the UDF sample where σ_{t} ≡ σ(2R_{e})^{10}.
4.3.2. The j − M − (V/σ) global fit
We use a Bayesian algorithm to fit the 5 parameters (α, β, γ, a, b) using a Studentt likelihood^{11} with the NoUTurn Sampler (NUTS) from the PYMC3 python package (Salvatier et al. 2016). In Fig. 11, we show the resulting fit to the j − M − (V/σ_{t}) relation using Eq. (9) whose fitted parameters are listed in Table 1. From this fit, one sees that the relation is a nonlinear function of M_{⋆}, as the relation flattens toward low masses given that the quadratic factor α is significantly non zero (α = 0.04 ± 0.01). Specifically, the slope of the sAM relation goes from s ≈ 0.5 at log M_{⋆}/M_{⊙} = 10.5, consistent with numerous results at z = 0 (e.g., Obreschkow & Glazebrook 2014; Cortese et al. 2016; Posti et al. 2018a; Fall & Romanowsky 2018) or highredshifts (e.g., Swinbank et al. 2017; Marasco et al. 2019), to s ≈ 0.25 in the low mass regime at log M_{⋆} = 8. In the remainder of this paper, we use f(log M) as short hand for this nonlinear mass dependence in Eq. (9) for the remainder of this work.
This nonlinear j − M relation is expected from Eq. (7) where the stellar mass fraction f_{m} plays a significant role in this mass regime. Together, this result indicates the angular momentum retention fraction f_{j} at z ∼ 1 is not far from being constant with mass (see Posti et al. 2018b, 2019, for discussions on this at z = 0).
4.3.3. The j − M redshift evolution
In regards to the redshift dependence of the sAM sequence, the inset in Fig. 11(right) shows the redshift evolution of the zeropoint ZP_{10} (Eq. (9)) of the relation fitted over our redshift range 0.4 < z < 1.4 for M_{⋆} = 10^{10} M_{⊙} and V/σ_{t} = 5. The red circle shows the RF12 z = 0 zeropoint, and the squares show the Swinbank et al. (2017) results^{12}. The dashed line shows the (1 + z)^{−1/2} expected evolution (Eq. (7)). This panel shows that the sAM relation (at fixed V/σ_{t}) is evolving as
with (95%) which is consistent with the theoretical evolution (1 + z)^{−1/2} from Eq. (7). This redshift dependence is not consistent with the slope of −1 found by Swinbank et al. (2017) to z = 1.5, but it should be noted that they did not perform a global fit to their data and that their binned data (squares in the inset of Fig. 11(right) are in fact consistent with a slope of −0.5 up to z = 1.0, with the z = 1.5 binned data not matching the (1 + z)^{−0.5} trend is most likely due to incompleteness issues at these redshifts.
While Swinbank et al. (2017) and Harrison et al. (2017) found some evidence for a redshift evolution of the j − M relation with a combined dataset of ∼1000 SFGs, in contrast, Marasco et al. (2019) found no evidence for a redshift evolution of the relation on a subset (made of 17 z = 1 disk galaxies with M_{⋆} > 10^{9.5} M_{⊙}) of the KROSS and KMOS3D survey. This discrepancy implies no evolution of both the TFR relation and of the sizemass relation, in apparent contradiction to several studies (e.g., Straatman et al. 2017; Übler et al. 2017; van der Wel et al. 2014). This discrepancy could originate from the small sample size, from the different methodologies involved, and/or could be resolved under our (V/σ_{t}) framework.
5. Comparisons to SFGs in TNG50
In order to compare our results on the sAM relation presented in Sect. 4.3 to TNG50, we selected 18000 SFGs at z = 0.5, 1.0 and 1.5 from TNG50. We matched the properties of TNG50 SFGs to the UDF SFGs as follows. We first match the SFR distribution by imposing a log(SFR/[M_{⊙} yr^{−1}]) > −1.5 as well as impose a similar observational limit of , and then match both the redshift and stellar mass distributions of the UDF sample. The resulting sample is made of 5180 SFGs.
Figure 12(left) shows the SFRweighted sAM as a function of M_{⋆} for SFGs in TNG50 with M_{⋆} > 10^{8} M_{⊙}, colorcoded with (V/σ_{t}). Fig. 12(right) shows the relation between the normalized sAM and (V/σ_{t})_{toz1}, that is V/σ_{t} corrected for its redshift evolution to z = 1 using the scaling V/σ_{t} ∝ (1 + z)^{0.87} found in Fig. D.1. In this figure, we show only 1000 SFGs for clarity. Figure 12 shows that SFGs in TNG50 follow similar trends with M_{⋆} and with V/σ_{t} as in the UDF sample. In particular, the relation is also a strong function of V/σ_{t}.
To quantify potential differences, we applied the same parametric fit as for the UDF sample on a randomly selected subset of 250 galaxies in order to have statistical errors similar to the UDF sample, namely using Eq. (9) with the redshift evolution of V/σ_{t} (Eq. (10)). The results are listed in Table 1 with and without the correction for the redshift evolution of V/σ_{t}. One sees that the V/σ_{t} dependence is stronger in TNG50 than in the UDF sample with γ ∼ 1.3, and the redshift evolution for the TNG50 SFGs is consistent with the (1 + z)^{−0.5} expectation, with a = −0.45 ± 0.25, against for the UDF SFGs.
We end this section by noting that we found no strong dependence of the sAM with age (formation redshift) or with morphological indicators (Tacchella et al. 2019), such as GiniM20, and bulgetototal ratio B/T.
6. The j − M − (V/σ) relation: implications
As shown in Figs. 11 and 12, the sAM sequence of SFGs in the UDF and in TNG50 represents a strong function of the galaxy dynamical state as characterized by V/σ_{t}, where rotation dominated SFGs with V/σ_{t} ∼ 5 and dispersion dominated SFGs with V/σ_{t} = 0.1 are offset by ∼1dex. This result reinforces the preliminary results of Contini et al. (2016) on only 27 intermediate redshift SFGs, and extends the similar trend present in Burkert et al. (2016) for massive (with M_{⋆} > 10^{10} M_{⊙}) SFGs.
This j − M − (V/σ_{t}) relation has important consequences. If SFGs form with high sAM, on the rotationdominated j − M sequence, this j − M − (V/σ_{t}) relation indicates that SFGs may experience a dynamical transformation to lower their sAMs prior to becoming passive, and remain on the earlytype j − M sequence. This might also lead to a morphological transformation as the j − M residuals also correlate with Sérsic n (Harrison et al. 2017), but this trend is not found in the UDF data. Indeed, the vast majority of our SFGs have Sérsic indices n ≲ 1, thus they have not undergone the morphological transformation related to bulge formation. This indicates that the kinematical transformation (from high V/σ_{t} to low V/σ_{t}) occurs before the morphological transformation toward bulge formation.
This scenario of kinematical transformation toward low j with smaller sizes ^{13} and larger dispersions is in agreement with the recent finding from the CANDELS survey (Osborne et al. 2020) that SFGs contains ∼15% of compact objects (‘blue nuggets’), indicating that compaction precedes quenching.
In local galaxies, Oh et al. (2020) separated the bulge and disk (stellar) kinematics of galaxies in the SAMI sample and found that the bulge (disk) component have low (high) V/σ, and follows the FaberJackson (TullyFisher) relation, respectively. We hypothesize that our 0.4 < z < 1.4 SFGs are their progenitors, and z = 1 SFGs with low V/σ will evolve into bulge dominated systems at z = 0. This is consistent with the TNG results of Tacchella et al. (2019), who found that the morphology of galaxies is set during their starforming phase.
We return to Eq. (8) to gain insight on our results and in particular to address whether our relatively broad sAM sequence is consistent with the expected distribution of halo spins λ. One might expect the distribution of halo spins to be reflected in the scatter of the sAM sequence if each SFG retains a similar amount of sAM (a constant f_{j} = 1), as in semianalytical models. However, in hydrodynamical simulations the situation is more complex (Genel et al. 2015; Jiang et al. 2019) given the potential large AM exchanges.
Figure 13 shows the normalization of Eq. (8), namely , as a function of V/σ (corrected again to z = 1) for the MUSE sample (left) and for the TNG50 sample (right), respectively. The left yaxis is found from the “fiducial” fits shown in Figs. 11 and 12, which are listed in Table 1. The histograms in Fig. 13 show the distributions of Eq. (8) recast as λf_{j} (right yaxis). For comparison, the solid line shows the lognormal distribution P(λ) expected for darkmatter halos (Bett et al. 2007) assuming a retention factor f_{j} of unity.
Fig. 13.
Left: sAM sequence normalization as a function of (V/σ_{t})_{toz1} for the UDF mosaic sample of SFGs where the color representes the SFRweighted log R_{1/2,[O II]}. The left yaxis is found from the fiducial fits shown in Figs. 11–12 (see Table 1). The right yaxis shows the corresponding effective spin λ × f_{j} according to Eq. (8) using the zeropoing of Eq. (8). The solid histograms show the distributions from the data. The solid line shows the lognormal distribution P(λ) expected for halos (Bett et al. 2007). Right: same for the matched TNG50 sample of SFGs. 
This figure shows that the distribution of SFGs in sAM is consistent with the expectation from darkmatter halos in TNG50, but is broader in our observations. The trend with V/σ is a strong function of galaxy sizes (SFRweighted) in TNG50, but not in the MUSE sample. Another difference is the paucity of dispersion dominated galaxies in TNG50 (with 15% having V/σ < 3) compared to the MUSE sample (with 50% having V/σ < 3). Interestingly, the MUSE sample restricted to S/N_{eff} > 5 appears to show a more bimodal distributions in λf_{j}, consistent with the z = 0 RF12 results where rotation dominated galaxies have retained most of their halo sAM, while dispersion dominated systems have lost around 80% of their sAM.
7. Conclusions
The main purpose of this paper was to develop a robust methodology to measure j_{tot} from the angular momentum profiles j(< r) of galaxies in the lowmass regime (down to M_{⋆} = 10^{8} M_{⊙}), even when one is severely limited by S/N, surfacebrightness limits and/or by the spatial resolution of groundbased instruments. The methodology consists of a twostep process (Sect. 2.3) and is based on the fitting algorithm GAlPAK^{3D} (Bouché et al. 2015). The algorithm adjusts the morphological and dynamical parameters of a 3dimensional disk model against the IFU data directly.
Using mock observations of SFGs taken from the TNG50 simulations (Nelson et al. 2019; Pillepich et al. 2019) with M_{⋆} = 10^{8.5−10.5} M_{⊙}, we have shown that our 3D methodology is able to recover the angular momentum profiles of SFGs with little bias (< 5%) (Fig. 6) provided the S/N and/or surface brightness is sufficient (see Sect. 3.4). In particular, we find that we are able to estimate the j profiles of lowmass SFGs down to 10^{8} M_{⊙} with correspondingly small sizes. Specifically, galaxies must have S/N_{eff} ≳ 3, where S/N_{eff} = S/N_{max} (R_{e}/R_{psf}) and S/N_{max} is the S/N in the brightest pixel, R_{e} is the halflight radius and R_{psf} is the PSF radius.
Using our 3D methodology on a sample of 182 SFGs (selected solely from their [O II] emission) with M_{⋆} > 10^{8} M_{⊙} at 0.4 < z < 1.5 from the 9 arcmin^{2} MUSE observations of Bacon et al. (2017), we have found the following:
− Our methodology allows us to determine galaxy sizes of unresolved SFGs with down to R_{1/2,[O II]} ≃ 1 kpc or 0.2″ (Fig. 10) thanks to the builtin use of the PSF in our modeling. This opens a new parameter space (Fig. 10, right) between 1 < R_{e} < 3 kpc (and 10^{8} < M_{⋆} < 10^{9.5} M_{⊙}) unexplored so far with IFS surveys beyond the local universe (e.g., Burkert et al. 2016; Swinbank et al. 2017; Harrison et al. 2017).
− We find that the relation for SFGs is a nonlinear relation (Fig. 11, left) with at M_{⋆} < 10^{9.5} M_{⊙} and at M_{⋆} > 10^{9.5} M_{⊙}, similarly to RF12. This quadratic j − M relation appears to follow the expectations from simple disk formation models with the observed stellar mass fractions f_{m} (Moster et al. 2010; Behroozi et al. 2013). The zeropoint of the sAM sequence at z = 1 is found to be for rotation dominated log M_{⋆} = 10 SFGs both in the UDF and TNG50 samples.
− The relation is a strong function of the dynamical status of SFGs quantified by V/σ_{t}, leading to a sequence (Fig. 11, right), where rotation dominated SFGs (with V/σ_{t} > 1–5) follow the j − M relation of local disks, while dispersion dominated galaxies (with V/σ_{t} = 0.1) follow the j − M relation of early types. This confirms the preliminary results of Contini et al. (2016) on a small sample and extends the results of Burkert et al. (2016) to the relatively unexplored regime below 10^{10} M_{⊙} down to 10^{8} M_{⊙}.
− This j − M − (V/σ_{t}) relation is reminiscent of the z = 0 results where the j − M relation is found to be a function of the bulgetototal B/T ratio in local galaxies (e.g., Obreschkow & Glazebrook 2014; Fall & Romanowsky 2018) and is in relatively good qualitative agreement with theoretical framework on AM loss increasing with B/T in the EAGLE simulations (Zavala et al. 2016; Lagos et al. 2017) and with the results of Genel et al. (2015) from the Illustris simulations. Consequently, SFGs experience a dynamical transformation (from the disk to early type sequence of the j − M relation) and lower their sAMs prior to becoming passive, on the earlytype j − M sequence.
− The ZP of the j − M − (V/σ_{t}) relation in the UDF sample (Fig. 11, right) evolves as ∝(1 + z)^{a} with (2σ), which is consistent with the expectation of (1 + z)^{−0.5}, when taking into account the redshift evolution of V/σ_{t} (∝(1 + z)^{−0.5}) for our definition of σ_{t} = σ(2R_{e}).
− The ZP of the j − M − (V/σ_{t}) relation in the TNG50 sample (Fig. 12, right) evolves as ∝(1 + z)^{−0.45 ± 0.25}, which is consistent with the expectation of (1 + z)^{−0.5}, after taking into account the redshift evolution of V/σ_{t}, which goes as ∝(1 + z)^{−0.85} for our definition of σ_{t} = σ(2R_{e}). The j − M relation in the TNG50 sample is steeper with V/σ_{t}. This is likely to a combination of the difficulty in measuring σ_{t} in observations down to 10 km s^{−1} and to the physics detail in TNG50.
− The galaxy λ spin parameters (from the scatter of the j − M − (V/σ) relation) is broader in the MUSE sample than in TNG50 (Fig. 13), and also correlates better with galaxy sizes (SFRweighted) in TNG50 than in the UDF observations.
− The MUSE sample restricted to S/N_{eff} > 5 appears (Fig. 13, left) to show a more bimodal distributions in sAM, consistent with the z = 0 RF12 results where rotation dominated galaxies have retained most of their halo sAM, while dispersion dominated systems have lost around 80% of their sAM.
We have shown that both in the UDF MUSE mosaic and in TNG50, SFGs follow a j − M − (V/σ) fundamental relation, and have highlighed the importance of taking into account the redshift evolution of V/σ in the redshift evolution of the sAM sequence j − M.
Our results demonstrate that it is possible to study the kinematics of large samples of SFGs blindly (i.e., without preselection on morphology, nor prior HST continuum detection) with little biases including the lowS/N regime of 3–5. Our 3D modeling technique opens the lowmass regime, M_{⋆} < 10^{9} M_{⊙}, at intermediate redshifts 0.5 < z < 1.5, which had been unexplored so far.
See Marasco et al. (2019) for details on the conversion between j_{gas} and j_{⋆}.
This ratio is slightly massdependent beyond 10^{10} M_{⊙} reaching 1.2–1.3 at 10^{11} M_{⊙} in the KMOS3D survey (Wilman et al. 2020).
Available at http://galpak3d.univlyon1.fr.
Others uses j(< R_{e}) (Cortese et al. 2016) or even j(< R_{max}) where R_{max} is the last radius where Σ_{⋆}(r) can be measured (e.g., Posti et al. 2018a; Marasco et al. 2019; Mancera Piña et al. 2021), which is prone to biases.
This approach follows in spirit the calculations in Appendix A of Pillepich et al. (2019). However, it is important to note that Pillepich et al. (2019), in method (ii) of their appendix, have inadvertently added σ_{th} ≈ 15.7 km s^{−1} to their lineofsight velocity dispersions, while for a T = 10^{4} K ideal monoatomic Hydrogen gas that value corresponds to the 3D velocity dispersion, as implemented in this work.
When using the alternative definition of σ_{t} ≡ σ_{0} (discussed in Sect. 2.3) used in Übler et al. (2019), we find that q ≃ −0.7 to −0.8 very similar to the value of −0.8 we derived from the Übler et al. (2019) data (their Fig. 6).
For Swinbank et al. (2017), we assumed a mean V/σ_{t} = 2.5 for their V/σ_{t} > 1 sample according to their sample properties. For Romanowsky & Fall (2012), we assumed a V/σ = 5 for their diskdominated sample and shifted their j_{⋆} zeropoint to our SFRweighted by −0.25 dex.
Acknowledgments
This work made use of the following open source software: matplotlib (Hunter 2007), NumPy (Van Der Walt et al. 2011), SciPy (Jones et al. 2001), Astropy (Astropy Collaboration 2018). We thank the referee for his/her constructive comments that improved the quality of the paper. We thank L. Posti for insightful discussions that led to an improved manuscript. This study is based on observations collected at the European Southern Observatory under ESO programs 097.A0138(A), 097.A0144(A), 0100.A0089(A). This work has been carried out thanks to the support of the ANR 3DGasFlows (ANR17CE310017), the OCEVU Labex (ANR11LABX0060). B. E. acknowledges financial support from the program National Cosmology et Galaxies (PNCG) of CNRS/INSU with INP and IN2P3, cofunded by CEA and CNES. R. B. acknowledges support from the ERC advanced grant 339659MUSICOS. S. G., through the Flatiron Institute, is supported by the Simons Foundation. The primary TNG simulations were realized with compute time granted by the Gauss Centre for Super computing (GCS): TNG50 under GCS LargeScale Project GCS DWAR (2016; PIs Nelson/Pillepich), and TNG100 and TNG300 under GCSILLU (2014; PI Springel) on the GCS share of the supercomputer Hazel Hen at the High Performance Computing Center Stuttgart (HLRS). GCS is the alliance of the three national supercomputing centres HLRS (Universität Stuttgart), JSC (Forschungszentrum Jülich), and LRZ (Bayerische Akademie der Wissenschaften), funded by the German Federal Ministry of Education and Research (BMBF) and the German State Ministries for Research of BadenWürttemberg (MWK), Bayern (StMWFK) and NordrheinWestfalen (MIWF).
References
 AbrilMelgarejo, V., Epinat, B., Mercier, W., et al. 2021, A&A, 647, A152 [Google Scholar]
 Alcorn, L. Y., Tran, K.V., Glazebrook, K., et al. 2018, ApJ, 858, 47 [Google Scholar]
 Astropy Collaboration (PriceWhelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
 Bacchini, C., Fraternali, F., Iorio, G., & Pezzulli, G. 2019, A&A, 622, A64 [Google Scholar]
 Bacon, R., Accardo, M., Adjali, L., et al. 2010, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, SPIE Conf. Ser., 7735, 8 [Google Scholar]
 Bacon, R., Conseil, S., Mary, D., et al. 2017, A&A, 608, A1 [Google Scholar]
 Barcons, X., Lanzetta, K. M., & Webb, J. K. 1995, Nature, 376, 321 [Google Scholar]
 Behroozi, P. S., Wechsler, R. H., & Conroy, C. 2013, ApJ, 770, 57 [Google Scholar]
 Bett, P., Eke, V., Frenk, C. S., et al. 2007, MNRAS, 376, 215 [Google Scholar]
 Bett, P., Eke, V., Frenk, C. S., Jenkins, A., & Okamoto, T. 2010, MNRAS, 404, 1137 [NASA ADS] [Google Scholar]
 Bezanson, R., van der Wel, A., Pacifici, C., et al. 2018, ApJ, 858, 60 [Google Scholar]
 Boogaard, L. A., Brinchmann, J., Bouché, N., et al. 2018, A&A, 619, A27 [Google Scholar]
 Bouché, N., Cresci, G., Davies, R., et al. 2007, ApJ, 671, 303 [Google Scholar]
 Bouché, N., Murphy, M. T., Kacprzak, G. G., et al. 2013, Science, 341, 50 [Google Scholar]
 Bouché, N., Carfantan, H., Schroetter, I., MichelDansac, L., & Contini, T. 2015, AJ, 150, 92 [Google Scholar]
 Bouché, N., Finley, H., Schroetter, I., et al. 2016, ApJ, 820, 121 [Google Scholar]
 Bryan, G. L., & Norman, M. L. 1998, ApJ, 495, 80 [Google Scholar]
 Burkert, A., Genzel, R., Bouché, N., et al. 2010, ApJ, 725, 2324 [Google Scholar]
 Burkert, A., Förster Schreiber, N. M., Genzel, R., et al. 2016, ApJ, 826, 214 [Google Scholar]
 Butler, K. M., Obreschkow, D., & Oh, S.H. 2017, ApJ, 834, L4 [Google Scholar]
 Catinella, B., Giovanelli, R., & Haynes, M. P. 2006, ApJ, 640, 751 [Google Scholar]
 Ceverino, D., Primack, J., Dekel, A., & Kassin, S. A. 2017, MNRAS, 467, 2664 [Google Scholar]
 Chabrier, G. 2003, PASP, 115, 763 [Google Scholar]
 Cole, J., Bezanson, R., van der Wel, A., et al. 2020, ApJ, 890, L25 [Google Scholar]
 Contini, T., Epinat, B., Bouché, N., et al. 2016, A&A, 591, A49 [Google Scholar]
 Cortese, L., Fogarty, L. M. R., Bekki, K., et al. 2016, MNRAS, 463, 170 [Google Scholar]
 Cresci, G., Hicks, E. K. S., Genzel, R., et al. 2009, ApJ, 697, 115 [Google Scholar]
 Daddi, E., Elbaz, D., Walter, F., et al. 2010, ApJ, 714, L118 [Google Scholar]
 Dalcanton, J. J., & Stilp, A. M. 2010, ApJ, 721, 547 [Google Scholar]
 Dalcanton, J. J., Spergel, D. N., & Summers, F. J. 1997, ApJ, 482, 659 [Google Scholar]
 Danovich, M., Dekel, A., Hahn, O., Ceverino, D., & Primack, J. 2015, MNRAS, 449, 2087 [Google Scholar]
 Davies, R., Förster Schreiber, N. M., Cresci, G., et al. 2011, ApJ, 741, 69 [Google Scholar]
 DeFelippis, D., Genel, S., Bryan, G. L., & Fall, S. M. 2017, ApJ, 841, 16 [Google Scholar]
 DeFelippis, D., Genel, S., Bryan, G. L., et al. 2020, ApJ, 895, 17 [Google Scholar]
 DeFelippis, D., Bouché, N. F., Genel, S., et al. 2021, ApJ, submitted, [arXiv:2102.08383] [Google Scholar]
 Di Teodoro, E. M., & Fraternali, F. 2015, MNRAS, 451, 3021 [Google Scholar]
 Di Teodoro, E. M., Fraternali, F., & Miller, S. H. 2016, A&A, 594, A77 [Google Scholar]
 Dutton, A. A., & van den Bosch, F. C. 2009, MNRAS, 396, 141 [Google Scholar]
 Dutton, A. A., van den Bosch, F. C., & Dekel, A. 2010, MNRAS, 405, 1690 [NASA ADS] [Google Scholar]
 Efstathiou, G., & Jones, B. J. T. 1979, MNRAS, 186, 133 [Google Scholar]
 ElBadry, K., Quataert, E., Wetzel, A., et al. 2018, MNRAS, 473, 1930 [Google Scholar]
 Emsellem, E., Cappellari, M., Krajnović, D., et al. 2007, MNRAS, 379, 401 [Google Scholar]
 Epinat, B., Tasca, L., Amram, P., et al. 2012, A&A, 539, A92 [Google Scholar]
 Fall, S. M. 1983, IAU, 100, 391 [Google Scholar]
 Fall, S. M., & Efstathiou, G. 1980, MNRAS, 193, 189 [Google Scholar]
 Fall, S. M., & Romanowsky, A. J. 2013, ApJ, 769, L26 [Google Scholar]
 Fall, S. M., & Romanowsky, A. J. 2018, ApJ, 868, 133 [Google Scholar]
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
 Förster Schreiber, N. M., Genzel, R., Lehnert, M. D., et al. 2006, ApJ, 645, 1062 [Google Scholar]
 Genel, S., Vogelsberger, M., Springel, V., et al. 2014, MNRAS, 445, 175 [Google Scholar]
 Genel, S., Fall, S. M., Hernquist, L., et al. 2015, ApJ, 804, L40 [Google Scholar]
 Genel, S., Nelson, D., Pillepich, A., et al. 2018, MNRAS, 474, 3976 [Google Scholar]
 Genzel, R., Burkert, A., Bouché, N., et al. 2008, ApJ, 687, 59 [Google Scholar]
 Genzel, R., Tacconi, L. J., GraciaCarpio, J., et al. 2010, MNRAS, 407, 2091 [Google Scholar]
 Genzel, R., Schreiber, N. M. F., Übler, H., et al. 2017, Nature, 543, 397 [Google Scholar]
 Genzel, R., Price, S. H., Übler, H., et al. 2020, ApJ, 902, 98 [Google Scholar]
 Geweke, J. 1992, In Bayesian Statistics (University Press), 169 [Google Scholar]
 Gillman, S., Tiley, A. L., Swinbank, A. M., et al. 2020, MNRAS, 492, 1492 [Google Scholar]
 Graham, A. W., Driver, S. P., Petrosian, V., et al. 2005, AJ, 130, 1535 [Google Scholar]
 Grand, R. J. J., van de Voort, F., Zjupa, J., et al. 2019, MNRAS, 490, 4786 [Google Scholar]
 Guérou, A., Krajnović, D., Epinat, B., et al. 2017, A&A, 608, A5 [Google Scholar]
 Harrison, C. M., Johnson, H. L., Swinbank, A. M., et al. 2017, MNRAS, 467, 1965 [Google Scholar]
 Ho, S. H., Martin, C. L., Kacprzak, G. G., & Churchill, C. W. 2017, ApJ, 835, 267 [Google Scholar]
 Ho, S. H., Martin, C. L., & Turner, M. L. 2019, ApJ, 875, 54 [Google Scholar]
 Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [Google Scholar]
 Inami, H., Bacon, R., Brinchmann, J., et al. 2017, A&A, 608, A2 [Google Scholar]
 Jiang, F., Dekel, A., Kneller, O., et al. 2019, MNRAS, 488, 4801 [Google Scholar]
 Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open Source Scientific Tools for Python [Google Scholar]
 Kassin, S. A., Weiner, B. J., Faber, S. M., et al. 2007, ApJ, 660, L35 [Google Scholar]
 Kassin, S. A., Weiner, B. J., Faber, S. M., et al. 2012, ApJ, 758, 106 [NASA ADS] [CrossRef] [Google Scholar]
 Kennicutt, R. C. 1998, ARA&A, 36, 189 [Google Scholar]
 Kretschmer, M., Agertz, O., & Teyssier, R. 2020, MNRAS, 497, 4346 [Google Scholar]
 Kriek, M., van Dokkum, P. G., Franx, M., Illingworth, G. D., & Magee, D. K. 2009, ApJ, 705, L71 [Google Scholar]
 Lagos, C., Theuns, T., Stevens, A. R. H., et al. 2017, MNRAS, 3870, 3850 [Google Scholar]
 Lang, P., Meidt, S. E., Rosolowsky, E., et al. 2020, ApJ, 897, 122 [Google Scholar]
 Lapi, A., Salucci, P., & Danese, L. 2018, ApJ, 859, 2 [Google Scholar]
 Leroy, A. K., Walter, F., Brinks, E., et al. 2008, AJ, 136, 2782 [Google Scholar]
 Li, J., Obreschkow, D., Lagos, C., et al. 2020, MNRAS, 493, 5024 [Google Scholar]
 Lopez, S., Tejos, N., Ledoux, C., et al. 2018, Nature, 554, 493 [Google Scholar]
 Macciò, A. V., Dutton, A. A., van den Bosch, F. C., et al. 2007, MNRAS, 378, 55 [Google Scholar]
 Mancera Piña, P. E., Posti, L., Fraternali, F., Adams, E. A. K., & Oosterloo, T. 2021, A&A, 647, A76 [CrossRef] [EDP Sciences] [Google Scholar]
 Marasco, A., Fraternali, F., Posti, L., et al. 2019, A&A, 621, L6 [Google Scholar]
 Marinacci, F., Vogelsberger, M., Pakmor, R., et al. 2018, MNRAS, 480, 5113 [Google Scholar]
 Martin, C. L., Ho, S. H., Kacprzak, G. G., & Churchill, C. W. 2019, ApJ, 878, 84 [Google Scholar]
 Martinsson, T. P. K., Verheijen, M. A. W., Westfall, K. B., et al. 2013, A&A, 557, A130 [Google Scholar]
 Meurer, G. R., Carignan, C., Beaulieu, S. F., & Freeman, K. C. 1996, AJ, 111, 1551 [Google Scholar]
 Mitchell, P. D., Lacey, C. G., Lagos, C. D. P., et al. 2018, MNRAS, 474, 492 [Google Scholar]
 Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 [Google Scholar]
 Moster, B. P., Somerville, R. S., Maulbetsch, C., et al. 2010, ApJ, 710, 903 [Google Scholar]
 Mowla, L., van der Wel, A., van Dokkum, P., & Miller, T. B. 2019, ApJ, 872, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Murugeshan, C., Kilborn, V., Jarrett, T., et al. 2020, MNRAS, 496, 2516 [Google Scholar]
 Naiman, J. P., Pillepich, A., Springel, V., et al. 2018, MNRAS, 477, 1206 [Google Scholar]
 Navarro, J. F., & Steinmetz, M. 2000, ApJ, 538, 477 [Google Scholar]
 Nelson, E. J., van Dokkum, P. G., Förster Schreiber, N. M., et al. 2016, ApJ, 828, 27 [Google Scholar]
 Nelson, D., Pillepich, A., Springel, V., et al. 2018, MNRAS, 475, 624 [Google Scholar]
 Nelson, D., Pillepich, A., Springel, V., et al. 2019, MNRAS, 490, 3234 [Google Scholar]
 Obreschkow, D., & Glazebrook, K. 2014, ApJ, 784, 26 [Google Scholar]
 Obreschkow, D., Glazebrook, K., Bassett, R., et al. 2015, ApJ, 815, 97 [Google Scholar]
 Obreschkow, D., Glazebrook, K., Kilborn, V., & Lutz, K. 2016, ApJ, 824, L26 [Google Scholar]
 Oh, S., Colless, M., Barsanti, S., et al. 2020, MNRAS, 495, 4638 [Google Scholar]
 Osborne, C., Salim, S., Damjanov, I., et al. 2020, ApJ, 902, 77 [Google Scholar]
 Peebles, P. J. E. 1969, ApJ, 155, 393 [Google Scholar]
 Persic, M., Salucci, P., & Stel, F. 1996, MNRAS, 281, 27 [Google Scholar]
 Pillepich, A., Nelson, D., Hernquist, L., et al. 2018a, MNRAS, 475, 648 [Google Scholar]
 Pillepich, A., Springel, V., Nelson, D., et al. 2018b, MNRAS, 473, 4077 [Google Scholar]
 Pillepich, A., Nelson, D., Springel, V., et al. 2019, MNRAS, 490, 3196 [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [Google Scholar]
 Posti, L., Pezzulli, G., Fraternali, F., & Di Teodoro, E. M. 2018a, MNRAS, 475, 232 [Google Scholar]
 Posti, L., Fraternali, F., Di Teodoro, E. M., & Pezzulli, G. 2018b, A&A, 612, L6 [Google Scholar]
 Posti, L., Marasco, A., Fraternali, F., & Famaey, B. 2019, A&A, 629, A59 [Google Scholar]
 Prugniel, P., & Simien, F. 1997, A&A, 321, 111 [NASA ADS] [Google Scholar]
 Refregier, A., Kacprzak, T., Amara, A., Bridle, S., & Rowe, B. 2012, MNRAS, 425, 1951 [Google Scholar]
 Renzini, A. 2020, MNRAS, 495, L42 [Google Scholar]
 Rizzo, F., Fraternali, F., & Iorio, G. 2018, MNRAS, 476, 2137 [Google Scholar]
 Romanowsky, A. J., & Fall, S. M. 2012, ApJS, 203, 17 [Google Scholar]
 Romeo, A. B. 2020, MNRAS, 491, 4843 [Google Scholar]
 Romeo, A. B., & Mogotsi, K. M. 2018, MNRAS, 480, L23 [Google Scholar]
 Salvatier, J., Wiecki, T. V., & Fonnesbeck, C. 2016, PeerJ Comput. Sci., 2 [Google Scholar]
 Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521 [Google Scholar]
 Sérsic, J. L. 1963, Boletin de la Asociacion Argentina de Astronomia La Plata Argentina, 6, 41 [Google Scholar]
 Sharma, S., Steinmetz, M., & BlandHawthorn, J. 2012, ApJ, 750, 107 [Google Scholar]
 Sharples, R., Bender, R., Agudo Berbel, A., et al. 2013, Messenger, 151, 21 [Google Scholar]
 Somerville, R. S., Behroozi, P., Pandya, V., et al. 2018, MNRAS, 473, 2714 [Google Scholar]
 Springel, V. 2010, MNRAS, 401, 791 [Google Scholar]
 Springel, V., & Hernquist, L. 2003, MNRAS, 339, 312 [Google Scholar]
 Springel, V., Pakmor, R., Pillepich, A., et al. 2018, MNRAS, 475, 676 [Google Scholar]
 Steidel, C. C., Kollmeier, J. A., Shapley, A. E., et al. 2002, ApJ, 570, 526 [Google Scholar]
 Stevens, A. R. H., Croton, D. J., & Mutch, S. J. 2016, MNRAS, 461, 859 [Google Scholar]
 Stewart, K. R., Brooks, A. M., Bullock, J. S., et al. 2013, ApJ, 769, 74 [Google Scholar]
 Stewart, K. R., Maller, A. H., Oñorbe, J., et al. 2017, ApJ, 843, 47 [Google Scholar]
 Stott, J. P., Swinbank, A. M., Johnson, H. L., et al. 2016, MNRAS, 457, 1888 [Google Scholar]
 Straatman, C. M. S., Glazebrook, K., Kacprzak, G. G., et al. 2017, ApJ, 839, 57 [Google Scholar]
 Sweet, S. M., Fisher, D., Glazebrook, K., et al. 2018, ApJ, 860, 37 [Google Scholar]
 Swinbank, A. M., Harrison, C. M., Trayford, J., et al. 2017, MNRAS, 467, 3140 [NASA ADS] [Google Scholar]
 Tacchella, S., Diemer, B., Hernquist, L., et al. 2019, MNRAS, 487, 5416 [Google Scholar]
 Teklu, A. F., Remus, R.S., Dolag, K., et al. 2015, ApJ, 812, 29 [Google Scholar]
 Tiley, A. L., Bureau, M., Cortese, L., et al. 2019, MNRAS, 482, 2166 [Google Scholar]
 Übler, H., Förster Schreiber, N. M., Genzel, R., et al. 2017, ApJ, 842, 121 [Google Scholar]
 Übler, H., Genzel, R., Wisnioski, E., et al. 2019, ApJ, 880, 48 [Google Scholar]
 van den Bosch, F. C., Abel, T., & Hernquist, L. 2003, MNRAS, 346, 177 [Google Scholar]
 Van Der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22 [Google Scholar]
 van der Wel, A., Chang, Y.Y., Bell, E. F., et al. 2014, ApJ, 792, L6 [Google Scholar]
 van de Sande, J., Kriek, M., Franx, M., et al. 2013, ApJ, 771, 85 [Google Scholar]
 Vogelsberger, M., Sijacki, D., Kereš, D., Springel, V., & Hernquist, L. 2012, MNRAS, 425, 3024 [Google Scholar]
 Vogelsberger, M., Genel, S., Sijacki, D., et al. 2013, MNRAS, 436, 3031 [Google Scholar]
 Vogelsberger, M., Genel, S., Springel, V., et al. 2014a, Nature, 509, 177 [Google Scholar]
 Vogelsberger, M., Genel, S., Springel, V., et al. 2014b, MNRAS, 444, 1518 [Google Scholar]
 Vogelsberger, M., Marinacci, F., Torrey, P., & Puchwein, E. 2020, Nat. Rev. Phys., 2, 42 [Google Scholar]
 Weinberger, R., Springel, V., Hernquist, L., et al. 2017, MNRAS, 465, 3291 [Google Scholar]
 Whitaker, K. E., Franx, M., Leja, J., et al. 2014, ApJ, 795, 104 [Google Scholar]
 White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 [Google Scholar]
 Wilman, D. J., Fossati, M., Mendel, J. T., et al. 2020, ApJ, 892, 1 [Google Scholar]
 Wisnioski, E., Förster Schreiber, N. M., Wuyts, S., et al. 2015, ApJ, 799, 209 [Google Scholar]
 Wisnioski, E., Förster Schreiber, N. M., Fossati, M., et al. 2019, ApJ, 886, 124 [Google Scholar]
 Zabl, J., Bouché, N. F., Schroetter, I., et al. 2019, MNRAS, 485, 1961 [Google Scholar]
 Zavala, J., Frenk, C. S., Bower, R., et al. 2016, MNRAS, 460, 4466 [Google Scholar]
Appendix A: Analytical derivation
For a Sérsic profile I(r)∝exp[−b_{n} (R/R_{e})^{1/n}], we can calculate the coefficient k_{n} in Eq. (3) assuming a constant flat rotation curve such as V(R) = V_{c} as in RF12. Using Eq. (2) and by substituting by , we find
which means that the k_{n} coefficient in Eq. (3) is
where Γ is the gamma function. As discussed in Graham et al. (2005), the Sérsic coefficient b_{n} is given from the solution of the following equation
where Γ is the (complete) gamma function and γ is the incomplete gamma function. Note, the Sérsic coefficient b_{n} can be also approximated by b_{n} ≈ 1.9992 n − 0.327 as discussed in Graham et al. (2005) (see also Prugniel & Simien 1997).
Appendix B: Selection function
Due to the instrumental resolution of R ∼ 3000 (∼ 100 km s^{−1} or ∼ 50 km s^{−1} per pixel) for MUSE, we are limited to a few tens of km s^{−1} in either V_{max} or σ_{t}. In order to quantify the impact of this limit to the j − M − (V/σ) relation, , we assume a constant floor K for , which becomes
after a change of variables with and K=r^{2}. This implies that the j − M − (V/σ) relation becomes
Finally, the last term can be expressed as a function of M_{⋆} using the S05M_{⋆} relation (e.g., Kassin et al. 2007; Übler et al. 2017; Tiley et al. 2019), namely log S05 ≃ 0.3 log M_{⋆,9} + 1.5, and .
Appendix C: Kinematics
Fig. C.1 shows examples of kinematic properties for UDF mosaic galaxies analyzed in § 4.3.
Fig. C.1.
Morphology and kinematics of examples from the UDFmosaic sample. The panels show the flux map, the SB profile Σ_{SFR}, the velocity field, the rotation curve and the SFRweighted profiles, respectively. The 1D profiles are determined along the pseudolong slits represented by the gray straight lines (black solid lines), along with the 95% confidence intervals. The vertical dotted lines represent 2 × R_{e}. In the last column, the horizontal blue lines show the total modeled angular momentum j_{3D}, while the horizontal red lines show the total using the Eq. (3) approximation. 
Fig. C.2 shows examples of kinematic properties for TNG50 galaxies analyzed in § 3.4.
Fig. C.2.
Kinematic properties of some of the TNG50 galaxies. The panels show the flux map, the SB profile Σ_{SFR}, the velocity field, the rotation curve and the SFRweighted profiles, respectively. The 1D profiles are determined along the pseudolong slits represented by the gray straight lines (black solid lines), along with the 95% confidence intervals. The vertical dotted lines represent 2 × R_{e}. The 1D red solid lines represent the true profile determined directly from the TNG50 data. In the last column, the horizontal blue lines show the total modeled angular momentum j_{3D}, while the horizontal red lines show the total using the Eq. (3) formula. This approximation can lead to an overestimation of the total angular momentum by 10–30% depending on the steepness of the v(r) profile. 
Appendix D: Redshift evolution of V/σ
In order to estimate the redshift evolution q of V/σ_{t}, we minimized the V/σ_{t} distribution with respect to q. Fig. D.1 shows the results of this excercise for the TNG50 (UDF) sample in the left (right) panel, respectively. This figure shows that V/σ_{t} evolves as (1 + z)^{−0.8} in TNG50 and as (1 + z)^{−0.5} in the UDF sample.
Fig. D.1.
Left:V/σ standard distribution as a function of the redshift evolution q (Eq. (10)). One sees that q ∼ −0.8 minimizes the distribution across all redshifts. Right: Same for the UDF sample. The redshift evolution for V/σ here is q ≃ −0.5. 
All Tables
All Figures
Fig. 1.
Systematics at play in measuring j(< r). Left: sAM profile j(< r) normalized by the RF12 approximation (Eq. (3)) as a function of radius for disks with Sérsic n = 1 and a tanh rotation curve. The solid blue horizontal line represents the Eq. (3) approximation. The solid curves represent the sAM for various rotation curves with various turnover radii R_{t} where the opacity of the line is inversely proportional to R_{t}/R_{e}. The vertical dotted line represents 2R_{e}. For galaxies with slowly rising rotation curves, with larger R_{t}, the RF12 approximation leads to a large overestimation of the total angular momentum (represented by the horizontal dashed lines), by 10–20%. Right: relative difference for the k_{n} coefficient between the k_{n} formula from RF12 and Eq. (4) as a function of Sérsic index n. One sees that Eq. (3) under or overestimate k_{n} by up to ∼10% for 0.5 < n < 4. 

In the text 
Fig. 2.
Schematic description of our 3D methodology using a 2step process. First, we determine the morphological parameters by fitting a 2D Sérsic model for the flux profile Σ(r) on a [O II] image from MUSE. Then, we fit the MUSE [O II] data with 3D model with the Sérsic n fixed to 0.5, 1, 2 or 4 based on n_{2d} and priors on the halflight radius R_{e} and inclination i from the results of the 2D analysis. We run the MCMC algorithm in two steps, first with Cauchy sampling which reaches convergence more efficiently, then with EMCEE which better sample the posteriors. 

In the text 
Fig. 3.
Sizes and SFRs of the 158 galaxies selected from the overall TNG50 galaxy population. Left: sizemass relation for the selected TNG50 for the mocks where the sizes are halfSFR sizes (calculated in faceon projection and without any dust attenuation). The open (solid) squares represent the z = 0.5 (z = 1) galaxies where the sizes for the z = 0.5 galaxies are adjusted to z = 1 according to the (1 + z)^{−0.5} evolution found by van der Wel et al. (2014). The background gray points and contours represent the sizemass distribution for the z = 1 TNG50 galaxy population. The solid line (gray band) represents the sizemass relation for late types from Mowla et al. (2019) (van der Wel et al. 2014), respectively. Right: SFRmass or “mainsequence” relation for the selected TNG50 galaxies (red dots, where SFRs for the z = 0.5 galaxies are adjusted to z = 1 according to the (1 + z)^{2.2} evolution found by Whitaker et al. 2014). The background gray points and contours represent the SFRmass distribution for the z = 1 TNG50 galaxy population. For comparison, we show with the gray band the z = 1 mainsequence from the MUSE UDF analysis of Boogaard et al. (2018). 

In the text 
Fig. 4.
Visual representations of three example TNG50 starforming galaxies from the set analyzed in this work. The two columns on the left represent stellar light (dustfree color composite of the (SDSS)rg(Johnson)B bands) and the two on the right the SFR in the gas phase. Each panel is 30 kpc on a side. 

In the text 
Fig. 5.
Example of GAlPAK^{3D} fits for three SFGs from TNG50 shown in Fig. 4. The columns show the flux map, the flux profile Σ(r), the intrinsic velocity field, the deprojected velocity profile and the deprojected angular momentum profile, from left to right respectively. The 1D profiles are computed along the pseudoslit shown in the first and third columns by the gray lines. In the second column, the SB profile for a single orientation is shown (black points). In the second, fourth and fifth columns, the solid gray curves show the model fits for the various inclinations and orientations. The red solid lines represent the true profiles determined directly from the TNG50 data. The vertical dotted lines represent twice the halflight radius (R_{1/2,SFR}). The horizontal blue lines show the total modeled angular momentum j_{3D}, while the horizontal red lines show the total j from RF12 (Eq. (3)), which can lead to overestimation of the total angular momentum by 10–20%. 

In the text 
Fig. 6.
Reliability of our sAM estimator. Left: bottom panel: total SFRweighted sAM, (in km s^{−1} kpc), as a function of the true angular momentum (SFRweighted) for the 790 mock galaxies. Top panel: relative error ( as a function of the effective S/N (see text) along with the histogram of . The points are colored according to the stellar mass in both panels. The red circles represent the median whose errorbar represent 1σ. Right: relative error ( for each mass bin of 1 dex centered on log M_{⋆} = 8.5, 9.5, 10.5, from top to bottom. The histograms are weighted by the errors. Left (right) panels: distributions of , in three mass bins, resulting from the 3D method (the Eq. (3) approximation), respectively. The red solid line represents the continuous density estimation of the distributions using a Gaussian kernel with bandwidth of ∼0.2. This shows that, for massive galaxies (M_{⋆} > 10^{10} M_{⊙}), the Eq. (3) approximation is less biased, while for galaxies with M_{⋆} < 10^{10} M_{⊙}, the 3D method is less biased. 

In the text 
Fig. 7.
Relation between the stellar sAM j_{⋆} and for z = 1 TNG50 SFGs. The solid points with errorbars represent the median difference between the two as a function of the stellar angular momentum, where the error bars represent the standard deviation. Thus, the SFRweighted sAM, , is higher by ≈ 0.25 dex (dotted line) compared to j_{⋆} over a wide range of sAM from 50 to 2000 km s^{−1} kpc. A second order polynomial fit is shown with the dotdashed black line. 

In the text 
Fig. 8.
Left: following AbrilMelgarejo et al. (2021), we show the S/N of the total [O II] flux as a function of the resolved ratio (R_{e, hst}/PSF FWHM) for the sample of SFGs. The gray bands represent the selection of resolved galaxies in AbrilMelgarejo et al. (2021), consisting of S/N > 30 and R_{e, hst}/FWHM > 0.5. Right: comparison between the total [O II] S/N and the effective S/N (Eq. (6)), colorcoded with R_{e}. In both panels, the red points are galaxies with S/N_{eff} < 3.0, while the blue points correspond to galasxes with S/N_{eff} > 3.0. 

In the text 
Fig. 9.
Comparison of selection techniques. Left: sizemass relation for the SFGs in the UDF mosaic sample (blue circles). The sizes are adjusted to z = 1 according to the (1 + z)^{−0.5} evolution found by van der Wel et al. (2014). The solid line (gray band) represent the sizemass relation of Mowla et al. (2019) (van der Wel et al. 2014). Right: mainsequence relation for the SFGs in the UDF mosaic sample. The SFRs are adjusted to z = 1 according to the (1 + z)^{1.74} evolution of Boogaard et al. (2018) whose z = 1 main sequence is represented by the gray band. In both panels, the points are colorcoded according to effective S/N. For comparison we show the properties of the mock 158 galaxies (Sect. 3.2) as red circles. 

In the text 
Fig. 10.
Left: comparison between [O II]based sizes R_{e} from MUSE (at FWHM ≈ 0.7″ resolution) versus Hband continuum R_{e, ⋆} from HST/F160W at ≈0.15″ resolution, for the full sample of 494 SFGs, where the points are colorcoded according to Eq. (6). The top panel shows the relative error δR_{e}/R_{e} as a function of the effective S/N (see text) and the side histogram shows its distribution. The gray errorbars represent the 3σ_{t} statistical uncertainties. This shows that our MUSE[O II] sizes are of sufficient accuracy down to ≈1 kpc, provided that the effective S/N is ⪆3. Right: sizes R_{e} (from HST/F160W) as a function of redshift for 182 galaxies in the MUSE UDFmosaic field whose effective S/N is sufficient. The gray region represents the region probed by Swinbank et al. (2017). The dashed line represents ∝(1 + z)^{1}. The vertical error bars are 3σ. 

In the text 
Fig. 11.
Specific angular momentum for the 182 galaxies in the MUSE UDF mosaic fields with a max S/N ≳ 3. Left: total specific angular momentum (SFR weighted) as a function of stellar mass M_{⋆}. The points are colorcoded as a function of V/σ_{t}. The rightaxis shows the corresponding j_{⋆} according to the 0.25dex offset found in Fig. 7. The blue (red) dashed line show the sequence for late (early) type galaxies at z = 0 from Fall & Romanowsky (2018, FR18). The dotted line (solid squares) show the z = 0 sAM j_{⋆} − M_{⋆} sequence from Posti et al. (2018a) (Lapi et al. 2018), respectively. The dotted line show the sequence for z = 0 disk galaxies of Posti et al. (2018a). The solid black line shows the expected j_{⋆} − M_{⋆} relation from Eq. (8). The solid blue and red lines show the fitted fiducial model from Eq. (9) and the gray band represents the 95% confidence interval. Right: normalized AM sequence as a function of (V/σ_{t})_{toz1} (normalized to z = 1 see text), i.e. after taking into account for the redshift evolution of σ_{t}(z), where the points are colorcoded with redshift. In the inset, we show the redshift evolution of the sAM relation at M_{⋆} = 10^{10} M_{⊙} (solid line) at fixed (V/σ_{t})_{toz1} (=5) along with the 95% predictive interval (gray band), with the z = 0 ZP from Fall & Romanowsky (2018) (red square), the Swinbank et al. (2017) binned data (see text), and the expected (1 + z)^{−0.5} evolution (dotted line) from Eq. (8). 

In the text 
Fig. 12.
Same as Fig. 11, but for TNG50 SFGs. 

In the text 
Fig. 13.
Left: sAM sequence normalization as a function of (V/σ_{t})_{toz1} for the UDF mosaic sample of SFGs where the color representes the SFRweighted log R_{1/2,[O II]}. The left yaxis is found from the fiducial fits shown in Figs. 11–12 (see Table 1). The right yaxis shows the corresponding effective spin λ × f_{j} according to Eq. (8) using the zeropoing of Eq. (8). The solid histograms show the distributions from the data. The solid line shows the lognormal distribution P(λ) expected for halos (Bett et al. 2007). Right: same for the matched TNG50 sample of SFGs. 

In the text 
Fig. C.1.
Morphology and kinematics of examples from the UDFmosaic sample. The panels show the flux map, the SB profile Σ_{SFR}, the velocity field, the rotation curve and the SFRweighted profiles, respectively. The 1D profiles are determined along the pseudolong slits represented by the gray straight lines (black solid lines), along with the 95% confidence intervals. The vertical dotted lines represent 2 × R_{e}. In the last column, the horizontal blue lines show the total modeled angular momentum j_{3D}, while the horizontal red lines show the total using the Eq. (3) approximation. 

In the text 
Fig. C.2.
Kinematic properties of some of the TNG50 galaxies. The panels show the flux map, the SB profile Σ_{SFR}, the velocity field, the rotation curve and the SFRweighted profiles, respectively. The 1D profiles are determined along the pseudolong slits represented by the gray straight lines (black solid lines), along with the 95% confidence intervals. The vertical dotted lines represent 2 × R_{e}. The 1D red solid lines represent the true profile determined directly from the TNG50 data. In the last column, the horizontal blue lines show the total modeled angular momentum j_{3D}, while the horizontal red lines show the total using the Eq. (3) formula. This approximation can lead to an overestimation of the total angular momentum by 10–30% depending on the steepness of the v(r) profile. 

In the text 
Fig. D.1.
Left:V/σ standard distribution as a function of the redshift evolution q (Eq. (10)). One sees that q ∼ −0.8 minimizes the distribution across all redshifts. Right: Same for the UDF sample. The redshift evolution for V/σ here is q ≃ −0.5. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.