A simple approach to the supernova progenitor–explosion connection

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parametrized approaches and find good agreement with semi-empirical measures for the ‘explodability’ of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

Progenitor-Dependent Explosion Dynamics in Self-Consistent, Axisymmetric Simulations of Neutrino-Driven Core-Collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M ⊙, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

The dynamics of neutrino-driven supernova explosions after shock revival in 2D and 3D

We study the growth of the explosion energy after shock revival in neutrino-driven explosions in two and three dimensions (2D/3D) using multi-group neutrino hydrodynamics simulations of an 11.2 M⊙ star. The 3D model shows a faster and steadier growth of the explosion energy and already shows signs of subsiding accretion after one second. By contrast, the growth of the explosion energy in 2D is unsteady, and accretion lasts for several seconds as confirmed by additional long-time simulations of stars of similar masses. Appreciable explosion energies can still be reached, albeit at the expense of rather high neutron star masses. In 2D, the binding energy at the gain radius is larger because the strong excitation of downward-propagating g modes removes energy from the freshly accreted material in the downflows. Consequently, the mass outflow rate is considerably lower in 2D than in 3D. This is only partially compensated by additional heating by outward-propagating acoustic waves in 2D. Moreover, the mass outflow rate in 2D is reduced because much of the neutrino energy deposition occurs in downflows or bubbles confined by secondary shocks without driving outflows. Episodic constriction of outflows and vertical mixing of colder shocked material and hot, neutrino-heated ejecta due to Rayleigh–Taylor instability further hamper the growth of the explosion energy in 2D. Further simulations will be necessary to determine whether these effects are generic over a wider range of supernova progenitors.

Outer crust of a cold non-accreting magnetar

The outer-crust structure and composition of a cold, non-accreting magnetar are studied. We model the outer crust to be made of fully equilibrated matter where ionized nuclei form a Coulomb crystal embedded in an electron gas. The main effects of the strong magnetic field are those of quantizing the electron motion in Landau levels and of modifying the nuclear single-particle levels producing, on average, an increased binding of nucleons in nuclei present in the Coulomb lattice. The effect of a homogeneous and constant magnetic field on nuclear masses has been predicted by using a covariant density functional in which induced currents and axial deformation due to the presence of a magnetic field that breaks time-reversal symmetry have been included self-consistently in the nucleon and meson equations of motion. Although not yet observed, for Ba 1016 G both effects contribute to produce different compositions - odd-mass nuclei are frequently predicted - and to increase the neutron-drip pressure as compared to a typical neutron star. Specifically, in such a regime, the magnetic-field effects on nuclei favor the appearance of heavier nuclei at low pressures. As B increases, such heavier nuclei are also preferred up to larger pressures. For the most extreme magnetic field considered, B=1018 G, and for the models studied, almost the whole outer crust is made of 4092Zr52.

Triton-He-3 relative and differential flows as probes of the nuclear symmetry energy at supra-saturation densities

Using a transport model coupled with a phase-space coalescence afterburner, we study the triton-He-3 (t-He-3) ratio with both relative and differential transverse flows in semicentral Sn-132 + Sn-124 reactions at a beam energy of 400 MeV/nucleon. The neutron-proton ratios with relative and differential flows are also discussed as a reference. We find that similar to the neutron-proton pairs, the t-He-3 pairs also carry interesting information regarding the density dependence of the nuclear symmetry energy. Moreover, the nuclear symmetry energy affects more strongly the t-He-3 relative and differential flows than the pi(-)/pi(+) ratio in the same reaction. The t-He-3 relative flow can be used as a particularly powerful probe of the high-density behavior of the nuclear symmetry energy.

Isospin dependence of S-1(0) proton and neutron superfluidity in asymmetric nuclear matter

