Proca Condensate: Concepts & Applications

• Proca condensate is a versatile term describing various constructions of massive vector fields with a physical longitudinal polarization across different physical settings.
• It encompasses distinct interpretations such as renormalized local vacuum bilinears, emergent effective fields, and self-gravitating Bose–Einstein condensates.
• Applications range from probing Casimir effects and holographic orders to understanding black-hole vacua and compact Proca stars.

“Proca condensate” is not a single universal object in current literature. The phrase is used for several distinct constructions built from a massive vector field $A_\mu$ obeying the Proca equation

$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$

together with the Lorenz-type constraint $\nabla_\mu A^\mu=0$, which is a consequence of the massive field equations rather than a gauge choice. Depending on context, the term denotes a renormalized local vacuum bilinear, an electric condensed phase whose infrared excitation is a Proca field, a self-gravitating Bose–Einstein condensate of complex massive vector bosons, a holographic or cosmological vector background, or an emergent or effective massive vector sector in analogue or non-Abelian settings (Saharian et al., 9 Jul 2025, Guimaraes et al., 2012, Herdeiro et al., 2020, Arias et al., 2016, Pastor-Marcos et al., 1 Apr 2026).

1. Terminological scope and common structure

Across the literature, the same phrase labels physically different quantities and phases.

Usage	Condensate object	Representative reference
Casimir vacuum density	$\langle F_{\mu\nu}F^{\mu\nu}\rangle$	(Saharian et al., 9 Jul 2025)
Curved-spacetime vacuum polarization	$\langle A_\mu A^\mu\rangle_\Psi$	(Prakash et al., 26 May 2026)
Electric condensed phase	Proca theory as low-energy effective theory of an electric condensate	(Guimaraes et al., 2012)
Self-gravitating vector soliton	Macroscopic Bose–Einstein condensate of complex massive vector bosons	(Herdeiro et al., 2020, García et al., 2016)
Holographic order parameter	Nonzero normalizable mode of a charged bulk Proca field	(Arias et al., 2016)
Cosmological vector condensate	Background expectation value $X_0=\langle X\rangle$, $X=\frac{Z_A}{2}A_\mu A^\mu$	(Pastor-Marcos et al., 1 Apr 2026)
Emergent analogue field	Spin-nematic mode cast into a Proca field on an acoustic spacetime	(Brunner et al., 16 Jun 2025)
Non-Abelian effective medium	Classical Proca sector sourced by a gluon condensate	(Dzhunushaliev et al., 2024)

The common kinematic element is the massive spin-1 field with a physical longitudinal polarization. This is decisive in several of the cited applications. In the Casimir problem, the longitudinal mode is constrained by PMC plates but not by PEC plates; in Schwarzschild quantization, it survives as a genuine even-parity branch with no Maxwell analogue; in generalized Proca theories, the entire interaction program is organized so that only the three desired polarizations propagate (Saharian et al., 9 Jul 2025, Prakash et al., 26 May 2026, Heisenberg, 2017).

This multiplicity of meanings is the first point of orientation. In some subfields, “condensate” means a local quadratic observable. In others, it means a many-body or mean-field phase. In yet others, it means a background value of the invariant $A_\mu A^\mu$ or an emergent effective field. A frequent misconception is therefore to assume that every “Proca condensate” is a symmetry-breaking order parameter. That is explicitly false in the Casimir and black-hole-vacuum settings (Saharian et al., 9 Jul 2025, Prakash et al., 26 May 2026).

2. Local vacuum condensates and vacuum polarization

In the Casimir problem for a Proca field between parallel plates, the condensate is defined precisely as the renormalized local invariant

$$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$

Here $\langle\cdots\rangle_0$ is the boundary-free Minkowski contribution, subtracted before the coincidence limit. In this usage, the “Proca condensate” is not a symmetry-breaking order parameter but a local vacuum expectation value of a quadratic field operator in a Casimir background, and the paper explicitly identifies $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$0 as “the condensate,” analogous to the gluon condensate in QCD (Saharian et al., 9 Jul 2025).

