Skip to content
AI.info

Research

Hellinger Multimodal Variational Autoencoders

Hellinger Multimodal Variational Autoencoders Overview Research area: Multimodal generative modeling — specifically multimodal variational autoencoders (VAEs) for weakly supervised learning across mul

arXiv
2601.06572
Published
2026-01-10
Authors
Huyen Vo, Isabel Valera

AI summary

Hellinger Multimodal Variational Autoencoders

Overview

Research area: Multimodal generative modeling — specifically multimodal variational autoencoders (VAEs) for weakly supervised learning across multiple data types.

Technical level: Intermediate. The paper's core derivation draws on information geometry (α-divergences, Hölder pooling, moment matching), but its argument and empirical results can be followed with a working knowledge of VAEs and Gaussian distributions.

Scope: The paper reinterprets how multimodal VAEs combine per-modality inference distributions, derives a new aggregation rule from a symmetric divergence, and evaluates the resulting model (HELVAE) against existing aggregation-based methods on three benchmarks.

What This Paper Is About

Multimodal VAEs must approximate a single joint posterior over a shared latent variable from several unimodal inference distributions ("experts"). The dominant strategies — product of experts (PoE) and mixture of experts (MoE) — each have a structural flaw: PoE can collapse or be biased toward overly sharp experts, while MoE cannot produce a posterior sharper than its components and typically requires sub-sampling of modality subsets during training, which imposes an upper bound on the multimodal ELBO and hurts generative quality. The paper reframes this aggregation step as probabilistic opinion pooling, selects the unique symmetric member of the α-divergence family (α = 0.5, proportional to the squared Hellinger distance), and derives a closed-form moment-matching approximation for Gaussian experts that removes the need for sub-sampling.

Key Contributions

  1. An opinion-pooling view of multimodal VAEs. The authors show that PoE and MoE are special cases of Hölder pooling — PoE corresponds to log-linear pooling in the limit α → 0 and MoE to linear pooling at α = 1 — placing existing aggregation schemes within a single parameterized family derived from minimizing α-divergence.

  2. A moment-matching approximation of Hölder pooling at α = 0.5 ("Hellinger aggregation"). Because α = 0.5 is the unique symmetric case of the α-divergence family (and the Hellinger distance is symmetric and bounded between 0 and 1), the pooled density can be projected onto a single diagonal Gaussian via matched first and second moments. The derivation yields explicit formulas for the pooled mean and variance in terms of pairwise auxiliary parameters and Bhattacharyya coefficients, with uniform weights across modalities as in PoE and MoE.

  3. HELVAE, a multimodal VAE that avoids sub-sampling. Hellinger aggregation defines the model, which trains on all modalities at once rather than random subsets. The authors report that it learns more expressive latent representations as additional modalities are observed and achieves better trade-offs between generative coherence and quality than existing methods.

  4. MoHELVAE, a mixture variant. Following MoPoE, the authors construct a mixture over all 2^M − 1 non-empty modality subsets, where each subset posterior is a Hellinger aggregation ("mixture of HELVAEs").

