Horizon Run 5 (HR5) is a cosmological hydrodynamics simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1 kpc. This enormous dynamic range allows us to simultaneously capture the physics of the cosmic web on very large scales and account for the formation and evolution of dwarf galaxies on much smaller scales. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of 1049x114x114 Mpc^3. The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, reionization, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and active galactic nuclei (AGN) feedback in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OpenMP version of the RAMSES code, specifically targeted for modern many-core many thread parallel architectures. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder to analyse the large outputs of HR5. The simulation successfully reproduces many observations, such as the cosmic star formation history, connectivity of galaxy distribution and stellar mass functions. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one should carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe.

We present the variability study of the lowest-luminosity Seyfert 1 galaxy NGC 4395 based on the photometric monitoring campaigns in 2017 and 2018. Using 22 ground-based and space telescopes, we monitored NGC 4395 with a $\sim$5 minute cadence during a period of 10 days and obtained light curves in the UV, V, J, H, and K/Ks bands as well as the H$\alpha$ narrow-band. The RMS variability is $\sim$0.13 mag on \emph{Swift}-UVM2 and V filter light curves, decreasing down to $\sim$0.01 mag on K filter. After correcting for continuum contribution to the H$\alpha$ narrow-band, we measured the time lag of the H$\alpha$ emission line with respect to the V-band continuum as ${55}^{+27}_{-31}$ to ${122}^{+33}_{-67}$ min. in 2017 and ${49}^{+15}_{-14}$ to ${83}^{+13}_{-14}$ min. in 2018, depending on the assumption on the continuum variability amplitude in the H$\alpha$ narrow-band. We obtained no reliable measurements for the continuum-to-continuum lag between UV and V bands and among near-IR bands, due to the large flux uncertainty of UV observations and the limited time baseline. We determined the AGN monochromatic luminosity at 5100\AA\ $\lambda L_\lambda = \left(5.75\pm0.40\right)\times 10^{39}\,\mathrm{erg\,s^{-1}}$, after subtracting the contribution of the nuclear star cluster. While the optical luminosity of NGC 4395 is two orders of magnitude lower than that of other reverberation-mapped AGNs, NGC 4395 follows the size-luminosity relation, albeit with an offset of 0.48 dex ($\geq$2.5$\sigma$) from the previous best-fit relation of Bentz et al. (2013).

Using high resolution hydrodynamical cosmological simulations, we study the gas accretion history of low mass halos located in a field-like, low density environment. We track their evolution individually from the early, pre-reionization era, through reionization, and beyond until $z=0$. Before reionization, low mass halos accrete cool cosmic web gas at a very rapid rate, often reaching the highest gas mass they will ever have. But when reionization occurs, we see that almost all halos lose significant quantities of their gas content, although some respond less quickly than others. We find that the response rate is influenced by halo mass first, and secondarily by their internal gas density at the epoch of reionization. Reionization also fully ionises the cosmic web gas by z$\sim$6. As a result, the lowest mass halos (M$\sim$10$^6~h^{-1}$M$_\odot$ at $z=6$) can never again re-accrete gas from the cosmic web, and by $z\sim5$ have lost all their internal gas to ionisation, resulting in a halt in star formation at this epoch. However, more massive halos can recover from their gas mass loss, and re-accrete ionised cosmic web gas. We find the efficiency of this re-accretion is a function of halo mass first, followed by local surrounding gas density. Halos that are closer to the cosmic web structure can accrete denser gas more rapidly. We find that our lower mass halos have a sweet spot for rapid, dense gas accretion at distances of roughly 1-5 virial radii from the most massive halos in our sample ($>$10$^8~h^{-1}$M$_\odot$), as these tend to be embedded deeply within the cosmic web.

The disruption of substructure in galaxy clusters likely plays an important role in shaping the cluster population as a significant fraction of cluster galaxies today have spent time in a previous host system, and thus may have been pre-processed. Once inside the cluster, group galaxies face the combined environmental effects from group and cluster - so called 'post-processing'. We investigate these concepts, by tracking the evolution of satellites and their hosts after entering the cluster and find that tidal forces during their first pericentric passage are very efficient at breaking up groups, preferentially removing satellites at larger distances from their hosts. 92.2% of satellites whose host has passed pericentre will leave their host by $z=0$, typically no later than half a Gyr after pericentric passage. We find satellites leave with high velocities, and quickly separate to large distances from their hosts, making their identification within the cluster population challenging. Those few satellites ($\sim$7.8%) that remain bound to their hosts after a pericentric passage are typically found close to their host centres. This implies that substructure seen in clusters today is very likely on first infall into the cluster, and yet to pass pericentre. This is even more likely if the substructure is extended, with satellites beyond R$_{200}$ of their host. We find the cluster dominates the tidal mass loss and destruction of satellites, and is responsible for rapidly halting the accretion of new satellites onto hosts once they reach 0.6-0.8 R$_{200}$ radii from the cluster.

