Implications of WMAP Observations On the Population III Star Formation Processes

In an earlier paper (Cen 2003) we pointed out the strong likelihood for the universal reionization to occur twice, giving rise to a larger Thomson optical depth. Here we perform a more focused analysis of the Thomson optical depth in light of the WMAP observations. While the current uncertainties on the observed Thomson optical depth are still relatively large, with tau_e=0.17+-0.04 (68%) (Kogut et al 2003), important implications on Pop III star formation processes at high redshift can already be inferred. We are able to draw four conclusions: (1) in the absence of a top-heavy initial stellar mass function (IMF) for Pop III metal-free stars and without a dramatic upturn in the star formation efficiency and ionizing photon escape fraction at high redshift (z>6), we find tau_e =< 0.09; (2) with a top-heavy IMF for the Pop III metal-free stars and plausible star formation efficiency and ionizing photon escape fraction, it is expected that tau_e =< 0.12; (3) it is possible to reach tau_e = 0.15, if the metal enrichment efficiency of the intergalactic medium by Pop III stars is very low thus Pop III era is prolonged; (4) to reach tau_e >= 0.17 requires either of the following two conditions: the cosmological model power index n is positively tilted to n >= 1.03, Pop III star formation in minihalos with molecular hydrogen cooling has an efficiency c_*(H_2,III)>0.01 (with ionizing photon escape fraction greater than 30%). If the current observed value of Thomson optical depth withstands future data, we will have strong observational evidence that Pop III stars are massive and their formation efficiency may be much higher than current theoretical works suggest. Alternatively, there may be unknown, non-stellar ionizing sources at very high redshift.Comment: a numerical error corrected, conclusions strengthened, submitted to ApJ Letters, 13 page

Simulation of stellar instabilities with vastly different timescales using domain decomposition

Strange mode instabilities in the envelopes of massive stars lead to shock waves, which can oscillate on a much shorter timescale than that associated with the primary instability. The phenomenon is studied by direct numerical simulation using a, with respect to time, implicit Lagrangian scheme, which allows for the variation by several orders of magnitude of the dependent variables. The timestep for the simulation of the system is reduced appreciably by the shock oscillations and prevents its long term study. A procedure based on domain decomposition is proposed to surmount the difficulty of vastly different timescales in various regions of the stellar envelope and thus to enable the desired long term simulations. Criteria for domain decomposition are derived and the proper treatment of the resulting inner boundaries is discussed. Tests of the approach are presented and its viability is demonstrated by application to a model for the star P Cygni. In this investigation primarily the feasibility of domain decomposition for the problem considered is studied. We intend to use the results as the basis of an extension to two dimensional simulations.Comment: 15 pages, 10 figures, published in MNRA

Galactic Twins of the Ring Nebula Around SN1987A and a Possible LBV-like Phase for Sk-69 202

Some core-collapse supernovae show clear signs of interaction with dense circumstellar material that often appears to be non-spherical. Circumstellar nebulae around supernova progenitors provide clues to the origin of that asymmetry in immediate pre-supernova evolution. Here I discuss outstanding questions about the formation of the ring nebula around SN1987A and some implications of similar ring nebulae around Galactic B supergiants. Several clues hint that SN1987A's nebula may have been ejected in an LBV-like event, rather than through interacting winds in a transition from a red supergiant to a blue supergiant.Comment: 2 pages, to appear in procedings of "Massive stars: fundamental parameters and circumstellar interactions", conference in honor of Virpi Niemela's 70th birthda

The Structure of the Homunculus. III. Forming a Disk and Bipolar Lobes in a Rotating Surface Explosion

We present a semi-analytic model for shaping the nebula around eta Carinae that accounts for the simultaneous production of bipolar lobes and an equatorial disk through a rotating surface explosion. Material is launched normal to the surface of an oblate rotating star with an initial kick velocity that scales approximately with the local escape speed. Thereafter, ejecta follow ballistic orbital trajectories, feeling only a central force corresponding to a radiatively reduced gravity. Our model is conceptually similar to the wind-compressed disk model of Bjorkman & Cassinelli, but we modify it to an explosion instead of a steady line-driven wind, we include a rotationally-distorted star, and we treat the dynamics somewhat differently. Continuum-driving avoids the disk inhibition that normally operates in line-driven winds. Our model provides a simple method by which rotating hot stars can simultaneously produce intrinsically bipolar and equatorial mass ejections, without an aspherical environment or magnetic fields. Although motivated by eta Carinae, the model may have generic application to other LBVs, B[e] stars, or SN1987A's nebula. When near-Eddington radiative driving is less influential, our model generalizes to produce bipolar morphologies without disks, as seen in many PNe.Comment: ApJ accepted, 9 page

Long-range interactions in the effective low energy Hamiltonian of Sr2IrO4: a core level resonant inelastic x-ray scattering study

We have investigated the electronic structure of Sr2IrO4 using core level resonant inelastic x-ray scattering. The experimental spectra can be well reproduced using ab initio density functional theory based multiplet ligand field theory calculations, thereby validating these calculations. We found that the low-energy, effective Ir t2g orbitals are practically degenerate in energy. We uncovered that covalency in Sr2IrO4, and generally in iridates, is very large with substantial oxygen ligand hole character in the Ir t2g Wannier orbitals. This has far reaching consequences, as not only the onsite crystal-field energies are determined by the long range crystal-structure, but, more significantly, magnetic exchange interactions will have long range distance dependent anisotropies in the spin direction. These findings set constraints and show pathways for the design of d^5 materials that can host compass-like magnetic interactions

A Multi-scale Approach for Simulations of Kelvin Probe Force Microscopy with Atomic Resolution

