Optical observations of a sample of 12 $gamma$-ray bright blazars from four optical data archives, AAVSO, SMARTS, Catalina, and Steward Observatory, are compiled to create densely sampled light curves spanning more than a decade. As a part of the bla
zar multi-wavelength studies, several methods of analyses, e. g., flux distribution and RMS-flux relation, are performed on the observations with an aim to compare the results with the similar ones in the gama-ray band presented in Bhatta & Dhital 2020. It is found that, similar to $gamma$-ray band, blazars display significant variability in the optical band that can be characterized with log-normal flux distribution and a power-law dependence of RMS on flux. It could be an indication of possible inherent linear RMS-flux relation, yet the scatter in the data does not allow to rule out other possibilities. When comparing variability properties in the two bands, the blazars in the gama-rays are found to exhibit stronger variability with steeper possible linear RMS-flux relation and the flux distribution that is more skewed towards higher fluxes. The cross-correlation study shows that except for the source 3C 273, the overall optical and the $gamma$-ray emission in the sources are highly correlated, suggesting a co-spatial existence of the particles responsible for both the optical and $gamma$-ray emission. Moreover, the sources S5 0716+714, Mrk 421, Mrk 501, PKS 1424-418 and PKS 2155-304 revealed possible evidence for quasi-periodic oscillations in the optical emission with the characteristic timescales, which are comparable to those in the $gamma$-ray band detected in our previous work.
We investigate the thermal emission and extinction from dust associated with the nearby superluminous supernova (SLSN) 2018bsz. Our dataset has daily cadence and simultaneous optical and near-infrared coverage up to ~ 100 days, together with late tim
e (+ 1.7 yr) MIR observations. At 230 days after light curve peak the SN is not detected in the optical, but shows a surprisingly strong near-infrared excess, with r - J > 3 mag and r - Ks > 5 mag. The time evolution of the infrared light curve enables us to investigate if the mid-infrared emission is from newly formed dust inside the SN ejecta, from a pre-existing circumstellar envelope, or interstellar material heated by the radiation from the SN. We find the latter two scenarios can be ruled out, and a scenario where new dust is forming in the SN ejecta at epochs > 200 days can self-consistently reproduce the evolution of the SN flux. We can fit the spectral energy distribution well at +230 d with 5 x 10^-4 solar mass of carbon dust, increasing over the following several hundred days to 10^-2 solar mass by +535 d. SN 2018bsz is the first SLSN showing evidence for dust formation within the SN ejecta, and appears to form ten times more dust than normal core-collapse SNe at similar epochs. Together with their preference for low mass, low metallicity host galaxies, we suggest that SLSNe may be a significant contributor to dust formation in the early Universe.
We perform a new dark matter hot spot analysis using ten years of public IceCube data. In this analysis we assume dark matter self-annihilates to neutrino pairs and treat the production sites as discrete point sources. For neutrino telescopes these s
ites will appear as hot spots in the sky, possibly outshining other standard model neutrino sources. Comparing to galactic center analyses, we show that this approach is a powerful tool and capable of setting the highest neutrino detector limits for dark matter masses between 10 TeV and 100 PeV. This is due to the inclusion of spatial information in addition to the typically used energy deposition in the analysis.
With the remarkable advent of gravitational-wave astronomy, we have shed light on previously shrouded events: compact binary coalescences. Neutron stars are promising (and confirmed) sources of gravitational radiation and it proves timely to consider
the ways in which these stars can be deformed. Gravitational waves provide a unique window through which to examine neutron-star interiors and learn more about the equation of state of ultra-dense nuclear matter. In this work, we study two relevant scenarios for gravitational-wave emission: neutron stars that host (non-axially symmetric) mountains and neutron stars deformed by the tidal field of a binary partner. Although they have yet to be seen with gravitational waves, rotating neutron stars have long been considered potential sources. By considering the observed spin distribution of accreting neutron stars with a phenomenological model for the spin evolution, we find evidence for gravitational radiation in these systems. We study how mountains are modelled in both Newtonian and relativistic gravity and introduce a new scheme to resolve issues with previous approaches to this problem. The crucial component of this scheme is the deforming force that gives the star its non-spherical shape. We find that the force (which is a proxy for the stars formation history), as well as the equation of state, plays a pivotal role in supporting the mountains. Considering a scenario that has been observed with gravitational waves, we calculate the structure of tidally deformed neutron stars, focusing on the impact of the crust. We find that the effect on the tidal deformability is negligible, but the crust will remain largely intact up until merger.
Lorentz Invariance Violation in Quantum Gravity (QG) models or a non-zero photon mass, $m_gamma$, would lead to an energy-dependent propagation speed for photons, such that photons of different energies from a distant source would arrive at different
times, even if they were emitted simultaneously. By developing source-by-source, Monte Carlo-based forward models for such time delays from Gamma Ray Bursts, and marginalising over empirical noise models describing other contributions to the time delay, we derive constraints on $m_gamma$ and the QG length scale, $ell_{rm QG}$, using spectral lag data from the BATSE satellite. We find $m_gamma < 4.0 times 10^{-5} , h , {rm eV}/c^2$ and $ell_{rm QG} < 5.3 times 10^{-18} , h , {rm , GeV^{-1}}$ at 95% confidence, and demonstrate that these constraints are robust to the choice of noise model. The QG constraint is among the tightest from studies which consider multiple Gamma Ray Bursts and the constraint on $m_gamma$, although weaker than from using radio data, provides an independent constraint which is less sensitive to the effects of dispersion by electrons.
Cosmic ray propagation is determined by the properties of interstellar turbulence. The multiphase nature of interstellar medium (ISM) and diversity of driving mechanisms give rise to spatial variation of turbulence properties. Meanwhile, precision as
troparticle experiments pose challenges to the conventional picture of homogeneous and isotropic transport of cosmic rays (CRs). We are beginning a new chapter of CR propagation research when studies of particle transport and interstellar turbulence confront each other. Here we review our recent developement on understandings of magnetohydrodynamic (MHD) turbulence and its connection to the fundamental processes governing cosmic ray propagation, different regimes of particle transport, that are augmented with observational discovery and analysis from multi-wavelength observations.
Star-forming regions have been proposed as potential Galactic cosmic-ray accelerators for decades. Cosmic ray acceleration can be probed through observations of gamma-rays produced in inelastic proton-proton collisions, at GeV and TeV energies. We an
alyze more than 11 years of Fermi-LAT data from the direction of Westerlund 2, one of the most massive and best-studied star-forming regions in our Galaxy. The spectral and morphological characteristics of the LAT source agree with the ones in the TeV regime (HESS J1023-575), allowing the description of the gamma-ray source from a few hundreds of MeV to a few tens of TeVs. We will present the results and discuss the implications of the identification with the stellar cluster and the radiation mechanism involved.
Ultrahigh-energy photons up to 1.4 peta-electronvolts have been observed by new cosmic-ray telescope in China -- a hint that Lorentz invariance might break down at the Planck-scale level.
GRB 060505 was the first well-known nearby (at redshift 0.089) hybrid gamma-ray burst that has a duration longer than 2 seconds but without the association of a supernova down to very stringent limits. The prompt $gamma-$ray flash lasting $sim 4$ sec
could consist of an intrinsic short burst and its tail emission, but the sizable temporal lag ($sim 0.35$ sec) as well as the environment properties led to the widely-accepted classification of a long duration gamma-ray burst originated from the collapse of a massive star. Here for the $ first$ time we report the convincing evidence for a thermal-like optical radiation component in the spectral energy distribution of the early afterglow emission. In comparison to AT2017gfo, the thermal radiation is $sim 2$ times brighter and the temperature is comparable at similar epochs. The optical decline is much quicker than that in X-rays, which is also at odds with the fireball afterglow model but quite natural for the presence of a blue kilonova. Our finding reveals a neutron star merger origin of the hybrid GRB 060505 and strongly supports the theoretical speculation that some binary neutron stars can merge ultra-quickly (within $sim 1$ Myr) after their formation when the surrounding region is still highly star-forming and the metallicity remains low. Gravitational wave and electromagnetic jointed observations are expected to confirm such scenarios in the near future.
We extend previous work on gamma-ray burst (GRB) afterglows involving hot thermal electrons at the base of a shock-accelerated tail. Using a physically-motivated electron distribution based on first-principles simulations, we compute broadband emissi
on from radio to TeV gamma-rays. For the first time, we present the effects of a thermal distribution of electrons on synchrotron self-Compton (SSC) emission. The presence of thermal electrons causes temporal and spectral structure across the entire observable afterglow, which is substantively different from models that assume a pure power-law distribution for the electrons. We show that early-time TeV emission is enhanced by more than an order of magnitude for our fiducial parameters, with a time-varying spectral index that does not occur for a pure power law of electrons. We further show that the X-ray closure relations take a very different, also time-dependent, form when thermal electrons are present; the shape traced out by the X-ray afterglows is a qualitative match to observations of the traditional decay phase.
