CAVITY: Void Galaxy Evolution Survey

• CAVITY is a legacy integral-field spectroscopy survey that targets void galaxies to probe environmental influences on galaxy evolution.
• It uses coordinated optical, radio, and millimeter observations to derive spatially resolved star-formation histories, metallicity gradients, and structural properties.
• The survey’s design enables detailed comparisons of internal processes versus external effects, challenging ΛCDM predictions in low-density cosmic voids.

The Calar Alto Void Integral-field Treasury surveY (CAVITY) is a legacy integral-field spectroscopy (IFS) project dedicated to the systematic characterization of galaxies residing within the most underdense regions of the local Universe—the cosmic voids. By leveraging multi-wavelength spatially resolved spectroscopy and coordinated optical, radio, and millimeter follow-up, CAVITY provides an unprecedented benchmark for constraining the impact of low-density environments on galaxy mass assembly, star-formation histories, chemical enrichment, and structural evolution. The survey targets a statistically significant sample of several hundred void galaxies, supporting detailed analysis of internal processes and large-scale environmental influences on galaxy evolution (García-Benito et al., 2024, Pérez et al., 2024, Azevedo et al., 12 Jan 2026, Conrado et al., 2024).

1. Scientific Motivation and Survey Objectives

CAVITY is designed to address key outstanding questions in extragalactic astrophysics associated with environmental regulation of galaxy evolution. Void galaxies, relatively isolated from the strong tidal forces and frequent interactions that prevail in clusters and filaments, represent an ideal laboratory to disentangle internal processes—including secular evolution and feedback—from environmental effects such as mergers and ram-pressure stripping (García-Benito et al., 2024, Pérez et al., 2024).

The primary scientific objectives include:

• Mapping spatially resolved star-formation histories, stellar populations, ionized gas properties, and dark matter content in void galaxies.
• Testing inside-out versus outside-in stellar mass assembly scenarios within minimal external perturbation environments.
• Comparing radial gradients in stellar age, metallicity, and specific star-formation rate (sSFR) between void galaxies and analogues in higher-density regimes.
• Assessing how baryonic and dark matter relations and scaling laws (e.g., mass–metallicity relation, Tully–Fisher relation) behave at the extremes of cosmic density.
• Constraining the efficiency and timescales of internal quenching mechanisms (e.g., AGN feedback) in the absence of frequent external interactions.

These objectives collectively aim to challenge and refine theoretical predictions of the ΛCDM paradigm, especially regarding the role of environment at small (sub-megaparsec) scales (Pérez et al., 2024, Conrado et al., 2024).

2. Survey Design, Sample Selection, and Observational Strategy

CAVITY systematically draws its targets from the SDSS DR7 void catalogue using the Pan et al. (2012) algorithm. The parent sample is defined as galaxies occupying underdense regions: $0.005 < z < 0.050$, within $0.8\,R_{\rm void}$ of void centers, absolute $r$-band magnitudes $-21.5 \leq M_r \leq -17.0$, and stellar masses $10^{8.5}$–$10^{11}\,M_\odot$ (García-Benito et al., 2024, Pérez et al., 2024, Conrado et al., 2024). To ensure high-quality data, galaxies are selected for intermediate inclinations (20–70°), well-determined morphological types (E, S0, Sa, Sb, Sc, Sd using the Domínguez-Sánchez et al. 2018 scheme), and sufficient surface brightness.

The final CAVITY sample covers 15 representative voids with at least 20 members each; cross-survey overlaps enable comparison with MaNGA, CALIFA, and WALLABY. The first data release (DR1) provides 100 science-grade data cubes, with full completion targeting $\sim$300 galaxies (García-Benito et al., 2024).

Observations are carried out on the 3.5 m telescope at Calar Alto Observatory using the PMAS/PPak IFU. The instrument configuration employs V500 (R~850, 3745–7500 Å, FWHM ~6 Å) for all galaxies and V1200 (R~1650, 3400–4750 Å) for brighter targets, using dithered exposures to achieve complete spatial coverage and minimize vignetting (García-Benito et al., 2024, Pérez et al., 2024). On-source exposure times are typically 1.5–3 h per galaxy.

3. Data Processing, Calibration, and Data Products

The reduction pipeline, adapted from the CALIFA framework and written in Python3/Cython, performs the following steps (García-Benito et al., 2024):

