We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolution of AKARI compared to the previous IRAS survey, and reliable identification from WISE MIR pointing, the sample is unique in identifying optically-faint (i$\sim$20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum-SED decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ($M_{\rm star}\sim10^{11}$-$10^{12}$ M$_{\odot}$), associated with intense star formation activities (SFR $\sim$ 200-2000 M$_{\odot}$ yr$^{-1}$). The sample covers a range of AGN bolometric luminosity of $10^{10}$-$10^{13}$ L$_{\odot}$, and the outflow velocity measured from [OIII] 5007A line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500-2000 km s$^{-1}$. However, the co-existence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving stellar mass and mass fraction of young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that SF-dominated ULIRGs show smaller stellar mass and younger stellar populations compared to AGN-dominated ULIRGs.

We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply the Lyman break technique to a large search for HzRGs. For an HzRG candidate among pre-selected $r$-band dropouts with a radio detection, a follow-up optical spectroscopy with GMOS/Gemini has been performed. The obtained spectrum presents a clear Ly$\alpha$ emission line redshifted to $z=4.72$. The SED fitting analysis with the rest-frame UV and optical photometries suggests the massive nature of this HzRG with $\log{M_*/M_{\odot}} = 11.4$. The small equivalent width of Ly$\alpha$ and the moderately red UV colors indicate its dusty host galaxy, implying a chemically evolved and dusty system. The radio spectral index does not meet a criterion for an ultra-steep spectrum: $\alpha^{325}_{1400}$ of $-1.1$ and $\alpha^{150}_{1400}$ of $-0.9$, demonstrating that the HSC-SSP survey compensates for a sub-population of HzRGs which are missed in surveys focusing on an ultra-steep spectral index.

We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray spectra and those from the infrared ones. The torus angular widths determined from the infrared spectra ($\sigma_{\mathrm{IR}}$) are systematically larger than those from the X-ray data ($\sigma_{\mathrm{X}}$); the difference ($\sigma_{\mathrm{IR}}-\sigma_{\mathrm{X}}$) correlates with the inclination angle determined from the X-ray spectrum. These results can be explained by the contribution from dusty polar outflows to the observed infrared flux, which becomes more significant at higher inclinations (more edge-on views). The ratio of the hydrogen column density and V-band extinction in the line of sight absorber shows large scatter ($\simeq$1 dex) around the Galactic value, suggesting that a significant fraction of AGNs have dust-rich circumnuclear environments.

Theory predicts that a supermassive black hole binary (SMBHB) could be observed as a luminous active galactic nucleus (AGN) that periodically varies on the order of its orbital timescale. In X-rays, periodic variations could be caused by mechanisms including relativistic Doppler boosting and shocks. Here we present the first systematic search for periodic AGNs using $941$ hard X-ray light curves (14-195 keV) from the first 105 months of the Swift Burst Alert Telescope (BAT) survey (2004-2013). We do not find evidence for periodic AGNs in Swift-BAT, including the previously reported SMBHB candidate MCG+11$-$11$-$032. We find that the null detection is consistent with the combination of the upper-limit binary population in AGNs in our adopted model, their expected periodic variability amplitudes, and the BAT survey characteristics. We have also investigated the detectability of SMBHBs against normal AGN X-ray variability in the context of the eROSITA survey. Under our assumptions of a binary population and the periodic signals they produce which have long periods of hundreds of days, up to $13$% true periodic binaries can be robustly distinguished from normal variable AGNs with the ideal uniform sampling. However, we demonstrate that realistic eROSITA sampling is likely to be insensitive to long-period binaries because longer observing gaps reduce their detectability. In contrast, large observing gaps do not diminish the prospect of detecting binaries of short, few-day periods, as 19% can be successfully recovered, the vast majority of which can be identified by the first half of the survey.

The vast majority of supermassive black holes (SMBHs) in the local universe exhibit levels of activity much lower than those expected from gas supplying rates onto the galactic nuclei, and only a small fraction of silent SMBHs can turn into active galactic nuclei. Revisiting observational data of very nearby SMBHs whose gravitational spheres of influence are spatially reached by the Chandra X-ray satellite, we find that the level of BH activity drastically increases from the quiescent phase when the inflow rate outside of the BH influence radius is higher than 0.1% of the Eddington accretion rate. We also show that the relation between the nuclear luminosity and gas accretion rate from the BH influence radius measured from X-ray observations is well described by the universal state transition of accreting SMBHs, as predicted by recent hydrodynamical simulations with radiative cooling and BH feedback. After the state transition, young massive stars should form naturally in the nucleus, as observed in the case of the nearest SMBH, Sagittarius A$^\ast$, which is currently quiescent but was recently active.

The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows, is not thoroughly studied. Here we report the first detection of polarized thermal emission by means of magnetically aligned dust grains in the dusty torus of NGC 1068 using ALMA Cycle 4 polarimetric dust continuum observations ($0.07"$, $4.2$ pc; 348.5 GHz, $860$ $\mu$m). The polarized torus has an asymmetric variation across the equatorial axis with a peak polarization of $3.7\pm0.5$\% and position angle of $109\pm2^{\circ}$ (B-vector) at $\sim8$ pc east from the core. We compute synthetic polarimetric observations of magnetically aligned dust grains assuming a toroidal magnetic field and homogeneous grain alignment. We conclude that the measured 860 $\mu$m continuum polarization arises from magnetically aligned dust grains in an optically thin region of the torus. The asymmetric polarization across the equatorial axis of the torus arises from 1) an inhomogeneous optical depth, and 2) a variation of the velocity dispersion, i.e. variation of the magnetic field turbulence at sub-pc scales, from the eastern to the western region of the torus. These observations and modeling constrain the torus properties beyond spectral energy distribution results. This study strongly supports that magnetic fields up to a few pc contribute to the accretion flow onto the active nuclei.

Observed extinction curves of active galactic nuclei (AGNs) are significantly different from those observed in the Milky Way. The observations require preferential removal of small grains at the AGN environment; however, the physics for this remains unclear. In this paper, we propose that dust destruction by charging, or Coulomb explosion, may be responsible for AGN extinction curves. Harsh AGN radiation makes a dust grain highly charged through photoelectric emission, and grain fission via the Coulomb explosion occurs when the electrostatic tensile stress of a charge grain exceeds its tensile strength. We show that the Coulomb explosion can preferentially remove both small silicate and graphite grains and successfully reproduce both flat extinction curves and the absence of 2175\AA~bump.

Recent mid-infrared high spatial resolution observations have revealed that active galactic nuclei (AGNs) may host a polar dust region with the size of several pc, and such dust may be carried by radiation from the central engine. The polar dust emission often exhibits very weak or absence of the silicate 10-$\mu$m emission feature. A possible explanation is that the polar dust is dominated by micron-sized large grains because these grains do not show the silicate feature, while it remains unclear how large grains are preferentially supplied to the polar region. We here propose a new scenario describing the prevalence of large grains at the polar region. We show that grains are accelerated to the hypersonic drift velocity by the radiation pressure from AGN, and the hypersonic drift results in dust destruction via kinetic sputtering. Sputtering destroys small grains faster than the large ones, and thus larger grains will be preferentially blown over longer distance. Although the hypersonic drift, or kinetic sputtering, tends to be suppressed for very small grains due to the Coulomb drag, they might also be disrupted by Coulomb explosion. Removal of small grains and/or survival of large grains may explain the lack of a silicate 10-$\mu$m emission feature in polar dust emission.

We study the obscuring structure of circumnuclear disks (CNDs) by considering supernova (SN) feedbacks from nuclear starburst and the effect of anisotropic radiative pressure from AGNs. We suppose that the mass accretion onto a central supermassive black hole (SMBH) is triggered by SN-driven turbulence within CNDs, and we explore how the structures of CNDs depend on the BH mass ($M_{\rm BH}$) and AGN luminosity ($L_{\rm AGN}$). We find that the obscuring fraction ($f_{\rm obs}$) peaks at $\sim10\%$ of the Eddington luminosity ($L_{\rm Edd}$), and its maximal value is $f_\mathrm{obs} \sim 0.6$ for less massive SMBHs (e.g., $M_{\rm BH} < 10^{8}M_{\odot}$). This is because the scale height of CNDs is determined by the SN-driven accretion for a smaller $L_{\rm AGN}$, while the dusty molecular gas in CNDs is blown away by the radiation pressure from AGNs beyond the critical luminosity. On the other hand, for massive SMBHs (e.g., $M_{\rm BH} > 10^{8}M_{\odot}$), $f_{\rm obs}$ is always smaller than $0.2$, and it is almost independent of $L_{\rm AGN}$ because the scale height of CNDs is mainly controlled by the maximal star-formation efficiency ($C_{\rm *, max}$) in CNDs. By comparison with the obscuring fractions suggested from the mid-infrared observations of nearby AGNs, the SN plus radiative feedback model with $C_{\rm *, max} = 10^{-7}\, {\rm yr}^{-1}$ well reproduces the observations for $M_\mathrm{BH} = 10^8 M_\odot$. We also find that the intense starburst or the existence of dust-free absorbers inside CNDs are necessary, to explain X-ray observations.

Kuiper belt objects (KBOs) are thought to be the remnant of the early solar system, and their size distribution provides an opportunity to explore the formation and evolution of the outer solar system. In particular, the size distribution of kilometre-sized (radius = 1-10 km) KBO represents a signature of initial planetesimal sizes when planets form. These kilometre-sized KBOs are extremely faint, and it is impossible to detect them directly. Instead, monitoring of stellar occultation events is one possible way to discover these small KBOs. Hitherto, however, there has been no observational evidence for the occultation events by KBOs with radii of 1-10 km. Here we report the first detection of a single occultation event candidate by a KBO with a radius of $\sim$1.3 km, which is simultaneously provided by two low-cost small telescopes coupled with commercial CMOS cameras. From this detection, we conclude that a surface number density of KBOs with radii exceeding $\sim 1.2$ km is $\sim 6 \times 10^5 \ {\rm deg^{-2}}$. This surface number density favours a theoretical size distribution model with an excess signature at a radius of 1-2 km. If this is a true detection, this implies that planetesimals before their runaway growth phase grow into kilometre-sized objects in the primordial outer solar system and remain as a major population of the present-day Kuiper belt.

Recent active galactic nucleus (AGN) and quasar surveys have revealed a population showing rapid AGN luminosity variability by a factor of $\sim10$. Here we present the most drastic AGN luminosity decline by a factor of $\gtrsim 10^{3}$ constrained by a NuSTAR X-ray observation of the nearby galaxy Arp 187, which is a promising "dead" quasar whose current activity seems quiet but whose past activity of $L_\mathrm{bol} \sim 10^{46}$ erg s$^{-1}$ is still observable at a large scale by its light echo. The obtained upper bound of the X-ray luminosity is $\log (L_{\rm 2-10 keV}/{\rm erg} {\rm s}^{-1}) < 41.2$, corresponding to $\log (L_\mathrm{bol}/{\rm erg} {\rm s}^{-1}) < 42.5$, indicating an inactive central engine. Even if a putative torus model with $N_\mathrm{H} \sim 1.5 \times 10^{24}$ cm$^{-2}$ is assumed, the strong upper-bound still holds with $\log (L_{\rm 2-10 keV}/{\rm erg} {\rm s}^{-1}) < 41.8$ or $\log (L_\mathrm{bol}/{\rm erg} {\rm s}^{-1}) < 43.1$. Given the expected size of the narrow line region, this luminosity decrease by a factor of $\gtrsim 10^3$ must have occurred within $\lesssim 10^4$ yr. This extremely rapid luminosity/accretion shutdown is puzzling and it requires one burst-like accretion mechanism producing a clear outer boundary for an accretion disk. We raise two possible scenarios realizing such an accretion mechanism: a mass accretion 1) by the tidal disruption of a molecular cloud and/or 2) by the gas depletion as a result of vigorous nuclear starformation after rapid mass inflow to the central engine.