Active galactic nuclei (AGNs) show a correlation between the size of the broad line region (BLR) and the monochromatic continuum luminosity at 5100 \AA, allowing black hole mass estimation based on single-epoch spectra. However, the validity of the correlation is yet to be clearly tested for high-luminosity AGNs. We present the first reverberation-mapping results of the Seoul National University AGN monitoring program (SAMP), which is designed to focus on luminous AGNs for probing the high end of the size-luminosity relation. We report time lag measurements of two AGNs, namely, 2MASS J10261389+5237510 and SDSS J161911.24+501109.2, using the light curves obtained over a $\sim$1000 day period with an average cadence of $\sim$10 and $\sim$20 days, respectively for photometry and spectroscopy monitoring. Based on a cross-correlation analysis and H$\beta$ line width measurements, we determine the H$\beta$ lag as $41.8^{+4.9}_{-6.0}$ and $52.6^{+17.6}_{-14.7}$ days in the observed-frame, and black hole mass as $3.65^{+0.49}_{-0.57} \times 10^7 M_{\odot}$ and $23.02^{+7.81}_{-6.56} \times 10^7 M_{\odot}$, respectively for 2MASS J1026 and SDSS J1619.

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 use the IllustrisTNG cosmological hydrodynamical simulation to study the evolution of star formation rate (SFR)-density relation over cosmic time. We construct several samples of galaxies at different redshifts from $z=2.0$ to $z=0.0$, which have the same comoving number density. The SFR of galaxies decreases with local density at $z=0.0$, but its dependence on local density becomes weaker with redshift. At $z\gtrsim1.0$, the SFR of galaxies increases with local density (reversal of the SFR-density relation), and its dependence becomes stronger with redshift. This change of SFR-density relation with redshift still remains even when fixing the stellar masses of galaxies. The dependence of SFR on the distance to a galaxy cluster also shows a change with redshift in a way similar to the case based on local density, but the reversal happens at a higher redshift, $z\sim1.5$, in clusters. On the other hand, the molecular gas fraction always decreases with local density regardless of redshift at $z=0.0-2.0$ even though the dependence becomes weaker when we fix the stellar mass. Our study demonstrates that the observed reversal of the SFR-density relation at $z\gtrsim1.0$ can be successfully reproduced in cosmological simulations. Our results are consistent with the idea that massive, star-forming galaxies are strongly clustered at high redshifts, forming larger structures. These galaxies then consume their gas faster than those in low-density regions through frequent interactions with other galaxies, ending up being quiescent in the local universe.

We present a spatially-resolved analysis of ionized and molecular gas in a nearby Seyfert 2 galaxy NGC 5728, using the VLT/MUSE and ALMA data. We find ionized gas outflows out to ~kpc scales, which encounter the star formation ring at 1 kpc radius. The star formation rate of the encountering region is significantly high (~1.8 M$_{\rm sol}/yr/kpc^2$) compared to other regions in the ring. In contrast, the CO (2-1) emission is significantly weaker by a factor of ~3.5, indicating very high star formation efficiency. These results support the positive feedback scenario that the AGN-driven outflows compress the ISM in the ring, enhancing the star formation activity. In addition, we detect outflow regions outside of spiral arms, in which gas is likely to be removed from the spiral arms and no clear sign of star formation is detected. The overall impact of AGN outflows on the global star formation in NGC 5728 is limited, suggesting the feedback of the low-luminosity AGN is insignificant.

While the reverberation mapping technique is the best available method for measuring black hole mass in active galactic nuclei (AGN) beyond the local volume, this method has been mainly applied to relatively low-to-moderate luminosity AGNs at low redshift. We present the strategy of the Seoul National University AGN Monitoring Project, which aims at measuring the time delay of the \Hb\ line emission with respect to AGN continuum, using a sample of relatively high luminosity AGNs out to z $\sim$0.5. We present simulated cross correlation results based on a number of mock light curves, in order to optimally determine monitoring duration and cadence. We describe our campaign strategy based on the simulation results and the availability of observing facilities. We present the sample selection, and the properties of the selected 100 AGNs, including the optical luminosity, expected time lag, black hole mass, and Eddington ratio.

We trace the cosmological origin of satellites around a dwarf galaxy using a very high resolution (12~pc/h) cosmological hydrodynamic zoom simulation. To realistically describe the formation and evolution of small-mass stellar satellites, our model includes a full baryonic physics treatment including a recipe for UV self-shielding. We find that the mini-halos form objects resembling dwarf galaxies. Despite our high resolution, none of our 'mini-halos' form objects that might be considered globular cluster-like. Instead such objects are formed in the host galaxy's gaseous disk. We investigate why some form stars more efficiently than others. We find that the majority of their star forming gas is accreted after reionisation, thus the survival of a mini-halo's gas to reionisation is not an important factor. Instead, the key factor seems to be the ability for a mini-halo to cool its recently accreted gas, which is more efficient in more massive halos. We find halos in denser environments suffer more mergers, enabling them to grow their mass such that cooling of accreted gas can occur efficiently. Although the host galaxy is only a dwarf galaxy itself, we find that ram pressure is an efficient means by which accreted mini-halos lose their gas content, both by interacting with hot halo gas but also in direct collisions with the gas disk of the host. The satellites are also disrupted by the tidal forces near the center of the host galaxy. Compared to the disrupted satellites, surviving satellites are relatively more massive, but tend to infall later into the host galaxy, thus reducing the time they are subjected to destructive environmental mechanisms and dynamical friction. In summary, our results suggest that the characteristics of satellites are mainly determined by their ability to efficiently cool gas that is accreted in the redshift range z=3-5, prior to their infall into the host galaxy.

