We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

Space asteroseismology reveals that stellar structure and evolution models of intermediate- and high-mass stars are in need of improvement in terms of angular momentum and chemical element transport. We aim to probe the interior structure of a hot massive star in the core-hydrogen burning phase of its evolution. We analyse CoRoT space photometry, Gaia DR2 space astrometry, and high-resolution high signal-to-noise HERMES and HARPS time-series spectroscopy of the slowly rotating B2V star HD 170580. From the time-series spectroscopy we derive $vsini=4\pm 2\,$km\,s$^{-1}$, where the uncertainty results from the complex pulsational line-profile variability. We detect 42 frequencies with amplitude above 5$\sigma$. Among these we identify 5 rotationally split triplets and 1 quintuplet. Asteroseismic modelling based on CoRoT, Gaia DR2 and spectroscopic data leads to a star of 8M$_\odot$ near core-hydrogen exhaustion and an extended overshoot zone. The detected low-order pressure-mode frequencies cannot be fit within the uncertainties of the CoRoT data by models without atomic diffusion. Irrespective of this limitation, the low-order gravity modes reveal HD 170580 to be a slow rotator with an average rotation period between 73 and 98 d and a hint of small differential rotation. Future Gaia DR3 data taking into account the multiplicity of the star, along with long-term TESS photometry would allow to put better observational constraints on the asteroseismic models of this blue evolved massive star. Improved modelling with atomic diffusion, including radiative levitation, is needed to achieve compliance with the low helium surface abundance of the star. This poses immense computational challenges but is required to derive the interior rotation and mixing profiles of this star.

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction.

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact.

We report on the discovery of a planetary companion candidate with a minimum mass Msini = 4.6 M_J orbiting the K2 III giant star HD 175370 (KIC 007940959). This star was a target in our program to search for planets around a sample of 95 giant stars observed with Kepler. This detection was made possible using precise stellar radial velocity measurements of HD 175370 taken over five years and four months using the coude echelle spectrograph of the 2-m Alfred Jensch Telescope and the fibre-fed echelle spectrograph HERMES of the 1.2-m Mercator Telescope. Our radial velocity measurements reveal a periodic (349.5 days) variation with a semi-amplitude K = 133 m/s, superimposed on a long-term trend. A low-mass stellar companion with an orbital period of ~88 years in a highly eccentric orbit and a planet in a Keplerian orbit with an eccentricity e = 0.22 are the most plausible explanation of the radial velocity variations. However, we cannot exclude the existence of stellar envelope pulsations as a cause for the low-amplitude radial velocity variations and only future continued monitoring of this system may answer this uncertainty. From Kepler photometry we find that HD 175370 is most likely a low-mass red-giant branch or asymptotic-giant branch star.

Star clusters are known as superb tools for understanding stellar evolution. In a quest for understanding the physical origin of magnetism and chemical peculiarity in about 7% of the massive main-sequence stars, we analysed two of the ten brightest members of the ~10 Myr old Galactic open cluster NGC 2264, the early B-dwarfs HD47887 and HD47777. We find accurate rotation periods of 1.95 and 2.64 days, respectively, from MOST photometry. We obtained ESPaDOnS spectropolarimetric observations, through which we determined stellar parameters, detailed chemical surface abundances, projected rotational velocities, and the inclination angles of the rotation axis. Because we found only small (<5 km/s) radial velocity variations, most likely caused by spots, we can rule out that HD47887 and HD47777 are close binaries. Finally, using the least-squares deconvolution technique, we found that both stars possess a large-scale magnetic field with an average longitudinal field strength of about 400 G. From a simultaneous fit of the stellar parameters we determine the evolutionary masses of HD47887 and HD47777 to be 9.4+/-0.7 M0 and 7.6+/-0.5 M0. Interestingly, HD47777 shows a remarkable helium underabundance, typical of helium-weak chemically peculiar stars, while the abundances of HD47887 are normal, which might imply that diffusion is operating in the lower mass star but not in the slightly more massive one. Furthermore, we argue that the rather slow rotation, as well as the lack of nitrogen enrichment in both stars, can be consistent with both the fossil and the binary hypothesis for the origin of the magnetic field. However, the presence of two magnetic and apparently single stars near the top of the cluster mass-function may speak in favour of the latter.

