Search

X-ray polarization is a unique new probe of the particle acceleration in astrophysical jets made possible through the Imaging X-ray Polarimetry Explorer. Here we report on the first dense X-ray polarization monitoring campaign on the blazar Mrk 421. Our observations were accompanied by an even denser radio and optical polarization campaign. We find significant short-timescale variability in both X-ray polarization degree and angle, including a $\sim90^\circ$ angle rotation about the jet axis. We attribute this to random variations of the magnetic field, consistent with the presence of turbulence but also unlikely to be explained by turbulence alone. At the same time, the degree of lower-energy polarization is significantly lower and shows no more than mild variability. Our campaign provides further evidence for a scenario in which energy-stratified shock-acceleration of relativistic electrons, combined with a turbulent magnetic field, is responsible for optical to X-ray synchrotron emission in blazar jets., Comment: 23 pages, including 8 pages of appendices. 12 figures, 3 tables. Submitted to ApJ

We present multiwavelength polarization measurements of the luminous blazar Mrk~501 over a 14-month period. The 2--8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100-ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optical-infrared polarization measurements were also available in the B, V, R, I, and J bands, as were radio polarization measurements from 4.85 GHz to 225.5 GHz. Among the first five IXPE observations, we did not find significant variability in the X-ray polarization degree and angle with IXPE. However, the most recent sixth observation found an elevated polarization degree at $>3\sigma$ above the average of the other five observations. The optical and radio measurements show no apparent correlations with the X-ray polarization properties. Throughout the six IXPE observations, the X-ray polarization degree remained higher than, or similar to, the R-band optical polarization degree, which remained higher than the radio value. This is consistent with the energy-stratified shock scenario proposed to explain the first two IXPE observations, in which the polarized X-ray, optical, and radio emission arises from different regions., Comment: Accepted for publication in ApJ

Variability can be the pathway to understanding the physical processes in astrophysical jets, however, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, >24 hour polarization light curves of blazars using telescopes distributed across the globe and the rotation of the Earth to avoid the rising Sun. Our campaign involved 16 telescopes in Asia, Europe, and North America. We observed BL Lacertae and CGRaBS J0211+1051 for a combined 685 telescope hours. We find large variations in the polarization degree and angle for both sources in sub-hour timescales as well as a ~180 degree rotation of the polarization angle in CGRaBS J0211+1051 in less than two days. We compared our high-cadence observations to Particle-In-Cell magnetic reconnection and turbulent plasma simulations. We find that although the state of the art simulation frameworks can produce a large fraction of the polarization properties, they do not account for the entirety of the observed polarization behavior in blazar jets., Comment: 20 pages, 15 figures, 2 tables, accepted for publication in A&A. The data used in the paper are available here: https://doi.org/10.7910/DVN/IETSXS

We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the highest X-ray polarization degree detected for an HSP blazar thus far, (30.7$\pm$2.0)%, which dropped to (15.3$\pm$2.1)% during the second half ($T_2$). The X-ray polarization angle remained stable during the IXPE pointing at 129.4$^\circ$$\pm$1.8$^\circ$ and 125.4$^\circ$$\pm$3.9$^\circ$ during $T_1$ and $T_2$, respectively. Meanwhile, the optical polarization degree remained stable during the IXPE pointing, with average host-galaxy-corrected values of (4.3$\pm$0.7)% and (3.8$\pm$0.9)% during the $T_1$ and $T_2$, respectively. During the IXPE pointing, the optical polarization angle changed achromatically from $\sim$140$^\circ$ to $\sim$90$^\circ$ and back to $\sim$130$^\circ$. Despite several attempts, we only detected (99.7% conf.) the radio polarization once (during $T_2$, at 225.5 GHz): with degree (1.7$\pm$0.4)% and angle 112.5$^\circ$$\pm$5.5$^\circ$. The direction of the broad pc-scale jet is rather ambiguous and has been found to point to the east and south at different epochs; however, on larger scales (> 1.5 pc) the jet points toward the southeast ($\sim$135$^\circ$), similar to all of the MW polarization angles. Moreover, the X-ray to optical polarization degree ratios of $\sim$7 and $\sim$4 during $T_1$ and $T_2$, respectively, are similar to previous IXPE results for several HSP blazars. These findings, combined with the lack of correlation of temporal variability between the MW polarization properties, agree with an energy-stratified shock-acceleration scenario in HSP blazars., Comment: 17 pages, 10 figures, 4 tables, Accepted for publication in A&A

Observations of linear polarization in the 2-8 keV energy range with the Imaging X-ray Polarimetry Explorer (IXPE) explore the magnetic field geometry and dynamics of the regions generating non-thermal radiation in relativistic jets of blazars. These jets, particularly in blazars whose spectral energy distribution peaks at X-ray energies, emit X-rays via synchrotron radiation from high-energy particles within the jet. IXPE observations of the X-ray selected BL Lac-type blazar 1ES 1959+650 in 2022 May 3-4 showed a significant linear polarization degree of $\Pi_\mathrm{x} = 8.0\% \pm 2.3\%$ at an electric-vector position angle $\psi_\mathrm{x} = 123^\circ \pm 8^\circ$. However, in 2022 June 9-12, only an upper limit of $\Pi_\mathrm{x} \leq 5.1\%$ could be derived (at the 99% confidence level). The degree of optical polarization at that time $\Pi_\mathrm{O} \sim 5\%$ is comparable to the X-ray measurement. We investigate possible scenarios for these findings, including temporal and geometrical depolarization effects. Unlike some other X-ray selected BL Lac objects, there is no significant chromatic dependence of the measured polarization in 1ES 1959+650, and its low X-ray polarization may be attributed to turbulence in the jet flow with dynamical timescales shorter than 1 day., Comment: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal

We present X-ray polarimetry observations from the Imaging X-ray Polarimetry Explorer (IXPE) of three low spectral peak and one intermediate spectral peak blazars, namely 3C 273, 3C 279, 3C 454.3, and S5 0716+714. For none of these objects was IXPE able to detect X-ray polarization at the 3$\sigma$ level. However, we placed upper limits on the polarization degree at $\sim$10-30\%. The undetected polarizations favor models where the X-ray band is dominated by unpolarized photons upscattered by relativistic electrons in the jets of blazars, although hadronic models are not completely eliminated. We discuss the X-ray polarization upper limits in the context of our contemporaneous multiwavelength polarization campaigns., Comment: 19 pages, 10 figures, 6 tables; submitted to the Astrophysical Journal