The geometry is $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$1-dimensional Minkowski spacetime with plates at $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$2 and $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$3. The analysis gives explicit two-point functions for the vector potential and field tensor, as well as the induced vacuum expectation values of $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$4, $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$5, $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$6, and $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$7. For PMC conditions, the condensate is positive; for PEC conditions, it is negative for $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$8 and vanishes for $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$9. After Minkowski subtraction it is finite for interior points $\nabla_\mu A^\mu=0$0, but it diverges near a plate with the standard surface behavior. In $\nabla_\mu A^\mu=0$1, the interpretation

$\nabla_\mu A^\mu=0$2

shows that the condensate measures the imbalance between magnetic and electric vacuum fluctuations induced by the boundaries (Saharian et al., 9 Jul 2025).

A structurally important result is that all observables determined purely by $\nabla_\mu A^\mu=0$3, including $\nabla_\mu A^\mu=0$4, have a smooth massless limit. By contrast, the zero-mass limit of $\nabla_\mu A^\mu=0$5 is subtle for PMC conditions because the energy-momentum tensor contains $\nabla_\mu A^\mu=0$6, which retains sensitivity to the longitudinal polarization mode even as $\nabla_\mu A^\mu=0$7. This sharply separates the field-strength condensate from vector-potential bilinears (Saharian et al., 9 Jul 2025).

A different local usage appears in canonical quantization of the Proca field on Schwarzschild spacetime. There the condensate is defined as

$\nabla_\mu A^\mu=0$8

This is explicitly a scalar contraction of the vector-potential two-point function, not a field-strength invariant. Because the Proca field is massive and the longitudinal mode is physical, the object is not affected by Maxwell-type gauge ambiguity in the same way as $\nabla_\mu A^\mu=0$9 for a massless gauge field. The paper does not carry out a full Hadamard renormalization, but evaluates finite differences such as $\langle F_{\mu\nu}F^{\mu\nu}\rangle$0 by subtracting the Boulware correlator (Prakash et al., 26 May 2026).

In that Schwarzschild setting, the condensate is strongly vacuum-dependent. The Boulware state is singular on the future horizon; the Unruh state is regular on $\langle F_{\mu\nu}F^{\mu\nu}\rangle$1 and naturally tied to outgoing Hawking flux; the Hartle–Hawking state is thermal in both in and up sectors. The numerics show that the condensate becomes significant near the boundary of the future horizon. A notable polarization effect is that the monopole mode does not contribute to the horizon condensate, even though the longitudinal Proca sector has no Maxwell analogue and contributes to the asymptotic Hawking spectrum (Prakash et al., 26 May 2026).

Taken together, these two lines of work show that “Proca condensate” can mean a renormalized local probe of vacuum polarization, either as $\langle F_{\mu\nu}F^{\mu\nu}\rangle$2 in a boundary-value problem or as $\langle F_{\mu\nu}F^{\mu\nu}\rangle$3 in a black-hole vacuum-state problem. The choice of bilinear is therefore not universal.

3. Condensed phases, order parameters, and effective massive-vector descriptions

A very different meaning arises in the generalized Julia–Toulouse approach to electric condensation. In that framework, Proca theory is not introduced as an explicit mass deformation of electrodynamics; it is derived as the low-energy effective theory of an electric condensate. Starting from Maxwell theory with diluted electric charges and external monopoles, one passes to the dual description and replaces the regular non-minimal combination by a collective condensate field,

$\langle F_{\mu\nu}F^{\mu\nu}\rangle$4

After dualizing back, the infrared theory becomes a Proca theory in which the massive vector mode is the propagating excitation of the electric condensed phase (Guimaraes et al., 2012).

In this usage, the mass parameter $\langle F_{\mu\nu}F^{\mu\nu}\rangle$5 is the scale of the condensate and produces a Meissner effect. Magnetic flux is then confined into vortices, and in the presence of monopoles the physical object is not a Dirac string but the Dirac brane invariant

$\langle F_{\mu\nu}F^{\mu\nu}\rangle$6

Integrating out the vector field yields an effective monopole–antimonopole potential with both a Yukawa piece and a linear confining term,

$\langle F_{\mu\nu}F^{\mu\nu}\rangle$7

Here the “Proca condensate” is a phase of condensed electric charges whose long-distance excitations are described by a massive vector field (Guimaraes et al., 2012).

In holography, the phrase refers to spontaneous condensation of a charged massive vector field in AdS. The model contains a Maxwell field $\langle F_{\mu\nu}F^{\mu\nu}\rangle$8 and a complex Proca field $\langle F_{\mu\nu}F^{\mu\nu}\rangle$9 with

