In this paper, we present a new calculation of composition-dependent radiative cooling and heating curves of low-density gas, intended primarily for use in numerical simulations of galaxy formation and evolution. These curves depend on only five parameters: temperature, density, redshift, [Fe/H], and [Mg/Fe]. They are easily tabulated and can be efficiently interpolated during a simulation. The ionization equilibrium of 14 key elements is determined for temperatures between 10K and 10^9K and densities up to 100 amu/cm^3 taking into account collisional and radiative ionization, by the cosmic UV background and an interstellar radiation field, and by charge-transfer reactions. These elements, ranging from H to Ni, are the ones most abundantly produced and/or released by SNIa, SNII, and intermediate-mass stars. Self-shielding of the gas at high densities by neutral Hydrogen is taken into account in an approximate way by exponentially suppressing the H-ionizing part of the cosmic UV background for HI densities above a threshold density of n_HI,crit=0.007 cm^-3. We discuss how the ionization equilibrium, and the cooling and heating curves depend on the physical properties of the gas. The main advantage of the work presented here is that, within the confines of a well-defined chemical evolution model and adopting the ionization equilibrium approximation, it provides accurate cooling and heating curves for a wide range of physical and chemical gas properties, including the effects of self-shielding. The latter is key to resolving the formation of cold, neutral, high-density clouds suitable for star formation in galaxy simulations.

Simulation codes for galaxy formation and evolution take on board as many physical processes as possible beyond the standard gravitational and hydrodynamical physics. Most of this extra physics takes place below the resolution level of the simulations and is added in a sub-grid fashion. However, these sub-grid processes affect the macroscopic hydrodynamical properties of the gas and thus couple to the on-grid physics that is explicitly integrated during the simulation. In this paper, we focus on the link between partial ionization and the hydrodynamical equations. We show that the energy stored in ions and free electrons constitutes a potential energy term which breaks the linear dependence of the internal energy on temperature. Correctly taking into account ionization hence requires modifying both the equation of state and the energy-temperature relation. We implemented these changes in the cosmological simulation code Gadget2. As an example of the effects of these changes, we study the propagation of Sedov-Taylor shock waves through an ionizing medium. This serves as a proxy for the absorption of supernova feedback energy by the interstellar medium. Depending on the density and temperature of the surrounding gas, we find that up to 50% of the feedback energy is spent ionizing the gas rather than heating it. Thus, it can be expected that properly taking into account ionization effects in galaxy evolution simulations will drastically reduce the effects of thermal feedback. To the best of our knowledge, this potential energy term is not used in current simulations of galaxy formation and evolution.

We present a novel method to study the dynamics of bulk fermion systems such as the neutron-star crust. By introducing periodic boundary conditions into Fermionic Molecular Dynamics, it becomes possible to examine the long-range many-body correlations induced by antisymmetrisation in bulk fermion systems. The presented technique treats the spins and the fermionic nature of the nucleons explicitly and permits investigating the dynamics of the system. Despite the increased complexity related to the periodic boundary conditions, the proposed formalism remains computationally feasible.

Recently, we suggested that low-energy beta-beam neutrinos can be very useful for the study of supernova neutrino interactions. In this paper, we examine the use of a such experiment for the analysis of a supernova neutrino signal. Since supernova neutrinos are oscillating, it is very likely that the terrestrial spectrum of supernova neutrinos of a given flavor will not be the same as the energy distribution with which these neutrinos were first emitted. We demonstrate the efficacy of the proposed method for untangling multiple neutrino spectra. This is an essential feature of any model aiming at gaining information about the supernova mechanism, probing proto-neutron star physics, and understanding supernova nucleosynthesis, such as the neutrino process and the r-process. We also consider the efficacy of different experimental approaches including measurements at multiple beam energies and detector configurations.

Only weakly interacting, neutrinos are the principal messengers reaching us from the center of a supernova. Terrestrial neutrino telescopes, such as SNO and SuperKamiokande, can provide precious information about the processes in the core of a collapsing and exploding star. But the information about the supernova that a neutrino detector can supply, is restricted by the fact that little experimental data on the neutrino-nucleus cross sections exists and by the uncertainties in theoretical calculations. In this letter, we propose a novel procedure that determines the response of a target nucleus in a supernova-neutrino detector directly, by using low-energy beta-beams. We show that fitting 'synthetic' spectra, constructed by taking linear combinations of beta-beam spectra, to the original supernova-neutrino spectra reproduces the folded differential cross sections very accurately. Comparing the response in a terrestrial detector to these synthetic responses provides a direct way to determine the main parameters of the supernova-neutrino energy-distribution.

We study neutrino-induced nucleon knockout from nuclei. Expressions for the induced polarization are derived within the framework of the independent-nucleon model and the non-relativistic plane-wave approximation. Large dissimilarities in the nucleon polarization asymmetries are observed between neutrino- and antineutrino-induced processes. These asymmetries represent a potential way to distinguish between neutrinos and antineutrinos in neutral-current neutrino-scattering on nuclei. We discuss astrophysical applications of these polarization asymmetries. Our findings are illustrated for neutrino scattering on $^{16}$O and $^{208}$Pb.