We suggest in this note that "spider" systems are the naturally expected progenitors of the highest neutron star masses, and possibly low-mass black holes, based on their long-term evolutionary features and actual mass measurements.

An approach to estimate the spatial curvature $\Omega_k$ from data independently of dynamical models is suggested, through kinematic parametrizations: at first for comoving distance ($D_{C}(z)$) with third degree polynomial; at second for Hubble parameter ($H(z)$) with a second degree polynomial and at third for deceleration parameter ($q(z)$) with first order polynomial. All these parametrizations were done as function of redshift $z$. We used SNe Ia dataset from Pantheon compilation with 1048 apparent magnitudes estimated in the range $0.01<z<2.3$ with systematic and statistical errors and a compilation of 51 estimated $H(z)$ data. The spatial curvature found for $H(z)$ parametrization was $\Omega_{k}=-0.03^{+0.21+0.46}_{-0.24-0.45}$. Comoving distance $D_{C}$ parametrization has provided $\Omega_{k}=0.11^{+0.21+0.48}_{-0.24-0.44}$ and the parametrization for deceleration parameter $q(z)$ resulted in $\Omega_{k}=-0.05^{+0.21+0.48}_{-0.25-0.45}$. The three results are self-consistent and fully compatible with the spatially flat Universe as predicted by many inflation models and data from CMB. However our results do not discard the recent possibility of a closed Universe with $-0.095<\Omega_k<-0.007$ found by Planck Legacy 2018 (PL2018). This type of analysis may be interesting as it avoids any bias because it does not depend on assumptions about the matter content for estimating $\Omega_k$.

This paper aims to put constraints on the transition redshift $z_t$, which determines the onset of cosmic acceleration, in cosmological-model independent frameworks. In order to do that, we use the non-parametric Gaussian Process method with $H(z)$ and SNe Ia data. The deceleration parameter reconstruction from $H(z)$ data yields $z_t=0.59^{+0.12}_{-0.11}$. The reconstruction from SNe Ia data assumes spatial flatness and yields $z_t=0.683^{+0.11}_{-0.082}$. These results were found with a Gaussian kernel and we show that they are consistent with two other kernel choices.

We compile 41 $H(z)$ data from literature and use them to constrain O$\Lambda$CDM and flat $\Lambda$CDM parameters. We show that the available $H(z)$ suffers from uncertainties overestimation and propose a Bayesian method to reduce them. As a result of this method, using $H(z)$ only, we find, in the context of O$\Lambda$CDM, $H_0=69.5\pm2.5\mathrm{\,km\,s^{-1}Mpc^{-1}}$, $\Omega_m=0.242\pm0.036$ and $\Omega_\Lambda=0.68\pm0.14$. In the context of flat $\Lambda$CDM model, we have found $H_0=70.4\pm1.2\mathrm{\,km\,s^{-1}Mpc^{-1}}$ and $\Omega_m=0.256\pm0.014$. This corresponds to an uncertainty reduction of up to 30\% when compared to the uncorrected analysis in both cases.

The SN1987A in the Giant Magellanic Cloud was an amazing and extraordinary event because it was detected in real time for different neutrinos experiments ($\nu$s) around the world. Approximate $\sim25$ events were observed in three different experiments: Kamiokande II (KII) $\sim 12$, Irvine-Michigan-Brookhaven (IMB) $\sim 8$ e Baksan $\sim 5$, plus a contrived burst at Mont Blanc (Liquid Scintillator Detector - LSD) later dismissed because of energetic requirements (Aglietta et al. 1988). The neutrinos have an important play role into the neutron star newborn: at the moment when the supernova explodes the compact object remnant is freezing by neutrinos ($\sim99\%$ energy is lost in the few seconds of the explosion). The work is motivated by neutrinos' event in relation arrival times where there is a temporal gap between set of events ($\sim6\mbox{s}$). The first part of dataset came from the ordinary mechanism of freezing and the second part suggests different mechanism of neutrinos production. We tested two models of cooling for neutrinos from SN1987A: 1st an exponential cooling is an ordinary model of cooling and 2nd a two-step temperature model that it considers two bursts separated with temporal gap. Our analysis was done with Bayesian tools ({\it Bayesian Information Criterion} - BIC) The result showed strong evidence in favor of a two-step model against one single exponential cooling ($\ln\mbox{B}_{ij} > 5.0$), and suggests the existence of two neutrino bursts at the moment the neutron star was born.

Creation of Cold Dark Matter (CCDM), in the context of Einstein Field Equations, leads to negative creation pressure, which can be used to explain the accelerated expansion of the Universe. In this work we tested six different spatially flat models for matter creation using statistical tools, at light of SN Ia data: Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC) and Bayesian Evidence (BE). These approaches allow to compare models considering goodness of fit and number of free parameters, penalizing excess of complexity. We find that JO model is slightly favoured over LJO/$\Lambda$CDM model, however, neither of these, nor $\Gamma=3\alpha H_0$ model can be discarded from the current analysis. Three other scenarios are discarded either from poor fitting, either from excess of free parameters.

We present in this article an overview of the problem of neutron star masses. After a brief appraisal of the methods employed to determine the masses of neutron stars in binary systems, the existing sample of measured masses is presented, with a highlight on some very well-determined cases. We discuss the analysis made to uncover the underlying distribution and a few robust results that stand out from them. The issues related to some particular groups of neutron stars originated from different channels of stellar evolution are shown. Our conclusions are that last century's paradigm that there a single, $1.4 M_{\odot}$ scale is too simple. A bimodal or even more complex distribution is actually present. It is confirmed that some neutron stars have masses of $\sim 2 M_{\odot}$, and, while there is still no firm conclusion on the maximum and minimum values produced in nature, the field has entered a mature stage in which all these and related questions can soon be given an answer.

The distribution of masses for neutron stars is analyzed using the Bayesian statistical inference, evaluating the likelihood of proposed gaussian peaks by using fifty-four measured points obtained in a variety of systems. The results strongly suggest the existence of a bimodal distribution of the masses, with the first peak around $1.37 {M_{\odot}}$, and a much wider second peak at $1.73 {M_{\odot}}$. The results support earlier views related to the different evolutionary histories of the members for the first two peaks, which produces a natural separation (even if no attempt to "label" the systems has been made here), and argues against the single-mass scale viewpoint. The bimodal distribution can also accommodate the recent findings of $\sim M_{\odot}$ masses quite naturally. Finally, we explore the existence of a subgroup around $1.25 {M_{\odot}}$, finding weak, if any, evidence for it. This recently claimed low-mass subgroup, possibly related to $O-Mg-Ne$ core collapse events, has a monotonically decreasing likelihood and does not stand out clearly from the rest of the sample.