Forming stars on a viscous time-scale: the key to exponential stellar profiles in disc galaxies?

We argue for implementing star formation on a viscous time-scale in hydrodynamical simulations of disc galaxy formation and evolution. Modelling two-dimensional isolated disc galaxies with the Bhatnagar–Gross–Krook (BGK) hydrocode, we verify the analytic claim of various authors that if the characteristic time-scale for star formation is equal to the viscous time-scale in discs, the resulting stellar profile is exponential on several scalelengths whatever the initial gas and dark matter profile. This casts new light on both numerical and semi-analytical disc formation simulations that either (a) commence star formation in an already exponential gaseous disc, (b) begin a disc simulation with conditions known to lead to an exponential, i.e. the collapse of a spherically symmetric nearly uniform sphere of gas in solid-body rotation under the assumption of specific angular momentum conservation, or (c) in simulations performed in a hierarchical context, tune their feedback processes to delay disc formation until the dark matter haloes are slowly evolving and without much substructure so that the gas has the chance to collapse under conditions known to give exponentials. In such models, star formation follows a Schmidt-like law, which for lack of a suitable time-scale, resorts to an efficiency parameter. With star formation prescribed on a viscous time-scale, however, we find gas and star fractions after ∼12 Gyr that are consistent with observations without having to invoke a 'fudge factor' for star formation. Our results strongly suggest that despite our gap in understanding the exact link between star formation and viscosity, the viscous time-scale is indeed the natural time-scale for star formation.     

Meridional circulation in the radiation zones of rotating stars: Origins,behaviors and consequences on stellar evolution

Stellar radiation zones are the seat of meridional currents. This circulation has a strong impact on the transport of angular momentum and the mixing of chemicals that modify the evolution of stars. First, we recall in details the dynamical processes that are taking place in differentially rotating stellar radiation zones and the assumptions which are adopted for their modelling in stellar evolution. Then, we present our new results of numerical simulations which allow us to follow in 2D the secular hydrodynamics of rotating stars, assuming that anisotropic turbulence enforces a shellular rotation law and taking into account the transport of angular momentum by internal gravity waves. The different behaviors of the meridional circulation in function of the type of stars which is studied are discussed with their physical origin and their consequences on the transport of angular momentum and of chemicals. Finally, we show how this work is leading to a dynamical vision of the evolution of rotating stars from their birth to their death. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Metallicities and ages of stellar populations at a high Galactic latitude field

We present an analysis of UBVRI data from the selected area SA 141. By applying recalibrated methods of measuring ultraviolet excess (UVX), we approximate abundances and absolute magnitudes for 368 stars over 1.3 deg² out to distances over 10 kpc. With the density distribution constrained from our previous photometric parallax investigations and with sufficient accounting for the metallicity bias in the UVX method, we are able to compare the vertical abundance distribution to those measured in previous studies. We find that the abundance distribution has an underlying uniform component consistent with previous spectroscopic results that posit a monometallic thick disc and halo with abundances of  [Fe/H]=−0.8  and −1.4, respectively. However, there are a number of outlying data points that may indicate contamination by more metal-rich halo streams. The absence of vertical abundance gradients in the Galactic stellar populations and the possible presence of interloping halo streams would be consistent with expectations from merger models of Galaxy formation. We find that our UVX method has limited sensitivity in exploring the metallicity distribution of the distant Galactic halo, owing to the poor constraint on the UBV properties of very metal-poor stars. The derivation of metallicities from broad-band UBV photometry remains fundamentally sound for the exploration of the halo but is in need of both improved calibration and superior data.     

Transport of angular momentum and chemical species by anisotropic mixing in stellar radiative interiors

