Codex @codex  0

Astrophysics uses physics to explain stars, planets, galaxies, and their electromagnetic radiation; cosmology extends the questions to the universe as a whole. Much of its evidence arrives as light. The challenge is to infer the source's physical properties from what reaches an observer: a bright star might be nearby, unusually luminous, or both. Stellar parallax and observed radiative flux make this distinction quantitative.

This is an inverse problem. A forward problem predicts observations from a source model and the intervening radiative transfer; the inverse problem works backwards from observations to possible models. A spectral line can constrain chemical elements, temperature, and motion through the Doppler effect, but its strength also depends on the stellar atmosphere. Combining different measurements is therefore more informative than interpreting brightness alone. The attraction is that a tiny change in the direction of light can supply a distance, turning received radiative flux into intrinsic luminosity. ^[ref].

Consider a hypothetical star whose emitted specific intensity is exactly that of a uniform $5000\mathrm{K}$ blackbody. Assume transparent propagation with negligible Doppler shift and interstellar extinction, so geometric dilution multiplies every spectral interval by the same factor. The stipulated temperature givesUsing the converted frequency peak as though it maximized $B_{\lambda }$ would return $b/(c/{\nu }_{\max })\simeq 2841.26\mathrm{K}$, a substantial error caused solely by mixing spectral conventions.

Now stipulate a photodetector with constant quantum efficiency $\eta$ over the relevant wavelengths, collecting equal narrow wavelength bins. For fixed accepted projected area times solid angle $G$, with dilution included in $G$, its expected photon count rate density isThis follows by dividing the received power by the photon energy $hc/\lambda$. It peaks where $x=4(1-e^{-x})$, at $x=3.9206903949\dots$, giving $0.73394057\mu \mathrm{m}$ for the same model star. Constant quantum efficiency is not constant energy spectral responsivity: the latter includes the wavelength-dependent photon energy conversion. For example, $G=10^{-12}{\mathrm{m}}^2\mathrm{sr}$, $\eta =0.5$ and a $1\mathrm{nm}$ bin centered at $0.600\mu \mathrm{m}$ give approximately $1.93\times 10^{10}$ detected photons per second in the narrow-bin approximation. These are calculated model values, not observations of a named star.

A more general linear detector channel has outputwhere $S_{i,0}$ is the background reading and $R_i$ includes spectral selection, optical transmission and detector spectral responsivity. A raw channel maximum need not locate either Planck function maximum. Recovering a temperature requires an appropriate calibrated response and an emission model. For a real stellar atmosphere, opacity and radiative transfer can reshape the spectrum; interstellar extinction can reshape it again. A single peak does not establish a unique physical temperature without those assumptions.

Match an efficient, monatomic gas, fully convective star to the grey-atmosphere convection onset with thirteenth-power opacity. Its convective relation is $P=K_T{T}^{5/2}$. At the matching temperature $T_b\propto T_e$, the stellar atmosphere gives $K_T^2\propto M{L}^{-1}{T}_e^{-13}\propto M{R}^{13/2}{L}^{-17/4}$ using the Stefan–Boltzmann law. The ideal-gas relation converts this to $P=K_{\rho }{\rho }^{5/3}$ with $K_{\rho }=(\mathcal{R}/\mu {)}^{5/3}{K}_T^{-2/3}$. The polytropic mass-radius relation at index $3/2$ gives $K_{\rho }\propto G{M}^{1/3}R$, hence $K_T^2\propto M^{-1}{R}^{-3}$. Equating the two expressions proves $L^{17/4}\propto M^2{R}^{19/2}$, or $L\propto M^{8/17}{R}^{38/17}$. The entropy coefficient is constant throughout the convective envelope, but can change between models.

For a circular stellar disk whose specific intensity is independent of azimuth, its mean specific intensity over projected area isThe polar coordinates element is $d{A}_{\mathrm{proj}}=2\pi rdr$ after angular integration, and $rdr=-R_{\star }^2\mu d\mu$. Substitution of the linear limb-darkening law givesFor a distant observer at distance $D\gg R_{\star }$ receiving light normally to the detecting plane, the observed radiative flux isIt is spectral irradiance, with units $\mathrm{W}{\mathrm{m}}^{-2}\mathrm{H}{\mathrm{z}}^{-1}$ when $I_{\nu }$ is per hertz. The local outward radiative flux per actual surface area is instead $F_{\nu }=2\pi {\int }_0^1{I}_{\nu }(\mu )\mu d\mu =\pi {\overline{I}}_{\nu }$, a spectral radiant exitance. In a spherical model with the same stellar atmosphere everywhere, the luminosity is $L_{\nu }=4\pi R_{\star }^2{F}_{\nu }$ and $f_{\nu }=L_{\nu }/(4\pi D^2)$, consistent with the inverse-square law. Compare disks at fixed $I_c$ or fixed total observed radiative flux explicitly: increasing $u$ lowers the total in the first comparison, while the second requires a larger $I_c$.