Kinematic studies have produced accurate measurements of the total dark matter mass and mean dark matter density within the optical extent of galaxies, for large samples of objects. Here we consider theoretical predictions for the latter quantity, $\bar{\rho}_{dm}$, measured within the isophotal radius $R_{23.5}$, for isolated halos with universal density profiles. Through a combination of empirical scaling relations, we show that $\bar{\rho}_{dm}$ is expected to depend weakly on halo mass and strongly on redshift. When galaxy halos fall into larger groups or clusters they become tidally stripped, reducing their total dark matter mass, but this process is expected to preserve central density until an object is close to disruption. We confirm this with collisonless simulations of cluster formation, finding that subhalos have values of $\bar{\rho}_{dm}$ close to the analytic predictions. This suggests that $\bar{\rho}_{dm}$ may be a useful indicator of infall redshift onto the cluster. We test this hypothesis with data from the SHIVir survey, which covers a reasonable fraction of the Virgo cluster. We find that galaxies with high $\bar{\rho}_{dm}$ do indeed trace the densest regions of the cluster, with a few notable exceptions. Samples selected by environment have higher densities at a significance of 3.5-4$\sigma$, while samples selected by density are more clustered at 3-3.5$\sigma$ significance. We conclude that halo density can be a powerful tracer of the assembly history of clusters and their member galaxies.

The motions of gas and stars in the nuclei of nearby large galaxies have demonstrated that massive black holes are common and that their masses strongly correlate with the stellar velocity dispersion $\sigma_{\star}$ of the bulge. This correlation suggests that massive black holes and galaxies influence each other's growth. Dynamical measurements are less reliable when the sphere of influence is unresolved, thus it remains unknown whether this correlation exists in galaxies much smaller than the Milky Way, as well as what fraction of these galaxies have central black holes. Light echoes from photoionized clouds around accreting black holes, in combination with the velocity of these clouds, yield a direct mass measurement that circumvents this difficulty. Here we report an exceptionally low reverberation delay of $83\pm14$ minutes between variability in the accretion disk and high velocity H$\alpha$ emission from the nucleus of the bulgeless dwarf galaxy NGC~4395. Combined with the H$\alpha$ line-of-sight velocity dispersion $\sigma_{\rm line}=426\pm1$~km~s$^{-1}$, this lag determines a mass of about 10,000~$M_{\odot}$ for the black hole. This mass is among the smallest central black hole masses reported, near the low end of expected masses for heavy "seeds", and the best direct mass measurement for a galaxy of this size. Despite the lack of a bulge, NGC~4395 is consistent with the $M_{\rm BH} - \sigma_{\star}$ relation when $\sigma_{\star}$ is measured from the central region. This indicates that the relation need not originate from hierarchical galaxy assembly nor from black hole feedback.

AGN-driven outflows are considered as one of the processes driving the co-evolution of supermassive black holes with their host galaxies. We present integral field spectroscopy of six Type 2 AGNs at z < 0.1, which are selected as AGNs without strong outflows based on the kinematics of [Oiii] gas. Using spatially resolved data, we investigate the ionized gas kinematics and photoionization properties in comparison with AGNs with strong outflows. We find significant differences between the kinematics of ionized gas and stars for two AGNs, which indicates the presence of AGN-driven outflows. Nevertheless, the low velocity and velocity dispersion of ionized gas indicate relatively weak outflows in these AGNs. Our results highlight the importance of spatially-resolved observation in investigating gas kinematics and identifying the signatures of AGN-driven outflows. While it is unclear what determines the occurrence of outflows, we discuss the conditions and detectability of AGN-driven outflows based on a larger sample of AGNs with and without outflows, suggesting the importance of gas content in the host galaxies.

The FeII/MgII line flux ratio has been used to investigate the chemical evolution of high redshift active galactic nuclei (AGNs). No strong evolution has been found out to z~6, implying that the type 1a supernova activity has been already occurred in early universe. However, the trend of no evolution can be caused by the sample selection bias since the previous studies have utilized mostly very luminous AGNs, which may be already chemically matured at the observed redshift. As motivated by the previously reported correlation between AGN luminosity and metallicity, we investigate the FeII/MgII flux ratio over a large dynamic range of luminosity, by adding a new sample of 12 quasars at z~3, of which the lower luminosity limit is more than 1 dex smaller than that of the previously studied high-z quasars. Based on the Gemini/GNIRS observations, we find that the seven low-luminosity quasars with a mean bolometric luminosity log L_bol~46.5 \pm 0.2 has an average FeII/MgII ratio of 0.68 \pm 0.11 dex. This ratio is comparable to that of high-luminosity quasars (log L_bol~47.3 \pm 0.3) in our sample (i.e., FeII/MgII ratio of 0.59 \pm 0.15 dex) and that of the previously studied high-luminosity quasars at higher redshifts. One possible scenario is that the low-luminosity quasars in our sample are still relatively luminous and already chemically matured. To search for chemically-young AGNs, and fully understand the chemical evolution based on the FeII/MgII flux ratio, it is required to investigate much lower-luminosity AGNs.