We investigate the S-1(0) neutron and proton superfluidity in isospin-asymmetric nuclear matter. We have concentrated on the isospin dependence of the pairing gaps and the effect of a microscopic three-body force. It is found that as the isospin asymmetry goes higher, the neutron S-1(0) superfluid phase shrinks gradually to a smaller density domain, whereas the proton one extends rapidly to a much wider density domain. The three-body force turns out to weaken the neutron S-1(0) superfluidity slightly, but it suppresses strongly the proton S-1(0) superfluidity at high densities in nuclear matter with large isospin asymmetry.

Studies of the star-forming structures in the dense interstellar medium: a view by dust extinction

New stars in galaxies form in dense, molecular clouds of the interstellar medium. Measuring how the mass is distributed in these clouds is of crucial importance for the current theories of star formation. This is because several open issues in them, such as the strength of different mechanism regulating star formation and the origin of stellar masses, can be addressed using detailed information on the cloud structure. Unfortunately, quantifying the mass distribution in molecular clouds accurately over a wide spatial and dynamical range is a fundamental problem in the modern astrophysics. This thesis presents studies examining the structure of dense molecular clouds and the distribution of mass in them, with the emphasis on nearby clouds that are sites of low-mass star formation. In particular, this thesis concentrates on investigating the mass distributions using the near infrared dust extinction mapping technique. In this technique, the gas column densities towards molecular clouds are determined by examining radiation from the stars that shine through the clouds. In addition, the thesis examines the feasibility of using a similar technique to derive the masses of molecular clouds in nearby external galaxies. The papers presented in this thesis demonstrate how the near infrared dust extinction mapping technique can be used to extract detailed information on the mass distribution in nearby molecular clouds. Furthermore, such information is used to examine characteristics crucial for the star formation in the clouds. Regarding the use of extinction mapping technique in nearby galaxies, the papers of this thesis show that deriving the masses of molecular clouds using the technique suffers from strong biases. However, it is shown that some structural properties can still be examined with the technique.

A Triggering Mechanism for Enhanced Star-Formation in Colliding Galaxies

We propose a physical mechanism to explain the origin of the intense burst of massive-star formation seen in colliding/merging, gas-rich, field spiral galaxies. We explicitly take account of the different parameters for the two main mass components, H-2 and H I, of the interstellar medium within a galaxy and follow their consequent different evolution during a collision between two galaxies. We also note that, in a typical spiral galaxy-like our galaxy, the Giant Molecular Clouds (GMCs) are in a near-virial equilibrium and form the current sites of massive-star formation, but have a low star formation rate. We show that this star formation rate is increased following a collision between galaxies. During a typical collision between two field spiral galaxies, the H I clouds from the two galaxies undergo collisions at a relative velocity of approximately 300 km s-1. However, the GMCs, with their smaller volume filling factor, do not collide. The collisions among the H I clouds from the two galaxies lead to the formation of a hot, ionized, high-pressure remnant gas. The over-pressure due to this hot gas causes a radiative shock compression of the outer layers of a preexisting GMC in the overlapping wedge region. This makes these layers gravitationally unstable, thus triggering a burst of massive-star formation in the initially barely stable GMCs.The resulting value of the typical IR luminosity from the young, massive stars from a pair of colliding galaxies is estimated to be approximately 2 x 10(11) L., in agreement with the observed values. In our model, the massive-star formation occurs in situ in the overlapping regions of a pair of colliding galaxies. We can thus explain the origin of enhanced star formation over an extended, central area approximately several kiloparsecs in size, as seen in typical colliding galaxies, and also the origin of starbursts in extranuclear regions of disk overlap as seen in Arp 299 (NGC 3690/IC 694) and in Arp 244 (NGC 4038/39). Whether the IR emission from the central region or that from the surrounding extranuclear galactic disk dominates depends on the geometry and the epoch of the collision and on the initial radial gas distribution in the two galaxies. In general, the central starburst would be stronger than that in the disks, due to the higher preexisting gas densities in the central region. The burst of star formation is expected to last over a galactic gas disk crossing time approximately 4 x 10(7) yr. We can also explain the simultaneous existence of nearly normal CO galaxy luminosities and shocked H-2 gas, as seen in colliding field galaxies.This is a minimal model, in that the only necessary condition for it to work is that there should be a sufficient overlap between the spatial gas distributions of the colliding galaxy pair.

