Star Formation History and Extinction in the central kpc of M82-like Starbursts

We report on the star formation histories and extinction in the central kpc region of a sample of starburst galaxies that have similar far infrared (FIR), 10 micron and K-band luminosities as those of the archetype starburst M82. Our study is based on new optical spectra and previously published K-band photometric data, both sampling the same area around the nucleus. Model starburst spectra were synthesized as a combination of stellar populations of distinct ages formed over the Hubble time, and were fitted to the observed optical spectra and K-band flux. The model is able to reproduce simultaneously the equivalent widths of emission and absorption lines, the continuum fluxes between 3500-7000 Ang, the K-band and the FIR flux. We require a minimum of 3 populations -- (1) a young population of age < 8 Myr, with its corresponding nebular emission, (2) an intermediate-age population (age < 500 Myr), and (3) an old population that forms part of the underlying disk or/and bulge population. The contribution of the old population to the K-band luminosity depends on the birthrate parameter and remains above 60% in the majority of the sample galaxies. Even in the blue band, the intermediate age and old populations contribute more than 40% of the total flux in all the cases. A relatively high contribution from the old stars to the K-band nuclear flux is also apparent from the strength of the 4000 Ang break and the CaII K line. The extinction of the old population is found to be around half of that of the young population. The contribution to the continuum from the relatively old stars has the effect of diluting the emission equivalent widths below the values expected for young bursts. The mean dilution factors are found to be 5 and 3 for the Halpha and Hbeta lines respectively.Comment: 20 pages, uses emulateapj.cls. Scheduled to appear in ApJ Jan 1, 200

The boundary Riemann solver coming from the real vanishing viscosity approximation

We study a family of initial boundary value problems associated to mixed hyperbolic-parabolic systems: v^{\epsilon} _t + A (v^{\epsilon}, \epsilon v^{\epsilon}_x ) v^{\epsilon}_x = \epsilon B (v^{\epsilon} ) v^{\epsilon}_{xx} The conservative case is, in particular, included in the previous formulation. We suppose that the solutions $v^{\epsilon}$ to these problems converge to a unique limit. Also, it is assumed smallness of the total variation and other technical hypotheses and it is provided a complete characterization of the limit. The most interesting points are the following two. First, the boundary characteristic case is considered, i.e. one eigenvalue of $A$ can be $0$. Second, we take into account the possibility that $B$ is not invertible. To deal with this case, we take as hypotheses conditions that were introduced by Kawashima and Shizuta relying on physically meaningful examples. We also introduce a new condition of block linear degeneracy. We prove that, if it is not satisfied, then pathological behaviours may occur.Comment: 84 pages, 6 figures. Text changes in Sections 1 and 3.2.3. Added Section 3.1.2. Minor changes in other section

CAVITATION EROSION WEAR OF METALLIC SPECIMENS USING THE NEW COMPACT ROTATING DISK DEVICE

A more compact version of the rotating disk device (where a disk with cavitation inducers and specimens fixed on it rotates in water to provide cavitating flow) is used here to study the failure by cavitation in metallic specimens. The damage in the specimens is measured by mass loss versus exposition time to cavitating flow diagrams. Images of the test specimens obtained by electronic microscope are also shown and discussed. After 25 hours working in cavitating conditions, aluminum, brass and cast iron are weared by the cavitation phenomenon, resulting in pitting formation and mass loss. After each 5 hours operating in cavitating conditions, the specimens are cleaned by ultrasound, dried and weighted in a digital balance to obtain the mass loss in the process by comparing to its initial weight. All specimens, aluminum, brass and cast iron, were eroded by cavitation although in the aluminum the mass loss was lesser than the expected. The brass specimen was eroded as expected, mas no mass loss could be attained for the cast iron due to specimen oxidation

Boundedness and Stability of Impulsively Perturbed Systems in a Banach Space

Consider a linear impulsive equation in a Banach space $$\dot{x}(t)+A(t)x(t) = f(t), ~t \geq 0,$$ $$x(\tau_i +0)= B_i x(\tau_i -0) + \alpha_i,$$ with $\lim_{i \rightarrow \infty} \tau_i = \infty$. Suppose each solution of the corresponding semi-homogeneous equation $$\dot{x}(t)+A(t)x(t) = 0,$$ (2) is bounded for any bounded sequence $\{ \alpha_i \}$. The conditions are determined ensuring (a) the solution of the corresponding homogeneous equation has an exponential estimate; (b) each solution of (1),(2) is bounded on the half-line for any bounded $f$ and bounded sequence $\{ \alpha_i \}$ ; (c) $\lim_{t \rightarrow \infty}x(t)=0$ for any $f, \alpha_i$ tending to zero; (d) exponential estimate of $f$ implies a similar estimate for $x$.Comment: 19 pages, LaTex-fil

Cepheid Mass-loss and the Pulsation -- Evolutionary Mass Discrepancy

I investigate the discrepancy between the evolution and pulsation masses for Cepheid variables. A number of recent works have proposed that non-canonical mass-loss can account for the mass discrepancy. This mass-loss would be such that a 5Mo star loses approximately 20% of its mass by arriving at the Cepheid instability strip; a 14Mo star, none. Such findings would pose a serious challenge to our understanding of mass-loss. I revisit these results in light of the Padova stellar evolutionary models and find evolutionary masses are ($17\pm5$)% greater than pulsation masses for Cepheids between 5<M/Mo<14. I find that mild internal mixing in the main-sequence progenitor of the Cepheid are able to account for this mass discrepancy.Comment: 15 pages, 3 figures, ApJ accepte