We present a multi wavelength analysis of 28 of the most luminous low-redshift narrow-line, ultra-hard X-ray selected active galactic nuclei (AGN) drawn from the 70 month Swift/BAT all-sky survey, with bolometric luminosities of log(L_bol/erg/s) > 45.25. The broad goal of our study is to determine whether these objects have any distinctive properties, potentially setting them aside from lower-luminosity obscured AGN in the local Universe. Our analysis relies on the first data release of the BAT AGN Spectroscopic Survey (BASS/DR1) and on dedicated observations with the VLT, Palomar, and Keck observatories. We find that the vast majority of our sources agree with commonly used AGN selection criteria which are based on emission line ratios and on mid-infrared colours. Our AGN are predominantly hosted in massive galaxies (9.8 < log(M_*/M_sun) < 11.7); based on visual inspection of archival optical images, they appear to be mostly ellipticals. Otherwise, they do not have distinctive properties. Their radio luminosities, determined from publicly available survey data, show a large spread of almost 4 orders of magnitude - much broader than what is found for lower X-ray luminosity obscured AGN in BASS. Moreover, our sample shows no preferred combination of black hole masses (M_BH) and/or Eddington ratio (lambda_Edd), covering 7.5 < log(M_BH/M_sun) < 10.3 and 0.01 < lambda_Edd < 1. Based on the distribution of our sources in the lambda_Edd-N_H plane, we conclude that our sample is consistent with a scenario where the amount of obscuring material along the line of sight is determined by radiation pressure exerted by the AGN on the dusty circumnuclear gas.

We systematically investigate the mid-infrared (MIR; $\lambda>3 ~\mu$m) time variability of uniformly selected $\sim800$ massive young stellar objects (MYSOs) from the Red MSX Source (RMS) survey. Out of the 806 sources, we obtain reliable 9-year-long MIR magnitude variability data of 331 sources at the 3.4~$\mu$m (W1) and 4.6~$\mu$m (W2) bands by cross-matching the MYSO positions with ALLWISE and NEOWISE catalogs. After applying the variability selections using ALLWISE data, we identify 5 MIR-variable candidates. The light curves show various classes, with the periodic, plateau-like, and dipper features. Out of the obtained two color-magnitude diagram of W1 and W1$-$W2, one shows "bluer when brighter and redder when fainter" trends in variability, suggesting change in extinction or accretion rate. Finally, our results show that G335.9960$-$00.8532 (hereafter, G335) has a periodic light curve, with a $\approx 690$-day cycle. Spectral energy density model fitting results indicate that G335 is a relatively evolved MYSO; thus, we may be witnessing the very early stages of a hyper- or ultra-compact HII region, a key source for understanding MYSO evolution.

We explore X-ray evidence for the presence of active galactic nuclei (AGNs) in the two most actively star-forming Green Pea galaxies (GPs), SDSS J0749+3337 and SDSS J0822+2241, which have star-formation rates (SFRs) of $123~M_\odot$ yr$^{-1}$ and $78~M_\odot$ yr$^{-1}$, respectively. The GPs have red mid-infrared (MIR) spectral energy distributions and higher 22 $\mu$m luminosities than expected from a proxy of the SFR (H$\alpha$ luminosity), consistent with hosting AGNs with 2-10 keV luminosities of $\sim10^{44}$ erg s$^{-1}$. We thus obtain and analyze the first hard ($>$ 10 keV) X-ray data observed with NuSTAR and archival XMM-Newton data below 10 keV. From the NuSTAR $\approx$20 ksec data, however, we find no significant hard X-ray emission. By contrast, soft X-ray emission with 0.5--8 keV luminosities of $\approx10^{42}$ erg s$^{-1}$ is significantly detected in both targets, which can be explained only by star formation (SF). A possible reason for the lack of clear evidence is that a putative AGN torus absorbs most of the X-ray emission. Applying a smooth-density AGN torus model, we determine minimum hydrogen column densities along the equatorial plane ($N_{\rm H}^{\rm eq}$) consistent with the non-detection. The results indicate $N_{\rm H}^{\rm eq} \gtrsim 2\times10^{24}$ cm$^{-2}$ for SDSS J0749+3337 and $N_{\rm H}^{\rm eq} \gtrsim 5\times10^{24}$ cm$^{-2}$ for SDSS J0822+2241. Therefore, the GPs may host such heavily obscured AGNs. Otherwise, no AGN exists and the MIR emission is ascribed to SF. Active SF in low-mass galaxies is indeed suggested to reproduce red MIR colors. This would imply that diagnostics based on MIR photometry data alone may misidentify such galaxies as AGNs.

We study the properties of rotating accretion flows onto supermassive black holes (SMBHs) using axisymmetric two-dimensional hydrodynamical simulations with radiative cooling and BH feedback. The simulations resolve the accretion dynamics of gas outside from the BH influence radius through an inner accretion disk. For lower Bondi accretion rates in units of the Eddington rate ($\dot{M}_{\rm B}\ll 10^{-3}~\dot{M}_{\rm Edd}$), the BH feeding is suppressed due to turbulent motion by several orders of magnitudes from the Bondi rate. Thus, the radiative luminosity results in as low as $\sim 10^{-10}-10^{-7}~L_{\rm Edd}$, where $L_{\rm Edd}$ is the Eddington luminosity. For higher rates of $\dot{M}_{\rm B}> 10^{-3}~\dot{M}_{\rm Edd}$, the optically-thin accreting gas cools via free-free emission and forms a geometrically-thin disk, which feeds the BH efficiently and increases the radiative luminosity to $> 10^{-3}~L_{\rm Edd}$. The transitional behavior of accreting BHs in galactic nuclei from radiatively inefficient phases to cold disk accretion naturally explains (1) the reason for the offset between the observed luminosities and theoretical predictions for nearby quiescent SMBHs, and (2) the conditions to fuel gas into the nuclear SMBH. In addition, the cold disk formed in galactic nuclei tends to be gravitationally unstable and leads to star formation when the Bondi rate is as high as $ \dot{M}_{\rm B} > 10^{-2}~M_\odot~{\rm yr}^{-1}$. This is a plausible explanation of the correlation observed between star formation rates and BH feeding rates in Seyfert galaxies.

We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We derived their color excess ($E (B-V)_{*}$), stellar mass, star formation rate (SFR), IR luminosity, the ratio of IR and radio luminosity ($q_{\rm IR}$), and radio spectral index ($\alpha_{\rm radio}$) that are derived from the SED fitting with CIGALE. We also estimated Eddington ratio based on stellar mass and integration of the best-fit SEDs of AGN component. We found that $E (B-V)_{*}$, SFR, and IR luminosity clearly depend on redshift while stellar mass, $q_{\rm IR}$, and $\alpha_{\rm radio}$ do not significantly depend on redshift. Since optically-faint ($i_{\rm AB} \geq 21.3$) RGs that are newly discovered by our RG survey tend to be high redshift, they tend to not only have a large dust extinction and low stellar mass but also have high SFR and AGN luminosity, high IR luminosity, and high Eddington ratio compared to optically-bright ones. The physical properties of a fraction of RGs in our sample seem to differ from a classical view of RGs with massive stellar mass, low SFR, and low Eddington ratio, demonstrating that our RG survey with HSC and FIRST provides us curious RGs among entire RG population.

In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$\mu$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a ULIRG at z = 0.49 in the spectroscopic follow-up observation, indicates signatures of an extremely strong outflow in its emission line profiles. Its [OIII] 5007\AA\ emission line shows FWHM of 1830 km s$^{-1}$ and velocity shift of -770 km s$^{-1}$ in relative to the stellar absorption lines. Furthermore, low-ionization [OII] 3726\AA\ 3729\AA\ doublet also shows large FWHM of 910 km s$^{-1}$ and velocity shift of -380 km s$^{-1}$. After the removal of an unresolved nuclear component, the long-slit spectroscopy 2D image possibly shows that the outflow extends to radius of 4 kpc. The mass outflow and energy ejection rates are estimated to be 500 M$_{\odot}$ yr$^{-1}$ and $4\times10^{44}$ erg s$^{-1}$, respectively, which imply that the outflow is among the most powerful ones observed in ULIRGs and QSOs at 0.3 < z < 1.6. The co-existence of the strong outflow and intense star formation (star formation rate of 990 M$_{\odot}$ yr$^{-1}$) indicates that the feedback of the strong outflow has not severely affect the star-forming region of the galaxy.

We present arcsecond-scale observations of the active galactic nuclei (AGNs) of seven nearby Seyfert galaxies observed from the Stratospheric Observatory For Infrared Astronomy (SOFIA) using the 31.5 and 37.1 um filters of the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST). We isolate unresolved emission from the torus and find extended diffuse emission in six 37.1 um residual images in our sample. Using Spitzer/IRS spectra, we determine the dominant mid-infrared (MIR) extended emission source and attribute it to dust in the narrow line region (NLR) or star formation. We compare the optical NLR and radio jet axes to the extended 37.1 um emission and find coincident axes for three sources. For those AGNs with extended emission coincident with the optical axis, we find that spatial scales of the residual images are consistent with 0.1 - 1 kpc scale distances to which dust can be heated by the AGN. Using previously published subarcsecond 1 - 20 um imaging and spectroscopic data along with our new observations, we construct broadband spectral energy distributions (SEDs) of the AGNs at wavelengths 1 - 40 um. We find that three AGNs in our sample tentatively show a turnover in the SED between 30 - 40 um. Using results from Clumpy torus models and the Bayesian inference tool BayesClumpy, we find that the posterior outputs for AGNs with MIR turnover revealed by SOFIA/FORCAST have smaller uncertainties than AGNs that do not show a turnover.

Arp 187 is one of the fading active galactic nuclei (AGN), whose AGN activity is currently decreasing in luminosity. We investigate the observational signatures of AGN in Arp 187, which trace various physical scales from less than 0.1 pc to the nearly 10 kpc, to estimate the longterm luminosity change over $10^{4}$ years. The VLA 5 GHz, 8 GHz, and the ALMA 133 GHz images reveal bimodal jet lobes with $\sim$5 kpc size and the absence of the central radio-core. The 6dF optical spectrum shows that Arp 187 hosts narrow line region with the estimated size of $\sim$1 kpc, and the line strengths give the AGN luminosity of $L_{\rm bol}=1.5 \times 10^{46}$ erg s$^{-1}$. On the other hand, the current AGN activity estimated from the AGN torus emission gives the upper bound of $L_{\rm bol} < 2.2 \times 10^{43}$ erg s$^{-1}$. The absence of the radio-core gives the more strict upper bound of the current AGN luminosity of $L_{\rm bol} < 8.0 \times 10^{40}$ erg s$^{-1}$, suggesting that the central engine is already quenched. These multi-wavelength signatures indicate that Arp 187 hosts a "dying" AGN: the central engine is already dead, but the large scale AGN indicators are still observable as the remnant of the past AGN activity. The central engine has experienced the drastic luminosity decline by a factor of $\sim10^{3-5}$ fainter within $\sim10^{4}$ years, which is roughly consistent with the viscous timescale of the inner part of the accretion disk within $\sim$500 years.

We quantify the luminosity contribution of active galactic nuclei (AGN) to the 12 $\mu$m, mid-infrared (MIR; 5-38 $\mu$m), and the total IR (5-1000 $\mu$m) emission in the local AGN detected in the all-sky 70-month Swift/Burst Alert Telescope (BAT) ultra hard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into AGN and starburst components using AGN torus and star-forming galaxy templates. This enables us to recover the AGN torus emission also for low-luminosity end, down to $\log (L_{14-150}/{\rm erg}~{\rm s}^{-1}) \simeq 41$, which typically have significant host galaxy contamination. We find that the luminosity contribution of the AGN to the 12 $\mu$m, the MIR, and the total IR band is an increasing function of the 14-150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 $\mu$m. The obtained total IR AGN luminosity through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median of the dust covering factor is always smaller than that of the X-ray obscuration fraction above the AGN bolometric luminosity of $\log (L_{\rm bol}/{\rm erg}~{\rm s}^{-1}) \simeq 42.5$. Considering that X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, it indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission.

Gravitational waves (GWs) in the nano-hertz band are great tools for understanding the cosmological evolution of supermassive black holes (SMBHs) in galactic nuclei. We consider SMBH binaries in high-$z$ ultra-luminous infrared galaxies (ULIRGs) as sources of a stochastic GW background (GWB). ULIRGs are likely associated with gas-rich galaxy mergers containing SMBHs that possibly occur at most once in the life of galaxies, unlike multiple dry mergers at low redshift. Adopting a well-established sample of ULIRGs, we study the properties of the GWB due to coalescing binary SMBHs in these galaxies. Since the ULIRG population peaks at $z>1.5$, the amplitude of the GWB is not affected even if BH mergers are delayed by as long as $\sim $ 10 Gyrs. Despite the rarity of the high-$z$ ULIRGs, we find a tension with the upper limits from Pulsar Timing Array (PTA) experiments. This result suggests that if a fraction $f_{\rm m,gal}$ of ULIRGs are associated with SMBH binaries, then no more than $20 f_{\rm m,gal}(\lambda_{\rm Edd}/0.3)^{5/3}(t_{\rm life}/30~{\rm Myr})~\%$ of the binary SMBHs in ULIRGs can merge within a Hubble time, for plausible values of the Eddington ratio of ULIRGs ($\lambda_{\rm Edd}$) and their lifetime ($t_{\rm life}$).