Outflows in active galactic nuclei (AGNs) are crucial to understand in investigating the co-evolution of supermassive black holes (SMBHs) and their host galaxies since outflows may play an important role as an AGN feedback mechanism. Based on the archival UV spectra obtained with HST and IUE, we investigate outflows in the broad-line region (BLR) in low-redshift AGNs (z < 0.4) through the detailed analysis of the velocity profile of the CIV emission line. We find a dependence of the outflow strength on the Eddington ratio and the BLR metallicity in our low-redshift AGN sample, which is consistent with the earlier results obtained for high-redshift quasars. These results suggest that the BLR outflows, gas accretion onto SMBH, and past star-formation activity in the host galaxies are physically related in low-redshift AGNs as in powerful high-redshift quasars.

As an evolutionary phase of galaxies, active galactic nuclei (AGN) over a large range of redshift have been utilized for understanding cosmic evolution. In particular, the population and evolution of AGNs have been investigated through the study of the cosmic X-ray background in various fields. As one of the deep fields observed by Chandra with a total of 2.8 Msec exposures, Abell 133 is a special region to investigate AGNs, providing a testbed for probing the environmental effect on AGN trigger, since cluster environments can be different from field environments. The achieved flux limit of data at 50% completeness level of 6.95x10^(-16), 1.43x10^(-16), and 1.57x10^(-15) ergs cm^(-2) in 0.5-8, 0.5-2, and 2-8 keV. Using the wavdetect and no-source binomial probability (i.e., p<0.007), we analyze the combined \textit{Chandra} image, detecting 1617 (in 0.5-8 keV), 1324 (in 0.5-2 keV), and 1028 (in 2-8 keV), X-ray point sources in the Abell 133 region. Here, we present the X-ray point source catalogue with the source fluxes, which can be combined with multiwavelength data for future works. We find that the number count distribution of the X-ray point sources is well reproduced with a broken power-law model while the best-fit model parameters are sensitive to the fitting range of the number count distribution. Finally, we find an excess of the number density (decrease of AGN fraction) at the central region of the cluster, which reflects the effect of dense environments on AGN trigger, as similarly found in studies of other galaxy clusters.

Existence of hypothetical fourth neutrino, so-called sterile neutrino, is one of open issues in the particle and neutrino physics. This fourth neutrino is a candidate for explaining some anomalies reported in LSND, MiniBoone, reactor experiments, and gallium experiments. To search for the existence of the sterile neutrino, we report detailed analysis of a feasible experiment for short baseline electron antineutrino (${\bar{\nu}}_{e}$) disappearance study, in which a ${\bar{\nu}}_{e}$ source from $^{8}$Li generator is considered under non-accelerator system. For $^{8}$Li production, we suggest to use $^{252}$Cf source as an intense neutron emitter, by which one can produce $^{8}$Li isotope through $^{7}$Li(n,$\gamma$)$^{8}$Li reaction, effectively. Using the $^{8}$Li generator, one does not need any accelerator or reactor facilities because the generator can be placed on any present and/or planned neutrino detectors as closely as possible. For the effect of the possible sterile neutrinos, we estimate expected neutrino flux and event rates from the neutrino source scheme, and show neutrino disappearance features and possible reaction rate changes by the sterile neutrino using the spectral shape analysis.

Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is established using a parametrisation-independent extraction of the bias function. This average bias model is then combined with the dark matter PDF, neglecting any scatter around it: it nevertheless yields an excellent model for densities-in-cells statistics of mass tracers that is parametrised in terms of the underlying dark matter variance and three bias parameters. The procedure is validated on measurements of both the one and two point statistics of subhalo densities in the state-of-the-art Horizon Run 4 simulation showing excellent agreement for measured dark matter variance and bias parameters. Finally, it is demonstrated that this formalism allows for a joint estimation of the nonlinear dark matter variance and the bias parameters using solely the statistics of subhaloes. Having verified that galaxy counts in hydrodynamical simulations sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work provides important steps towards making theoretical predictions for density-in-cells statistics applicable to upcoming galaxy surveys like Euclid or WFIRST.

We have measured the probability distribution function (PDF) of cosmic matter density field from a suite of N-body simulations. We propose the generalized normal distribution of version 2 (Nv2) as an alternative fitting formula to the well-known log-normal distribution. We find that Nv2 provides significantly better fit than the log-normal distribution for all smoothing radii (2, 5, 10, 25 [Mpc/h]) that we studied. The improvement is substantial in the underdense regions. The development of non- Gaissianities in the cosmic matter density field is captured by continuous evolution of the skewness and shifts parameters of the Nv2 distribution. We present the redshift evolution of these parameters for aforementioned smoothing radii and various background cosmology models. All the PDFs measured from large and high-resolution N-body simulations that we use in this study can be obtained from a Web site at https://astro.kias.re.kr/jhshin.