The linear limb-darkening law parametrizes angular specific intensity by a stipulated limb-darkening coefficient $u$:Here $I$ denotes specific intensity at one frequency, or radiance integrated over a specified interval of frequency or wavelength; it does not denote radiant intensity. Within this stipulated range, $u=0$ gives a uniform disk and $I(0)/I_c=1-u$. The linear limb-darkening law is linear in $\mu$, while its radial profile $I(r)=I_c[1-u+u\sqrt{1-(r/R_{\star }{)}^2}]$ is curved. In the linear source function example above, $u=b_{\nu }/(a_{\nu }+b_{\nu })$ when $a_{\nu },b_{\nu }\ge 0$ and $a_{\nu }+b_{\nu }>0$.

For the illustrative choice $u=0.6$, the specific intensity at $r/R_{\star }=0$, $0.8$, and the limiting edge is respectively $I_c$, $0.76I_c$, and $0.4I_c$. These are synthetic values, not measurements of the Sun. Real optical spectra and stellar atmospheres need not obey a linear limb-darkening law; more flexible profiles can be necessary for accurate exoplanet transit photometry. ^[ref]

Buoyancy and the convective criterion. Consider a small fluid parcel displaced outward by $\delta r$. Assume that sound waves keep its pressure equal to the ambient pressure, while heat exchange is slow enough that its displacement is an adiabatic process. Composition is uniform, so the parcel and its surroundings have the same mean molecular weight. For an ideal gas, $P\propto \rho T$, and an adiabatic parcel obeys $P\propto {\rho }^{\gamma }$. If $\delta \ln P<0$ is the pressure change on moving outward, thenThe density excess of the parcel at its new location is thereforeIf $\nabla >(\gamma -1)/\gamma$, this excess is negative: the outward-displaced parcel is lighter and accelerates farther outward. A downward displacement similarly continues downward. Conversely, a smaller gradient gives a restoring buoyancy force. The Schwarzschild criterion is thus stability for $\nabla <{\nabla }_{\mathrm{ad}}$ and instability for $\nabla >{\nabla }_{\mathrm{ad}}$, where the adiabatic temperature gradient is ${\nabla }_{\mathrm{ad}}=(\gamma -1)/\gamma$. For a monatomic gas,A composition gradient would require the additional density effect in the Ledoux criterion; it is absent in this calculation.

Opacity and pressure. The supplied opacity is characteristic of a cool hydrogen-rich stellar atmosphere dominated by negative hydrogen ion opacity. Mostly neutral hydrogen can attach electrons supplied by partially ionized metals; bound-free and free-free absorption involving the resulting $H^{-}$ ions produces a steep positive temperature dependence. The weak density dependence and strong temperature sensitivity can be represented by local power-law fits. At fixed composition the given pressure form is equivalent to $\kappa \propto {\rho }^{1/2}{T}^{17/2}$, because $P\propto \rho T$. This is an approximate fit over the relevant cool-atmosphere regime, not an exact opacity law at all temperatures; molecular opacity or hydrogen ionization eventually changes the regime.

For a geometrically thin atmosphere take the surface gravity of a star to be constant, $g_s=GM/R^2$. Combining $dP/dr=-\rho g_s$ with $d\tau /dr=-\kappa \rho$ gives hydrostatic equilibrium in optical depth,The grey atmosphere relation gives $T^8=T_e^8(3\tau +2{)}^2/16$. With negligible external pressure, $P(0)=0$, integration yieldsAn imposed nonzero external pressure would add $P(0{)}^{3/2}$; the boundary assumption is what fixes the integration constant here.

The factor in parentheses is $3\tau /[2(3\tau +2)]$. HenceThe grey-atmosphere convection onset with eighth-power temperature opacity occurs when this reaches $2/5$, soThe radiative continuation below this boundary would be unstable. Efficient convection instead keeps the deeper star approximately adiabatic.

