The Milky Way, like other disc galaxies, underwent violent mergers and accretion of smaller satellite galaxies in its early history. Thanks to $\mathit Gaia$-DR2 and spectroscopic surveys, the stellar remnants of such mergers have been identified. The chronological dating of such events is crucial to uncover the formation and evolution of the Galaxy at high redshift, but it has so far been challenging owing to difficulties in obtaining precise ages for these oldest stars. Here we combine asteroseismology -- the study of stellar oscillations -- with kinematics and chemical abundances, to estimate precise stellar ages ($\sim$ 11%) for a sample of stars observed by the $\mathit{Kepler}$ space mission. Crucially, this sample includes not only some of the oldest stars that were formed inside the Galaxy, but also stars formed externally and subsequently accreted onto the Milky Way. Leveraging this unprecedented temporal resolution in age, we provide overwhelming evidence in favour of models in which the Galaxy had already formed a significant population of its stars (which now reside mainly in its thick disc) before the in-fall of the satellite galaxy $\mathit{Gaia}$-$\mathit{Enceladus}$ 10 billions years ago.

Context. The CoRoT and Kepler missions have paved the way for synergies between exoplanetology and asteroseismology. The use of seismic data helps providing stringent constraints on the stellar properties which directly impact the results of planetary studies. Amongst the most interesting planetary systems discovered by Kepler, Kepler-444 is unique by the quality of its seismic and classical stellar constraints. Its magnitude, age and the presence of 5 small-sized planets orbiting this target makes it an exceptional testbed for exoplanetology. Aims. We aim at providing a detailed characterization of Kepler-444, focusing on the dependency of the results on variations of key ingredients of the theoretical stellar models. This thorough study will serve as a basis for future investigations of the planetary evolution of the system orbiting Kepler-444. Methods. We use local and global minimization techniques to study the internal structure of the exoplanet-host star Kepler-444. We combine seismic observations from the Kepler mission, Gaia DR2 data and revised spectroscopic parameters to precisely constrain its internal structure and evolution. Results. We provide updated robust and precise determinations of the fundamental parameters of Kepler-444 and demonstrate that this low-mass star bore a convective core during a significant portion of its life on the main-sequence. Using seismic data, we are able to estimate the lifetime of the convective core to approximately 8 Gyr out of the 11 Gyr of the evolution of Kepler-444. The revised stellar parameters found by our thorough study are M = 0.754 +- 0.03M_Sun , R = 0.753 +- 0.01R_Sun , Age = 11 +- 1 Gy.

Models of core accretion assume that in the radiative zones of accreting gas envelopes, radiation diffuses. But super-Earths/sub-Neptunes (1-4$R_\oplus$, 2-20$M_\oplus$) point to formation conditions that are optically thin: their modest gas masses are accreted from short-lived and gas-poor nebulae reminiscent of the transparent cavities of transitional disks. Planetary atmospheres born in such environments can be optically thin to both incident starlight and internally generated thermal radiation. We construct time-dependent models of such atmospheres, showing that super-Earths/sub-Neptunes can accrete their $\sim$1%-by-mass gas envelopes, and super-puffs/sub-Saturns their $\sim$20%-by-mass envelopes, over a wide range of nebular depletion histories requiring no fine tuning. Although nascent atmospheres can exhibit stratospheric temperature inversions effected by atomic Fe and various oxides that absorb strongly at visible wavelengths, the rate of gas accretion remains controlled by the radiative-convective boundary (rcb) at much greater pressures. For dusty envelopes, the temperature at the rcb $T_{\rm rcb} \simeq 2500$ K is still set by ${\rm H}_2$ dissociation; for dust-depleted envelopes, $T_{\rm rcb}$ tracks the temperature of the visible or thermal photosphere, whichever is deeper, out to at least $\sim$5 AU. The rate of envelope growth remains largely unchanged between the old radiative diffusion models and the new optically thin models, reinforcing how robustly super-Earths form as part of the endgame chapter in disk evolution.

We present a new set of horizontal-branch (HB) models computed with the MESA stellar evolution code. The models adopt $\alpha$-enhanced \cite{ags09} metals mixtures and include the gravitational settling of He. They are used in our HB population synthesis tool to generate theoretical distributions of HB stars in order to describe the multiple stellar populations in the globular clusters 47Tuc, M3, and M13. The observed HB in 47Tuc is reproduced very well by our simulations for [Fe/H] $= -0.70$ and [$\alpha$/Fe] $= +0.4$ if the initial helium mass fraction varies by $\Delta Y_0 \sim 0.03$ and approximately 21%, 37%, and 42% of the stars have $Y_0 = 0.257$, 0.270, and 0.287, respectively. These simulations yield $(m-M)_V = 13.27$, implying an age near 13.0 Gyr. In the case of M3 and M13, our synthetic HBs for [Fe/H] $= -1.55$ and [$\alpha$/Fe] $= 0.4$ match the observed ones quite well if M3 has $\Delta Y_0 \sim 0.01$ and $(m-M)_V = 15.02$, resulting in an age of 12.6 Gyr, whereas M13 has $\Delta Y_0 \sim 0.08$ and $(m-M)_V = 14.42$, implying an age of 12.9 Gyr. Mass loss during giant-branch evolution and $\Delta Y_0$ appear to be the primary second parameters for M3 and M13. New observations for 7 of the 9 known RR Lyrae in M13 are also reported. Surprisingly, periods predicted for the $c$-type variables tend to be too high (by up to $\sim 0.1$~d).