Temperature profiles of accretion discs around rapidly rotating strange stars in general relativity: A comparison with neutron stars

We compute the temperature profiles of accretion discs around rapidly rotating strange stars, using constant gravitational mass equilibrium sequences of these objects, considering the full effect of general relativity. Beyond a certain critical value of stellar angular momentum (J), we observe the radius ( $r_{\rm orb}$) of the innermost stable circular orbit (ISCO) to increase with J (a property seen neither in rotating black holes nor in rotating neutron stars). The reason for this is traced to the crucial dependence of ${\rm d}r_{\rm orb}/{\rm d}J$ on the rate of change of the radial gradient of the Keplerian angular velocity at $r_{\rm orb}$ with respect to J. The structure parameters and temperature profiles obtained are compared with those of neutron stars, as an attempt to provide signatures for distinguishing between the two. We show that when the full gamut of strange star equation of state models, with varying degrees of stiffness are considered, there exists a substantial overlap in properties of both neutron stars and strange stars. However, applying accretion disc model constraints to rule out stiff strange star equation of state models, we notice that neutron stars and strange stars exclusively occupy certain parameter spaces. This result implies the possibility of distinguishing these objects from each other by sensitive observations through future X-ray detectors.

The intergalactic and circumgalactic medium surrounding star-forming galaxies at redshifts 2 < z < 3

We present measurements of the spatial distribution, kinematics, and physical properties of gas in the circumgalactic medium (CGM) of 2.0<z<2.8 UV color-selected galaxies as well as within the 2<z<3 intergalactic medium (IGM). These measurements are derived from Voigt profile decomposition of the full Lyα and Lyβ forest in 15 high-resolution, high signal-to-noise ratio QSO spectra resulting in a catalog of ∼6000 HI absorbers.

Chapter 2 of this thesis focuses on HI surrounding high-z star-forming galaxies drawn from the Keck Baryonic Structure Survey (KBSS). The KBSS is a unique spectroscopic survey of the distant universe designed to explore the details of the connection between galaxies and intergalactic baryons within the same survey volumes. The KBSS combines high-quality background QSO spectroscopy with large densely-sampled galaxy redshift surveys to probe the CGM at scales of ∼50 kpc to a few Mpc. Based on these data, Chapter 2 presents the first quantitative measurements of the distribution, column density, kinematics, and absorber line widths of neutral hydrogen surrounding high-z star-forming galaxies.

Chapter 3 focuses on the thermal properties of the diffuse IGM. This analysis relies on measurements of the ∼6000 absorber line widths to constrain the thermal and turbulent velocities of absorbing "clouds." A positive correlation between the column density of HI and the minimum line width is recovered and implies a temperature-density relation within the low-density IGM for which higher-density regions are hotter, as is predicted by simple theoretical arguments.

Chapter 4 presents new measurements of the opacity of the IGM and CGM to hydrogen-ionizing photons. The chapter begins with a revised measurement of the HI column density distribution based on this new absorption line catalog that, due to the inclusion of high-order Lyman lines, provides the first statistically robust measurement of the frequency of absorbers with HI column densities 14 ≲ log(N_HI/cm^-2) ≲ 17.2. Also presented are the first measurements of the column density distribution of HI within the CGM (50 <d < 300 pkpc) of high-z galaxies. These distributions are used to calculate the total opacity of the IGM and IGM+CGM and to revise previous measurements of the mean free path of hydrogen-ionizing photons within the IGM. This chapter also considers the effect of the surrounding CGM on the transmission of ionizing photons out of the sites of active star-formation and into the IGM.

This thesis concludes with a brief discussion of work in progress focused on understanding the distribution of metals within the CGM of KBSS galaxies. Appendix B discusses my contributions to the MOSFIRE instrumentation project.