$\langle A_\mu A^\mu\rangle_\Psi$0

on an AdS-Schwarzschild background. The full ansatz

$\langle A_\mu A^\mu\rangle_\Psi$1

supports an $\langle A_\mu A^\mu\rangle_\Psi$2-wave phase with $\langle A_\mu A^\mu\rangle_\Psi$3, a $\langle A_\mu A^\mu\rangle_\Psi$4-wave phase with $\langle A_\mu A^\mu\rangle_\Psi$5, and an $\langle A_\mu A^\mu\rangle_\Psi$6-wave coexistence phase with all fields nontrivial (Arias et al., 2016).

The condensates are read from the normalizable coefficients of the near-boundary expansions of $\langle A_\mu A^\mu\rangle_\Psi$7 and $\langle A_\mu A^\mu\rangle_\Psi$8, with source-free conditions $\langle A_\mu A^\mu\rangle_\Psi$9 and $X_0=\langle X\rangle$0. For $X_0=\langle X\rangle$1, the $X_0=\langle X\rangle$2-wave instability appears first, with $X_0=\langle X\rangle$3. For $X_0=\langle X\rangle$4, the $X_0=\langle X\rangle$5-wave phase appears first at $X_0=\langle X\rangle$6, and an $X_0=\langle X\rangle$7 branch emerges continuously at $X_0=\langle X\rangle$8. In the coexistence phase, a spontaneous current $X_0=\langle X\rangle$9 turns on because the mixed bulk couplings source $X=\frac{Z_A}{2}A_\mu A^\mu$0 once both $X=\frac{Z_A}{2}A_\mu A^\mu$1 and $X=\frac{Z_A}{2}A_\mu A^\mu$2 are nonzero. In this setting, both the $X=\frac{Z_A}{2}A_\mu A^\mu$3- and $X=\frac{Z_A}{2}A_\mu A^\mu$4-wave order parameters are Proca condensates, since they are different components of the same bulk Proca field (Arias et al., 2016).

A third usage arises in generalized Proca cosmology and asymptotic safety. There the condensate is encoded in the invariant

$X=\frac{Z_A}{2}A_\mu A^\mu$5

motivated by the fact that an isotropic FLRW background can arise naturally for a purely timelike condensate of the Proca field. In the functional-renormalization-group truncation, the effective action is expanded around $X=\frac{Z_A}{2}A_\mu A^\mu$6, and the condensate is related to fluctuation couplings by

$X=\frac{Z_A}{2}A_\mu A^\mu$7

The same work identifies several fixed points, including an interacting Proca$X=\frac{Z_A}{2}A_\mu A^\mu$8 point with four relevant directions in the strict Proca limit, providing evidence for non-perturbative renormalisability of vector-tensor theories. The paper is explicit, however, that $X=\frac{Z_A}{2}A_\mu A^\mu$9 at the fluctuation level and that kinetic terms for the condensate variable $A_\mu A^\mu$0 are neglected, so the condensate is only partially dynamical in the truncation (Pastor-Marcos et al., 1 Apr 2026).

These condensed-phase usages are constrained by the broader generalized-Proca consistency program. A single-field generalized Proca theory admits a finite family of derivative self-interactions $A_\mu A^\mu$1 through $A_\mu A^\mu$2, constructed so that only three polarizations propagate; in multifield settings, secondary-constraint conditions are essential, and many rotationally symmetric multi-Proca interactions suggested previously propagate ghosts (Heisenberg, 2017, Díez et al., 2019). A plausible implication is that not every proposed vector order parameter in cosmology or holography can be embedded into a healthy Proca theory without checking the constraint structure.

4. Self-gravitating and solitonic Proca condensates

In gravitational physics, “Proca condensate” most often means a self-gravitating Bose–Einstein condensate of complex massive vector bosons. Proca stars are everywhere regular, asymptotically flat self-gravitating solitons described by the Einstein–complex–Proca system. With the stationary ansatz

$A_\mu A^\mu$3

