Spectral Core: A Multi-Domain Overview

• Spectral core is a polysemous technical label across various fields, defining central structures inferred from spectral or frequency-domain data.
• In VLBI studies, it denotes the frequency-dependent τ=1 surface whose position shifts reveal key jet parameters like magnetic fields and particle densities.
• In network science and spectroscopy, spectral core methods optimize core–periphery detection and predict core-level spectra through nonlinear eigenproblems and state-specific protocols.

“Spectral core” is not a single standardized technical term across the arXiv literature. It denotes several distinct objects that are all defined by combining a notion of core with a spectral or frequency-domain characterization. In VLBI studies of relativistic jets, the core is the synchrotron self-absorbed $\tau=1$ surface whose position and spectrum vary with observing frequency and with flaring state (Fromm et al., 2013). In network science, the spectral core is the set of nodes with the largest entries in a nonlinear eigenvector that scores core–periphery structure (Tudisco et al., 2018). In radio-source population studies, the spectral core is the compact, relativistically boosted flat-spectrum central component whose prominence is quantified by the core-dominance parameter (Pei et al., 2019). In molecular electronic-structure theory, “Spectral-Core” is used for a protocol based on Restricted Open-Shell Kohn-Sham theory and Square Gradient Minimization for predicting core-level spectra (Hait et al., 2019). These usages are conceptually related by their emphasis on extracting a central structure from spectral observables, but they are not interchangeable.

1. Terminological scope

The term appears in at least four research contexts, each with its own ontology, observables, and mathematical machinery.

Domain	Meaning of “spectral core”	Representative source
Relativistic-jet VLBI	Frequency-dependent $\tau=1$ radio core and its spectral turnover structure	(Fromm et al., 2013)
Network science	High-coreness nodes from a nonlinear spectral fixed point	(Tudisco et al., 2018)
Radio-source statistics	Doppler-boosted flat-spectrum compact core in two-component AGN models	(Pei et al., 2019)
Molecular spectroscopy	ROKS/SGM protocol for core-excited spectra at DFT cost	(Hait et al., 2019)

In the jet literature, “core” is explicitly not a fixed physical knot but the surface along the jet where the synchrotron optical depth reaches unity, $\tau=1$; at frequencies below this surface the jet is opaque, and above it the jet is transparent (Fromm et al., 2013). In the network literature, the term refers to a continuous coreness score whose largest entries define the network’s spectral core (Tudisco et al., 2018). In AGN population studies, it refers to the compact radio component with approximately flat spectral index, revealed most clearly in high-$R$ sources (Pei et al., 2019). In molecular spectroscopy, the phrase labels a computational protocol rather than a physical substructure (Hait et al., 2019).

A common misconception is that “spectral core” names a unified theory. The literature instead shows a family of domain-specific usages. A plausible implication is that the phrase is best treated as a polysemous technical label rather than as a single canonical concept.

2. Spectral core in relativistic-jet VLBI

In VLBI images of relativistic jets, the observed radio core is the $\tau=1$ surface. Because synchrotron self-absorption depends on frequency, this surface moves upstream at higher frequencies and downstream at lower frequencies. Quantitatively, the distance of the core from the jet base at observing frequency $\nu$ follows

$r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$

In the simplest conical-jet, equipartition case one finds $k_r=1$; if $B(r)\propto r^{-b}$ and particle density $N(r)\propto r^{-n}$, then

$\tau=1$0

where $\tau=1$1 is the optically thin spectral index defined by $\tau=1$2 (Fromm et al., 2013).

CTA 102 provides a detailed case study. During the 2006 flare, analysis of eight multi-frequency VLBI observations detected a radio component ejected around 2005.9 and a stationary feature around $\tau=1$3 mas away from the core; the overall picture was interpreted as an over-pressured jet perturbed by traveling shock waves (Fromm et al., 2013). A related multi-epoch VLBA analysis spanning May 2005 to April 2007 reported that the position of the jet core is proportional to $\tau=1$4 with some temporal variations, suggesting possible equipartition between magnetic field energy and particle kinetic energy densities at the most compact regions (Fromm et al., 2013).

The core-shift law is not immutable. In CTA 102 during the 2006 flare, the eight-epoch average gave

$\tau=1$5

from fits of

$\tau=1$6

The slight departure from $\tau=1$7 signals a shift away from pure equipartition, indicating that the balance between magnetic and particle energy densities is evolving during the flare (Fromm et al., 2013).