We conducted a polarimetry campaign from radio to X-ray wavelengths of the high-synchrotron-peak (HSP) blazar Mrk 421, including Imaging X-ray Polarimetry Explorer (IXPE) measurements on 2022 December 6-8. We detected X-ray polarization of Mrk 421 with a degree of $\Pi_{\rm X}$=14$\pm$1$\%$ and an electric-vector position angle $\psi_{\rm X}$=107$\pm$3$^{\circ}$ in the 2-8 keV band. From the time variability analysis, we find a significant episodic variation in $\psi_{\rm X}$. During 7 months from the first IXPE pointing of Mrk 421 in 2022 May, $\psi_{\rm X}$ varied across the range of 0$^{\circ}$ to 180$^{\circ}$, while $\Pi_{\rm X}$ maintained similar values within $\sim$10-15$\%$. Furthermore, a swing in $\psi_{\rm X}$ in 2022 June was accompanied by simultaneous spectral variations. The results of the multiwavelength polarimetry show that the X-ray polarization degree was generally $\sim$2-3 times greater than that at longer wavelengths, while the polarization angle fluctuated. Additionally, based on radio, infrared, and optical polarimetry, we find that rotation of $\psi$ occurred in the opposite direction with respect to the rotation of $\psi_{\rm X}$ over longer timescales at similar epochs. The polarization behavior observed across multiple wavelengths is consistent with previous IXPE findings for HSP blazars. This result favors the energy-stratified shock model developed to explain variable emission in relativistic jets. The accompanying spectral variation during the $\psi_{\rm X}$ rotation can be explained by a fluctuation in the physical conditions, e.g., in the energy distribution of relativistic electrons. The opposite rotation direction of $\psi$ between the X-ray and longer-wavelength polarization accentuates the conclusion that the X-ray emitting region is spatially separated from that at longer wavelengths., Comment: 17 pages, 13 figures, 4 tables; Accepted for publication in A&A

The highly variable blazar OJ~287 is commonly discussed as an example of a binary black hole system. The 130 year long optical light curve is well explained by a model where the central body is a massive black hole of 18.35$\times$10$^9$ solar mass that supports a thin accretion disc. The secondary black hole of 0.15$\times$10$^9$ solar mass impacts the disc twice during its 12 year orbit, and causes observable flares. Recently, it has been argued that an accretion disc with a typical AGN accretion rate and above mentioned central body mass should be at least six magnitudes brighter than OJ~287's host galaxy and would therefore be observationally excluded. Based on the observations of OJ~287's radio jet, detailed in Marscher and Jorstad (2011), and up-to-date accretion disc models of Azadi et al. (2022), we show that the V-band magnitude of the accretion disc is unlikely to exceed the host galaxy brightness by more than one magnitude, and could well be fainter than the host. This is because accretion power is necessary to launch the jet as well as to create electromagnetic radiation, distributed across many wavelengths, and not concentrated especially on the optical V-band. Further, we note that the claimed V-band concentration of accretion power leads to serious problems while interpreting observations of other Active Galactic Nuclei. Therefore, we infer that the mass of the primary black hole and its accretion rate do not need to be smaller than what is determined in the standard model for OJ~287.

We present a summary of the results of the OJ 287 observational campaign, which was carried out during the 2021/2022 observational season. This season is special in the binary model because the major axis of the precessing binary happens to lie almost exactly in the plane of the accretion disc of the primary. This leads to pairs of almost identical impacts between the secondary black hole and the accretion disk in 2005 and 2022. In 2005, a special flare called "blue flash" was observed 35 days after the disk impact, which should have also been verifiable in 2022. We did observe a similar flash and were able to obtain more details of its properties. We describe this in the framework of expanding cloud models. In addition, we were able to identify the flare arising exactly at the time of the disc crossing from its photo-polarimetric and gamma-ray properties. This is an important identification, as it directly confirms the orbit model. Moreover, we saw a huge flare that lasted only one day. We may understand this as the lighting up of the jet of the secondary black hole when its Roche lobe is suddenly flooded by the gas from the primary disk. Therefore, this may be the first time we directly observed the secondary black hole in the OJ 287 binary system., Comment: 24 pages, 13 figures, published in Galaxies as part of the Special Issue Distant Glowing Objects: Quest for Quasars

We report the first $> 99\%$ confidence detection of X-ray polarization in BL Lacertae. During a recent X-ray/$\gamma$-ray outburst, a 287 ksec observation (2022 November 27-30) was taken using the Imaging X-ray Polarimetry Explorer ({\it IXPE}), together with contemporaneous multiwavelength observations from the Neil Gehrels {\it Swift} observatory and {\it XMM-Newton} in soft X-rays (0.3--10~keV), {\it NuSTAR} in hard X-rays (3--70~keV), and optical polarization from the Calar Alto, and Perkins Telescope observatories. Our contemporaneous X-ray data suggest that the {\it IXPE} energy band is at the crossover between the low- and high-frequency blazar emission humps. The source displays significant variability during the observation, and we measure polarization in three separate time bins. Contemporaneous X-ray spectra allow us to determine the relative contribution from each emission hump. We find $>99\%$ confidence X-ray polarization $\Pi_{2-4{\rm keV}} = 21.7^{+5.6}_{-7.9}\%$ and electric vector polarization angle $\psi_{2-4{\rm keV}} = -28.7 \pm 8.7^{\circ}$ in the time bin with highest estimated synchrotron flux contribution. We discuss possible implications of our observations, including previous {\it IXPE} BL Lacertae pointings, tentatively concluding that synchrotron self-Compton emission dominates over hadronic emission processes during the observed epochs., Comment: Accepted to ApJL. 16 pages, 8 figures

The magnetic field conditions in astrophysical relativistic jets can be probed by multiwavelength polarimetry, which has been recently extended to X-rays. For example, one can track how the magnetic field changes in the flow of the radiating particles by observing rotations of the electric vector position angle $\Psi$. Here we report the discovery of a $\Psi_{\mathrm x}$ rotation in the X-ray band in the blazar Mrk 421 at an average flux state. Across the 5 days of Imaging X-ray Polarimetry Explorer (IXPE) observations of 4-6 and 7-9 June 2022, $\Psi_{\mathrm x}$ rotated in total by $\geq360^\circ$. Over the two respective date ranges, we find constant, within uncertainties, rotation rates ($80 \pm 9$ and $91 \pm 8 ^\circ/\rm day$) and polarization degrees ($\Pi_{\mathrm x}=10\%\pm1\%$). Simulations of a random walk of the polarization vector indicate that it is unlikely that such rotation(s) are produced by a stochastic process. The X-ray emitting site does not completely overlap the radio/infrared/optical emission sites, as no similar rotation of $\Psi$ was observed in quasi-simultaneous data at longer wavelengths. We propose that the observed rotation was caused by a helical magnetic structure in the jet, illuminated in the X-rays by a localized shock propagating along this helix. The optically emitting region likely lies in a sheath surrounding an inner spine where the X-ray radiation is released.