Dwarf spheroidal galaxies (dSphs) are promising targets for the gamma-ray dark matter (DM) search. In particular, DM annihilation signal is expected to be strong in some of the recently discovered nearby ultra-faint dSphs, which potentially give stringent constraints on the $\mathcal{O}(1)$ TeV WIMP DM. However, various non-negligible systematic uncertainties complicate the estimation of the astrophysical factors relevant for the DM search in these objects. Among them, the effects of foreground stars particularly attract attention because the contamination is unavoidable even for the future kinematical survey. In this article, we assess the effects of the foreground contamination on the astrophysical $J$-factor estimation by generating mock samples of stars in the four ultra-faint dSphs and using a model of future spectrographs. We investigate various data cuts to optimize the quality of the data and apply a likelihood analysis which takes member and foreground stellar distributions into account. We show that the foreground star contaminations in the signal region (the region of interest) and their statistical uncertainty can be estimated by interpolating the foreground star distribution in the control region where the foreground stars dominate the member stars. Such regions can be secured at future spectroscopic observations utilizing a multiple object spectrograph with a large field of view; e.g. the Prime Focus Spectrograph mounted on Subaru Telescope. The above estimation has several advantages: The data-driven estimation of the contamination makes the analysis of the astrophysical factor stable against the complicated foreground distribution. Besides, foreground contamination effect is considered in the likelihood analysis.

We combine new (NGC 1275, NGC 4151, and NGC 5506) and previously published (Cygnus A, Mrk 231, and NGC 1068) sub-arcsecond resolution mid-infrared (MIR; 8-13 $\mu$m) imaging- and spectro-polarimetric observations of six Seyfert galaxies using CanariCam on the 10.4-m Gran Telescopio CANARIAS. These observations reveal a diverse set of physical processes responsible for the nuclear polarization, and permit characterization of the origin of the MIR nuclear polarimetric signature of active galactic nuclei (AGN). For all radio quiet objects, we found that the nuclear polarization is low (<1 per cent), and the degree of polarization is often a few per cent over extended regions of the host galaxy where we have sensitivity to detect such extended emission (i.e., NGC 1068 and NGC 4151). We suggest that the higher degree of polarization previously found in lower resolution data arises only on the larger-than-nuclear scales. Only the radio-loud Cygnus A exhibits significant nuclear polarization ($\sim$11 per cent), attributable to synchrotron emission from the pc-scale jet close to the core. We present polarization models that suggest that the MIR nuclear polarization for highly obscured objects arises from a self-absorbed MIR polarized clumpy torus and/or dichroism from the host galaxy, while for unabsorbed cores, MIR polarization arises from dust scattering in the torus and/or surrounding nuclear dust.

We present observations of NGC 1068 covering the $19.7-53.0$ $\mu$m wavelength range using FORCAST and HAWC+ onboard SOFIA. Using these observations, high-angular resolution infrared (IR) and sub-mm observations, we find an observational turn-over of the torus emission in the $30-40$ $\mu$m wavelength range with a characteristic temperature of $70-100$ K. This component is clearly different from the diffuse extended emission in the narrow line and star formation regions at 10-100 $\mu$m within the central 700 pc. We compute $2.2-432$ $\mu$m 2D images using the best inferred \textsc{clumpy} torus model based on several nuclear spectral energy distribution (SED) coverages. We find that when $1-20$ $\mu$m SED is used, the inferred result gives a small torus size ($<4$ pc radius) and a steep radial dust distribution. The computed torus using the $1-432$ $\mu$m SED provides comparable torus sizes, $5.1^{+0.4}_{-0.4}$ pc radius, and morphology to the recently resolved 432 $\mu$m ALMA observations. This result indicates that the $1-20$ $\mu$m wavelength range is not able to probe the full extent of the torus. The characterization of the turn-over emission of the torus using the $30-60$ $\mu$m wavelength range is sensitive to the detection of cold dust in the torus. The morphology of the dust emission in our 2D image at 432 $\mu$m is spatially coincident with the cloud distribution, while the morphology of the emission in the $1-20$ $\mu$m wavelength range shows an elongated morphology perpendicular to the cloud distribution. We find that our 2D \textsc{clumpy} torus image at 12 $\mu$m can produce comparable results to those observed using IR interferometry.

Hard X-ray ($\geq 10$ keV) observations of Active Galactic Nuclei (AGN) can shed light on some of the most obscured episodes of accretion onto supermassive black holes. The 70-month Swift/BAT all-sky survey, which probes the 14-195 keV energy range, has currently detected 838 AGN. We report here on the broad-band X-ray (0.3-150 keV) characteristics of these AGN, obtained by combining XMM-Newton, Swift/XRT, ASCA, Chandra, and Suzaku observations in the soft X-ray band ($\leq 10$ keV) with 70-month averaged Swift/BAT data. The non-blazar AGN of our sample are almost equally divided into unobscured ($N_{\rm H}< 10^{22}\rm cm^{-2}$) and obscured ($N_{\rm H}\geq 10^{22}\rm cm^{-2}$) AGN, and their Swift/BAT continuum is systematically steeper than the 0.3-10 keV emission, which suggests that the presence of a high-energy cutoff is almost ubiquitous. We discuss the main X-ray spectral parameters obtained, such as the photon index, the reflection parameter, the energy of the cutoff, neutral and ionized absorbers, and the soft excess for both obscured and unobscured AGN.

We present the cross-correlation between 151 luminous quasars ($M_{ \mathrm{UV}} < -26$) and 179 protocluster candidates at $z \sim 3.8$, extracted from the Wide imaging survey ($ \sim 121~ $deg$^2$) performed with a part of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We find that only two out of 151 quasars reside in regions that are more overdense compared to the average field at $ > 4 \sigma $. The distributions of the distance between quasars and the nearest protoclusters and the significance of the overdensity at the position of quasars are statistically identical to those found for $g$-dropout galaxies, suggesting that quasars tend to reside in almost the same environment as star-forming galaxies at this redshift. Using stacking analysis, we find that the average density of $g$-dropout galaxies around quasars is slightly higher than that around $g$-dropout galaxies on $1.0 - 2.5$ pMpc scales, while at $ < 0.5$ pMpc that around quasars tends to be lower. We also find that quasars with higher UV-luminosity or with more massive black holes tend to avoid the most overdense regions, and that the quasar near zone sizes are anti-correlated with overdensity. These findings are consistent with a scenario in which the luminous quasar at $z \sim4 $ resides in structures that are less massive than those expected for the progenitors of today's rich clusters of galaxies, and possibly that luminous quasars may be suppressing star formation in their close vicinity.

