Inferring hemispheric asymmetries of stellar active regions through the information content of astrometric signals

1. Problem Statement

Stellar surface mapping (reconstructing the distribution of starspots, faculae, and magnetic activity) is a fundamental inverse problem in astrophysics. Current techniques, primarily photometry (light curves) and Doppler imaging, suffer from severe degeneracies that prevent unique surface reconstruction:

• Photometry: Integrates flux over the visible disk. It is fundamentally insensitive to odd-degree spherical harmonic modes ($\ell \geq 3$), which encode north-south asymmetries. It also suffers from a "north-south" degeneracy where spots mirrored across the equator produce identical light curves.
• Doppler Imaging: Suffers from latitude smearing near the equator and cannot distinguish the sign of latitude (north vs. south) for features mirrored across the equator.
• General Limitation: As spatial resolution increases (higher spherical harmonic degree $\ell$), the fraction of recoverable modes decreases asymptotically because the "null space" (unobservable modes) grows quadratically while the observable rank grows only linearly.

The paper addresses the need for a complementary technique to break these degeneracies, specifically focusing on high-precision astrometry (measuring the photocentre position of a star).

2. Methodology

The authors extend the theoretical framework established by Luger et al. (2021) for photometric signals to include astrometric signals.

A. Mathematical Framework

• Spherical Harmonic Decomposition: The stellar surface brightness $I(\theta, \phi)$ is decomposed into spherical harmonics $Y_\ell^m$.
• Linear Observation Model: The mapping from surface coefficients to observations is modeled as a linear system $\mathbf{y} = \mathbf{A}\mathbf{c}$, where $\mathbf{A}$ is the design matrix.
• Selection Rules (Analytical Derivation):
  • Photometry: The authors re-derive that the photometric operator annihilates all odd-degree modes ($\ell \geq 3$) due to symmetry.
  • Astrometry: They derive the astrometric operator (flux-weighted centroid position). Crucially, they prove analytically that the position-weighting ($x, y$ coordinates) breaks the symmetry of odd modes. Specifically, they demonstrate that astrometry is sensitive to odd-degree spherical harmonics ($\ell \geq 3$), which are invisible to photometry.
• Fisher Information Matrix (FIM): To quantify the constraining power of observations, the authors use the FIM ($\mathbf{F} = \mathbf{A}^T\mathbf{A}$). They calculate the posterior shrinkage factor ($S$), which measures how much the uncertainty of a mode coefficient is reduced by the data (where $S=1$ is fully constrained and $S=0$ is unconstrained).

B. Numerical Simulation

• Discretization: The stellar surface is discretized using a Fibonacci spiral packing to ensure uniform point distribution without area weighting.
• Simulation: They simulate light curves and photocentre trajectories ($X_c, Y_c$) for various viewing inclinations ($i$) and rotation phases.
• Rank Analysis: Using QR decomposition with column pivoting, they determine the numerical rank and nullity of the design matrices for photometry, astrometry (X and Y components), and combined observations.
• Reconstruction Tests: They test surface recovery on synthetic maps containing:
  1. A single spot.
  2. Multiple random spots.
  3. A complex text pattern ("GREETINGS EARTHLINGS").
  4. Scenarios with single vs. multiple viewing inclinations.

3. Key Contributions

1. Analytical Proof of Complementarity: The paper provides a rigorous proof that astrometry accesses a distinct set of spherical harmonic modes (specifically odd-$\ell$ modes) that lie entirely in the photometric null space. This allows for the detection of hemispheric asymmetries (e.g., active belts in one hemisphere vs. the other).
2. Quantification of Information Content: They demonstrate that while the fraction of recoverable modes still decreases with increasing spatial resolution ($\ell$), the combined rank of photometry and astrometry grows faster than either technique alone. The two techniques are not redundant; they provide orthogonal constraints.
3. Reframing Astrometric Jitter: The authors argue that "astrometric jitter" (noise in exoplanet searches caused by stellar activity) should be treated as a signal encoding surface structure. This is particularly valuable for missions like Gaia and TOLIMAN, where this jitter is a byproduct of planet hunting.
4. Inclination Dependence: They map the constraining power of specific modes against stellar inclination, showing that constraints weaken for pole-on views but that ensemble observations (many stars at random inclinations) can statistically recover population-level properties (e.g., typical spot latitudes).

4. Key Results

• Mode Recovery:
  • Photometry: Only recovers even-degree modes ($\ell$).
  • Astrometry: Recovers both even and odd-degree modes.
  • Combined: Recovers the union of modes, significantly increasing the total number of accessible degrees of freedom.
• Breaking Degeneracies:
  • North-South: Astrometry breaks the north-south degeneracy inherent in photometry because odd-$\ell$ modes encode the asymmetry.
  • Bright/Dark: While astrometry alone has a bright/dark spot degeneracy (a dark spot at location A vs. a bright spot at location B), combining it with photometry (which constrains total flux deficit) resolves this ambiguity.
• Surface Reconstruction:
  • Single Spot: Combined data localizes the spot in both latitude and longitude with high fidelity, whereas individual methods suffer from smearing or ambiguity.
  • Complex Structures: For the "GREETINGS EARTHLINGS" text test, neither photometry nor astrometry alone could recover the fine structure. However, the combined data successfully reconstructed the message, demonstrating that the information content is sufficient for complex mapping if the signal-to-noise ratio is high enough.
• Inclination Effects: Recovery is best for equator-on views ($i \approx 90^\circ$) and degrades for pole-on views, though the combined method still outperforms single-method approaches at all inclinations.

5. Significance and Applications

• Exoplanet Science: Understanding stellar surfaces is critical for removing stellar activity signals that mimic or mask exoplanet signals (especially Earth-analogs). By characterizing the surface via astrometry, the "jitter" can be modeled and subtracted, improving planet detection limits.
• Stellar Physics: The ability to detect odd-$\ell$ modes allows astronomers to study hemispheric asymmetries in magnetic dynamos, persistent active latitudes, and polar spots, which are impossible to constrain with photometry alone.
• Mission Design:
  • Immediate Application: Gaia (specifically the upcoming DR4 release) can map evolved giant stars, which have large radii and large, long-lived spots, producing detectable microarcsecond signals.
  • Future Missions: Next-generation missions targeting sub-microarcsecond precision (e.g., TOLIMAN, Theia, MASS, CHES, ARMADA) will be capable of mapping Sun-like stars. The paper suggests that these missions will naturally provide surface maps as a byproduct of their exoplanet searches.
• Methodological Shift: The paper advocates for treating astrometric jitter not as noise to be minimized, but as a rich data source for stellar surface physics, encouraging the joint analysis of photometric and astrometric time series.

In conclusion, Deagan & Montet demonstrate that high-precision astrometry is a powerful, complementary tool to photometry that breaks fundamental degeneracies in stellar surface mapping, enabling the recovery of hemispheric asymmetries and improving the characterization of both stars and their planetary systems.