Supermassive binary black holes (SMBBHs) in galactic centers may radiate gravitational wave (GW) in the nano-Hertz frequency band, which are expected to be detected by pulsar timing arrays (PTAs) in the near future. GW signals from individual SMBBHs at cosmic distances, if detected by PTAs, are potentially powerful standard sirens that can be used to independently measure distances and thus put constraints on cosmological parameters. In this paper, we investigate the constraint that may be obtained on the equation of state ($w$) of dark energy by using those SMBBHs, expected to be detected by the PTAs in the Square Kilometre Array (SKA) era. By considering both the currently available SMBBH candidates and mock SMBBHs in the universe resulting from a simple galaxy major merger model, we find that $\sim 200$ to $3000$ SMBBHs with chirp mass $>10^9M_\odot$ are expected to be detected with signal-to-noise ratio $>10$ by SKA-PTA with conservative and optimistic settings and they can be used to put a constraint on $w$ to an uncertainty of $\Delta w\sim 0.02-0.1$. If further information on the mass and mass ratio of those SMBBHs can be provided by electromagnetic observations (e.g., chirp mass uncertainty $\lesssim 50\%$), the constraint may be further improved to $\lesssim 0.01$ level, as many more SMBBHs will be detected by SKA-PTA with relatively better distance measurements and can be used as the standard sirens.

(abridged) The far-UV wavelength range (912-2000A) provides access to atomic and molecular transitions of many species the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium, within phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way thanks to absorption features appearing in the UV spectra of hot stars and yielding fundamental insights into the composition and physical characteristics of all phases of the ISM along with the processes that influence them. However, we have yet to design a spectrometer able to observe the full UV domain at resolving power R>10^5 with a signal-to-noise ratio SNR>500. Such a resolution is necessary to resolve lines from both the cold molecular hydrogen and the warm metal ions with a turbulent velocity of about 1 km s-1, and to differentiate distinct velocity components. Future UV spectroscopic studies of the Milky Way ISM must revolutionize our understanding of the ISM as a dynamical, unstable, and magnetized medium, and rise to the challenge brought forward by current theories. Another interesting prospect is to transpose the same level of details that has been reached for the Milky Way to the ISM in external galaxies, in particular in metal-poor galaxies, where the ISM chemical composition, physical conditions, and topology change dramatically, with significant consequences on the star-formation properties. Finally, we need to be able to perform statistical analyses of background quasar lines of sight intersecting the CGM of galaxies at various redshifts and to comprehend the role of gas exchanges and flows for galaxy evolution.