The distance dependence and atomic-scale contrast observed in nominal contact potential difference (CPD) signals recorded by KPFM on surfaces of insulating and semiconducting samples, have stimulated theoretical attempts to explain such effects. We attack this problem in two steps. First, the electrostatics of the macroscopic tip-cantilever-sample system is treated by a finite-difference method on an adjustable nonuniform mesh. Then the resulting electric field under the tip apex is inserted into a series of atomistic wavelet-based density functional theory (DFT) calculations. Results are shown for a realistic neutral but reactive silicon nano-scale tip interacting with a NaCl(001) sample. Bias-dependent forces and resulting atomic displacements are computed to within an unprecedented accuracy. Theoretical expressions for amplitude modulation (AM) and frequency modulation (FM) KPFM signals and for the corresponding local contact potential differences (LCPD) are obtained by combining the macroscopic and atomistic contributions to the electrostatic force component generated at the voltage modulation frequency, and evaluated for several tip oscillation amplitudes A up to 10 nm. Being essentially constant over a few Volts, the slope of atomistic force versus bias is the basic quantity which determines variations of the atomic-scale LCPD contrast. Already above A = 0.1 nm, the LCPD contrasts in both modes exhibit almost the same spatial dependence as the slope. In the AM mode, this contrast is approximately proportional to $A^{-1/2}$, but remains much weaker than the contrast in the FM mode, which drops somewhat faster as A is increased. These trends are a consequence of the macroscopic contributions to the KPFM signal, which are stronger in the AM-mode and especially important if the sample is an insulator even at sub-nanometer separations where atomic-scale contrast appears.Comment: 19 pages, 13 figure

On the stability of very massive primordial stars

The stability of metal-free very massive stars ($Z$ = 0; M = 120 - 500 \msol) is analyzed and compared with metal-enriched stars. Such zero-metal stars are unstable to nuclear-powered radial pulsations on the main sequence, but the growth time scale for these instabilities is much longer than for their metal-rich counterparts. Since they stabilize quickly after evolving off the ZAMS, the pulsation may not have sufficient time to drive appreciable mass loss in Z = 0 stars. For reasonable assumptions regarding the efficiency of converting pulsational energy into mass loss, we find that, even for the larger masses considered, the star may die without losing a large fraction of its mass. We find a transition between the $\epsilon$- and $\kappa$-mechanisms for pulsational instability at Z\sim 2\E{-4} - 2\E{-3}. For the most metal-rich stars, the $\kappa$-mechanism yields much shorter $e$-folding times, indicating the presence of a strong instability. We thus stress the fundamental difference of the stability and late stages of evolution between very massive stars born in the early universe and those that might be born today.Comment: 7 pages, 5 figures. Minor changes, more results given in Table 1, accepted for publication in Ap

A Chandra X-ray Study of NGC 1068: II. The Luminous X-ray Source Population

We present an analysis of the compact X-ray source population in the Seyfert~2 galaxy NGC 1068, imaged with Chandra. We find a total of 84 compact sources, of which 66 are projected onto the galactic disk of NGC 1068. Spectra of the brightest sources have been modeled with both multi-color disk blackbody and power-law models. The power-law model provides the better description of the spectrum for most of these sources. Five sources have 0.4-8 keV intrinsic luminosities greater than 10^{39} erg/s, assuming that their emission is isotropic and that they are associated with NGC 1068. We refer to these sources as Intermediate Luminosity X-ray Objects (IXOs). If these five sources are X-ray binaries accreting with luminosities that are both sub-Eddington and isotropic, then the implied source masses are >7 solar masses, and so they are inferred to be black holes. The brightest source has a much harder spectrum (Gamma = 0.9\pm0.1) than that found in Galactic black hole candidates and other IXOs. It also shows large-amplitude variability on both short-term and long-term timescales. The ratio of the number of sources with luminosities greater than 2.1 x 10^{38} erg/s in the 0.4-8 keV band to the rate of massive star formation is the same, to within a factor of two, for NGC 1068, the Antennae, NGC 5194 (the main galaxy in M51), and the Circinus galaxy. This suggests that the rate of production of X-ray binaries per massive star is approximately the same for galaxies with currently active star formation, including ``starbursts''.Comment: 33 pages, 10 figures. To appear in The Astrophysical Journal, v591 n1, July 1, 2003 issu

On the rotational dynamics of magnetically threaded disks around neutron stars

We investigate the rotational dynamics of disk accretion around a strongly magnetized neutron star with an aligned dipole field. The magnetospheric field is assumed to thread the disk plasma both inside and outside the corotation radius. As a result of disk-star interaction, the magnetic torque on the disk affects the structure of accretion flow to yield the observed spin- up or spin- down rates for a source of given fastness, magnetic field strength, and mass accretion rate. Within the model we obtain a prescription for the dynamical viscosity of such magnetically modified solutions for a Keplerian disk. We then use this prescription to find a model solution for the rotation rate profile throughout the entire disk, including the non-Keplerian inner disk. We find that the non-Keplerian angular velocity transition region is not necessarily narrow for a source of given spin state. The boundary layer approximation, as in the standard magnetically threaded disk model, holds only in the case of dynamical viscosity decreasing all the way to the innermost edge of the disk. These results are applied to several observed disk-fed X-ray pulsars that have exhibited quasi-periodic oscillations (QPOs). The QPO frequencies provide a constraint on the fastness parameter and enable one to determine uniquely the width of the angular velocity transition zone for each source within model assumptions. We discuss the implications of these results on the value of the critical fastness parameter for a magnetized star in spin equilibrium. Applications of our model are also made with relevant parameters from recent numerical simulations of quasi-stationary disk - magnetized star interactions