Thermal Deformations in Super-Eddington Magnetized Neutron Stars: Implications for Continuous Gravitational-Wave Detectability

Basic Information

• Paper ID: 2510.10443
• Title: Thermal Deformations in Super-Eddington Magnetized Neutron Stars: Implications for Continuous Gravitational-Wave Detectability
• Authors: Hong-Bo Li, Yacheng Kang, Ren-Xin Xu (Peking University)
• Classification: astro-ph.HE (High Energy Astrophysics), gr-qc (General Relativity and Quantum Cosmology)
• Publication Date: October 14, 2025 (Draft version)
• Paper Link: https://arxiv.org/abs/2510.10443

Abstract

Rapidly rotating neutron stars represent promising targets for continuous gravitational-wave searches with current and next-generation ground-based gravitational-wave detectors. This paper presents the first investigation of thermal deformations in super-Eddington accreting magnetized neutron stars, wherein magnetic field-induced anisotropic heat conduction leads to asymmetric crustal temperature distributions. The study calculates the resulting mass quadrupole moments and estimates the associated continuous gravitational-wave strain amplitudes. Results demonstrate that Galactic magnetized neutron stars undergoing super-Eddington columnar accretion can emit detectable continuous gravitational waves in forthcoming observatories. Assuming 2-year coherent integration, the Einstein Telescope and Cosmic Explorer can detect such continuous gravitational-wave signals from rapidly rotating neutron stars with spin periods P ≲ 20 ms, while LIGO O5 operations may detect systems with P ≲ 6 ms.

Research Background and Motivation

Problem Definition

The core scientific question addressed by this research is: Can thermal deformations produced in super-Eddington magnetized neutron stars during columnar accretion serve as detectable continuous gravitational-wave sources?

Research Significance

1. Astrophysical Implications: The discovery of ultraluminous X-ray sources (ULXs) challenges conventional accretion theory, particularly the existence of pulsating ULXs (ULXPs) indicates that neutron stars can achieve super-Eddington luminosities
2. Gravitational-Wave Astronomy: Continuous gravitational-wave detection represents an important frontier in current gravitational-wave astronomy, and identifying new gravitational-wave source types is of considerable significance
3. Multi-Messenger Astrophysics: Combining accretion physics with gravitational-wave astronomy provides novel pathways for probing neutron star internal structure

Limitations of Existing Approaches

1. Electromagnetic Observation Constraints: Electromagnetic observations alone cannot fully determine the intrinsic properties of super-Eddington magnetized neutron stars
2. Theoretical Model Uncertainties: The physical origin of extreme luminosities in ULXPs remains controversial, involving a broad range of magnetic field strengths
3. Lack of Systematic Studies: Previous work lacks quantitative investigation of thermal deformations in super-Eddington magnetized neutron stars

Research Motivation

Based on the theory that rapidly rotating neutron stars can generate continuous gravitational waves through mass quadrupole moments, the authors propose that magnetic field-guided columnar accretion may induce non-axisymmetric crustal deformations in rapidly rotating neutron stars, thereby producing significant continuous gravitational-wave signals.

Core Contributions

1. First Systematic Investigation: First quantitative study of thermal deformation mechanisms in columnar accretion onto super-Eddington magnetized neutron stars
2. Physical Model Development: Establishes a temperature asymmetry model accounting for magnetic field-modulated anisotropic heat conduction
3. Gravitational-Wave Predictions: Calculates corresponding mass quadrupole moments and assesses continuous gravitational-wave detectability
4. Observational Predictions: Predicts that rapidly rotating super-Eddington magnetized neutron stars in the Galaxy may constitute a new class of continuous gravitational-wave sources
5. Detector Sensitivity Analysis: Systematically evaluates detection capabilities of LIGO, Einstein Telescope, and Cosmic Explorer for such sources

Methodology Details

Task Definition

Input: Physical parameters of super-Eddington magnetized neutron stars (mass, radius, magnetic field strength, accretion rate, spin period) Output: Mass quadrupole moment Q₂₂ and continuous gravitational-wave strain amplitude h₀ produced by thermal deformation Constraints: Magnetic field strength B ≲ 10¹³ G, steady-state thermal equilibrium assumption

Model Architecture

1. Accretion Column Structure Model

Based on radiation pressure-dominated shock theory, the governing equations are:

dε/dr = (3/4)(ε/rD - GM/r²)     (1)
dρ/dr = -3ρ/4ε(ε/rD + 3GM/r²)  (2)

where ε is the specific radiation energy, ρ is mass density, and rD = 1.4×10⁵(Ṁ/10¹⁷ g s⁻¹) cm is the photon diffusion length.

2. Heat Conduction Model

Heat flux is described using Fourier's law:

F = -κ∇T                        (3)

where κ is the thermal conductivity, and the presence of magnetic fields renders the thermal conductivity anisotropic.

3. Temperature Perturbation Calculation

Spherical harmonic expansion is employed:

δT(r,θ,φ) = Σ δTℓm(r)Yℓm(θ,φ)   (4)

with emphasis on the quadrupole component (ℓ=m=2).

4. Crustal Deformation Calculation

Mass quadrupole moment expression:

Q₂₂ = ∫ (∂lnρ/∂lnT)|p (δT/T) ρr⁴ dr   (5)

Technical Innovations

1. Multi-Physics Coupling: First unified modeling of super-Eddington accretion, magnetic field effects, heat conduction, and structural deformation
2. Self-Consistent Numerical Method: Employs the self-consistent numerical method developed by Li et al. (2025) to calculate mass quadrupole moments
3. Parametric Description: Proposes simplified expression for ellipticity: ε ≃ 2×10⁻⁷B₁₂
4. Systematic Analysis: Complete physical chain from accretion column structure to gravitational-wave detection

Experimental Setup

Physical Parameters

• Neutron Star Models: M = 1.4M☉, R = 10⁶ cm (Section 2); M ≃ 1.84M☉, R ≃ 1.17×10⁶ cm (temperature perturbation calculations, using DH equation of state)
• Magnetic Field Strength: B = 10¹² - 10¹³ G
• Accretion Rate: Ṁ = 10¹⁷ - 10¹⁹ g s⁻¹
• Spin Period: P = 3 - 20 ms

Numerical Methods

• Boundary Conditions: Column top determined by free-fall region conditions, bottom determined by gas pressure-magnetic pressure balance
• Integration Strategy: Integration from initial guess at rS inward to stellar surface
• Convergence Criterion: Self-consistency of bottom boundary conditions

Comparison Benchmarks

• Known ULXPs: 12 observed pulsating ultraluminous X-ray sources listed in Table 1
• Typical LMXB Systems: Temperature perturbations compared with low-mass X-ray binary systems
• Gravitational-Wave Detector Sensitivities: LIGO O5, Einstein Telescope, Cosmic Explorer

Experimental Results

Main Results

1. Accretion Column Structure (Figure 1)

• Extremely high accretion rates produce taller accretion columns
• Column base temperature shows weak dependence on accretion rate, with characteristic temperature ~4×10⁹ K
• Significantly exceeds typical LMXB system temperatures

2. Temperature Perturbations (Figure 2)

• Relative temperature perturbation δT/T increases with magnetic field strength
• At B = 10¹² G: δT/T ≈ 0.02
• At B = 10¹³ G: δT/T ≈ 0.14
• Up to 4 orders of magnitude higher than typical LMXB systems

3. Mass Quadrupole Moments

• B = 10¹² G: Q₂₂ ≃ 0.4×10³⁸ g cm²
• B = 10¹³ G: Q₂₂ ≃ 1.2×10³⁸ g cm²
• Scales linearly with magnetic field strength

4. Gravitational-Wave Detectability (Figure 3)

Known ULXPs: Undetectable by current detectors due to large distances and insufficient spin frequencies

Predicted Galactic NS ULXs:

• Einstein Telescope & Cosmic Explorer: Can detect systems with P ≲ 20 ms
• LIGO O5: Can detect systems with P ≲ 6 ms
• PSR J1928+1815-like systems (P = 10.55 ms): Detectable by next-generation detectors

Key Findings

1. Temperature Enhancement Effect: The combination of high temperatures from super-Eddington accretion and strong magnetic fields produces significantly enhanced temperature perturbations
2. Ellipticity Parametrization: ε ≃ 2×10⁻⁷B₁₂ provides simplified detectability assessment
3. Frequency Dependence: Detectability strongly depends on spin period, with millisecond pulsars being most promising
4. Importance of Galactic Sources: Distance is the critical limiting factor for detection