We present deep $J$ and $H$-band images in the extended Great Observatories Origins Deep Survey-North (GOODS-N) field covering an area of 0.22 $\rm{deg}^{2}$. The observations were taken using WIRCam on the 3.6-m Canada France Hawaii Telescope (CFHT). Together with the reprocessed $K_{\rm s}$-band image, the $5\sigma$ limiting AB magnitudes (in 2" diameter apertures) are 24.7, 24.2, and 24.4 AB mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. We also release a multi-band photometry and photometric redshift catalog containing 93598 sources. For non-X-ray sources, we obtained a photometric redshift accuracy $\sigma_{\mathrm{NMAD}}=0.036$ with an outlier fraction $\eta = 7.3\%$. For X-ray sources, which are mainly active galactic nuclei (AGNs), we cross-matched our catalog with the updated 2M-CDFN X-ray catalog from Xue et al. (2016) and found that 658 out of 683 X-ray sources have counterparts. $GALEX$ UV data are included in the photometric redshift computation for the X-ray sources to give $\sigma_{\mathrm{NMAD}} = 0.040$ with $\eta=10.5\%$. Our approach yields more accurate photometric redshift estimates compared to previous works in this field. In particular, by adopting AGN-galaxy hybrid templates, our approach delivers photometric redshifts for the X-ray counterparts with fewer outliers compared to the 3D-HST catalog, which fit these sources with galaxy-only templates.

The third-generation gravitational wave detectors are under development by operating the detector in cryogenic temperature to reduce the thermal noise. Silicon and sapphire are promising candidate materials for the test masses and suspension elements due to their remarkable mechanical and thermal properties at cryogenic temperature. Here we present the performances of the cryogenic thermal cycling and strength testing on hydroxide catalysis bonding between sapphire and silicon. Our results suggest that although these two materials have very different coefficients of thermal expansion, but if the flatness and the thermally grown $\mathrm{SiO_2}$ oxidation layer on the silicon surface are controlled well, the bonded samples can still survive thermal cycling from room temperature to 5.5 K. A breaking strength of 3.6$\pm 0.6$ MPa is measured for the bonds between sapphire and silicon with a 190 nm silicon oxidation thickness after cooling cycle. We construct a hybrid sapphire-silicon Fabry-Perot cavity with the developing bonding technique in our lab. The measurement results reveal that the cavity can survive repeated thermal cycling while maintaining a good finesse.

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

The Prime Focus Spectrograph (PFS) is a new optical/near-infrared multi-fiber spectrograph design for the prime focus of the 8.2m Subaru telescope. PFS will cover 1.3 degree diameter field with 2394 fibers to complement the imaging capability of Hyper SuprimeCam (HSC). The prime focus unit of PFS called Prime Focus Instrument (PFI) provides the interface with the top structure of Subaru telescope and also accommodates the optical bench in which Cobra fiber positioners are located. In addition, the acquisition and guiding cameras (AGCs), the optical fiber positioner system, the cable wrapper, the fiducial fibers, illuminator, and viewer, the field element, and the telemetry system are located inside the PFI. The mechanical structure of the PFI was designed with special care such that its deflections sufficiently match those of the HSC's Wide Field Corrector (WFC) so the fibers will stay on targets over the course of the observations within the required accuracy. In this report, the latest status of PFI development will be given including the performance of PFI components, the setup and performance of the integration and testing equipment.

Supermassive binary black holes (BBHs) are unavoidable products of galaxy mergers and are expected to exist in the cores of many quasars. Great effort has been made during the past several decades to search for BBHs among quasars; however, observational evidence for BBHs remains elusive and ambiguous, which is difficult to reconcile with theoretical expectations. In this paper, we show that the distinct optical-to-UV spectrum of Mrk 231 can be well interpreted as emission from accretion flows onto a BBH, with a semimajor axis of ~590AU and an orbital period of ~1.2 year. The flat optical and UV continua are mainly emitted from a circumbinary disk and a mini-disk around the secondary black hole (BH), respectively; and the observed sharp drop off and flux deficit at wavelength lambda ~ 4000-2500 Angstrom is due to a gap (or hole) opened by the secondary BH migrating within the circumbinary disk. If confirmed by future observations, this BBH will provide a unique laboratory to study the interplay between BBHs and accretion flows onto them. Our result also demonstrates a new method to find sub-parsec scale BBHs by searching for deficits in the optical-to-UV continuum among the spectra of quasars.

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multifiber spectrograph with 2394 science fibers distributed across a 1.3-deg diameter field of view at the Subaru 8.2-m telescope. The wide wavelength coverage from 0.38 {\mu}m to 1.26 {\mu}m, with a resolving power of 3000, simultaneously strengthens its ability to target three main survey programs: cosmology, galactic archaeology and galaxy/AGN evolution. A medium resolution mode with a resolving power of 5000 for 0.71 {\mu}m to 0.89 {\mu}m will also be available by simply exchanging dispersers. We highlight some of the technological aspects of the design. To transform the telescope focal ratio, a broad-band coated microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of the cable system, optimizing overall throughput; a fiber with low focal ratio degradation is selected for the fiber-positioner and fiber-slit components, minimizing the effects of fiber movements and fiber bending. Fiber positioning will be performed by a positioner consisting of two stages of piezo-electric rotary motors. The positions of these motors are measured by taking an image of artificially back-illuminated fibers with the metrology camera located in the Cassegrain container; the fibers are placed in the proper location by iteratively measuring and then adjusting the positions of the motors. Target light reaches one of the four identical fast-Schmidt spectrograph modules, each with three arms. The PFS project has passed several project-wide design reviews and is now in the construction phase.

The Prime Focus Spectrograph (PFS) is a new optical/near-infrared multi-fiber spectrograph designed for the prime focus of the 8.2m Subaru telescope. The metrology camera system of PFS serves as the optical encoder of the COBRA fiber motors for the configuring of fibers. The 380mm diameter aperture metrology camera will locate at the Cassegrain focus of Subaru telescope to cover the whole focal plane with one 50M pixel Canon CMOS sensor. The metrology camera is designed to provide the fiber position information within 5{\mu}m error over the 45cm focal plane. The positions of all fibers can be obtained within 1s after the exposure is finished. This enables the overall fiber configuration to be less than 2 minutes.

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multi-fiber spectrograph with 2394 science fibers, which are distributed in 1.3 degree diameter field of view at Subaru 8.2-meter telescope. The simultaneous wide wavelength coverage from 0.38 um to 1.26 um, with the resolving power of 3000, strengthens its ability to target three main survey programs: cosmology, Galactic archaeology, and galaxy/AGN evolution. A medium resolution mode with resolving power of 5000 for 0.71 um to 0.89 um also will be available by simply exchanging dispersers. PFS takes the role for the spectroscopic part of the Subaru Measurement of Images and Redshifts project, while Hyper Suprime-Cam works on the imaging part. To transform the telescope plus WFC focal ratio, a 3-mm thick broad-band coated glass-molded microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of cable system, while one with a better FRD performance is selected for the fiber-positioner and fiber-slit components, given the more frequent fiber movements and tightly curved structure. Each Fiber positioner consists of two stages of piezo-electric rotary motors. Its engineering model has been produced and tested. Fiber positioning will be performed iteratively by taking an image of artificially back-illuminated fibers with the Metrology camera located in the Cassegrain container. The camera is carefully designed so that fiber position measurements are unaffected by small amounts of high special-frequency inaccuracies in WFC lens surface shapes. Target light carried through the fiber system reaches one of four identical fast-Schmidt spectrograph modules, each with three arms. Prototype VPH gratings have been optically tested. CCD production is complete, with standard fully-depleted CCDs for red arms and more-challenging thinner fully-depleted CCDs with blue-optimized coating for blue arms.

Sub-parsec binary massive black holes (BBHs) are long anticipated to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circum-binary disk, each component of the BBH is surrounded by a small accretion disk, and a gap is opened by the secondary component in between the circum-binary disk and the two small disks. Assuming such a BBH structure, we generate mock microlensing light curves for some QSO systems that host BBHs with typical physical parameters. We show that microlensing light curves of a BBH QSO system at the infrared-optical-UV bands can be significantly different from those of corresponding QSO system with a single massive black hole (MBH), mainly because of the existence of the gap and the rotation of the BBH (and its associated small disks) around the center of mass. We estimate the half-light radii of the emission region at different wavelengths from mock light curves and find that the obtained half-light radius vs. wavelength relations of BBH QSO systems can be much flatter than those of single MBH QSO systems at a wavelength range determined by the BBH parameters, such as the total mass, mass ratio, separation, accretion rates, etc. The difference is primarily due to the existence of the gap. Such unique features on the light curves and half-light radius-wavelength relations of BBH QSO systems can be used to select and probe sub-parsec BBHs in a large number of lensed QSOs to be discovered by current and future surveys, including the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Large Synoptic Survey telescope (LSST) and Euclid.

We report JVLA 8-10 GHz ($\lambda$=3.0-3.7 cm) monitoring observations toward the YSO cluster R Coronae Australis (R\,CrA), taken in 2012, from March 15 to September 12. These observations were planned to measure the radio flux variabilities in timescales from 0.5 hours to several days, to tens of days, and up to $\sim$200 days. We found that among the YSOs detectable in individual epochs, in general, the most reddened objects in the \textit{Spitzer} observations show the highest mean 3.5 cm Stokes \textit{I} emission, and the lowest fractional variabilities on $<$200-day timescales. The brightest radio flux emitters in our observations are the two reddest sources IRS7W and IRS7E. In addition, by comparing with observations taken in 1996-1998 and 2005, we found that the radio fluxes of these two sources have increased by a factor $\sim$1.5. The mean 3.5-cm fluxes of the three Class I/II sources IRSI, IRS2, and IRS6 appear to be correlated with their accretion rates derived by a previous near infrared line survey. The weakly accreting Class I/II YSOs, or those in later evolutionary stages, present radio flux variability on $<$0.5-hour timescales. Some YSOs were detected only during occasional flaring events. The source R\,CrA went below our detection limit during a few fading events.

We present new near-IR H2, CO J=2-1, and CO J = 3-2 observations to study outflows in the massive star forming region IRAS 05358+3543. The Canada-France-Hawaii Telescope H2 images and James Clerk Maxwell Telescope CO data cubes of the IRAS 05358 region reveal several new outflows, most of which emerge from the dense cluster of sub-mm cores associated with the Sh 2-233IR NE cluster to the northeast of IRAS 05358. We used Apache Point Observatory (APO) JHK spectra to determine line of sight velocities of the outflowing material. Analysis of archival VLA cm continuum data and previously published VLBI observations reveal a massive star binary as a probable source of one or two of the outflows. We have identified probable sources for 6 outflows and candidate counterflows for 7 out of a total of 11 seen to be originating from the IRAS 05358 clusters. We classify the clumps within Sh 2-233IR NE as an early protocluster and Sh 2-233IR SW as a young cluster, and conclude that the outflow energy injection rate approximately matches the turbulent decay rate in Sh 2-233IR NE.

The Prime Focus Spectrograph (PFS) is a new multi-fiber spectrograph on Subaru telescope. PFS will cover around 1.4 degree diameter field with ~2400 fibers. To ensure precise positioning of the fibers, a metrology camera is designed to provide the fiber position information within 5 {\mu}m error. The final positioning accuracy of PFS is targeted to be better than 10 {\mu}m. The metrology camera will locate at the Cassegrain focus of Subaru telescope to cover the whole focal plane. The PFS metrology camera will also serve for the existing multi-fiber infrared spectrograph FMOS.

The Prime Focus Spectrograph (PFS) of the Subaru Measurement of Images and Redshifts (SuMIRe) project has been endorsed by Japanese community as one of the main future instruments of the Subaru 8.2-meter telescope at Mauna Kea, Hawaii. This optical/near-infrared multi-fiber spectrograph targets cosmology with galaxy surveys, Galactic archaeology, and studies of galaxy/AGN evolution. Taking advantage of Subaru's wide field of view, which is further extended with the recently completed Wide Field Corrector, PFS will enable us to carry out multi-fiber spectroscopy of 2400 targets within 1.3 degree diameter. A microlens is attached at each fiber entrance for F-ratio transformation into a larger one so that difficulties of spectrograph design are eased. Fibers are accurately placed onto target positions by positioners, each of which consists of two stages of piezo-electric rotary motors, through iterations by using back-illuminated fiber position measurements with a wide-field metrology camera. Fibers then carry light to a set of four identical fast-Schmidt spectrographs with three color arms each: the wavelength ranges from 0.38 {\mu}m to 1.3 {\mu}m will be simultaneously observed with an average resolving power of 3000. Before and during the era of extremely large telescopes, PFS will provide the unique capability of obtaining spectra of 2400 cosmological/astrophysical targets simultaneously with an 8-10 meter class telescope. The PFS collaboration, led by IPMU, consists of USP/LNA in Brazil, Caltech/JPL, Princeton, & JHU in USA, LAM in France, ASIAA in Taiwan, and NAOJ/Subaru.

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

FS Tau B is one of the few T Tauri stars that possess a jet and a counterjet as well as an optically-visible cavity wall. We obtained images and spectra of its jet-cavity system in the near-infrared H and K bands using Subaru/IRCS and detected the jet and the counterjet in the [Fe II] 1.644 \mu m line for the first time. Within the inner 2" the blueshifted jet is brighter, whereas beyond ~ 5" the redshifted counterjet dominates the [Fe II] emission. The innermost blueshifted knot is spectrally resolved to have a large line width of ~ 110 km/s, while the innermost redshifted knot appears spectrally unresolved. The velocity ratio of the jet to the counterjet is ~ 1.34, which suggests that FS Tau B is driving an asymmetric jet, similar to those found in several T Tauri Stars. Combining with optical observations in the literature, we showed that the blueshifted jet has lower density and higher excitation than the redshifted counterjet. We suggest that the asymmetry in brightness and velocity is the manifestation of a bipolar outflow driving at different mass-loss rates, while maintaining balance of linear momentum. A full explanation to the asymmetry in the FS Tau B system awaits detail modeling and further investigation of the kinematic structure of the wind-associated cavity walls.

It has been suggested that the high metallicity generally observed in active galactic nuclei (AGNs) and quasars originates from ongoing star formation in the self-gravitating part of accretion disks around the supermassive black holes. We designate this region as the star forming (SF) disk, in which metals are produced from supernova explosions (SNexp) while at the same time inflows are driven by SNexp-excited turbulent viscosity to accrete onto the SMBHs. In this paper, an equation of metallicity governed by SNexp and radial advection is established to describe the metal distribution and evolution in the SF disk. We find that the metal abundance is enriched at different rates at different positions in the disk, and that a metallicity gradient is set up that evolves for steady-state AGNs. Metallicity as an integrated physical parameter can be used as a probe of the SF disk age during one episode of SMBH activity. In the SF disk, evaporation of molecular clouds heated by SNexp blast waves unavoidably forms hot gas. This heating is eventually balanced by the cooling of the hot gas, but we show that the hot gas will escape from the SF disk before being cooled, and diffuse into the BLRs forming with a typical rate of $\sim 1\sunmyr$. The diffusion of hot gas from a SF disk depends on ongoing star formation, leading to the metallicity gradients in BLR observed in AGNs. We discuss this and other observable consequences of this scenario.

We observed the Sh 2-233IR (S233IR) region with better sensitivity in near-infrared than previous studies for this region. By applying statistical subtraction of the back- ground stars, we identified member sources and derived the age and mass of three distinguishable sub-groups in this region: Sh 2-233IR NE, Sh 2-233IR SW, and the "distributed stars" over the whole cloud. Star formation may be occurring sequentially with a relatively small age difference (\sim 0.2-0.3 Myrs) between subclusters. We found that the slopes for initial mass function (Gamma \sim -0.5) of two subclusters are flatter than that of Salpeter, which suggests that more massive stars were preferentially formed in those clusters compared to other Galactic star-forming regions. These subclusters may not result from the overall collapse of the whole cloud, but have formed by triggering before the previous star formation activities disturbed the natal molecular cloud. Addi- tionally, high star formation efficiency (&40%) of the subclusters may also suggest that stars form very efficiently in the center of NE.

We present evidence for unusually strong Halpha emission in galaxies with spectroscopic redshifts in the range of 3.8<z<5.0, over the Great Observatories Origins Deep Survey (GOODS) fields. Among 74 galaxies detected in the Spitzer IRAC 3.6 and 4.5um bands, more than 70% of the galaxies show clear excess at 3.6um compared to the expected flux density from stellar continuum only. We provide evidence that this 3.6um excess is due to Halpha emission redshifted into the 3.6um band, and classify these 3.6um excess galaxies to be Halpha emitter (HAE) candidates. The selection of HAE candidates using an excess in broad-band filters is sensitive to objects whose rest-frame H$\alpha$ equivalent width is larger than 350A, star formation rates of 20-500 Msun/yr. The Halpha-to-UV luminosity ratio of HAEs is on average larger than that of local starbursts. Possible reasons for such strong Halpha emission in these galaxies include different dust extinction properties, young stellar population age, continuous star formation history, low metallicity, and top-heavy stellar initial mass function. Although the correlation between UV slope beta and L(Ha)/L(UV) raises the possibility that HAEs prefer a dust extinction curve which is steeper in the UV, the most dominant factor that results in strong Halpha emission appears to be star formation history. The Halpha equivalent widths of HAEs are large despite their relatively old stellar population ages constrained by SED fitting, suggesting that at least 60% of HAEs produce stars at a constant rate. Under the assumption that the cold gas supply is sustained, HAEs are able to produce ~50% of the observed stellar mass density at z~3 that is caught in massive (M*>10^11 Msun) galaxies. This 'strong Halpha phase' of star formation plays a dominant role in early galaxy growth, being a feasible progenitors of massive red galaxies at lower redshifts.

(Abridged) High redshift galaxies are undergoing intensive evolution of dynamical structure and morphologies. We incorporate the feedback into the dynamical equations through mass dropout and angular momentum transportation driven by the SNexp-excited turbulent viscosity. We numerically solve the equations and show that there can be intensive evolution of structure of the gaseous disk. Secular evolution of the disk shows interesting characteristics that are 1) high viscosity excited by SNexp can efficiently transport the gas from 10kpc to $\sim 1$kpc forming a stellar disk whereas a stellar ring forms for the case with low viscosity; 2) starbursts trigger SMBH activity with a lag $\sim 10^8$yr depending on star formation rates, prompting the joint evolution of SMBHs and bulges; 3) the velocity dispersion is as high as $\sim 100~\kms$ in the gaseous disk. In order to compare the present models with the observed dynamical structure and images, we use the incident continuum from the simple stellar synthesis (GALAXEV) and CLOUDY to calculate emission line ratios of H$\alpha$, H$\beta$, $\OIII$ and $\NII$, and H$\alpha$ brightness of gas photoionized by young massive stars formed on the disks. The models can produce the main features of emission from star forming galaxies and the observed relation between turbulent velocity and the H$\alpha$ brightness. We successfully apply the present model to BX 389 and BX 482 observed in SINS high$-z$ sample, which are bulge and disk-dominated, respectively. High viscosity excited by SNexp is able to efficiently transport the gas into a bulge to maintain high star formation rates, or, to form a stellar ring close enough to the bulge so that it immigrates into the bulge of its host galaxy. This leads to a fast growing bulge. Implications and future work of the present models have been extensively discussed for galaxy formation.

