Infrared imaging of the colliding-wind binary Apep has revealed a spectacular dust plume with complicated internal dynamics that challenges standard colliding-wind binary physics. Such challenges can be potentially resolved if a rapidly-rotating Wolf-Rayet star is located at the heart of the system, implicating Apep as a Galactic progenitor system to long-duration gamma-ray bursts. One of the difficulties in interpreting the dynamics of Apep is that the spectral composition of the stars in the system was unclear. Here we present visual to near-infrared spectra that demonstrate that the central component of Apep is composed of two classical Wolf-Rayet stars of carbon- (WC8) and nitrogen-sequence (WN4-6b) subtypes. We argue that such an assignment represents the strongest case of a classical WR+WR binary system in the Milky Way. The terminal line-of-sight wind velocities of the WC8 and WN4-6b stars are measured to be $2100 \pm 200$ and $3500 \pm 100$ km s$^{-1}$, respectively. If the mass-loss rate of the two stars are typical for their spectral class, the momentum ratio of the colliding winds is expected to be $\approx$ 0.4. Since the expansion velocity of the dust plume is significantly smaller than either of the measured terminal velocities, we explore the suggestion that one of the Wolf-Rayet winds is anisotropic. We can recover a shock-compressed wind velocity consistent with the observed dust expansion velocity if the WC8 star produces a significantly slow equatorial wind with a velocity of $\approx$530 km s$^{-1}$. Such slow wind speeds can be driven by near-critical rotation of a Wolf-Rayet star.

Because the binary neutron star (BNS) merger GW170817 and its electromagnetic (EM) counterparts were successfully detected, there is a chance to explore the joint effect of the host galaxy and our Milky Way on the weak equivalence principle (WEP) test. In this paper, through using the Navarro-Frenk-White profile and the Herquist profile, we present a new analytic model to calculate the galactic gravitational potential, in which the possible locations of the source by the observed angle offset and the second SN (SN2) kick is accounted for. We show that the upper limit of $\Delta \gamma$ is $10^{-9}$ for the comparison of GW170817/GRB 170817A, and is $10^{-4}$ for the comparison of GW170817/AT2017gfo, which are tighter by one or two orders of magnitude than these of only by the measured Milky Way gravitational potential in literatures. We demonstrate that the WEP test could be strengthened by the contribution of the host galaxy NGC 4993 to the Shapiro time delay. Meanwhile, we also find the large natal kick can produce a maximum deviation of about $20\%$ on the results with a typical kick velocity $400\sim 500 km s^{-1}$. Finally, we study and analyse the impact from the halo mass of NGC 4993 with a typical 0.2 dex uncertainty, and find the upper limit of $\Delta \gamma$ of the maximum mass $10^{12.4}h^{-1} M_{\odot}$ is nearly two times tighter than that of the minimum mass $10^{12.0}h^{-1} M_{\odot}$. It suggests the constraint on WEP test would be also enhanced when the halo mass increases in the host galaxy NGC 4993.

We present ALMA observations and multiwavelength spectral energy distribution (SED) analysis in a WISE-selected, hyperluminous dust-obscured quasar W0533-3401 at $z=2.9$. We derive its physical properties of each component, such as molecular gas, stars, dust and the central supermassive black hole (SMBH). Both the dust continuum at 3 mm and the CO(3-2) line are detected. The derived molecular gas mass $M_{\rm gas}=8.4\times10^{10}\ M_\odot$ and its fraction $f_{\rm gas}=0.7$ suggest that W0533-3401 is gas-rich. The star formation rate (SFR) has been estimated to be $\sim3000-7000\ M_\odot$ yr$^{-1}$ by using different methods. The high values of SFR and specific SFR suggest that W0533-3401 is a maximum-starburst. The corresponding gas depletion timescales are very short ($t_{\rm depl}\sim12-28$ Myr). The CO(3-2) emission line is marginally resolved and has a velocity gradient, which is possibly due to a rotating gas disk, gas outflow or merger. Finally, we infer the black hole mass growth rate of W0533-3401 (${\dot{M}}_{\rm BH}$ = 49 $M_\odot$ yr$^{-1}$), which suggests a rapid growth of the central SMBH. The observed black hole to stellar mass ratio $M_{\rm BH}/M_\star$ of W0533-3401, which is dependent on the adopted Eddington ratio, is over one order of magnitude higher than the local value, and is evolving towards the evolutionary trend of unobscured quasars. Our results are consistent with the scenario that \wobs, with both a gas-rich maximum-starburst and a rapid black hole growth, is experiencing a short transition phase towards an unobscured quasar.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