Main Findings

  • A structural advantage on controlled Gaussian experts. In an illustration with three Gaussian experts — expert 1 with mean (0, 0) and variance (0.5, 0.5), expert 2 with mean (1.0, 0.2) and variance (0.6, 0.6), and a sharper, distant expert 3 with mean (4, 0) and variance (0.2, 0.2) — PoE shifts toward the sharpest expert, MoE spreads mass across all experts and misses the local agreement between experts 1 and 2, while Hölder pooling at α = 0.5 places more mass around the first two experts and yields a sharper distribution with reduced variance relative to MoE. The cross terms of the form sqrt(q_i(z) q_j(z)) are described as inducing soft dependencies between experts by amplifying regions of mutual support.

  • Better behavior as agreeing experts accumulate. In a synthetic setup where "good" experts are N((0,0), 0.5I), "bad" experts are N((4,4), 0.2I), and the true distribution is N(0, I), PoE shows high negative log-likelihood and a low Bhattacharyya coefficient because it is biased toward the sharper bad expert. Hölder (α = 0.5) and Hellinger yield lower NLL, sharper distributions, and stronger overlap with the true distribution than MoE as the number of good experts increases.

  • PolyMNIST trade-offs. Across β ∈ {1, 2.5, 5, 10}, MVAE gives tight log-likelihood estimates but poor coherence; mixture-based models (MMVAE, MoPoE, MWBVAE, CoDEVAE) reach relatively high coherence but not the best overall; MMVAE+ and HELVAE both achieve better trade-offs, with HELVAE attaining higher coherence in both conditional and unconditional settings while MMVAE+ provides a tighter log-likelihood bound (attributed to its modality-specific variables).

  • Bimodal CelebA (Table 1). HELVAE achieves the highest conditional coherence (0.386 ± 0.003, rank 1) with second-best conditional FID (80.17 ± 1.52, rank 2) and second-best unconditional FID (81.89 ± 1.56, rank 2). MMVAE+ attains the highest unconditional coherence (0.258 ± 0.024) but high FID (101.85 ± 5.13). MVAE has the best FID on both settings (75.62 ± 3.01 conditional, 78.64 ± 2.73 unconditional) at a large cost in coherence (0.315 ± 0.013 conditional, 0.199 ± 0.021 unconditional).

  • CUB Image-Captions (Table 2). For caption-to-image generation, HELVAE ties CoDEVAE for highest conditional coherence (both 0.750; HELVAE ± 0.105, CoDEVAE ± 0.050) while achieving lower conditional FID than other models apart from MVAE: HELVAE 157.56 ± 0.95 (rank 2) versus CoDEVAE 175.97 ± 0.30 (rank 4), MMVAE+ 164.94 ± 1.50 (rank 3), and MVAE 148.69 ± 1.97 (rank 1) at coherence 0.233 ± 0.038 (rank 7).

  • Latent representation quality and synergy (Table 3). On CelebA, HELVAE ties CoDEVAE for the highest image accuracy (0.365), ranks fifth on text (0.400), and reaches the second-best joint score (0.408), improving over its own single-modality performance. MWBVAE has the best joint score (0.410); MMVAE+ has the best text score (0.461) but the weakest joint (0.346). The authors state that MWBVAE, CoDEVAE, and MMVAE+ do not benefit from additional modalities in the same way.

  • Multimodal synergy on PolyMNIST (Figure 6). At β = 2.5, MoPoE and CoDEVAE show decreased generative quality as more modalities are observed, and MMVAE+ shows limited benefit; HELVAE's FID steadily decreases as additional modalities are provided.

  • Robustness to β (Figure 7). Across β ∈ {3e-4, 3e-3, 3e-2, 3e-1, 1, 2.5, 3, 5, 10}, MMVAE+ reaches the highest log-likelihood but its coherence drops and FID rises sharply at small and large β. HELVAE achieves the highest coherence, competitive FID, and high log-likelihood across a broad range of β. MMVAE, MoPoE, MWBVAE, and CoDEVAE remain stable but with lower coherence and higher FID than HELVAE.

  • MoHELVAE improves coherence but slightly worsens FID (Table 4). On CelebA, MoHELVAE beats HELVAE on conditional coherence (0.394 ± 0.011 vs 0.386 ± 0.003), unconditional coherence (0.229 ± 0.020 vs 0.204 ± 0.013), and latent representation accuracy on image (0.379 ± 0.013 vs 0.365 ± 0.004), text (0.441 ± 0.019 vs 0.400 ± 0.004), and image-plus-text (0.430 ± 0.018 vs 0.408 ± 0.005). Its FID is slightly higher (82.28 ± 4.72 conditional, 83.14 ± 4.87 unconditional). Representation accuracy drops when combining image and text relative to text alone, a pattern also observed in other mixture-based models on this dataset.

  • Computational cost (Table 5). HELVAE has the shortest average batch training time on all three benchmarks: 0.0895 s on PolyMNIST (M = 5), 0.0863 s on CUB (M = 2), and 0.5272 s on CelebA (M = 2). For comparison, MVAE is 0.1887 / 0.2184 / 1.0294 s, MMVAE+ 0.2884 / 0.2270 / 0.9916 s, and MoHELVAE 0.5918 s on CelebA. The paper notes that CoDEVAE, MoPoE, and MWBVAE incur a computational cost of O(2^M − 1) in the number of modalities M.

Methodology in Plain English

Each modality gets its own encoder that outputs a Gaussian distribution over a shared latent variable. The question is how to merge those per-modality Gaussians into one. The authors treat the merge as an optimization problem: pick the pooled density that minimizes a weighted average of a discrepancy measure to each expert. Choosing the α-divergence family as that measure yields Hölder pooling; the standard PoE and MoE fall out as the limiting cases α → 0 and α = 1.

The interesting choice is α = 0.5, the only member of the family that is symmetric. At this value the pooling rule becomes a squared weighted sum of the square roots of the experts' densities (Equation 2), which expands into individual expert terms plus cross terms between every pair of experts — those cross terms are what let the aggregation reward agreement between experts instead of treating them as independent. To keep the result usable inside a VAE, the authors project this pooled density back onto a diagonal Gaussian by matching its first and second moments. Working through the algebra for diagonal-covariance Gaussian experts produces pairwise "auxiliary" means and variances plus Bhattacharyya coefficients that measure pairwise overlap; these give closed-form expressions for the pooled mean and variance. Weights across modalities are set uniformly, following PoE and MoE.

