SciELO - Scientific Electronic Library Online

vol.55 número1Extensive photometry of V1838 Aql during the 2013 superoutburstOptimal trajectories to Kuiper Belt Objects índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados




Links relacionados

  • No hay artículos similaresSimilares en SciELO


Revista mexicana de astronomía y astrofísica

versión impresa ISSN 0185-1101

Rev. mex. astron. astrofis vol.55 no.1 Ciudad de México abr. 2019



Asymmetric shapes of radio recombination lines from ionized stellar winds

R. Ignace1 

1 Department of Physics & Astronomy, East Tennessee State University, USA.


Recombination line profile shapes are derived for ionized spherical stellar winds at radio wavelengths. It is assumed that the wind is optically thick owing to free-free opacity. Emission lines of arbitrary optical depth are obtained assuming that the free-free photosphere forms in the outer, constant expansion portion of the wind. Previous works have derived analytic results for isothermal winds when the line and continuum source functions are equal. Here, semi-analytic results are derived for unequal source functions to reveal that line shapes can be asymmetric about line center. A parameter study is presented and applications discussed.

Key Words: line; profiles; radiative transfer; stars; early-type; stars; winds; outflows; radio lines; stars


Se calcula la forma de los perfiles en radiofrecuencia de las líneas de recombinación para vientos estelares esféricos ionizados. Se supone que el viento es óptica-mente grueso debido a la opacidad libre-libre. Se obtienen líneas de emisión para profundidades ópticas arbitrarias, suponiendo que la fotósfera libre-libre se forma en la parte externa del viento, la cual se expande de manera constante. En trabajos anteriores se habían obtenido resultados analíticos para vientos isotérmicos, en los que las funciones fuente para las líneas y el continuo se suponían iguales. En este artículo obtenemos resultados semi-analíticos cundo las funciones fuente no son iguales. Las líneas resultantes pueden mostrar formas asimétricas. Se presenta un estudio de los parámetros y se discuten algunas aplicaciones.

1. Introduction

Radio astronomy has long proven to be an important window into the study of stellar astrophysics, and stellar outflows have been no exception (e.g., Dulk 1995; Güdel 2002; Kurt et al. 2002). For stellar winds a key driver has been the prospect of measuring wind mass-loss rates, M˙, from the excess infrared (IR) and radio continuum emission relative to the stellar atmosphere (e.g., Panagia & Felli 1975; Wright & Barlow 1975). Numerous studies have focused on determining M values based on this approach (e.g., Abbott et al. 1980; Abbot, Bieging, & Churchwell 1981; Abbott et al. 1986; Bieging, Abbott, & Churchwell 1989; Leitherer, Chapman, & Koribalski 1995).

One of the main results from a consideration of free-free excesses formed in the wind is that the spectral energy distribution (SED) at long wavelengths will have a power-law slope with flux fνλ-0.6. However, this outcome depends on several assumptions: isothermal, spherical symmetry, large optical depth, negligible contribution from the stellar atmosphere, and constant outflow speed. Cassinelli & Hartmann (1977) explored the effects of different power laws for the wind density and temperature distributions to relate the SED power-law slope to these influences. Schmid-Burgk (1982) showed that such SED slopes persist even for axisymmetric stellar envelopes, as long as the same power-law relations are adopted. The main difference is that flux levels are modified, which would have implications for inferring M˙ values.

Of greater relevance in recent decades has been the abundance of evidence for clumping in massive star winds. In this regard the literature is voluminous, and there has even been a conference to focus on the topic (Hamann, Feldmeier, & Oskinova 2008). The line-driven winds of massive stars (Castor, Abbot, & Klein 1975; Friend & Abbott 1986; Pauldrach, Puls, & Kudritzski 1986) are known to be subject to an instability (e.g., Lucy & White 1980; Owocki, Castor, & Rybicki 1988). This instability produces shocks in the flow and is a natural culprit for stochastic wind clumping. The clumping is well-known to affect the long wavelength emission because of the density-square dependence of the free-free emissivity. In the presence of clumping, the radio emission is overly bright for a given value of M˙ as compared to a smooth (i.e., unclumped) wind with the same mass loss. Neglecting the clumping leads to overestimates of M˙, scaling as the square root of the clumping factor, or inverse to the square root of the volume filling factor of clumps. These factors will be defined precisely in the following section.