We have shown that Lyman-$\alpha$ blobs (LABs) may still exist even at $z\sim0.3$, about 7 billion years later than most other LABs known (Schirmer et al. 2016). Their luminous Ly$\alpha$ and [OIII] emitters at $z\sim0.3$ offer new insights into the ionization mechanism. This paper focuses on the two X-ray brightest LABs at $z\sim0.3$, SDSS J0113$+$0106 (J0113) and SDSS J1155$-$0147 (J1155), comparable in size and luminosity to `B1', one of the best-studied LABs at $z \gtrsim$ 2. Our NuSTAR hard X-ray (3--30 keV) observations reveal powerful active galactic nuclei (AGN) with $L_{2-10{\;\rm keV}}=(0.5$--$3)\times10^{44}$ erg cm$^{-2}$ s$^{-1}$. J0113 also faded by a factor of $\sim 5$ between 2014 and 2016, emphasizing that variable AGN may cause apparent ionization deficits in LABs. Joint spectral analyses including Chandra data constrain column densities of $N_{\rm H}=5.1^{+3.1}_{-3.3}\times10^{23}$ cm$^{-2}$ (J0113) and $N_{\rm H}=6.0^{+1.4}_{-1.1}\times10^{22}$ cm$^{-2}$ (J1155). J0113 is likely buried in a torus with a narrow ionization cone, but ionizing radiation is also leaking in other directions as revealed by our Gemini/GMOS 3D spectroscopy. The latter shows a bipolar outflow over $10$ kpc, with a peculiar velocity profile that is best explained by AGN flickering. X-ray analysis of J1155 reveals a weakly absorbed AGN that may ionize over a wide solid angle, consistent with our 3D spectra. Extinction corrected [OIII] log-luminosities are high, $\sim43.6$. The velocity dispersions are low, $\sim100$--$150$ km s$^{-1}$, even at the AGN positions. We argue that this is a combination of high extinction hiding the turbulent gas, and previous outflows that have cleared the escape paths for their successors.

We present the first catalog and data release of the Swift-BAT AGN Spectroscopic Survey (BASS). We analyze optical spectra of the majority of AGN (77%, 641/836) detected based on their 14-195 keV emission in the 70-month Swift BAT all-sky catalog. This includes redshift determination, absorption and emission line measurements, and black hole mass and accretion rate estimates for the majority of obscured and un-obscured AGN (74%, 473/641) with 340 measured for the first time. With ~90% of sources at z<0.2, the survey represents a significant census of hard-X-ray selected AGN in the local universe. In this first catalog paper, we describe the spectroscopic observations and datasets, and our initial spectral analysis. The FWHM of the emission lines show broad agreement with the X-ray obscuration (~94%), such that Sy 1-1.8 have NH<10^21.9 cm^-2, and Seyfert 2, have NH>10^21.9 cm^-2. Seyfert 1.9 show a range of column densities. Compared to narrow line AGN in the SDSS, the X-ray selected AGN have a larger fraction of dusty host galaxies suggesting these types of AGN are missed in optical surveys. Using the most sensitive [OIII]/Hbeta and [NII]/Halpha emission line diagnostic, about half of the sources are classified as Seyferts, ~15% reside in dusty galaxies that lack an Hbeta detection, but for which the line upper limits imply either a Seyfert or LINER, ~15% are in galaxies with weak or no emission lines despite high quality spectra, and a few percent each are LINERS, composite galaxies, HII regions, or in known beamed AGN.

We estimate the dust torus cooling timescale once the active galactic nucleus (AGN) is quenched. In a clumpy torus system, once the incoming photons are suppressed, the cooling timescale of one clump from $T_{\rm dust}=1000$ K to several $10$ K is less than $10$ years, indicating that the dust torus cooling time is mainly governed by the light crossing time of the torus from the central engine. After considering the light crossing time of the torus, the AGN torus emission at $12~\mu$m becomes over two orders of magnitude fainter within $100$ years after the quenching. We also propose that those "dying" AGN could be found using the AGN indicators with different physical scale $R$ such as $12~\mu$m band luminosity tracing AGN torus ($R \sim 10$ pc) and the optical [OIII]$\lambda5007$ emission line narrow line regions ($R=10^{2-4}$ pc).

Organized Autotelescopes for Serendipitous Event Survey (OASES) is an optical observation project that aims to detect and investigate stellar occultation events by kilometer-sized trans-Neptunian objects (TNOs). In this project, multiple low-cost observation systems for wide-field and high-speed photometry were developed in order to detect rare and short-timescale stellar occultation events. The observation system consists of commercial off-the-shelf $0.28 \ {\rm m}$ aperture $f/1.58$ optics providing a $2.3 \times 1.8$ square-degree field of view. A commercial CMOS camera is coupled to the optics to obtain full-frame imaging with a frame rate greater than $10 \ {\rm Hz}$. As of September 2016, this project exploits two observation systems, which are installed on Miyako Island, Okinawa, Japan. Recent improvements in CMOS technology in terms of high-speed imaging and low readout noise mean that the observation systems are capable of monitoring $\sim 2000$ stars in the Galactic plane simultaneously with magnitudes down to ${\rm V} \sim 13.0$, providing $\sim 20\%$ photometric precision in light curves with a sampling cadence of $15.4 \ {\rm Hz}$. This number of monitored stars is larger than for any other existing instruments for coordinated occultation surveys. In addition, a precise time synchronization method needed for simultaneous occultation detection is developed using faint meteors. The two OASES observation systems are executing coordinated monitoring observations of a dense stellar field in order to detect occultations by kilometer-sized TNOs for the first time.

The gamma-ray observation of the dwarf spheroidal galaxies (dSphs) is a promising approach to search for the dark matter annihilation (or decay) signal. The dSphs are the nearby satellite galaxies with a clean environment and dense dark matter halo so that they give stringent constraints on the 1 TeV dark matter. However, recent studies have revealed that current estimation of astrophysical factors relevant for the dark matter searches are not conservative, where the various non-negligible systematic uncertainties are not taken into account. Among them, the effect of foreground stars on the astrophysical factors has not been paid much attention, which becomes more important for deeper and wider stellar surveys in the future. In this article, we assess the effects of the foreground contamination by generating the mock samples of stars and using a model of future spectrographs. We investigate various data cuts to optimize the quality of the data and find that the cuts on the velocity and surface gravity can efficiently eliminate the contamination. We also propose a new likelihood function which includes the foreground distribution function. We apply this likelihood function to the fit of the three types of the mock data (Ursa Minor, Draco with large dark matter halo, and Draco with small halo) and three cases of the observation. The likelihood successfully reproduces the input $J$-factor value while the fit without considering the foreground distribution gives large deviation from the input value by a factor of three.

Recent quasar surveys have revealed that supermassive black holes (SMBHs) rarely exceed a mass of $M_{\rm BH} \sim {\rm a~few}\times10^{10}~M_{\odot}$ during the entire cosmic history. It has been argued that quenching of the BH growth is caused by a transition of a nuclear accretion disk into an advection dominated accretion flow, with which strong outflows and/or jets are likely to be associated. We investigate a relation between the maximum mass of SMBHs and the radio-loudness of quasars with a well-defined sample of $\sim 10^5$ quasars at a redshift range of $0<z<2$, obtained from the Sloan Digital Sky Surveys DR7 catalog. We find that the number fraction of the radio-loud (RL) quasars increases above a threshold of $M_{\rm BH} \simeq 10^{9.5}~M_{\odot}$, independent of their redshifts. Moreover, the number fraction of RL quasars with lower Eddington ratios (out of the whole RL quasars), indicating lower accretion rates, increases above the critical BH mass. These observational trends can be natural consequences of the proposed scenario of suppressing BH growth around the apparent maximum mass of $\sim 10^{10}~M_{\odot}$. The ongoing VLA Sky Survey in radio will allow us to estimate of the exact number fraction of RL quasars more precisely, which gives further insights to understand quenching processes for BH growth.

We present the results of clustering analyses of Lyman break galaxies (LBGs) at $z\sim3$, $4$, and $5$ using the final data release of the Canada--France--Hawaii Telescope Legacy Survey (CFHTLS). Deep- and wide-field images of the CFHTLS Deep Survey enable us to obtain sufficiently accurate two-point angular correlation functions to apply a halo occupation distribution analysis. Mean halo masses, calculated as $\langle M_{h} \rangle=10^{11.7}-10^{12.8}h^{-1}M_{\odot}$, increase with stellar-mass limit of LBGs. The threshold halo mass to have a central galaxy follows the same increasing trend with the low-$z$ results, whereas the threshold halo mass to have a satellite galaxy shows higher values at $z=3-5$ than $z=0.5-1.5$ over the entire stellar mass range. Satellite fractions of dropout galaxies, even at less massive haloes, are found to drop sharply from $z=2$ down to less than $0.04$ at $z=3-5$. These results suggest that satellite galaxies form inefficiently within dark haloes at $z=3-5$ even for less massive satellites with $M_{\star}<10^{10}M_{\odot}$. We compute stellar-to-halo mass ratios (SHMRs) assuming a main sequence of galaxies, which is found to provide consistent SHMRs with those derived from a spectral energy distribution fitting method. The observed SHMRs are in good agreement with the model predictions based on the abundance-matching method within $1\sigma$ confidence intervals. We derive observationally, for the first time, $M_{{\rm h}}^{{\rm pivot}}$, which is the halo mass at a peak in the star-formation efficiency, at $3<z<5$, and it shows a little increasing trend with cosmic time at $z>3$. In addition, $M_{{\rm h}}^{{\rm pivot}}$ and its normalization are found to be almost unchanged during $0<z<5$. Our study shows an observational evidence that galaxy formation is ubiquitously most efficient near a halo mass of $M_{{\rm h}}\sim10^{12}M_{\odot}$ over cosmic time.

We present mid-infrared (MIR, 7.5-13.5 $\mu$m) imaging and spectroscopy observations obtained with the CanariCam (CC) instrument on the 10.4m Gran Telescopio CANARIAS for a sample of 20 nearby, MIR bright and X-ray luminous QSOs. We find that for the majority of QSOs the MIR emission is unresolved at angular scales nearly 0.3 arcsec, corresponding to physical scales $<600$ pc. We find that the higher-spatial resolution CC spectra have similar shapes to those obtained with Spitzer/IRS, and hence we can assume that the spectra are not heavily contaminated by extended emission in the host galaxy. We thus take advantage of the higher signal to noise Spitzer/IRS spectra, as a fair representation of the nuclear emission, to decompose it into a combination of active galactic nuclei (AGN), polycyclic aromatic hydrocarbon (PAH) and stellar components. In most cases the AGN is the dominant component, with a median contribution of 85 per cent of the continuum light at MIR (5-15 $\mu$m) within the IRS slit. This IR AGN emission is well reproduced by clumpy torus models. We find evidence for significant differences in the parameters that describe the dusty tori of QSOs when compared with the same parameters of Seyfert 1 and 2 nuclei. In particular, we find a lower number of clouds ($N_{0}<12$), steeper radial distribution of clouds ($q=1.5-3.0$), and clouds that are less optically thick ($\tau_{V}<100$) than in Seyfert 1, which could be attributed to dusty structures that have been partially evaporated and piled up by the higher radiation field in QSOs. We find that the combination of the angular width $\sigma_{torus}$, viewing angle $i$, and number of clouds along the equatorial line $N_{0}$, produces large escape probabilities ($P_{esc} > 2$ per cent) and low geometrical covering factors ($f_{2}<0.6$), as expected for AGN with broad lines in their optical spectra.

We systematically investigate the near- (NIR) to far-infrared (FIR) photometric properties of a nearly complete sample of local active galactic nuclei (AGN) detected in the Swift/Burst Alert Telescope (BAT) all-sky ultra hard X-ray (14-195 keV) survey. Out of 606 non-blazar AGN in the Swift/BAT 70-month catalog at high galactic latitude of $|b|>10^{\circ}$, we obtain IR photometric data of 604 objects by cross-matching the AGN positions with catalogs from the WISE, AKARI, IRAS, and Herschel infrared observatories. We find a good correlation between the ultra-hard X-ray and mid-IR (MIR) luminosities over five orders of magnitude ($41 < \log (L_{14-195}/{\rm erg}~{\rm s}^{-1})< 46$). Informed by previous measures of the intrinsic spectral energy distribution of AGN, we find FIR pure-AGN candidates whose FIR emission is thought to be AGN-dominated with low starformation activity. We demonstrate that the dust covering factor decreases with the bolometric AGN luminosity, confirming the luminosity-dependent unified scheme. We also show that the completeness of the WISE color-color cut in selecting Swift/BAT AGN increases strongly with 14-195 keV luminosity.

We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $\mu$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $\mu$m, 10.3 $\mu$m, and 11.6 $\mu$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $\mu$m total and polarized nuclear spectral energy distribution (SED). The total flux SED in the central 400 pc is explained as the combination of 1) a hot (731 $\pm$ 4 K) dusty structure, directly irradiated by the central engine, which is at 1.6 $\pm$ 0.1 pc away and attributed to be in the pc-scale polar region, 2) an optically-thick, smooth and disk-like dusty structure (`torus') with an inclination of 48 $\pm$ 23$^{\circ}$ surrounding the central engine, and 3) an extinguished (A$_{\mbox{V}} =$ 36 $\pm$ 5 mag) starburst component. The polarized SED decreases from 0.77 $\pm$ 0.14 per cent at 1.2 $\mu$m to 0.31 $\pm$ 0.15 per cent at 11.6 $\mu$m and follows a power-law function, $\lambda^{\sim0.57}$. The polarization angle remains constant ($\sim$108$^{\circ}$) in the 1$-$12 $\mu$m wavelength range. The dominant polarization mechanism is explained as scattering off hot dust grains in the pc-scale polar regions.

We investigate the evolutionary connection between local IR-bright galaxies ($\log L_{\rm IR}\ge 11.4\,L_\odot$) and quasars. We use high angular resolution ($\sim$ 0.3-0.4 arcsec $\sim$ few hundred parsecs) $8-13\,\mu$m ground-based spectroscopy to disentangle the AGN mid-IR properties from those of star formation. The comparison between the nuclear $11.3\,\mu$m PAH feature emission and that measured with Spitzer/IRS indicates that the star formation is extended over a few kpc in the IR-bright galaxies. The AGN contribution to the total IR luminosity of IR-bright galaxies is lower than in quasars. Although the dust distribution is predicted to change as IR-bright galaxies evolve to IR-bright quasars and then to optical quasars, we show that the AGN mid-IR emission of all the quasars in our sample is not significantly different. In contrast, the nuclear emission of IR-bright galaxies with low AGN contributions appears more heavily embedded in dust although there is no clear trend with the interaction stage or projected nuclear separation. This suggests that the changes in the distribution of the nuclear obscuring material may be taking place rapidly and at different interaction stages washing out the evidence of an evolutionary path. When compared to normal AGN, the nuclear star formation activity of quasars appears to be dimming whereas it is enhanced in some IR-bright nuclei, suggesting that the latter are in an earlier star-formation dominated phase.

We present subarcsecond resolution mid-infrared (MIR) images obtained with 8-10 m-class ground-based telescopes of a complete volume-limited (DL<40 Mpc) sample of 24 Seyfert galaxies selected from the Swift/BAT nine month catalog. We use those MIR images to study the nuclear and circumnuclear emission of the galaxies. Using different methods to classify the MIR morphologies on scales of ~400 pc, we find that the majority of the galaxies (75-83%) are extended or possibly extended and 17-25% are point-like. This extended emission is compact and it has low surface brightness compared with the nuclear emission, and it represents, on average, ~30% of the total MIR emission of the galaxies in the sample. We find that the galaxies whose circumnuclear MIR emission is dominated by star formation show more extended emission (650+-700 pc) than AGN-dominated systems (300+-100 pc). In general, the galaxies with point-like MIR morphologies are face-on or moderately inclined (b/a~0.4-1.0), and we do not find significant differences between the morphologies of Sy1 and Sy2. We used the nuclear and circumnuclear fluxes to investigate their correlation with different AGN and SF activity indicators. We find that the nuclear MIR emission (the inner ~70 pc) is strongly correlated with the X-ray emission (the harder the X-rays the better the correlation) and with the [O IV] lambda 25.89 micron emission line, indicating that it is AGN-dominated. We find the same results, although with more scatter, for the circumnuclear emission, which indicates that the AGN dominates the MIR emission in the inner ~400 pc of the galaxies, with some contribution from star formation.

The dwarf spheroidal galaxies (dSphs) in the Milky Way are the primary targets in the indirect searches for particle dark matter. To set robust constraints on candidate of dark matter particles, understanding the dark halo structure of these systems is of substantial importance. In this paper, we first evaluate the astrophysical factors for dark matter annihilation and decay for 24 dSphs, with taking into account a non-spherical dark halo, using generalized axisymmetric mass models based on axisymmetric Jeans equations. First, from a fitting analysis of the most recent kinematic data available, our axisymmetric mass models are a much better fit than previous spherical ones, thus, our work should be the most realistic and reliable estimator for astrophysical factors. Secondly, we find that among analysed dSphs, the ultra-faint dwarf galaxies Triangulum II and Ursa Major II are the most promising but large uncertain targets for dark matter annihilation while the classical dSph Draco is the most robust and detectable target for dark matter decay. It is also found that the non-sphericity of luminous and dark components influences the estimate of astrophysical factors, even though these factors largely depend on the sample size, the prior range of parameters and the spatial extent of the dark halo. Moreover, owing to these effects, the constraints on the dark matter annihilation cross-section are more conservative than those of previous spherical works. These results are important for optimizing and designing dark matter searches in current and future multi-messenger observations by space and ground-based telescopes.