Related Work

Continuous Gravitational-Wave Theory

• Elastic Strain Mechanisms: Haskell et al. (2006), Johnson-McDaniel & Owen (2013)
• Magnetic Field Deformations: Melatos & Payne (2005), Haskell et al. (2008)
• r-Mode Instability: Andersson et al. (1999), Levin (1999)

ULX and ULXP Research

• Observational Discoveries: Discovery of M82 X-2 (Bachetti et al. 2014) initiated ULXP research
• Theoretical Models: Multiple mechanisms from geometric beaming to magnetic field suppression of radiation pressure
• Population Synthesis: Predicts 7-20 NS ULXs in the Galaxy

Heat Conduction and Magnetic Field Effects

• Anisotropic Heat Conduction: Geppert et al. (2004, 2006), Page et al. (2007)
• Temperature Perturbation Methods: Spherical harmonic expansion method of Osborne & Jones (2020)

Conclusions and Discussion

Main Conclusions

1. New Gravitational-Wave Source Type: Super-Eddington magnetized neutron stars may represent a new class of continuous gravitational-wave sources
2. Detection Prospects: Next-generation gravitational-wave detectors possess the capability to detect rapidly rotating systems in the Galaxy
3. Multi-Messenger Opportunities: Provides unique opportunities for probing neutron star crustal structure and connecting accretion physics with gravitational-wave astronomy

Limitations

1. Magnetic Field Strength Constraints: Method applicable for B ≲ 10¹³ G
2. Steady-State Assumption: Assumes the system is in thermal equilibrium
3. Simplified Magnetic Field Configuration: Does not account for complex non-dipolar magnetic field structures
4. Limited Observational Sample: No rapidly rotating NS ULX has yet been discovered in the Galaxy

Future Directions

1. More Realistic Magnetic Field Models: Including complex internal magnetic field configurations
2. Improved Microphysics: Updated microscopic physical models for accreting crusts
3. Optimized Search Strategies: Optimization of search algorithms for future gravitational-wave detectors
4. Population Studies: Further exploration of occurrence rates of rapidly rotating super-Eddington magnetized neutron stars

In-Depth Evaluation

Strengths

1. High Originality: First systematic investigation of continuous gravitational-wave emission from this class of objects
2. Physical Completeness: Complete physical chain from microscopic heat conduction to macroscopic gravitational-wave radiation
3. Numerical Rigor: Employs self-consistent numerical methods and reasonable physical assumptions
4. Observational Relevance: Directly interfaces with capabilities of current and future gravitational-wave detectors
5. Predictive Value: Provides clear targets and expectations for future observations

Weaknesses

1. Model Simplifications: Simplified magnetic field configuration and heat conduction may affect quantitative prediction accuracy
2. Parameter Uncertainties: Intrinsic parameters of ULXPs (particularly magnetic field strength) exhibit substantial uncertainties
3. Lack of Statistical Analysis: Insufficient consideration of statistical impacts of parameter uncertainties
4. Observational Verification Challenges: Predicted Galactic sources have not yet been observed

Impact

1. Interdisciplinary Connection: Successfully bridges high-energy astrophysics and gravitational-wave astronomy
2. Observational Guidance: Provides theoretical guidance for gravitational-wave detector search strategies
3. Theoretical Advancement: Advances understanding of neutron star physics in extreme environments
4. Technical Inspiration: Numerical methods applicable to other types of neutron star deformation studies

Applicable Scenarios

1. Gravitational-Wave Data Analysis: Provides new candidate source types for continuous gravitational-wave searches
2. Multi-Messenger Astrophysics: Joint analysis combining X-ray and gravitational-wave observations
3. Neutron Star Physics Research: Probes neutron star internal structure in extreme magnetic field environments
4. Population Synthesis Models: Improves theoretical models of binary evolution and neutron star formation

References

The paper cites extensive relevant literature spanning ULX observations, neutron star physics, gravitational-wave theory, and numerical methods, reflecting the comprehensive and cutting-edge nature of the research. Key references include Bachetti et al. (2014) on the discovery of the first ULXP, Osborne & Jones (2020) on temperature perturbation methods, and the authors' previous work Li et al. (2025) on self-consistent numerical methods.