Experiments use three benchmarks: PolyMNIST (a synthetic dataset with five modalities, each sample made of MNIST digits of the same label patched onto random crops from five distinct background images), CUB Image-Captions (bird images paired with textual descriptions), and bimodal CelebA (face images paired with text describing facial attributes). Models are compared with MVAE, MMVAE, MoPoE, MMVAE+, MWBVAE, and CoDEVAE. Evaluation covers generative coherence, generative quality via FID, log-likelihood estimation as a lower bound on posterior approximation tightness, and a linear classifier on the learned latent representation. All results are averaged over three seeds, with each model reported at its best β ∈ {1, 2.5, 5} for CelebA and CUB (and β ∈ {1, 2.5, 5, 10} in the PolyMNIST trade-off analysis).

Why This Matters

The paper's central insight is that the two dominant multimodal VAE objectives are two ends of one spectrum, and that the symmetric midpoint of that spectrum has properties neither endpoint has. That reframing matters because it gives a principled route to aggregation rules that both sharpen the joint posterior and preserve the experts — a tension the authors note was previously documented as a hard trade-off between generative coherence and quality.

Real-world applications:

  • Cross-modal retrieval and generation with partially annotated data, such as generating images from captions or captions from images when only weak supervision is available.
  • Medical or scientific datasets where different measurements (imaging, text reports, sensor streams) describe the same subject but are rarely all available for every subject.
  • Robust inference under missing modalities, since a model that trains on all modalities at once can still condition on any subset at test time.
  • Representation learning for downstream classifiers, where shared latent spaces that improve when more modalities are observed are more useful for tasks like attribute prediction.

Industry relevance: The reported batch training times are shorter than the compared models on all three benchmarks, and HELVAE avoids the O(2^M − 1) subset enumeration that the paper attributes to MoPoE, MWBVAE, and CoDEVAE. The authors also highlight that HELVAE achieves results comparable to MMVAE+ without requiring modality-specific latent variables, which is relevant for teams that want strong multimodal generation with a simpler model.

Future Directions

  • Negative pooling weights. The authors observe that the pooled density at α = 0.5 takes the form of a squared mixture model, related to recent squared subtractive mixture formulations that showed increased expressiveness for density estimation. They suggest that exploring negative pooling weights is a promising direction.
  • Beyond Gaussian, diagonal-covariance experts. The derivation assumes multimodal Gaussians with diagonal covariance, and the pooled covariance is simplified by discarding off-diagonal terms; relaxing either assumption would extend applicability.
  • Learning rather than fixing the subset weights. The full derivation is given under subset weights, but the authors set them uniformly as in PoE and MoE because determining them is nontrivial in multimodal VAEs. How best to set these weights remains open.
  • Bridging the coherence/quality gap that remains. HELVAE does not dominate every metric: MVAE retains the best FID, MMVAE+ retains the tightest log-likelihood bound (attributed to its modality-specific variables), and the authors note that HELVAE images are less diverse because of the lack of a modality-specific latent space. Closing these gaps while keeping the sub-sampling-free training scheme is a natural next step.

Target Audience

Researchers and practitioners working on multimodal generative models, variational inference, or weakly supervised learning who want a principled alternative to product-of-experts and mixture-of-experts aggregation. It is also relevant to readers interested in probabilistic opinion pooling and α-divergences as tools for machine learning, and to engineers who need a multimodal VAE that trains on all modalities simultaneously with low per-batch cost.

Authors’ abstract

Multimodal variational autoencoders (VAEs) are widely used for weakly supervised generative learning with multiple modalities. Predominant methods aggregate unimodal inference distributions using either a product of experts (PoE), a mixture of experts (MoE), or their combinations to approximate the joint posterior. In this work, we revisit multimodal inference through the lens of probabilistic opinion pooling, an optimization-based approach. We start from Hölder pooling with $α=0.5$, which corresponds to the unique symmetric member of the $α\text{-divergence}$ family, and derive a moment-matching approximation, termed Hellinger. We then leverage such an approximation to propose HELVAE, a multimodal VAE that avoids sub-sampling, yielding an efficient yet effective model that: (i) learns more expressive latent representations as additional modalities are observed; and (ii) empirically achieves better trade-offs between generative coherence and quality, outperforming state-of-the-art multimodal VAE models.

Read the original paper