Self-gravitating accretion disks collapse to star-forming(SF) regions extending to the inner edge of the dusty torus in active galactic nuclei (AGNs). A full set of equations including feedback of star formation is given to describe the dynamics of the regions. We explore the role of supernovae explosion (SNexp), acting to excite turbulent viscosity, in the transportation of angular momentum in the regions within 1pc scale. We find that accretion disks with typical rates in AGNs can be driven by SNexp in the regions and metals are produced spontaneously. The present model predicts a metallicity--luminosity relationship consistent with that observed in AGNs. As relics of SF regions, a ring (or belt) consisting of old stars remains for every episode of supermassive black hole activity. We suggest that multiple stellar rings with random directions interact and form a nuclear star cluster after episodes driven by star formation.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

We present a near-infrared extinction study of nine dense cores at evolutionary stages between starless to Class I. Our results show that the density structure of all but one observed cores can be modeled with a single power law rho \propto r^p between ~ 0.2R-R of the cores. The starless cores in our sample show two different types of density structures, one follows p ~ -1.0 and the other follows p ~ -2.5, while the protostellar cores all have p ~ -2.5. The similarity between the prestellar cores with p ~ -2.5 and protostellar cores implies that those prestellar cores could be evolving towards the protostellar stage. The slope of p ~ -2.5 is steeper than that of an singular isothermal sphere, which may be interpreted with the evolutionary model of cores with finite mass.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

