Simulation of turbulent convection in a slowly rotating red giant star

The first 3-D non-linear hydrodynamical simulation of the inner convective envelope of a rotating low mass red giant star is presented. This simulation, computed with the ASH code, aims at understanding the redistribution of angular momentum and heat in extended convection zones. The convection patterns achieved in the simulation consist of few broad and warm upflows surrounded by a network of cool downflows. This asymmetry between up and downflows leads to a strong downward kinetic energy flux, that must be compensated by an overluminous enthalpy flux in order to carry outward the total luminosity of the star. The influence of rotation on turbulent convection results in the establishment of largescale mean flows: a strong radial differential rotation and a single cell poleward meridional circulation per hemisphere. A detailed analysis of angular momentum redistribution reveals that the meridional circulation transports angular momentum outward in the radial direction and poleward in the latitudinal direction, with the Reynolds stresses acting in the opposite direction. This simulation indicates that the classical hypothesis of mixing length theory and solid-body rotation in the envelope of red giants assumed in 1-D stellar evolution models are unlikely to be realized and thus should be reconsidered. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling photospheric magnetoconvection

The increasing power of computers makes it possible to model the non-linear interaction between magnetic fields and convection at the surfaces of solar-type stars in ever greater detail. We present the results of idealized numerical experiments on two-dimensional magnetoconvection in a fully compressible perfect gas. We first vary the aspect ratio λ of the computational box and show that the system runs through a sequence of convective patterns, and that it is only for a sufficiently wide box (λ ≥ 6) that the flow becomes insensitive to further increases in λ. Next, setting λ = 6, we decrease the field strength from a value strong enough to halt convection and find transitions to small-scale steady convection, next to spatially modulated oscillations (first periodic, then chaotic) and then to a new regime of flux separation, with regions of strong field (where convection is almost completely suppressed) separated by broad convective plumes. We also explore the effects of altering the boundary conditions and show that this sequence of transitions is robust. Finally, we relate these model calculations to recent high-resolution observations of solar magnetoconvection, in plage regions as well as in light bridges and the umbrae of sunspots.     

Convective intensification of magnetic fields in the quiet Sun

Kilogauss-strength magnetic fields are often observed in intergranular lanes at the photosphere in the quiet Sun. Such fields are stronger than the equipartition field B _e, corresponding to a magnetic energy density that matches the kinetic energy density of photospheric convection, and comparable with the field B _p that exerts a magnetic pressure equal to the ambient gas pressure. We present an idealized numerical model of three-dimensional compressible magnetoconvection at the photosphere, for a range of values of the magnetic Reynolds number. In the absence of a magnetic field, the convection is highly supercritical and characterized by a pattern of vigorous, time-dependent, 'granular' motions. When a weak magnetic field is imposed upon the convection, magnetic flux is swept into the convective downflows where it forms localized concentrations. Unless this process is significantly inhibited by magnetic diffusion, the resulting fields are often much greater than B _e and the high magnetic pressure in these flux elements leads to their being partially evacuated. Some of these flux elements contains ultraintense magnetic fields that are significantly greater than B _p. Such fields are contained by a combination of the thermal pressure of the gas and the dynamic pressure of the convective motion, and they are constantly evolving. These ultraintense fields develop owing to non-linear interactions between magnetic fields and convection; they cannot be explained in terms of 'convective collapse' within a thin flux tube that remains in overall pressure equilibrium with its surroundings.     

The structure of the solar convective overshooting zone