Our goal is to perform in-depth ensemble asteroseismology of the young open cluster NGC2244 with the 2-wheel Kepler mission. While the nominal Kepler mission already implied a revolution in stellar physics for solar-type stars and red giants, it was not possible to perform asteroseismic studies of massive OB stars because such targets were carefully avoided in the FoV in order not to disturb the exoplanet hunting. Now is an excellent time to fill this hole in mission capacity and to focus on the metal factories of the Universe, for which stellar evolution theory is least adequate. Our white paper aims to remedy major shortcomings in the theory of stellar structure and evolution of the most massive stars by focusing on a large ensemble of stars in a carefully selected young open cluster. Cluster asteroseismology of very young stars such as those of NGC2244 has the major advantage that all cluster stars have similar age, distance and initial chemical composition, implying drastic restrictions for the stellar modeling compared to asteroseismology of single isolated stars with very different ages and metallicities. Our study requires long-term photometric measurements of stars with visual magnitude ranging from 6.5 to 15 in a large FoV with a precision better than 30 ppm for the brightest cluster members (magnitude below 9) up to 500 ppm for the fainter ones, which is well achievable with 2-Wheel Kepler, in combination with high-precision high-resolution spectroscopy and spectro-polarimetry of the brightest pulsating cluster members. These ground-based spectroscopic data will be assembled with the HERMES and CORALIE spectrographs (twin 1.2m Mercator and Euler telescopes, La Palma, Canary Islands and La Silla, Chile), as well as with the spectro-polarimetric NARVAL instrument (2m BLT at the Pic du Midi, French Pyrenees), to which we have guaranteed access.

We present the CoRoT light curve of the bright B2.5V star HD 48977 observed during a short run of the mission in 2008, as well as a high-resolution spectrum gathered with the HERMES spectrograph at the Mercator telescope. We use several time series analysis tools to explore the nature of the variations present in the light curve. We perform a detailed analysis of the spectrum of the star to determine its fundamental parameters and its element abundances. We find a large number of high-order g-modes, and one rotationally induced frequency. We find stable low-amplitude frequencies in the p-mode regime as well. We conclude that HD 48977 is a new Slowly Pulsating B star with fundamental parameters found to be Teff = 20000 $\pm$ 1000 K and log(g)=4.2 $/pm$ 0.1. The element abundances are similar to those found for other B stars in the solar neighbourhood. HD 48977 was observed during a short run of the CoRoT satellite implying that the frequency precision is insufficient to perform asteroseismic modelling of the star. Nevertheless, we show that a longer time series of this star would be promising for such modelling. Our present study contributes to a detailed mapping of the instability strips of B stars in view of the dominance of g-mode pulsations in the star, several of which occur in the gravito-inertial regime.

We present modelling of the $\beta\,$Cep star HD\,180642 based on its observational properties deduced from CoRoT and ground-based photometry as well as from time-resolved spectroscopy. We investigate whether present-day state-of-the-art models are able to explain the full seismic behaviour of this star, which has extended observational constraints for this type of pulsator. We constructed a dedicated database of stellar models and their oscillation modes tuned to fit the dominant radial mode frequency of HD\,180642, by means of varying the hydrogen content, metallicity, mass, age, and core overshooting parameter. We compared the seismic properties of these models with those observed. We find models that are able to explain the numerous observed oscillation properties of the star, for a narrow range in mass of 11.4--11.8\,M$_\odot$ and no or very mild overshooting (with up to 0.05 local pressure scale heights), except for an excitation problem of the $\ell=3$, p$_1$ mode. We deduce a rotation period of about 13\,d, which is fully compatible with recent magnetic field measurements. The seismic models do not support the earlier claim of solar-like oscillations in the star. We instead ascribe the power excess at high frequency to non-linear resonant mode coupling between the high-amplitude radial fundamental mode and several of the low-order pressure modes. We report a discrepancy between the seismic and spectroscopic gravity at the $2.5\sigma$ level.