A direct consequence of hierarchical galaxy formation is the existence of dual supermassive black holes (SMBHs), which may be preferentially triggered as active galactic nuclei (AGN) during galaxy mergers. Despite decades of searching, however, dual AGNs are extremely rare, and most have been discovered serendipitously. Using the all-sky WISE survey, we identified a population of over 100 morphologically identified interacting galaxies or mergers that display red mid-infrared colors often associated in extragalactic sources with powerful AGNs. The vast majority of these advanced mergers are optically classified as star-forming galaxies suggesting that they may represent an obscured population of AGNs that cannot be found through optical studies. In this work, we present Chandra/ACIS observations and near-infrared spectra with the Large Binocular Telescope of six advanced mergers with projected pair separations less than ~ 10 kpc. The combined X-ray, near-infrared, and mid-infrared properties of these mergers provide confirmation that four out of the six mergers host at least one AGN, with four of the mergers possibly hosting dual AGNs with projected separations less than ~10 kpc, despite showing no firm evidence for AGNs based on optical spectroscopic studies. Our results demonstrate that 1) optical studies miss a significant fraction of single and dual AGNs in advanced mergers, and 2) mid-infrared pre-selection is extremely effective in identifying dual AGN candidates in late-stage mergers. Our multi-wavelength observations suggest that the buried AGNs in these mergers are highly absorbed, with intrinsic column densities in excess of N_H >10^24cm^-2, consistent with hydrodynamic simulations.

The equilibrium-diffusion approximation (EDA) is used to describe certain radiation-hydrodynamic (RH) environments. When this is done the RH equations reduce to a simplified set of equations. The EDA can be derived by asymptotically analyzing the full set of RH equations in the equilibrium-diffusion limit. We derive the EDA this way and show that it and the associated set of simplified equations are both first-order accurate with transport corrections occurring at second order. Having established the EDA's first-order accuracy we then analyze the grey nonequilibrium-diffusion approximation and the grey Eddington approximation and show that they both preserve this first-order accuracy. Further, these approximations preserve the EDA's first-order accuracy when made in either the comoving-frame (CMF) or the lab-frame (LF). While analyzing the Eddington approximation, we found that the CMF and LF radiation-source equations are equivalent when neglecting ${\cal O}(\beta^2)$ terms and compared in the LF. Of course, the radiation pressures are not equivalent. It is expected that simplified physical models and numerical discretizations of the RH equations that do not preserve this first-order accuracy will not retain the correct equilibrium-diffusion solutions. As a practical example, we show that nonequilibrium-diffusion radiative-shock solutions devolve to equilibrium-diffusion solutions when the asymptotic parameter is small.

The individual characteristics of C, N, O, Na, Mg, Al, Si, and Fe on isochrones have been investigated in this study. Stellar models have been constructed for various mixtures in which the content of each element is changed up to the extreme value reported in recent studies, and the changes in isochrone shape have been analyzed for the various mixtures. To express the abundance variation of different elements with a single parameter, we have focused on the relative changes in the total number of metal ions. A review of the shape changes revealed that Na, Mg, and Al work the same way in stellar models, as the well-known fact that C, N, and O have the same reactions in the stellar interior. In addition, it was found that in high-metallicity conditions the influence of Si and Fe on the red giant branch becomes smaller than that of Na, Mg, and Al closer to the tip. Furthermore, the influence of Fe on the main sequence is larger than that of Na, Mg, Al, and even Si.

We present stellar atmosphere and evolution models of main sequence stars in two stellar populations of the Galactic globular cluster NGC 6752. These populations represent the two extremes of light-element abundance variations in the cluster. NGC 6752 is a benchmark cluster in the study of multiple stellar populations because of the rich array of spectroscopic abundances and panchromatic Hubble Space Telescope photometry. The spectroscopic abundances are used to compute stellar atmosphere and evolution models. The synthetic spectra for the two populations show significant differences in the ultraviolet and, for the coolest temperatures, in the near-infrared. The stellar evolution models exhibit insignificant differences in the H-R diagram except on the lower main sequence. The appearance of multiple sequences in the colour-magnitude diagrams (CMDs) of NGC 6752 is almost exclusively due to spectral effects caused by the abundance variations. The models reproduce the observed splitting and/or broadening of sequences in a range of CMDs. The ultraviolet CMDs are sensitive to variations in carbon, nitrogen, and oxygen but the models are not reliable enough to directly estimate abundance variations from photometry. On the other hand, the widening of the lower main sequence in the near-infrared CMD, driven by oxygen-variation via the water molecule, is well-described by the models and can be used to estimate the range of oxygen present in a cluster from photometry. We confirm that it is possible to use multiband photometry to estimate helium variations among the different populations, with the caveat that the estimated amount of helium-enhancement is model-dependent.