Clumping affects any density-square emissivity, including recombination lines. Clumping has been incorporated into several detailed complex numerical codes for modeling massive star atmospheres and their winds, such as CMFGEN (Hillier & Miller 1999) and PoWR (Hamann, Gräfener, & Liermann 2006). An important distinction for clumping is between “macroclumping” and “microclumping”. The former leads to modifications of observables that can depend on the shape of the clump and is sometimes synonymous with a “porosity” treatment. The latter occurs when clumps are all optically thin, so that the radiative transfer does not depend on details of clump morphology. Consequently, microclumping can be handled in terms of a scale parameter, and in fact does not alter the SED slope relative to an unclumped wind (Nugis, Crowther, & Willis 1998). Ignace (2016a) considered the impact of macroclumping vs microclumping for ionized winds at long wavelengths.

This contribution is concerned with modeling a radio recombination line (RRL) profile shape that also includes continuum free-free opacity. The problem has been addressed many times before. Rodríguez (1982) derived the line profile shape for this case, with the interest of supplementing the use of the continuum to obtain M with line broadening formed in the same spatial locale to obtain the wind terminal speed v. Hillier, Jones, & Hyland (1983) did so as well. Ignace (2009) repeated the derivation, and expanded the consideration to include line blends. All of these treatments assume that the source function for the line and continuum is the same, as given by the Planck function for an isothermal wind. Using a numerical radiative transfer calculation, Viner, Vallee, & Hughes (1979) showed that an asymmetric line shape can result when the line and continuum source functions are unequal. Here, this result is explored further through analytic derivations. § 2 introduces the model assumptions and presents a derivation for the line shape. Unlike most previous treatments, the derivation also allows for a power-law distribution of microclumping in the wind. § 3 provides for a parameter study for line profile shapes. § 4 discusses relevant applications for various astrophysical sources.

2. Radio recombination line modeling

Various authors have addressed the relevance of non-LTE effects for interpreting observed RRLs. A discussion of progress on the topic can be found in Gordon & Sorochenko (2002). Relevant to wind-broadened emission lines, Viner et al. undertook a calculation of departure coefficients for studies of H II regions. As previously noted, they allowed for spherical outflow and found that line shapes can be asymmetric. Peters, Longmore, & Dullemond (2012) conducted a similar study for H II regions, and elaborated further on line asymmetry for an outflow. However, neither Viner et al. nor Peters et al. explored the possibility of analytic solutions for the radiative transfer. Here, the approach largely follows Ignace (2009), but relaxing the assumption that the line and continuum source functions are equal. The primary assumptions of the model are as follows:

  1. The wind is spherically symmetric in time average.

  2. The wind is optically thick to free-free opacity. The line can be thin or thick.

  3. While the line and continuum source functions may not be equal, they are taken as constant with radius.

  4. Microclumping is included in the treatment, specifically as a power-law distribution1 with radius. Clumping in massive star winds is both predicted and measured to vary with radius (e.g., Runacres & Owocki 2002; Blomme et al. 2002, 2003; Puls et al. 2006).

2.1. Wind Parameters

Spherical symmetry requires the wind to have a strictly radial wind velocity and density. Being optically thick to free-free opacity, only the large radius flow at constant expansion will be considered. The wind terminal speed is represented by v. The wind also has a time-average mass-loss rate of M˙, and the star has radius R.

Microclumping is represented as

ρ2=Dcl ρ2, (1)

where x represents spatial averaging. On the right-hand side, the average density is given by the smooth wind relation for spherical symmetry, with

ρ=M˙4πR*2vR*r2ρ0R*r2. (2)

