High-Redshift Galaxies in Cold Dark Matter Models

We use hydrodynamic cosmological simulations to predict the star formation properties of high-redshift galaxies (z=2-6) in five variants of the inflationary cold dark matter scenario, paying particular attention to z=3, the redshift of the largest "Lyman-break galaxy" (LBG) samples. Because we link the star formation timescale to the local gas density, the rate at which a galaxy forms stars is governed mainly by the rate at which it accretes cooled gas from the surrounding medium. At z=3, star formation in most of the simulated galaxies is steady on 200 Myr timescales, and the instantaneous star formation rate (SFR) is correlated with total stellar mass. However, there is enough scatter in this correlation that a sample selected above a given SFR threshold may contain galaxies with a fairly wide range of masses. The redshift history and global density of star formation in the simulations depend mainly on the amplitude of mass fluctuations in the underlying cosmological model. The three models whose mass fluctuation amplitudes agree with recent analyses of the Lyman-alpha forest also reproduce the observed luminosity function of LBGs reasonably well, though the dynamic range of the comparison is small and the theoretical and observational uncertainties are large. The models with higher and lower amplitudes appear to predict too much and too little star formation, respectively, though they are not clearly ruled out. The intermediate amplitude models predict SFR ~ 30-40 Msun/yr for galaxies with a surface density 1 per arcmin^2 per unit redshift at z=3. They predict much higher surface densities at lower SFR, and significant numbers of galaxies with SFR > 10 Msun/yr at z >= 5.Comment: Submitted to ApJ. 31 pages including 10 ps figures. Full resolution version of Fig 2 available at http://www.astronomy.ohio-state.edu/~dhw/Sph/zgal.fig2.ps.g

Solving One Dimensional Scalar Conservation Laws by Particle Management

We present a meshfree numerical solver for scalar conservation laws in one space dimension. Points representing the solution are moved according to their characteristic velocities. Particle interaction is resolved by purely local particle management. Since no global remeshing is required, shocks stay sharp and propagate at the correct speed, while rarefaction waves are created where appropriate. The method is TVD, entropy decreasing, exactly conservative, and has no numerical dissipation. Difficulties involving transonic points do not occur, however inflection points of the flux function pose a slight challenge, which can be overcome by a special treatment. Away from shocks the method is second order accurate, while shocks are resolved with first order accuracy. A postprocessing step can recover the second order accuracy. The method is compared to CLAWPACK in test cases and is found to yield an increase in accuracy for comparable resolutions.Comment: 15 pages, 6 figures. Submitted to proceedings of the Fourth International Workshop Meshfree Methods for Partial Differential Equation

Tidal spin-up of stars in dense stellar cusps around massive black holes

We show that main-sequence stars in dense stellar cusps around massive black holes are likely to rotate at a significant fraction of the centrifugal breakup velocity due to spin-up by hyperbolic tidal encounters. We use realistic stellar structure models to calculate analytically the tidal spin-up in soft encounters, and extrapolate these results to close and penetrating collisions using smoothed particle hydrodynamics simulations. We find that the spin-up falls off only slowly with distance from the black hole because the increased tidal coupling in slower collisions at larger distances compensates for the decrease in the stellar density. We apply our results to the stars near the massive black hole in the Galactic Center. Over their lifetime, ~1 Msol main sequence stars in the inner 0.3 pc of the Galactic Center are spun-up on average to ~10%--30% of the centrifugal breakup limit. Such rotation is ~20--60 times higher than is usual for such stars and may affect their subsequent evolution and their observed properties.Comment: 25 pages, 7 figures. Submitted to Ap

COSMOS: A Hybrid N-Body/Hydrodynamics Code for Cosmological Problems

We describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and individual galaxies. The code, named COSMOS, possesses several new features which distinguish it from other PM-PPM codes. In particular, to solve the Poisson equation we have written a new multigrid solver which can determine the gravitational potential of isolated matter distributions and which properly takes into account the finite-volume discretization required by PPM. All components of the code are constructed to work with a nonuniform mesh, preserving second-order spatial differences. The PPM code uses vacuum boundary conditions for isolated problems, preventing inflows when appropriate. The PM code uses a second-order variable-timestep time integration scheme. Radiative cooling and cosmological expansion terms are included. COSMOS has been implemented for parallel computers using the Parallel Virtual Machine (PVM) library, and it features a modular design which simplifies the addition of new physics and the configuration of the code for different types of problems. We discuss the equations solved by COSMOS and describe the algorithms used, with emphasis on these features. We also discuss the results of tests we have performed to establish that COSMOS works and to determine its range of validity.Comment: 43 pages, 14 figures, submitted to ApJS and revised according to referee's comment

The Chemical Compositions of the Type II Cepheids -- The BL Her and W Vir Variables

Abundance analyses from high-resolution optical spectra are presented for 19 Type II Cepheids in the Galactic field. The sample includes both short-period (BL Her) and long-period (W Vir) stars. This is the first extensive abundance analysis of these variables. The C, N, and O abundances with similar spreads for the BL Her and W Vir show evidence for an atmosphere contaminated with $3\alpha$-process and CN-cycling products. A notable anomaly of the BL Her stars is an overabundance of Na by a factor of about five relative to their presumed initial abundances. This overabundance is not seen in the W Vir stars. The abundance anomalies running from mild to extreme in W Vir stars but not seen in the BL Her stars are attributed to dust-gas separation that provides an atmosphere deficient in elements of high condensation temperature, notably Al, Ca, Sc, Ti, and $s$-process elements. Such anomalies have previously been seen among RV Tau stars which represent a long-period extension of the variability enjoyed by the Type II Cepheids. Comments are offered on how the contrasting abundance anomalies of BL Her and W Vir stars may be explained in terms of the stars' evolution from the blue horizontal branch.Comment: 41 pages including 11 figures and 4 tables; Accepted for publication in Ap

Numerical Simulations of Globular Cluster Formation

We examine various physical processes associated with the formation of globular clusters by using the three-dimensional Smoothed Particle Hydrodynamics (SPH) code. Our code includes radiative cooling of gases, star formation, energy feedback from stars including stellar winds and supernovae, and chemical enrichment by stars. We assume that, in the collapsing galaxy, isothermal cold clouds form through thermal condensations and become proto-globular clouds. We calculate the size of proto-globular clouds by solving the linearized equations for perturbation. We compute the evolution of the inner region of the proto-cloud with our SPH code for various initial radius and initial composition of gases. When the initial gases contain no heavy elements, the evolution of proto-clouds sensitively depends on the initial radius. For a smaller initial radius, the initial star burst is so intense that the subsequent star formation occurs in the central regions to form a dense star cluster as massive as the globular cluster. When the initial gases contain some heavy elements, the metallicity of gases affects the evolution and the final stellar mass. If the initial radius of the proto-globular clouds was relatively large, the formation of a star cluster as massive as the globular clusters requires the initial metallicity as high as [Fe/H] $\geq -2$. The self-enrichment of heavy elements in the star cluster does not occur in all cases.Comment: Accpeted for publication in the ApJ. Correctiong errors in Table