Fitting the multi-wavelength spectral energy distributions (SEDs) of galaxies is a widely used technique to extract information about the physical properties of galaxies. However, a major difficulty lies in the numerous uncertainties regarding almost all ingredients of the SED modeling of galaxies. The Bayesian methods provide a consistent conceptual basis for dealing with the problem of inference with many uncertainties. While the Bayesian parameter estimation method have become quite popular in the field of SED fitting of galaxies, the Bayesian model comparison method, which is based on the same Bayes' rule, is still not widely used in this field. With the application of Bayesian model comparison method in a series of papers, we show that the results obtained with Bayesian model comparison are understandable in the context of stellar/galaxy physics. These results indicate that Bayesian model comparison is a reliable and very powerful method for the SED fitting of galaxies.

When modeling and interpreting the spectral energy distributions (SEDs) of galaxies, the simple stellar population (SSP) model, star formation history (SFH) and dust attenuation law (DAL) are three of the most important components. However, each of them carries significant uncertainties which have seriously limited our ability to reliably recover the physical properties of galaxies from the analysis of their SEDs. In this paper, we present a Bayesian framework to deal with these uncertain components simultaneously. Based on the Bayesian evidence, a quantitative implement of the principle of Occam's razor, the method allows a more objective and quantitative discrimination among the different assumptions about these uncertain components. With a Ks-selected sample of 5467 low-redshift (mostly with $z\lesssim 1$) galaxies in the COSMOS/UltraVISTA field and classified into passively evolving galaxies (PEGs) and star-forming galaxies (SFGs) with UVJ diagram, we present a Bayesian discrimination of a set of 16 SSP models from five research groups (BC03 and CB07, M05, GALEV, Yunnan-II, BPASS V2.0), five forms of SFH (Burst, Constant, Exp-dec, Exp-inc, Delayed-$\tau$), and four kinds of DAL (Calzetti law, MW, LMC, SMC). We show that the results obtained with the method are either obvious or understandable in the context of galaxy physics. We conclude that the Bayesian model comparison method, especially that for a sample of galaxies, is very useful for discriminating the different assumptions in the SED modeling of galaxies. The new version of the BayeSED code, which is used in this work, is publicly available at https://bitbucket.org/hanyk/bayesed/.

We present a multiwavelength study of an atypical submillimeter galaxy in the GOODS-North field, with the aim to understand its physical properties of stellar and dust emission, as well as the central AGN activity. Although it is shown that the source is likely an extremely dusty galaxy at high redshift, its exact position of submillimeter emission is unknown. With the new NOEMA interferometric imaging, we confirm that the source is a unique dusty galaxy. It has no obvious counterpart in the optical and even NIR images observed with HST at lambda~<1.4um. Photometric-redshift analyses from both stellar and dust SED suggest it to likely be at z~>4, though a lower redshift at z~>3.1 cannot be fully ruled out (at 90% confidence interval). Explaining its unusual optical-to-NIR properties requires an old stellar population (~0.67 Gyr), coexisting with a very dusty ongoing starburst component. The latter is contributing to the FIR emission, with its rest-frame UV and optical light being largely obscured along our line of sight. If the observed fluxes at the rest-frame optical/NIR wavelengths were mainly contributed by old stars, a total stellar mass of ~3.5x10^11Msun would be obtained. An X-ray spectral analysis suggests that this galaxy harbors a heavily obscured AGN with N_H=3.3x10^23 cm^-2 and an intrinsic 2-10 keV luminosity of L_X~2.6x10^44 erg/s, which places this object among distant type 2 quasars. The radio emission of the source is extremely bright, which is an order of magnitude higher than the star-formation-powered emission, making it one of the most distant radio-luminous dusty galaxies. The combined characteristics of the galaxy suggest that the source appears to have been caught in a rare but critical transition stage in the evolution of submillimeter galaxies, where we are witnessing the birth of a young AGN and possibly the earliest stage of its jet formation and feedback.