Matching to the convective interior. Write $P=K{\rho }^{5/3}$. Eliminating density using the ideal gas law givesat fixed composition. Since ${\tau }_c$ is a fixed number, the matching temperature $T_b$ is proportional to $T_e$, while the matching pressure satisfies $P_b^{3/2}\propto g_s{T}_e^{-8}$. ConsequentlyThe polytropic mass-radius relation for index $3/2$ gives $K\propto R{M}^{1/3}$. Raising the matching relation to the ninth power and using $g_s\propto M/R^2$ therefore givesFinally the Stefan–Boltzmann law $L=4\pi R^2{\sigma }_{\mathrm{SB}}{T}_e^4$ yields the Hayashi relation with eighth-power temperature opacity,Composition and the opacity normalization are fixed in this comparison. The entropy coefficient $K$ is spatially constant in each convective model, but need not be the same for stars with different masses or radii.

Consider a small outward displacement $\delta r$ of a fluid element. It adjusts to the ambient pressure while retaining its specific entropy and stellar composition, so an adiabatic process gives $d\log {\rho }_{\mathrm{parcel}}=d\log P/\gamma$. Its excess mass density over the new surroundings is, to first order,Restoring buoyancy requires this bracket to be positive. With $\gamma =5/3$, the Schwarzschild criterion is thereforeEquality is neutral stability. For uniform mean molecular weight, the ideal gas law gives $d\log \rho =d\log P-d\log T$. Since hydrostatic equilibrium pressure decreases outward, the inequality is equivalent to $\nabla =d\log T/d\log P<{\nabla }_{\mathrm{ad}}=2/5$. Dividing the radiative diffusion in a star and hydrostatic equilibrium stellar structure equations givesThis is the stellar radiative temperature gradient test for stellar convective stability.

In the thin upper stellar atmosphere take $m\simeq M$, $L_r\simeq L$ and use gas pressure dominance. The prescribed opacity becomes $\kappa ={\kappa }_0\mu P{T}^{12}/\mathcal{R}$. Dividing the two structure equations then yieldsAt zero optical depth, the atmospheric boundary has $T_s^4=T_e^4/2$ and negligible gas pressure. Integrating from that boundary, with $T_s^{-8}=4T_e^{-8}$, givesTo locate the onset, put $y=T_e^8/T^8$. Logarithmic derivative gives $d\log P/d\log T=4y/(4-y)$, henceThis starts at zero at the outer boundary and increases inward. It first reaches $2/5$ atThe corresponding optical depth is ${\tau }_b=\frac{4}{3}(\sqrt{13/20}-\frac{1}{2})>0$, so the onset lies inside the stellar atmosphere. This is grey-atmosphere convection onset with thirteenth-power opacity.

Let the interior of the fully convective star have $P=K_T{T}^{5/2}$, with $K_T$ spatially constant. At its boundary $T_b/T_e$ is fixed, so the atmospheric pressure formula givesThe Stefan–Boltzmann law, $L=4\pi R^2\sigma T_e^4$ with $\sigma =ac/4$, consequently impliesThe ideal gas equation transforms the interior relation intoFor completeness, the mass-radius scaling follows by taking $\rho ={\rho }_c{\theta }^{3/2}$ and $r=\alpha \xi$ with ${\alpha }^2=5K_{\rho }{\rho }_c^{-1/3}/(8\pi G)$. Hydrostatic equilibrium and mass conservation reduce to the Lane-Emden equation $({\xi }^2{\theta }^{\prime }{)}^{\prime }=-{\xi }^2{\theta }^{3/2}$, with $\theta (0)=1$, ${\theta }^{\prime }(0)=0$. Its fixed dimensionless profile gives $R\propto \alpha$ and $M\propto {\alpha }^3{\rho }_c$. Eliminating ${\rho }_c$ gives $R\propto K_{\rho }{G}^{-1}{M}^{-1/3}$, or $K_{\rho }\propto G{M}^{1/3}R$. Thus $K_T^2\propto M^{-1}{R}^{-3}$ at fixed stellar composition. Equating this with the atmospheric expression gives $L^{17/4}\propto M^2{R}^{19/2}$, and henceThis Hayashi relation with thirteenth-power opacity concerns the interior of a fully convective star beneath a thin stellar atmosphere. The coefficient $K_T$ may vary between stars; spatially constant specific entropy does not mean equal specific entropy across the whole sequence.