Diffusion of atoms can be important during quiescent phases of stellar evolution. Particularly in the very thin inert envelopes of subdwarf B stars, diffusive movements will considerably change the envelope structure and the surface abundances on a short timescale. Also, the subdwarfs will inherit the effects of diffusion in their direct progenitors, namely giants near the tip of the red giant branch. This will influence the global evolution and the pulsational properties of subdwarf B stars. We investigate the impact of gravitational settling, thermal diffusion and concentration diffusion on the evolution and pulsations of subdwarf B stars. Our diffusive stellar models are compared with models evolved without diffusion. We constructed subdwarf B models with a mass of 0.465 Msun from a 1 and 3 Msun ZAMS progenitor. The low mass star ignited helium in an energetic flash, while the intermediate mass star started helium fusion gently. For each progenitor type we computed series with and without atomic diffusion. Atomic diffusion in red giants causes the helium core mass at the onset of helium ignition to be larger. We find an increase of 0.0015 Msun for the 1 Msun model and 0.0036 Msun for the 3 Msun model. The effects on the red giant surface abundances are small after the first dredge up. The evolutionary tracks of the diffusive subdwarf B models are shifted to lower surface gravities and effective temperatures due to outward diffusion of hydrogen. This affects both the frequencies of the excited modes and the overall frequency spectrum. Especially the structure and pulsations of the post-non-degenerate sdB star are drastically altered, proving that atomic diffusion cannot be ignored in these stars.

We search for new variable B-type pulsators in the CoRoT data assembled primarily for planet detection, as part of CoRoT's Additional Programme. We aim to explore the properties of newly discovered B-type pulsators from the uninterrupted CoRoT space-based photometry and to compare them with known members of the Beta Cep and slowly pulsating B star (SPB) classes. We developed automated data analysis tools that include algorithms for jump correction, light-curve detrending, frequency detection, frequency combination search, and for frequency and period spacing searches. Besides numerous new, classical, slowly pulsating B stars, we find evidence for a new class of low-amplitude B-type pulsators between the SPB and Delta Sct instability strips, with a very broad range of frequencies and low amplitudes, as well as several slowly pulsating B stars with residual excess power at frequencies typically a factor three above their expected g-mode frequencies. The frequency data we obtained for numerous new B-type pulsators represent an appropriate starting point for further theoretical analyses of these stars, once their effective temperature, gravity, rotation velocity, and abundances will be derived spectroscopically in the framework of an ongoing FLAMES survey at the VLT.

Seismology of stars provides insight into the physical mechanisms taking place in their interior, with modes of oscillation probing different layers. Low-amplitude acoustic oscillations excited by turbulent convection were detected four decades ago in the Sun and more recently in low-mass main-sequence stars. Using data gathered by the Convection Rotation and Planetary Transits mission, we report here on the detection of solar-like oscillations in a massive star, V1449 Aql, which is a known large-amplitude (b Cephei) pulsator.

We present the results of a spectroscopic multisite campaign for the beta Cephei star 12 (DD) Lacertae. Our study is based on more than thousand high-resolution high S/N spectra gathered with 8 different telescopes in a time span of 11 months. In addition we make use of numerous archival spectroscopic measurements. We confirm 10 independent frequencies recently discovered from photometry, as well as harmonics and combination frequencies. In particular, the SPB-like g-mode with frequency 0.3428 1/d reported before is detected in our spectroscopy. We identify the four main modes as (l1,m1) = (1, 1), (l2,m2) = (0, 0), (l3,m3) = (1, 0) and (l4,m4) = (2, 1) for f1 = 5.178964 1/d, f2 = 5.334224 1/d, f3 = 5.066316 1/d and f4 = 5.490133 1/d, respectively. Our seismic modelling shows that f2 is likely the radial first overtone and that the core overshooting parameter alpha_ov is lower than 0.4 local pressure scale heights.

