--- title: Dark Photon Solitons url: https://www.emergentmind.com/topics/dark-photon-solitons type: topic --- # Dark Photon Solitons Dark photon solitons are localized, self-gravitating, nonrelativistic configurations of a massive spin-1 field, usually described as a dark photon or Proca boson. In the dark-matter context they appear as coherent ground-state clumps of wave-like vector dark matter, and the literature refers to closely related objects as “dark photon stars,” “Proca stars,” and “vector solitons” [2203.10100][2301.11470][2406.05031]. Across these formulations, the common structure is a bound condensate whose leading dynamics are governed by a Schrödinger–Poisson system, with polarization, spin density, self-interactions, and nonminimal gravitational interactions controlling the detailed profile, stability, and phenomenology [2406.05031]. ## 1. Terminology and scope The term *dark photon soliton* denotes an astrophysical or cosmological soliton built from a massive vector field in the dark sector. In the minimal setup the field is a Proca boson with mass \(m\), and in the weak-field, nonrelativistic regime its bound states are described by three complex mode functions corresponding to the vector components [2203.10100]. When the emphasis is on compact-object phenomenology, the same solutions are often called “dark photon stars” or “Proca stars”; when the emphasis is on wave mechanics and EFT, the name “vector soliton” is common [2203.10100][2301.11470]. A useful terminological distinction is that this subject is not the same as the nonlinear-optics literature on “dark” photonic solitons. “Dark topological valley Hall edge solitons,” for example, are intensity-dip edge excitations in a photonic topological insulator formed at domain walls between honeycomb lattices with broken inversion symmetry, described by a \((2+1)\)-dimensional nonlinear Schrödinger equation in waveguide arrays [2206.14460]. This is a distinct usage of *dark*: there it refers to a dark soliton on a finite optical background, not to a dark-sector vector field. Within dark-matter theory, the broader unifying point is that solitons are a generic prediction for ultralight dark matter. A “unified view” has shown that real or complex, scalar or vector dark matter can exhibit universal nonrelativistic soliton properties such as conserved charges, mass–radius relations, stability criteria, and approximate profiles, while vector-specific effects arise from macroscopic spin density and polarization-dependent self-interactions [2406.05031]. ## 2. Field-theoretic description and nonrelativistic limit In a curved FRW spacetime, the minimal dark-photon setup starts from the Proca action \[ S = \int d^4x \sqrt{-g}\left[-\frac14 F_{\mu\nu}F^{\mu\nu} + \frac12 m^2 A_\mu A^\mu \right], \] with \(F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu\). Variation yields the Proca equation \(\nabla_\nu F^{\nu\mu}+m^2A^\mu=0\) together with the constraint \(\nabla_\mu A^\mu=0\) [2203.10100]. In the nonrelativistic, weak-field limit one decomposes the spatial components as \[ A_i(x,t)=\mathrm{Re}\!\left[\frac{\psi_i(x,t)}{\sqrt{2m^2a^3}}e^{-imt}\right], \] and introduces a Newtonian potential \(\Phi\) satisfying \(\nabla^2\Phi=4\pi G(\rho-\langle\rho\rangle)\) [2203.10100]. The resulting dynamics are a vector Schrödinger–Poisson system, \[ i\partial_t\psi_i=\left[-\frac{\nabla^2}{2m}+m\Phi\right]\psi_i,\qquad \nabla^2\Phi=4\pi G\,\frac{\sum_i|\psi_i|^2-\langle|\psi_i|^2\rangle}{a}, \] which is formally identical to three copies of a scalar Schrödinger–Poisson problem after eliminating \(A_0\), but with an intrinsic spin \[ S_p=\frac{i}{m}\epsilon_{pqr}\int \psi_q\psi_r^*\,d^3x \] that can take any value \(0\le |S|\le N\hbar\), where \(N\) is the particle number [2203.10100]. The ground-state soliton has zero orbital angular momentum but can carry macroscopic spin. A more general complex-vector theory augments the Proca action by a quartic self-interaction and nonminimal couplings to curvature, \[ -\frac{\lambda}{4}(X_\mu^*X^\mu)^2-\xi_1 R X_\mu^*X^\mu-\xi_2 R^{\mu\nu}X_\mu^*X_\nu. \] In the nonrelativistic limit this produces a vector Gross–Pitaevskii–Poisson system with a spin density \[ \mathcal S_i=i\varepsilon_{ijk}\psi_j\psi_k^* \] and a gradient-dependent term proportional to \(\xi\,\nabla^2\rho\), where \(\xi\equiv \xi_1+\tfrac12\xi_2\) [2406.05031]. The associated conserved quantities include particle number \(N=\int d^3x\,|\psi|^2\), total spin \(S_i=\int d^3x\,\mathcal S_i\), gravitational mass \(M=\int d^3x\,\rho_\xi\), and an energy functional containing a spin-spin contribution proportional to \(-\mathcal S^2\) [2406.05031]. This is the precise sense in which dark photon solitons generalize scalar Schrödinger–Poisson solitons. ## 3. Formation in the early universe and dark-matter substructure One concrete formation channel begins during inflation. If a massive dark photon is present while \(H_I\gg m\), the longitudinal mode behaves like a light scalar and acquires a nearly scale-invariant spectrum at horizon exit, \[ P_{A_L}(k)\big|_{a=k/H_I}=\left(\frac{kH_I}{2\pi m}\right)^2. \] After reheating and radiation domination, the perturbations redshift differently above and below the scale \(k_\star\) defined by \(H(a_\star)=m\) [2203.10100]. In parallel, quantum pressure enters the Euler equation through the Bohm potential and defines a comoving quantum Jeans scale \[ k_J(\rho)\equiv a\,(16\pi G\rho m^2)^{1/4}. \] A key result is the parametric coincidence \(k_J(\langle\rho\rangle)\approx 1.9\,k_\star\) at matter–radiation equality, so perturbations at \(k\sim k_\star\) first become nonlinear just as quantum pressure becomes subdominant [2203.10100]. This mechanism yields a rich hierarchy of substructure rather than a smooth dark-matter fluid. A substantial fraction of the dark matter collapses into gravitationally bound solitons, which are fully quantum coherent objects [2203.10100]. The characteristic masses are \[ 10^{-16} M_\odot \left(\frac{10^{-5}\,\mathrm{eV}}{m}\right)^{3/2}, \] and the central densities are typically a factor \(10^6\) larger than the local background dark matter density, with the detailed study reporting a range of \(10^2\)–\(10^8\) times the ambient density depending on the population considered [2203.10100]. Numerical simulations find formation with \(M\simeq c_M M_J(a)\), \(c_M\approx 0.45\), and today’s most abundant solitons near \(M\sim 10^{-16}M_\odot(10^{-5}\,\mathrm{eV}/m)^{3/2}\) [2203.10100]. The same simulations indicate that a comparable fraction of the energy density is initially stored in, and subsequently radiated from, long-lived quasi-normal modes, and that solitons are surrounded by characteristic “fuzzy” halos in which wave effects are enhanced relative to virialized dark matter expectations [2203.10100]. The soliton mass function \(df_s/d\log M\) peaks at \(M\sim 0.1\,M_J^{eq}\) and integrates to \(f_s\gtrsim 0.05\) of the total dark matter, while lower-density compact halos of mass \(\sim 10^3 M_J^{eq}\) are also produced at larger scales [2203.10100]. A plausible implication is that dark photon solitons should be treated not as rare exotica but as a potentially generic component of vector dark matter substructure whenever the production history populates the relevant wave modes. ## 4. Internal structure, polarization, and stability For isolated bound states, a general polarized vector-soliton ansatz may be written as \[ X_i(t,\mathbf x)=\frac12 e^{-i\omega t}\sum_{a=1}^3 c^{(a)}X(r)\,\hat e_i^{(a)}+\text{c.c.},\qquad \sum_a |c^{(a)}|^2=1, \] with \(\omega=m-\mu\), \(X(r)\) a real spherically symmetric profile, and \(c^{(a)}\) specifying the polarization or spin state [2301.11470]. Empirically, the profile is fitted by \[ X(r)\simeq \frac{X_0}{[1+0.077\,(\mu/m)\,(mr)^2]^4},\qquad X_0\simeq 2.04\,\frac{\mu}{m}, \] and the mass, binding energy, and radius scale as \[ M\simeq 62.3\,\frac{M_{\rm Pl}^2}{m}\left(\frac{\mu}{m}\right)^{1/2},\quad |E|\simeq 20.8\,\frac{M_{\rm Pl}^2}{m}\left(\frac{\mu}{m}\right)^{3/2},\quad R\simeq 3.16\,\frac1m\left(\frac{\mu}{m}\right)^{-1/2} \] [2301.11470]. The choice of \(c^{(a)}\) distinguishes, for example, linear polarization from a purely circular soliton with \(c^{(+)}=1\). The cosmological analysis of dark photon stars gives a complementary structural description. Stationary, minimal-energy solitons obey \[ M\approx \frac{2\alpha}{Gm},\qquad R\approx \frac{1.9}{\alpha m},\qquad MR\approx \frac{3.9}{Gm^2}, \] so \(R\propto 1/M\) and the central density scales as \(\rho_s\propto G^3m^6M^4\) [2203.10100]. The physical radius is \[ R_s\approx 5\times10^3\,\mathrm{km}\, \left(\frac{m}{10^{-5}\,\mathrm{eV}}\right)^{-1/2} \left(\frac{M}{10^{-23}M_\odot}\right), \] and the density profile is well matched by \[ \rho_{\rm soliton}(r)=\frac{\rho_s}{[1+(0.228\,mr)^2]^4}. \] Outside the core, the surrounding fuzzy halo is described by an NFW form with \(r_0\approx 1.56\,\lambda_J(\rho_s)\) and \(\rho_0\approx 0.042\,\rho_s\) [2203.10100]. Self-interactions and nonminimal gravitational interactions change these relations in polarization-dependent ways. For a stationary polarized ansatz \[ \psi_i(t,\mathbf x)=e_i(\sigma)\,f(r)e^{i\mu t}, \] with \(e_i(\sigma)\in\{(0,0,1),(1,\pm i,0)/\sqrt2\}\), the radial profile equation contains \[ \lambda_\sigma \equiv \left(1-\frac{\sigma^2}{3}\right)\lambda, \] so the effective quartic self-interaction depends on polarization [2406.05031]. In the minimal case one obtains the approximate analytic relations \[ f(r)\approx \frac{f_0^2}{(1+0.0529\,f_0^2 r^2)^4},\qquad M\simeq 25.9\,f_0,\qquad R\simeq \frac{3.80}{f_0}, \] while in the strong repulsive Thomas–Fermi limit with \(\xi=0\), \[ R\simeq \frac{\pi}{2}\sqrt{\lambda_\sigma} \] [2406.05031]. Stability is controlled by the mass–radius curve. The linear criterion is \[ \frac{dM}{dR}<0, \] so the negative-slope branch is stable and the turnover point \(dM/dR=0\) marks the onset of instability [2406.05031]. Attractive self-interactions \(\lambda_\sigma<0\) yield a maximum mass and minimum radius, while \(\xi<0\) produces an analogous turnover from nonminimal gravity. For \(\xi>0\), regularity imposes an absolute amplitude bound, \[ f_0<\frac{1}{2|\xi|}, \] because the effective inertia denominator \(1-4\xi^2 f^2(r)\) must remain positive [2406.05031]. This bound caps the central amplitude and therefore constrains the high-mass tail of the soliton distribution. ## 5. Couplings to photons and radiative processes A central question is whether dark photon solitons emit observable electromagnetic radiation. In vacuo, gauge–kinetic mixing \[ \mathcal L_{\rm int}^{(4)}\supset \frac{\varepsilon}{2}F_{\mu\nu}X^{\mu\nu} \] does **not** source photon emission from a coherent dark photon condensate [2301.11470]. The leading interactions that do so in the isolated-soliton problem are instead dimension-6 operators \[ \mathcal L_{\rm int}^{(6)}=\sum_{i=1}^5 g_i^2 O_i, \] with five independent gauge- and Lorentz-invariant structures built from \(F_{\mu\nu}\), \(\tilde F_{\mu\nu}\), and \(X_\mu\) [2301.11470]. In the time-periodic background of a vector soliton, the photon mode equations acquire periodic coefficients, and Floquet analysis shows exponential growth for modes near the resonance band \(k\simeq \omega\simeq m\), \[ A_j(t)\sim e^{\mu_{\mathbf k}t},\qquad \Re\{\mu_{\mathbf k}\}>0. \] For \(O_{1\text{–}4}\), the maximal growth rate in the linear-polarization case is \[ \mu_{\max}^{(\mathrm{lin})}\simeq \frac12 g^2X_0^2 m, \] whereas for circular polarization only \(O_3\) and \(O_4\) resonate and \(O_5\) is suppressed at \(O(g^4)\) [2301.11470]. The emitted radiation carries operator- and polarization-specific signatures. The central frequency is \[ \nu_0=\frac{m}{2\pi}\approx 200\,\mathrm{MHz}\left(\frac{m}{10^{-6}\,\mathrm{eV}}\right), \] with bandwidth \[ \Delta \nu \sim \frac{g^2X_0^2 m}{2\pi}. \] For \(O_1\) and \(O_2\), the emission is isotropic and unpolarized; for \(O_3\) and \(O_4\) with a linearly polarized soliton along \(\hat z\), it peaks in the equatorial plane \(\theta=\pi/2\), with \(O_3\) linearly polarized parallel to \(X\) and \(O_4\) polarized in the azimuthal \(\hat\phi\) direction; for a circular soliton in the \(xy\) plane, only \(O_3\) and \(O_4\) resonate, producing circularly polarized radiation along \(\pm \hat z\) with the same handedness as the soliton [2301.11470]. The resonance condition is \(\mu_{\max}^{(\mathrm{hom})}R\gtrsim 1\), equivalently \(g^2X_0^2\gtrsim (\mu/m)^{1/2}\) [2301.11470]. A distinct radiative channel appears when external electromagnetic fields or charge densities are present. In that case, a dimension-6 dipole operator \[ \mathcal L_{\rm int}^{(6)}=-\frac{1}{4\Lambda^2}(A'_\mu A'^\mu)F_{\alpha\beta}\tilde F^{\alpha\beta} \] allows linearly polarized solitons to induce an oscillating charge/current at \(2\omega\), while gauge kinetic mixing implies \({\cal L}\supset -e\varepsilon J_\mu A'^\mu\), so electrons in the soliton background experience an oscillating Lorentz force and produce a current at \(\omega\simeq m\) [2601.05351]. For both mechanisms, plasma effects are crucial: in the resonant limits \(\omega_p\to 2m\) for the dipole channel and \(\omega_p\to m\) for the mixing channel, the exponential suppression from finite soliton size is removed and replaced by a milder power-law behavior [2601.05351]. This suggests that compact magnetospheres and dense plasma environments can be more important observationally than empty space. ## 6. Observational signatures, survival, and open distinctions The survival of dark photon solitons depends both on gravitational disruption and on radiative decay. From the gravitational side, a clump of mean density \(\langle\rho\rangle\) survives tidal forces from a host of mean density \(\langle\rho_{\rm host}\rangle\) if \(\langle\rho\rangle\gtrsim 0.05\,\langle\rho_{\rm host}\rangle\). In the Milky Way, where \(\langle\rho_{\rm host}\rangle\sim 10^{-5}\,\mathrm{eV}^4\), the soliton cores with \(\langle\rho\rangle\sim 0.1\)–\(100\,\mathrm{eV}^4\) are robust, and fuzzy halos survive down to local densities \(\rho(r)\gtrsim 10^{-2}\,\mathrm{eV}^4\) [2203.10100]. The same study argues that, at minimum, the solitons are likely to survive to the present day without being tidally disrupted [2203.10100]. Radiative stability is more model-dependent. If the parametric-resonance channel is open, an isolated soliton decays on a timescale \[ \tau\sim \frac{2}{g^2X_0^2 m}\ll t_0, \] so long-lived present-day objects require \(g^2X_0^2<(\mu/m)^{1/2}\) and therefore lie below the critical mass \[ M