A common misunderstanding is to identify the VLBI core with a stationary emission knot. The cited analyses explicitly treat it as a dynamic opacity surface whose position and spectrum encode local jet physics, especially during flaring states (Fromm et al., 2013).

3. Core-shift, spectral turnover, and intrinsic jet parameters

Measurement of the jet spectral core requires map co-registration, because standard VLBI calibration loses absolute sky coordinates. The standard procedure is to cross-correlate optically thin jet regions—bright, extended features whose positions do not depend on frequency—to determine relative shifts, using a two-dimensional, intensity-weighted cross-correlation (Fromm et al., 2013). After alignment, the fitted position of the $\tau=1$8 core in each map is measured relative to a reference frequency and fitted with the core-shift model above.

After core-shift correction, all frequency maps are convolved to a common beam and pixel grid. At each pixel along the jet ridge-line, the spectrum is fitted with a homogeneous synchrotron self-absorbed model,

$\tau=1$9

with

$\tau=1$0

The fit returns the turnover frequency $\tau=1$1 and turnover flux $\tau=1$2; the optically thick spectral index is fixed at $\tau=1$3 for a homogeneous sphere, whereas the optically thin index is fitted in the high-frequency tail and must satisfy $\tau=1$4 (Fromm et al., 2013).

These observables can be converted into intrinsic parameters. From the core-shift normalization $\tau=1$5 and $\tau=1$6, and assuming a lower cutoff Lorentz factor $\tau=1$7, the magnetic field and particle density follow

$\tau=1$8

In CTA 102 the $\tau=1$9 GHz core sits at $R$0 pc from the jet base. Assuming equipartition, $R$1, and $R$2 gives

$R$3

(Fromm et al., 2013). A related analysis reported $R$4 G and $R$5–$R$6, depending on model assumptions (Fromm et al., 2013).

Along the jet, pixel-by-pixel spectral fits in CTA 102 yielded $R$7 decreasing from $R$8 mG at $R$9 mas down to $\tau=1$0 mG at $\tau=1$1 mas, and $\tau=1$2 dropping from $\tau=1$3 to $\tau=1$4 (Fromm et al., 2013). The magnetization $\tau=1$5 calculated downstream falls from $\tau=1$6 at $\tau=1$7 mas to $\tau=1$8 at $\tau=1$9 mas, with local bumps up to $\nu$0 at recollimation sites. Downstream, local peaks in turnover parameters occur at $\nu$1 mas, $\nu$2 mas, $\nu$3 mas, and $\nu$4 mas (Fromm et al., 2013).

The physical interpretation advanced for the 2006 flare is a shock–shock interaction in an over-pressured jet. As a traveling shock passes the core region, both $\nu$5 and $\nu$6 rise and then decline as the disturbance propagates downstream; the core shift temporarily flattens because particle energy density exceeds magnetic energy density (Fromm et al., 2013). The same line of work also states that the source kinematics together with the spectral and structural variations can be described by helical motions in an over-pressured jet (Fromm et al., 2013).

4. Opacity, viewing geometry, and spectral flattening

The frequency dependence of the spectral core can also be derived from explicit synchrotron radiative-transfer models. For a perfect conical jet with half-opening angle $\nu$7, viewing angle $\nu$8, magnetic field $\nu$9, and electron distribution $r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$0 with $r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$1, the synchrotron emissivity and absorption take the form

$r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$2

$r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$3

with $r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$4 for a purely longitudinal field and $r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$5 for a purely toroidal field (Sharma et al., 2022).

The optical depth to a point $r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$6 is written

$r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$7

and the core is defined by $r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$8 (Sharma et al., 2022). In the finite-opening-angle geometry,

$r_{\rm core}(\nu)\propto \nu^{-1/k_r}.$9

For moderate $k_r=1$0, this reduces to the familiar power law $k_r=1$1, with the equivalent expression

$k_r=1$2

This is the same functional dependence used in observational core-shift studies (Sharma et al., 2022).

The same model links core displacement to spectral flattening. In the classical Blandford–Konigl picture,

$k_r=1$3

with

$k_r=1$4

The paper concludes that a smaller jet inclination angle or a higher electron density causes the jet core position to move downstream of the jet and that this displacement gives rise to spectral flattening (Sharma et al., 2022).

