---
title: 'ALMA-CRISTAL: ISM Mapping at Cosmic Dawn'
url: https://www.emergentmind.com/topics/alma-cristal-program
type: topic
---

# ALMA-CRISTAL: ISM Mapping at Cosmic Dawn

ALMA-CRISTAL—“[CII] Resolved ISM in STar-forming galaxies with ALMA”—is an ALMA Large Program (ID 2021.1.00280.L) designed to spatially resolve the [CII] \(158\,\mu\mathrm{m}\) line and far-infrared dust continuum in typical main-sequence star-forming galaxies at \(4<z<6\), and to register those cold-ISM tracers against rest-UV and rest-optical structure from HST and JWST. The program targets the interval when the Universe was \(\sim 1\) Gyr old and reionization had just ended, with the explicit aim of building a kiloparsec-scale census of ISM structure, dust-obscured star formation, kinematics, outflows, and circumgalactic gas in representative massive galaxies rather than compact starburst-selected systems [2505.06340, 2507.11600].

## 1. Program definition and scientific rationale

CRISTAL is motivated by the need to move beyond small, heterogeneous, and often starburst-biased high-redshift samples. Its central tracer, the [CII] line at rest frequency \(\nu_{\rm rest}=1900.537\) GHz, is a dominant coolant of the multiphase ISM and can arise in photodissociation regions around young star clusters as well as in diffuse neutral and ionized gas. In the CRISTAL redshift range, the line is shifted into ALMA Band 7 and remains bright enough to trace gas kinematics and morphology to and beyond effective radii [2408.03374, 2507.11600].

Program papers organize the science around four primary goals: kinematics, outflows, spatial distribution of the ISM and star formation, and ISM conditions. The survey was designed to couple ALMA Band 7 spectroscopy and continuum imaging with HST/WFC3 rest-UV morphology and, in later work, JWST/NIRCam and JWST/NIRSpec observations. This combination enables direct comparison of stellar light, unobscured star formation, dust-embedded star formation, and cool gas geometry on \(\sim\)kpc scales [2505.06340, 2510.16111].

A defining feature of CRISTAL is its emphasis on typical main-sequence galaxies. The program explicitly contrasts with earlier ALMA work dominated by massive dusty star-forming galaxies and quasars. This shift in target class is essential because the representative galaxy population at \(z\sim 5\) governs most stellar-mass build-up, yet is difficult to characterize in the far-infrared and submillimeter regime. CRISTAL therefore functions both as a resolved ISM survey and as a calibration dataset for interpreting JWST-era studies of early galaxy assembly [2311.17671, 2505.06340].

## 2. Survey design, target selection, and sample growth

The initial targeted sample consisted of 19 main-sequence galaxies drawn from ALPINE. Selection required robust spectroscopic redshifts at \(4\le z\le 6\), [CII] detections in ALPINE, stellar masses \(\log(M_\star/M_\odot)\gtrsim 9.5\), HST imaging, and specific star formation rates within a factor of three of the star-forming main sequence. Several papers also note the practical role of HST/WFC3 UV coverage in ensuring that stellar and UV morphologies could be compared directly to ALMA maps [2505.06340, 2311.17671].

As the program developed, the working sample expanded through deeper CRISTAL detections, archival ALMA data, and pilot targets. Different program papers report different accounting schemes tied to their scientific scope: the resolved SED study spans 25 fields containing 36 galaxies, the outflow analysis refers to 37 galaxies including serendips and resolved multiples, and the survey overview reports 39 main-sequence galaxies after including newly detected galaxies in CRISTAL fields, archival data, and pilot study targets [2409.10961, 2504.17877, 2505.06340].