Canonical grids of stellar evolutionary sequences have been computed for the helium mass-fraction abundances $Y = 0.25$, $0.29$, and $0.33$, and for iron abundances that vary from $-2.4$ to $+0.4$ (in 0.2 dex increments) when [$\alpha$/Fe] $= +0.4$, or for the ranges $-2.0 \le$ [Fe/H] $\le +0.6$, $-1.8 \le$ [Fe/H] $\le +0.6$ when [$\alpha$/Fe] $= 0.0$ and $-0.4$, respectively. The grids, which consist of tracks for masses from $0.12 {{\cal M}_\odot}$ to $1.1$-$1.5 {{\cal M}_\odot}$ (depending on the metallicity) are based on up-to-date physics, including the gravitational settling of helium (but not metals diffusion). Interpolation software is provided to generate isochrones for arbitrary ages between $\approx 5$ and $15$ Gyr and any values of $Y$, [$\alpha$/Fe], and [Fe/H] within the aforementioned ranges. Comparisons of isochrones with published color-magnitude diagrams (CMDs) for the open clusters M67 ([Fe/H] $\approx 0.0$) and NGC 6791 ([Fe/H] $\approx 0.3$) and for four of the metal-poor globular clusters (47 Tuc, M3, M5, and M92) indicate that the models for the observed metallicities do a reasonably good job of reproducing the locations and slopes of the cluster main sequences and giant branches. The same conclusion is reached from a consideration of plots of nearby subdwarfs that have accurate Hipparcos parallaxes and metallicities in the range $-2.0 \le$ [Fe/H] $\le -1.0$ on various CMDs and on the ($\log T_{\rm eff},\,M_V$)-diagram. A relatively hot temperature scale similar to that derived in recent calibrations of the infrared flux method is favored by both the isochrones and the adopted color transformations, which are based on the latest MARCS model atmospheres.

Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite; the protostar is born. This sequence of events entail many physical processes, of which radiative transfer is of paramount importance. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency. In this paper, we follow-up on a previous paper on the collapse and formation of Larson's first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson's second core. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and 10 solar mass cloud cores. We find that, for a same central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass and temperature of the second cores are independent of initial cloud mass, size and temperature. Our simulations support the idea of a standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar masses.

We have computed low-mass stellar models and synthetic spectra for an initial chemical composition that includes the full C-N, O-Na, and Mg-Al abundance anticorrelations observed in second generation stars belonging to a number of massive Galactic globular clusters. This investigation extends a previous study that has addressed the effect of only the C-N and O-Na anticorrelations, seen in all globulars observed to date. We find that the impact of Mg-Al abundance variations at fixed [Fe/H] and Helium abundance is negligible on stellar models and isochrones (from the main sequence to the tip of the red giant branch) and bolometric corrections, when compared to the effect of C-N and O-Na variations. We identify a spectral feature at 490-520 nm, for low-mass main sequence stars, caused by MgH molecular bands. This feature has a vanishingly small effect on bolometric corrections for Johnson and Stroemgren filters that cover that spectral range. However, specific narrow-band filters able to target this wavelength region can be powerful tools to investigate the Mg-poor unevolved stars and highlight possible splittings of the MS due to variations of Mg abundances.

[Abbreviated] We have investigated the color-magnitude diagram of Omega Centauri and find that the blue main sequence (bMS) can be reproduced only by models that have a of helium abundance in the range Y=0.35-$0.40. To explain the faint subgiant branch of the reddest stars ("MS-a/RG-a" sequence), isochrones for the observed metallicity ([Fe/H]\approx0.7) appear to require both a high age (~13Gyr) and enhanced CNO abundances ([CNO/Fe]\approx0.9$). Y~0.35 must also be assumed in order to counteract the effects of high CNO on turnoff colors, and thereby to obtain a good fit to the relatively blue turnoff of this stellar population. This suggest a short chemical evolution period of time (<1Gyr) for Omega Cen. Our intermediate-mass (super-)AGB models are able to reproduce the high helium abundances, along with [N/Fe]~2 and substantial O depletions if uncertainties in the treatment of convection are fully taken into account. These abundance features distinguish the bMS stars from the dominant [Fe/H] $\approx1.7$ population. The most massive super-AGB stellar models (M_zams>=6.8M_sun, M_He,core>=1.245M_sun) predict too large N-enhancements, which limits their role in contributing to the extreme populations. We show quantitatively that highly He- and N-enriched AGB ejecta have particularly efficient cooling properties. Based on these results and on the reconstruction of the orbit of Omega Cen with respect to the Milky Way we propose the galactic plane passage gas purging scenario for the chemical evolution of this cluster. Our model addresses the formation and properties of the bMS population (including their central location in the cluster). We follow our model descriptively through four passage events, which could explain not only some key properties of the bMS, but also of the MS-a/RGB-a and the s-enriched stars.

