[Shortened] CONTEXT: [...] AIMS: To identify and verify the BdHNe I properties in the additional sources GRB160509A, GRB160625B and GRB1340427A, and compare and contrast the results with the ones of a BdHN II source GRB180728A. We have also identified in all four sources, following the analysis GRB 130427A in the companion paper, the GeV radiation during and following the UPE phase. Also in all the four sources, we describe the spectral properties of their afterglow emission, including the mass estimate of the $\nu$NS, following the results presented in the companion paper. METHODS: [...] RESULTS: The results of the spectral analysis have validated the common properties in all BdHNe I: the three Episodes as well as the self-similar structures and the associated power-laws in the UPE phase. The profound similarities of the results have made a significant step forward in the taxonomy of GRBs and in evidencing a standard composition of the BdHN I. This opens the opportunity of a vaster inquire of the astrophysical nature of their components in the population synthesis approach: e.g., the BH formation in all BdHN I occurs due to accretion of the SN ejecta in a tight binary system with a neutron star companion which reaches its critical mass, leading to the formation of the BH. The SN-rise in all five BdHNe are compare and contrasted. CONCLUSIONS: The most far reaching discovery of self-similarities and power-laws here extensively confirmed, thanks also to the conclusions presented in the companion papers, leads to the existence of a discrete quantized repetitive polarized emission, both in the GeV and MeV observed by {\it Fermi}-GBM and {\it Fermi}-LAT, on a timescale as short as $10^{-14}$s. These results open new paths in the discovery of fundamental physical laws.

Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimeter broadband flux to total infrared luminosity (L_IR), and hence star formation rate, in high-z galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g., redshift and star formation rate. In this paper, we analyse the dust temperatures of massive (M* > 10^10 Msun) z=2-6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. At z~2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature does not strongly evolve with redshift over z=2-6 at fixed IR luminosity but is tightly correlated with L_IR at fixed z. The mass-weighted temperature is important for accurately estimating the total dust mass. We also analyse an 'equivalent' dust temperature for converting (sub)millimeter flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature ('warmer dust') do not necessarily have higher mass-weighted temperatures. A 'two-phase' picture for interstellar dust can explain the different scaling relations of the various dust temperatures.

