In the formation of massive stars (MS), the classical prediction of the radiation pressure problem requires that stars above ≈20 M⊙ cannot form because their own luminosity would halt further accretion, remains a central puzzle in MS formation. We present a unified analytic model that synthesizes four key mechanisms: anisotropic radiation escape (flashlight effect), diffusive flux transport in optically thick regions, temperature‑dependent dust opacity, and a power‑law density profile in the infalling envelope. The resulting force balance yields a maximum stellar mass scaling \ (M_*=100 f () ^3{4} (Ṁ10^{-3 M_ yr^-1}) ^4{13} (R10 ~cm² ~g^{-1}) ^-4{13} M_, \), where f (θ) ≈ 0. 2 is the geometric reduction factor in a disk geometry, is the mass accretion rate, and is the Rosseland mean opacity, represents the temperature power-law index. The weak exponent is the geometric reduction factor in a disk geometry, \ (Ṁ\) is the mass accretion rate, and is the Rosseland mean opacity, represents the temperature power-law index. The weak exponent \ (113\) ≈ 0. 308 implies that increasing \ (Ṁ\) by a factor of 10 raises M* by only a factor of ≈ 2. 03, while doubling reduces \ (Ṁ\) by ≈ 1. 24. Numerical evaluation for f (θ) = 0. 2 and Tsub = 1500 K (Tsub is the sublimation temperature) gives \ (Ṁ\) ≈3. 6 M⊙ at \ (Ṁ\) =10-6M⊙ yr-1, KR =10 cm2 g−1; rising to \ (Ṁ\ = 10-2M⊙ yr-1, KR =2 cm2 g−1. These values span the observed stellar initial mass function from intermediate‑mass to the most MSs known and favor the formation of MS in the early universe, due to low metallicity. The model demonstrates that no single effect overcomes the radiation barrier; instead, multiplicative suppression through anisotropy, high accretion rates, and reduced dust opacity (e. g. , in low‑metallicity environments) allows gravity to dominate. Observational tests include resolved imaging of polar cavities in high‑mass protostars and measurements of dust opacities in star‑forming clouds. The derived scaling provides a physically motivated, predictive criterion for the maximum stellar mass achievable under given infall and opacity conditions. The limitation of the ADR-RB lies in its simplified 1D formulation, which neglects magnetohydrodynamic effects, dynamic dust opacity evolution, and multidimensional processes such as stellar multiplicity and Rayleigh–Taylor instabilities that critically regulate angular momentum transport, accretion flow, and anisotropic radiation escape.
Onuchukwu et al. (Thu,) studied this question.