Recent work has shown that most globular clusters have at least two chemically distinct components, as well as cluster-to-cluster differences in the mean [O/Fe], [Mg/Fe], and [Si/Fe] ratios at similar [Fe/H] values. In order to investigate the implications of variations in the abundances of these and other metals for H-R diagrams and predicted ages, grids of evolutionary sequences have been computed for scaled solar and enhanced alpha-element mixtures, and for mixtures in which the assumed [m/Fe] value for each of the metals C, N, O, Ne, Na, Mg, Si, S, Ca, and Ti has been increased, in turn, by 0.4 dex at constant [Fe/H]. These tracks, together with isochrones for ages from 6 to 14 Gyr, have been computed for -3.0 < [Fe/H] < -0.6, with helium abundances Y = 0.25, 0.29, and 0.33 at each [Fe/H] value, using upgraded versions of the Victoria stellar structure program and the Regina interpolation code, respectively. Turnoff luminosity versus age relations from isochrones are found to depend almost entirely on the importance of the CNO-cycle, and thereby mainly on the abundance of oxygen. Since C, N, and O, as well as Ne and S, do not contribute significantly to the opacities at low temperatures and densities, variations in their abundances do not impact the Teff scale of red giants. The latter is a strong function of the abundances of only Mg and Si (and Fe, possibly to a lesser extent), because they are so abundant and because they are strong sources of opacity at low temperatures. For these reasons, Mg and Si also have important effects on the temperatures of main-sequence stars. Due to their low abundances, Na, Ca, and Ti are of little consequence for stellar models. The effects of varying the adopted solar metals mix and the helium abundance at a fixed [Fe/H] are also briefly discussed.

We construct updated solar models with different sets of solar abundances, including the most recent determinations by Asplund et al. (2009). The latter work predicts a larger ($\sim 10%$) solar metallicity compared to previous measurements by the same authors but significantly lower ($\sim 25%$) than the recommended value from a decade ago by Grevesse & Sauval (1998). We compare the results of our models with determinations of the solar structure inferred through helioseismology measurements. The model that uses the most recent solar abundance determinations predicts the base of the solar convective envelope to be located at $R_{\rm CZ}= 0.724{\rm R_\odot}$ and a surface helium mass fraction of $Y_{\rm surf}=0.231$. These results are in conflict with helioseismology data ($R_{\rm CZ}= 0.713\pm0.001{\rm R_\odot}$ and $Y_{\rm surf}=0.2485\pm0.0035$) at 5$-\sigma$ and 11$-\sigma$ levels respectively. Using the new solar abundances, we calculate the magnitude by which radiative opacities should be modified in order to restore agreement with helioseismology. We find that a maximum change of $\sim 15%$ at the base of the convective zone is required with a smooth decrease towards the core, where the change needed is $\sim 5%$. The required change at the base of the convective envelope is about half the value estimated previously. We also present the solar neutrino fluxes predicted by the new models. The most important changes brought about by the new solar abundances are the increase by $\sim 10%$ in the predicted $^{13}$N and $^{15}$O fluxes that arise mostly due to the increase in the C and N abundances in the newly determined solar composition.

We present a new grid of stellar model calculations for stars on the Asymptotic Giant Branch between 1.0 and 6.0 M_sun. Our grid consists of 5 chemical mixtures between Z=0.0005 and Z=0.04, with both solar-like and $\alpha$-element enhanced metal ratios. We treat consistently the carbon-enhancement of the stellar envelopes by using opacity tables with varying C/O-ratio and by employing theoretical mass loss rates for carbon stars. The low temperature opacities have been calculated specifically for this project. For oxygen stars we use an empirical mass loss formalism. The third dredge-up is naturally obtained by including convective overshooting. Our models reach effective temperatures in agreement with earlier synthetic models, which included approximative carbon-enriched molecular opacities and show good agreement with empirically determined carbon-star lifetimes. A fraction of the models could be followed into the post-AGB phase, for which we provide models in a mass range supplementing previous post-AGB calculations. Our grid constitutes the most extensive set of AGB-models, calculated with the latest physical input data and treating carbon-enhancement due to the third dredge-up most consistently.

Among the newly discovered features of multiple stellar populations in Globular Clusters, the cluster NGC 1851 harbours a double subgiant branch, that can be explained in terms of two stellar generations, only slightly differing in age, the younger one having an increased total C+N+O abundance. Thanks to this difference in the chemistry, a fit can be made to the subgiant branches, roughly consistent with the C+N+O abundance variations already discovered two decades ago, and confirmed by recent spectroscopic data. We compute theoretical isochrones for the main sequence turnoff, by adopting four chemical mixtures for the opacities and nuclear reaction rates. The standard mixture has Z=10$^{-3}$ and [$\alpha$/Fe]=0.4, the others have C+N+O respectively equal to 2, 3 and 5 times the standard mixture, according to the element abundance distribution described in the text. We compare tracks and isochrones, and show how the results depend on the total CNO abundance. We notice that different initial CNO abundances between two clusters, otherwise similar in metallicity and age, may lead to differences in the turnoff morphology that can be easily attributed to an age difference. We simulate the main sequence and subgiant branch data for NGC 1851 and show that an increase of C+N+O by a factor $\sim$3 best reproduces the shift between the subgiant branches. We compare the main sequence width in the color m$_{F336W}$-m$_{F814W}$ with models, and find that the maximum helium abundance compatible with the data is Y$\simeq$0.29. We consider the result in the framework of the formation of the second stellar generation in globular clusters, for the bulk of which we estimate a helium abundance of Y$\simlt 0.26$.