The survey overview gives the broadest census: 39 main-sequence galaxies with median \(z\approx 5.1\), median stellar mass \(M_\star \approx 10^{10.1}\,M_\odot\), and median \(\mathrm{SFR}\approx 58\,M_\odot\,\mathrm{yr}^{-1}\). [CII] is detected in 38/39 systems, with CRISTAL-18 as the non-detection, and dust continuum is detected in 24/39. The kinematic analysis narrows this to 32 targets with sufficient signal-to-noise ratio and spatial resolution, while the dedicated [CII]-size study reports structural parameters for \(N=34\) [CII]-emitting components after excluding the non-detection and treating interacting systems component-wise [2505.06340, 2507.11600, 2408.03374].

This sample construction matters methodologically. The survey was designed to be representative of the massive main-sequence population at \(4<z<6\), not of extreme compact starbursts. At the same time, several papers note that pre-selection on ALPINE [CII] detections, HST visibility, and UV brightness biases the survey against the most heavily obscured systems and toward well-characterized COSMOS and GOODS-S fields with rich ancillary data [2507.11600, 2505.06340].

## 3. Observational strategy and analysis framework

CRISTAL’s ALMA setup was optimized for both uv-coverage and spatial dynamic range. Targets were observed in at least two 12 m-array configurations: a compact configuration sensitive to large scales and an extended configuration sensitive to sub-kpc structure. Typical final [CII] cubes have synthesized beams of \(\theta_{\rm beam}\approx 0.3\)–\(0.5\) arcsec, corresponding to \(\sim 2\)–3 kpc at \(z\approx 5\), while higher-resolution products reach \(\sim 0.2\)–\(0.3\) arcsec in favorable cases. Continuum subtraction is performed in the uv-plane with CASA/uvcontsub, and imaging generally uses natural weighting for sensitivity, with Briggs products used where finer spatial detail is required [2408.03374, 2505.06340].

A recurrent methodological choice is direct visibility fitting. In the [CII]-size analysis, emission is modeled with UVMULTIFIT as a 2D exponential disk with parameters for position, total line flux, major-axis FWHM, axis ratio \(q\), and position angle. The surface-brightness law is
\[
I(r)=I_0\exp(-r/R_s),
\]
with
\[
R_e=1.678\,R_s,
\qquad
R_e=1.21\,\mathrm{FWHM}\,\sqrt{q}.
\]
This uv-plane approach is used to avoid CLEAN-mask choices, weighting trade-offs, and PSF deconvolution uncertainties. Rest-UV sizes are measured with GALFIT assuming \(n=1\), and FIR sizes are fit as circular exponentials in the uv-plane [2408.03374].

The kinematic papers adopt a separate forward-modeling framework. DysmalPy constructs full 3D model cubes including projection, beam convolution, and spectral broadening, and fits 1D velocity and dispersion profiles extracted along the kinematic major axis. The baryonic model comprises an exponential disk and, where needed, a compact de Vaucouleurs bulge; the dark matter halo follows an NFW profile. Pressure support is included through
\[
V_{\rm rot}^2(R)=V_{\rm circ}^2(R)-3.36\,\sigma_0^2\,(R/R_e),
\]
and inclinations are inferred from axis ratios using an intrinsic disk thickness \(\epsilon_0=0.25\) [2507.11600].

CRISTAL also introduced resolved UV-to-FIR SED modeling at \(\sim 0.5\)–1 kpc resolution. The 14-galaxy study uses PSF-matched HST, JWST, and ALMA data with magphys, while the JWST PRIMER study of four galaxies uses BAGPIPES pixel-by-pixel fits to NIRCam photometry. These analyses quantify \(\Sigma_\star\), \(\Sigma_{\rm SFR}\), \(L_{\rm dust}\), \(A_V\), and related maps, and show explicitly that ALMA continuum data reduce uncertainties in \(L_{\rm dust}\), \(A_V\), and SFR by breaking age–dust degeneracies [2409.10961, 2409.10963].

## 4. Spatial distribution of gas, dust, and star formation