the explicit time dependence cancels out of the geometry because the stress-energy depends on bilinears such as $A_\mu A^\mu$4 and $A_\mu A^\mu$5. The family of solutions forms a spiral in the $A_\mu A^\mu$6-versus-$A_\mu A^\mu$7 plane, beginning at the Newtonian limit $A_\mu A^\mu$8, $A_\mu A^\mu$9; the branch from the maximum mass back to the Newtonian limit is perturbatively stable, while beyond the maximum mass the solutions are unstable (Herdeiro et al., 2020).

This self-gravitating usage extends beyond the minimally coupled, asymptotically flat case. Charged Proca stars, Proca Q-balls, Proca Q-stars, and their charged counterparts arise when one combines a complex vector field with a gauged $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$0 sector and, optionally, a self-interaction potential

$$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$1

In flat space, non-gravitating Proca Q-balls exist in a bounded frequency interval $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$2; with $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$3, charged solutions exist up to about $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$4. With gravity included and $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$5, the maximal masses and minimal radii reported are $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$6, $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$7 for PS; $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$8, $$\langle F_{\mu\nu}F^{\mu\nu}\rangle = g^{\mu\rho}g^{\nu\sigma} \lim_{x'\to x} \left[ \langle F_{\mu\nu}F'_{\rho\sigma}\rangle - \langle F_{\mu\nu}F'_{\rho\sigma}\rangle_0 \right].$$9 for CPS; $\langle\cdots\rangle_0$0, $\langle\cdots\rangle_0$1 for IPS; and $\langle\cdots\rangle_0$2, $\langle\cdots\rangle_0$3 for CIPS. The corresponding maximum compactnesses are $\langle\cdots\rangle_0$4, $\langle\cdots\rangle_0$5, $\langle\cdots\rangle_0$6, and $\langle\cdots\rangle_0$7 (García et al., 2016).

Self-interactions introduce an EFT subtlety that has become central in the Proca-star literature. For a quartic self-interaction

$\langle\cdots\rangle_0$8

the effective metric for longitudinal perturbations can become singular, with the radial component $\langle\cdots\rangle_0$9 vanishing. In the EFT, this is the point where the radial field equation becomes singular. The UV-complete analysis with an auxiliary heavy scalar $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$00 shows, however, that Proca-star solutions continue to exist beyond the EFT threshold: $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$01 signals breakdown of the low-energy description rather than a fundamental pathology, while the EFT ceases to be trustworthy before a ghost region associated with $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$02 is reached (Aoki et al., 2022).

The same condensate idea survives on nontrivial topology. In asymptotically AdS Ellis wormholes, a complex Proca field forms a wormhole-supported condensate with harmonic dependence

$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$03

These solutions are classified by reflection parity across the throat. In the symmetric AdS case, the asymptotic expansion has no $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$04 term, so the mass vanishes while the Noether charge remains finite. As the cosmological constant decreases, the familiar $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$05-$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$06 spiral gradually disappears, the condensate becomes more sharply localized near the throat, and the geometry can approach a black-bounce-like configuration, either at the throat or on both sides depending on the Proca parity class (Li et al., 25 Mar 2025).

This self-gravitating literature is the one in which “condensate” is closest to the standard many-body usage. Even there, however, the precise realization is relativistic and field-theoretic: the condensate is a coherent classical configuration of a complex massive vector field, supported by gravity, self-interactions, or background geometry.

5. Response, stability, mergers, and cloud dynamics

Once Proca condensates are viewed as compact or quasi-bound objects, the relevant questions become tidal response, nonlinear stability, merger phenomenology, and environmental sensitivity.

For spherical Proca stars, the nonlinear stability picture tracks the turning-point structure of the equilibrium branch. Fully nonlinear evolutions confirm that the separation between stable and unstable configurations occurs at the solution with maximal ADM mass. Depending on the sign of the binding energy and on the perturbation, unstable stars have three possible fates: migration to the stable branch, total dispersion, or collapse to a Schwarzschild black hole. In the collapse channel, a long-lived exterior Proca remnant—a “Proca wig” composed of quasi-bound states—can remain outside the horizon, with a lifetime that scales inversely with the Proca mass (Sanchis-Gual et al., 2017).

Binary dynamics sharpen the distinction between vector condensates and vacuum black holes. In head-on collisions of equal-mass Proca stars, low-compactness configurations can leave a stable Proca-star remnant, whereas more compact stars form a transient hypermassive Proca star that later decays into a black hole, often temporarily surrounded by Proca quasi-bound states. In orbital mergers, the most compact binaries produce a Kerr black hole with a transient Proca remnant, while less compact binaries can form a massive Proca star with angular momentum, though out of equilibrium. The waveform can show delayed collapse, a two-stage burst, and late-time distortions associated with the external Proca cloud, all absent in ordinary vacuum black-hole binaries (Sanchis-Gual et al., 2018).

The linear tidal response of spherically symmetric Proca stars provides another diagnostic. The electric-type quadrupolar Love number is positive and the magnetic-type one is negative. For the 28 background solutions studied, compactness ranges from $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$07 to $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$08, with $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$09 roughly $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$10 to $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$11 and $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$12 roughly $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$13 to $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$14. At fixed compactness, the electric and magnetic Love numbers are closer in magnitude than in scalar boson stars, with $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$15 but only slightly so (Herdeiro et al., 2020).

A related but distinct dynamical object is the Proca cloud produced by black-hole superradiance. On a fixed Kerr–de Sitter background, the cloud is a quasibound accumulation of a massive vector field rather than a fully backreacting self-gravitating star. The principal result is that parameter choices producing growth at $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$16 can become decaying states when $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$17. For example, at $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$18, modes with $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$19 and $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$20 have positive $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$21 at $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$22 but negative $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$23 at $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$24. This indicates that cosmological expansion can quench part of the superradiant instability window (Fell et al., 2024).

A common misconception is to identify all massive-vector condensates in gravity with the same object. The literature distinguishes at least three: self-gravitating Proca stars, post-collapse quasi-bound “wigs,” and black-hole superradiant clouds. They are connected, but not interchangeable.

6. Emergent, non-Abelian, and effective-medium realizations

The term also appears in analogue and effective-field settings where the Proca field is not fundamental.

In a spin-1 Bose–Einstein condensate, excitations around the polar phase can be reorganized into emergent relativistic fields on acoustic spacetimes. The spin-nematic rotation mode $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$25 obeys, in the hydrodynamic limit,

$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$26

After decomposing the planar mode into longitudinal and transverse parts, the transverse sector can be cast into a Proca equation on the spin-wave acoustic metric,

$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$27

This emergent Proca field is valid on length scales larger than the spin-healing length. By tuning the quadratic Zeeman coefficient $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$28, one can realize an FLRW-type metric with scale factor

$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$29

opening a pathway toward quantum simulation of cosmological particle production of Proca quanta via quenches or magnetic-field ramps (Brunner et al., 16 Jun 2025).

In non-Abelian Proca theory with external sources, the condensate language shifts again. The theory

$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$30

admits chromoelectric flux tubes. The paper then proposes an interpretation in which the almost-classical $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$31 components satisfy

$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$32

while the coset components satisfy

$$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$33

In that scenario, the purely quantum sector forms a gluon condensate that both sources the classical Proca sector and can generate an effective Proca mass through nonperturbative correlators. The resulting flux tube has a longitudinal chromoelectric field with nonlinear and gradient components, together with a transverse chromomagnetic field (Dzhunushaliev et al., 2024).

A related SU(3) Proca–Higgs model constructs cylindrically symmetric tubes carrying either longitudinal color electric flux or energy flux and momentum. The authors state that the existence of such tubes depends crucially on the presence of the Higgs field and that there are no such solutions without it. These topologically trivial tubes demonstrate the dual Meissner effect, in the sense that the electric field is pushed out by the Higgs scalar field (Dzhunushaliev et al., 2021).

The conceptual boundary is therefore sharp. In the spinor-BEC case, the Proca field is emergent. In the non-Abelian QCD-inspired case, the Proca field is an effective classical sector interacting with a gluon condensate. In neither case is the term restricted to a simple many-body occupation number.

A plausible synthesis is that the phrase “Proca condensate” is best understood as a family resemblance term rather than a single invariant concept. What unifies the usages is the centrality of a massive vector sector with a physical longitudinal mode. What differentiates them is the object being called a condensate: a local bilinear such as $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$34 or $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$35, an electric condensed phase whose infrared excitation is Proca-like, a self-gravitating Bose–Einstein condensate of complex vector bosons, a background value of $$\nabla_\nu F^{\mu\nu}-m^2A^\mu=0, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu,$$36, an emergent analogue vector field, or an effective classical field supported by a gluon condensate.