Motivated by Genzel et al.'s observations of high-redshift star-forming galaxies, containing clumpy and turbulent rings or disks, we build a set of equations describing the dynamical evolution of gaseous disks with inclusion of star formation and its feedback. Transport of angular momentum is due to "turbulent" viscosity induced by supernova explosions in the star formation region. Analytical solutions of the equations are found for the initial cases of a gaseous ring and the integrated form for a gaseous disk, respectively. For a ring with enough low viscosity, it evolves in a slow processes of gaseous diffusion and star formation near the initial radius. For a high viscosity, the ring rapidly diffuses in the early phase. The diffusion drives the ring into a region with a low viscosity and start the second phase undergoing pile-up of gas at a radius following the decreased viscosity torque. The third is a sharply deceasing phase because of star formation consumption of gas and efficient transportation of gas inward forming a stellar disk. We apply the model to two $z\sim 2$ galaxies BX 482 and BzK 6004, and find that they are undergoing a decline in their star formation activity.

The growth of supermassive black holes (BHs) located at the centers of their host galaxies comes mainly from accretion of gas, but how to fuel them remains an outstanding unsolved problem in quasar evolution. This issue can be elucidated by quantifying the radiative efficiency parameter ($\eta$) as a function of redshift, which also provides constraints on the average spin of the BHs and its possible evolution with time. We derive a formalism to link $\eta$ with the luminosity density, BH mass density, and duty cycle of quasars, quantities we can estimate from existing quasar and galaxy survey data. We find that $\eta$ has a strong cosmological evolution: at z~2, $\eta \approx 0.3$, and by $z\approx 0$ it has decreased by an order of magnitude, to $\eta\approx 0.03$. We interpret this trend as evolution in BH spin, and we appeal to episodic, random accretion as the mechanism for reducing the spin. The observation that the fraction of radio-loud quasars decreases with increasing redshift is inconsistent with the popular notion that BH spin is a critical factor for generating strong radio jets. In agreement with previous studies, we show that the derived history of BH accretion closely follows the cosmic history of star formation, consistent with other evidence that BHs and their host galaxies coevolve.