We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of halos up to M_halo~10^12 M_sun (M_star~10^10.5 M_sun) at z>=5 from the Feedback in Realistic Environments project. We post-process our simulations with a three-dimensional Monte Carlo dust radiative transfer code to study dust extinction, dust emission, and dust temperature within these simulated z>=5 galaxies. Our sample forms a tight correlation between infrared excess (IRX=F_IR/F_UV) and ultraviolet (UV)-continuum slope (beta_UV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX-beta_UV relation is mainly determined by the shape of the extinction curve and is independent of its normalization (set by the dust-to-gas ratio). The bolometric IR luminosity (L_IR) correlates with the intrinsic UV luminosity and the star formation rate (SFR) averaged over the past 10 Myr. We predict that at a given L_IR, the peak wavelength of the dust spectral energy distributions for z>=5 galaxies is smaller by a factor of 2 (due to higher dust temperatures on average) than at z=0. The higher dust temperatures are driven by higher specific SFRs and SFR surface densities with increasing redshift. We derive the galaxy UV luminosity functions (LFs) at z=5-10 from our simulations and confirm that a heavy attenuation is required to reproduce the observed bright-end UVLFs. We also predict the IRLFs and UV luminosity densities at z=5-10. We discuss the implications of our results on current and future observations probing dust attenuation and emission in z>=5 galaxies.

Recent long wavelength observations on the thermal dust continuum suggest that the Rayleigh-Jeans (RJ) tail can be used as a time-efficient quantitative probe of the dust and ISM mass in high-$z$ galaxies. We use high-resolution cosmological simulations from the Feedback in Realistic Environment (FIRE) project to analyze the dust emission of $M_*>10^{10}\;M_{\odot}$ galaxies at $z=2-4$. Our simulations (MassiveFIRE) explicitly include various forms of stellar feedback, and they produce the stellar masses and star formation rates of high-$z$ galaxies in agreement with observations. Using radiative transfer modelling, we show that sub-millimeter (sub-mm) luminosity and molecular ISM mass are tightly correlated and that the overall normalization is in quantitative agreement with observations. Notably, sub-mm luminosity traces molecular ISM mass even during starburst episodes as dust mass and mass-weighted temperature evolve only moderately between $z=4$ and $z=2$, including during starbursts. Our finding supports the empirical approach of using broadband sub-mm flux as a proxy for molecular gas content in high-$z$ galaxies. We thus expect single-band sub-mm observations with ALMA to dramatically increase the sample size of high-$z$ galaxies with reliable ISM masses in the near future.

The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback.

Two CCD epochs of light minimum and a complete R light curve of SS Ari are presented. The light curve obtained in 2007 was analyzed with the 2003 version of the W-D code. It is shown that SS Ari is a shallow contact binary system with a mass ratio $q=3.25$ and a degree of contact factor f=9.4(\pm0.8%). A period investigation based on all available data shows that there may exist two distinct solutions about the assumed third body. One, assuming eccentric orbit of the third body and constant orbital period of the eclipsing pair results in a massive third body with $M_3=1.73M_{\odot}$ and P_3=87.0$yr. On the contrary, assuming continuous period changes of the eclipsing pair the orbital period of tertiary is 37.75yr and its mass is about $0.278M_{\odot}$. Both of the cases suggest the presence of an unseen third component in the system.

We obtained complete $R$ and $I$ light curves of GSC 02393-00680 in 2008 and analyzed them with the 2003 version of the W-D code. It is shown that GSC 02393-00680 is a W-type shallow contact binary system with a high mass ratio $q=1.600$ and a degree of contact factor $f=5.0%(\pm1.3%)$. It will be a good example to check up on the TRO theory. A period investigation based on all available data suggests that the system has a small-amplitude period oscillation ($A_3=0.^{d}0030$; $T_3=1.92$years). This may indicate it has a moderate mass close third body, which is similar to XY Leo.

The CCD photometric data of the EW-type binary, II CMa, which is a contact star in the field of the middle-aged open cluster Berkeley 33, are presented. The complete R light curve was obtained. In the present paper, using the five CCD epochs of light minimum (three of them are calculated from Mazur et al. (1993)'s data and two from our new data), the orbital period P was revised to 0.22919704 days. The complete R light curve was analyzed by using the 2003 version of W-D (Wilson-Devinney) program. It is found that this is a contact system with a mass ratio $q=0.9$ and a contact factor $f=4.1%$. The high mass ratio ($q=0.9$) and the low contact factor ($f=4.1%$) indicate that the system just evolved into the marginal contact stage.

TX Cnc is a member of the young open cluster NGC 2632. In the present paper, four CCD epochs of light minimum and a complete V light curve of TX Cnc are presented. A period investigation based on all available photoelectric or CCD data showed that it is found to be superimposed on a long-term increase ($dP/dt=+3.97\times{10^{-8}}$\,days/year), and a weak evidence suggests that it includes a small-amplitude period oscillation ($A_3=0.^{d}0028$; $T_3=26.6\,years$). The light curves in the V band obtained in 2004 were analyzed with the 2003 version of the W-D code. It was shown that TX Cnc is an overcontact binary system with a degree of contact factor $f=24.8%(\pm0.9%)$. The absolute parameters of the system were calculated: $M_1=1.319\pm0.007M_{\odot}$, $M_2=0.600\pm0.01M_{\odot}$; $R_1=1.28\pm0.19R_{\odot}$, $R_2=0.91\pm0.13R_{\odot}$. TX Cnc may be on the TRO-controlled stage of the evolutionary scheme proposed by Qian (2001a, b; 2003a), and may contains an invisible tertiary component ($m_3\approx0.097M_{\odot}$). If this is true, the tertiary component has played an important role in the formation and evolution of TX Cnc by removing angular momentum from the central system(Pribulla & Rucinski, 2006). In this way the contact binary configuration can be formed in the short life time of a young open cluster via AML.

This paper analyzes the first secured four color light curves of V396 Mon using the 2003 version of the WD code. It is confirmed that V396 Mon is a shallow W-type contact binary system with a mass ratio $q=2.554(\pm0.004)$ and a degree of contact factor $f=18.9%(\pm1.2%)$. A period investigation based on all available data shows that the period of the system includes a long-term decrease ($dP/dt=-8.57\times{10^{-8}}$ days/year) and an oscillation ($A_3=0.^{d}0160$; $T_3=42.4\,years$). They are caused by angular momentum loss (AML) and light-time effect, respectively. The suspect third body perhaps is a small M-type star (about 0.31 solar mass). Though some proofs show that this system has strong magnetic activity, through analyzing we found that the Applegate mechanism cannot explain the periodic changes. This binary is an especially important system according to Qian's statistics of contact binaries as its mass ratio lies near the proposed pivot point about which the physical structure of contact binaries supposedly oscillate.