Nburst: Integral-Field Spectrum Fitting

• Nburst is a full-spectrum stellar-population framework that fits integral-field spectra using single-stellar-population templates convolved with line-of-sight velocity distributions.
• It employs a joint fitting technique that simultaneously optimizes stellar age, metallicity, kinematics, and nebular emission through chi-square minimization.
• The derived mass-to-light ratios and stellar masses from Nburst are crucial for interpreting galaxy dynamics and interactions in complex environments like Abell 2142.

Searching arXiv for papers on Nburst to ground the article in published sources. Nburst is a full-spectrum stellar-population fitting framework used to model integral-field spectra with one or more single-stellar-population templates convolved with line-of-sight velocity distributions, while simultaneously accounting for emission lines and continuum-shape mismatches. In the Abell 2142 infalling-group study, it is the key tool for decomposing multiple galaxies superposed along the line of sight in MaNGA integral-field spectroscopy and for deriving approximate single-stellar-population ages and metallicities; those parameters are then converted into SDSS $r$-band mass-to-light ratios and stellar masses that feed directly into the system’s dynamical interpretation (Shaji et al., 10 Jul 2025).

1. Conceptual role and fitting formalism

In the usage documented for the Abell 2142 system, Nburst follows a standard full-spectrum philosophy. Each spaxel spectrum is fit by a model built from one or more SSP templates, convolved with a LOSVD, with nebular emission and continuum residuals handled in the same optimization. For the one-component analysis, each spaxel is represented by a single SSP template from the e-MILES library, and the fit searches over stellar age, metallicity, and kinematics simultaneously (Shaji et al., 10 Jul 2025).

The fit is expressed as a minimization of

$$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$

Here, $p_\lambda$ masks the good wavelength pixels and $\Delta F_\lambda$ is the spectral uncertainty. The stellar spectrum is interpolated on the SSP grid at each trial pair $(T_{\rm SSP}, [Z/H]_{\rm SSP})$, then convolved with a Gaussian LOSVD. This makes Nburst, in this context, a joint stellar-population and kinematic inference engine rather than a line-index method or a narrow-feature fitting procedure.

A central methodological point is that the stellar-population parameters are not treated as isolated descriptive quantities. The ages and metallicities returned by Nburst are subsequently transformed into mass-to-light ratios and stellar masses, so the outputs of the spectral fit become inputs to the astrophysical interpretation of the group as a dynamical system.

2. Spectral assumptions and instrumental treatment

Several assumptions define the specific Nburst implementation used in the Abell 2142 analysis. The stellar continuum in each fitted component is approximated by a single SSP from the e-MILES library. For the one-component fit, the LOSVD is initially taken to be a pure Gaussian. Strong nebular lines are included as additive components with Gaussian kinematics, but they are fit independently from the stellar kinematics, and their fluxes enter linearly at each nonlinear iteration (Shaji et al., 10 Jul 2025).

Continuum-shape mismatches are absorbed with a 19th-degree multiplicative Legendre polynomial. The reported choice is roughly “1 degree per 150–200 Å” across the full spectral range. In practice, this polynomial is used to absorb residual calibration errors and dust reddening rather than to encode the stellar-population signal itself.

The instrumental setup is also explicit. The MaNGA cube spans $3600$–$10300$ Å at $R \sim 2000$, and the analysis uses that broad range in a true full-spectrum sense rather than restricting the fit to selected diagnostics. Both the SSP templates and the emission-line templates are pre-convolved with the wavelength-dependent MaNGA instrumental resolution before fitting. This makes the comparison between model and data resolution-consistent across the full bandpass.

These assumptions impose a well-defined interpretation on the returned parameters. The ages and metallicities are approximate SSP-equivalent quantities for each recovered stellar component, not nonparametric star-formation histories. That limitation is methodological rather than incidental, because the decomposition problem is already highly degenerate in the presence of overlapping galaxies.

3. Two-component decomposition in overlapping systems

The most technically distinctive use of Nburst in the Abell 2142 study is the two-component analysis applied to overlapping galaxies. The system contains a merging pair, $a_1$ and $a_2$, together with additional projected contributions from galaxies $$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$0 and $$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$1. In such spaxels, the observed spectrum is modeled as the sum of two independent stellar populations, each with its own LOSVD (Shaji et al., 10 Jul 2025).

The authors do not run this two-component fit blindly on the full cube. Instead, they use a constrained workflow designed to reduce degeneracy:

1. approximate kinematics are first recovered from a non-parametric LOSVD analysis;
2. each galaxy’s rotation field is modeled with a simple circular rotation law;
3. those velocities are used as starting values;
4. only a few high-S/N spaxels per configuration are fit, typically 3–4 spaxels per case;
5. fits are repeated with different combinations of components to avoid false $$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$2 minima.

This restricted strategy is motivated by the coarse spatial resolution of MaNGA relative to the sizes of the galaxies, which means that many spaxels contain mixed light from multiple systems. In that regime, unconstrained multi-component fitting would be especially vulnerable to local minima and component swapping. The two-component Nburst analysis is therefore not a brute-force decomposition of the entire field, but a targeted inference procedure anchored by external kinematic information.