For emissivity jνρ2, the emission is enhanced above the smooth wind by the clumping factor D cl. It is common to represent the clumping in the wind with a volume filling factor, fV=Dcl-1. Both approaches are used in the literature (c.f., clumping factor: Hamann & Koesterke 1998 or Ignace, Quigley, & Cassinelli 2003; volume filling factor: Abbott et al. 1981 or Dessart et al. 2000).

For this study the clumping factor is allowed to vary with radius as a power law, with

Dcl r-m. (3)

The case of m = 0 is for clumping that is constant throughout the flow; m > 0 implies that clumping declines with radius; m < 0 is the opposite case. (Note that some care must be taken with use of the power law for clumping, since D c ≥ 1.)

2.2. Line and Continuum Opacities

The free-free opacity, κν, is given by

kv ρ= Kff ni ne, (4)

where ni=ρ/μimH is the number density of ions, with μi the mean molecular weight per free ion; ne=ρ/μemH is the number density of electrons, with μe the mean molecular weight per free electron; m H is the mass of hydrogen; and (Cox 2001)

Kff=3.692×108(1-e-hv/kTc)Zi2gvTC-1/2v-3. (5)

In the preceding equation, h is Planck’s constant, k is the Boltzmann constant, T C is the temperature of the gas appropriate for the continuum emission, Z i is the root mean square ion charge, ν is the frequency, and g ν is the Gaunt factor.

Figure 1 shows the geometry for evaluating the optical depth τ along a ray. Cylindrical coordinates for the observer are (p, α, z), with the observer located at great distance along the +z-axis. Spherical observer coordinates are (r, θ, α), with r 2 = p 2 + z 2. The continuum optical depth along a ray of fixed impact parameter, p, is

τC=τCλr~ρ~2Dclr~ dz~, (6)

Fig. 1 Coordinate definitions used in the derivation for the line and continuum emission from a spherical ionized wind. See text for explanation. The color figure can be viewed online. 

where x~ signifies a normalized parameter, in this case ρ~=ρ/ρ0 and lengths are relative to R, and the optical depth scaling is

τC=KffR*ρ02μi μe mH2. (7)

At long wavelengths that are the focus of this paper, τC gv λ2 for the Rayleigh-Jeans limit, and gvλ0.1.

The line opacity is somewhat similar to that of the continuum in the sense that there is a dependence on the square of the density for recombination. Assuming that the wind speed is highly supersonic, the line optical depth can be approximated from Sobolev theory (Sobolev 1960). The line optical depth becomes

τL=κLρλ(v/r)(1-μ2), (8)

where κLρ Dcl ρ2 F (TL), for F(T L ) a function of temperature appropriate for the line emission, and μ=cos θ.

For the case of constant expansion of the wind at v, the line-of-sight velocity shift due to the Doppler effect is vz=-v μ µ. It is convenient to introduce a normalized velocity shift with

wz=vz/v=-μ. (9)

Also note that p=r sin θ. Then the line optical depth becomes

τL=τLp~-3-m (sin θ)1+m, (10)

where the power-law dependence of D cl with radius has been substituted into the expression, along with ρ̃2=r̃-4, and τL is the optical depth scaling for the line. Casting the line optical depth in terms of p and θ will prove useful for solving the radiative transfer problems in the following sections.

2.3. Solution for the Case of S L = S C

When the line and continuum source functions are the same, let S0=SL=SC. At wavelength λ, just outside the maximum velocity shift of the line, the flux of continuum emission is given by

fC=2πR*2d20Iνp̃dp̃, (11)

where Iν is the emergent intensity as given by

Iν=S01-e-τC(p̃). (12)

When the wind is optically thick, such that the excess emission from the wind greatly exceeds the attentuated stellar emisison through the wind, the radiative transfer has a well-known solution when there is no clumping (Panagia & Felli 1975; Wright & Barlow 1975). When constant clumping is present (m=0), the spectral energy distribution is unchanged, and the flux is simply enhanced above that of a smooth wind (Nugis et al. 1998).