The bright blazar OJ~287 routinely parades high brightness bremsstrahlung flares, which are explained as being a result of a secondary supermassive black hole (SMBH) impacting the accretion disc of a more massive primary SMBH in a binary system. The accretion disc is not rigid but rather bends in a calculable way due to the tidal influence of the secondary. Below we refer to this phenomenon as a variable disc level. We begin by showing that these flares occur at times predicted by a simple analytical formula, based on general relativity inspired modified Kepler equation, which explains impact flares since 1888. The 2022 impact flare, namely flare number 26, is rather peculiar as it breaks the typical pattern of two impact flares per 12-year cycle. This is the third bremsstrahlung flare of the current cycle that follows the already observed 2015 and 2019 impact flares from OJ~287. It turns out that the arrival epoch of flare number 26 is sensitive to the level of primary SMBH's accretion disc relative to its mean level in our model. We incorporate these tidally induced changes in the level of the accretion disc to infer that the thermal flare should have occurred during July-August 2022, when it was not possible to observe it from the Earth. Thereafter, we explore possible observational evidence for certain pre-flare activity by employing spectral and polarimetric data from our campaigns in 2004/05 and 2021/22. We point out theoretical and observational implications of two observed mini-flares during January-February 2022., Comment: 29 pages, 6 figures, 1 table. arXiv admin note: text overlap with arXiv:2209.08360

Blazars are a class of jet-dominated active galactic nuclei with a typical double-humped spectral energy distribution. It is of common consensus the Synchrotron emission to be responsible for the low frequency peak, while the origin of the high frequency hump is still debated. The analysis of X-rays and their polarization can provide a valuable tool to understand the physical mechanisms responsible for the origin of high-energy emission of blazars. We report the first observations of BL Lacertae performed with the Imaging X-ray Polarimetry Explorer ({IXPE}), from which an upper limit to the polarization degree $\Pi_X<$12.6\% was found in the 2-8 keV band. We contemporaneously measured the polarization in radio, infrared, and optical wavelengths. Our multiwavelength polarization analysis disfavors a significant contribution of proton synchrotron radiation to the X-ray emission at these epochs. Instead, it supports a leptonic origin for the X-ray emission in BL Lac., Comment: 17 pages, 5 figures, accepted for publication in ApJL

In the first days of WZ Sge-type dwarf nova (DN) outbursts, the 2:1 resonance induces a spiral arm structure in the accretion disk, which is observed as early superhumps in optical light curves. This paper reports our optical observations of an eclipsing WZ Sge-type DN PNV J00444033+4113068 during its 2021 superoutburst with the 3.8m Seimei telescope and through VSNET collaboration. The eclipse analysis gave its orbital period as 0.055425534(1) d. Our observations confirmed early superhumps with an amplitude of 0.7 mag, the largest amplitude among known WZ Sge-type DNe. More interestingly, its early superhumps became the reddest around their secondary minimum, whereas other WZ Sge-type DNe show the reddest color around the early superhump maximum. The spectrum around the peak of the outburst showed the double-peaked emission lines of He II 4686\AA~ and H$\alpha$ with a peak separation of $\ge 700$ km/s, supporting a very high-inclination system. With the early superhump mapping, the unique profile and color of the early superhump of PNV J00444033+4113068 are successfully reproduced by the accretion disk with vertically extended double arm structure. Therefore, the large amplitude and unique color behavior of the early superhumps in PNV J00444033+4113068 can be explained by the 2:1 resonance model along with other WZ Sge-type DNe., Comment: 8 pages, 4 figures, 1 table, accepted for publication in PASJ

Simultaneous X-ray and optical observations of black hole X-ray binaries have shown that the light curves contain multiple correlated and anti-correlated variation components when the objects are in the hard state. In the case of the black hole X-ray binary, GX 339-4, the cross correlation function (CCF) of the light curves suggests a positive correlation with an optical lag of 0.15 s and anti-correlations with an optical lag of 1 s and X-ray lag of 4 s. This indicates the two light curves have some common signal components with different delays. In this study, we extracted and reconstructed those signal components from the data for GX 339-4. The results confirmed that correlation and anti-correlation with the optical lag are two common components. However, we found that the reconstructed light curve for the anti-correlated component indicates a positively correlated variation with an X-ray lag of ~ +1 s. In addition, the CCF for this signal component shows anti-correlations not only with the optical lag, but also with the X-ray lag, which is consistent with the CCF for the data. Therefore, our results suggest that the combination of the two positively correlated components, that is, the X-ray preceding signal with the 0.15-s optical lag and the optical preceding signal with the 1-s X-ray lag, can make the observed CCF without anti-correlated signals. The optical preceding signal may be caused by synchrotron emission in a magnetically dominated accretion flow or in a jet, while further study is required to understand the mechanism of the X-ray time lag.

Blazars can be divided into two subtypes, flat spectrum radio quasars (FSRQs) and BL Lac objects, which have been distinguished phenomenologically by the strength of their optical emission lines, while their physical nature and relationship are still not fully understood. In this paper, we focus on the differences in their variability. We characterize the blazar variability using the Ornstein-Uhlenbeck (OU) process, and investigate the features that are discriminative for the two subtypes. We used optical photometric and polarimetric data obtained with the 1.5-m Kanata telescope for 2008-2014. We found that four features, namely the variation amplitude, characteristic timescale, and non-stationarity of the variability obtained from the light curves and the median of the degree of polarization (PD), are essential for distinguishing between FSRQs and BL Lac objects. FSRQs are characterized by rare and large flares, while the variability of BL Lac objects can be reproduced with a stationary OU process with relatively small amplitudes. The characteristics of the variability are governed not by the differences in the jet structure between the subtypes, but by the peak frequency of the synchrotron emission. This implies that the nature of the variation in the jets is common in FSRQs and BL Lac objects. We found that BL Lac objects tend to have high PD medians, which suggests that they have a stable polarization component. FSRQs have no such component, possibly because of a strong Compton cooling effect in sub-pc scale jets., Comment: 13 pages, 10 figures, accepted for publication in PASJ

We present optical and near-infrared observations of SN~Ib~2019ehk. We show that it evolved to a Ca-rich transient according to its spectral properties and evolution in late phases. It, however, shows a few distinguishable properties from the canonical Ca-rich transients: a short-duration first peak in the light curve, high peak luminosity, and association with a star-forming environment. Indeed, some of these features are shared with iPTF14gqr and iPTF16hgs, which are candidates for a special class of core-collapse SNe (CCSNe): the so-called ultra-stripped envelope SNe, i.e., a relatively low-mass He (or C+O) star explosion in a binary as a precursor of double neutron star binaries. The estimated ejecta mass ($0.43 M_\odot$) and explosion energy ($1.7 \times 10^{50} $~erg) are consistent with this scenario. The analysis of the first peak suggests existence of dense circumstellar material in the vicinity of the progenitor, implying a CCSN origin. Based on these analyses, we suggest SN 2019ehk is another candidate for an ultra-stripped envelope SN. These ultra-stripped envelope SN candidates seem to form a subpopulation among Ca-rich transients, associated with young population. We propose that the key to distinguishing this population is the early first peak in their light curves., Comment: submitted to ApJ on Mar 27, 2020

