9 research outputs found
Implementing efficient concerted rotations using Mathematica and C code
In this article we demonstrate a general and efficient metaprogramming implementation of concerted rotations using Mathematica. Concerted rotations allow the movement of a fixed portion of a polymer backbone with fixed bending angles, like a protein, while maintaining the correct geometry of the backbone and the initial and final points of the portion fixed. Our implementation uses Mathematica to generate a C code which is then wrapped in a library by a Python script. The user can modify the Mathematica notebook to generate a set of concerted rotations suited for a particular backbone geometry, without having to write the C code himself. The resulting code is highly optimized, performing on the order of thousands of operations per second
Sub-surface convection zones in hot massive stars and their observable consequences
We study the convection zones in the outer envelope of hot massive stars
which are caused by opacity peaks associated with iron and helium ionization.
We determine the occurrence and properties of these convection zones as
function of the stellar parameters. We then confront our results with
observations of OB stars. A stellar evolution code is used to compute a grid of
massive star models at different metallicities. In these models, the mixing
length theory is used to characterize the envelope convection zones. We find
the iron convection zone (FeCZ) to be more prominent for lower surface gravity,
higher luminosity and higher initial metallicity. It is absent for luminosities
below about 10^{3.2}\Lsun, 10^{3.9}\Lsun, and \Lsun$ for the
Galaxy, LMC and SMC, respectively. We map the strength of the FeCZ on the
Hertzsprung-Russell diagram for three metallicities, and compare this with the
occurrence of observational phenomena in O stars: microturbulence, non-radial
pulsations, wind clumping, and line profile variability. The confirmation of
all three trends for the FeCZ as function of stellar parameters by empirical
microturbulent velocities argues for a physical connection between
sub-photospheric convective motions and small scale stochastic velocities in
the photosphere of O- and B-type stars. We further suggest that clumping in the
inner parts of the winds of OB stars could be caused by the same mechanism, and
that magnetic fields produced in the FeCZ could appear at the surface of OB
stars as diagnosed by discrete absorption components in ultraviolet absorption
lines.Comment: Accepted for publication in Astronomy and Astrophysic