Neural audio codecs based on Residual Vector Quantization (RVQ) have become the dominant discrete representation for token-based general audio generation, yet resynthesizing high-quality audio from coarse codec tokens remains an open problem and bounds the fidelity of every system that generates them. Prior work has framed resynthesis as a choice between discrete token prediction and continuous regression. We argue that this dichotomy is incomplete and introduce geometric iterative retrieval, a paradigm that uses the RVQ layer hierarchy itself as a natural iterative decomposition in continuous codebook space. Rather than classifying over discrete vocabularies or regressing to a single target vector, our method performs contrastive retrieval in the codebook's geometric space. We evaluate our method on codec restoration tasks across speech and music, and show improvements over both single-pass token prediction and one-step regression baselines.
Figures & tables
Single-step
Iterative
Discrete
(degenerate)
Token prediction
Continuous
Regression
Diffusion / Retrieval (ours)
Table 1: Codec resynthesis methods organized along two axes, the prediction space (discrete vs. continuous) and the refinement strategy (single-step vs. iterative). Prior work occupies three of the four cells. Our method shares the continuous-iterative cell with diffusion, but iterates along the codec’s own RVQ hierarchy rather than along an externally imposed noise schedule, giving each step a meaning grounded in the codec rather than in the sampler.
Figure 1: Combined view of geometric iterative retrieval. (a) During training, the model predicts layers 2,…,D from the masked input stack in a single forward pass, supervised by a CLIP-style contrastive loss per predicted layer. (b) At inference, layer k+1 is predicted from layers 1,…,k , quantized by nearest-neighbor lookup against the layer- k+1 codebook, and concatenated for the next step. (c) Architecture: per-position codebook embeddings are aggregated by self-attention into one hidden state, processed by the DeBERTa-v3 encoder, and projected by Wout into the codebook space. Orange = discrete tokens, blue = continuous latents.
Method
LSD ↓
SI-SDR ↑
FAD ↓
Naive
10.82±0.08
−1.88±0.24
0.99
CE
12.42±0.20
−6.93±0.30
2.46
OSR
11.03±0.07
−0.22±0.27
0.61
MSE
12.86±0.22
−7.83±0.30
3.08
Ours
10.61±0.07
−0.63±0.28
0.94
Table 2: Codec restoration results on the combined val split (1500 clips). Lower is better for LSD and FAD, higher is better for SI-SDR. Best per column in bold , second-best underlined . LSD and SI-SDR are reported as mean ± half-width of the 95 % normal CI over the 1500 clips. FAD is a single set-level number per cell and admits no per-sample CI. OSR is one-step, layer-wise regression under cosine-similarity loss. MSE is the same backbone trained with MSE against the residual codebook latent.
K
LSD ↓
SI-SDR ↑
FAD ↓
1 (input only)
10.81
−1.88
0.99
2
10.26 ( −5.1% )
−0.34 ( +82% )
0.43 ( −57% )
3
10.20 ( −5.6% )
−0.22 ( +88% )
0.43 ( −57% )
4
10.23 ( −5.4% )
−0.30 ( +84% )
0.56 ( −43% )
5
10.34 ( −4.3% )
−0.41 ( +78% )
0.69 ( −30% )
6
10.41 ( −3.7% )
−0.50 ( +73% )
0.77 ( −22% )
Table 3: Reconstruction quality as a function of the number of predicted residual layers K (layer 1 is the ground-truth input, layers 2,…,K are our predictions). All three objective metrics improve from K=1 through K=3 and then degrade monotonically, because each additional residual layer brings progressively less spectral signal while still injecting prediction noise of roughly the same magnitude. Beyond K=3 , the noise outweighs the signal. Parenthesized values report the percentage change relative to the K=1 input-only baseline. Best per column in bold .
Condition
Mean (95% CI)
Δ
W/L
Naive
23.8 [21.0, 26.7]
+31.4
19 / 148
CE
0 9.1 [7.5, 10.7]
+45.5
0 6 / 164
MSE
0 7.6 [5.7, 9.6]
+46.8
0 7 / 164
OSR
43.4 [39.7, 47.1]
+11.2
48 / 117
Ours (K=3)
49.3 [45.7, 52.8]
+5.3
64 / 0 99
Ours
54.0 [50.1, 57.9]
—
—
Table 4: Double-blind human evaluation (19 raters × 9 clips). Higher is better. Δ is the mean per-trial difference (Ours − baseline). W ins/ L osses count the trials in which the baseline scored above or below Ours (ties omitted). All five pairwise differences favor our method and are significant under a Wilcoxon signed-rank test (Naive, CE, MSE with ∣z∣>9.7 , p<10−21 ; OSR with ∣z∣=5.5 , p=3.4×10−8 ; Ours (K=3) with ∣z∣=3.9 , p=1.1×10−4 ).
Variant
LSD ↓
SI-SDR ↑
FAD ↓
Ours
10.61±0.07
−0.63±0.28
0.94
Full-residual
12.29±0.22
−4.22±0.32
1.14
No-contrastive
10.75±0.07
−1.17±0.31
0.99
Additive
11.74±0.09
−1.24±0.30
0.96
Table 5: Ablation results on the combined validation split. Each variant perturbs one component, the prediction target (per-layer ek vs. cumulative residual ∑j=k+1Dej ), the loss (contrastive retrieval vs. cosine regression to the codebook vector), or the layer-aggregation mechanism (self-attention vs. additive summation). Best per column in bold , second-best underlined . LSD and SI-SDR are reported as mean ± half-width of the 95 % normal CI over the 1500 clips. FAD is a single set-level number per cell.
Neural audio codecs with residual vector quantization (RVQ) normally treat all frequencies uniformly, so their codebooks become spectrally entangled. Truncating stages then removes an unpredictable mix of frequencies. Parallel band decomposition addresses this by splitting audio into independent bands, but fragments the latent space and loses cross-frequency coherence. We introduce HARP (Harmonic-Aware Residual Partitioning), a training strategy that partitions RVQ stages into frequency-ordered groups where each group refines its target band while the decoder retains access to all lower frequencies. Overtones are reconstructed in the context of their fundamentals, preserving coherence that parallel methods lose. HARP requires no architectural changes; it only modifies the training loss, leaving inference identical to standard RVQ. On speech, music, and general audio, HARP outperforms both standard RVQ and parallel decomposition. MUSHRA listening tests also show perceptual improvements.
Qiaoyu Yang, Lixing He, Binyue Deng +1
Georgia Institute of Technology, Atlanta, United States · The Chinese University of Hong Kong, Hong Kong, China · Tencent Music Entertainment, Shenzhen, China
Neural audio codecs are widely used as tokenizers for spoken language models, but they are optimized for waveform reconstruction rather than autoregressive prediction. This mismatch injects acoustically driven uncertainty into the discrete token space and increases language-model perplexity. We propose \ours, which augments codec training with language-model-facing objectives while keeping both codec and LLM architectures unchanged. \ours introduces (i) future token prediction with Medusa-style multi-step heads to encourage multi-step predictability, and (ii) semantic alignment that matches audio and text representations via a memory-bank contrastive loss. A differentiable Gumbel bridge enables end-to-end gradients from these objectives to the codec encoder. On SALMon speech coherence, token LMs trained on \ours reach 61.6% accuracy (+12.1 points over AUV) while reducing perplexity 35. On Codec-SUPERB-tiny, \ours improves speech Mel distance by 5.0% over AUV while simultaneously achieving the learnability gains, demonstrating that reconstruction fidelity and token predictability can be improved together.
Ho-Lam Chung, Yiming Chen, Hung-yi Lee
Graduate Institute of Communication Engineering, National Taiwan University · ASUS Intelligent Cloud Services · NTU Artificial Intelligence Center of Research Excellence (NTU AI-CoRE)
Neural audio codecs compress waveforms into compact discrete tokens that underpin speech language models, real-time communication, and large-scale audio storage. Almost every dominant design, including residual vector quantization, finite scalar quantization, and single-codebook variants, follows the VQ-VAE template by partitioning the encoder latent through a learned codebook or a fixed scalar grid. We ask whether this partition is necessary. We introduce GS-Codec, a neural speech codec whose bottleneck is a parametric signal decomposition rather than a quantizer. We adapt Gaussian splatting from 3D scene reconstruction to one-dimensional latents. An inner optimization loop fits each encoder segment as a weighted sum of 1D Gaussian primitives. The decoder then reconstructs the waveform from the rendered sum. To avoid the cost of this iterative inner loop at inference time, we additionally train a lightweight GS Predictor Net that regresses the primitive parameters in a single forward pass. The encoder and decoder are trained end-to-end through the inner loop, with no quantizer anywhere in the training pipeline: the bottleneck is the decomposition itself, and scalar quantization is applied only post-training to the fitted parameters. Rather than relying on discrete codebook stages for bitrate control, our representation exposes a fine-grained rate-quality tradeoff: a single trained checkpoint supports post-training bitrate control by varying the number of primitives and the per-parameter bit depth, with no retraining required. GS-Codec matches or exceeds well-established open-source codecs such as EnCodec and DAC on speaker similarity (SIM), intelligibility (STOI), and perceptual quality (UTMOS) at comparable bitrates, while achieving comparable semantic performance (WER). Code and audio samples are available at https://ronaluf.github.io/gs-codec/
Ron Aluf, Alon Canfi, Eliya Nachmani
School of Electrical and Computer Engineering Ben-Gurion University of the Negev