Tracing the star formation history in circumnuclear regions (CNRs) is a key step towards understanding the starburst-AGN connection. However, bright nuclei outshining the entire host galaxy prevent the analysis of the stellar populations of CNRs around type-I AGNs. Obscuration of the nuclei by the central torus provides an unique opportunity to study the stellar populations of AGN host galaxies. We assemble a sample of 10, 848 type-II AGNs with a redshift range of $0.03\le z\le 0.08$ from the Sloan Digital Sky Survey's Data Release 4, and measure the mean specific star formation rates (SSFRs) over the past 100Myr in the central $\sim1-2$ kpc . We find a tight correlation between the Eddington ratio ($\lambda$) of the central black hole (BH) and the mean SSFR, strongly implying that supernova explosions (SNexp) play a role in the transportation of gas to galactic centers. We outline a model for this connection by accounting for the role of SNexp in the dynamics of CNRs. In our model, the viscosity of turbulence excited by SNexp is enhanced, and thus angular momentum can be efficiently transported, driving inflows towards galactic centers. Our model explains the observed relation $\lambda \propto \rm SSFR^{1.5-2.0}$, suggesting that AGN are triggered by SNexp in CNRs.

We use deep far-IR, submm, radio and X-ray imaging and mid-IR spectroscopy to explore the nature of a sample of Spitzer-selected dust-obscured galaxies (DOGs) in GOODS-N. A sample of 79 galaxies satisfy the criteria R-[24]>14 (Vega) down to S24>100 microJy. Twelve of these galaxies have IRS spectra available which we use to measure redshifts and classify these objects as being dominated by star formation or active galactic nuclei (AGN) activity in the mid-IR. The IRS spectra and Spitzer photometric redshifts confirm that the DOGs lie in a tight redshift distribution around z~2. Based on mid-IR colors, 80% of DOGs are likely dominated by star formation; the stacked X-ray emission from this sub-sample of DOGs is also consistent with star formation. Since only a small number of DOGs are individually detected at far-IR and submm wavelengths, we use a stacking analysis to determine the average flux from these objects and plot a composite IR (8-1000 microns) spectral energy distribution (SED). The average luminosity of these star forming DOGs is LIR~1e12 Lsun. We compare the average star forming DOG to the average bright (S850>5 mJy) submillimeter galaxy (SMG); the S24>100 microJy DOGs are 3 times more numerous but 8 times less luminous in the IR. The far-IR SED shape of DOGs is similar to that of SMGs (average dust temperature of around 30 K) but DOGs have a higher mid-IR to far-IR flux ratio. The average star formation-dominated DOG has a star formation rate of 200 Msun/yr which, given their space density, amounts to a contribution of 0.01 Msun/yr/Mpc3 (or 5-10%) to the star formation rate density at z~2. We use the composite SED to predict the average flux of DOGs in future Herschel/PACS 100 micron and SCUBA-2 450 micron surveys and show that the majority of them will be detected.