To trace the source of the unique 13, 19.5, and 28 $\mu$m emission features in the spectra of oxygen-rich circumstellar shells around AGB stars, we have compared dust extinction spectra obtained by aerosol measurements. We have measured the extinction spectra for 19 oxide powder samples of eight different types, such as Ti-compounds (TiO, TiO$_2$, Ti$_2$O$_3$, Ti$_3$O$_5$, Al$_2$TiO$_5$, CaTiO$_3$), $\alpha$-, $\gamma$-, $\chi$-$\delta$-$\kappa$-Al$_2$O$_3$, and MgAl$_2$O$_4$ in the infrared region (10 - 50 $\mu$m) paying special attention to the morphological (size, shape, and agglomeration) effects and the differences in crystal structure. Anatase (TiO$_2$) particles with rounded edges are the possible 13, 19.5 and 28 $\mu$m band carriers as the main contributor in the spectra of AGB stars, and spherically shaped nano-sized spinel and Al$_2$TiO$_5$ dust grains are possibly associated with the anatase, enhancing the prominence of the 13 $\mu$m feature and providing additional features at 28 $\mu$m. The extinction data sets obtained by the aerosol and CsI pellet measurements have been made available for public use at http://elbe.astro.uni-jena.de

We present a new grid of stellar models and isochrones for old stellar populations, covering a large range of [Fe/H] values, for an heavy element mixture characterized by CNONa abundance anticorrelations as observed in Galactic globular cluster stars. The effect of this metal abundance pattern on the evolutionary properties of low mass stars, from the main sequence to the horizontal branch phase is analyzed. We perform comparisons between these new models, and our reference alpha-enhanced calculations, and discuss briefly implications for CMDs showing multiple main sequence or subgiant branches. A brief qualitative discussion of the effect of CN abundances on color-T_{eff} transformations is also presented, highlighting the need to determine theoretical color transformations for the appropriate metal mixture, if one wants to interpret observations in the Stroemgren system, or broadband filters blueward of the Johnson V-band.

The first paper in this series explored the effects of altering the chemical mixture of the stellar population on an element by element basis on stellar evolutionary tracks and isochrones to the end of the red giant branch. This paper extends the discussion by incorporating the fully consistent synthetic stellar spectra with those isochrone models in predicting integrated colors, Lick indices, and synthetic spectra. Older populations display element ratio effects in their spectra at higher amplitude than younger populations. In addition, spectral effects in the photospheres of stars tend to dominate over effects from isochrone temperatures and lifetimes, but, further, the isochrone-based effects that are present tend to fall along the age-metallicity degeneracy vector, while the direct stellar spectral effects usually show considerable orthogonality.

Observations suggest that a relationship exists between the driving mechanism of roAp star pulsations and the heavy element distribution in these stars. We attempt to study the effects of local and global metallicity variations on the excitation mechanism of high order p-modes in A star models. We developed stellar evolutionary models to describe magnetic A stars with different global metallicity or local metal accumulation profiles. These models were computed with CLES ("Code Li\`egeois d'\'evolution stellaire"), and the stability of our models was assessed with the non-adiabatic oscillation code MAD. Our models reproduce the blue edge of the roAp star instability strip, but generate a red edge hotter than the observed one, regardless of metallicity. Surprisingly, we find that an increase in opacity inside the driving region can produce a lower amount of driving, which we refer to as the "inverse $\kappa$-mechanism".

The ever-expanding depth and quality of photometric and spectroscopic observations of stellar populations increase the need for theoretical models in regions of age-composition parameter space that are largely unexplored at present. Stellar evolution models that employ the most advanced physics and cover a wide range of compositions are needed to extract the most information from current observations of both resolved and unresolved stellar populations. The Dartmouth Stellar Evolution Database is a collection of stellar evolution tracks and isochrones that spans a range of [Fe/H] from -2.5 to +0.5, [alpha/Fe] from -0.2 to +0.8 (for [Fe/H] <=0) or +0.2 (for [Fe/H] >0), and initial He mass fractions from Y=0.245 to 0.40. Stellar evolution tracks were computed for masses between 0.1 and 4 Msun, allowing isochrones to be generated for ages as young as 250 Myr. For the range in masses where the core He flash occurs, separate He-burning tracks were computed starting from the zero age horizontal branch. The tracks and isochrones have been transformed to the observational plane in a variety of photometric systems including standard UBV(RI)c, Stromgren uvby, SDSS ugriz, 2MASS JHKs, and HST ACS-WFC and WFPC2. The Dartmouth Stellar Evolution Database is accessible through a website at http://stellar.dartmouth.edu/~models/ where all tracks, isochrones, and additional files can be downloaded.