We present 31.5 micron imaging photometry of 11 nearby Seyfert galaxies observed from the Stratospheric Observatory For Infrared Astronomy (SOFIA) using the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST). We tentatively detect extended 31 micron emission for the first time in our sample. In combination with this new data set, subarcsecond resolution 1-18 micron imaging and 7.5-13 micron spectroscopic observations were used to compute the nuclear spectral energy distribution (SED) of each galaxy. We found that the turnover of the torus emission does not occur at wavelengths <31.5 micron, which we interpret as a lower-limit for the wavelength of peak emission. We used CLUMPY torus models to fit the nuclear infrared (IR) SED and infer trends in the physical parameters of the AGN torus for the galaxies in the sample. Including the 31.5 micron nuclear flux in the SED 1) reduces the number of clumpy torus models compatible with the data, and 2) modifies the model output for the outer radial extent of the torus for 10 of the 11 objects. Specifically, six (60%) objects show a decrease in radial extent while four (40%) show an increase. We find torus outer radii ranging from <1pc to 8.4 pc

We present far-infrared (FIR) $70-500\,\mu$m imaging observations obtained with Herschel/PACS and SPIRE of 33 nearby (median distance of 30 Mpc) Seyfert galaxies from the Revised Shapley-Ames (RSA) catalogue. We obtain the FIR nuclear ($r=1\,$kpc and $r=2\,$kpc) and integrated spectral energy distributions (SEDs). We estimate the unresolved nuclear emission at 70 $\mu$m and we fit the nuclear and integrated FIR SEDs with a grey body model. We find that the integrated FIR emission of the RSA Seyferts in our sample is dominated by emission from the host galaxy, with dust properties similar to those of normal galaxies (non AGN). We use four criteria to select galaxies whose nuclear $70\,\mu$m emission has a significant AGN contribution: (1) elevated 70/160 $\mu$m flux ratios, (2)spatially resolved, high dust temperature gradient, (3) $70\,\mu$m excess emission with respect to the fit of the FIR SEDs with a grey body, and (4) excess of nuclear SFR obtained from $70\,\mu$m over SFR from mid-infrared indicators. 16 galaxies (48 per cent of the initial sample) satisfy at least one of these conditions, whereas 10 satisfy half or more. After careful examination of these, we select six bona fide candidates (18 per cent of the initial sample) and estimate that $\sim 40-70$ per cent of their nuclear ($r=1-2\,$kpc) $70\,\mu$m emission is contributed by dust heated by the AGN.

We study the correlation between the [O III]$\lambda 5007$ and X-ray luminosities of local Active Galactic Nuclei (AGNs), using a complete, hard X-ray ($>10$ keV) selected sample in the Swift/BAT 9-month catalog. From our optical spectroscopic observations at the South African Astronomical Observatory and the literature, a catalog of [O III]$\lambda 5007$ line flux for all 103 AGNs at Galactic latitudes of $|b|>15^\circ$ is complied. Significant correlations with intrinsic X-ray luminosity ($L_{\rm X}$) are found both for observed ($L_{\rm [O~III]}$) and extinction-corrected ($L_{\rm [O~III]}^{\rm cor}$) luminosities, separately for X-ray unabsorbed and absorbed AGNs. We obtain the regression form of $L_{\rm [O~III]}$ $\propto L_{\rm 2-10\; keV}^{1.18\pm0.07}$ and $L_{\rm [O~III]}^{\rm cor}$ $\propto L_{\rm 2-10\; keV}^{1.16\pm0.09}$ from the whole sample. The absorbed AGNs with low ($<$0.5\%) scattering fractions in soft X-rays show on average smaller $L_{\rm [O~III]}/L_{\rm X}$ and $L_{\rm [O~III]}^{\rm cor}/L_{\rm X}$ ratios than the other absorbed AGNs, while those in edge-on host galaxies do not. These results suggest that a significant fraction of this population are buried in tori with small opening angles. By using these $L_{\rm [O~III]}$ vs. $L_{\rm X}$ correlations, the X-ray luminosity function of local AGNs (including Compton thick AGNs) in a standard population synthesis model gives much better agreement with the [O III]$\lambda 5007$ luminosity function derived from the Sloan Digital Sky Survey than previously reported. This confirms that hard X-ray observations are a very powerful tool to find AGNs with high completeness.

We investigate optical, infrared, and radio active galactic nucleus (AGN) signs in the merger remnant Arp 187, which hosts luminous jets launched in the order of $10^5$ yr ago but whose present-day AGN activity is still unknown. We find AGN signs from the optical BPT diagram and infrared [OIV]25.89 $\mu$m line, originating from the narrow line regions of AGN. On the other hand, Spitzer/IRS show the host galaxy dominated spectra, suggesting that the thermal emission from the AGN torus is considerably small or already diminished. Combining the black hole mass, the upper limit of radio luminosity of the core, and the fundamental plane of the black hole enable us to estimate X-ray luminosity, which gives $<10^{40}$ erg s$^{-1}$. Those results suggest that the AGN activity of Arp 187 has already been quenched, but the narrow line region is still alive owing to the time delay of emission from the past AGN activity.

We investigate the relationship between X-ray and optical line emission in 340 nearby AGN selected above 10 keV using Swift BAT. We find a weak correlation between the extinction corrected [O III] and hard X-ray luminosity (14-195 keV) with a [OIII] large scatter (R_Pear = 0.64, sigma = 0.62 dex) and a similarly large scatter with the intrinsic 2-10 keV to [O III] luminosities (RPear=0.63, sigma = 0.63 dex). Correlations of the hard X-ray fluxes with the fluxes of high-ionization narrow lines ([O III], He II, [Ne III] and [Ne V]) are not significantly better than with the low ionization lines (Halpha, [SII]). Factors like obscuration or physical slit size are not found to be a significant part of the large scatter. In contrast, the optical emission lines show much better correlations with each other (sigma = 0.3 dex) than with the X-ray flux. The inherent large scatter questions the common usage of narrow emission lines as AGN bolometric luminosity indicators and suggests that other issues such as geometrical differences in the scattering of the ionized gas or long term AGN variability are important.

We present results from the fitting of infrared (IR) spectral energy distributions of 21 active galactic nuclei (AGN) with clumpy torus models. We compiled high spatial resolution ($\sim 0.3$--$0.7$ arcsec) mid-IR $N$-band spectroscopy, $Q$-band imaging and nuclear near- and mid-IR photometry from the literature. Combining these nuclear near- and mid-IR observations, far-IR photometry and clumpy torus models, enables us to put constraints on the torus properties and geometry. We divide the sample into three types according to the broad line region (BLR) properties; type-1s, type-2s with scattered or hidden broad line region (HBLR) previously observed, and type-2s without any published HBLR signature (NHBLR). Comparing the torus model parameters gives us the first quantitative torus geometrical view for each subgroup. We find that NHBLR AGN have smaller torus opening angles and larger covering factors than those of HBLR AGN. This suggests that the chance to observe scattered (polarized) flux from the BLR in NHBLR could be reduced by the dual effects of (a) less scattering medium due to the reduced scattering volume given the small torus opening angle and (b) the increased torus obscuration between the observer and the scattering region. These effects give a reasonable explanation for the lack of observed HBLR in some type-2 AGN.

We discuss the indirect detection of the wino dark matter utilizing gamma-ray observations of dwarf spheroidal galaxies (dSphs). After carefully reviewing current limits with particular attention to astrophysical uncertainties, we show prospects of the wino mass limit in future gamma-ray observation by the Fermi-LAT and the GAMMA-400 telescopes. We find that the improvement of the so-called $J$-factor of both the classical and the ultra-faint dSphs will play a crucial role to cover whole mass range of the wino dark matter. For example, with $\delta (\log_{10}J) = 0.1$ for both the classical and the ultra-faint dSphs, whole wino dark matter mass range can be covered by 15 years and 10 years data at the Fermi-LAT and GAMMA-400 telescopes, respectively.

We present the result of a systematic infrared 2.5-5 um spectroscopic study of 22 nearby infrared galaxies over a wide infrared luminosity range (10 < log(L_IR / Lsun) < 13) obtained from AKARI Infrared Camera (IRC). The unique band of the AKARI IRC spectroscopy enables us to access both the 3.3 um polycyclic aromatic hydrocarbon (PAH) emission feature from star forming activity and the continuum of torus-dust emission heated by an active galactic nucleus (AGN). Applying our AGN diagnostics to the AKARI spectra, we discover 14 buried AGNs. The large fraction of buried AGNs suggests that AGN activity behind the dust is almost ubiquitous in ultra-/luminous infrared galaxies (U/LIRGs). We also find that both the fraction and energy contribution of buried AGNs increase with infrared luminosity from 10 < log(L_IR / Lsun) < 13, including normal infrared galaxies with log (L_IR / Lsun) < 11. The energy contribution from AGNs in the total infrared luminosity is only ~7% in LIRGs and ~20% in ULIRGs, suggesting that the majority of the infrared luminosity originates from starburst activity. Using the PAH emission, we investigate the luminosity relation between star formation and AGN. We find that these infrared galaxies exhibit higher star formation rates than optically selected Seyfert galaxies with the same AGN luminosities, implying that infrared galaxies could be an early evolutionary phase of AGN.

We explore the relationships between the 3.3 {\mu}m polycyclic aromatic hydrocarbon (PAH) feature and active galactic nucleus (AGN) properties of a sample of 54 hard X-ray selected bright AGNs, including both Seyfert 1 and Seyfert 2 type objects, using the InfraRed Camera (IRC) on board the infrared astronomical satellite AKARI. The sample is selected from the 9-month Swift/BAT survey in the 14-195 keV band and all of them have measured X-ray spectra at $E \lesssim 10$ keV. These X-ray spectra provide measurements of the neutral hydrogen column density ($N_{\rm H}$) towards the AGNs. We use the 3.3 {\mu}m PAH luminosity ($L_{\rm 3.3{\mu}m}$) as a proxy for star formation activity and hard X-ray luminosity ($L_{\rm 14-195keV}$) as an indicator of the AGN activity. We search for possible difference of star-formation activity between type 1 (un-absorbed) and type 2 (absorbed) AGNs. We have made several statistical analyses taking the upper-limits of the PAH lines into account utilizing survival analysis methods. The results of our $\log(L_{\rm 14-195keV})$ versus $\log(L_{\rm 3.3{\mu}m})$ regression shows a positive correlation and the slope for the type 1/unobscured AGNs is steeper than that of type 2/obscured AGNs at a $3\sigma$ level. Also our analysis show that the circum-nuclear star-formation is more enhanced in type 2/absorbed AGNs than type 1/un-absorbed AGNs for low X-ray luminosity/low Eddington ratio AGNs, while there is no significant dependence of star-formation activities on the AGN type in the high X-ray luminosities/Eddington ratios.