Neutron source from thin foil confined by two laser pulses

Neutron production from a thin deuterium-tritium (D-T) foil irradiated by two intense femtosecond laser pulses from opposite sides with zero phase difference is studied analytically and numerically. For the interaction of a laser pulse of amplitude a = 7, focal area 100 mu m(2) and areal density 4.4 x 10(18) cm(-2) with a D-T plasma foil, about 1.17 x 10(21) neutron s(-1) can be obtained, much more than from other methods. The profiles of the ion and electron densities are also calculated.

An experimental study of a stellar neutron source: ^(13)C(α, n)^(16)O

Described in this thesis are measurements made of the thick-target neutron yield from the reaction ¹³C(α, n)¹⁶O. The yield was determined for laboratory bombarding energies between 0.475 and 0.700 MeV, using a stilbene crystal neutron detector and pulse-shape discrimination to eliminate gamma rays. Stellar temperatures between 2.5 and 4.5 x 10^{8 o}K are involved in this energy region. From the neutron yield was extracted the astrophysical cross-section factor S(E), which was found to fit a linear function: S(E) = [(5.48 ± 1.77) + (12.05 ± 3.91)E] x 10⁵ MeV-barns, center-of-mass system. The stellar rate of the ¹³C(α, n)¹⁶O reaction if calculated, and discussed with reference to helium burning and neutron production in the core of a giant star.

Results are also presented of measurements carried out on the reaction ⁹Be(α, n)¹²C, taken with a thin Be target. The bombarding energy-range covered was from 0.340 to 0.680 MeV, with excitation curves for the ground- and first excited-state neutrons being reported. Some angular distributions were also measured. Resonances were found at bombarding energies of E_LAB = 0.520 MeV (E_CM = 0.360 MeV, Γ ~ 55 keV CM, ωγ = 3.79 eV CM) and E_LAB = 0.600 MeV (E_CM = 0.415 MeV, Γ ˂ 4 keV CM, ωγ = 0.88 eV CM). The astrophysical rate of the ⁹Be(α, n)¹²C reaction due to these resonances is calculated.

Systematic study of the pi(-)/pi(+) ratio in heavy-ion collisions with the same neutron/proton ratio but different masses

A systematic study of the pi(-)/pi(+) ratio in heavy-ion collisions with the same neutron/proton ratio but different masses can help single out effects of the nuclear mean field on pion production. Based on simulations using the IBUU04 transport model, it is found that the pi(-)/pi(+) ratio in head-on collisions of Ca-48 + Ca-48, Sn-124 + Sn-124, and Au-197 + Au-197 at beam energies from 0.25 to 0.6 GeV/nucleon increases with increasing the system size or decreasing the beam energies. A comprehensive analysis of the dynamical isospin fractionation and the pi(-)/pi(+) ratio as well as their time evolution and spatial distributions demonstrates clearly that the pi(-)/pi(+) ratio is an effective probe of the high-density behavior of the nuclear symmetry energy.

Proton and neutron S-1(0) superfluidity in asymmetric nuclear matter

The proton and neutron S-1(0), pairing gaps and their isospin dependence in isospin asymmetric nuclear matter have been studied by the isospin dependent Brueckner-Hartree-Fock approach and the BCS theory. We have focused on investigating and discussing the effect of three-body force. The calculated results indicate that as the isospin asymmetry increases, the density range of the S-1(0) neutron superfluidity is narrowed slightly and the maximum value of the neutron pairing gap increases 9 while the density domain for the proton superfluidity enlarges rapidly and the peak value of the proton gap decreases remarkably. The three-body force turns out to affect only weakly the neutron S-1(0) superfluidity and its isospin dependence, i. e., it leads to a small reduction of the neutron S-1(0) paring gap. However, the three-body force not only reduces largely the strength of the proton S-1(0) gaps at high densities in highly asymmetric nuclear matter but also suppresses strongly the density domain for the proton S-1(0) superfluidity phase.