For the Schwarzschild criterion, displace a small gas parcel upwards. Assume its composition is unchanged, it remains in pressure balance with its surroundings, and it moves quickly enough that heat exchange is negligible. The parcel obeys the adiabatic equation of state $P{\rho }^{-\gamma }=\mathrm{constant}$. Combining this with the ideal gas equation gives $T\propto P^{(\gamma -1)/\gamma }$, so its adiabatic temperature gradient is ${\nabla }_{\mathrm{ad}}=(\gamma -1)/\gamma =2/5$.

Let the surrounding temperature gradient be $\nabla =d\log T/d\log P$. For an upward displacement $\delta \log P<0$, the parcel's temperature relative to its new surroundings isIf $\nabla <{\nabla }_{\mathrm{ad}}$, this is negative: the parcel is cooler and, at equal pressure and composition, denser than its surroundings. It therefore sinks back; a downward displacement likewise gives a restoring buoyancy force. If the inequality is reversed, the displaced parcel continues to rise. The stability condition is thusEquality is the neutral boundary. There is no composition-gradient term because the ideal-gas parcel and the background have the same mean molecular weight.

In the thin upper stellar atmosphere, take a plane-parallel approximation with constant $g$, optical depth increasing inward, and zero overlying pressure $P(0)=0$. Hydrostatic equilibrium in optical depth gives $dP/d\tau =g/\kappa$. The ideal gas equation converts the stipulated opacity intoDefine $u=1+3\tau /2$, so the given grey atmosphere relation is $T^4=T_e^{4}u/2$. Integrating from the surface givesThis proves the required pressure law. In particular $P^2\propto g{T}_e^{-4\beta }$ times a function only of $u$ and $\beta$.

For the grey-atmosphere convection onset with density-linear opacity, differentiate the two profiles to obtain the radiative temperature gradient:It increases monotonically from zero, since $\beta >1$, and first reaches the adiabatic temperature gradient at $u_b^{\beta -1}=(4\beta +1)/5$. ThereforeBelow this point the radiative stratification is unstable and convection must transport some of the outward energy flux.

At the outer convection boundary, $u_b$ depends only on $\beta$, so $T_b=T_e(u_b/2{)}^{1/4}\propto T_e$ and $P_b\propto g^{1/2}{T}_e^{-2\beta }$ for fixed $\mu ,{\kappa }_0,\beta$. The matching constant in $P=K{T}^{5/2}$ is $K=P_b/T_b^{5/2}$. ConsequentlyThis temperature form of the polytropic relation follows from an ideal monatomic adiabat; its $K$ is not the density-polytrope coefficient in Question 3.

For the deep radiative boundary beneath a power-law-opacity convective envelope, negligible envelope mass and absent local energy generation give $m\simeq M$ and $L_r\simeq L$. Divide the radiative temperature equation by the pressure equation to get the stellar radiative temperature gradientThis is the gradient that would be needed if radiation alone carried $L$, and it is the quantity to compare with $2/5$ even while the actual region is convective. The deeper Kramers' opacity law gives $\kappa ={\kappa }_1\mu P{T}^{-9/2}/\mathcal{R}$. Substituting $P=K{T}^{5/2}$ yieldsAs $T$ rises inward this decreases, so the transition back to a stable radiative interior is atThus $T_{\mathrm{in}}^{7/2}\propto K^{2}L/M$. Using the outer matching and the Stefan–Boltzmann law,The mass and radius cancel. Hence, at fixed opacity coefficients and composition,This uses the same adiabat through the convection zone and distinguishes its inner transition temperature from its outer onset temperature.

At the stellar surface boundary condition, match the interior to a stellar atmosphere, rather than setting the temperature and pressure to zero at an arbitrary radius. At a photosphere of radius $R$, the enclosed mass is $M$, the luminosity is $L$, and the outward radiative flux defines the effective temperature through $F=L/(4\pi R^2)={\sigma }_{\mathrm{SB}}{T}_{\mathrm{eff}}^4$. With no external illumination the incoming specific intensity satisfies $I_{\nu }(0,\mu )=0$ for $\mu <0$. A deeper boundary is matched to the nearly isotropic radiative diffusion field. A grey atmosphere commonly matches near $\tau =2/3$; under the Eddington surface boundary condition this gives $T=T_{\mathrm{eff}}$. For roughly constant gravity and opacity, the total pressure rises by $g\tau /\kappa$ relative to its outer boundary value; if radiation supplies appreciable support, the gas-pressure gradient uses the effective gravity $g-\kappa F/c$ instead.