We present the first models allowing one to explore in a consistent way the influence of changes in the alpha-element-to-iron abundance ratio on the high-resolution spectral properties of evolving stellar populations. The models cover the wavelength range 300-1340nm at a resolution of FWHM=1AA, for metallicities in the range 0.005<=Z<=0.048 and stellar population ages 3 to 14 Gyr. These models are based on a recent library of synthetic stellar spectra and a new library of stellar evolutionary tracks, both computed for three different [Fe/H] (-0.5,0.0 and 0.2) and two different [alpha/Fe] (0.0 and 0.4). We expect our fully synthetic models to be primarily useful for evaluating the differential effect of changes in the alpha/Fe ratio on spectral properties such as broad-band colours and narrow spectral features. In addition, we assess the accuracy of absolute model predictions in two ways: first, by comparing the predictions of models for scaled-solar metal abundances [alpha/Fe]=0.0) to those of existing models based on libraries of observed stellar spectra; and secondly, by comparing the predictions of models for alpha-enhanced metal abundances ([alpha/Fe]=0.4) to observed spectra of massive early-type galaxies in the SDSS-DR4. We find that our models predict accurate strengths for those spectral indices that are strongly sensitive to the abundances of Fe and alpha elements. The predictions are less reliable for the strengths of other spectral features, such as those dominated by the abundances of C and N, as expected from the fact that the models do not yet allow one to explore the influence of these elements in an independent way. We conclude that our models are a powerful tool for extracting new information about the chemical properties of galaxies for which high-quality spectra have been gathered by modern surveys.

This paper is a demonstration of how the WKB approximation can be used to help solve the linearised 3D MHD equations. Using Charpit's Method and a Runge-Kutta numerical scheme, we have demonstrated this technique for a potential 3D magnetic null point, ${\bf{B}}=(x,\epsilon y -(\epsilon +1)z)$. Under our cold plasma assumption, we have considered two types of wave propagation: fast magnetoacoustic and Alfv\'en waves. We find that the fast magnetoacoustic wave experiences refraction towards the magnetic null point, and that the effect of this refraction depends upon the Alfv\'en speed profile. The wave, and thus the wave energy, accumulates at the null point. We have found that current build up is exponential and the exponent is dependent upon $\epsilon$. Thus, for the fast wave there is preferential heating at the null point. For the Alfv\'en wave, we find that the wave propagates along the fieldlines. For an Alfv\'en wave generated along the fan-plane, the wave accumulates along the spine. For an Alfv\'en wave generated across the spine, the value of $\epsilon$ determines where the wave accumulation will occur: fan-plane ($\epsilon=1$), along the $x-$axis ($0<\epsilon <1$) or along the $y-$axis ($\epsilon>1$). We have shown analytically that currents build up exponentially, leading to preferential heating in these areas. The work described here highlights the importance of understanding the magnetic topology of the coronal magnetic field for the location of wave heating.

The mass at which a transition is made between stars that have radiative or convective cores throughout the core H-burning phase is a fairly sensitive function of Z (particularly the CNO abundances). As a consequence, the ~4 Gyr, open cluster M67 provides a constraint on Z_\odot (and the solar heavy-element mixture) because (i) high-resolution spectroscopy indicates that this system has virtually the same metal abundances as the Sun, and (ii) its turnoff stars have masses just above the lower limit for sustained core convection on the main sequence. In this study, evolutionary tracks and isochrones using the latest MARCS model atmospheres as boundary conditions have been computed for 0.6-1.4 solar masses on the assumption of a metals mix (implying Z_\odot = 0.0125) based on the solar abundances derived by M. Asplund and collaborators using 3-D model atmospheres. These calculations do not predict a turnoff gap where one is observed in M67. No such difficulty is found if the analysis uses isochrones for Z_\odot = 0.0165, assuming the Grevesse & Sauval (1998) mix of heavy elements. Our findings, like the inferences from helioseismology, indicate a problem with the Asplund et al. abundances. However, it is possible that low-Z models with diffusive processes taken into account will be less problematic.

Integrated light from distant galaxies is often compared to stellar population models via the equivalent widths of spectral features--spectral indices--whose strengths rely on the abundances of one or more elements. Such comparisons hinge not only on the overall metal abundance but also on relative abundances. Studies have examined the influence of individual elements on synthetic spectra but little has been done to address similar issues in the stellar evolution models that underlie most stellar population models. Stellar evolution models will primarily be influenced by changes in opacities. In order to explore this issue in detail, twelve sets of stellar evolution tracks and isochrones have been created at constant heavy element mass fraction Z that self-consistently account for varying heavy element mixtures. These sets include scaled-solar, alpha-enhanced, and individual cases where the elements C, N, O, Ne, Mg, Si, S, Ca, Ti, and Fe have been enhanced above their scaled-solar values. The variations that arise between scaled-solar and the other cases are examined with respect to the H-R diagram and main sequence lifetimes.

The ACS Survey of Galactic Globular Clusters, an HST Treasury Project, will deliver high quality, homogeneous photometry of 65 globular clusters. This paper introduces a new collection of stellar evolution tracks and isochrones suitable for analyzing the ACS Survey data. Stellar evolution models were computed at [Fe/H]= -2.5, -2.0, -1.5, -1.0, -0.5, and 0; [alpha/Fe]= -0.2, 0, 0.2, 0.4, 0.6, and 0.8; and three initial He abundances for masses from 0.1 to 1.8 Msun and ages from 2 to 15 Gyr. Each isochrone spans a wide range in luminosity from Mv~14 up to the tip of the red giant branch. These are complemented by a set of He-burning tracks that extend from the zero age horizontal branch to the onset of thermal pulsations on the asymptotic giant branch. In addition, a set of computer programs are provided that make it possible to interpolate the isochrones in [Fe/H], generate luminosity functions from the isochrones, and create synthetic horizontal branch models. The tracks and isochrones have been converted to the observational plane with two different color-Teff transformations, one synthetic and one semi-empirical, in ground-based B, V, and I, and F606W and F814W for both ACS-WFC and WFPC2 systems. All models and programs presented in this paper are available from http://stellar.dartmouth.edu/~models/