We present hydrodynamic simulations of gas clouds inflowing from the disk to a few hundred parsec region of the Milky Way. A gravitational potential is generated to include realistic Galactic structures by using thousands of multipole expansions that describe 6.4 million stellar particles of a self-consistent Galaxy simulation. We find that a hybrid multipole expansion model, with two different basis sets and a thick disk correction, accurately reproduces the overall structures of the Milky Way. Through non-axisymmetric Galactic structures of an elongated bar and spiral arms, gas clouds in the disk inflow to the nuclear region and form a central molecular zone (CMZ)-like nuclear ring. We find that the size of the nuclear ring evolves into ~240 pc at T~1500 Myr, regardless of the initial size. For most simulation runs, the rate of gas inflow to the nuclear region is equilibrated to ~0.02 M_sun/yr. The nuclear ring is off-centered, relative to the Galactic center, by the lopsided central mass distribution of the Galaxy model, and thus an asymmetric mass distribution of the nuclear ring arises accordingly. The vertical asymmetry of the the Galaxy model also causes the nuclear ring to be tilted along the Galactic plane. During the first ~100 Myr, the vertical frequency of the gas motion is twice that of the orbital frequency, thus the projected nuclear ring shows a twisted, infinity-like shape.

To investigate the existence of sterile neutrino, we propose a new neutrino production method using $^{13}$C beams and a $^{9}$Be target for short-baseline electron antineutrino (${\bar{\nu}}_{e}$) disappearance study. The production of secondary unstable isotopes which can emit neutrinos from the $^{13}$C + $^{9}$Be reaction is calculated with three different nucleus-nucleus (AA) reaction models. Different isotope yields are obtained using these models, but the results of the neutrino flux are found to have unanimous similarities. This feature gives an opportunity to study neutrino oscillation through shape analysis. In this work, expected neutrino flux and event rates are discussed in detail through intensive simulation of the light ion collision reaction and the neutrino flux from the beta decay of unstable isotopes followed by this collision. Together with the reactor and accelerator anomalies, the present proposed ${\bar{\nu}}_{e}$ source is shown to be a practically alternative test of the existence of the $\Delta m^{2}$ $\sim$ 1 eV$^{2}$ scale sterile neutrino.

We perform a comprehensive study of X-ray cavities using a large sample of X-ray targets selected from the Chandra archive. The sample is selected to cover a large dynamic range including galaxy clusters, groups, and individual galaxies. Using $\beta$-modeling and unsharp masking techniques, we investigate the presence of X-ray cavities for 133 targets that have sufficient X-ray photons for analysis. We detect 148 X-ray cavities from 69 targets and measure their properties, including cavity size, angle, and distance from the center of the diffuse X-ray gas. We confirm the strong correlation between cavity size and distance from the X-ray center similar to previous studies (i.e., Birzan et al. 2004; Diehl et al. 2008; Dong et al. 2010). We find that the detection rates of X-ray cavities are similar among galaxy clusters, groups and individual galaxies, suggesting that the formation mechanism of X-ray cavities is independent of environment.

The Quintuplet and Arches clusters, which were formed in the harsh environment of the Galactic Center (GC) a few million years ago, have been excellent targets for studying the effects of a star-forming environment on the initial mass function (IMF). In order to estimate the shape of the low-end IMF of the Arches cluster, Shin & Kim devised a novel photometric method that utilizes pixel intensity histograms (PIHs) of the observed images. Here, we apply the PIH method to the Quintuplet cluster and estimate the shape of its low-end IMF below the magnitude of completeness limit as set by conventional photometry. We found that the low-end IMF of the Quintuplet is consistent with that found for the Arches cluster--Kroupa MF, with a significant number of low-mass stars below 1 $M_{\odot}$. We conclude that the most likely IMFs of the Quintuplet and the Arches clusters are not too different from the IMFs found in the Galactic disc. We also find that the observed PIHs and stellar number density profiles of both clusters are best reproduced when the clusters are assumed to be at three-dimensional distances of approximately 100 pc from the GC.

Production of a neutrino source through proton induced reaction is studied by using the particle transport code, GEANT4. Unstable isotope such as $^{27}$Si can be produced when $^{27}$Al target is bombarded by 15 MeV energetic proton beams. Through the beta decay process of the unstable isotope, a new electron-neutrino source in the 1.0 $\sim$ 5.0 MeV energy range is obtained. Proton induced reactions are simulated with JENDL High Energy File 2007 (JENDL/HE-2007) data and other nuclear data. For radioactive decay processes, we use "G4RadioactiveDecay" model based on the Evaluated Nuclear Structure Data File (ENSDF). We suggest target systems required for future's solar neutrino experiments, in particular, for the vacuum-matter transition region. As for the detection system of the new neutrino source, we evaluate reaction rates for available radiochemical detectors and LENA type scintillator detector. Possibility of detecting sterile neutrinos is also discussed.

The initial mass function (IMF) of the Arches cluster, which was formed a few million years ago in the harsh environment of the Galactic center (GC), has long been a target of interest to those who study the GC and the theory of star formation. The distinct star-forming conditions in the GC might have caused the cluster to have a shallower slope or an elevated lower mass cutoff in its IMF. But its mass function has been revealed only down to 1-2 Msun (the lower limit of resolved stars), and the low- end mass function of the Arches is still unknown. To estimate the unresolved part of the Arches mass function, we have devised a novel photometric method that involves the histogram of pixel intensities in the observed image, which contains information on the unresolved, faint stars. By comparing the pixel intensity histograms (PIHs) of numerous artificial images constructed from model IMFs with the observed PIH, we find that the best-fit model IMF for the Arches cluster has a cutoff mass less than or similar to 0.1 Msun and a shape very close to that of the Kroupa MF. Our findings imply that the IMF of the Arches cluster is similar to those found in the Galactic disk.