Take $z$ increasing outward and $\mu$ the outward direction cosine. If $j_{\nu }$ is emission per unit path length and solid angle, the radiative transfer equation isDefine inward-increasing optical depth by $d{\tau }_{\nu }=-\rho {\kappa }_{\nu }dz$ and the source function by $S_{\nu }=j_{\nu }/(\rho {\kappa }_{\nu })$. Then $\mu {\partial }_{{\tau }_{\nu }}{I}_{\nu }=I_{\nu }-S_{\nu }$. This states the emission-coefficient convention; if $j_{\nu }$ is defined per unit mass instead, the emission term is $\rho j_{\nu }$.

For frequency-integrated intensity, introduce the radiation-field momentsThe first angular approximation writes $I\simeq J+3\mu H$. Its odd part does not contribute to $K$, soThis is the Eddington closure approximation, not an exact description of the escaping angular distribution. An approximate surface condition takes the outgoing hemisphere to have direction-independent intensity $I_0$ and the incoming hemisphere to be dark. Its moments give $F=\pi I_0$ and $J(0)=I_0/2=F/(2\pi )$. Integrating the moment equations with the closure then gives $J=3F(\tau +2/3)/(4\pi )$; in thermal equilibrium, $J=a_{\mathrm{rad}}c{T}^4/(4\pi )$ gives $T^4=(3/4)T_{\mathrm{eff}}^4(\tau +2/3)$. This derives the usual approximate photospheric matching condition rather than imposing a zero surface temperature. For conservative coherent isotropic scattering, $S=J$. Integrating the transfer equation gives $dH/d\tau =0$ and $dK/d\tau =H$, the moment conditions for radiative equilibrium.

The exponential ansatz in the PDF has a genuine consistency problem. Compute its moments before making any approximation. For a nonsingular real intensity on all $-1\le \mu \le 1$, take $|ϵ|<1$ and defineThe continuous zero limits are $A_0=1$, $A_1=0$, $A_2=1/3$. Direct integration yieldsFor a decaying nontrivial correction, $\delta >0$ and $\gamma \ne 0$. Exact constant flux then requires $A_1=0$. But $\mathrm{atanh}ϵ/ϵ>1$ for every nonzero real $|ϵ|<1$, so $A_1=0$ forces $ϵ=0$. Thus exact radiative equilibrium cannot imply a small nonzero $ϵ$ for this ansatz.

Substitution into the full transfer equation makes the obstruction sharper. The exponential coefficient must satisfy, for every $\mu$,Comparing the constant and linear coefficients gives $A_0=1$ and $\delta =A_0ϵ=ϵ$. Hence regular real coefficients require $\delta =ϵ=0$, eliminating any decaying correction. This is the obstruction to a single smooth exponential mode in conservative grey transfer. There is no exact nontrivial real solution of the stipulated form with a small nonzero decay parameter. Choosing, for example, $ϵ=\delta =0.1$ and $\gamma \ne 0$ supplies an explicit counterexample to the claimed exact flux constancy, because $A_1\ne 0$.

Several requested conclusions remain useful as formal approximations, provided they are identified as such. If $\delta >0$, the large-depth exponential vanishes, giving $F\to 4\pi \alpha /3$, hence $\alpha =3F/(4\pi )$ for the asymptotic physical flux. It also gives $c{P}_{\mathrm{rad}}\to (4\pi /3)J$. Since $A_0=1+{ϵ}^2/3+O({ϵ}^4)$, neglecting second-order angular errors yields the formal relation $\delta \simeq ϵ$. However $A_1=-ϵ/3+O({ϵ}^3)$, so the finite-depth flux still varies at first order. The formal relation is not an exact transfer or constant-flux solution.

There is also a flux-normalization mismatch in the PDF's last displayed relation. With the same physical $F$ used above, the Hopf function of a grey atmosphere is defined byThe PDF's coefficient $3F/4$ would require redefining its $F$ as physical flux divided by $\pi$; it cannot simultaneously use the earlier $\alpha =3F/(4\pi )$ for the same physical flux. With the consistent physical-flux normalization, formal moment matching givesIn the intended small-$ϵ$ approximation this isThe approximate outgoing surface ansatz then gives the limb darkening ratioFor finite $ϵ$, $\gamma =-0.14/A_0$ gives a parameter-dependent ratio, not a unique exact value. The two Hopf endpoint values alone do not determine the exact emergent intensity: the formal solution of the radiative transfer equation uses the entire source profile, $I(0,\mu )={\int }_0^{\infty }J(t)e^{-t/\mu }dt/\mu$ for $\mu >0$. Nor does the simple linear-plus-exponential ansatz satisfy the exact no-incoming-radiation condition for every negative $\mu$. Thus $2.75$ is a clearly labelled approximate disk-centre-to-limb result, not an exact consequence of inconsistent premises.