We report on the variations of the physical parameters of the jet observed in the blazar Mrk 421, and discuss the origin of X-ray flares in the jet, based on the analysis of the several spectral energy distributions (SEDs). The SEDs are modeled using the one-zone synchrotron self-Compton (SSC) model and its parameters determined using a Markov chain Monte Carlo method. The lack of data at TeV energies means many of the parameters cannot be uniquely determined and are correlated. These are studied in detail. We found that the optimal solution can be uniquely determined only when we apply a constraint to one of four parameters: the magnetic field (B), Doppler factor, size of the emitting region, and normalization factor of the electron energy distribution. We used 31 sets of SED from 2009 to 2014 with optical-UV data observed with UVOT/Swift and the Kanata telescope, X-ray data with XRT/Swift, and gamma-ray data with the Fermi Large Area Telescope (LAT). The result of our SED analysis suggests that, in the X-ray faint state, the emission occurs in a relatively small area (~ 10^16 cm) with relatively strong magnetic field (B~10^-1 G). The X-ray bright state shows a tendency opposite to that of the faint state, that is, a large emitting area (~10^18 cm), probably in the downstream of the jet and weak magnetic field (B~10^-3 G). The high X-ray flux was due to an increase in the maximum energy of electrons. On the other hand, the presence of two kinds of emitting areas implies that the one-zone model is unsuitable to reproduce, at least a part of the observed SEDs., Comment: 20 pages, 16 figures, accepted for publication in PASJ

We present optical and near-infrared observations of the rapidly evolving supernova (SN) 2017czd that shows hydrogen features. The optical light curves exhibit a short plateau phase ($\sim 13$ days in the $R$-band) followed by a rapid decline by $4.5$ mag in $\sim 20 \mathrm{days}$ after the plateau. The decline rate is larger than those of any standard SNe, and close to those of rapidly evolving transients. The peak absolute magnitude is $-16.8$ mag in the $V$-band, which is within the observed range for SNe IIP and rapidly evolving transients. The spectra of SN 2017czd clearly show the hydrogen features and resemble those of SNe IIP at first. The H$\alpha$ line, however, does not evolve much with time and it becomes similar to those in SNe IIb at decline phase. We calculate the synthetic light curves using a SN IIb progenitor which has 16 M$_{\odot}$ at the zero-age main sequence and evolves in a binary system. The model with a low explosion energy ($5\times 10^{50}$ erg) and a low ${}^{56}$Ni mass ($0.003 \mathrm{M}_{\odot}$) can reproduce the short plateau phase as well as the sudden drop of the light curve as observed in SN 2017czd. We conclude that SN 2017czd might be the first identified weak explosion from a SN IIb progenitor. We suggest that some rapidly evolving transients can be explained by such a weak explosion of the progenitors with little hydrogen-rich envelope., Comment: 11 pages, 15 figures. Accepted for publication in ApJ

We present optical and near-infrared observations of the nearby Type Iax supernova (SN) 2014dt from 14 to 410 days after the maximum light. The velocities of the iron absorption lines in the early phase indicated that SN 2014dt showed slower expansion than the well-observed Type Iax SNe 2002cx, 2005hk and 2012Z. In the late phase, the evolution of the light curve and that of the spectra were considerably slower. The spectral energy distribution kept roughly the same shape after ~100 days, and the bolometric light curve flattened during the same period. These observations suggest the existence of an optically thick component that almost fully trapped the {\gamma}-ray energy from 56 Co decay. These findings are consistent with the predictions of the weak deflagration model, leaving a bound white dwarf remnant after the explosion., Comment: 39 pages, 14 figures, accepted for publication in PASJ

We present optical and near-infrared observations of a low-luminosity Type IIP supernova (SN) 2016bkv from the initial rising phase to the plateau phase. Our observations show that the end of the plateau is extended to $\gtrsim 140$ days since the explosion, indicating that this SN takes one of the longest time to finish the plateau phase. among Type IIP SNe (SNe IIP), including low-luminosity (LL) SNe IIP. The line velocities of various ions at the middle of the plateau phase are as low as 1,000--1,500 km s$^{-1}$, which is the lowest even among LL SNe IIP. These measurements imply that the ejecta mass in SN 2016bkv is larger than that of the well-studied LL IIP SN 2003Z. In the early phase, SN 2016bkv shows a strong bump in the light curve. In addition, the optical spectra in this bump phase exhibit a blue continuum accompanied with a narrow H$\alpha$ emission line. These features indicate an interaction between the SN ejecta and the circumstellar matter (CSM) as in SNe IIn. Assuming the ejecta-CSM interaction scenario, the mass loss rate is estimated to be $\sim 1.7 \times 10^{-2} M_{\odot}$ yr$^{-1}$ within a few years before the SN explosion. This is comparable to or even larger than the largest mass loss rate observed for the Galactic red supergiants ($\sim 10^{-3} M_{\odot}$ yr$^{-1}$ for VY CMa). We suggest that the progenitor star of SN 2016bkv experienced a violent mass loss just before the SN explosion., Comment: 9 pages, 12 figures. Accepted for publication in ApJ

The first detected gravitational wave from a neutron star merger was GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the $z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder for later epochs; the $z-H$ color changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The rapid evolution of both the optical brightness and the color are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that $r$-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an additional energy source., Comment: 7 pages, 5 figures. Published in PASJ (Publications of the Astronomical Society of Japan)

We present the results of the detailed analysis of an optical imaging survey conducted using the Subaru / Hyper Suprime-Cam (HSC), which aims to identify an optical counterpart to the gravitational wave event GW151226. In half a night, the $i$- and $z$-band imaging survey by HSC covers 63.5deg$^2$ of the error region, which contains about 7\% of the LIGO localization probability, and the same field is observed in three different epochs. The detectable magnitude of the candidates in a differenced image is evaluated as $i \sim 23.2$ mag for the requirement of at least two 5$\sigma$ detections, and 1744 candidates are discovered. Assuming a kilonova as an optical counterpart, we compared the optical properties of the candidates with model predictions. A red and rapidly declining light curve condition enables the discrimination of a kilonova from other transients, and a small number of candidates satisfy this condition. The presence of stellar-like counterparts in the reference frame suggests that the surviving candidates are likely to be flare stars. The fact that most of those candidates are in galactic plane, $|b|<5^{\circ}$, supports this interpretation. We also checked whether the candidates are associated with the nearby GLADE galaxies, which reduces the number of contaminants even with a looser color cut. When a better probability map (with localization accuracy of $\sim50{\rm deg}^2$) is available, kilonova searches of up to approximately $200$ Mpc will become feasible by conducting immediate follow-up observations with an interval of 3--6 days., Comment: 13 pages, 10 figures, Submitted to PASJ on 12 Aug 2017