The information on stellar parameters and on the stellar interior we can get by studying pulsating stars depends crucially on the available observational constraints: both seismic constraints precision and number of detected modes, identification, nature of the modes) and "classical" observations (photospheric abundances, effective temperature, luminosity, surface gravity). We consider the case of beta Cephei pulsators and, with the aim of estimating quantitatively how the available observational constraints determine the type and precision of our inferences, we set the stage for Hare&Hound exercises. In this contribution we present preliminary results for one simple case, where we assume as "observed" frequencies a subset of frequencies of a model and then evaluate a seismic merit function on a dense and extensive grid of models of B-type stars. We also compare the behaviour of chi^2 surfaces obtained with and without mode identification.

The Liege Oscillation code can be used as a stand-alone program or as a library of subroutines that the user calls from a Fortran main program of his own to compute radial and non-radial adiabatic oscillations of stellar models. We describe the variables and the equations used by the program and the methods used to solve them. A brief account is given of the use and the output of the program.

Cles is an evolution code recently developed to produce stellar models meeting the specific requirements of studies in asteroseismology. It offers the users a lot of choices in the input physics they want in their models and its versatility allows them to tailor the code to their needs and implement easily new features. We describe the features implemented in the current version of the code and the techniques used to solve the equations of stellar structure and evolution. A brief account is given of the use of the program and of a solar calibration realized with it.

We present a seismic study of the beta Cephei star theta Ophiuchi. Our analysis is based on the observation of one radial mode, one rotationally split l = 1 triplet and three components of a rotationally split l = 2 quintuplet for which the m-values were well identified by spectroscopy. We identify the radial mode as fundamental, the triplet as p_1 and the quintuplet as g_1. Our NLTE abundance analysis results in a metallicity and CNO abundances in full agreement with the most recent updated solar values. With X \in [0.71,0.7211] and Z \in [0.009,0.015], and using the Asplund et al. (2005) mixture but with a Ne abundance about 0.3 dex larger (Cunha et al. 2006), the matching of the three independent modes, enables us to deduce constrained ranges for the mass (M = 8.2 +/- 0.3 Msun) and central hydrogen abundance (X_c = 0.38 +/- 0.02) of theta Oph and to prove the occurrence of core overshooting (alpha_ov = 0.44 +/- 0.07). We also derive an equatorial rotation velocity of 29 +/- 7 km/s. Moreover, we show that the observed non-equidistance of the l=1 triplet can be reproduced by second order effects of rotation. Finally, we show that the observed rotational splitting of two modes cannot rule out a rigid rotation model.

Models of alpha Cen A & B have been computed using the masses determined by Pourbaix et al. (2002) and the data derived from the spectroscopic analysis of Neuforge and Magain (1997). The seismological data obtained by Bouchy and Carrier (2001, 2002) do help improve our knowledge of the evolutionary status of the system. All the constraints are satisfied with a model which gives an age of about 6 Gyr for the binary.

We present an improved version of the method of photometric mode identification of Heynderickx et al. (1994). Our new version is based on the inclusion of precise non-adiabatic eigenfunctions determined in the outer stellar atmosphere according to the formalism recently proposed by Dupret et al.(2002). Our improved photometric mode identification technique is therefore no longer dependent on ad hoc parameters for the non-adiabatic effects. It contains the complete physical conditions of the outer atmosphere of the star, provided that rotation does not play a key role. We apply our improved method to the two slowly pulsating B stars HD 74560 and HD 138764 and to the beta Cephei star EN (16) Lac. Besides identifying the degree l of the pulsating stars, our method is also a tool for improving the knowledge of stellar interiors and atmospheres, by imposing constraints on parameters such as the metallicity and the mixing-length parameter alpha (a procedure we label non-adiabatic asteroseismology).