The most systematic morphological result is that [CII] emission is more extended than both stellar UV light and, to a lesser degree, FIR dust emission. The dedicated size census finds \(R_{e,[\mathrm{CII}]}\) in the range \(0.5\)–\(3.5\) kpc, with average \(\langle R_{e,[\mathrm{CII}]}\rangle = 1.90\pm 0.77\) kpc. The mean size ratios are \(\langle R_{e,[\mathrm{CII}]}/R_{e,\mathrm{UV}}\rangle = 2.90\pm 1.40\) and \(\langle R_{e,[\mathrm{CII}]}/R_{e,\mathrm{FIR}}\rangle = 1.54\pm 0.63\). Median [CII]-to-UV peak offsets are only \(\sim 0.11\) arcsec, implying co-located centers but more extended [CII] outskirts [2408.03374].

A key negative result is that the achieved data do not require a bright secondary [CII] halo component. Residual visibilities show no statistically significant excess beyond an exponential disk, and stacked short-baseline residuals remain below \(\lesssim 1.2\sigma\). In the stacked fit of single galaxies, a one-component model gives \(R_{e,[\mathrm{CII}]}=1.97\pm 0.10\) kpc, while a two-component “galaxy + halo” fit yields a very faint, poorly constrained extended component with flux fraction \(8.7\pm 13.4\%\) and \(R_{e,\mathrm{halo}}=9.17\pm 28.2\) kpc, with reduced \(\chi^2\) indistinguishable from the single-disk fit [2408.03374].

The dust-continuum program reaches a complementary conclusion. In 26 galaxies, 19 are individually detected in rest-frame \(158\,\mu\mathrm{m}\) continuum, with \(\log_{10}L_{\rm IR}[L_\odot]=10.9\)–12.4. Dust effective radii are on average \(\sim 1.5\) kpc and are \(\sim 1.8\) times larger than the rest-UV. The median infrared surface density is \(\Sigma_{\rm IR}\sim 2.0\times 10^{10}\,L_\odot\,\mathrm{kpc}^{-2}\), about an order of magnitude below compact DSFG starbursts, supporting a picture of disk-wide rather than nuclear star formation in typical \(z=4\)–6 galaxies [2311.17671].

Resolved SED work further shows that CRISTAL galaxies obey a resolved star-forming main sequence at \(z\sim 5\). Combining all 14 galaxies yields a best-fit slope \(0.85\pm 0.02\) with median scatter \(\simeq 0.1\) dex, and the summed stellar masses and SFRs from spatial bins agree with integrated SED fits, implying no systematic outshining bias in these data. The same study finds that less obscured star formation is located where \(\Sigma_\star\) and \(\Sigma_{\rm SFR}\) are highest, typically in central regions, and explicitly speculates that feedback may be driving gas and dust out of those regions [2409.10961].

The program also ties [CII] extent to ISM phase balance. The [CII]-size paper reports a negative correlation between \(\Sigma_{[\mathrm{CII}]}\) and \(\mathrm{EW}_{\mathrm{Ly}\alpha}\) with Spearman \(\rho=-0.51\) and \(p=0.039\), a positive correlation between \(\Sigma_{[\mathrm{CII}]}/\Sigma_{\mathrm{IR}}\) and \(\mathrm{EW}_{\mathrm{Ly}\alpha}\) with \(\rho=0.56\) and \(p=0.090\), and a possible negative correlation between \(R_{e,[\mathrm{CII}]}/R_{e,\mathrm{UV}}\) and \(\mathrm{EW}_{\mathrm{Ly}\alpha}\). The interpretation advanced there is that inner [CII] primarily traces PDRs directly associated with star formation, while the extra spatial extent likely includes a substantial diffuse neutral atomic-gas contribution [2408.03374].

Resolved ISM-condition work on HZ10 extends this logic to component scales. Multi-band ALMA SED fitting gives a global \(T_{\rm SED}=46.7\pm 6.8\) K, while the more UV-obscured HZ10-W component has \(T_{\rm SED}=51.2\pm 13.1\) K and a lower [CII]/FIR ratio than HZ10-C and the Bridge. HZ10-C follows the local UV-starburst IRX–\(\beta_{\rm UV}\) relation, whereas HZ10-W and the Bridge are better described by dust-screen models with holes in front of a hard UV radiation field [2407.09681].