We develop a formalism to help calculate in quantum field theory the departures from the description of a system by classical field equations. We apply the formalism to a homogeneous condensate with attractive contact interactions and to a homogeneous self-gravitating condensate in critical expansion. In their classical descriptions, such condensates persist forever. We show that in their quantum description, parametric resonance causes quanta to jump in pairs out of the condensate into all modes with wavector less than some critical value. We calculate in each case the time scale over which the homogeneous condensate is depleted, and after which a classical description is invalid. We argue that the duration of classicality of inhomogeneous condensates is shorter than that of homogeneous condensates.

Dark matter axions and other highly degenerate bosonic fluids are commonly described by classical field equations. In a recent paper \cite{BECprop} we calculated the duration of classicality of homogeneous condensates with attractive contact interactions and of self-gravitating homogeneous condensates in critical expansion. According to their classical equations of motion, such condensates persist forever. In their quantum evolution parametric resonance causes quanta to jump in pairs out of the condensate into all modes with wavevector less than some critical value. We estimated in each case the time scale over which the condensate is depleted and after which a classical description is invalid.

Broad-lined type Ic supernovae (SNe Ic-BL) are a subclass of rare core collapse SNe whose energy source is debated in the literature. Recently a series of investigations on SNe Ic-BL with the magnetar (plus 56Ni) model were carried out. Evidence for magnetar formation was found for the well-observed SNe Ic-BL 1998bw and 2002ap. In this paper we systematically study a large sample of SNe Ic-BL not associated with gamma-ray bursts. We use photospheric velocity data determined in a homogeneous way. We find that the magnetar+56Ni model provides a good description of the light curves and velocity evolution of our sample of SNe Ic-BL, although some SNe (not all) can also be described by the pure-magnetar model or by the two-component pure-56Ni model (3 out of 12 are unlikely explained by two-component model). In the magnetar+56Ni model, the amount of 56Ni required to explain their luminosity is significantly reduced, and the derived initial explosion energy is, in general, in accordance with neutrino heating. Some correlations between different physical parameters are evaluated and their implications regarding magnetic field amplification and the total energy reservoir are discussed.

Hot dust-obscured galaxies (Hot DOGs) are a new population recently discovered in the \wise All-Sky survey. Multiwavelength follow-up observations suggest that they are luminous, dust-obscured quasars at high redshift. Here we present the JCMT SCUBA-2 850 $\mu m$ follow-up observations of 10 Hot DOGs. Four out of ten Hot DOGs have been detected at $>3\sigma$ level. Based on the IR SED decomposition approach, we derive the IR luminosities of AGN torus and cold dust components. Hot DOGs in our sample are extremely luminous with most of them having $L_{\rm IR}^{\rm tot}>10^{14} L_\odot$. The torus emissions dominate the total IR energy output. However, the cold dust contribution is still non-negligible, with the fraction of the cold dust contribution to the total IR luminosity $(\sim 8-24\%)$ being dependent on the choice of torus model. The derived cold dust temperatures in Hot DOGs are comparable to those in UV bright quasars with similar IR luminosity, but much higher than those in SMGs. Higher dust temperatures in Hot DOGs may be due to the more intense radiation field caused by intense starburst and obscured AGN activities. Fourteen and five submillimeter serendipitous sources in the 10 SCUBA-2 fields around Hot DOGs have been detected at $>3\sigma$ and $>3.5\sigma$ levels, respectively. By estimating their cumulative number counts, we confirm the previous argument that Hot DOGs lie in dense environments. Our results support the scenario in which Hot DOGs are luminous, dust-obscured quasars lying in dense environments, and being in the transition phase between extreme starburst and UV-bright quasars.