Some recent observations of the abundances of s-process, r-process, and alpha elements in metal-poor stars have led to a new scenario for their formation. According to this scenario, these stars were born in a globular cluster and accreted the s-process enriched gas expelled by cluster stars of higher-mass, thereby modifying their surface abundances. Later on, these polluted stars evaporated from the globular cluster to constitute an important fraction of the current halo population. In addition, there are now many direct observations of abundance anomalies not only in globular cluster giant stars but also in subgiant and main-sequence stars. Accretion provides again a plausible explanation for (at least some of) these peculiarities. Here we investigate further the efficiency of the accretion scenario. We find that in concentrated clusters with large escape velocities, accretion is very efficient and can indeed lead to major modifications of the stellar surface abundances.

Based on a model of globular cluster self-enrichment published in a previous paper, we present an explanation for the metallicity gradient observed throughout the galactic Old Halo. Our self-enrichment model is based on the ability of globular cluster progenitor clouds to retain the ejecta of a first generation of Type II Supernovae. The key point is that this ability depends on the pressure exerted on the progenitor cloud by the surrounding protogalactic medium and therefore on the location of the cloud in the protoGalaxy. Since there is no significant (if any) metallicity gradient in the whole halo, we also present a review in favour of a galactic halo partly build via accretions and mergers of satellite systems. Some of them bear their own globular clusters and therefore ``contaminate'' the system of globular clusters formed ``in situ'', namely within the original potential well of the Galaxy. Therefore, the comparison between our self-enrichment model and the observational data should be limited to the genuine galactic globular clusters, the so-called Old Halo group.

We present a model of globular cluster self-enrichment. In the protogalaxy, cold and dense clouds embedded in the hot protogalactic medium are assumed to be the progenitors of galactic halo globular clusters. The massive stars of a first generation of metal-free stars, born in the central areas of the proto-globular cluster clouds, explode as Type II supernovae. The associated blast waves trigger the expansion of a supershell, sweeping all the material of the cloud, and the heavy elements released by these massive stars enrich the supershell. A second generation of stars is born in these compressed and enriched layers of gas. These stars can recollapse and form a globular cluster. This work aims at revising the most often encountered argument against self-enrichment, namely the presumed ability of a small number of supernovae to disrupt a proto-globular cluster cloud. We describe a model of the dynamics of the supershell and of its progressive chemical enrichment. We show that the minimal mass of the primordial cluster cloud required to avoid disruption by several tens of Type II supernovae is compatible with the masses usually assumed for proto-globular cluster clouds. Furthermore, the corresponding self-enrichment level is in agreement with halo globular cluster metallicities.

Accurate relative abundances have been obtained for a sample of 21 mildly metal-poor stars from the analysis of high resolution and high signal-to-noise spectra. With these accurate results, correlations between relative abundances have been searched for, with a special emphasis on the neutron capture elements. This analysis shows that the r elements are closely correlated to the alpha elements, which is in agreement with the generally accepted idea that the r-process takes place during the explosion of massive stars. The situation is more complex as far as the s elements are concerned. Their relation with the alpha elements is not linear. In a first group of stars, the relative abundance of the s elements increases only slightly with the alpha elements overabundance until the latter reaches a maximum value. For the second group, the s elements show a rather large range of enhancement and a constant (and maximum) value of the alpha elements overabundance. This peculiar behaviour leads us to distinguish between two sub-populations of metal-poor stars, namely PopIIa (first group) and PopIIb (second group). We suggest a scenario of formation of metal-poor stars based on two distinct phases of chemical enrichment, a first phase essentially consisting in supernova explosions of massive stars, and a second phase where the enrichment is provided by stellar winds from intermediate mass stars. More specifically, we assume that all thick disk and field halo stars were born in globular clusters, from which they escaped, either during an early disruption of the cluster (PopIIa) or, later, through an evaporation process (PopIIb).