Episodic activity of quasars is driving growth of supermassive black holes (SMBHs) via accretion of baryon gas. In this Letter, we develop a simple method to analyse the duty cycle of quasars up to redshift $z\sim 6$ universe from luminosity functions (LFs). We find that the duty cycle below redshift $z\sim 2$ follows the cosmic history of star formation rate (SFR) density. Beyond $z\sim 2$, the evolutionary trends of the duty cycle are just opposite to that of the cosmic SFR density history, implying the role of feedback from black hole activity. With the duty cycle, we get the net lifetime of quasars ($z\le 5$) about $\sim 10^9$yrs. Based on the local SMBHs, the mean mass of SMBHs is obtained at any redshifts and their seeds are of $10^5\sunm$ at the reionization epoch ($z_{\rm re}$) of the universe through the conservation of the black hole number density in comoving frame. We find that primordial black holes ($\sim 10^3\sunm$) are able to grow up to the seeds via a moderate super-Eddington accretion of $\sim 30$ times of the critical rate from $z=24$ to $z_{\rm re}$. Highly super-Eddington accretion onto the primordials is not necessary.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

Feedback from black hole activity is widely believed to play a key role in regulating star formation and black hole growth. A long-standing issue is the relation between the star formation and fueling the supermassive black holes in active galactic nuclei (AGNs). We compile a sample of 57 Seyfert galaxies to tackle this issue. We estimate the surface densities of gas and star formation rates in circumnuclear regions (CNRs). Comparing with the well-known Kennicutt-Schmidt (K-S) law, we find that the star formation rates in CNRs of most Seyfert galaxies are suppressed in this sample. Feedback is suggested to explain the suppressed star formation rates.

We have searched three hot molecular cores for submillimeter emission from the nucleic acid building-block pyrimidine. We obtain upper limits to the total pyrimidine (beam-averaged) column densities towards Sgr B2(N), Orion KL and W51 e1/e2 of 1.7E+14 cm^{-2}, 2.4E+14 cm^{-2} and 3.4E+14 cm^{-2}, respectively. The associated upper limits to the pyrimidine fractional abundances lie in the range (0.3-3)E-10. Implications of this result for interstellar organic chemistry, and for the prospects of detecting nitrogen heterocycles in general, are briefly discussed.