Broad-lined type Ic supernovae (SNe Ic-BL) are peculiar stellar explosions that distinguish themselves from ordinary SNe. Some SNe Ic-BL are associated with long-duration (\gtrsim 2 s) gamma-ray bursts (GRBs). Black holes and magnetars are two types of compact objects that are hypothesized to be central engines of GRBs. In spite of decades of investigations, no direct evidence for the formation of black holes or magnetars has been found for GRBs so far. Here we report the finding that the early peak (t \lesssim 50 days) and late-time (t \gtrsim 300 days) slow decay displayed in the light curves of both SNe 1998bw (associated with GRB 980425) and 2002ap (not GRB-associated) can be attributed to magnetar spin-down with initial rotation period P0 \sim 20 ms, while the intermediate-time (50 \lesssim t \lesssim 300 days) linear decline is caused by radioactive decay of 56Ni. The connection between the early peak and late-time slow decline in the light curves is unexpected in alternative models. We thus suggest that GRB 980425 and SN 2002ap were powered by magnetars.

Broad-lined type Ic supernovae (SNe Ic-BL) are of great importance because their association with long-duration gamma-ray bursts (LGRBs) holds the key to deciphering the central engine of LGRBs, which refrains from being unveiled despite decades of investigation. Among the two popularly hypothesized types of central engine, i.e., black holes and strongly magnetized neutron stars (magnetars), there is mounting evidence that the central engine of GRB-associated SNe (GRB-SNe) is rapidly rotating magnetars. Theoretical analysis also suggests that magnetars could be the central engine of SNe Ic-BL. What puzzled the researchers is the fact that light curve modeling indicates that as much as 0.2-0.5 solar mass of 56Ni was synthesized during the explosion of the SNe Ic-BL, which is unfortunately in direct conflict with current state-of-the-art understanding of magnetar-powered 56Ni synthesis. Here we propose a dynamic model of magnetar-powered SNe to take into account the acceleration of the ejecta by the magnetar, as well as the thermalization of the injected energy. Assuming that the SN kinetic energy comes exclusively from the magnetar acceleration, we find that although a major fraction of the rotational energy of the magnetar is to accelerate the SNe ejecta, a tiny fraction of this energy deposited as thermal energy of the ejecta is enough to reduce the needed 56Ni to 0.06 solar mass for both SN 1997ef and SN 2007ru. We therefore suggest that magnetars could power SNe Ic-BL both in aspects of energetics and of 56Ni synthesis.

Previous studies have shown that WISE-selected hyperluminous, hot dust-obscured galaxies (Hot DOGs) are powered by highly dust-obscured, possibly Compton-thick AGNs. High obscuration provides us a good chance to study the host morphology of the most luminous AGNs directly. We analyze the host morphology of 18 Hot DOGs at $z\sim3$ using Hubble Space Telescope/WFC3 imaging. We find that Hot DOGs have a high merger fraction ($62\pm 14 \%$). By fitting the surface brightness profiles, we find that the distribution of S\'ersic indices in our Hot DOG sample peaks around 2, which suggests that most of Hot DOGs have transforming morphologies. We also derive the AGN bolometric luminosity ($\sim10^{14}L_\odot$) of our Hot DOG sample by using IR SEDs decomposition. The derived merger fraction and AGN bolometric luminosity relation is well consistent with the variability-based model prediction (Hickox et al. 2014). Both the high merger fraction in IR-luminous AGN sample and relatively low merger fraction in UV/optical-selected, unobscured AGN sample can be expected in the merger-driven evolutionary model. Finally, we conclude that Hot DOGs are merger-driven and may represent a transit phase during the evolution of massive galaxies, transforming from the dusty starburst dominated phase to the unobscured QSO phase.