We report on the variation in the optical polarization of the blazar PKS 1749+096 observed in 2008--2015. The degree of polarization (PD) tends to increase in short flares having a time-scale of a few days. The object favors a polarization angle (PA) of $40^\circ$--$50^\circ$ at the flare maxima, which is close to the position angle of the jet ($20^\circ$--$40^\circ$). Three clear polarization rotations were detected in the negative PA direction associated with flares. In addition, a rapid and large decrease in the PA was observed in the other two flares, while another two flares showed no large PA variation. The light curve maxima of the flares possibly tend to lag behind the PD maxima and color-index minima. The PA became $-50^\circ$ to $-20^\circ$ in the decay phase of active states, which is almost perpendicular to the jet position angle. We propose a scenario to explain these observational features, where transverse shocks propagate along curved trajectories. The favored PA at the flare maxima suggests that the observed variations were governed by the variations in the Doppler factor, $\delta$. Based on this scenario, the minimum viewing angle of the source, $\theta_\mathrm{min}=4.8^\circ$--$6.6^\circ$, and the location of the source, $\Delta r\gtrsim 0.1$pc, from the central black hole were estimated. In addition, the acceleration of electrons by the shock and synchrotron cooling would have a time-scale similar to that of the change in $\delta$. The combined effect of the variation in $\delta$ and acceleration/cooling of electrons is probably responsible for the observed diversity of the polarization variations in the flares., Comment: 12 pages, 9 figures, accepted for publication in PASJ

We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~\mu$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $\alpha_{\nu}\sim-1.4$, where $f_{\nu}\propto\nu^{\alpha_{\nu}}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $\nu_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations., Comment: 22 pages, 7 figures, 2 tables. Accepted for publication in PASJ

We propose a K-sparse exhaustive search (ES-K) method and a K-sparse approximate exhaustive search method (AES-K) for selecting variables in linear regression. With these methods, K-sparse combinations of variables are tested exhaustively assuming that the optimal combination of explanatory variables is K-sparse. By collecting the results of exhaustively computing ES-K, various approximate methods for selecting sparse variables can be summarized as density of states. With this density of states, we can compare different methods for selecting sparse variables such as relaxation and sampling. For large problems where the combinatorial explosion of explanatory variables is crucial, the AES-K method enables density of states to be effectively reconstructed by using the replica-exchange Monte Carlo method and the multiple histogram method. Applying the ES-K and AES-K methods to type Ia supernova data, we confirmed the conventional understanding in astronomy when an appropriate K is given beforehand. However, we found the difficulty to determine K from the data. Using virtual measurement and analysis, we argue that this is caused by data shortage., Comment: 19pages, 3 figures

We present the early-phase spectra and the light curves of the broad-lined supernova (SN) 2016coi from $t=7$ to $67$ days after the estimated explosion date. This SN was initially reported as a broad-lined Type SN Ic (SN Ic-BL). However, we found that spectra up to $t=12$ days exhibited the He~{\sc i} $\lambda$5876, $\lambda$6678, and $\lambda$7065 absorption lines. We show that the smoothed and blueshifted spectra of normal SNe Ib are remarkably similar to the observed spectrum of SN 2016coi. The line velocities of SN 2016coi were similar to those of SNe Ic-BL and significantly faster than those of SNe Ib. Analyses of the line velocity and light curve suggest that the kinetic energy and the total ejecta mass of SN 2016coi are similar to those of SNe Ic-BL. Together with broad-lined SNe 2009bb and 2012ap for which the detection of He~{\sc i} were also reported, these SNe could be transitional objects between SNe Ic-BL and SNe Ib, and be classified as broad-lined Type `Ib' SNe (SNe `Ib'-BL). Our work demonstrates the diversity of the outermost layer in broad-lined SNe, which should be related to the variety of the evolutionary paths., Comment: 9 pages, 10 figures. Accepted for publication in ApJ

V404 Cygni is a well-known black hole binary candidate thought to have relativistic jets. It showed extreme outbursts on June 2015, characterized by a large amplitude and short time variation of flux in the radio, optical, and X-ray bands. Not only disk emission, but also synchrotron radiation from the relativistic jets were suggested by radio observations. However, it is difficult to measure the accurate spectral shape in the optical/near infrared band because there are uncertainties of interstellar extinction. To estimate the extinction value for V404 Cygni, we performed photopolarimetric and spectroscopic observations of V404 Cygni and nearby field stars. Here, we estimate the Galactic extinction using interstellar polarization based on the observation that the origin of the optical polarization is the interstellar medium, and investigate the properties of interstellar polarization around V404~Cygni. We found a good correlation between the color excess and polarization degree in the field stars. We also confirmed that the wavelength dependence of the polarization degree in the highly polarized field stars was similar to that of V404~Cygni. Using the highly polarized field stars, we estimated the color excess and the extinction, $E(B-V)=1.2 \pm 0.2$ and $3.0 < A(V) < 3.6$, respectively. A tendency for a bluer peak of polarization ($\lambda_{\rm max}<5500$ \AA) was commonly seen in the highly polarized field stars, suggesting that the dust grains toward this region are generally smaller than the Galactic average. The corrected spectral energy distribution of V404~Cygni in the near infrared (NIR) and optical bands in our results indicated a spectral break between $2.5 \times 10^{14}$ Hz and $3.7 \times 10^{14}$ Hz, which might be originated in the synchrotron self absorption., Comment: 11 pages, 6 figures, 3 tables. Accepted for publication in PASJ

We report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226., Comment: 12 pages, 6 figures, 4 tables, accepted for publication in PASJ

Cross-validation (CV) is a technique for evaluating the ability of statistical models/learning systems based on a given data set. Despite its wide applicability, the rather heavy computational cost can prevent its use as the system size grows. To resolve this difficulty in the case of Bayesian linear regression, we develop a formula for evaluating the leave-one-out CV error approximately without actually performing CV. The usefulness of the developed formula is tested by statistical mechanical analysis for a synthetic model. This is confirmed by application to a real-world supernova data set as well., Comment: 5 pages, 2 figures, invited paper for Allerton2016 conference

Blazars are highly variable active galactic nuclei which emit radiation at all wavelengths from radio to gamma-rays. Polarized radiation from blazars is one key piece of evidence for synchrotron radiation at low energies and it also varies dramatically. The polarization of blazars is of interest for understanding the origin, confinement, and propagation of jets. However, even though numerous measurements have been performed, the mechanisms behind jet creation, composition and variability are still debated. We performed simultaneous gamma-ray and optical photopolarimetry observations of 45 blazars between Jul. 2008 and Dec. 2014 to investigate the mechanisms of variability and search for a basic relation between the several subclasses of blazars. We identify a correlation between the maximum degree of optical linear polarization and the gamma-ray luminosity or the ratio of gamma-ray to optical fluxes. Since the maximum polarization degree depends on the condition of the magnetic field (chaotic or ordered), this result implies a systematic difference in the intrinsic alignment of magnetic fields in pc-scale relativistic jets between different types blazars (FSRQs vs. BL Lacs), and consequently between different types of radio galaxies (FR Is vs. FR IIs)., Comment: 17 pages, 30 figures, 6 tables. Accepted for publication in Astrophysical Journal

