Polarization from Rapidly Rotating Massive Stars

Basic Information

• Paper ID: 2510.12500
• Title: Polarization from Rapidly Rotating Massive Stars
• Authors: J. Patrick Harrington (University of Maryland), Richard Ignace (East Tennessee State University), K. G. Gayley (University of Iowa)
• Classification: astro-ph.SR (Solar and Stellar Astrophysics)
• Publication Date: October 14, 2025
• Paper Link: https://arxiv.org/abs/2510.12500

Abstract

This study diagnoses the rotational criticality of rapidly rotating massive stars by computing continuum polarization. The research combines the effects of rotational stellar distortion (oblateness) and stellar flux redistribution (gravity darkening), comparing polarization produced by the von Zeipel approximation and the Espinosa Lara & Rieutard (ELR) method. Results demonstrate that in the far-ultraviolet spectral region, photospheric polarization of B-type stars rises significantly, reaching approximately 1% maximum polarization at 140 nm for early main-sequence and subgiant stars rotating at 98% critical velocity. The predicted signal strengths indicate that future deep ultraviolet spectropolarimetric instruments will enable precise determination of near-critical rotation in B-type stars.

Research Background and Motivation

Core Questions

1. Diagnosis of stellar rotational criticality: How to determine the degree to which massive stars approach rotational breakup, which is crucial for understanding the formation mechanisms of Be stars (with disk-like emission lines) and Bn stars (with rotationally broadened absorption lines)
2. Angular momentum transfer mechanisms: Understanding how angular momentum transfer leads to the production of two important categories of rapidly rotating massive stars
3. Mass transfer conservative limits: The impact of rotation on mass transfer conservatism in binary star systems

Research Significance

• Massive stars commonly exist in binary or multiple star systems, with a substantial fraction undergoing mass transfer evolution
• Massive stars are typically rapid rotators (>100 km/s), with rotation inducing rotational mixing
• Be stars require near-critical rotation to form Keplerian disks, yet observations show approximately half of Be stars rotate below 70% of the breakup velocity

Limitations of Existing Methods

• Current uncertainties in rotational velocity determination are substantial
• Optical interferometric observations, while showing some stars near critical rotation (>90%), have limited samples
• Traditional methods struggle to precisely distinguish rotational state differences between Be and Bn stars

Core Contributions

1. Established a systematic theoretical model of continuum polarization for rapidly rotating massive stars, combining rotational distortion and gravity darkening effects
2. Compared two gravity darkening theories (von Zeipel and ELR), finding that the ELR method produces approximately 25% lower polarization than von Zeipel
3. Identified the far-ultraviolet spectral region as the optimal wavelength range for detecting rotational polarization, with strongest polarization signals in the 120-200 nm band
4. Provided detailed polarization predictions for different spectral-type B stars, including B0V, B1IV, B3V, and B8V
5. Established theoretical foundation for future ultraviolet spectropolarimetric observation missions, supporting scientific objectives of missions such as CASSTOR and Polstar

Methodology Details

Task Definition

Computing continuum line polarization of rapidly rotating massive stars at different rotation velocities (W = Ω/Ωc, where Ωc is the critical angular velocity) and observational inclination angles, with stellar fundamental parameters (mass, radius, temperature, surface gravity) as input and wavelength-dependent polarization degree as output.

Model Architecture

1. Stellar Shape Model (Roche Model)

Adopting the Roche model with concentrated central mass to describe rotating stellar shape:

x(θ*, W) = (3/W sin θ*) × cos{(1/3)[π + cos⁻¹(W sin θ*)]}

where x = R(θ,Ω)/Rp is the normalized radius and θ* is the colatitude.

2. Gravity Darkening Model

von Zeipel Law:

Teff ∝ g^(1/4)

ELR Method: Introducing correction factor F(Ω̃, θ*), accounting for Eddington-Sweet circulation effects:

Teff(θ) ∝ [FΩ(θ) g(θ)]^(1/4)

3. Polarization Calculation

Using TLUSTY model atmospheres to extract opacity and scattering information, computing Stokes parameters I and Q:

LI = 8π Rp² ∫∫ μIλ(θ*, μ) x² sin θ*/cos(θ* - δ) dφ* dθ*
LQ = 8π Rp² ∫∫ μQλ(θ*, μ) cos(2ξ) x² sin θ*/cos(θ* - δ) dφ* dθ*

Technical Innovations

1. Unified treatment of shape distortion and gravity darkening: Simultaneously accounting for two symmetry-breaking mechanisms induced by rotation
2. First systematic application of ELR gravity darkening: More consistent with observations compared to traditional von Zeipel method
3. High-precision numerical integration: Dividing stellar surface into 1°×1° grid for accurate computation
4. Multi-spectral-type coverage: Systematically studying stars from B0V to B8V at different evolutionary stages

Experimental Setup

Stellar Sample

Four representative B-type stars selected:

• B0V: Teff=29,000K, log g=3.9, M=17.5M⊙
• B1IV: Teff=25,000K, log g=3.6, M=16.3M⊙
• B3V: Teff=18,000K, log g=4.0, M=8.0M⊙
• B8V: Teff=12,000K, log g=4.0, M=3.8M⊙

Model Atmospheres

• Early-type stars: Using TLUSTY model atmospheres
• B8V stars: Due to extreme rotation causing excessively low equatorial temperature, using Kurucz model with synspec205 program

Parameter Ranges

• Rotation velocity: W = 0.6-0.99 (60%-99% critical velocity)
• Observational inclination: i = 20°-90° (10° intervals)
• Wavelength range: 120-800 nm, with emphasis on 120-300 nm ultraviolet region