Tables of mean opacities are often used to compute the transfer of radiation in a variety of astrophysical simulations from stellar evolution models to proto-planetary disks. Often tables, such as Ferguson et al. (2005), are computed with a predetermined set of physical assumptions that may or may not be valid for a specific application. This paper explores the effects of several assumptions of grain physics on the Rosseland mean opacity in an oxygen rich environment. We find that changing the distribution of grain sizes, either the power-law exponent or the shape of the distribution, has a marginal effect on the total mean opacity. We also explore the difference in the mean opacity between solid homogenous grains and grains that are porous or conglomorations of several species. Changing the amount of grain opacity included in the mean by assuming a grain-to-gas ratio significantly affects the mean opacity, but in a predictable way.

We are developing a brand new stellar population models with flexible chemistry (isochrones plus stellar colors and spectra) in order to set a new standard of completeness and excellence. Here we present preliminary results to assess the effects of stellar evolution models and stellar model atmosphere to the well-known Lick indices at constant heavy element mass fraction Z that self-consistently account for varying heavy element mixtures. We have enhanced chemical elements one by one. Our ultimate goal is to demonstrate 10% absolute mean ages for a sample of local galaxies derived from an integrated light spectrum.

Since the pioneering observations of Spite & Spite in 1982, the constant lithium abundance of metal-poor ([Fe/H]<-1.3) halo stars near the turnoff has been attributed a cosmological origin. Closer analysis revealed that the observed abundance lies at $\Delta$ Li~0.4 dex below the predictions of Big Bang Nucleosynthesis. The measurements of deuterium abundances on the lines of sight toward quasars and the recent data from the Wilkinson Microwave Anisotropy Probe independently confirmed this gap. We suggest here that part of the discrepancy is explained by the first generation of stars that efficiently depleted lithium. Assuming that the models for lithium evolution in the halo turnoff stars and $\Delta$ Li estimates are correct, we infer that between 1/3 and 1/2 of the baryonic matter of the early halo (~10^9 Mo) was processed through Population III stars. This new paradigm proposes a very economical solution to the lingering difficulty of understanding the properties of the Spite Plateau and its lack of star-to-star scatter down to [Fe/H]=-2.5. It is moreover in agreement both with the absence of lithium detection in the most metal-poor star presently known (HE 1327-2326), and also with new trends of the Plateau suggesting its low metallicity edge may be reached around [Fe/H]=-2.5. We discuss the turbulent mixing associated with enhanced supernovae explosions in the early interstellar medium in this picture. Finally we show that other chemical properties of the extremely metal-poor stars are in agreement with a significant Population III processing in the halo, provided these models include mass-loss and rotationally-induced mixing.

We investigate the influence of both a new generation of low-temperature opacities and of various amounts of alpha-element enhancements on stellar evolution models. New stellar models with two different alpha-element mixtures and two sets of appropriate opacity tables are computed and compared. The influence of the different mixtures as well as that of the improved generation of opacity tables is investigated. It is found that around solar metallicity the new opacity tables have a drastic influence on stellar temperatures, which is mainly an effect of the new low-temperature tables, and not of variations in alpha-element enhancement factors. The latter, however, influence stellar lifetimes via systematic opacity effects at core temperatures. We trace the reason for the low-temperature table changes to errors in the old tables. We conclude that variations in alpha-element abundance ratios affect the main-sequence properties of super-solar metallicity stars significantly. Red giant branch effective temperatures depend only slightly on the specific mixture. Our older low-temperature opacity tables were shown to be erroneous and should no longer be used for stellar models with near- or super-solar metallicity. Corrected tables have already been produced.

Self-enrichment processes occurring in the early stages of a globular cluster lifetime are generally invoked to explain the observed CNONaMgAl abundance anticorrelations within individual Galactic globulars.We have tested, with fully consistent stellar evolution calculations, if theoretical isochrones for stars born with the observed abundance anticorrelations satisfy the observational evidence that objects with different degrees of these anomalies lie on essentially identical sequences in the Color-Magnitude-Diagram (CMD). To this purpose, we have computed for the first time low-mass stellar models and isochrones with an initial metal mixture that includes the extreme values of the observed abundance anticorrelations, and varying initial He mass fractions. Comparisons with 'normal' alpha-enhanced isochrones and suitable Monte Carlo simulations that include photometric errors show that a significant broadening of the CMD sequences occurs only if the helium enhancement is extremely large (in this study, when Y=0.35) in the stars showing anomalous abundances. Stellar luminosity functions up to the Red Giant Branch tip are also very weakly affected, apart from - depending on the He content of the polluting material - the Red Giant Branch bump region. We also study the distribution of stars along the Zero Age Horizontal Branch, and derive general constraints on the relative location of objects with and without abundance anomalies along the observed horizontal branches of globular clusters.