Ignace (2009) also showed that an analytic solution can result with a power-law distribution in the clumping. The following integral relation will be found of general use in subsequent steps:

01-e-axβxdx=1βΓ2βa2/β, (13)

where Γ is the Gamma-function.

For the case at hand, the continuum optical depth is

τCp̃=-+τCdz~r~4+m (14)

=τCp~-3-m0π(sin θ)2+m dθ (15)

=τCp~-3-mGmπ, (16)

where the second line above uses a change of variable to θ, with tan θ=p/z, and

Gmθ=0θ(sin x)2+m dx (17)

The flux of continuum emission becomes

fC=2πR*2d2S013+mΓ23+m×[GmπτCλ]2/(3+m). (18)

In the Rayleigh-Jeans limit, the continuum flux will have a power-law slope of -2+4.2/(3+m), with S0 scaling as λ-2 for the Planck function. When m=0, the canonical slope of -0.6 results. Formally, the analytic solution of equation (18) requires that m>-1.

Within the line, the solution is really no more complicated. Again, Ignace (2009) showed that

fwz=2πR*2d2S00{1-exp[-(τCGm(π)+τL sin θ)p~-3-m]} p~dp~ (19)

While the argument of the exponential now has two terms, the form of the integral is just like that of the pure continuum. The analytic result is

fwz=2πR*2d2S013+mΓ23+m×[GmπτCλ+τL sin θ]2/(3+m). (20)

Note that sin θ=1-wz2. When m=0, the result of Rodríguez (1982) is recovered. The foremost outcome for equal line and continuum source functions is that regardless of the value of m, the line profile is always symmetric about line center. However, when the two source functions are not equal, the line shape will be asymmetric, as demonstrated in the next section.

2.4. Solution for the Case of SL ≠ SC

In the previous section, a relatively complicated radiative transfer problem for line and continuum was found to have an analytic solution. The simplifications required to obtain that solution were spherical symmetry, time-independent flow, and the limit of constant wind expansion. Variation in the clumping factor could be included if the variation could be treated as a power law. Especially important was that both the free-free and line opacities scaled as the square of density.

The final key assumption was that the line and continuum source functions were equal. However, this assumption can be relaxed to allow for unequal source functions (yet still constant throughout the flow at large radius). In this case the solution for the emergent intensity is more complicated, and becomes

Ivp~,wz=SCe-τW-e-τCe-τL+SC1-e-τW+SL1-e-τLe-τW. (21)

This expression has three terms. A ray at impact parameter p̃ intersects the conical isovelocity zone in the form of a ring. Considering just one point on this ring, corresponding to z̃ for a given velocity shift w z, we have two path segments and one point to consider for the accumulation of sinks and sources that contribute to the emergent intensity. The first term in the expression is for the continuum emission up to the point of interest, and then its attenuation by the line opacity at the point. The second term is for the continuum emission from the point of interest to the observer. The third term is the contribution by the line emission, as attenuated by the foreground continuum opacity. Thus as before, τC is the total continuum optical depth along the ray, and we also have τW as the continuum optical depth from the observer to the point of interest where the line emissivity contributes to the emission.

When SL=SC, terms involving τW cancel out. With unequal source functions, the dependence on τW persists. The emergent intensity now becomes:

Iν(p̃)=SC1-e-τC-τL-(SL-SC)e-τW-τL+(SL-SC)e-τW. (22)

For this expression the first term closely mimics the result from the preceding section when SL=SC. Thus the other two terms in the arrangement of equation (22) represent modifications when the source functions are unequal.

The flux still has an analytic solution. However, an additional standard integral relation is required, of the form

0x-βe-axdx=Γ(1-β)aβ-1. (23)

This relation can be applied to the solution for the flux by allowing x=p̃-3-m, for which p̃=x-1/(3+m). One also has pdp=-(3+m)x-βdx with β=(5+m)/(3+m).

The flux in the continuum, outside the velocities of the line, is the same as in the preceding section. However, within the line, the flux now becomes


[GmπτCλ+τL (sin θ)1+m]2/(3+m)-