We utilize a Bayesian approach to fit the observed mid-IR-to-submm/mm spectral energy distributions (SEDs) of 22 WISE-selected and submm-detected, hyperluminous hot dust-obscured galaxies. By adopting the Torus+GB model, we decompose the observed IR SEDs of Hot DOGs into torus and cold dust components. The main results are: 1) Hot DOGs in our submm-detected sample are hyperluminous, with torus emission dominating the IR energy output. However, cold dust emission is non-negligible, averagely contributing ~24% of total IR luminosity. 2) Compared to QSO and starburst SED templates, the median SED of Hot DOGs shows the highest luminosity ratio between mid-IR and submm at rest-frame, while it is very similar to that of QSOs at 10-50um suggesting that the heating sources of Hot DOGs should be buried AGNs. 3) Hot DOGs have both high dust temperatures ~73K and IR luminosity of cold dust. The T-L relation of Hot DOGs suggests that the increase in IR luminosity for Hot DOGs is mostly due to the increase of the dust temperature, rather than dust mass. Hot DOGs have lower dust masses than those of submillimeter galaxies (SMGs) and QSOs within the similar redshift range. Both high IR luminosity of cold dust and relatively low dust mass in Hot DOGs can be expected by their relatively high dust temperatures. 4) Hot DOGs have high dust covering factors, which deviate the previously proposed trend of the dust covering factor decreasing with increasing bolometric luminosity. Finally, we can reproduce the observed properties in Hot DOGs by employing a physical model of galaxy evolution. The result suggests that Hot DOGs may lie at or close to peaks of both star formation and black hole growth histories, and represent a transit phase during the evolution of massive galaxies, transforming from the dusty starburst dominated phase to the optically bright QSO phase. (abridged)

Millisecond magnetars can be formed via several channels: core-collapse of massive stars, accretion-induced collapse of white dwarfs (WDs), double WD mergers, double neutron star (NS) mergers, and WD-NS mergers. Because the mass of ejecta from these channels could be quite different, their light curves are also expected to be diverse. We evaluate the dynamic evolution of optical transients powered by millisecond magnetars. We find that the magnetar with short spin-down timescale converts its rotational energy mostly into the kinetic energy of the transient, while the energy of a magnetar with long spin-down timescale goes into radiation of the transient. This leads us to speculate that hypernovae could be powered by magnetars with short spin-down timescales. At late times the optical transients will gradually evolve into a nebular phase because of the photospheric recession. We treat the photosphere and nebula separately because their radiation mechanisms are different. In some cases the ejecta could be light enough that the magnetar can accelerate it to a relativistic speed. It is well known that the peak luminosity of a supernova (SN) occurs when the luminosity is equal to the instantaneous energy input rate, as shown by Arnett (1979). We show that photospheric recession and relativistic motion can modify this law. The photospheric recession always leads to a delay of the peak time $t_{\mathrm{pk}}$ relative to the time $t_{\times }$ at which the SN luminosity equals the instantaneous energy input rate. Relativistic motion, however, may change this result significantly.

We present a newly developed version of BayeSED, a general Bayesian approach to the spectral energy distribution (SED) fitting of galaxies. The new BayeSED code has been systematically tested on a mock sample of galaxies. The comparison between estimated and inputted value of the parameters show that BayeSED can recover the physical parameters of galaxies reasonably well. We then applied BayeSED to interpret the SEDs of a large Ks-selected sample of galaxies in the COSMOS/UltraVISTA field with stellar population synthesis models. With the new BayeSED code, a Bayesian model comparison of stellar population synthesis models has been done for the first time. We found that the model by Bruzual & Charlot (2003), statistically speaking, has larger Bayesian evidence than the model by Maraston (2005) for the Ks-selected sample. Besides, while setting the stellar metallicity as a free parameter obviously increases the Bayesian evidence of both models, varying the IMF has a notable effect only on the Maraston (2005) model. Meanwhile, the physical parameters estimated with BayeSED are found to be generally consistent with those obtained with the popular grid-based FAST code, while the former exhibits more natural distributions. Based on the estimated physical parameters of galaxies in the sample, we qualitatively classified the galaxies in the sample into five populations that may represent galaxies at different evolution stages or in different environments. We conclude that BayeSED could be a reliable and powerful tool for investigating the formation and evolution of galaxies from the rich multi-wavelength observations currently available. A binary version of MPI parallelized BayeSED code is publicly available at https://bitbucket.org/hanyk/bayesed.