Previous computations of low temperature Rosseland and Planck mean opacities from Alexander & Ferguson (1994) are updated and expanded. The new computations include a more complete equation of state with more grain species and updated optical constants. Grains are now explicitly included in thermal equilibrium in the equation of state calculation, which allows for a much wider range of grain compositions to be accurately included than was previously the case. The inclusion of high temperature condensates such as Al$_2$O$_3$ and CaTiO$_3$ significantly affects the total opacity over a narrow range of temperatures before the appearance of the first silicate grains. The new opacity tables are tabulated for temperatures ranging from 30000 K to 500 K with gas densities from 10$^{-4}$ g cm$^{-3}$ to 10$^{-19}$ g cm$^{-3}$. Comparisons with previous Rosseland mean opacity calculations are discussed. At high temperatures, the agreement with OPAL and Opacity Project is quite good. Comparisons at lower temperatures are more divergent as a result of differences in molecular and grain physics included in different calculations. The computation of Planck mean opacities performed with the opacity sampling method are shown to require a very large number of opacity sampling wavelength points; previously published results obtained with fewer wavelength points are shown to be significantly in error. Methods for requesting or obtaining the new tables are provided.

Important revisions of the solar model ingredients appear after 35 years of intense work which have led to an excellent agreement between solar models and solar neutrino detections. We first show that the updated CNO composition suppresses the anomalous position of the Sun in the known galactic enrichment. The following law: He/H= 0.075 + 44.6 O/H in fraction number is now compatible with all the indicators. We then examine the existing discrepancies between the standard model and solar - seismic and neutrino - observations and suggest some directions of investigation to solve them. We update our predicted neutrino fluxes using the recent composition, new nuclear reaction rates and seismic models as the most representative of the central plasma properties. We get 5.31 +- 0.6 10^6/cm^{2}/s for the total ^8B neutrinos, 66.5 SNU and 2.76 SNU for the gallium and chlorine detectors, all in remarquable agreement with the detected values including neutrino oscillations for the last two. We conclude that the acoustic modes and detected neutrinos see the same Sun, but that clear discrepancies in solar modelling encourage further observational and theoretical efforts.

We present the extension of our NextGen model atmosphere grid to the regime of giant stars. The input physics of the models presented here is nearly identical to the NextGen dwarf atmosphere models, however spherical geometry is used self-consistently in the model calculations (including the radiative transfer). We re-visit the discussion of the effects of spherical geometry on the structure of the atmospheres and the emitted spectra and discuss the results of NLTE calculations for a few selected models.

The launch of the Infrared Space Observatory and new atomic data have opened a window to the study of high ionization gas in active galactic nuclei (AGN). We present the results of a large number of photoionization simulations of the ``coronal line'' region in AGN, employing new atomic data from the Opacity and Iron Projects. Our grid of line emission spans 8 orders of magnitude in gas density and 14 orders of magnitude in ionizing flux in an effort to identify the optimal conditions in which these lines form. We show that coronal lines form at distances from just outside the broad line region to ~400L_{43.5}^{1/2} pc, in gas with ionization parameter -2.0 < log U(H) < 0.75, corresponding to gas densities of 10$^2$ to 10$^{8.5}$ cm$^{-3}$, with electron temperatures ~12,000K -- 150,000K. A large range of distances from the central source implies significant line width variation among the coronal lines. We identify several line ratios which could be used to measure relative abundances, and we use these to show that the coronal line gas is likely to be dust free.

The narrow emission line spectra of active galactic nuclei are not accurately described by simple photoionization models of single clouds. Recent Hubble Space Telescope images of Seyfert 2 galaxies show that these objects are rich with ionization cones, knots, filaments, and strands of ionized gas. Here we extend to the narrow line region the ``locally optimally emitting cloud'' (LOC) model, in which the observed spectra are predominantly determined by powerful selection effects. We present a large grid of photoionization models covering a wide range of physical conditions and show the optimal conditions for producing many of the strongest emission lines. We show that the integrated narrow line spectrum can be predicted by an integration of an ensemble of clouds, and we present these results in the form of diagnostic line ratio diagrams making comparisons with observations. We also predict key diagnostic line ratios as a function of distance from the ionizing source, and compare these to observations. The predicted radial dependence of the [O III]/[O II] ratio may be matched to the observed one in NGC4151, if the narrow line clouds see a more intense continuum than we see. The LOC scenario when coupled with a simple Keplerian gravitational velocity field will quite naturally predict the observed line width versus critical density relationship. The influence of dust within the ionized portion of the clouds is discussed and we show that the more neutral gas is likely to be dusty, although a high ionization dust-free region is most likely present too. This argues for a variety of NLR cloud origins.

Many large scale numerical simulations of astrophysical plasmas must also reproduce the hydrogen ionization and the resulting emission spectrum, in some cases quite accurately. We describe a compact model hydrogen atom that can be readily incorporated into such simulations. It reproduces the recombination efficiency and line spectrum predicted by much larger calculations for a very broad range of densities and temperatures. Uncertainties in hydrogen collision data are the largest source of differences between our compact atom and predictions of more extensive calculations, and underscore the need for accurate atomic data.