## 5. Kinematics, interactions, and circumgalactic structures

CRISTAL’s resolved kinematic analysis provides the first systematic kpc-scale characterization of typical massive main-sequence galaxies at \(4<z<6\). In the 32-galaxy kinematic sample, \(22\pm 7\%\) are classified as Best Disks, \(28\pm 8\%\) as Disks, and \(50\pm 9\%\) as Non-Disks, giving a total disk fraction of \(50\pm 9\%\). The disks have median intrinsic velocity dispersion \(\sigma_0\approx 70\,\mathrm{km\,s^{-1}}\), typical \(V_{\rm rot}(R_e)\approx 100\)–\(180\,\mathrm{km\,s^{-1}}\), and median \(V_{\rm rot}/\sigma_0\approx 2\), indicating highly turbulent, moderately rotation-supported systems. Roughly half the disks have \(f_{\rm DM}(<R_e)\lesssim 30\%\), and the median enclosed dark-matter fraction is \(\approx 18\%\) [2507.11600].

These disk properties are interpreted within a marginally Toomre-stable, gas-rich framework. Using CRISTAL medians \(\sigma_0/v_c\approx 0.35\)–0.4 and \(f_{\rm molgas}\approx 0.5\), the kinematic paper infers typical \(Q_{\rm gas}\approx 0.6\), with the sample spanning \(Q\approx [0.4,2.0]\). The elevated dispersions are described as consistent with dominant gravity driving, and the weak or absent correlation of \(\sigma_0\) with \(\Sigma_{\rm SFR}\) argues against feedback alone as the dominant driver of turbulence across the disk sample [2507.11600].

A recurrent lesson of the program is that apparently simple systems often resolve into interacting or multi-component structures at higher angular resolution. HZ10 (CRISTAL-22) at \(z=5.65\), studied with ALMA and JWST/NIRSpec, separates into at least three components—HZ10-C, HZ10-E, and HZ10-W—with projected separations of \(\sim 1.5\) and 4 kpc from the central component. The [CII] morphology and kinematics closely resemble the broad [OIII] \(\lambda 5007\) component, leading to dynamical scenarios ranging from a disturbed rotating disk plus companions to triple or quadruple merger configurations [2411.09033].

CRISTAL-05/HZ3 at \(z=5.54\) provides a second representative example. High-resolution [CII] maps show that a source previously treated as single is instead a close interacting pair, C05-NW and C05-SE, separated by \(\sim 0.4\) arcsec (\(\sim 2.4\) kpc), and embedded in extended carbon-enriched gas. The global [CII] radius is \(r_{[\mathrm{CII}]}\approx 4\,r_{\mathrm{UV}}\), significant [CII] emission extends to \(\sim 10\) kpc, and the dominant component is fit by a rotation-supported, baryon-dominated DysmalPy model with only a minor dark-matter contribution within the radial range probed. Spatially resolved [CII]/FIR ratios in the outskirts are consistent with merger-driven shocks or turbulence enhancing [CII] emission [2403.03379].

CRISTAL also resolves extreme circumgalactic structures. In the J1000+0234 system at \(z=4.54\), the program identified an elongated [CII]-emitting plume extending \(\sim 15\) kpc from the central DSFG, with no clear counterpart at other wavelengths. Its centroid velocity increases from \(\sim 180\) to \(\sim 400\,\mathrm{km\,s^{-1}}\) relative to systemic along the plume, while the linewidth narrows from \(\sim 400\)–600 to \(\sim 190\,\mathrm{km\,s^{-1}}\). The paper evaluates four possible origins—conical outflow, cold accretion stream, ram-pressure stripping, and gravitational interactions—and explicitly disfavors ram-pressure stripping while leaving the others viable [2401.04919].