We present high-angular (~0.4") resolution mid-infrared (MIR) polarimetric observations in the 8.7 ${\mu}$m and 11.6 ${\mu}$m filters of Cygnus A using CanariCam on the 10.4-m Gran Telescopio CANARIAS. A highly polarized nucleus is observed with a degree of polarization of 11${\pm}$3% and 12${\pm}$3% and position angle of polarization of 27${\pm}$8 degrees and 35${\pm}$8 degrees in a 0.38" (~380 pc) aperture for each filter. The observed rising of the polarized flux density with increasing wavelength is consistent with synchrotron radiation from the pc-scale jet close to the core of Cygnus A. Based on our polarization model, the synchrotron emission from the pc-scale jet is estimated to be 14% and 17% of the total flux density in the 8.7 ${\mu}$m and 11.6 ${\mu}$m filters, respectively. A blackbody component with a characteristic temperature of 220 K accounts for >75% of the observed MIR total flux density. The blackbody emission arises from a combination of (1) dust emission in the torus; and (2) diffuse dust emission around the nuclear region, but the contributions of the two components cannot be well constrained in these observations.

The narrow Fe K$\alpha$ line is one of the main signatures of the reprocessing of X-ray radiation from the material surrounding supermassive black holes, and it has been found to be omnipresent in the X-ray spectra of active galactic nuclei (AGN). In this work we study the characteristics of the narrow Fe K$\alpha$ line in different types of AGN. Using the results of a large Suzaku study we find that Seyfert 2s have on average lower Fe K$\alpha$ luminosities than Seyfert 1s for the same 10-50 keV continuum luminosity. Simulating dummy Sy1s and Sy2s populations using physical torus models of X-ray reflected emission, we find that this difference can be explained by means of different average inclination angles with respect to the torus, as predicted by the unified model. Alternative explanations include differences in the intensities of Compton humps or in the photon index distributions. We show that the ratio between the flux of the broad and narrow Fe K$\alpha$ line in the 6.35-6.45 keV range depends on the torus geometry considered, and is on average $<25\%$ and $<15\%$ for type I and type II AGN, respectively. We find evidence of absorption of the narrow Fe K$\alpha$ line flux in Compton-thick AGN, which suggests that part of the reflecting material is obscured. We estimate that on average in obscured AGN the reflected radiation from neutral material is seen through a column density which is 1/4 of that absorbing the primary X-ray emission. This should be taken into account in synthesis models of the CXB and when studying the luminosity function of heavily obscured AGN. We detect the first evidence of the X-ray Baldwin effect in Seyfert 2s, with the same slope as that found for Seyfert 1s, which suggests that the mechanism responsible for the decrease of the equivalent width with the continuum luminosity is the same in the two classes of objects.

The narrow component of the iron K$\alpha$ is an almost ubiquitous feature in the X-ray spectra of active galactic nuclei (AGN) and is believed to originate in neutral material, possibly located in the molecular torus. This would imply a tight connection between the Fe K$\alpha$ equivalent width (EW) and the physical properties of the torus. In a recent work we have shown that the decrease of the covering factor of the torus with the luminosity, as expected by luminosity-dependent unification models, would be able to explain the decrease of Fe K$\alpha$ EW with the luminosity (i.e., the X-ray Baldwin effect). Recent developments in the study of the mid-IR (MIR) spectrum of AGN allow important parameters of the torus to be deduced, such as its covering factor ($f_{\rm\,obs}$) and equatorial column density ($N_{\rm\,H}^{\rm\,T}$), by applying clumpy torus models. Using XMM-Newton/EPIC observations of a sample of 24 type-I AGN, we investigate the relation between the physical parameters of the torus obtained by recent MIR works and the properties of the Fe K$\alpha$ line. We correct the values of the Fe K$\alpha$ EW by taking the inclination angle, the photon index, the equatorial column density, and half-opening angle of the torus into account using a physical torus model of X-ray reprocessed radiation. We find that the relation between Fe K$\alpha$ EW and $f_{\rm\,obs}$ shows a slope that is consistent with the expected value, albeit with a low statistical significance. A trend that is consistent with the theoretical prediction is also found when comparing the Fe K$\alpha$ EW to $N_{\rm\,H}^{\rm\,T}$. Our work seems to confirm that the bulk of the narrow Fe K$\alpha$ line is produced by the same material responsible for the MIR emission.

The large value of the tensor-to-scalar ratio in the cosmic microwave background radiation reported by the BICEP2 collaboration gives strong impact on models of supersymmetry (SUSY). The large ratio indicates inflation with a high-energy scale and thus a high reheating temperature in general, and various SUSY models suffer from the serious gravitino and Polonyi problems. In this article, we discuss a class of the high-scale SUSY breaking models which are completely free from those problems. With especially focusing on the dark matter relic abundance, we examine how the BICEP2 result narrows down the parameter space of the models, assuming the simplest chaotic inflation model. We find that the mass of the dark matter is predicted to be less than about 1 TeV thanks to the non-thermal production in the early universe through the decay of abundant gravitinos produced after the reheating process. We also discuss implications in some details to dark matter searches at collider and indirect dark matter detection experiments.

X-ray surveys have revealed a new class of active galactic nuclei (AGN) with a very low observed fraction of scattered soft X-rays, f_scat < 0.5%. Based on X-ray modeling these "X-ray new-type", or low observed X-ray scattering (hereafter:"low-scattering") sources have been interpreted as deeply-buried AGN with a high covering factor of gas. In this paper we address the questions whether the host galaxies of low-scattering AGN may contribute to the observed X-ray properties, and whether we can find any direct evidence for high covering factors from the infrared (IR) emission. We find that X-ray low-scattering AGN are preferentially hosted by highly-inclined galaxies or merger systems as compared to other Seyfert galaxies, increasing the likelihood that the line-of-sight toward the AGN intersects with high columns of host-galactic gas and dust. Moreover, while a detailed analysis of the IR emission of low-scattering AGN ESO 103-G35 remains inconclusive, we do not find any indication of systematically higher dust covering factors in a sample of low-scattering AGN based on their IR emission. For ESO 103-G35, we constrained the temperature, mass and location of the IR emitting dust which is consistent with expectations for the dusty torus. However, a deep silicate absorption feature probably from much cooler dust suggests an additional screen absorber on larger scales within the host galaxy. Taking these findings together, we propose that the low f_scat observed in low-scattering AGN is not necessarily the result of circumnuclear dust but could originate from interference of host-galactic gas with a column density of the order of 10^22 cm^-2 with the line-of-sight. We discuss implications of this hypothesis for X-ray models, high-ionization emission lines, and observed star-formation activity in these objects.

We carried out optical high-dispersion spectroscopic monitoring of the Be disk in a Be/X-ray binary A 0535+262/V725 Tau from 2009 to 2012, covering two giant outbursts and several normal outbursts. This monitoring was performed in order to investigate variabilities of the Be disk due to the interaction with the neutron star in recent X-ray active phase from 2008 to 2011. Such variabilities give a clue to uncleared detailed mechanism for very bright X-ray outbursts, which are unique to some Be/X-ray binaries with relatively wide and eccentric orbit. In the previous letter (Moritani et al. 2011), a brief overview of line profile variabilities around the 2009 giant outburst was given and the possibility of a warped Be disk was discussed. In this paper, a full analysis of the Halpha line profiles as well as other line profiles is carried out. A bright blue component, or blue "shoulder", showing up after periastron indicates the presence of a dense gas stream toward the neutron star, which is associated with observed outbursts. We re-analyze the Halpha line profiles before 2009 (down to 2005) in order to investigate the variability of the the disk structure in the innermost region, which seems to have detached from the Be star surface by 2008. A redshifted enhanced component is remarkable in all emission lines observed around the 2009 giant outburst, occasionally forming a triple peak. These features indicate that the Be disk was warped in X-ray active phase. We estimate the position of the warped region from fitting the radial velocity of the redshifted enhanced component of Halpha, and find that it was very close to the periastron when two giant outbursts in 2009 and 2011 and a bright normal outburst in 2010 March occurred. These facts strongly suggest that the warped Be disk triggered these giant outbursts.

We present results of four-pointing Suzaku X-ray observations (total ~200 ks) of the intracluster medium (ICM) in the Abell 1835 galaxy cluster (kT ~ 8 keV, z = 0.253) out to the virial radius (r_vir ~ 2.9 Mpc) and beyond. Faint X-ray emission from the ICM out to r_vir is detected. The temperature gradually decreases with radius from ~8 keV in the inner region to ~2 keV at r_vir. The entropy profile is shown to flatten beyond r_500, in disagreement with the r_1.1 dependence predicted from the accretion shock heating model. The thermal pressure profile in the range 0.3r_500 < r < r_vir agrees well with that obtained from the stacked Sunyaev-Zel'dovich effect observations with the Planck satellite. The hydrostatic mass profile in the cluster outskirts (r_500 < r < r_vir) falls well short of the weak lensing one derived from Subaru/Suprime-Cam observations, showing an unphysical decrease with radius. The gas mass fraction at r_vir defined with the lensing total mass agrees with the cosmic baryon fraction from the WMAP 7-year data. All these results indicate, rather than the gas-clumping effect, that the bulk of the ICM in the cluster outskirts is far from hydrostatic equilibrium and infalling matter retained some of its kinetic energy. Finally, combining with our recent Suzaku and lensing analysis of Abell 1689, a cluster of similar mass, temperature, and redshift, we show that the cluster temperature distribution in the outskirts is significantly correlated with the galaxy density field in the surrounding large-scale environment at (1-2)r_vir.

We present results of the 2.5-5 {\mu}m spectroscopy of a sample of hard X-ray selected active galactic nuclei (AGNs) using the grism mode of the InfraRed Camera (IRC) on board the infrared astronomical satellite AKARI. The sample is selected from the 9-month Swift/BAT survey in the 14-195 keV band, which provides a fair sample of AGNs including highly absorbed ones. The 2.5-5 {\mu}m spectroscopy provide a strong diagnostic tool for the circumnuclear environment of AGNs through the continuum shapes and emission/absorption features such as the 3.3 {\mu}m polycyclic aromatic hydrocarbon (PAH) emission and the broad 3.1 {\mu}m H2O ice, 3.4 {\mu}m bare carbonaceous dust, 4.26 {\mu}m CO2 and 4.67 {\mu}m CO absorptions. As our first step, we use the 3.3 {\mu}m PAH emission as a proxy for the star-formation activity and searched for possible difference of star-formation activity between type 1 (unabsorbed) and type 2 (absorbed) AGNs. We found no significant dependence of the 3.3 {\mu}m PAH line luminosity, normalized by the black hole mass, on optical AGNs types or the X-ray measured column densities.

We report Suzaku observations of the northern half of the Hydra A cluster out to ~1.4 Mpc, reaching the virial radius. This is the first Suzaku observations of a medium-size (kT ~3 keV) cluster out to the virial radius. Two observations were conducted, north-west and north-east offsets, which continue in a filament direction and a void direction of the large-scale structure of the Universe, respectively. The X-ray emission and distribution of galaxies elongate in the filament direction. The temperature profiles in the two directions are mostly consistent with each other within the error bars and drop to 1.5 keV at 1.5 r_500. As observed by Suzaku in hot clusters, the entropy profile becomes flatter beyond r_500, in disagreement with the r^1.1 relationship that is expected from accretion shock heating models. When scaled with the average intracluster medium (ICM) temperature, the entropy profiles of clusters observed with Suzaku are universal and do not depend on system mass. The hydrostatic mass values in the void and filament directions are in good agreement, and the Navarro, Frenk, and White universal mass profile represents the hydrostatic mass distribution up to ~ 2 r_500. Beyond r_500, the ratio of gas mass to hydrostatic mass exceeds the result of the Wilkinson microwave anisotropy probe, and at r_100, these ratios in the filament and void directions reach 0.4 and 0.3, respectively. We discuss possible deviations from hydrostatic equilibrium at cluster outskirts. We derived radial profiles of the gasmass- to-light ratio and iron-mass-to-light ratio out to the virial radius. Within r_500, the iron-mass-to-light ratio of the Hydra A cluster was compared with those in other clusters observed with Suzaku.

We present initial results from the CO survey toward high redshift galaxies using the Nobeyama 45m telescope. Using the new wide bandwidth spectrometer equipped with a two-beam SIS receiver, we have robust new detections of three high redshift (z=1.6-3.4) submillimeter galaxies (SXDF 1100.001, SDP9, and SDP17), one tentative detection (SDSS J160705+533558), and one non-detection (COSMOS-AzTEC1). The galaxies observed during the commissioning phase are sources with known spectroscopic redshifts from previous optical or from wide-band submm spectroscopy. The derived molecular gas mass and line widths from Gaussian fits are ~10^11 Msun and 430-530 km/s, which are consistent with previous CO observations of distant submm galaxies and quasars. The spectrometer that allows a maximum of 32 GHz instantaneous bandwidth will provide new science capabilities at the Nobeyama 45m telescope, allowing us to determine redshifts of bright submm selected galaxies without any prior redshift information.

We investigate the mid- (MIR) to far-infrared (FIR) properties of a nearly complete sample of local Active Galactic Nuclei (AGNs) detected in the Swift/BAT all sky hard X-ray (14-195 keV) survey, based on the cross correlation with the AKARI infrared survey catalogs complemented by those with IRAS and WISE. Out of 135 non-blazer AGNs in the Swift/BAT 9 month catalog, we obtain the MIR photometric data for 128 sources either in the 9, 12, 18, 22, and/or 25 um band. We find good correlation between their hard X-ray and MIR luminosities over 3 orders of magnitude (42< log lambda L_{lambda}(9, 18 um)< 45), which is tighter than that with the FIR luminosities at 90 um. This suggests that thermal emission from hot dusts irradiated by the AGN emission dominate the MIR fluxes. Both X-ray unabsorbed and absorbed AGNs follow the same correlation, implying isotropic infrared emission, as expected in clumpy dust tori rather than homogeneous ones. We find excess signals around 9 um in the averaged infrared spectral energy distribution from heavy obscured "new type" AGNs with small scattering fractions in the X-ray spectra. This could be attributed to the PAH emission feature, suggesting that their host galaxies have strong starburst activities.

We combine data from two all-sky surveys in order to study the connection between the infrared and hard X-ray (>10keV) properties for local active galactic nuclei (AGN). The Swift/Burst Alert Telescope all-sky survey provides an unbiased, flux-limited selection of hard X-ray detected AGN. Cross-correlating the 22-month hard X-ray survey with the AKARI all-sky survey, we studied 158 AGN detected by the AKARI instruments. We find a strong correlation for most AGN between the infrared (9, 18, and 90 micron) and hard X-ray (14-195 keV) luminosities, and quantify the correlation for various subsamples of AGN. Partial correlation analysis confirms the intrinsic correlation after removing the redshift contribution. The correlation for radio galaxies has a slope and normalization identical to that for Seyfert 1s, implying similar hard X-ray/infrared emission processes in both. In contrast, Compton-thick sources show a large deficit in the hard X-ray band, because high gas column densities diminish even their hard X-ray luminosities. We propose two photometric diagnostics for source classification: one is an X-ray luminosity vs. infrared color diagram, in which type 1 radio-loud AGN are well isolated from the others in the sample. The other uses the X-ray vs. infrared color as a useful redshift-independent indicator for identifying Compton-thick AGN. Importantly, Compton-thick AGN and starburst galaxies in composite systems can also be differentiated in this plane based upon their hard X-ray fluxes and dust temperatures. This diagram may be useful as a new indicator to classify objects in new and upcoming surveys such as WISE and NuSTAR.

We report on high-dispersion optical spectroscopic observations of the Be/X-ray binary A0535+262/V725 Tau during the giant outburst in November/December 2009 and after it. The observed emission line profiles, reflecting the structure of the geometrically thin circumstellar envelope of the Be star (Be disk), show drastic variabilities and indicate the existence of a warped component. The enhanced blue shoulder seen after periastron passage implies the gas stream from a dense part of the Be disk to the neutron star.

Objects in the Edgeworth-Kuiper belt and the main asteroid belt should emit microwaves that may give rise to extra anisotropy signals in the multipole of the cosmic microwave background (CMB) experiment. Constraints are derived from the absence of positive detection of such anisotropies for ell < 50, giving the total mass of Edgeworth-Kuiper belt objects to be smaller than 0.2 earth mass. This limit is consistent with the mass extrapolated from the observable population with the size of a > 15 km, assuming that the small-object population follows the power law in size dN/da ~ a^{-q} with the canonical index expected for collisional equilibrium, q ~ 3.5, with which 23% of the mass is ascribed to objects smaller than are observationally accessible down to grains. A similar argument applied to the main asteroid belt indicates that the grain population should not increase faster than q ~ 3.6 towards smaller radii, if it follows the power law continued to observed asteroids with larger radii. It is underlined that both cases are at or only slightly above the limit that can be physically significant, implying the importance of tightening further the CMB anisotropy limit, which may be attained with the observation at higher radio frequencies.

We present the results of a ground-based, high spatial resolution infrared 18 micron imaging study of nearby luminous infrared galaxies (LIRGs), using the Subaru 8.2-m and Gemini South 8.1-m telescopes. The diffraction-limited images routinely achieved with these telescopes in the Q-band (17-23 micron) allow us to investigate the detailed spatial distribution of infrared emission in these LIRGs. We then investigate whether the emission surface brightnesses are modest, as observed in starbursts, or are so high that luminous active galactic nuclei (AGNs; high emission surface brightness energy sources) are indicated. The sample consists of 18 luminous buried AGN candidates and starburst-classified LIRGs identified in earlier infrared spectroscopy. We find that the infrared 18 micron emission from the buried AGN candidates is generally compact, and the estimated emission surface brightnesses are high, sometimes exceeding the maximum value observed in and theoretically predicted for a starburst phenomenon. The starburst-classified LIRGs usually display spatially extended 18 micron emission and the estimated emission surface brightnesses are modest, within the range sustained by a starburst phenomenon. The general agreement between infrared spectroscopic and imaging energy diagnostic methods suggests that both are useful tools for understanding the hidden energy sources of the dusty LIRG population.

We present the results of infrared L-band (3-4 micron) slit spectroscopy of 30 PG QSOs at z < 0.17, the representative sample of local high-luminosity, optically selected AGNs. The 3.3 micron polycyclic aromatic hydrocarbon (PAH) emission feature is used to probe nuclear (< a few kpc) starburst activity and to investigate the connections between AGNs and nuclear starbursts in PG QSOs. The 3.3 micron PAH emission is detected in the individual spectra of 5/30 of the observed PG QSOs. We construct a composite spectrum of PAH-undetected PG QSOs and discern the presence of the 3.3 micron PAH emission therein. We estimate the nuclear-starburst and AGN luminosities from the observed 3.3 micron PAH emission and 3.35 micron continuum luminosities, respectively, and find that the nuclear-starburst-to-AGN luminosity ratios in PG QSOs are similar to those of previously studied AGN populations with lower luminosities, suggesting that AGN-nuclear starburst connections are valid over the wide luminosity range of AGNs in the local universe. The observed nuclear-starburst-to-AGN luminosity ratios in PG QSOs with available supermassive black hole masses are comparable to a theoretical prediction based on the assumption that the growth of a supermassive black hole is controlled by starburst-induced turbulence.

We investigate constraints on the time variation of the fine structure constant between the recombination epoch and the present epoch, \Delta\alpha/\alpha \equiv (\alpha_{rec} - \alpha_{now})/\alpha_{now}, from cosmic microwave background (CMB) taking into account simultaneous variation of other physical constants, namely the electron mass m_{e} and the proton mass m_{p}. In other words, we consider the variation of Yukawa coupling and the QCD scale \Lambda_{QCD} in addition to the electromagnetic coupling. We clarify which parameters can be determined from CMB temperature anisotropy in terms of singular value decomposition. Assuming a relation among variations of coupling constants governed by a single scalar field (the dilaton), the 95% confidence level (C.L.) constraint on \Delta\alpha/\alpha is found to be -8.28 \times 10^{-3} < \Delta\alpha/\alpha < 1.81 \times 10^{-3}, which is tighter than the one obtained by considering only the change of \alpha and m_{e}. We also obtain the constraint on the time variation of the proton-to-electron mass ratio \mu \equiv m_{p}/m_{e} to be -0.52 < \Delta\mu/\mu < 0.17 (95% C.L.) under the same assumption. Finally, we also implement a forecast for constraints from the PLANCK survey.

Non-zero neutrino mass would affect the evolution of the Universe in observable ways, and a strong constraint on the mass can be achieved using combinations of cosmological data sets. We focus on the power spectrum of cosmic microwave background (CMB) anisotropies, the Hubble constant H_0, and the length scale for baryon acoustic oscillations (BAO) to investigate the constraint on the neutrino mass, m_nu. We analyze data from multiple existing CMB studies (WMAP5, ACBAR, CBI, BOOMERANG, and QUAD), recent measurement of H_0 (SHOES), with about two times lower uncertainty (5%) than previous estimates, and recent treatments of BAO from the Sloan Digital Sky Survey (SDSS). We obtained an upper limit of m_nu < 0.2eV (95% C.L.), for a flat LambdaCDM model. This is a 40% reduction in the limit derived from previous H_0 estimates and one-third lower than can be achieved with extant CMB and BAO data. We also analyze the impact of smaller uncertainty on measurements of H_0 as may be anticipated in the near term, in combination with CMB data from the Planck mission, and BAO data from the SDSS/BOSS program. We demonstrate the possibility of a 5 sigma detection for a fiducial neutrino mass of 0.1eV or a 95% upper limit of 0.04eV for a fiducial of m_nu = 0eV. These constraints are about 50% better than those achieved without external constraint. We further investigate the impact on modeling where the dark-energy equation of state is constant but not necessarily -1, or where a non-flat universe is allowed. In these cases, the next-generation accuracies of Planck, BOSS, and 1% measurement of H_0 would all be required to obtain the limit m_nu < 0.05 - 0.06eV (95% C.L.) for the fiducial of m_nu = 0eV. The independence of systematics argues for pursuit of both BAO and H_0 measurements.

We investigate the possibilities of constraining the light gravitino mass m_{3/2} from future cosmic microwave background (CMB) surveys. A model with light gravitino with the mass m_{3/2}<O(10) eV is of great interest since it is free from the cosmological gravitino problem and, in addition, can be compatible with many baryogenesis/leptogenesis scenarios such as the thermal leptogenesis. We show that the lensing of CMB anisotropies can be a good probe for m_{3/2} and obtain an expected constraint on m_{3/2} from precise measurements of lensing potential in the future CMB surveys, such as the PolarBeaR and CMBpol experiments. If the gravitino mass is m_{3/2}=1 eV, we will obtain the constraint for the gravitino mass as m_{3/2} < 3.2 eV (95% C.L.) for the case with Planck+PolarBeaR combined and m_{3/2}=1.04^{+0.22}_{-0.26} eV (68% C.L.) for CMBpol. The issue of Bayesian model selection is also discussed.

The 21-cm PDF (i.e., distribution of pixel brightness temperatures) is expected to be highly non-Gaussian during reionization and to provide important information on the distribution of density and ionization. We measure the 21-cm PDF as a function of redshift in a large simulation of cosmic reionization and propose a simple empirical fit. Guided by the simulated PDF, we then carry out a maximum likelihood analysis of the ability of upcoming experiments to measure the shape of the 21-cm PDF and derive from it the cosmic reionization history. Under the strongest assumptions, we find that upcoming experiments can measure the reionization history in the mid to late stages of reionization to 1-10% accuracy. Under a more flexible approach that allows for four free parameters at each redshift, a similar accuracy requires the lower noise levels of second-generation 21-cm experiments.

In this paper, we put constraints on neutrino properties such as mass $m_{\nu}$ and degeneracy parameters $\xi_i$ from WMAP5 data and light element abundances by using a Markov chain Monte Carlo (MCMC) approach. In order to take consistently into account the effects of the degeneracy parameters, we run the Big Bang Nucleosynthesis code for each value of $\xi_i$ and the other cosmological parameters to estimate the Helium abundance, which is then used to calculate CMB anisotropy spectra instead of treating it as a free parameter. We find that the constraint on $m_{\nu}$ is fairly robust and does not vary very much even if the lepton asymmetry is allowed, and is given by $\sum m_\nu < 1.3 \rm eV$ ($95 % \rm C.L.$).

We discuss how much we can probe the effective number of neutrino species N_nu with cosmic microwave background alone. Using the data of WMAP, ACBAR, CBI and BOOMERANG experiments, we obtain a constraint on the effective number of neutrino species as 0.96< N_nu <7.94 at 95% C.L. for a power-law LCDM flat universe model. The limit is improved to be 1.39 < N_nu < 6.38 at 95% C.L. if we assume that the baryon density, N_nu and the helium abundance are related by the big bang nucleosynthesis theory. We also provide a forecast for the PLANCK experiment using a Markov chain Monte Carlo approach. In addition to constraining N_nu, we investigate how the big bang nucleosynthesis relation affects the estimation for these parameters and other cosmological parameters.

We studied a constraint on the primordial helium abundance Y_p from current and future observations of CMB. Using the currently available data from WMAP, ACBAR, CBI and BOOMERANG, we obtained the constraint as Y_p = 0.25^{+0.10}_{-0.07} at 68% C.L. We also provide a forecast for the Planck experiment using the Markov chain Monte Carlo approach. In addition to forecasting the constraint on Y_p, we investigate how assumptions for Y_p affect constraints on the other cosmological parameters.

We study non-Gaussianity, the spectral index of primordial scalar fluctuations and tensor modes in models where fluctuations from the inflaton and the curvaton can both contribute to the present cosmic density fluctuations. Even though simple single-field inflation models generate only tiny non-Gaussianity, if we consider such a mixed scenario, large non-Gaussianity can be produced. Furthermore, we study the inflationary parameters such as the spectral index and the tensor-to-scalar ratio in this kind of models and discuss in what cases models predict the spectral index and tensor modes allowed by the current data while generating large non-Gaussianity, which may have many implications for model-buildings of the inflationary universe.

We investigate non-Gaussianity in the modulated reheating scenario where fluctuations of the decay rate of the inflaton generate adiabatic perturbations, paying particular attention to the non-linearity parameters $f_{\rm NL}, \tau_{\rm NL}$ and $g_{\rm NL}$ as well as the scalar spectral index and tensor-to-scalar ratio which characterize the nature of the primordial power spectrum. We also take into account the pre-existing adiabatic perturbations produced from the inflaton fluctuations. It has been known that the non-linearity between the curvature perturbations and the fluctuations of the decay rate can yield non-Gaussianity at the level of $f_{\rm NL} \sim \mathcal{O}(1)$, but we find that the non-linearity between the decay rate and the modulus field which determines the decay rate can generate much greater non-Gaussianity. We also discuss a consistency relation among non-linearity parameters which holds in the scenario and find that the modulated reheating yields a different one from that of the curvaton model. In particular, they both can yield a large positive $f_{\rm NL}$ but with a different sign of $g_{\rm NL}$. This provides a possibility to discriminate these two competitive models by looking at the sign of $g_{\rm NL}$. Furthermore, we work on some concrete inflation models and investigate in what cases models predict the spectral index and the tensor-to-scalar ratio allowed by the current data while generating large non-Gaussianity, which may have many implications for model-buildings of the inflationary universe.

We study how the determination of the Hubble constant from cosmological distance measures is affected by models of dark energy and vice versa. For this purpose, constraints on the Hubble constant and dark energy are investigated using the cosmological observations of cosmic microwave background, baryon acoustic oscillations and type Ia suprenovae. When one investigates dark energy, the Hubble constant is often a nuisance parameter, thus it is usually marginalized over. On the other hand, when one focuses on the Hubble constant, simple dark energy models such as a cosmological constant and a constant equation of state are usually assumed. Since we do not know the nature of dark energy yet, it is interesting to investigate the Hubble constant assuming some types of dark energy and see to what extent the constraint on the Hubble constant is affected by the assumption concerning dark energy. We show that the constraint on the Hubble constant is not affected much by the assumption for dark energy. We furthermore show that this holds true even if we remove the assumption that the universe is flat. We also discuss how the prior on the Hubble constant affects the constraints on dark energy and/or the curvature of the universe.

We discuss the possibilities of the simultaneous determination of the neutrino masses and the evolution of dark energy from future cosmological observations such as cosmic microwave background (CMB), large scale structure (LSS) and the cross correlation between them. Recently it has been discussed that there is a degeneracy between the neutrino masses and the equation of state for dark energy. It is also known that there are some degeneracies among the parameters describing the dark energy evolutions. We discuss the implications of these on the cross correlation of CMB with LSS in some details. Then we consider to what extent we can determine the neutrino masses and the dark energy evolution using the expected data from CMB, LSS and their cross correlation.

We show that the cosmic microwave background (CMB) data of WMAP can give subelectronvolt limit on the neutrino mass: m_nu < 0.63 eV (95% CL). We also investigate its degeneracy with other cosmological parameters. In particular, we show the Hubble constant derived from the WMAP data decreases considerably when the neutrino mass is a few times 0.1 eV.

We investigate the modifications to predictions for the abundances of light elements from standard Big-Bang nucleosynthesis when exotic late-decaying particles with lifetimes exceeding ~1 sec are prominent in the early Universe. Utilising a model-independent analysis of the properties of these long-lived particles, we identify the parameter space associated with models that are consistent with all observational data and hence resolve the much discussed discrepancies between observations and theoretical predictions for the abundances of Li^7 and Li^6.

We discuss constraints on the effective number of neutrino species Nnu from recent cosmological observations such as CMB, LSS, BBN, including our own analysis which uses the WMAP and the Luminous Red Galaxy power spectrum data. We also discuss their implications on some non-standard cosmological scenarios such as the low (MeV-scale) reheating temperature scenario and the scenario with decaying particles between BBN and structure formation.

We present models of decaying particles for increasing the effective number of neutrinos N_nu after big bang nucleosynthesis but before the structure formation begins. We point out that our scenario not only solves the discrepancy between the constraints on N_nu from these two epochs, but also provides a possible answer to deeper inconsistency in the estimation of the matter power spectrum amplitude at small scales, represented by sigma_8, between the WMAP and some small scale matter power measurements such as the Lyman-alpha forest and weak lensing. We consider (a) saxion decay into two axions; (b) gravitino decay into axino and axion; (c) Dirac right-handed sneutrino decay into gravitino and right-handed neutrino.

We derive constraint on the effective number of neutrino species N_nu from the cosmic microwave background power spectrum of the WMAP and galaxy clustering power spectrum of the SDSS luminous red galaxies (LRGs). Using these two latest data sets of CMB and galaxy clustering alone, we obtain the limit 0.9 < N_nu < 8.2 (95% C.L.) for the power-law LambdaCDM flat universe, with no external prior. The lower limit corresponds to the lower bound on the reheating temperature of the universe T_R > 2 MeV.

We investigate observational constraints on the curvature of the universe not restricting ourselves to a cosmological constant as dark energy, in particular allowing a dark energy equation of state to evolve with time in several ways. We use type Ia supernovae (SNeIa) data from the latest gold data set which includes 182 SNeIa, along with the CMB shift parameter and the baryon acoustic oscillation peak. We show quantitatively that the constraint on the curvature of the universe depends on dark energy model: some popular parametrizations give constraints closely around the flat universe at 5% level (2 sigma C.L.) whereas some parametrizations allow the universe to be as open as Omega_k ~ 0.2.

We investigate cosmological evolution and implications of cosmic strings with time-dependent tension. We derive basic equations of time development of the correlation length and the velocity of such strings, based on the one scale model. Then, we find that, in the case where the tension depends on some power of the cosmic time, cosmic strings with time-dependent tension goes into the scaling solution if the power is lower than a critical value. We also discuss cosmic microwave background anisotropy and matter power spectra produced by these strings. The constraints on their tensions from the Wilkinson microwave anisotropy probe (WMAP) three year data and Sloan digital sky survey (SDSS) data are also given.

We investigate how the nature of dark energy affects the determination of the curvature of the universe from recent observations. For this purpose, we consider the constraints on the matter and dark energy density using observations of type Ia supernovae, baryon acoustic oscillation peak and cosmic microwave background with several types of dark energy equation of state. Although it is usually said that the combination of current observations favors a flat universe, we found that a relatively large parameter space allows the universe to be open for a particular model of dark energy. We also discuss what kind of dark energy model or prior allow a non-flat universe.

We study constraints on time variation of the fine structure constant alpha from cosmic microwave background (CMB) taking into account simultaneous change in alpha and the electron mass m_e which might be implied in unification theories. We obtain the constraints -0.097 < Delta alpha/alpha < 0.034 at 95% C.L. using WMAP data only, and -0.042 < Delta alpha/alpha < 0.026 combining with the constraint on the Hubble parameter by the HST Hubble Key Project. These are improved by 15% compared with constraints assuming only alpha varies. We discuss other relations between variations in alpha and m_e but we do not find evidence for varying alpha.

We derive an upper limit on the neutrino mass from the WMAP three-year data alone by employing a deterministic minimisation method based on a grid search in multi-parameter space. Assuming the flat LCDM model with power-law adiabatic perturbations, we find the sum of m_nu < 2.0 eV in agreement with the result of the WMAP team. This result, the limit being nearly the same as that from the WMAP first-year data, means that the fundamental limit on the neutrino mass obtainable from the cosmic microwave background alone is already nearly met, as anticipated from the previous analysis. We also clarify the role of the polarisation data in deriving the limit on the neutrino mass.

We revisit the constraint on the primordial helium mass fraction Yp from observations of cosmic microwave background (CMB) alone. By minimizing chi square of recent CMB experiments over 6 other cosmological parameters, we obtained rather weak constraints as 0.17 < Yp < 0.52 at 1 sigma C.L. for a particular data set. We also study the future constraint on cosmological parameters when we take account of the prediction of the standard big bang nucleosynthesis (BBN) theory as a prior on the helium mass fraction where Yp can be fixed for a given energy density of baryon. We discuss the implications of the prediction of the standard BBN on the analysis of CMB.

We discuss the constraints on the time-varying equation of state for dark energy and the curvature of the universe using observations of type Ia supernovae from Riess et al. and the most recent Supernova Legacy Survey (SNLS), the baryon acoustic oscillation peak detected in the SDSS luminous red galaxy survey and cosmic microwave background. Due to the degeneracy among the parameters which describe the time dependence of the equation of state and the curvature of the universe, the constraints on them can be weakened when we try to constrain them simultaneously, in particular when we use a single observational data. However, we show that we can obtain relatively severe constraints when we use all data sets from observations above even if we consider the time-varying equation of state and do not assume a flat universe. We also found that the combined data set favors a flat universe even if we consider the time variation of dark energy equation of state.

We study how the oscillations of the neutrinos affect their thermalization process during the reheating period with temperature O(1) MeV in the early universe. We follow the evolution of the neutrino density matrices and investigate how the predictions of big bang nucleosynthesis vary with the reheating temperature. For the reheating temperature of several MeV, we find that including the oscillations makes different predictions, especially for $^4$He abundance. Also, the effects on the lower bound of the reheating temperature from cosmological observations are discussed.

It is shown that a subelectronvolt upper limit can be derived on the neutrino mass from the CMB data alone in the Lambda CDM model with the power-law adiabatic perturbations, without the aid of any other cosmological data. Assuming the flatness of the universe, the constraint we can derive from the current WMAP observations is \sum m_nu < 2.0 eV at the 95% confidence level for the sum over three species of neutrinos (m_nu < 0.66 eV for the degenerate neutrinos). This constraint modifies little even if we abandon the flatness assumption for the spatial curvature. We argue that it would be difficult to improve the limit much beyond \sum m_nu \lesssim 1.5 eV using only the CMB data, even if their statistics are substantially improved. However, a significant improvement of the limit is possible if an external input is introduced that constrains the Hubble constant from below. The parameter correlation and the mechanism of CMB perturbations that give rise to the limit on the neutrino mass are also elucidated.

Within the standard big bang nucleosynthesis (BBN) and cosmic microwave background (CMB) framework, the baryon density measured by the Wilkinson Microwave Anisotropy Probe (WMAP) or the primordial D abundance is much higher than the one measured by the $^4$He or $^7$Li abundances. To solve the discrepancy, we propose a scenario in which additional baryons appear after BBN. We show that simply adding the baryons can not be a solution but the existence of a large lepton asymmetry before BBN makes the scenario successful. These extra baryons and leptons, in addition to the initial baryons which exist before the BBN, can be all produced from $Q$-balls.

We calculate primordial abundances of light elements produced during big bang nucleosynthesis when the fine structure constant and/or the cosmic expansion rate take non-standard values. We compare them with the recent values of observed D, He4 and Li7 abundances, which show slight inconsistency among themselves in the standard big bang nucleosynthesis scenario. This inconsistency is not solved by considering either a varying fine structure constant or a non-standard expansion rate separately but solutions are found by their simultaneous existence.

We reconsider the evolution of strongly degenerate neutrinos in the early universe. Our chief concern is the validity of the entropy conservation after the neutrino annihilation process has frozen out (so that the establishment of chemical equilibrium is not trivial). We argue that the entropy indeed conserves because elastic scattering keeps the neutrino and antineutrino distribution functions in the equilibrium form and the sum of their chemical potential keeps zero even after the neutrino annihilation freeze-out. We also simulate the evolution of the degenerate neutrino spectrum to support the argument. We conclude that the change in the neutrino degeneracy parameter when the relativistic degrees of freedom in the universe decreases is calculated using the entropy conservation and the lepton number conservation without worrying about at what temperature the neutrino annihilation process freezes out.

We re-investigate the evolution of the strongly degenerate neutrinos in the early universe. With the larger degeneracy, the neutrino number freezes at higher temperatures because the neutrino annihilation rate decreases. We consider very large degeneracy so large that the neutrino number freezes before events in which the particle degrees of freedom in the universe decrease (e.g. the muon annihilation and the quark-hadron phase transition). In such a case, the degeneracy by the time of nucleosynthesis becomes smaller than the initial degeneracy. We calculate how much it decreases from the initial value on the basis of the conservation of the neutrino number and the total entropy. We found a large drop in the degeneracy but it is not large enough to affect the current constraints on the neutrino degeneracy from BBN and CMBR.