Photoionization by the high-redshift ultraviolet radiation background heats low density gas before it falls into dark matter potential wells, and it eliminates the neutral hydrogen and singly ionized helium that dominate cooling of primordial gas at temperatures of $10^4-10^5$K. We investigate the influence of photoionization on galaxy formation using high-resolution simulations with a 1-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code. We find that the presence of a photoionizing background suppresses the formation of galaxies with circular velocities $v_{circ} \simlt 30\,$\kms\ and substantially reduces the mass of cooled baryons in systems with circular velocities up to $v_{circ} \sim 50\,$\kms. Above $v_{circ} \sim 75\,$\kms, photoionization has no significant effect. Photoionization exerts its influence primarily by heating gas before collapse; the elimination of line cooling processes is less important. We discuss the implications of these results for hierarchical theories of galaxy formation.

We investigate the early formation of bound objects with masses comparable to the cosmological Jeans mass (10^5 solar masses). We follow the growth of isolated spherically symmetric density peaks starting from the linear perturbative regime. The initial parameters correspond to density peaks of various widths and heights in a Cold Dark Matter cosmology. We use a one-dimensional spherical Lagrangian hydrodynamics code to follow the dynamical, thermal, and non-equilibrium chemical evolution of the gas. The system includes a collisionless dark matter component and a baryonic component composed of the nine species H, H^-, H^+, He, He^+, He^{++}, H_2, H_2^+, and e^-. All relevant chemical reactions between these species and their cooling mechanisms are included in the calculations. We find that radiative cooling by H_2 affects the collapse dynamics of the gas only after it has already virialized and become part of the bound object. Therefore, radiative cooling is unlikely to have triggered the initial collapse of perturbations at redshifts z>10. Nevertheless, objects with baryonic masses well below the linear-theory Jeans mass (<10^3 solar masses) collapse due to shell crossing by the dark matter. Such objects could be the progenitors of a primordial population of high-mass stars in the intergalactic medium.

We describe an accurate, one-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code, designed to study the effects of radiative cooling and photo-ionization on the formation of protogalaxies. The code can treat an arbitrary number of fluid shells (representing baryons) and collisionless shells (representing cold dark matter). As a test of the code, we reproduce analytic solutions for the pulsation behavior of a polytrope and for the self-similar collapse of a spherically symmetric, cosmological perturbation. In this paper, we concentrate on the effects of radiative cooling, examining the ability of collapsing perturbations to cool within the age of the universe. In contrast to some studies based on order-of- magnitude estimates, we find that cooling arguments alone cannot explain the sharp upper cutoff observed in the galaxy luminosity function.

We have developed an accurate, one-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code, designed to study the effects of radiative cooling and photo-ionization on the formation of protogalaxies. We examine the ability of collapsing perturbations to cool within the age of the universe. In contrast to some studies based on order-of-magnitude estimates, we find that cooling arguments alone cannot explain the sharp upper cutoff observed in the galaxy luminosity function. We also look at the effect of a photoionizing background on the formation of low-mass galaxies.

We study the diffusion of helium and other heavy elements in the solar interior by solving exactly the set of flow equations developed by Burgers for a multi-component fluid, including the residual heat-flow terms. No approximation is made concerning the relative concentrations and no restriction is placed on the number of elements considered. We give improved diffusion velocities for hydrogen, helium, oxygen and iron, in the analytic form derived previously by Bahcall and Loeb. These expressions for the diffusion velocities are simple to program in stellar evolution codes and are expected to be accurate to $\sim 15\%$. Our complete treatment of element diffusion can be directly incorporated in a standard stellar evolution code by means of an exportable subroutine, but, for convenience, we also give simple analytical fits to our numerical results.