We present a detailed comparison between the 2-10 keV hard X-ray and infrared (IR) luminosity function (LF) of active galactic nuclei (AGN). The composite X-ray to IR spectral energy distributions (SEDs) of AGN used for connecting the hard X-ray LF (HXLF) and IR LF (IRLF) are modeled with a simple but well tested torus model based on the radiative transfer and photoionization code CLOUDY. Four observational determinations of the evolution of 2-10 keV HXLF and six evolution models of the obscured type-2 AGN fraction ($f_2$) have been considered. The 8.0 and 15 \micron LFs for the total, unobscured type-1 and obscured type-2 AGN are predicted from the HXLFs, and then compared with the measurements currently available. We find that the IRLFs predicted from HXLFs tend to underestimate the number of the most IR-luminous AGN. This is independent of the choices of HXLF and $f_2$, and even more obvious for the HXLFs recently measured. We show that the discrepancy between the HXLFs and IRLFs can be largely resolved when the anticorrelation between the UV to X-ray slope $\alpha_{\mathrm{ox}}$ and UV luminosity $L_{\rm UV}$ is appropriately considered. We also discuss other possible explanations for the discrepancy, such as the missing population of Compton-thick AGN and possible contribution of star-formation in the host to the mid-IR. Meanwhile, we find that the HXLFs and IRLFs of AGN can be more consistent with each other if the obscuration mechanisms of quasars and Seyferts are assumed to be different, corresponding to their different triggering and fueling mechanisms. More accurate measurements of the IRLFs of AGN, especially that determined at smaller redshift bins and more accurately separated to that for type-1 and type-2, are very helpful for clarifying these interesting issues.

We present BayeSED, a general purpose tool for doing Bayesian analysis of SEDs by using whatever pre-existing model SED libraries or their linear combinations. The artificial neural networks (ANNs), principal component analysis (PCA) and multimodal nested sampling (MultiNest) techniques are employed to allow a highly efficient sampling of posterior distribution and the calculation of Bayesian evidence. As a demonstration, we apply this tool to a sample of hyperluminous infrared galaxies (HLIRGs). The Bayesian evidences obtained for a pure Starburst, a pure AGN, and a linear combination of Starburst+AGN models show that the Starburst+AGN model have the highest evidence for all galaxies in this sample. The Bayesian evidences for the three models and the estimated contributions of starburst and AGN to infrared luminosity show that HLIRGs can be classified into two groups: one dominated by starburst and the other dominated by AGN. Other parameters and corresponding uncertainties about starburst and AGN are also estimated by using the model with the highest Bayesian evidence. We found that the starburst region of the HLIRGs dominated by starburst tends to be more compact and has a higher fraction of OB star than that of HLIRGs dominated by AGN. Meanwhile, the AGN torus of the HLIRGs dominated by AGN tend to be more dusty than that of HLIRGs dominated by starburst. These results are consistent with previous researches, but need to be tested further with larger samples. Overall, we believe that BayeSED could be a reliable and efficient tool for exploring the nature of complex systems such as dust-obscured starburst-AGN composite systems from decoding their SEDs.

Type Ia supernovae (SNe Ia) play an important role in diverse areas of astrophysics, from the chemical evolution of galaxies to observational cosmology. However, the nature of the progenitors of SNe Ia is still unclear. In this paper, according to a detailed binary population synthesis study, we obtained SN Ia birthrates and delay times from different progenitor models, and compared them with observations. We find that the Galactic SN Ia birthrate from the double-degenerate (DD) model is close to those inferred from observations, while the birthrate from the single-degenerate (SD) model accounts for only about 1/2-2/3 of the observations. If a single starburst is assumed, the distribution of the delay times of SNe Ia from the SD model is a weak bimodality, where the WD + He channel contributes to the SNe Ia with delay times shorter than 100Myr, and the WD + MS and WD + RG channels to those with age longer than 1Gyr.