A plausible implication is that, in this application, Nburst is most powerful when embedded in a broader inference pipeline rather than used as a standalone black-box fitter. The spectral model, the velocity-field model, and the spaxel selection strategy are tightly coupled.

4. Returned stellar-population parameters and derived masses

The two-component fits yield mean stellar-population parameters for each galaxy, which are then converted into SDSS $$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$3-band mass-to-light ratios and stellar masses (Shaji et al., 10 Jul 2025).

Galaxy	Age and metallicity	$$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$4 and stellar mass
$$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$5	age $$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$6 Gyr; $$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$7	$$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$8; $$\chi^2 = \sum_{N_\lambda} p_\lambda \, \frac{\left(F_\lambda - F_\lambda^\mathrm{model}\right)^2}{\Delta F_\lambda^2}.$$9
$p_\lambda$0	age $p_\lambda$1 Gyr; $p_\lambda$2	$p_\lambda$3; $p_\lambda$4
$p_\lambda$5	age $p_\lambda$6 Gyr; $p_\lambda$7	$p_\lambda$8; $p_\lambda$9
$\Delta F_\lambda$0	age $\Delta F_\lambda$1 Gyr; $\Delta F_\lambda$2	$\Delta F_\lambda$3; $\Delta F_\lambda$4

These outputs are central because they connect spectroscopy to dynamical interpretation. The Nburst-derived ages and metallicities provide the basis for estimating $\Delta F_\lambda$5, and the inferred $\Delta F_\lambda$6 values are combined with photometric decomposition to obtain stellar masses. In the paper’s analysis, Nburst is therefore not merely a classification tool; it supplies quantitative priors on the stellar content of each overlapping galaxy.

The distribution of fitted ages is also astrophysically informative. The particularly young fitted age of $\Delta F_\lambda$7, $\Delta F_\lambda$8 Gyr, is used together with the emission-line analysis to support a picture of ongoing or recent star formation in a dynamically disturbed system.

5. Astrophysical interpretation in the Abell 2142 infalling group

The Nburst results are integrated with stellar and gas kinematics, photometry, and emission-line diagnostics to characterize a compact, infalling, strongly interacting galaxy group projected on a filament about $\Delta F_\lambda$9 Mpc from the Abell 2142 center. The galaxies are described as perturbed systems with tidal tails and loops; they are mainly disks in rotation, although some regions show elevated dispersion characteristic of out-of-equilibrium gas (Shaji et al., 10 Jul 2025).

A major conclusion is that the group is not quenched. All galaxies show sustained star formation, with a global star formation rate of $(T_{\rm SSP}, [Z/H]_{\rm SSP})$0 after AGN correction. The Nburst-derived populations are relatively young overall, and especially young in the central merging galaxy $(T_{\rm SSP}, [Z/H]_{\rm SSP})$1. In context, these results support the interpretation that galaxy-galaxy interactions within the group are still enhancing star formation rather than suppressing it.

The stellar-population measurements also help distinguish between environmental mechanisms. The paper concludes that the spectacular trailing $(T_{\rm SSP}, [Z/H]_{\rm SSP})$2-kpc X-ray tail did not come from stripping of the individual galaxy disks’ own gas, but from the hot intra-group medium that was present before infall. Nburst contributes indirectly to that conclusion because the fitted populations, along with the derived stellar masses, support a picture of gas-rich, star-forming disks that are perturbed and interacting but not yet stripped to the degree required to account for the entire X-ray feature.

This suggests that, in this application, Nburst functions as part of a pre-processing diagnosis. Its outputs support an interpretation in which internal group interactions remain dynamically and star-formation relevant even as the larger cluster environment has already altered the surrounding hot gas reservoir.

6. Scope, ambiguity, and cross-domain usage

The term “Nburst” is not uniform across the supplied arXiv literature. In the Abell 2142 study, it denotes a spectral-fitting and decomposition framework for stellar populations in integral-field spectroscopy (Shaji et al., 10 Jul 2025). In a separate all-optical satellite-networking study, by contrast, “Nburst” is not explicitly defined as a named symbol in the text; the relevant quantity there is the burst-size threshold $(T_{\rm SSP}, [Z/H]_{\rm SSP})$3, meaning the amount of user traffic assembled into one optical burst switching burst at the optical ground station before transmission (Roethig et al., 25 Jun 2026).

A further source of ambiguity is lexical rather than conceptual. The neutron-star literature discusses thermonuclear burst oscillations, which are highly asymmetric brightness patterns on the burning surface layers of accreting neutron stars and are used as probes of spin, dense matter, and surface dynamics; this is unrelated to the stellar-population decomposition role of Nburst in extragalactic spectroscopy (Watts, 2012).

This cross-domain variation suggests that “Nburst” should be interpreted contextually rather than treated as a universal notation. In the specific astrophysical usage documented here, its distinguishing features are full-spectrum fitting, SSP-based decomposition, LOSVD modeling, and the derivation of SSP-equivalent ages, metallicities, mass-to-light ratios, and stellar masses from mixed integral-field spectra.