How black holes accrete surrounding matter is a fundamental, yet unsolved question in astrophysics. It is generally believed that matter is absorbed into black holes via accretion disks, the state of which depends primarily on the mass-accretion rate. When this rate approaches the critical rate (the Eddington limit), thermal instability is supposed to occur in the inner disc, causing repetitive patterns of large-amplitude X-ray variability (oscillations) on timescales of minutes to hours. In fact, such oscillations have been observed only in sources with a high mass accretion rate, such as GRS 1915+105. These large-amplitude, relatively slow timescale, phenomena are thought to have physical origins distinct from X-ray or optical variations with small amplitudes and fast ($\lesssim$10 sec) timescales often observed in other black hole binaries (e.g., XTE J1118+480 and GX 339-4). Here we report an extensive multi-colour optical photometric data set of V404 Cygni, an X-ray transient source containing a black hole of nine solar masses (and a conpanion star) at a distance of 2.4 kiloparsecs. Our data show that optical oscillations on timescales of 100 seconds to 2.5 hours can occur at mass-accretion rates more than ten times lower than previously thought. This suggests that the accretion rate is not the critical parameter for inducing inner-disc instabilities. Instead, we propose that a long orbital period is a key condition for these large-amplitude oscillations, because the outer part of the large disc in binaries with long orbital periods will have surface densities too low to maintain sustained mass accretion to the inner part of the disc. The lack of sustained accretion -- not the actual rate -- would then be the critical factor causing large-amplitude oscillations in long-period systems., Comment: Published in Nature on January 7th, 2016

From 2013 April to 2014 April, we performed an X-ray and optical simultaneous monitoring of the type 1.5 Seyfert galaxy NGC 3516. It employed Suzaku, and 5 Japanese ground-based telescopes, the Pirka, Kiso Schmidt, Nayuta, MITSuME, and the Kanata telescopes. The Suzaku observations were conducted seven times with various intervals ranging from days, weeks, to months, with an exposure of $\sim50$ ksec each. The optical $B$-band observations not only covered those of Suzaku almost simultaneously, but also followed the source as frequently as possible. As a result, NGC 3516 was found in its faint phase with the 2-10 keV flux of $0.21-2.70 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$. The 2-45 keV X-ray spectra were composed of a dominant variable hard power-law continuum with a photon index of $\sim1.7$, and a non-relativistic reflection component with a prominent Fe-K$\alpha$ emission line. Producing the $B$-band light curve by differential image photometry, we found that the $B$-band flux changed by $\sim2.7 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$, which is comparable to the X-ray variation, and detected a significant flux correlation between the hard power-law component in X-rays and the $B$-band radiation, for the first time in NGC 3516. By examining their correlation, we found that the X-ray flux preceded that of $B$ band by $2.0^{+0.7}_{-0.6}$ days ($1\sigma$ error). Although this result supports the X-ray reprocessing model, the derived lag is too large to be explained by the standard view which assumes a "lamppost"-type X-ray illuminator located near a standard accretion disk. Our results are better explained by assuming a hot accretion flow and a truncated disk., Comment: 16 pages, 12 figures, accepted for publication in The Astrophysical Journal

We present our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We carried out a wide-field imaging survey with Kiso Wide Field Camera (KWFC) on the 1.05-m Kiso Schmidt telescope in Japan and a galaxy-targeted survey with Tripole5 on the B&C 61-cm telescope in New Zealand. Approximately 24 deg2 regions in total were surveyed in i-band with KWFC and 18 nearby galaxies were observed with Tripole5 in g-, r-, and i-bands 4-12 days after the gravitational wave detection. Median 5-sigma depths are i~18.9 mag for the KWFC data and g~18.9 mag, r~18.7 mag, and i~18.3 mag for the Tripole5 data. Probability for a counterpart to be in the observed area is 1.2% in the initial skymap and 0.1% in the final skymap. We do not find any transient source associated to an external galaxy with spatial offset from its center, which is consistent with the local supernova rate. We summarize future prospects and ongoing efforts to pin down electromagnetic counterparts of binary black hole mergers as well as neutron star mergers., Comment: submitted to PASJ, 7 pages, 1 figure, 4 tables

We present extensively dense observations of the super-Chandrasekhar supernova (SC SN) candidate SN 2012dn from $-11$ to $+140$ days after the date of its $B$-band maximum in the optical and near-infrared (NIR) wavelengths conducted through the OISTER ToO program. The NIR light curves and color evolutions up to 35 days after the $B$-band maximum provided an excellent match with those of another SC SN 2009dc, providing a further support to the nature of SN 2012dn as a SC SN. We found that SN 2012dn exhibited strong excesses in the NIR wavelengths from $30$ days after the $B$-band maximum. The $H$ and $K_{s}$-band light curves exhibited much later maximum dates at $40$ and $70$ days after the $B$-band maximum, respectively, compared with those of normal SNe Ia. The $H$ and $K_{s}$-band light curves subtracted by those of SN 2009dc displayed plateaued evolutions, indicating a NIR echo from the surrounding dust. The distance to the inner boundary of the dust shell is limited to be $4.8 - 6.4\times10^{-2}$ pc. No emission lines were found in its early phase spectrum, suggesting that the ejecta-CSM interaction could not occur. On the other hand, we found no signature that strongly supports the scenario of dust formation. The mass loss rate of the pre-explosion system is estimated to be $10^{-6}-10^{-5}$ M$_{\odot}$ yr$^{-1}$, assuming that the wind velocity of the system is $10-100$ km~s$^{-1}$, which suggests that the progenitor of SN 2012dn could be a recurrent nova system. We conclude that the progenitor of this SC SN could be explained by the single-degenerate scenario., Comment: 22 pages, 10 pages, 8 tables, accepted for publication in PASJ. Reference corrected

The origin of short timescale (weeks/months) variability of AGN, whether due to intrinsic disc variations or reprocessing of X-ray emission by a surrounding accretion disc, has been a puzzle for many years. However recently a number of observational programmes, particularly of NGC5548 with Swift, have shown that the UV/optical variations lag behind the X-ray variations in a manner strongly supportive of X-ray reprocessing. Somewhat surprisingly the implied size of the accretion disc is ~3x greater than expected from a standard, smooth, Shakura-Sunyaev thin disc model. Although the difference may be explained by a clumpy accretion disc, it is not clear whether the difference will occur in all AGN or whether it may change as, eg, a function of black hole mass, accretion rate or disc temperature. Measurements of interband lags for most AGN require long timescale monitoring, which is hard to arrange. However for low mass (<1 million solar mass) AGN, the combination of XMM-Newton EPIC (X-rays) with the optical monitor in fast readout mode allows an X-ray/UV-optical lag to be measured within a single long observation. Here we summarise previous related observations and report on XMM-Newton observations of NGC4395 (mass ~100x lower and accretion rate ~20x lower than for NGC5548). We find that the UVW1 lags the X-rays by ~470s. Simultaneous observations at 6 different ground based observatories also allowed the g-band lag (~800s) to be measured. These observations are in agreement with X-ray reprocessing but initial analysis suggests that, for NGC4395, they do not differ markedly from the predictions of the standard thin disc model., Comment: 6 pages, 16 figures, accepted for publication in Astronomische Nachrichten