• CCD Preprocessing: Bias subtraction, cosmic-ray removal via PyCosmic, multi-amplifier frame combination.
• Fiber Extraction and Tracing: Flexure correction, optimal extraction of 331–382 science fibers per exposure, stray-light subtraction.
• Wavelength Calibration: HeHgCd arc lamps and correction for flexure, resampling spectra to 2 Å px⁻¹ linear grids, homogenization to the nominal FWHM.
• Flat-Field and Throughput: Twilight flats determine fiber transmission; fibers below 70% relative throughput are flagged.
• Flux Calibration: Master sensitivity curves (per García-Benito et al. 2015), nightly standards, atmospheric and telluric correction, absolute photometric scaling to SDSS $g, r$ 30″ diameter apertures (targeted $\leq$5% absolute accuracy, $\Delta$($g-r$)$\leq$0.05 mag).
• Sky Subtraction: Median combination of the 30 faintest of 36 sky fibers.
• Cube Construction: Dithered row-stacked spectra re-scaled, stacked, extinction-corrected, and sampled onto a $1''\times 1''$ grid. Differential atmospheric refraction (DAR) correction and WCS registration are applied.
• Quality Control and Masking: Vignetting flagged, cosmic rays and CCD defects masked (with BADPIX and ERRWEIGHT extensions), error propagation maintained throughout.

The principal science product is a FITS datacube with five extensions: flux, error, error weights, bad-pixel mask, and fiber coverage per spaxel. Metadata include calibration frames, observing conditions, and QA tables. Ancillary data products comprise master catalogs of physical and environmental galaxy properties, PPAK field overlays, and online cross-match capability (Simbad/VizieR). Access is provided via a web portal and TAP/ADQL interface (García-Benito et al., 2024).

CAVITY+ (the project extension) incorporates deep INT $g,r$ imaging, IRAM 30m and ALMA molecular gas (CO) data, and GBT/HI spectroscopy for subsets of the sample, allowing for a holistic baryon census in voids (Pérez et al., 2024).

4. Analysis Methods and Measurement Approaches

Data processing for stellar population analysis utilizes full spectral fitting codes, predominately STARLIGHT (for large samples) and FADO (for detailed studies, especially those requiring nebular emission modeling) (Conrado et al., 2024, Azevedo et al., 12 Jan 2026).

Key analytical steps:

• Spatial Binning: Signal-to-noise optimized tessellation via Voronoi binning to S/N$\sim$20 or higher per spatial bin.
• Spectrum Fitting: Population synthesis using MILES-based SSPs (ages 1 Myr–14 Gyr, multiple metallicities); direct fitting for light and mass fraction population vectors; self-consistent dust attenuation estimation.
• Stellar and Nebular Diagnostics: Derivation of stellar mass, mass surface density, light-weighted age, star formation rate (SFR; e.g., $t_{\rm SF}=32$ Myr), and sSFR. Emission-line diagnostics applied for nebular properties and AGN identification (BPT, WHAN diagrams).
• Radial Profiles and Morphometry: Profiles of $\Sigma_*(R)$, $\langle \log t \rangle_L$, $\Sigma_{\rm SFR}(R)$, and sSFR($R$) are constructed by scaling to each galaxy's half-light radius (HLR). Morphological analysis includes CASGM parameterization (concentration, asymmetry, Gini, $M_{20}$) and quantification of tidal and asymmetric features (Azevedo et al., 12 Jan 2026).
• Error Analysis: Covariance introduced by cube interpolation is empirically calibrated; noise properties validated via residual analysis of spectral fits.

Spatial mapping of quantities (stellar age, metallicity, SFR) is further refined using Integrated Nested Laplace Approximation (INLA), providing physically smooth 2D reconstructions (Azevedo et al., 12 Jan 2026).

5. Key Empirical Results and Early Science

Major findings from initial data releases and pilot analyses include (Conrado et al., 2024, García-Benito et al., 2024, Azevedo et al., 12 Jan 2026):

