Long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak, an extension of the two-dimensional Dainotti relation. This 3D relation identifies a GRB fundamental plane whose existence we here confirm. We extend the original analysis with X-ray data untill July 2016 gathering 183 Swift GRBs with afterglow plateaus and known redshifts. We compare several GRB categories, such as shorts with extended emission, SEE, X-ray Flashes, GRBs associated with SNe, a sample of only long-duration GRBs (132), selected from the total sample by excluding GRBs of the previous categories, and the gold sample, composed by GRBs with light curves with good data coverage and relatively flat plateaus. The relation planes for each of these categories are not statistically different from the gold fundamental plane, with the exception of the SEE, which are hence identified as a physically distinct class of objects. The gold fundamental plane has an intrinsic scatter smaller than any plane derived from the other sample categories. Thus, the distance of any particular GRB category from this plane becomes a key parameter. Additionally, we tested this 3D relations by using GRBs observed at high energy, namely the peak luminosity values derived by the Fermi-Gamma Ray Burst Monitor (GBM). The 3D relation is also confirmed for GRBs observed by the GBM, thus showing its independence from the energy range. Furthermore, we computed the several category planes with T*a as a dependent parameter obtaining for each category smaller intrinsic scatters (reaching a reduction of $24\%$ for the long GRBs). The fundamental plane is independent from several prompt and afterglow parameters.

The analysis of 176 gamma ray burst (GRB) afterglow plateaus observed by Swift from GRBs with known redshifts revealed that the subsample of long GRBs associated with supernovae (LONG-SNe) - 19 events - presents a very high correlation coefficient between the luminosity at the end of the plateau phase La and the end time of the plateau T*a, hereafter Dainotti relation. Furthermore, these SNe Ib/c associated with GRBs also obey the peak-magnitude stretch relation, similar to that used to standardize the SNe Ia. We here investigate a category of GRBs with plateau and associated with SNe to compare the Dainotti correlation for this sample with the correlation for long GRBs for which no associated SN has been observed (hereafter LONG-NO-SNe, 128 GRBs) and to check whether there is a difference among these sub-samples. We first adopted a non-parametric statistical method to take redshift evolution into account and to check if and how this effect may steepen the slope for the LONG-NO-SNe sample. Therefore, removing selection bias is the first step before any comparison among samples observed at different redshifts could be properly performed. Then, we applied the T-student test to evaluate a statistical difference between the slopes of the two samples. We demonstrate that there is no evolution for the slope of the LONG-NO-SNe sample and no evolution is expected for the LONG-SNe sample at small redshifts. The difference between the slope of the LONG-NO-SNe and the slope of LONG-SNe with firm spectral detection of SN components, is statistically significant. This possibly suggests that, unlike LONG-NO-SNe, LONG-SNe with firm spectroscopic features of the associated SNe might not require a standard energy reservoir in the plateau phase. Therefore, this analysis may open new perspectives in future theoretical investigations of the GRBs with plateau emission and that are associated with SNe.

A class of long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak. We extended the original analysis with X-ray data from July 2014 to July 2016 achieving a total sample of 183 Swift GRBs with afterglow plateaus and known redshifts. We added the most recent GRBs to the previous 'gold sample' (now including 45 GRBs) and obtained a relation plane with intrinsic scatter compatible within one sigma with the previous result. We compared several GRB categories, such as short with extended emission, X-ray Flashes, GRBs associated with SNe, long-duration GRBs, and the gold sample, composed only by GRBs with light curves with good data coverage and relatively flat plateaus and evaluated their relation planes. We found that they are not statistically different from the fundamental plane derived from the gold sample and that the fundamental plane still has the smallest scatter. We compared the jet opening angles tabulated in literature with the angles derived using the Eiso-Egamma relation of the method in Pescalli et al. (2015) and calculated the relation plane for a sample of long GRBs accounting for the different jet opening angles. We observed that this correction does not significantly reduce the scatter. In an extended analysis, we found that the fundamental plane is independent from several prompt and afterglow parameters.

A class of long Gamma-Ray Bursts (GRBs) presenting light curves with an extended plateau phase in their X-ray afterglows obeys a correlation between the rest frame end time of the plateau, $T_a$, and its corresponding X-ray luminosity, $L_{a}$, Dainotti et al. (2008). In this work we perform an analysis of a total sample of 176 {\it Swift} GRBs with known redshifts, exhibiting afterglow plateaus. By adding a third parameter, that is the peak luminosity in the prompt emission, $L_{peak}$, we discover the existence of a new three parameter correlation. The scatter of data about this plane becomes smaller when a class-specific GRB sample is defined. This sample of 122 GRBs is selected from the total sample by excluding GRBs with associated Supernovae (SNe), X-ray flashes and short GRBs with extended emission. {\bf With this sample the three parameter correlation identifies a GRB `fundamental plane'}. Moreover, we further limit our analysis to GRBs with lightcurves having good data coverage and almost flat plateaus, 40 GRBs forming our `gold sample'. The intrinsic scatter, $\sigma_{int}=0.27 \pm 0.04$, for the three-parameter correlation for this last subclass is more than twice smaller than the value for the $L_{a}-T_a$ one, making this the tightest three parameter correlation involving the afterglow plateau phase. Finally, we also show that a slightly less tight correlation is present between $L_{peak}$ and a proxy for the total energy emitted during the plateau phase, $L_a T_a$, confirming the existence of an energy scaling between prompt and afterglow phases.

Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. In this work, we present proof-of-principle 3D Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). We perform benchmark studies of $^{16} \mathrm{O}$, $^{208} \mathrm{Pb}$ and $^{238} \mathrm{U}$ nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24\:fm with an average density of about $\rho = 0.05 \:\mathrm{fm}^{-3}$, proton fraction of $Y_p = 0.3$ and temperature of $T=0\:$MeV. Our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. In the present pasta calculations spin-orbit forces are not included but will be added in the future. Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.

