Context. In analyses of stellar spectra and colours, and for the analysis of integrated light from galaxies, a homogeneous grid of model atmospheres of late-type stars and corresponding flux spectra is needed. Aims. We construct an extensive grid of spherically-symmetric models (supplemented with plane-parallel ones for the highest surface gravities), built on up-to-date atomic and molecular data, and make it available for public use. This includes "CN-cycled" models with C/N = 4.07 (solar), 1.5 and 0.5, C/O ranging from 0.09 to (normally) 5.0 to also represent stars of spectral types R, S and N, and with 1.0 ≤ ξ t ≤ 5 km s −1 . We also list thermodynamic quantities (T , P g , P e , ρ, partial pressures of molecules, etc.) and provide them on the World Wide Web, as well as calculated fluxes in approximately 108 000 wavelength points. Underlying assumptions in addition to 1D stratification (spherical or plane-parallel) include hydrostatic equilibrium, mixing-length convection and local thermodynamic equilibrium. We discuss a number of general properties of the models, in particular in relation to the effects of changing abundances, of blanketing, and of sphericity. We illustrate positive and negative feedbacks between sphericity and molecular blanketing. We compare the models with those of other available grids and find excellent agreement with planeparallel models of Castelli & Kurucz (if convection is treated consistently) within the overlapping parameter range. Although there are considerable departures from the spherically-symmetric NextGen models, the agreement with more recent PHOENIX models is gratifying.Conclusions. The models of the grid show considerable regularities, but some interesting departures from general patterns occur for the coolest models due to the molecular opacities. We have tested a number of approximate "rules of thumb" concerning effects of blanketing and sphericity and often found them to be astonishingly accurate. Some interesting new phenomena have been discovered and explored, such as the intricate coupling between blanketing and sphericity, and the strong effects of carbon enhancement on metal-poor models. We give further details of line absorption data for molecules, as well as details of models and comparisons with observations in subsequent papers.

No abstract