UGC 3995: A Close Pair of Spiral Galaxies

UGC 3995 is a close pair of spiral galaxies whose eastern component hosts a Seyfert 2 nucleus. We present a detailed analysis of this system using long slit spectroscopy and narrow (\ha + \nii) as well as broad band (B, R) imaging and an archive WFPC2 image. The component galaxies reveal surprisingly small signs of interaction considering their spatial proximity and almost identical recession velocities, as the bright filament is probably an optical illusion due to the superposition of the bar of the Seyfert galaxy and of the spiral arms of the companion. The broad band morphology, a B--R color map, and a continuum-subtracted \ha + \nii image demonstrate that the western component UGC 3995B is in front of the Seyfert-hosting component UGC 3995A, partly obscuring its western side. The small radial velocity difference leaves the relative motion of the two galaxies largely unconstrained. The observed lack of major tidal deformations, along with some morphological peculiarities, suggests that the galaxies are proximate in space but may have recently approached each other on the plane of the sky. The geometry of the system and the radial velocity curve at P. A. = 106 suggest that the encounter may be retrograde or, alternatively, prograde before perigalacticon. The partial overlap of the two galaxies allows us to estimate the optical thickness of the disk of component B. We derive an extinction = 0.18 visual magnitudes in the infra-arms parts of the foreground galaxy disk, and >= 1-1.5 visual magnitudes in correspondence of the spiral arms.Comment: Accepted for publication in the Astronomical Journal (June 1999 issue

A posteriori error analysis and adaptive non-intrusive numerical schemes for systems of random conservation laws

In this article we consider one-dimensional random systems of hyperbolic conservation laws. We first establish existence and uniqueness of random entropy admissible solutions for initial value problems of conservation laws which involve random initial data and random flux functions. Based on these results we present an a posteriori error analysis for a numerical approximation of the random entropy admissible solution. For the stochastic discretization, we consider a non-intrusive approach, the Stochastic Collocation method. The spatio-temporal discretization relies on the Runge--Kutta Discontinuous Galerkin method. We derive the a posteriori estimator using continuous reconstructions of the discrete solution. Combined with the relative entropy stability framework this yields computable error bounds for the entire space-stochastic discretization error. The estimator admits a splitting into a stochastic and a deterministic (space-time) part, allowing for a novel residual-based space-stochastic adaptive mesh refinement algorithm. We conclude with various numerical examples investigating the scaling properties of the residuals and illustrating the efficiency of the proposed adaptive algorithm

Theoretical Modeling of Starburst Galaxies

We have modeled a large sample of infrared starburst galaxies using both the PEGASE v2.0 and STARBURST99 codes to generate the spectral energy distribution of the young star clusters. PEGASE utilizes the Padova group tracks while STARBURST99 uses the Geneva group tracks, allowing comparison between the two. We used our MAPPINGS III code to compute photoionization models which include a self-consistent treatment of dust physics and chemical depletion. We use the standard optical diagnostic diagrams as indicators of the hardness of the EUV radiation field in these galaxies. These diagnostic diagrams are most sensitive to the spectral index of the ionizing radiation field in the 1-4 Rydberg region. We find that warm infrared starburst galaxies contain a relatively hard EUV field in this region. The PEGASE ionizing stellar continuum is harder in the 1-4 Rydberg range than that of STARBURST99. As the spectrum in this regime is dominated by emission from Wolf-Rayet (W-R) stars, this difference is most likely due to the differences in stellar atmosphere models used for the W-R stars. We believe that the stellar atmospheres in STARBURST99 are more applicable to the starburst galaxies in our sample, however they do not produce the hard EUV field in the 1-4 Rydberg region required by our observations. The inclusion of continuum metal blanketing in the models may be one solution. Supernova remnant (SNR) shock modeling shows that the contribution by mechanical energy from SNRs to the photoionization models is << 20%. The models presented here are used to derive a new theoretical classification scheme for starbursts and AGN galaxies based on the optical diagnostic diagrams.Comment: 36 pages, 16 figures, to be published in ApJ, July 20, 200

A Wide Symbiotic Channel to Type Ia Supernovae

As a promising channel to Type Ia supernovae (SNe Ia), we have proposed a symbiotic binary system consisting of a white dwarf (WD) and a low mass red-giant (RG), where strong winds from the accreting WD play a key role to increase the WD mass to the Chandrasekhar mass limit. Here we propose two new evolutionary processes which make the symbiotic channel to SNe Ia much wider. (1) We first show that the WD + RG close binary can form from a wide binary even with such a large initial separation as $a_i \lesssim 40000 R_\odot$. Such a binary consists of an AGB star and a low mass main-sequence (MS) star, where the AGB star is undergoing superwind before becoming a WD. If the superwind at the end of AGB evolution is as fast as or slower than the orbital velocity, the wind outflowing from the system takes away the orbital angular momentum effectively. As a result the wide binary shrinks greatly to become a close binary. Therefore, the WD + RG binary can form from much wider binaries than our earlier estimate. (2) When the RG fills its inner critical Roche lobe, the WD undergoes rapid mass accretion and blows a strong optically thick wind. Our earlier analysis has shown that the mass transfer is stabilized by this wind only when the mass ratio of RG/WD is smaller than 1.15. Our new finding is that the WD wind can strip mass from the RG envelope, which could be efficient enough to stabilize the mass transfer even if the RG/WD mass ratio exceeds 1.15. With the above two new effects (1) and (2), the symbiotic channel can account for the inferred rate of SNe Ia in our Galaxy.Comment: 29 pages including 14 firgures, to be published in ApJ, 521, No.