Neutron stars provide an excellent laboratory for physics under the most extreme conditions. Up to now, models of axisymmetric, stationary, differentially rotating neutron stars were constructed under the strong assumption of barotropicity, where a one-to-one relation between all thermodynamic quantities exists. This implies that the specific angular momentum of a matter element depends only on its angular velocity. The physical conditions in the early stages of neutron stars, however, are determined by their violent birth processes, typically a supernova or in some cases the merger of two neutron stars, and detailed numerical models show that the resulting stars are by no means barotropic. Here, we construct models for stationary, differentially rotating, non-barotropic neutron stars, where the equation of state and the specific angular momentum depend on more than one independent variable. We show that the potential formulation of the relativistic Euler equation can be extended to the non-barotropic case, which, to the best of our knowledge, is a new result even for the Newtonian case. We implement the new method into the XNS code and construct equilibrium configurations for non-barotropic equations of state. We scrutinize the resulting configurations by evolving them dynamically with the numerical relativity code BAM, thereby demonstrating that the new method indeed produces stationary, differentially rotating, non-barotropic neutron star configurations.

Short gamma-ray bursts (sGRBs) show a large diversity in their properties. This suggests that the observed phenomenon can be caused by different "central engines" or that the engine produces a variety of outcomes depending on its parameters, or possibly both. The most popular engine scenario, the merger of two neutron stars, has received support from the recent Fermi and INTEGRAL detection of a burst of gamma rays (GRB170817A) following the neutron star merger GW170817, but at the moment it is not clear how peculiar this event potentially was. Several sGRBs engine models involve the collapse of a supramassive neutron star that produces a black hole plus an accretion disk. We study this scenario for a variety of equations of states both via angular momentum considerations based on equilibrium models and via fully dynamical Numerical Relativity simulations. We obtain a broader range of disk forming configurations than earlier studies but we agree with the latter that none of these configurations is likely to produce a phenomenon that would be classified as an sGRB.

A proto-neutron star (PNS) is the first phase of life of a neutron star, and is likely to origin from a core-collapse supernova. After about 200 ms from core-collapse, the PNS evolution may be modeled as a sequence of quasi-stationary configurations. These configurations depend on the PNS thermodynamic profiles, whose evolution largely depends upon the neutrino diffusion. We developed a new PNS evolutionary code that solves by iteration the neutrino number and energy transport equations together with the relativistic stellar structure equations assuming spherical symmetry. The neutrino cross sections are determined consistently with the underlying equation of state (EoS). To include the EoSs in the evolution, we devised and tested a new fitting formula for the interacting part of the baryon free-energy, valid at finite temperature and arbitrary degeneracy. Using our code, we provide estimates for the neutrino signal in the Super-Kamiokande III detector and the frequencies of the gravitational waves due to stellar oscillations, for three stellar masses and three nucleonic EoSs. For the first time we evolve a PNS with a nuclear many-body theory EoS in a consistent way, that is, we take into account realistic nuclear interactions in the computation of the neutrino cross sections. By including rotation in an effective way, we have also determined the time variation of the rotation frequency due to PNS contraction and neutrino angular momentum loss, and the gravitational wave signal due to rotation. We find that the mass shedding limit restricts the initial angular momentum. Consequently, the final rotation frequency has to be smaller than about 300 Hz for a PNS of about 1.6 solar masses whose EoS is described by the GM3 mean-field model.

In a core-collapse supernova, a huge amount of energy is released in the Kelvin-Helmholtz phase subsequent to the explosion, when the proto-neutron star cools and deleptonizes as it loses neutrinos. Most of this energy is emitted through neutrinos, but a fraction of it can be released through gravitational waves. We model the evolution of a proto-neutron star in the Kelvin-Helmholtz phase using a general relativistic numerical code, and a recently proposed finite temperature, many-body equation of state; from this we consistently compute the diffusion coefficients driving the evolution. To include the many-body equation of state, we develop a new fitting formula for the high density baryon free energy at finite temperature and intermediate proton fraction. We estimate the emitted neutrino signal, assessing its detectability by present terrestrial detectors, and we determine the frequencies and damping times of the quasi-normal modes which would characterize the gravitational wave signal emitted in this stage.

We study the evolution of the rotation rate of a proto-neutron star, born in a core-collapse supernova, in the first seconds of its life. During this phase, the star evolution can be described as a sequence of stationary configurations, which we determine by solving the neutrino transport and the stellar structure equations in general relativity. We include in our model the angular momentum loss due to neutrino emission. We find that the requirement of a rotation rate not exceeding the mass-shedding limit at the beginning of the evolution implies a strict bound on the rotation rate at later times. Moreover, assuming that the proto-neutron star is born with a finite ellipticity, we determine the emitted gravitational wave signal, and estimate its detectability by present and future ground-based interferometric detectors.

Galaxy intrinsic alignment can be a severe source of error in weak-lensing studies. The problem has been widely studied by numerical simulations and with heuristic models, but without a clear theoretical justification of its origin and amplitude. In particular, it is still unclear whether intrinsic alignment of galaxies is dominated by formation and accretion processes or by the effects of the instantaneous tidal field acting upon them. We investigate this question by developing a simple model of intrinsic alignment for elliptical galaxies, based on the instantaneous tidal field. Making use of the galaxy stellar distribution function, we estimate the intrinsic alignment signal and find that although it has the expected dependence on the tidal field, it is too weak to account for the observed signal. This is an indirect validation of the standard view that intrinsic alignment is caused by formation and/or accretion processes.