Evaluation Metrics

• Polarization degree: pλ = |LQ/LI| × 100%
• Wavelength dependence: Polarization variation with wavelength
• Angular dependence: Polarization variation with observational inclination

Experimental Results

Main Results

1. Polarization Amplitude

• Maximum polarization: B1IV star reaches ~1% at W=0.98, i=90°, λ=140nm
• Significant polarization threshold: W>0.8 produces >0.1% polarization at high inclination angles
• Low inclination detection: Even at i=40°, W=0.98 models produce >0.1% polarization below 200nm

2. Wavelength Dependence

• Far-ultraviolet strongest: 120-200 nm band shows strongest polarization signal, rapidly decreasing toward longer wavelengths
• Optical band weak: Visible light polarization <0.05%, difficult to observe
• Steep rising feature: Rise toward FUV band steeper than Wien approximation predictions

3. Spectral Type Differences

• B0V: Relatively weak polarization due to short wavelength proximity to Wien peak and high surface gravity
• B1IV: Strongest polarization, optimal combination of temperature and gravity
• B3V: Moderate polarization, Wien peak located in UV band
• B8V: Reaches 0.75% maximum polarization at W=0.98

von Zeipel vs ELR Comparison

• ELR polarization systematically ~25% lower than von Zeipel
• Similar wavelength dependence between two methods
• ELR method more consistent with interferometric temperature contrasts

Relative Contributions of Shape Distortion and Gravity Darkening

• Approximately 50% each: Two effects contribute roughly equally to total polarization
• Shape effect slightly dominant under ELR
• Gravity darkening effect slightly stronger under von Zeipel

Related Work

Historical Development

• Harrington & Collins (1968): First calculation of intrinsic polarization of rapidly rotating early-type stars
• Collins et al. (1991): Established classical theoretical framework for rotating star polarization
• Espinosa Lara & Rieutord (2011): Proposed improved gravity darkening theory

Recent Applications

• Regulus: Cotton et al. (2017)
• ζ Pup: Bailey et al. (2024a)
• ε Sgr: Bailey et al. (2024b)
• ζ Aql: Howarth et al. (2023)

Future Missions

• Arago Mission Concept: Muslimov & Neiner (2023)
• Polstar Mission: Drake et al. (2025)
• CASSTOR: Neiner et al. (2025)
• POLLUX Instrument: Girardot et al. (2024a)

Conclusions and Discussion

Main Conclusions

1. FUV polarization is the optimal tool for rotation diagnostics: 120-200 nm band shows strongest and most characteristic polarization signals
2. Near-critical rotation is detectable: Rapidly rotating stars with W>0.8 can be reliably detected with appropriate instruments
3. ELR method is more realistic: Although polarization is slightly lower than von Zeipel, it is more consistent with observational constraints
4. Future observation prospects are promising: Planned UV polarimetric missions have capability to achieve scientific objectives

Limitations

1. Model assumptions: Adopting rigid body rotation and Roche approximation; actual stars may exhibit differential rotation
2. Atmospheric model limitations: Applicability of model atmospheres under extreme rotation conditions requires verification
3. Disk effects not considered: Circumstellar disks of Be stars produce additional polarization signals
4. Limited sample: Only four representative spectral types considered

Future Directions

1. Extension to cooler B-type stars: Investigating richer population of disk-free Bn stars
2. Combined disk polarization modeling: Addressing joint polarization of Be star photosphere and disk
3. Surface feature considerations: Effects of non-axisymmetric structures such as starspots
4. Observational verification: Validating theoretical predictions using future UV polarimetric missions

In-Depth Evaluation

Strengths

1. Theoretical completeness: Systematically addresses all major physical mechanisms of rotating star polarization
2. Advanced methodology: Employs latest ELR gravity darkening theory, enhancing physical realism of model
3. Forward-looking application: Provides important theoretical support for forthcoming space UV polarimetric observations
4. Result reliability: Consistent with historical calculations and recent observations, enhancing credibility

Weaknesses

1. Physical simplification: Rigid body rotation assumption may not apply to all rapidly rotating stars
2. Limited observational verification: Lacks large-sample observational data to validate theoretical predictions
3. Environmental effects neglected: Does not account for binary interactions, stellar winds, and other effects on polarization
4. Computational complexity: High-precision calculations require substantial computational resources, limiting parameter space exploration

Impact

1. Scientific value: Provides important observational diagnostic tool for understanding massive star rotational evolution
2. Technical advancement: Drives development needs for UV polarimetric observation technology
3. Application potential: Can be directly applied to interpretation and analysis of future observational data
4. Theoretical refinement: Advances theoretical understanding of rotational effects in stellar physics

Applicable Scenarios

1. Rotational velocity determination of rapidly rotating B-type stars
2. Be/Bn star classification and evolutionary studies
3. Binary mass transfer efficiency assessment
4. Scientific planning for future UV polarimetric observation missions

References

This paper cites extensive related research, primarily including:

• Theoretical foundations of rotating stars: von Zeipel (1924), Tassoul (1978), Maeder (2009)
• Polarization theory development: Harrington & Collins (1968), Collins et al. (1991)
• Gravity darkening improvements: Espinosa Lara & Rieutord (2011)
• Recent observational applications: Cotton et al. (2017), Bailey et al. (2024a,b), Howarth et al. (2023)
• Future mission planning: Drake et al. (2025), Neiner et al. (2025), Girardot et al. (2024a)

This research provides a solid theoretical foundation for polarimetric observations of rapidly rotating massive stars, with important significance for understanding stellar rotational evolution and guiding future observations. As space UV polarimetric observation capabilities develop, these theoretical predictions will receive direct observational verification.