We present infrared multi-epoch observations of the dust forming nova V1280 Sco over $\sim$2000 days from the outburst. The temporal evolution of the infrared spectral energy distributions at 1272, 1616 and 1947 days can be explained by the emissions produced by amorphous carbon dust of mass (6.6--8.7)$\times$10$^{-8}$M$_{\odot}$ with a representative grain size of 0.01$~\mu$m and astronomical silicate dust of mass (3.4--4.3)$\times$10$^{-7}$M$_{\odot}$ with a representative grain size of 0.3--0.5$~\mu$m. Both of these dust species travel farther away from the white dwarf without an apparent mass evolution throughout those later epochs. The dust formation scenario around V1280 Sco suggested from our analyses is that the amorphous carbon dust is formed in the nova ejecta followed by the formation of silicate dust in the expanding nova ejecta or as a result of the interaction between the nova wind and the circumstellar medium., Comment: 47 pages, 19 figures. Accepted for publication in the Astrophysical Journal on December 2, 2015

(abridged) We report on multi-wavelength measurements of Swift J1753.5-0127 in the hard state at L=2.7e36 erg/s (assuming d=3 kpc) in 2014. The radio emission is optically thick synchrotron, presumably from a compact jet. We take advantage of the low extinction and model the near-IR to UV emission with a multi-temperature disk model. Assuming a BH mass of M_BH=5 Msun and a system inclination of 40 deg, the fits imply an inner radius for the disk of Rin/Rg>212 d_3 (5Msun/M_BH). The outer radius is R_out/R_g=90,000 d_3 (5Msun/M_BH), which corresponds to 6.6e10 d_3 cm, consistent with the expected size of the disk. The 0.5-240 keV spectrum measured by Swift/XRT, Suzaku, and NuSTAR is relatively well characterized by a power-law with a photon index of Gamma=1.722+/-0.003, but a significant improvement is seen when a second continuum component is added. Reflection is a possibility, but no iron line is detected, implying a low iron abundance. We are able to fit the entire SED with a multi-temperature disk component, a Comptonization component, and a broken power-law, representing the emission from the compact jet. The broken power-law cannot significantly contribute to the soft X-ray emission, and this may be related to why Swift J1753.5-0127 is an outlier in the radio/X-ray correlation. The broken power-law might dominate above 20 keV, which would constrain the break frequency to be between 2.4e10 Hz and 3.6e12 Hz. Although the fits to the full SED do not include significant thermal emission in the X-ray band, previous observations have consistently seen such a component, and we find that there is evidence at the 3.1-sigma level for a disk-blackbody component with a temperature of 150(+30)(-20) eV and an inner radius of 5-14 R_g. If this component is real, it might imply the presence of an inner accretion disk in addition to the strongly truncated (R_in>212 R_g) disk., Comment: 18 pages, 9 figures, accepted by ApJ

We discuss what is an appropriate set of explanatory variables in order to predict the absolute magnitude at the maximum of Type Ia supernovae. In order to have a good prediction, the error for future data, which is called the "generalization error," should be small. We use cross-validation in order to control the generalization error and LASSO-type estimator in order to choose the set of variables. This approach can be used even in the case that the number of samples is smaller than the number of candidate variables. We studied the Berkeley supernova database with our approach. Candidates of the explanatory variables include normalized spectral data, variables about lines, and previously proposed flux-ratios, as well as the color and light-curve widths. As a result, we confirmed the past understanding about Type Ia supernova: i) The absolute magnitude at maximum depends on the color and light-curve width. ii) The light-curve width depends on the strength of Si II. Recent studies have suggested to add more variables in order to explain the absolute magnitude. However, our analysis does not support to add any other variables in order to have a better generalization error., Comment: 9 pages, 7 figures, accepted for publication in PASJ

Broadband spectrum of AGN consists of multiple components such as jet emission and accretion disk emission. Temporal correlation study is useful to understand emission components and their physical origins. We have performed optical monitoring using Kanata telescope for 4 radio galaxies and 6 radio-loud Narrow-Line Seyfert 1 (RL-NLSy1): 2 gamma-ray-loud RL-NLSy1s, 1H 0323+342 and PMN J0948+0022, and 4 gamma-ray-quiet RL-NLSy1s. From these results, it is suggested that RL-NLSy1s show a disk-dominant phase and a jet-dominant phase in the optical band, but it is not well correlated with brightness., Comment: 5 pages, 15 figures, 2014 Fermi Symposium proceedings - eConf C14102.1

We report on long-term multi-wavelength monitoring of blazar Mrk~421 from 2010 to 2011. The source exhibited extreme X-ray flares in 2010. Our research group performed optical photopolarimetric follow-up observations using the Kanata telescope. In 2010, the variability in the X-ray band was significant, while the optical and ultraviolet (UV) flux decreased gradually. Polarization properties also exhibited unique variability in 2010, suggesting the presence of systematic component of polarization and magnetic field alignment for the emergence of a new polarized emission region. In contrast, in 2011 the variability in the X-ray band was smaller, and the variability in the optical and UV bands was larger, than in 2010. To explore the reasons for these differences, spectral fitting analysis was performed via simple synchrotron-self Compton modelling; the results revealed different behaviors in terms of spectral evolution between these periods, suggesting different variability mechanisms between 2010 and 2011. In 2010, the radiation was likely the result of energy injection into the emitting regions with an aligned magnetic field. In contrast, in 2011 the superposition of different emission regions may have contributed to the low degree of observed polarization. It also implies that high-energy electron which were not accelerated to ultra-relativistic velocities were injected in 2011., Comment: 8 pages, 5 figures, 1 table. Accepted by PASJ

Blazars are thought to possess a relativistic jet that is pointing toward the direction of the Earth and the elect of relativistic beaming enhances its apparent brightness. They radiate in all wavebands from the radio to the gamma-ray bands via the synchrotron and the inverse Compton scattering process. Numerous observations are performed but the mechanism of variability, creation and composition of jets are still controversial. We performed multi-wavelength monitoring with optical polarization for 3C 66A, Mrk 421, CTA 102 and PMN J0948+0022 to investigate the mechanisms of variability and research the emission region in the relativistic jets. Consequently, an emergence of new emission component in flaring state is suggested in each object. The most significant aspect of these results is its wide range of sizes of emission regions from $10^{14}-10^{16}$ cm, which implies the model with a number of independent emission regions with variety sizes and randomly orientation. The "shock-in-jet" scenario can explain high PD and direction of PA in each objects. It might reflect the common mechanism of flares in the relativistic jets., Comment: 2014 Fermi Symposium proceedings - eConf C14102.1

