Assuming that the differential rotation of the massive neutron star (NS) formed in the double NS (DNS) mergers has been effectively terminated by the magnetic braking and a uniform rotation has been subsequently established (i.e., a supramassive NS is formed), we analytically derive in this work an approximated expression for the critical total gravitational mass ($M_{\rm tot,c}$) to form supramassive NS (SMNS) in the DNS mergers, benefited from some equation of state (EoS) insensitive relationships. The maximum gravitational mass of the nonrotating NSs ($M_{\rm TOV}$) as well as the dimensionless angular momentum of the remnant ($j$) play the dominant roles in modifying $M_{\rm tot,c}$, while the radius and mass differences of the premerger NSs do not. The GW170817/GRB 170817A/AT2017gfo observations have provided so far the best opportunity to quantitatively evaluate $M_{\rm TOV}$. Supposing the central engine for GRB 170817A is a black hole quickly formed in the collapse of an SMNS, we find $M_{\rm TOV}=2.13^{+0.09}_{-0.08}M_\odot$ (68.3% credibility interval, including also the uncertainties of the EoS insensitive relationships), which is consistent with the constraints set by current NS mass measurements.

In this work we parameterize the Equation of State of dense neutron star (NS) matter with four pressure parameters of $\{\hat{p}_1, \hat{p}_2, \hat{p}_3, \hat{p}_4\}$ and then set the combined constraints with the data of GW 170817 and the data of 6 Low Mass X-ray Binaries (LMXBs) with thermonuclear burst or alternatively the symmetry energy of the nuclear interaction. We find that the nuclear data effectively narrow down the possible range of $\hat{p}_1$, the gravitational wave data plays the leading role in bounding $\hat{p}_2$, and the LMXB data as well as the lower bound on maximal gravitational mass of non-rotating NSs govern the constraints on $\hat{p}_3$ and $\hat{p}_4$. Using posterior samples of pressure parameters and some universal relations, we further investigate how the current data sets can advance our understanding of tidal deformability ($\Lambda$), moment of inertia ($I$) and binding energy ($BE$) of NSs. For a canonical mass of $1.4M_\odot$, we have $I_{1.4} = {1.43}^{+0.30}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 390_{-210}^{+280}$ , $R_{1.4} = 11.8_{-0.7}^{+1.2}~{\rm km}$ and $BE_{1.4} = {0.16}^{+0.01}_{-0.02} M_{\odot}$ if the constraints from the nuclear data and the gravitational wave data have been jointly applied. For the joint analysis of gravitational wave data and the LMXB data, we have $I_{1.4} = {1.28}^{+0.15}_{-0.08} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 220_{-90}^{+90}$, $R_{1.4} = 11.1_{-0.6}^{+0.7}~{\rm km}$ and $BE_{1.4} = {0.18}^{+0.01}_{-0.01} M_{\odot}$. These results suggest that the current constraints on $\Lambda$ and $R$ still suffer from significant systematic uncertainties, while $I_{1.4}$ and $BE_{1.4}$ are better constrained.

As revealed recently by the modeling of the multi-wavelength data of the emission following GW170817/GRB 170817A, there was an off-axis energetic relativistic outflow component launched by this historic double neutron star merger event. In this work we use the results of these modeling to examine the energy extraction process of the central engine. We show that the magnetic process (i.e., the Blandford-Znajek mechanism) is favored, while the neutrino process usually requires a too massive accretion disk if the duration of the central engine activity is comparable to the observed $T_{90}$ of GRB 170817A, unless the timescale of the central engine activity is less than $\sim$ 0.2s. We propose that the GRB observations are helpful to constrain the combined tidal parameter $\tilde{\Lambda}$, and by adopting the accretion disk mass distribution estimated in BZ mechanism, the $90\%$ credible interval of $\tilde{\Lambda}$ for the progenitor of GW170817 is inferred as $309-954$.