In this paper we describe the FIT\textit{spec} code, a data mining tool for the automatic fitting of synthetic stellar spectra. The program uses a database of 27\,000 {\sc cmfgen} models of stellar atmospheres arranged in a six-dimensional (6D) space, where each dimension corresponds to one model parameter. From these models a library of 2\,835\,000 synthetic spectra were generated covering the ultraviolet, optical, and infrared region of the electromagnetic spectrum. Using FIT\textit{spec} we adjust the effective temperature and the surface gravity. From the 6D array we also get the luminosity, the metallicity, and three parameters for the stellar wind: the terminal velocity ($v_\infty$), the $\beta$ exponent of the velocity law, and the clumping filling factor ($F_{\rm cl}$). Finally, the projected rotational velocity ($v\cdot\sin i$) can be obtained from the library of stellar spectra. Validation of the algorithm was performed by analyzing the spectra of a sample of eight O-type stars taken from the {\sc iacob} spectroscopic survey of Northern Galactic OB stars. The spectral lines used for the adjustment of the analyzed stars are reproduced with good accuracy. In particular, the effective temperatures calculated with the FIT\textit{spec} are in good agreement with those derived from spectral type and other calibrations for the same stars. The stellar luminosities and projected rotational velocities are also in good agreement with previous quantitative spectroscopic analyses in the literature. An important advantage of FIT\textit{spec} over traditional codes is that the time required for spectral analyses is reduced from months to a few hours.

We present a multi-wavelength (X-ray to optical) analysis, based on non-local thermodynamic equilibrium photospheric+wind models, of the B0 Ia-supergiant: $\epsilon$~Ori. The aim is to test the consistency of physical parameters, such as the mass-loss rate and CNO abundances, derived from different spectral bands. The derived mass-loss rate is $\dot{M}/\sqrt{f_\infty}\sim$1.6$\times$10$^{-6}$ M$_\odot$ yr$^{-1}$ where $f_\infty$ is the volume filling factor. However, the S IV $\lambda\lambda$1062,1073 profiles are too strong in the models; to fit the observed profiles it is necessary to use $f_\infty<$0.01. This value is a factor of 5 to 10 lower than inferred from other diagnostics, and implies $\dot{M} \lesssim1 \times 10^{-7}$ M$_\odot$ yr$^{-1}$. The discrepancy could be related to porosity-vorosity effects or a problem with the ionization of sulfur in the wind. To fit the UV profiles of N V and O VI it was necessary to include emission from an interclump medium with a density contrast ($\rho_{cl}/\rho_{ICM}$) of $\sim$100. X-ray emission in H-He like and Fe L lines was modeled using four plasma components located within the wind. We derive plasma temperatures from $1 \times 10^{6}$ to $7\times 10^{6}$ K, with lower temperatures starting in the outer regions (R$_0\sim$3-6 R$_*$), and a hot component starting closer to the star (R$_0\lesssim$2.9 R$_*$). From X-ray line profiles we infer $\dot{M} <\, 4.9\times10^{-7}$ M$_\odot$ yr$^{-1}$. The X-ray spectrum ($\geq$0.1 kev) yields an X-ray luminosity $L_{\rm X}\sim 2.0\times10^{-7} L_{\rm bol}$, consistent with the superion line profiles. X-ray abundances are in agreement with those derived from the UV and optical analysis: $\epsilon$ Ori is slightly enhanced in nitrogen and depleted in carbon and oxygen, evidence for CNO processed material.

We present a method for computing the net transmission of X-rays emitted by shock-heated plasma distributed throughout a partially optically thick stellar wind from a massive star. We find the transmission by an exact integration of the formal solution, assuming that the emitting plasma and absorbing plasma are mixed at a constant mass ratio above some minimum radius, below which there is assumed to be no emission. This model is more realistic than either the slab absorption associated with a corona at the base of the wind or the exospheric approximation that assumes that all observed X-rays are emitted without attenuation from above the radius of optical depth unity. Our model is implemented in XSPEC as a pre-calculated table that can be coupled to a user-defined table of the wavelength dependent wind opacity. We provide a default wind opacity model that is more representative of real wind opacities than the commonly used neutral interstellar medium (ISM) tabulation. Preliminary modeling of \textit{Chandra} grating data indicates that the X-ray hardness trend of OB stars with spectral subtype can largely be understood as a wind absorption effect.

We fit every emission line in the high-resolution Chandra grating spectrum of zeta Pup with an empirical line profile model that accounts for the effects of Doppler broadening and attenuation by the bulk wind. For each of sixteen lines or line complexes that can be reliably measured, we determine a best-fitting fiducial optical depth, tau_* = kappa*Mdot/4{pi}R_{\ast}v_{\infty}, and place confidence limits on this parameter. These sixteen lines include seven that have not previously been reported on in the literature. The extended wavelength range of these lines allows us to infer, for the first time, a clear increase in tau_* with line wavelength, as expected from the wavelength increase of bound-free absorption opacity. The small overall values of tau_*, reflected in the rather modest asymmetry in the line profiles, can moreover all be fit simultaneously by simply assuming a moderate mass-loss rate of 3.5 \pm 0.3 \times 10^{-6} Msun/yr, without any need to invoke porosity effects in the wind. The quoted uncertainty is statistical, but the largest source of uncertainty in the derived mass-loss rate is due to the uncertainty in the elemental abundances of zeta Pup, which affects the continuum opacity of the wind, and which we estimate to be a factor of two. Even so, the mass-loss rate we find is significantly below the most recent smooth-wind H-alpha mass-loss rate determinations for zeta Pup, but is in line with newer determinations that account for small-scale wind clumping. If zeta Pup is representative of other massive stars, these results will have important implications for stellar and galactic evolution.