The nonlinear evolution of magnetized Keplerian shear ows is simulated in a local, three-dimensional model, including the e ects of compressibility and strati cation. Supersonic ows are initially generated by the Balbus-Hawley magnetic shear instability. The resulting ows regenerate a turbulent magnetic eld which, in turn, reinforces the turbulence. Thus, the system acts like a dynamo that generates its own turbulence. However, unlike usual dynamos, the magnetic energy exceeds the kinetic energy of the turbulence by a factor of 3{10. By assuming the eld to be vertical on the outer (upper and lower) surfaces we do not constrain the horizontal magnetic ux. Indeed, a large scale toroidal magnetic eld is generated, mostly in the form of toroidal ux tubes with lengths comparable to the toroidal extent of the box. This large scale eld is mainly of Present address: Nordita, Blegdamsvej 17, DK-2100 Copenhagen , Denmark y The National Center for Atmospheric Research is sponsored by the National Science Foundation 1 even (i.e. quadrupolar) parity with respect to the midplane and changes direction on a timescale of about 30 orbits, in a possibly cyclic manner. The e ective Shakura-Sunyaev alpha viscosity parameter is between 0.001 and 0.005, and the contribution from the Maxwell stress is about 3-7 times larger than the contribution from the Reynolds stress.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.Comment: 15 pages, 13 figures, to be submitted to MNRAS. For high-resolution figures, please see http://www-theorie.physik.unizh.ch/~agertz

Numerical simulations provide information on solar convection not available by direct observation. We present results of simulations of near surface solar convection with realistic physics : an equation of state including ionization and three-dimensional, LTE radiative transfer using a four-bin opacity distribution function. Solar convection is driven by radiative cooling in the surface thermal boundary layer, producing the familiar granulation pattern. In the interior of granules, warm plasma ascends with B10% ionized hydrogen. As it approaches and passes through the optical surface, the plasma cools, recombines, and loses entropy. It then turns over and converges into the dark intergranular lanes and further into the vertices between granulation cells. These vertices feed turbulent downdrafts below the solar surface, which are the sites of buoyancy work that drives the convection. Only a tiny fraction of the Ñuid ascending at depth reaches the surface to cool, lose entropy, and form the cores of these downdrafts. Granules evolve by pushing out against and being pushed in by their neighboring granules, and by being split by overlying Ñuid that cools and is pulled down by gravity. Convective energy transport properties that are closely related to integral constraints such as conservation of energy and mass are exceedingly robust. Other properties, which are less tightly constrained and/or involve higher order moments or derivatives, are found to depend more sensitively on the numerical resolution. At the highest numerical resolution, excellent agreement between simulated convection properties and observations is found. In interpreting observations it is crucial to remember that surfaces of constant optical depth are corrugated. The surface of unit optical depth in the continuum is higher above granules and lower in the intergranular lanes, while the surface of optical depth unity in a spectral line is corrugated in ways that are inÑuenced by both thermal and Doppler e †ects.

Transient heating events that formed calcium-aluminum-rich inclusions (CAIs) and chondrules are fundamental processes in the evolution of the solar protoplanetary disk, but their chronology is not understood. Using U-corrected Pb-Pb dating, we determined absolute ages of individual CAIs and chondrules from primitive meteorites. CAIs define a brief formation interval corresponding to an age of 4567.30 ± 0.16 million years (My), whereas chondrule ages range from 4567.32 ± 0.42 to 4564.71 ± 0.30 My. These data refute the long-held view of an age gap between CAIs and chondrules and, instead, indicate that chondrule formation started contemporaneously with CAIs and lasted ~3 My. This time scale is similar to disk lifetimes inferred from astronomical observations, suggesting that the formation of CAIs and chondrules reflects a process intrinsically linked to the secular evolution of accretionary disks.

We propose that the stellar initial mass function (IMF) is universal in the sense that its functional form arises as a consequence of the statistics of random supersonic flows.A model is developed for the origin of the stellar IMF, that contains a dependence on the average physical parameters (temperature, density, velocity dispersion) of the large scale site of star formation. The model is based on recent numerical experiments of highly supersonic random flows that have a strong observational counterpart.It is shown that a Miller-Scalo like IMF is naturally produced by the model for the typical physical conditions in molecular clouds. A more "massive" IMF in star bursts is also predicted.

This work presents a new physical model of the star formation rate (SFR), verified with an unprecedented set of large numerical simulations of driven, supersonic, self-gravitating, magnetohydrodynamic (MHD) turbulence, where collapsing cores are captured with accreting sink particles. The model depends on the relative importance of gravitational, turbulent, magnetic, and thermal energies, expressed through the virial parameter, α vir , the rms sonic Mach number, M S,0 , and the ratio of mean gas pressure to mean magnetic pressure, β 0 . The SFR is predicted to decrease with increasing α vir (stronger turbulence relative to gravity), to increase with increasing M S,0 (for constant values of α vir ), and to depend weakly on β 0 for values typical of star forming regions (M S,0 ≈ 4-20 and β 0 ≈ 1-20). In the unrealistic limit of β 0 → ∞, that is in the complete absence of a magnetic field, the SFR increases approximately by a factor of three, which shows the importance of magnetic fields in the star formation process, even when they are relatively weak (super-Alfvénic turbulence). In this non-magnetized limit, our definition of the critical density for star formation has the same dependence on α vir , and almost the same dependence on M S,0 , as in the model of Krumholz and McKee, although our physical derivation does not rely on the concepts of local turbulent pressure and sonic scale. However, our model predicts a different dependence of the SFR on α vir and M S,0 than the model of Krumholz and McKee. The star-formation simulations used to test the model result in an approximately constant SFR, after an initial transient phase. Both the value of the SFR and its dependence on the virial parameter found in the simulations are shown to agree very well with the theoretical predictions. A physical model of the SFR is needed for a realistic implementation of the star formation feedback in simulations of galaxy formation, and to retrieve the correct morphological and chemical evolution of galaxies. The new star formation law derived in this paper is suitable for such applications.