We have developed a parallel cosmological hydrodynamic simulation code designed for the study of formation and evolution of cosmological structures. The gravitational force is calculated using the TreePM method and the hydrodynamics is implemented based on the smoothed particle hydrodynamics. The initial displacement and velocity of simulation particles are calculated according to second-order linear perturbation theory using the power spectra of dark matter and baryonic matter. The initial background temperature is given by Recfast and the temperature uctuations at the initial particle position are determined by the adiabatic model. We use a time-limiter scheme over the individual time steps to capture shock-fronts and to ease the time-step tension between the shock and preshock particles. We also include the astrophysical gas processes of radiative heating/cooling, star formation, metal enrichment, and supernova feedback. We have tested the code in several standard cases such as one-dimensional Riemann problems, Kelvin-Helmholtz, and Sedov blast wave instability. Star formation on the galactic disk is investigated to check whether the Schmidt-Kennicutt relation is properly recovered. We also study global star formation histories at different simulation resolutions and compare them with observations.

Large-scale structure of the universe is a useful cosmological probe of the primordial non-Gaussianity and the expansion history of the universe because its topology does not change with time in the linear regime in the standard paradigm of structure formation. However, when the topology of iso-density contour surfaces is measured from an observational data, many systematic effects are introduced due to the finite size of pixels used to define the density field, non-linear gravitational evolution, redshift-space distortion, shot noise (discrete sampling), and bias in the distribution of the density field tracers. We study the various systematic effects on the genus curve to a great accuracy by using the Horizon Run 2 simulation of a {\Lambda}CDM cosmology. We numerically measure the genus curve from the gravitationally evolved matter and dark matter halo density fields. It is found that all the non-Gaussian deviations due to the systematic effects can be modeled by using a few low-order Hermite polynomials from H0 to H4. We compare our results with the analytic theories whenever possible, and find many new terms in the Hermite series that are making significant contributions to the non-Gaussian deviations. In particular, it is found that the amplitude drop of the genus curve due to the non-linear gravitational evolution can be accurately modeled by two terms H0 and H2 with coefficients both proportional to \sigma_0^2, the mean-square density fluctuation.

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

We search for correlations between positions of extragalactic objects and arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) with primary energy $E \ge 40$ EeV as observed by the surface detector array of the Telescope Array (TA) experiment during the first 40 months of operation. We examined several public astronomical object catalogs, including the Veron-Cetty and Veron catalog of active galactic nuclei. We counted the number of TA events correlated with objects in each catalog as a function of three parameters: the maximum angular separation between a TA event and an object, the minimum energy of the events, and the maximum redshift of the objects. We determine combinations of these parameters which maximize the correlations, and calculate the chance probabilities of having the same levels of correlations from an isotropic distribution of UHECR arrival directions. No statistically significant correlations are found when penalties for scanning over the above parameters and for searching in several catalogs are taken into account.

We present the single-epoch black hole mass estimators based on the CIV (1549 A) broad emission line, using the updated sample of the reverberation-mapped AGNs and high-quality UV spectra. By performing multi-component spectral fitting analysis, we measure the CIV line widths (FWHM_CIV) and line dispersion (sigma_CIV) and the continuum luminosity at 1350 A (L_1350) to calibrate the CIV-based mass estimators. By comparing with the Hbeta reverberation-based masses, we provide new mass estimators with the best-fit relationships, i.e., M_BH \propto L_1350 ^ (0.50+-0.07) sigma_CIV ^2 and M_BH \propto L_1350 ^ (0.52+-0.09) FWHM_CIV ^ (0.56+-0.48). The new CIV-based mass estimators show significant mass-dependent systematic difference compared to the estimators commonly used in the literature. Using the published Sloan Digital Sky Survey QSO catalog, we show that the black hole mass of high-redshift QSOs decreases on average by ~0.25 dex if our recipe is adopted.

We measure the spectrum of cosmic rays with energies greater than $10^{18.2}$ eV with the Fluorescence Detectors (FDs) and the Surface Detectors (SDs) of the Telescope Array Experiment using the data taken in our first 2.3-year observation from May 27 2008 to September 7 2010. A hybrid air shower reconstruction technique is employed to improve accuracies in determination of arrival directions and primary energies of cosmic rays using both FD and SD data. The energy spectrum presented here is in agreement with our previously published spectra and the HiRes results.

We present a measurement of the energy spectrum of ultra-high-energy cosmic rays performed by the Telescope Array experiment using monocular observations from its two new FADC-based fluorescence detectors. After a short description of the experiment, we describe the data analysis and event reconstruction procedures. Since the aperture of the experiment must be calculated by Monte Carlo simulation, we describe this calculation and the comparisons of simulated and real data used to verify the validity of the aperture calculation. Finally, we present the energy spectrum calculated from the merged monocular data sets of the two FADC-based detectors, and also the combination of this merged spectrum with an independent, previously published monocular spectrum measurement performed by Telescope Array's third fluorescence detector (Abu-Zayyad {\it et al.}, {Astropart. Phys.} 39 (2012), 109). This combined spectrum corroborates the recently published Telescope Array surface detector spectrum (Abu-Zayyad {\it et al.}, {Astrophys. Journ.} 768 (2013), L1) with independent systematic uncertainties.