Using a non-local theory of convection, we calculated the structure of the solar convection zone, paying special attention to the detailed structure of the lower overshooting zone. Our results show that an extended transition zone exists near the bottom of the convection zone, where the temperature gradient turns smoothly from adiabatic in the convection zone to radiative in solar interior. A super-radiative temperature region is found in the overshooting zone under the solar convection zone, where     ,     ,     and     . The extension of the super-radiative region (defined by     l is about 0.63  H _P (0.053 R_⊙). A careful comparison of the distribution of adiabatic sound speed and density with the local one is carried out. It is found, strikingly, that the distribution of adiabatic sound speed and density of our model is roughly consistent with the results of reversion from solar oscillation observations.     

The boundary layer of VW Hydri in quiescence

In this short paper, we suggest that the missing boundary layer luminosity of dwarf novae in quiescence is released mainly in the ultraviolet (UV) as the second component commonly identified in the far-UV as the 'accretion belt'. We present the well-studied SU UMa-type system VW Hyi in detail as a prototype for such a scenario. We consider detailed multiwavelength observations and in particular the recent Far Ultraviolet Spectroscopic Explorer ( FUSE ) observations of VW Hyi in quiescence, which confirm the presence of a second component (the 'accretion belt') in the far-UV spectrum. The temperature  (≈ 50 000 K)  and rotational velocity  (≈ 3000 km s ⁻¹)  of this second component are entirely consistent with the optically thick region  (τ≈ 1)  located just at the outer edge of optically thin boundary layer in the simulations of Popham. This second component contributes about 20 per cent of the far-UV flux. Using geometrical assumptions and taking into account the X-ray luminosity, we find that the total boundary layer luminosity sums up to   L _BL= (0.53 ± 0.25) L _disc  , while the theory (Kluźniak) predicts, for the rotation rate of VW Hyi's white dwarf,   L _BL≈ (0.76 ± 0.03) L _disc  . About one-fifth of the boundary layer energy is emitted in the X-ray and the remaining is emitted in the UV. This scenario is consistent with the recent simultaneous X-ray and UV observations of VW Hyi by Pandel, Córdova & Howell, from which we deduce here that the viscosity in the boundary layer region must be of the order of  ν≈ 10¹³–10¹⁴ cm² s ⁻¹  , depending on the white dwarf mass and the size of the boundary layer.     

Realistic simulations of stellar surface convection with ANTARES: I. Boundary conditions and model relaxation

We have implemented open boundary conditions into the ANTARES code to increase the realism of our simulations of stellar surface convection. Even though we greatly benefit from the high accuracy of our fifth order numerical scheme (WENO5) the broader stencils needed for the numerical scheme complicate the implementation of boundary conditions. We show that the effective temperature of a numerical simulation cannot be changed by corrections at the lower boundary since the thermal stratification does only change on the Kelvin–Helmholtz time scale. Except for very shallow models this time scale cannot be covered by multidimensional simulations due to the enormous computational requirements. We demonstrate to what extent numerical simulations of stellar surface convection are sensitive to the initial conditions and the boundary conditions. An ill-conceived choice of parameters for the boundary conditions can have a severe impact. Numerical simulations of stellar surface convection will only be (physically) meaningful and realistic if the initial model, the extent and position of the simulation box and the parameters from the boundary conditions are chosen adequately.     

Evolutionary calculations of carbon dredge-up in helium envelope white dwarfs

We investigate the evolution of cooling helium atmosphere white dwarfs using a full evolutionary code, specifically developed to follow the effects of element diffusion and gravitational settling on white dwarf cooling. The major difference between this work and previous work is that we use more recent opacity data from the OPAL project. Since, in general, these opacities are higher than those available 10 years ago, at a given effective temperature, convection zones go deeper than in models with older opacity data. Thus convective dredge-up of observationally detectable carbon in helium atmosphere white dwarfs can occur for thicker helium layers than found by Pelletier et al. We find that the range of observed C to He ratios in different DQ white dwarfs of similar effective temperature is well explained by a range of initial helium layer mass between 10⁻³ and 10⁻² M⊙, in good agreement with stellar evolution theory, assuming a typical white dwarf mass of 0.6 M⊙. We also predict that oxygen will be present in DQ white dwarf atmospheres in detectable amounts if the helium layer mass is near the lower limit compatible with stellar evolution theory. Determination of the oxygen abundance has the potential of providing information on the profile of oxygen in the core and hence on the important ¹²C(α,γ)¹⁶O reaction rate.     

Spatially resolved X-ray spectroscopy of the core of the Centaurus cluster

We present Chandra data from a 31.7-ks observation of the Centaurus cluster, using the ACIS-S detector. Images of the X-ray emission show a plume-like feature at the centre of the cluster, of extent 60 arcsec (20 kpc in projection). The feature has the same metallicity as gas at a similar radius, but is cooler. Using adaptive binning, we generate temperature, abundance and absorption maps of the cluster core. The radial abundance profile shows that the previously known, steep abundance gradient peaks with a metallicity of  1.3–1.8 Z_⊙  at a radius of about 45 arcsec (15 kpc), before falling back to 0.4 Z_⊙ at the centre of the cluster. A radial temperature profile shows that the temperature decreases inwards. We determine the spatial distributions of each of two temperature components, where applicable. The radiative cooling time of the cooler component within the inner 10 arcsec (3 kpc) is less than  2×10⁷ yr  . X-ray holes in the image coincident with the radio lobes are seen, as well as two outer sharp temperature drops, or cold fronts. The origin of the plume is unclear. The existence of the strong abundance gradient is a strong constraint on extensive convection or gas motion driven by a central radio source.     

High‐frequency waves and granulation

The theory of periodical shear flow is applied for the exploration of the effect of solar granulation on highfrequency waves in the solar photosphere. It is shown that upgoing and downgoing waves are trapped in intergranular spaces and granules, respectively. Upgoing waves in fast downflows are unstable. The theory is in a good agreement with observations.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

Non-adiabatic oscillations of red giants

Using our non-local time-dependent theory of convection, the linear non-adiabatic oscillations of 10 evolutionary model series with masses of  1–3 M_⊙  are calculated. The results show that there is a red giant instability strip in the lower temperature side of the Hertzsprung–Russell diagram which goes along the sequences of the red giant branch and the asymptotic giant branch. For red giants of lower luminosities, pulsation instability is found at high order overtones; the lower order modes from the fundamental to the second overtone are stable. Towards higher luminosity and lower effective temperature, instability moves to lower order modes, and the amplitude growth rate of oscillations also grows. At the high luminosity end of the strip, the fundamental and the first overtone become unstable, while all the modes above the fourth order become stable. The excitation mechanisms have been studied in detail. It is found that turbulent pressure plays a key role for excitation of red variables. The frozen convection approximation is unavailable for the low temperature stars with extended convective envelopes. In any case, this approximation can explain neither the red edge of the Cepheid instability strip, nor the blue edge of the pulsating red giant instability strip. An analytic expression of a pulsation constant as a function of stellar mass, luminosity and effective temperature is presented from this work.     

恒星内部的贯穿对流(英文)

贯穿对流是恒星演化理论长期未获很好解决的理论问题。尽管它已有近五十年的研究历史 ,但至今人们对贯穿对流仍存在很多的错误的理解。究其原因是因为人们对贯穿对流的理解依然停留在唯象的混合长理论图象上。根据一种非局部对流理论 ,我们揭示了贯穿对流区的结构。在太阳下部的贯穿对流区 ,温度梯度是亚绝热的 ,但是是超辐射的。本文还对恒星对流理论给出了一个简短的评述。     

X-ray and ultraviolet observations of the dwarf nova VW Hyi in quiescence

We present an analysis of X-ray and ultraviolet (UV) data of the dwarf nova VW Hyi that were obtained with XMM–Newton during the quiescent state. The X-ray spectrum indicates the presence of an optically thin plasma in the boundary layer that cools as it settles on to the white dwarf. The plasma has a continuous temperature distribution that is well described by a power law or a cooling flow model with a maximum temperature of 6–8 keV. We estimate from the X-ray spectrum a boundary layer luminosity of  8 × 10³⁰ erg s^-1  , which is only 20 per cent of the disc luminosity. The rate of accretion on to the white dwarf is  5 × 10⁻¹² M_⊙ yr⁻¹  , about half of the rate in the disc. From the high-resolution X-ray spectra, we estimate that the X-ray emitting part of the boundary layer is rotating with a velocity of 540 km s⁻¹, which is close to the rotation velocity of the white dwarf but is significantly smaller than the Keplerian velocity. We detect a 60-s quasi-periodic oscillation of the X-ray flux, which is likely to be due to the rotation of the boundary layer. The X-ray and the UV flux show strong variability on a time-scale of ∼1500 s. We find that the variability in the two bands is correlated and that the X-ray fluctuations are delayed by ∼100 s. The correlation indicates that the variable UV flux is emitted near the transition region between the disc and the boundary layer and that accretion rate fluctuations in this region are propagated to the X-ray emitting part of the boundary layer within ∼100 s. An orbital modulation of the X-ray flux suggests that the inner accretion disc is tilted with respect to the orbital plane. The elemental abundances in the boundary layer are close to their solar values.     

Widespread Occurrence of Trenching Patterns in the Granulation Field: Evidence for Roll Convection?

Time-averaged series of granulation images are analysed using COLIBRI, a purpose-adapted version of a code originally developed to detect straight or curvilinear features in aerospace images. The image-processing algorithm utilises a nonparametric statistical criterion that identifies a straight-line segment as a linear feature (lineament) if the photospheric brightness at a certain distance from this line on both sides is stochastically lower or higher than at the line itself. Curvilinear features can be detected as chains of lineaments, using a modified criterion. Once the input parameters used by the algorithm are properly adjusted, the algorithm highlights “ridges” and “trenches” in the relief of the brightness field, drawing white and dark lanes. The most remarkable property of the trenching patterns is a nearly universally present parallelism of ridges and trenches. Since the material upflows are brighter than the downflows, the alternating, parallel light and dark lanes should reflect the presence of roll convection in the subphotospheric layers. If the numerous images processed by us are representative, the patterns revealed suggest a widespread occurrence of roll convection in the outer solar convection zone. In particular, the roll systems could form the fine structure of larger scale, supergranular and/or mesogranular convection flows. Granules appear to be overheated blobs of material that could develop into convection rolls owing to instabilities of roll motion.     

Statistical theory of thermal instability

A new statistical approach is presented to study the thermal instability of an optically thin unmagnetized plasma. In the framework of this approach the time evolution of the mass distribution function over temperature φ( T ) is calculated. Function φ( T ) characterizes the statistical properties of the multiphase medium of arbitrarily spaced three-dimensional structure of arbitrary (small or large) temperature perturbations. We construct our theory under the isobarical condition ( P  = constant over space), which is satisfied in the short-wavelength limit of the perturbations. The developed theory is illustrated for the case of the thermal instability of a slowly expanding interstellar cloud (smooth scenario). Numerical solutions of equations of the statistical theory are constructed and compared with hydrodynamical solutions. The results of both approaches are identical in the short-wavelength range when the isobarity condition is satisfied. Also the limits of applicability of the statistical theory are estimated. The possible evolution of the initial spectrum of perturbations is discussed. The proposed theory and numerical models can be relevant to the formation of the two-phase medium in the ∼ 1 pc region around quasars. Then small warm ( T  ≃ 10⁴  K ) clouds are formed as the result of thermal instability in an expanded gas fragment, which is a product of either star–star or star–accretion disc collision.     

Magnetic draping of merging cores and radio bubbles in clusters of galaxies

Sharp fronts observed by the Chandra satellite between dense cool cluster cores moving with near-sonic velocity through the hotter intergalactic gas, require strong suppression of thermal conductivity across the boundary. This may be due to magnetic fields tangential to the contact surface separating the two plasma components. We point out that a super-Alfvenic motion of a plasma cloud (a core of a merging galaxy) through a weakly magnetized intercluster medium leads to 'magnetic draping', formation of a thin, strongly magnetized boundary layer with a tangential magnetic field. For supersonic cloud motion,   M _s≥ 1  , magnetic field inside the layer reaches near-equipartition values with thermal pressure. Typical thickness of the layer is  ∼ L / M ²_A≪ L   , where L is the size of the obstacle (plasma cloud) moving with Alfvén Mach number   M _A≫ 1  . To a various degree, magnetic draping occurs for both subsonic and supersonic flows, random and ordered magnetic fields and it does not require plasma compressibility. The strongly magnetized layer will thermally isolate the two media and may contribute to the Kelvin–Helmholtz stability of the interface. Similar effects occur for radio bubbles, quasi-spherical expanding cavities blown up by active galactic nucleus jets; in this case, the thickness of the external magnetized layer is smaller,  ∼ L / M ³_A≪ L   .     

A model of supernova feedback in galaxy formation

A model of supernova feedback during disc galaxy formation is developed. The model incorporates infall of cooling gas from a halo, and outflow of hot gas from a multiphase interstellar medium (ISM). The star formation rate is determined by balancing the energy dissipated in collisions between cold gas clouds with that supplied by supernovae in a disc marginally unstable to axisymmetric instabilities. Hot gas is created by thermal evaporation of cold gas clouds in supernova remnants, and criteria are derived to estimate the characteristic temperature and density of the hot component and hence the net mass outflow rate. A number of refinements of the model are investigated, including a simple model of a galactic fountain, the response of the cold component to the pressure of the hot gas, pressure-induced star formation and chemical evolution. The main conclusion of this paper is that low rates of star formation can expel a large fraction of the gas from a dwarf galaxy. For example, a galaxy with circular speed 50 km s¹ can expel 6080 per cent of its gas over a time-scale of 1 Gyr, with a star formation rate that never exceeds 0.1 M yr¹. Effective feedback can therefore take place in a quiescent mode and does not require strong bursts of star formation. Even a large galaxy, such as the Milky Way, might have lost as much as 20 per cent of its mass in a supernova-driven wind. The models developed here suggest that dwarf galaxies at high redshifts will have low average star formation rates and may contain extended gaseous discs of largely unprocessed gas. Such extended gaseous discs might explain the numbers, metallicities and metallicity dispersions of damped Lyman systems.     

The entropy and energy of intergalactic gas in galaxy clusters

Studies of the X-ray surface brightness profiles of clusters, coupled with theoretical considerations, suggest that the breaking of self-similarity in the hot gas results from an 'entropy floor', established by some heating process, which affects the structure of the intracluster gas strongly in lower-mass systems. By fitting analytical models for the radial variation in gas density and temperature to X-ray spectral images from the ROSAT PSPC and ASCA GIS, we have derived gas entropy profiles for 20 galaxy clusters and groups. We show that, when these profiles are scaled such that they should lie on top of one another in the case of self-similarity, the lowest-mass systems have higher-scaled entropy profiles than more massive systems. This appears to be due to a baseline entropy of depending on the extent to which shocks have been suppressed in low-mass systems. The extra entropy may be present in all systems, but is detectable only in poor clusters, where it is significant compared with the entropy generated by gravitational collapse. This excess entropy appears to be distributed uniformly with radius outside the central cooling regions.
We determine the energy associated with this entropy floor, by studying the net reduction in binding energy of the gas in low-mass systems, and find that it corresponds to a pre-heating temperature of 0.3 keV. Since the relationship between entropy and energy injection depends upon gas density, we are able to combine the excesses of 70140 keV cm² and 0.3 keV to derive the typical electron density of the gas into which the energy was injected. The resulting value of implies that the heating must have happened prior to cluster collapse but after a redshift z 710. The energy requirement is well matched to the energy from supernova explosions responsible for the metals which now pollute the intracluster gas.