• Structural Properties: Void galaxies have marginally larger HLRs (+0.05 dex in $\log\mathrm{HLR}$) and systematically lower stellar mass surface densities ($\Delta\log\Sigma_*\approx-0.17$ dex).
• Stellar Populations: All morphological types in voids are younger by $\Delta\langle\log t\rangle_L\approx-0.19$ dex. This age difference is most pronounced in the outer discs ($R>1\,\mathrm{HLR}$), consistent with more gradual evolution.
• Star Formation: Global SFR and sSFR are mildly enhanced in voids (+0.15 to +0.17 dex), statistically significant for early-type spirals (Sa). Outer discs of void spirals have higher sSFR, often remaining star-forming where filament analogues have quenched (sSFR$(R)>0.1\rm\,Gyr^{-1}$).
• Radial Trends: For low-mass galaxies $(M_*<10^{10.5}\,M_\odot)$, environmental effects (gas accretion modes, reduced merger rates) dominate the persistence of star-forming, extended discs; for higher masses, internal processes (e.g., mass quenching via AGN) dominate, and environmental effects are diminished.
• Interaction Case Studies: Triplet systems in voids (e.g., CAVITY5273X and VGS31) exhibit diverse assembly histories: ongoing star formation, tidal features, and AGN activity are evident, indicating that local interactions and filamentary inflows can drive significant evolutionary episodes even in very low-density environments (Azevedo et al., 12 Jan 2026).
• Mass-Metallicity Relation: Void galaxies adhere to a tight $M_*$–O/H relation ($$12+\log(\mathrm{O/H})=8.75+0.18\log(M_*/10^9\,M_\odot)$$), though interacting systems or those with ongoing accretion can deviate by $\sim$0.2–0.3 dex (Azevedo et al., 12 Jan 2026).
• Morphological Disturbance: A significant fraction ($\sim$26%) of observed void galaxies display low-surface-brightness tidal features, with central metallicity depressions linked to recent interactions.

These empirical results establish the distinctiveness of void galaxy evolution and provide rigorous comparison points for theoretical modeling (García-Benito et al., 2024, Conrado et al., 2024).

6. Data Access, Ancillary Products, and Limitations

CAVITY data releases are delivered through a dedicated web portal and ADQL/TAP service, offering download links, quality assurance tables, and cross-referenced value-added catalogs (García-Benito et al., 2024). The five-extension FITS datacubes include flux, error, error-weighting, bad-pixel flags, and fiber coverage layers.

Ancillary products encompass master catalogs (with $\sim$30 columns per galaxy: SDSS identifiers, photometry, stellar/gas properties, environmental metrics), footprint overlays, and multi-wavelength follow-up (HI, CO, deep optical imaging) via the CAVITY+ extension (Pérez et al., 2024).

Key data characteristics:

• Spatial resolution: $2.5''$ (seeing + fiber), corresponding to $\sim$1–2 kpc at $z\sim0.03$.
• Flux calibration: absolute accuracy $\lesssim$5%, robust across all fields.
• S/N: $>50$ per Å in central regions; $>3$ in outskirts (with binning).
• Limitations: DR1 under-samples the faint ($M_r>-18$ mag) and low-mass ($M_*<10^9\,M_\odot$) regime; small galaxies may show low-SB dithering patterns, but central regions are unaffected.

Table: CAVITY Data Structure (DR1)

Extension	Content	Units / Description
PRIMARY	Flux density	$10^{-16}$ erg s$^{-1}$ cm$^{-2}$ Å$^{-1}$
ERROR	1$\sigma$ uncertainty	same as PRIMARY
ERRWEIGHT	Error-weight factors	dimensionless
BADPIX	Bad-pixel flags (1=bad, 0=good)	per spaxel, per wavelength
FIBCOVER	Number of fibers contributing	integer (1–3)

7. Scientific Impact and Outlook

CAVITY, with its extension CAVITY+, is poised to have a transformative impact on the study of environmental influences on galaxy evolution. The survey directly addresses the long-standing question of how large-scale structure modulates star formation, chemical enrichment, and morphological transformation. Its legacy sample is designed for direct comparison with CALIFA, MaNGA, SAMI, and future HI/CO surveys, providing powerful leverage for modeling feedback and assembly in isolated versus grouped regimes (García-Benito et al., 2024, Conrado et al., 2024, Pérez et al., 2024).

A plausible implication is that void galaxies evolve more gradually and with less efficient quenching, especially in their outer regions, than their wall and filament counterparts; however, local interactions and filamentary accretion still play a dynamic role, refuting a simplistic picture of passivity. The spatially resolved, multi-wavelength dataset forms a comprehensive reference for confronting galaxy formation models under the extremes of cosmic density.