Thermal emission from bow shocks II. 3D magnetohydrodynamic models of ζ Ophiuchi

The nearby, massive, runaway star zeta Ophiuchi has a large bow shock detected in optical and infrared light and, uniquely among runaway O stars, diffuse X-ray emission detected from the shocked stellar wind. Here we make the first detailed computational investigation of the bow shock of zeta Ophiuchi, to test whether a simple model of the bow shock can explain the observed nebula, and to compare the detected X-ray emission with simulated emission maps. We reanalysed archival Chandra observations of the thermal diffuse X-ray emission from the shocked wind region of the bow shock, finding total unabsorbed X-ray flux in the 0.3-2 keV band corresponding to a diffuse X-ray luminosity of L-X = 2.33(-1.54)(+1.12)( )x 10(29) erg s(-1), consistent with previous work. The diffuse X-ray emission arises from the region between the star and the bow shock. Three-dimensional magnetohydrodyanmic simulations were used to model the interaction of the star's wind with a uniform interstellar medium (ISM) using a range of stellar and ISM parameters motivated by observational constraints. Synthetic infrared, H alpha, soft X-ray, emission measure, and radio 6 GHz emission maps were generated from three simulations, for comparison with the relevant observations. Simulations where the space velocity of zeta Ophiuchi has a significant radial velocity produce infrared emission maps with the opening angle of the bow shock in better agreement with observations than for the case where motion is fully in the plane of the sky. All three simulations presented here have X-ray emission fainter than observed, in contrast to results for NGC 7635. The simulation with the highest pressure has the closest match to X-ray observations, with a flux level within a factor of two of the observational lower limit, and emission weighted temperature of log(10)(T-A/K) = 6.4, although the morphology of the diffuse emission appears somewhat different. The observed X-ray emission is from a filled bubble that is brightest near the star, whereas simulations predict brightening towards the contact discontinuity as density increases. This first numerical study of the bow shock and wind bubble around zeta Ophiuchi uses a relatively simple model of a uniform ISM and ideal-magnetohydrodynamics, and can be used as a basis for comparing results from models incorporating more physical processes, or higher resolution simulations that may show more turbulent mixing.