3-D Radiative Transfer Modeling of Structured Winds in Massive Hot Stars with Wind3D A. Lobel1, J. A. Toala´2, and R. Blomme1 1 Royal Observatory of Belgium, Ringlaan 3, B-1180, Brussels, Belgium 2 Universidad Nacional Autonoma´ de Mexico,´ Centro de Radioastronom´ıa y Astrof´ısica, Campus Morelia, Michoacan,´ Mexico´ Abstract: We develop 3-D models of the structured winds of massive hot stars with the WIND3D radiative transfer (RT) code. We investigate the physical properties of large-scale structures observed in the wind of the B-type supergiant HD 64760 with detailed line profile fits to Discrete Absorption Components (DACs) and rotational modulations observed with IUE in Si IV λ1395. We develop parameterized input models for WIND3D with large-scale equatorial wind density- and velocity- structures, or so-called ‘Co-rotating Interaction Regions’ (CIRs) and ‘Rotational Modulation Regions’ (RMRs). The parameterized models offer important advantages for high-performance RT calculations over ab-initio hy- drodynamic input models. The acceleration of the input model calculations permits us to simulate and investi- gate a wide variety of physical conditions in the extended winds of massive hot stars. The new modeling method is very flexible for constraining the dynamic and geometric wind properties of RMRs in HD 64760. We compute that the modulations are produced by a regular pattern of radial density enhancements that protrude almost linearly into the equatorial wind. We find that the modulations are caused by narrow ‘spoke-like’ wind regions. We present a hydrodynamic model showing that the linearly shaped radial wind pattern can be caused by mechanical wave action at the base of the stellar wind from the blue supergiant. 1 Introduction Accurate mass-loss rates of massive hot stars determined from quantitative spectroscopy are important for understanding the physical properties of the radiative wind driving mechanism that is influenced arXiv:1012.5221v1 [astro-ph.SR] 23 Dec 2010 by dynamic structures on both large and small length scales in the wind. Rotational modulations and Discrete Absorption Components are important tracers of large-scale structures in the highly supersonic winds of these stars. DACs are recurring absorption features observed in UV resonance lines of many OB-stars. They drift bluewards in the absorption portion of P-Cygni profiles. DACs are caused by spiral-shaped density- and velocity structures winding up in the plane of the equator over several tens of stellar radii (e.g., Cranmer & Owocki 1996). Lobel & Blomme (2008) demonstrated with 3-D RT modeling (combined with hydrodynamic simulations) of the detailed DAC evolution in HD 64760 (B0.5 Ib) that these wind spirals are extended density waves emerging from two bright equatorial spots that rotate five times slower than the stellar surface. Detailed hydrodynamic models of its structured wind with these large density waves (or CIRs) reveal only a very small increase of less than 1 % above the smooth (symmetric) wind mass-loss rate. “The multi-wavelength view of hot, massive stars”; 39th Liege` Int. Astroph. Coll., 12-16 July 2010 Figure 2: The curvature, incidence angle on the Figure 1: Dynamic spectrum of Si IV λ1395 ob- surface, and opening angle of the RMRs shown served during 15.5 d in HD 64760 (right-hand in the right-hand panels determine the accelera- panel), compared to 3-D radiative transfer cal- tion (upper panels), inclination (middle panels), culations (left-hand panel) with a parameterized and duration (lower panels), respectively, of the structured wind model. Horizontal absorptions RT calculated modulations (left-hand panels). are the modulations we model in this paper. The bright surface spots produce the large-scale CIRs in the wind. The density enhancements and velocity plateaus of the CIRs yield migrating DACs with a recurrence period of 10.3 d in the UV line profiles of HD 64760. The modulations, on the other hand, show much shorter periods of ∼1.2 d and reveal a time-evolution that substantially differs from the rather slowly shifting DACs. The modula- tions are nearly-flat absorption features observed for only ∼0.5 d to 0.75 d with radial velocities that −1 −1 can range from ∼0 km s to ∼v1=1560 km s (Massa et al. 1995; Prinja 1998). They can intersect the slowly drifting DACs and sometimes reveal a remarkable bow shape (Fullerton et al. 2006) with broad flux minima around ∼930 km s−1. In this paper we present a semi-empiric model of the large- scale wind structures based on detailed RT fits to the time-evolution of the modulations observed in HD 64760. We utilize the WIND3D code for performing non-LTE RT calculations (Sect. 2) of im- portant diagnostic spectral lines, such as Si IV λ1395. We develop parameterized 3-D input models for Wind3D in Sect. 3 because they offer important advantages for high-performance RT calculations over ab-initio hydrodynamic input models. The acceleration of the input model calculations permits us to model and investigate a much broader range of 3-D physical conditions in the wind. 2 Parametrization of Wind3D Input Models Wind3D computes the 3-D non-LTE transport of radiation for the 2-level atom in optically thick resonance lines formed in scattering-dominated extended stellar winds. It lambda-iterates the line source function on a Cartesian grid of 713 points. The radiative transfer equation is solved in parallel over a mesh of 7013 voxels using the converged line source function. Typical computation times with WIND3D are ∼5 h for iterating the line source function, and ∼12 h for calculating the dynamic line profile variability at 90 viewing angles using 16 CPUs. We develop a new software module in WIND3D for semi-empiric modeling of large-scale equa- torial wind density- and velocity-structures. The computer code integrates the time-independent mo- mentum balance equation of radiation-driven rotating winds following Castor, Abbott & Klein (CAK, 1975). First, we parametrize the large-scale density structures in a stationary model of the equatorial wind with: q 0 0 0R − (x − x )2 + y2 111 m ? 0 ρ(x; y; ) = ρ0(r) @1 + A(r) sin @ n + 2π n fc @ AAA ; (1) R? − Rmod 2 2 2 where r =x +y , Rcr ≤ r ≤ Rmod, −π ≤ ≤ π, and Rmod is the outer radius of the wind model. x0=R?sin( 0), where 0 is the incidence angle of the equatorial wind structures on the stellar surface, and 0 ≤ 0 ≤ π/2. ρ0(r) is the density of the smooth (unperturbed) equatorial wind. n denotes the number of large-scale wind structures with enhanced density in one hemisphere, and m is an even exponent of the sine function that sets the opening angle or the tangential width of the wind structures. In wind regions where the sine function vanishes the local wind density equals the smooth wind density. The term in Eq. (1) with fc, where 0 ≤ fc ≤ 1, determines the curvature of the wind structures between R? and Rmod. If fc=1, they turn completely around the star over 2 π from the surface radius R? to Rmod. In case fc=0, the structures do not curve at all and stay strictly radial (linear) in the wind, ordered in equal sectors around the star. We parametrize the density constrast of the large-scale wind structures with the function A in Eq. (1): r 2 r 2 ! A(r) = c1 exp − − c2 + c3 + c4 ; (2) R? R? where c1;2;3;4 are constants. They determine the detailed density profile compared to ρ0(r). We vary the four constants until the parameterized density profile best fits the density contrast of the CIRs in the hydrodynamic wind model of HD 64760. Next, with the parameterized radial wind density-structure we compute the corresponding radial wind velocity-structure for the α parameter of CAK theory set equal to 1/2 in hot stars with 20 kK ≤ Teff ≤ 30 kK. We solve the CAK momentum balance equation in 1-D for many angles , and include the centrifugal force in the equatorial plane due to the photospheric rotation velocity vrot. With the conservation of angular momentum the equatorial tangential wind velocity is vrotR?=r, which becomes negligibly small in the highly supersonic outer regions of the parameterized wind models. We integrate the momentum equation together with the mass continuity equation towards the surface from a typical outer (reference) radius rref ≥30 R? with the boundary condition v(r = rref ) = v1, down to the critical point Rcr of the wind. The value of the mass-loss rate for calculating the CAK line force is computed with Eq. (5) of Lobel & Blomme (2008), but which typically increases the structured wind mass-loss rate by less than ∼1 % in the best-fit 2-spot model. We find that the radial wind velocity-structure computed with the CAK wind-momentum equation integration method, using parameterized stationary wind density models, compares very closely to the hydrodynamic wind velocity model of HD 64760 (Lobel & Toala´ 2009). 3 Radiative Transfer Modeling of Rotational Modulations Figure 1 compares the dynamic spectrum of Si IV observed over 15.5 d in HD 64760 (right-hand panel) and computed with WIND3D (left-hand panel) using a parameterized structured wind model. The slowly bluewards migrating upper and lower DACs result from two spiral-shaped CIRs shown in the right-hand panels of Fig. 2. The nearly horizontal absorption in the modulations is due to Figure 3: The upper panels show two modula- Figure 4: Parameterized model of the equatorial tions computed with Wind3D (left-hand panel) wind-density structure of HD 64760 computed and observed over ∼3.1 d (right-hand panel) in with 3-D RT fits to the modulations in Fig.