Small levels of turbulence can be present in stellar radiative interiors due to, e.g., the instability of rotational shear. In this paper we estimate turbulent transport coefficients for stably stratified rotating stellar radiation zones. Stable stratification induces strong anisotropy with a very small ratio of radial‐to‐horizontal turbulence intensities. Angular momentum is transported mainly due to the correlation between azimuthal and radial turbulent motions induced by the Coriolis force. This non‐diffusive transport known as the Λ‐effect has outward direction in radius and is much more efficient compared to the effect of radial eddy viscosity. Chemical species are transported by small radial diffusion only. This result is confirmed using direct numerical simulations combined with the test‐scalar method. As a consequence of the non‐diffusive transport of angular momentum, the estimated characteristic time of rotational coupling (≲100 Myr) between radiative core and convective envelope in young solar‐type stars is much shorter compared to the time‐scale of Lithium depletion (∼1 Gyr) (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hierarchical star formation: stars and stellar clusters in the Gould Belt

We perform a study of the spatial and kinematical distribution of young open clusters in the solar neighbourhood, discerning between bound clusters and transient stellar condensations within our sample. Then, we discriminate between Gould Belt (GB) and local Galactic disc (LGD) members, using our previous estimate of the structural parameters of both systems obtained from a sample of O-B6 Hipparcos stars. Single membership probabilities of the clusters are also calculated in the separation process. Using this classified sample, we analyse the spatial structure and the kinematic behaviour of the cluster system in the GB. The two star formation regions that dominate and give the GB its characteristic-inclined shape show a striking difference in their content of star clusters: while Ori OB1 is richly populated by open clusters, not a single one can be found within the boundaries of Sco OB2. This is mirrored in the velocity space, translating again into an abundance of clusters in the region of the kinematic space populated by the members of Ori OB1, and a marginal number of them associated with Sco OB2. We interpret all these differences by characterizing the Orion region as a cluster complex typically surrounded by a stellar halo, and the Sco-Cen region as an OB association in the outskirts of the complex. In the light of these results, we study the nature of the GB with respect to the optical segment of the Orion Arm, and we propose that the different content of star clusters, the different heights over the Galactic plane and the different residual velocities of Ori OB1 and Sco OB2 can be explained in terms of their relative position to the density maximum of the Local Arm in the solar neighbourhood. Although morphologically intriguing, the GB appears to be the result of our local and biased view of a larger star cluster complex in the Local Arm, that could be explained by the internal dynamics of the Galactic disc.     

The ages of pre-main-sequence stars

The position of pre-main-sequence or protostars in the Hertzsprung–Russell diagram is often used to determine their mass and age by comparison with pre-main-sequence evolution tracks. On the assumption that the stellar models are accurate, we demonstrate that, if the metallicity is known, the mass obtained is a good estimate. However, the age determination can be very misleading, because it is significantly (generally different by a factor of 2 to 5) dependent on the accretion rate and, for ages less than about 10⁶ yr, the initial state of the star. We present a number of accreting protostellar tracks that can be used to determine age if the initial conditions can be determined and the underlying accretion rate has been constant in the past. Because of the balance established between the Kelvin–Helmholtz, contraction time-scale and the accretion time-scale, a pre-main-sequence star remembers its accretion history. Knowledge of the current accretion rate, together with an HR-diagram position, gives information about the rate of accretion in the past, but does not necessarily improve any age estimate. We do not claim that ages obtained by comparison with these particular accreting tracks are likely to be any more reliable than those from comparisons with non-accreting tracks. Instead, we stress the unreliability of any such comparisons, and use the disparities between various tracks to estimate the likely errors in age and mass estimates. We also show how a set of coeval accreting objects do not appear coeval when compared with non-accreting tracks. Instead, accreting pre-main-sequence stars of around a solar mass are likely to appear older than those of either smaller or larger mass.     

Two populations of metal-free stars in the early Universe

We construct star formation histories at redshifts z ≳ 5 for two physically distinct populations of primordial, metal-free stars, motivated by theoretical and observational arguments that have hinted towards the existence of an intermediate stellar generation between Population III and Population I/II. Taking into account the cosmological parameters as recently revised by the Wilkinson Microwave Anisotropy Probe after three years of operation, we determine self-consistent reionization histories and discuss the resulting chemical enrichment from these early stellar generations. We find that the bulk of ionizing photons and heavy elements produced at high redshifts must have originated in Population II.5 stars, which formed out of primordial gas in haloes with virial temperatures ≳10⁴ K, and had typical masses ≳10 M_⊙. Classical Population III stars, formed in minihaloes and having masses ≳100 M_⊙, on the other hand, had only a minor impact on reionization and early metal enrichment. Specifically, we conclude that only ≃10 per cent by mass of metal-free star formation went into Population III.     

Binary formation in stellar clusters

We consider how the tidal potential of a stellar cluster or a dense molecular cloud affects the fragmentation of gravitationally unstable molecular cloud cores. We find that molecular cloud cores which would collapse to form a single star in the absence of tidal shear, can be forced to fragment if they are subjected to tides. This may enhance the frequency of binaries in star-forming regions such as Ophiuchus and the frequency of binaries with separations ≲100 au in the Orion Trapezium Cluster. We also find that clouds which collapse to form binary systems in the absence of a tidal potential will form bound binary systems if exposed to weak tidal shear. However, if the tidal shear is sufficiently strong, even though the cloud still collapses to form two fragments, the fragments are pulled apart while they are forming by the tidal shear and two single stars are formed. This sets an upper limit for the separation of binaries that form near dense molecular clouds or in stellar clusters.     

Rapid differential rotation in intermediate mass stars

ZAMS models with internal differential rotation are computed for 6.5 M_⊙ and 10 M_⊙ models with the 2D evolution code ROTORC. Two applications of these models, their pulsational frequencies and deduced locations in the HR diagram, are discussed in order to show how the internal differential rotation could be investigated. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

High-resolution simulations of stellar collisions between equal-mass main-sequence stars in globular clusters

We performed high-resolution simulations of two stellar collisions relevant for stars in globular clusters. We considered one head-on collision and one off-axis collision between two 0.6-M_⊙ main-sequence stars. We show that a resolution of about 100 000 particles is sufficient for most studies of the structure and evolution of blue stragglers. We demonstrate conclusively that collision products between main-sequence stars in globular clusters do not have surface convection zones larger than 0.004 M_⊙ after the collision, nor do they develop convection zones during the 'pre-main-sequence' thermal relaxation phase of their post-collision evolution. Therefore, any mechanism which requires a surface convection zone (i.e. chemical mixing or angular momentum loss via a magnetic wind) cannot operate in these stars. We show that no disc of material surrounding the collision product is produced in off-axis collisions. The lack of both a convection zone and a disc proves a continuing problem for the angular momentum evolution of blue stragglers in globular clusters.     

NLTE models of line-driven stellar winds – III. Influence of X-ray radiation on wind structure of O stars

We study the influence of X-rays on the wind structure of selected O stars. For this purpose we use our non-local thermodynamic equilibrium (NLTE) wind code with inclusion of additional artificial source of X-rays, assumed to originate in the wind shocks.
We show that the influence of shock X-ray emission on wind mass-loss rate is relatively small. Wind terminal velocity may be slightly influenced by the presence of strong X-ray sources, especially for stars cooler than   T _eff≲ 35 000 K  .
We discuss the origin of the   L _X/ L ∼ 10⁻⁷  relation. For stars with thick wind this relation can be explained assuming that the cooling time depends on wind density. Stars with optically thin winds exhibiting the 'weak wind problem' display enhanced X-ray emission which may be connected with large shock cooling length. We propose that this effect can explain the 'weak wind problem'.
Inclusion of X-rays leads to a better agreement of the model ionization structure with observations. However, we do not find any significant influence of X-rays on P  v ionization fraction implying that the presence of X-rays cannot explain the P  v problem.
We study the implications of modified ionization equilibrium due to shock emission on the line transfer in the X-ray region. We conclude that the X-ray line profiles of helium-like ions may be affected by the line absorption within the cool wind.     

High-velocity stars from decay of small stellar systems

In this paper we present numerical results on the decay of small stellar systems under different initial conditions (multiplicity 3 ≤  N  ≤ 10, and various mass spectra, initial velocities and initial configurations). The numerical treatment uses the CHAIN1 code (Mikkola &38; Aarseth). Particular attention is paid to the distribution of high-velocity escapers: we define these as stars with velocity above 30 km s⁻¹. These numerical experiments confirm that small N -body systems are dynamically unstable and produce cascades of escapers in the process of their decay. It is shown that the fraction of stars that escape from small dense stellar systems with an escape velocity greater than 30 km s⁻¹ is ∼1 per cent for all systems treated here. This relatively small fraction must be considered in relation to the rate of star formation in the Galaxy in small groups: this could explain some moderately high-velocity stars observed in the Galactic disc and possibly some young stars with relatively high metallicity in the thick disc.     

Physical properties of active stars and stellar systems

We study spectra of 18 active stars and stellar systems, drawn from lists of known active binaries or newly discovered extreme-ultraviolet and soft X-ray sources. Using spectral synthesis and template modelling, we measure (often for the first time) a variety of physical properties of the stars, including spectral types, luminosities and v  sin  i values. We report direct observation of previously undetected secondary components in three systems, and refine the orbital elements for two of these (HR 7275 and HD 217344).     

High-luminosity single carbon stars in stellar and galactic evolution

Summary In the solar neighborhood, approximately half of all intermediate mass main sequence stars with initially between 1 and about 5 Mbecome carbon stars with luminosities near 10⁴ L for typically less than 10⁶ years. These high luminosity carbon stars lose mass at rates nearly always in excess of 10^–7 M yr^–1 and sometimes in excess of 10^–5 M yr^–1. Locally, close to half of the mass returned into the interstellar medium by intermediate mass stars before they become white dwarfs is during the carbon star phase. A much greater fraction of lower metallicity stars become carbon-rich before they evolve into planetary nebulae than do higher metallicity stars; therefore, carbon stars are much more importan t in the outer than in the inner Galaxy.     

Non-emission-line young stars of intermediate mass

We present optical spectra of four intermediate-mass candidate young stellar objects that have often been classified as Herbig Ae/Be stars. Typical Herbig Ae/Be emission features are not present in the spectra of these stars. Three of them, HD 36917, HD 36982 and HD 37062, are members of the young Orion nebula cluster (ONC). This association constrains their ages to be ≲1 Myr. The lack of appreciable near-infrared excess in them suggests the absence of hot dust close to the central star. However, they do possess significant amounts of cold and extended dust as revealed by the large excess emission observed at far-infrared wavelengths. The fractional infrared luminosities  ( L _ir/ L _★)  and the dust masses computed from IRAS fluxes are systematically lower than those found for Herbig Ae/Be stars but higher than those for Vega-like stars. These stars may thus represent the youngest examples of the Vega phenomenon known so far. In contrast, the other star in our sample, HD 58647, is more likely to be a classical Be star, as is evident from the low   L _ir/ L _★  , the scarcity of circumstellar dust, the low polarization, the presence of H α emission and near-infrared excess, and the far-infrared spectral energy distribution consistent with free–free emission similar to other well-known classical Be stars.