The Telescope Array (TA) collaboration has measured the energy spectrum of ultra-high energy cosmic rays with primary energies above 1.6 x 10^(18) eV. This measurement is based upon four years of observation by the surface detector component of TA. The spectrum shows a dip at an energy of 4.6 x 10^(18) eV and a steepening at 5.4 x 10^(19) eV which is consistent with the expectation from the GZK cutoff. We present the results of a technique, new to the analysis of ultra-high energy cosmic ray surface detector data, that involves generating a complete simulation of ultra-high energy cosmic rays striking the TA surface detector. The procedure starts with shower simulations using the CORSIKA Monte Carlo program where we have solved the problems caused by use of the "thinning" approximation. This simulation method allows us to make an accurate calculation of the acceptance of the detector for the energies concerned.

Despite the importance of their size evolution in understanding the dynamical evolution of globular clusters (GCs) of the Milky Way, studies are rare that focus specifically on this issue. Based on the advanced, realistic Fokker-Planck (FP) approach, we predict theoretically the initial size distribution (SD) of the Galactic GCs along with their initial mass function and radial distribution. Over one thousand FP calculations in a wide parameter space have pinpointed the best-fit initial conditions for the SD, mass function, and radial distribution. Our best-fit model shows that the initial SD of the Galactic GCs is of larger dispersion than today's SD, and that typical projected half-light radius of the initial GCs is ~4.6 pc, which is 1.8 times larger than that of the present-day GCs (~2.5 pc). Their large size signifies greater susceptibility to the Galactic tides: the total mass of destroyed GCs reaches 3-5x10^8 M_sun$, several times larger than the previous estimates. Our result challenges a recent view that the Milky Way GCs were born compact on the sub-pc scale, and rather implies that (1) the initial GCs are generally larger than the typical size of the present-day GCs, (2) the initially large GCs mostly shrink and/or disrupt as a result of the galactic tides, and (3) the initially small GCs expand by two-body relaxation, and later shrink by the galactic tides.

We investigate the chemical properties of low-z QSOs, using archival UV spectra obtained with the HST and IUE for a sample of 70 Palomar-Green QSOs at z < 0.5. By utilizing the flux ratio of UV emission lines (i.e., NV /CIV, (SiIV+OIV])/CIV, and NV/HeII) as metallicity indicators, we compare broad-line region (BLR) gas metallicity with AGN properties, i.e., black hole mass, luminosity, and Eddington ratio. We find that BLR metallicity correlates with Eddington ratio while the dependency on black hole mass is much weaker. Although these trends of low-z AGNs appear to be different from those of high-z QSOs, the difference between low-z and high-z samples is partly caused by the limited dynamical range of the samples. We find that metal enrichment at the center of galaxies is closely connected to the accretion activity of black holes and that the scatter of metallicity correlations with black hole mass increases over cosmic time.

We study the anisotropy of Ultra-High Energy Cosmic Ray (UHECR) events collected by the Telescope Array (TA) detector in the first 40 months of operation. Following earlier studies, we examine event sets with energy thresholds of 10 EeV, 40 EeV, and 57 EeV. We find that the distributions of the events in right ascension and declination are compatible with an isotropic distribution in all three sets. We then compare with previously reported clustering of the UHECR events at small angular scales. No significant clustering is found in the TA data. We then check the events with E>57 EeV for correlations with nearby active galactic nuclei. No significant correlation is found. Finally, we examine all three sets for correlations with the large-scale structure of the Universe. We find that the two higher-energy sets are compatible with both an isotropic distribution and the hypothesis that UHECR sources follow the matter distribution of the Universe (the LSS hypothesis), while the event set with E>10 EeV is compatible with isotropy and is not compatible with the LSS hypothesis at 95% CL unless large deflection angles are also assumed. We show that accounting for UHECR deflections in a realistic model of the Galactic magnetic field can make this set compatible with the LSS hypothesis.

The Telescope Array (TA) experiment, located in the western desert of Utah,USA, is designed for observation of extensive air showers from extremely high energy cosmic rays. The experiment has a surface detector array surrounded by three fluorescence detectors to enable simultaneous detection of shower particles at ground level and fluorescence photons along the shower track. The TA surface detectors and fluorescence detectors started full hybrid observation in March, 2008. In this article we describe the design and technical features of the TA surface detector.

The Telescope Array's Middle Drum fluorescence detector was instrumented with telescopes refurbished from the High Resolution Fly's Eye's HiRes-1 site. The data observed by Middle Drum in monocular mode was analyzed via the HiRes-1 profile-constrained geometry reconstruction technique and utilized the same calibration techniques enabling a direct comparison of the energy spectra and energy scales between the two experiments. The spectrum measured using the Middle Drum telescopes is based on a three-year exposure collected between December 16, 2007 and December 16, 2010. The calculated difference between the spectrum of the Middle Drum observations and the published spectrum obtained by the data collected by the HiRes-1 site allows the HiRes-1 energy scale to be transferred to Middle Drum. The HiRes energy scale is applied to the entire Telescope Array by making a comparison between Middle Drum monocular events and hybrid events that triggered both Middle Drum and the Telescope Array's scintillator Ground Array.