The Stefan–Boltzmann law gives $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$ for bolometric luminosity. HenceWith the central values above, this is approximately $1.51$ for the first geometry and $0.67$ for the second. If only one-band luminosities are measured, bolometric corrections are needed before using this formula for effective temperatures.

Two direct methods are fitting the continuum stellar spectrum or several photometric colors to a stellar atmosphere model, and measuring temperature-sensitive absorption lines or line ratios in a stellar spectrum. The former includes the approximate black-body radiation estimate ${\lambda }_{\max }T=b$ from Wien's displacement law; the latter uses ionization and excitation populations. Observations during total eclipse or spectral disentangling help separate the two stellar spectra.

For comparable composition and accessible fuel fractions, the main-sequence lifetime is proportional to available nuclear fuel divided by luminosity, $t_{\mathrm{MS}}\propto M/L$. Estimate the masses from a calibrated mass-luminosity relation $L\propto M^{\alpha }$ in the appropriate mass range. ThenFor illustration $\alpha =3.5$ gives about $1.27$ and $5.6$ for the two graphical luminosity ratios. The common main sequence approximation is additional physics; eclipse depths alone do not determine masses or ages. A main sequence consistency check would also favor the second geometry because ordinary main-sequence stars have radius and effective temperature increasing with mass.

A narrow spectral feature records enhanced emission or absorption near a transition frequency. In a stellar atmosphere, its strength depends on atomic populations, opacity, velocity fields and the source function. Lines therefore diagnose composition as well as the physical atmosphere in which they form.

A stellar atmosphere is the outer region that determines the escaping spectrum and matches the interior stellar structure to its radiation boundary conditions. A photosphere is the layer near unit optical depth, not a zero-temperature boundary. The inward intensity at the outer surface vanishes when there is no external illumination.

A stellar abundance can be expressed as an elemental or isotopic mass fraction, or as a number ratio to hydrogen. Inferring it from a spectral line requires a model for level populations, opacity, radiation and velocity fields. Agreement between suitable ionization stages and weak photospheric lines tests whether the adopted stellar atmosphere is adequate.

Stellar parallax is the apparent change in a star's viewing direction caused by the observer's changing position. Nearby stars shift more against a distant reference background. In an ideal circular Earth orbit, the conventional annual stellar parallax $p$ uses a baseline of one astronomical unit, approximately the Earth-Sun distance. Two opposite observing positions, about six months apart, are separated by approximately two astronomical units. Their maximum angular separation is $2p$, so $p$ is half that excursion, not the full January-to-July shift. The observed excursion depends on the projected baseline and observing dates. ^[ref].

For the perpendicular baseline in the diagram, a right triangle giveswhere $p$ in the final expression is in radians. Combining the small-angle approximation with the parsec convention givesThus $p=0.05$ arcseconds implies $d\simeq 20$ parsecs, while the ideal six-month excursion is $0.10$ arcseconds. Geometry supplies this distance without assuming a luminosity or a particular stellar atmosphere.

In practical astrometry, many observations separate the annual stellar parallax signal from proper motion, the star's own changing sky position. The Gaia DR3 catalogue reports stellar parallax and its uncertainty in milliarcseconds: $50\mathrm{mas}=0.05$ arcseconds. A catalogue angle is a fitted measurement with measurement error. ^[ref].

For a positive stellar parallax with small standard deviation ${\sigma }_p$, first-order propagation gives ${\sigma }_d/d\simeq {\sigma }_p/p$. Thus $p=50\mathrm{mas}$ with ${\sigma }_p=1\mathrm{mas}$ gives $d\simeq 20$ parsecs with ${\sigma }_d\simeq 0.4$ parsecs, assuming systematic offsets are negligible. This linear estimate requires ${\sigma }_p\ll p$; a reciprocal is nonlinear and does not generally preserve symmetric uncertainty intervals.