[GmθτCλ+τL (sin θ)1+m]2/(3+m)+

δLCΓ-23+m[GmθτCλ]2/(3+m), (24)


δLC=SLSC-1, (25)


f0=2πR*2d2SC. (26)

It is frequently the case that line profile data are plotted as continuum normalized. The continuum-normalized emission line profile is given by

f(wz)fC=[1+tLCGmπsin θ)1+m+δLCγm×

GmθGmπ+tLCGmπsinθ1+m2/(3+m)-GmθGmπ2/(3+m), (27)


tLCλ=τL/τC, (28)


γm=-(3+m)Γ[-2/(3+m)]Γ[+2/(3+m)], (29)

where the negative anticipates that the Gamma function in the numerator is also negative.

Note the following special cases. Of course, where there is no line opacity, the radio SED will be a power law in wavelength with

fvBvτc2/(3+m)gv2/(3+m)λ-2(1+m)/(3+m), (30)

for SC = Bν (TC ). The opposite extreme occurs when the line opacity is significant, but the continuum is negligible. The emission line profile shape becomes

f(wz)SL(τL)2/(3+m) (1-wz2)1/(3+m). (31)

Note that the line shape is symmetric about line center in the limit of a strong line. Wavelength dependence pertinent to the specific line transition is implied through the factors SL and TL.

Equation (27) is the main result of this study. The first term represents a symmetric component to the emission line profile. The subsequent two terms contribute generally to asymmetric influences to the line in the form of G m(θ). These influences depend on the clumping power-law exponent m, on the ratio of the source functions S L /S C , and on the ratio of optical depths T L /T C . Note that if δLC>0, the line is in emission, whereas for δLC<0, the line is in absorption. Illustrative examples are given in the following section.

3. Results

Figures 2-4 provide illustrative results for line profile shapes. Figure 2 shows results for m = −0.5 (i.e., clumping that increases with radius); Figure 3 for m = 0 (i.e., clumping that is constant with radius); and Figure 4 for m = +1 (i.e., clumping that declines with radius). Each figure has 4 panels: upper left is for δ LC = −0.15, lower left is for δ LC = +0.23, upper right is for δ LC = +0.62, and lower right is for δ LC = +1.0 for Figures 2 and 3, but δ LC = −0.1, 0.4, 0.9, and 1.4 for Figure 4. Each panel has 5 line profiles, with tLC = 0.32, 0.56, 1.0, 1.8, and 3.2, from the weakest line to the strongest. The profiles are continuum normalized and plotted against velocity shift, wz. Note that each panel has a different ordinate scale.

Fig. 2 Continuum normalized emission line profiles for the case of m = −0.5. The 4 panels are: upper left for δLC = −0.15, lower left for δLC = +0.23, upper right for δLC = +0.62, and lower right for δLC = +1.0. In each panel, the 5 line profiles are shown for tLC = 0.32 (dot-dashed), 0.56 (long dashed), 1.0 (short dashed), 1.8 (dotted), and 3.2 (solid), from the weakest line to the strongest. Lines are plotted against normalized velocity shift. Note that each panel has a different scale. 

Fig. 3 As in Figure 2, except for m = 0. 

Fig. 4 As in Figure 2, except for m = 1 and with δLC = −0.1, 0.4, 0.9, and 1.4. 

When δLC<0, the line profile is actually in absorption. For m=+1, the line shape takes the appearance of a weak P Cygni line shape, with blueshifted net absorption and redshifted net emission.

For m=-0.5, the line profiles are more symmetric about line center as compared to either m=0.0 or m=+1.0. As m approaches -1, the line profile becomes perfectly symmetric, because the factor contributing to the line asymmetry cancels exactly when 2/(3+m)=1. As m increases, the line shapes become increasingly asymmetric with the line skewed preferentially toward blueshifted velocities. This is natural generally because the attentuation of line emission from the far-side hemisphere of the wind is greater than from the near-side hemisphere. When SL=SC, absorption is exactly compensated by emission, and no asymmetry in the line can result. For the given assumptions, the line asymmetry occurs only when the source functions are unequal.