Since 2007, the Telescope Array (TA) experiment, based in Utah, USA, has been observing ultra high energy cosmic rays to understand their origins. The experiment involves a surface detector (SD) array and three fluorescence detector (FD) stations. FD stations, installed surrounding the SD array, measure the air fluorescence light emitted from extensive air showers (EASs) for precise determination of their energies and species. The detectors employed at one of the three FD stations were relocated from the High Resolution Fly's Eye experiment. At the other two stations, newly designed detectors were constructed for the TA experiment. An FD consists of a primary mirror and a camera equipped with photomultiplier tubes. To obtain the EAS parameters with high accuracies, understanding the FD optical characteristics is important. In this paper, we report the characteristics and installation of new FDs and the performances of the FD components. The results of the monitored mirror reflectance during the observation time are also described in this report.

Evolution of the mass function (MF) and radial distribution (RD) of the Galactic globular cluster (GC) system is calculated using an advanced and a realistic Fokker-Planck (FP) model that considers dynamical friction, disc/bulge shocks and eccentric cluster orbits. We perform hundreds of FP calculations with different initial cluster conditions, and then search a wide-parameter space for the best-fitting initial GC MF and RD that evolves into the observed present-day Galactic GC MF and RD. By allowing both MF and RD of the initial GC system to vary, which is attempted for the first time in the present Letter, we find that our best-fitting models have a higher peak mass for a lognormal initial MF and a higher cut-off mass for a power-law initial MF than previous estimates, but our initial total masses in GCs, M_{T,i} = 1.5-1.8x10^8 Msun, are comparable to previous results. Significant findings include that our best-fitting lognormal MF shifts downward by 0.35 dex during the period of 13 Gyr, and that our power-law initial MF models well-fit the observed MF and RD only when the initial MF is truncated at >~10^5 Msun. We also find that our results are insensitive to the initial distribution of orbit eccentricity and inclination, but are rather sensitive to the initial concentration of the clusters and to how the initial tidal radius is defined. If the clusters are assumed to be formed at the apocentre while filling the tidal radius there, M_{T,i} can be as high as 6.9x10^8 Msun, which amounts to ~75 per cent of the current mass in the stellar halo.

The Fokker-Planck (FP) model is one of the commonly used methods for studies of the dynamical evolution of dense spherical stellar systems such as globular clusters and galactic nuclei. The FP model is numerically stable in most cases, but we find that it encounters numerical difficulties rather often when the effects of tidal shocks are included in two-dimensional (energy and angular momentum space) version of the FP model or when the initial condition is extreme (e.g., a very large cluster mass and a small cluster radius). To avoid such a problem, we have developed a new integration scheme for a two-dimensional FP equation by adopting an Alternating Direction Implicit (ADI) method given in the Douglas-Rachford split form. We find that our ADI method reduces the computing time by a factor of ~2 compared to the fully implicit method, and resolves problems of numerical instability.

The proton-proton fusion reaction, $pp\to de^+\nu$, is studied in pionless effective field theory (EFT) with di-baryon fields up to next-to leading order. With the aid of the di-baryon fields, the effective range corrections are naturally resummed up to the infinite order and thus the calculation is greatly simplified. Furthermore, the low-energy constant which appears in the axial-current-di-baryon-di-baryon contact vertex is fixed through the ratio of two- and one-body matrix elements which reproduces the tritium lifetime very precisely. As a result we can perform a parameter free calculation for the process. We compare our numerical result with those from the accurate potential model and previous pionless EFT calculations, and find a good agreement within the accuracy better than 1%.

We have measured the periods and light curves of 148 RR Lyrae variables from V=13.5 to 19.7 from the first 100 sq. degrees of the QUEST RR Lyrae survey. Approximately 55% of these stars belong to the clump of stars detected earlier by the Sloan Digital Sky Survey. According to our measurements, this feature has ~10 times the background density of halo stars, spans at least 37.5 deg by 3.5 deg in right ascension and declination (>=30 by >=3 kpc), lies ~50 kpc from the Sun, and has a depth along the line of sight of ~5 kpc (1 sigma). These properties are consistent with the recent models that suggest it is a tidal stream from the Sgr dSph galaxy. The mean period of the type ab variables, 0.58 d, is also consistent. In addition, we have found two smaller over-densities in the halo, one of which may be related to the globular cluster Pal 5.

We describe the discovery circumstances and photometric properties of 2000 EB173, now one of the brightest trans-Neptunian objects (TNOs) with opposition magnitude m_R=18.9 and also one of the largest Plutinos, found with the drift-scanning camera of the QUEST Collaboration, attached to the 1-m Schmidt telescope of the National Observatory of Venezuela. We measure B-V = 0.99 +/- 0.14 and V-R = 0.57 +/- 0.05, a red color observed for many fainter TNOs. At our magnitude limit m_R = 20.1 +/- 0.20, our single detection reveals a sky density of 0.015 (+0.034, -0.012) TNOs per deg^2 (the error bars are 68% confidence limits), consistent with fainter surveys showing a cumulative number proportional to 10^0.5m_R. Assuming an inclination distribution of TNOs with FWHM exceeding 30 deg, it is likely that one hundred to several hundred objects brighter than m_R=20.1 remain to be discovered.