Outflows are present in the program, but not as a dominant explanation for all extended [CII]. A stacking analysis of 15 disk-like galaxies, after excluding systems with disturbed kinematics or obvious interactions, finds only weak evidence for a broad [CII] component in the composite spectrum, and that evidence is largely driven by CRISTAL-02. Interpreting the modest residuals at \(v\sim 300\,\mathrm{km\,s^{-1}}\) as an outflow yields \(\dot{M}_{\rm out}=26\pm 11\,M_\odot\,\mathrm{yr}^{-1}\) and a cold atomic mass-loading factor \(\eta_m=0.49\pm 0.20\). The high-\(\Sigma_{\rm SFR}\) subsample shows moderate evidence for a broad component, whereas the low-\(\Sigma_{\rm SFR}\) subsample does not [2504.17877].

## 6. Program significance, limitations, and evolving legacy

At the program level, CRISTAL has established a quantitative picture of the cold ISM in typical galaxies during the first billion years. Its core empirical synthesis is that [CII] in these systems is usually well described by a single extended exponential disk rather than by a statistically required bright halo, that [CII] is systematically larger than rest-UV and moderately larger than FIR continuum, and that a substantial fraction of the emission likely arises outside the main star-forming radius in diffuse neutral gas. In parallel, the kinematic program shows that many \(z\sim 5\) main-sequence galaxies already host turbulent disks, often in interacting environments, with low dark-matter fractions inside \(R_e\) and substantial baryonic dominance on galaxy scales [2408.03374, 2507.11600].

The survey also demonstrates the value of strict multi-wavelength co-registration. The ALPINE-CRISTAL-JWST program combines CRISTAL ALMA maps with JWST/NIRSpec IFU and JWST/HST imaging for 18 representative galaxies at \(4.4\le z\le 5.7\), delivering co-spatial maps of H\(\alpha\), H\(\beta\), [OII], [OIII], [NII], [SII], stellar continuum, [CII], and dust on \(\sim 1\)–2 kpc scales. Early integrated results from that program show gas-phase metallicities \(12+\log(\mathrm{O/H})\approx 7.7\)–8.5, broad agreement among H\(\alpha\), UV+IR, [OII], and [CII]-based SFR indicators, and stacked H\(\alpha\) profiles that are more extended than rest-UV, while [CII] is more extended than H\(\alpha\) [2510.16111].

Several limitations are emphasized across the literature. Angular resolution remains finite, with many analyses operating at \(\theta_{\rm beam}\approx 0.3\)–0.5 arcsec and occasional need to fix the [CII] axis ratio \(q\). Sensitivity at the shortest baselines is reduced even when compact configurations are included. Morphological classification into single, pair, multiple-UV, or multiple-[CII] systems is explicitly described as crude, and JWST often reveals substructure in objects previously treated as single. Single-band dust measurements introduce uncertainty in \(L_{\rm IR}\), \(T_{\rm dust}\), and \(f_{\rm obs}\), while sample pre-selection on ALPINE detections and HST visibility biases the survey toward UV-bright systems [2408.03374, 2311.17671, 2505.06340].

Future work follows directly from these limitations. Program papers explicitly call for 3D analyses of Ly\(\alpha\) and H\(\alpha\) with JWST/NIRSpec IFUs, higher-resolution ALMA kinematics, direct [NII] \(205\,\mu\mathrm{m}\) measurements to quantify ionized-gas [CII] fractions, deeper and higher-frequency ALMA continuum to constrain dust temperatures, and broader multi-phase comparisons among [CII], H\(\alpha\), [OIII], and molecular gas tracers. In that sense, CRISTAL has become more than a single ALMA survey: it is the cold-ISM anchor of a wider observational program linking ALMA, HST, JWST, and IFU spectroscopy into a resolved benchmark for early galaxy assembly [2505.06340, 2510.16111].

Source: https://www.emergentmind.com/topics/alma-cristal-program