Parameter dependencies are explicit. As $k_r=1$5 decreases, the path length through the optically thick base increases, the $k_r=1$6 surface moves farther out, the core shift becomes larger, and $k_r=1$7 tends to be flatter or inverted. Larger $k_r=1$8 likewise produces stronger self-absorption and a stronger core shift. Parallel and perpendicular magnetic-field orientations produce different $k_r=1$9 values and therefore different core-shift slopes (Sharma et al., 2022). In the CTA 102 interpretation, the measured $B(r)\propto r^{-b}$0 during the flare is consistent with a temporarily particle-dominated state (Fromm et al., 2013).

5. Spectral core in network science

In network science, the spectral core is defined through a relaxed continuous optimization problem for core–periphery detection. For an undirected network with adjacency matrix $B(r)\propto r^{-b}$1, a natural discrete objective is

$B(r)\propto r^{-b}$2

where $B(r)\propto r^{-b}$3 is a ranking permutation. This objective is equivalent to the maximum-likelihood reordering under the logistic core–periphery random-graph model with

$B(r)\propto r^{-b}$4

(Tudisco et al., 2018).

The nonlinear spectral method replaces the discontinuous kernel $B(r)\propto r^{-b}$5 by the $B(r)\propto r^{-b}$6 surrogate

$B(r)\propto r^{-b}$7

and defines

$B(r)\propto r^{-b}$8

on the positive $B(r)\propto r^{-b}$9-sphere

$N(r)\propto r^{-n}$0

The relaxed logistic core–periphery problem is

$N(r)\propto r^{-n}$1

or equivalently $N(r)\propto r^{-n}$2 by $N(r)\propto r^{-n}$3-homogeneity (Tudisco et al., 2018).

Its gradient map is

$N(r)\propto r^{-n}$4

Karush–Kuhn–Tucker conditions yield the nonlinear eigenproblem

$N(r)\propto r^{-n}$5

With $N(r)\propto r^{-n}$6, defining

$N(r)\propto r^{-n}$7

one obtains a fixed-point characterization: $N(r)\propto r^{-n}$8 solves the relaxed problem if and only if $N(r)\propto r^{-n}$9 is a fixed point of $\tau=1$00 (Tudisco et al., 2018).

The algorithm iterates

$\tau=1$01

starting from $\tau=1$02, and outputs the positive score vector $\tau=1$03; nodes are then reordered by descending normalized score $\tau=1$04 (Tudisco et al., 2018). The theory is unusually strong. For any $\tau=1$05 and $\tau=1$06, the relaxed problem has a unique positive solution $\tau=1$07, the map $\tau=1$08 is a strict contraction in the Thompson metric

$\tau=1$09

with ratio $\tau=1$10, and the iteration converges globally at the geometric rate

$\tau=1$11

The fixed point $\tau=1$12 is the continuous coreness score, and nodes with large $\tau=1$13 form the network’s spectral core (Tudisco et al., 2018).

Computationally, if $\tau=1$14 is the number of nonzero edges, forming $\tau=1$15 costs $\tau=1$16 and the normalization step costs $\tau=1$17; on sparse networks with $\tau=1$18 each iteration is $\tau=1$19 (Tudisco et al., 2018). Empirically, on stochastic block-model benchmarks all methods performed similarly in one regime, whereas in another regime the nonlinear spectral method and Degree outperformed Sim-Ann; Sim-Ann was approximately $\tau=1$20 slower than the nonlinear spectral method. On the logistic core–periphery random model, the ordering quality satisfied NSM $\tau=1$21 Degree $\tau=1$22 Sim-Ann in maximizing likelihood. On real networks, reordered adjacency matrices showed sharper two-block patterns for NSM, and NSM produced the largest $\tau=1$23 in every case examined (Tudisco et al., 2018).

6. Multiplex and alternative spectral formulations

The single-layer nonlinear spectral core-periphery framework has been extended to multiplex networks. For a node-aligned multiplex network with nonnegative adjacency tensor

$\tau=1$24

the multiplex nonlinear objective is

$\tau=1$25

where $\tau=1$26 is the node-coreness vector and $\tau=1$27 is the layer-coreness vector (Bergermann et al., 2023). Constraining $\tau=1$28 and $\tau=1$29 yields a nonconvex homogeneous optimization problem with an alternating fixed-point iteration based on $\tau=1$30 and $\tau=1$31.

The convergence result is again explicit. If

$\tau=1$32

the iteration converges from any positive start to the unique global maximizer $\tau=1$33 at linear rate $\tau=1$34; for any $\tau=1$35, even if this condition fails, the algorithm monotonically increases the objective and converges to a stationary point (Bergermann et al., 2023). Per iteration, the work is $\tau=1$36 plus normalization costs $\tau=1$37 and $\tau=1$38, so runtime is linear in the total edge count. The method also introduces a multiplex QUBO score