The line profiles display some degeneracy between δ LC and t LC . As tLC becomes large, asymmetry in the line shape lessens, in the sense that the peak emission shifts closer to line center. A large line optical depth means that a (positive) δ LC has less influence on the line shape. Generally, t LC controls the degree of asymmetry in the line, and δ LC acts as an overall amplitude for the line emission (or line equivalent width).

4. Conclusions

The focus of this contribution has been to high-light the asymmetry of RRLs arising from a spherical wind using an analytic derivation. Previous analytic work produced symmetric line shapes. Numerical calculations have demonstrated that asymmetric lines can be produced. Here, with the assumption of constant but unequal line and continuum source functions, asymmetric line shapes are produced. The derivation allows for the presence of microclumping in the wind in terms of a power-law distribution (rising or declining with radius from the star). While the clumping distribution can impact line asymmetry, the line asymmetry results even with constant clumping, or no clumping whatsoever. Under the model assumptions, emission line asymmetry arises from the continuum opacity (specifically the appearance of the term with G m(θ) in equation 27) that absorbs the redshifted emission from the far hemisphere more than the blueshifted emission from the near hemisphere. The result is a line shape with blueshifted emission peak. (The same effect arises in X-ray lines; c.f., Ignace 2016b).

Fig. 5 The line profiles for the case m = 0.0 shown to highlight the evolution of line asymmetry as a function of tLC . The line profiles have been continuum subtracted and normalized to a peak emission of unity. The parame-ter δLC = 0.5 was held fixed. The line profiles are shown for log tLC = −0.5 (red), −0.125 (blue), +0.25 (green), +0.625 (magenta), and +1.0 (black). The color figure can be viewed online. 

RRLs are vigorously pursued as a diagnostic of source properties, from kinematics to geometrical aspects. Peters, Longmore, & Dullemond (2012) have made an in-depth study of various factors that affect the flux of line emission and the shape of the line profile, including line asymmetry. Observational motivation for understanding line asymmetry of RRLs includes some objects as the early-type binary MWC349, specifically the H76α line (Escalante et al. 1989). Understanding the line formation is key for distinguishing between radiative transfer effects for the line formation versus the influence of aspherical effects intrinsic to the source, such as binarity, or asymmetric mass-loss, or flow geometry. Applications for such effects include emission-line objects like MWC349 (e.g., LkHα 101; Thum et al. 2013), outflow from star-forming clumps (e.g., Kim et al. 2018), and planetary nebulae (e.g., Ershov & Berulis 1989; Sánchez Contreras et al. 2017). Analytic solutions are valuable to these studies in two main respects: (a) they allow for rapid evaluation of parameter space that can be honed with more detailed numerical calculations to fit the data, and (b) they are important for providing non-trivial benchmarks against which numerical codes can be tested.

The author expresses appreciation to an anony-mous referee for several helpful comments.


Abbott, D. C., Bieging, J. H., Churchwell, E., Cassinelli, J. P. 1980, ApJ, 238, 196 [ Links ]

Abbott, D. C., Bieging, J. H., & Churchwell E. B. 1981, ApJ, 250, 645 [ Links ]

Abbott, D. C., Bieging, J. H., Churchwell, E. B., & Torres, A. 1986, ApJ, 280, 671 [ Links ]

Bieging, J., Abbott, D., & Churchwell, E., 1989, ApJ, 340, 518 [ Links ]

Blomme, R., Prinja, R., Runacres, M., Colley, S. 2002, A&A, 382, 921 [ Links ]

Blomme, R., van de Steene, G., Prinja, R., Runacres, M., & Clark, J., 2003, A&A, 408, 715 [ Links ]

Cassinelli, J. P. & Hartmann, L. 1977, ApJ, 212, 488 [ Links ]

Castor, J., Abbott, D., & Klein, R. 1975, ApJ, 195, 157 [ Links ]