We study the dark energy equation of state as a function of redshift in a non-parametric way, without imposing any {\it a priori} $w(z)$ (ratio of pressure over energy density) functional form. As a check of the method, we test our scheme through the use of synthetic data sets produced from different input cosmological models which have the same relative errors and redshift distribution as the real data. Using the luminosity-time $L_{X}-T_{a}$ correlation for GRB X-ray afterglows (the Dainotti et al. correlation), we are able to utilize GRB sample from the {\it Swift} satellite as probes of the expansion history of the Universe out to $z \approx 10$. Within the assumption of a flat FLRW universe and combining SNeIa data with BAO constraints, the resulting maximum likelihood solutions are close to a constant $w=-1$. If one imposes the restriction of a constant $w$, we obtain $w=-0.99 \pm 0.06$ (consistent with a cosmological constant) with the present day Hubble constant as $H_{0}=70.0 \pm 0.6$ ${\rm km} \, {\rm s}^{-1} {\rm Mpc}^{-1}$ and density parameter as $\Omega_{\Lambda 0}=0.723 \pm 0.025$, while non-parametric $w(z)$ solutions give us a probability map which is centred at $H_{0}=70.04 \pm 1$ ${\rm km} \, {\rm s}^{-1} {\rm Mpc}^{-1}$ and $\Omega_{\Lambda 0}=0.724 \pm 0.03$. Our chosen GRB data sample with full correlation matrix allows us to estimate the amount, as well as quality (errors) of data, needed to constrain $w(z)$ in the redshift range extending an order of magnitude in beyond the farthest SNeIa measured.

The coefficients of diffusion, thermal conductivity, and shear viscosity are calculated for a system of non-relativistic particles interacting via a delta-shell potential V(r)=-v \delta(r-R) when the average distance between particles is smaller than R. The roles of resonances and long scattering lengths including the unitary limit are examined. Results for ratios of diffusion to viscosity and viscosity to entropy density are presented for varying scattering lengths.

Gravitational waves from the final stages of inspiralling binary neutron stars are expected to be one of the most important sources for ground-based gravitational wave detectors. The masses of the components are determinable from the orbital and chirp frequencies during the early part of the evolution during which tidal effects provide small correction; however, during this phase the signal is relatively clean. The accumulated phase shift due to tidal corrections is characterized by a single quantity, the Love number, which is sensitive to the compactness parameter M/R and the star's internal structure, and its determination could constrain the star's radius. We show that the Love number of normal neutron stars are much different from those of self-bound strange quark matter stars and could therefore provide an important way to distinguish between these two classes of stars.

The shear viscosity of the crust might have a damping effect on the amplitude of r-modes of rotating neutron stars. This damping has implications for the emission of gravitational waves. We calculate the contribution to the shear viscosity coming from the ions using both semi-analytical methods, that consider binary collisions, and Molecular Dynamics simulations. We compare these results with the contribution coming from electrons. We study how the shear viscosity depends on density for conditions of interest in neutron star envelopes and outer crusts. In the low density limit, we find good agreement between results of our molecular dynamics simulations and classical semi-analytic calculations.

A challenge to standard leptonic SSC models are so-called orphan TeV flares, i.e. enhanced very high energy (VHE) gamma-ray emission without any contemporaneous X-ray flaring activity, that have recently been observed in TeV-blazars (e.g., 1ES 1959+650). In order to explain the orphan TeV flare of 1ES 1959+650 observed in June 2002, the co-called hadronic synchrotron mirror model has been developed. Here, relativistic protons are proposed to exist in the jet, and interact with reflected electron-synchrotron radiation of the precursor SSC flare. If the reflector is located in the cloud region, time shifts of several days are possible between the precursor and the orphan flare. The external photons, blueshifted in the comoving jet frame, are able to excite the \Delta(1232)-resonance when interacting with protons of Lorentz factors \gamma'_p~10^3-10^4. The decay products of this resonance include charged pions which, upon decay, give rise to neutrino production during the orphan flare. In this paper we calculate the expected neutrino emission for the June 4, 2002, orphan TeV flare of 1ES 1959+650. We compare our results with the recent observations of AMANDA-II of a neutrino event in spatial and temporal coincidence with the orphan flare of this blazar. We find that the expected neutrino signal from the hadronic synchrotron mirror model is insufficient to explain the observed neutrino event from the direction of 1ES 1959+650.