$\tau=1$39

used to select the optimal binary core size (Bergermann et al., 2023).

Synthetic “one informative + one noise” experiments showed that the method can assign near-zero weight to a noise layer and recover the original core robustly; on seven real multiplex datasets, MP NSM with $\tau=1$40 achieved the highest multiplex QUBO scores in nearly all cases and produced visibly sharper L-shapes in reordered adjacency plots (Bergermann et al., 2023).

A different spectral tradition uses the random-walk Laplacian rather than a nonlinear fixed point. For an undirected graph with adjacency matrix $\tau=1$41 and diagonal degree matrix $\tau=1$42, the random-walk Laplacian is

$\tau=1$43

and coreness is inferred from the eigenvector associated with the largest eigenvalue of $\tau=1$44, or equivalently the second-largest eigenvector of $\tau=1$45 (Cucuringu et al., 2014). Cucuringu et al. define a coreness score from this eigenvector and classify vertices by maximizing the CP-density objective

$\tau=1$46

This method is $\tau=1$47 for eigenvector computation on sparse graphs and is described as scalable to networks with millions of edges (Cucuringu et al., 2014).

These two strands use different operators and guarantees. The nonlinear method is built around a max-type surrogate and nonlinear Perron–Frobenius theory, whereas the Laplacian approach uses linear diffusion modes. A common misconception is to treat them as equivalent “spectral core” algorithms. The literature presents them as distinct methodologies for the same mesoscale target.

7. Related uses in AGN population studies and molecular spectroscopy

In radio-source population studies, the “spectral core” denotes the compact, relativistically boosted jet base whose flat spectral signature is revealed in high-$\tau=1$48 objects (Pei et al., 2019). The central quantity is the core-dominance parameter

$\tau=1$49

with rest-frame correction

$\tau=1$50

and, in the two-component beaming model,

$\tau=1$51

with $\tau=1$52 for a continuous jet or $\tau=1$53 for a moving sphere (Pei et al., 2019). The radio spectral index is defined by $\tau=1$54, and the composite spectrum obeys

$\tau=1$55

For the sample of $\tau=1$56 AGNs, fitting this relation yielded $\tau=1$57, $\tau=1$58, $\tau=1$59, and $\tau=1$60 for the full sample (Pei et al., 2019). The reported sequence in beaming strength is BL Lac $\tau=1$61 FSRQ $\tau=1$62 Seyfert $\tau=1$63 Galaxy $\tau=1$64 FR I&II, whereas the sequence in spectral steepness is FR I&II $\tau=1$65 Seyfert $\tau=1$66 Galaxy $\tau=1$67 FSRQ $\tau=1$68 BL Lac (Pei et al., 2019). In this usage, the spectral core is an orientation- and beaming-sensitive compact component rather than a core-shift surface or a graph-theoretic subset.

In molecular electronic-structure theory, “Spectral-Core” refers to a state-specific protocol for core-level spectroscopy. The method combines Square Gradient Minimization (SGM), which optimizes excited-state orbitals by minimizing

$\tau=1$69

with the Gauss–Newton step

$\tau=1$70

and a Restricted Open-Shell Kohn-Sham energy functional for singlet core excitations,

$\tau=1$71

(Hait et al., 2019). The method is reported to predict the K edge of C, N, O, and F to a root mean squared error of $\tau=1$72 eV, with SCAN/ROKS and $\tau=1$73B97X-V/ROKS giving K-edge RMSE values of approximately $\tau=1$74–$\tau=1$75 eV; by contrast, TDDFT with the same functionals underestimates core excitations by $\tau=1$76–$\tau=1$77 eV and typically has greater than $\tau=1$78 eV error (Hait et al., 2019). Each SGM iteration costs approximately $\tau=1$79 to $\tau=1$80 a ground-state DFT orbital-gradient evaluation, and overall scaling remains that of ground-state DFT aside from the prefactor (Hait et al., 2019).

These uses broaden the semantic range of the term. In AGN surveys, the spectral core is a phenomenological radio component diagnosed statistically through $\tau=1$81 and $\tau=1$82 (Pei et al., 2019). In molecular spectroscopy, Spectral-Core is a computational protocol for core-level excitations (Hait et al., 2019). This suggests that “spectral core” functions less as a universal object name than as a reusable label for central structure inferred from spectral data or spectral optimization.