Cox, A. (ed.), 2000, Allen’s Astrophysical Quantities, (New York: AIP Press, Spring) [ Links ]

Dessart, L., Crowther, P., Hillier, D., et al. 2000, MNRAS, 315, 407 [ Links ]

Dulk, G. A. 1985, ARA&A, 23, 169 [ Links ]

Escalante, V., Rodríguez, L. F., Moran, J., & Cantó, J. 1989, RMxAA, 17, 11 [ Links ]

Friend, D. & Abbott, D. 1986, ApJ, 311, 701 [ Links ]

Gordon, M. & Sorochenko, R. 2002, Radio Recombination Lines. Their Physics and Astronomical Applications, (Kluwer: Dordrecht) [ Links ]

Güdel, M. 2002, ARA&A, 40, 217 [ Links ]

Hamann, W.-R. & Koesterke, L. 1998, A&A, 335, 1003 [ Links ]

Hamann, W.-R., Feldmeier, A., & Oskinova, L. (eds.) 2008, Clumping in Hot-Star Winds (Potsdam: Universitätsverlag Potsdam) [ Links ]

Hamann, W.-R., Gräfener, G., & Liermann, A. 2006, A&A, 457, 1015 [ Links ]

Hillier, D. I., Jones, T. J., & Hyland, A. R. 1983, ApJ, 271, 221 [ Links ]

Hillier, D. I. & Miller, D. L. 1999, ApJ, 519, 354 [ Links ]

Ignace, R. 2009, AN, 330, 717 [ Links ]

Ignace, R. 2016a, MNRAS, 457, 4123 [ Links ]

Ignace, R. 2016b, AdSpR, 58, 694 [ Links ]

Kim, W.-J., Urquhart, J., Wyrowski, F., Menten, K., & Csengeri, T. 2018, A&A, 616, A107 [ Links ]

Leitherer, C., Chapman, J. M., & Koribalski, B. 1995, ApJ, 450, 289 [ Links ]

Lucy, L. & White, R. 1980, ApJ, 241, 300 [ Links ]

Nugis, T., Crowther, P. A., & Willis A. J. 1998, A&A, 333, 956 [ Links ]

Owocki, S., Castor, J., & Rybicki, G. 1988, ApJ, 335, 914 [ Links ]

Panagia, N. & Felli, M. 1975, A&A, 39, 1 [ Links ]

Pauldrach, A., Puls, J., & Kudritzki, R. 1986, A&A, 164, 86 [ Links ]

Puls, J., Markova, N., Scuderi, S., et al. 2006, A&A, 454, 625 [ Links ]

Rodríguez, L. F., 1982, RMxAA, 5, 179 [ Links ]

Runacres, M. & Owocki, S. 2002, A&A, 381, 1015 [ Links ]

Sánchez Contreras, C., Báez-Rubio, A., Alcolea, J., Bujarrabal, V., & Martín-Pintado, J. 2017, A&A, 603, A67 [ Links ]

Sobolev, V. 1960, Moving Envelopes of Stars (Cambridge, MA: CUP) [ Links ]

Peters, P., Longmore, S., & Dullemond, C. 2012, MNRAS, 425, 2352 [ Links ]

Thum, C, Neri, R., Báez-Rubio, A., & Krips, M. 2013, A&A, 556, A129 [ Links ]

Weiler, K., Panagia, N., Monte, M., & Sramek, R. 2002, ARA&A, 40, 387 [ Links ]

Wright, A. E. & Barlow, M. J. 1975, MNRAS, 170, 41 [ Links ]

1The additional power-law distribution need not be attributed to clumping. It could be attributed to something else that modifies the density. However, it cannot be the velocity law, since that would lead to a different geometry for the isovelocity zones and would invalidate the derivation that follows. The inclusion of the additional power law follows the spirit of the approach in Cassinelli & Hartmann (1977).

Received: September 10, 2018; Accepted: October 30, 2018

Richard Ignace: Department of Physics & Astronomy, East Tennessee State University, Johnson City, TN 37614, USA (

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License