We have developed a new model of the Doppler tomography using total variation minimization (DTTVM). This method can reconstruct localized and non-axisymmetric profiles having sharp edges in the Doppler map. This characteristic is emphasized in the case that the number of the input data is small. We apply this model to real data of the dwarf novae, WZ Sge in superoutburst and TU Men in quiescence. We confirmed that DTTVM can reproduce the observed spectra with a high precision. Compared with the models based on the maximum entropy method, DTTVM provides the Doppler maps that little depend on the hyper-parameter and on the presence of the absorption core. We also introduce a cross-validation method to estimate reasonable values of a hyperparameter in the model by the data itself., Comment: 12 pages, 14 figures, PASJ in press

We performed optical simultaneous dual-band (SDSS $g'$- and $i'$- band) photometry and low-resolution spectroscopy for the WZ Sge-type dwarf nova EZ Lyn during its 2010 superoutburst. Dual-band photometry revealed that the $g'-i'$ color reddened with a decrease in brightness, during the main superoutburst and the following rebrightening phase, whereas the color became bluer with a further decrease in brightness during the slow, final decline phase. With a fit to our photometric results by a blackbody function, we estimated the disk radius ratio (ratio of the disk radius to the binary separation) and compared this with that of V455 And, a WZ Sge-type object that did not show any rebrightening in the 2007 superoutburst. The comparison revealed: (1) the disk radius ratio of EZ Lyn decreased more slowly than that of V455 And; and (2) the radius ratio of EZ Lyn at the end of the main superoutburst was larger than that of the V455 And. These results favor the mass reservoir model for the mechanism of rebrightening. During both the superoutburst plateau and subsequent rebrightening phase, H$\alpha$ and H$\beta$ lines were detected. The H$\alpha$ line showed a double-peak profile from which we estimated the disk radius ratio. The comparison of this ratio with that derived by photometry, indicates that the H$\alpha$ disk was larger than the photometric one, which suggests that the optically thin gas was extended to the outer region more than the optically thick gas disk and was possibly responsible for the rebrightening phenomenon. Time-series dual-band photometry during the main superoutburst revealed that color variations during the early superhump show roughly the same behavior as that of V455 And, whereas color variations during the ordinary superhump display clear anticorrelation with brightness, in contrast to that seen in the V455 And., Comment: 28 pages, 10 figures, 4 tables. Accepted for publication in PASJ

We propose a new technique to obtain super-resolution images with radio interferometer using sparse modeling. In standard radio interferometry, sampling of ($u$, $v$) is quite often incomplete and thus obtaining an image from observed visibilities becomes an underdetermined problem, and a technique so-called "zero-padding" is often used to fill up unsampled grids in ($u$, $v$) plane, resulting in image degradation by finite beam size as well as numerous side-lobes. In this paper we show that directly solving such an underdetermined problem based on sparse modeling (in this paper LASSO) avoids the above problems introduced by zero-padding, leading to super-resolution images in which structure finer than the standard beam size (diffraction limit) can be reproduced. We present results of one-dimensional and two-dimensional simulations of interferometric imaging, and discuss its implications to super-resolution imaging, particularly focusing on imaging of black hole shadows with millimeter VLBI., Comment: 12 pages, 8 figures, accepted for publication in PASJ

We present results of optical polarimetric and multi-band photometric observations for gamma-ray loud narrow-line Seyfert 1 galaxy 1H 0323+342. This object has been monitored by 1.5 m Kanata telescope since 2012 September but following a gamma-ray flux enhancement detected by Fermi-LAT on MJD 56483 (2013 July 10) dense follow-up was performed by ten 0.5-2.0 m telescopes in Japan over one week. The 2-year R_C-band light curve showed clear brightening corresponding to the gamma-ray flux increase and then decayed gradually. The high state as a whole lasted for ~20 days, during which we clearly detected optical polarization from this object. The polarization degree (PD) of the source increased from 0-1% in quiescence to ~3% at maximum and then declined to the quiescent level, with the duration of the enhancement of less than 10 days. The moderate PD around the peak allowed us to precisely measure the daily polarization angle (PA). As a result, we found that the daily PAs were almost constant and aligned to the jet axis, suggesting that the magnetic field direction at the emission region is transverse to the jet. This implies either a presence of helical/toroidal magnetic field or transverse magnetic field compressed by shock(s). We also found small-amplitude intra-night variability during the 2-hour continuous exposure on a single night. We discuss these findings based on the turbulent multi-zone model recently advocated by Marscher (2014). Optical to ultraviolet (UV) spectrum showed a rising shape in the higher frequency and the UV magnitude measured by Swift/UVOT was steady even during the flaring state, suggesting that thermal emission from accretion disk is dominant in that band., Comment: 7 pages, 4 figures, 1 table. Accepted by PASJ

The gamma-ray-detected blazar 3C 454.3 exhibits dramatic flux and polarization variations in the optical and near-infrared bands. In December 2010, the object emitted a very bright outburst. We monitored it for approximately four years (including the 2010 outburst) by optical and near-infrared photopolarimetry. During the 2010 outburst, the object emitted two rapid, redder brightenings, at which the polarization degrees (PDs) in both bands increased significantly and the bands exhibited a frequency-dependent polarization. The observed frequency-dependent polarization leads us to propose that the polarization vector is composed of two vectors. Therefore, we separate the observed polarization vectors into short and long-term components that we attribute to the emissions of the rapid brightenings and the outburst that varied the timescale of days and months, respectively. The estimated PD of the short-term component is greater than the maximum observed PD and is close to the theoretical maximum PD. We constrain the bulk Lorentz factors and inclination angles between the jet axis and the line of sight from the estimated PDs. In this case, the inclination angle of the emitting region of short-term component from the first rapid brightening should be equal to 90$^{\circ}$, because the estimated PD of the short-term component was approximately equal to the theoretical maximum PD. Thus, the Doppler factor at the emitting region of the first rapid brightening should be equal to the bulk Lorentz factor., Comment: 8 pages, 6 figures, accepted for publication in ApJ

We study the image reconstruction problem of a Compton camera which consists of semiconductor detectors. The image reconstruction is formulated as a statistical estimation problem. We employ a bin-mode estimation (BME) and extend an existing framework to a Compton camera with multiple scatterers and absorbers. Two estimation algorithms are proposed: an accelerated EM algorithm for the maximum likelihood estimation (MLE) and a modified EM algorithm for the maximum a posteriori (MAP) estimation. Numerical simulations demonstrate the potential of the proposed methods., Comment: 18 pages, 8 figures, submitted for publication in Nuclear Inst. and Methods in Physics Research, A