The annual modulation signal observed by the DAMA experiment is a long-standing question in the community of dark matter direct detection. This necessitates an independent verification of its existence using the same detection technique. The COSINE-100 experiment has been operating with 106~kg of low-background NaI(Tl) detectors providing interesting checks on the DAMA signal. However, due to higher backgrounds in the NaI(Tl) crystals used in COSINE-100 relative to those used for DAMA, it was difficult to reach final conclusions. Since the start of COSINE-100 data taking in 2016, we also have initiated a program to develop ultra-pure NaI(Tl) crystals for COSINE-200, the next phase of the experiment. The program includes efforts of raw powder purification, ultra-pure NaI(Tl) crystal growth, and detector assembly techniques. After extensive research and development of NaI(Tl) crystal growth, we have successfully grown a few small-size (0.61$-$0.78 kg) thallium-doped crystals with high radio-purity. A high light yield has been achieved by improvements of our detector assembly technique. Here we report the ultra-pure NaI(Tl) detector developments at the Institute for Basic Science, Korea. The technique developed here will be applied to the production of NaI(Tl) detectors for the COSINE-200 experiment.

NaI(Tl) crystals are used as particle detectors in a variety of rare-event search experiments because of their superb light emission quality. The crystal light yield is generally high, above 10 photoelectrons per keVee (electron-equivalent energy), and its emission spectrum is peaked around 400 nm, which matches well to the sensitive region of bialkali photocathode photomultiplier tubes. However, since NaI(Tl) crystals are hygroscopic, a sophisticated method of encapsulation has to be applied that prevents moisture from chemically attacking the crystal and thereby degrading the emission. In addition, operation with low energy thresholds, which is essential for a number of new phenomenon searches, is usually limited by the crystal light yield; in these cases higher light yields can translate into lower thresholds that improve the experimental sensitivity. Here we describe the development of an encapsulation technique that simplifies the overall design by attaching the photo sensors directly to the crystal so that light losses are minimized. The light yield of a NaI(Tl) crystal encapsulated with this technique was improved by more than 30%, and as many as 22 photoelectrons per keVee have been measured. Consequently, the energy threshold can be lowered and the energy resolution improved. Detectors with this higher light yield are sensitive to events with sub-keVee energies and well suited for low-mass dark matter particle searches and measurements of neutrino-nucleus coherent scattering.

The most significant and practical figure of merit in axion haloscope searches is the scanning rate, because of the unknown axion mass. Even under the best experimental parameters, the only way to improve the figure of merit is to increase the experimentally designed signal to noise ratio in the axion haloscope search analysis procedure. In this paper, we report an improved axion haloscope search analysis, developed using CAPP-8TB. By fully incorporating the known negative correlations induced by the background parametrizations, the figure of merit of the CAPP-8TB haloscope search was increased by about 17\% for the scanning rate, with about an 8\% improvement in the signal to noise ratio. In addition, the physics results can be retrieved directly from the standard Gaussian statistics of the data, at any confidence level.

This paper describes the liquid scintillator veto system for the COSINE-100 dark matter experiment and its performance. The COSINE-100 detector consists of 8 NaI(Tl) crystals immersed in 2,200~L of linear alkylbenzene based liquid scintillator. The liquid scintillator tags between 65 and 75\% of the internal $^{40}$K background in the 2--6~keV energy region. We also describe the background model for the liquid scintillator, which is primarily used to assess its energy calibration and threshold.

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

There is a long standing debate about whether or not the annual modulation signal reported by the DAMA/LIBRA collaboration is induced by Weakly Interacting Massive Particles~(WIMP) in the galaxy's dark matter halo scattering from nuclides in their NaI(Tl) crystal target/detector. This is because regions of WIMP-mass vs. WIMP-nucleon cross-section parameter space that can accommodate the DAMA/LIBRA-phase1 modulation signal in the context of the standard WIMP dark matter galactic halo and isospin-conserving~(canonical), spin-independent~(SI) WIMP-nucleon interactions have been excluded by many of other dark matter search experiments including COSINE-100, which uses the same NaI(Tl) target/detector material. Moreover, the recently released DAMA/LIBRA-phase2 results are inconsistent with an interpretation as WIMP-nuclide scattering via the canonical SI interaction and prefer, instead, isospin-violating or spin-dependent interactions. Dark matter interpretations of the DAMA/LIBRA signal are sensitive to the NaI(Tl) scintillation efficiency for nuclear recoils, which is characterized by so-called quenching factors~(QF), and the QF values used in previous studies differ significantly from recently reported measurements, which may have led to incorrect interpretations of the DAMA/LIBRA signal. In this article, the compatibility of the DAMA/LIBRA and COSINE-100 results, in light of the new QF measurements is examined for different possible types of WIMP-nucleon interactions. The resulting allowed parameter space regions associated with the DAMA/LIBRA signal are explicitly compared with 90\% confidence level upper limits from the initial 59.5~day COSINE-100 exposure. With the newly measured QF values, the allowed 3$\sigma$ regions from the DAMA/LIBRA data are still generally excluded by the COSINE-100 data.

COSINE-100 is a direct detection dark matter search experiment that uses a 106 kg array of eight NaI(Tl) crystals that are kept underground at the Yangyang Underground Laboratory to avoid cosmogenic activation of radioisotopes by cosmic rays. Even though the cosmogenic activity is declining with time, there are still significant background rates from the remnant nuclides. In this paper, we report measurements of cosmogenic isotope contaminations with less than one year half-lives that are based on extrapolations of the time dependent activities of their characteristic energy peaks to activity rates at the time the crystals were deployed underground. For longer-lived $^{109}$Cd ($T_{1/2}=1.6$ y) and $^{22}$Na ($T_{1/2}=2.6$ y), we investigate time correlations of characteristic $\gamma$/X-ray peaks. The inferred sea-level production rates are compared with caluclations based on the ACTIVIA and MENDL-2 model calculations and experimental data. For $^{3}$H, which has a long, 12.3 year half-life, we evaluated the activity levels from the exposure times and determined a cosmogenic activation rate that is consistent with other measurements.

Early type galaxies (ETGs) frequently show emission from warm ionized gas. These Low Ionization Emission Regions (LIERs) were originally attributed to a central, low-luminosity active galactic nuclei. However, the recent discovery of spatially-extended LIER emission suggests ionization by both a central source and an extended component that follows a stellar-like radial distribution. For passively-evolving galaxies with old stellar populations, hot post-Asymptotic Giant Branch (AGB) stars are the only viable extended source of ionizing photons. In this work, we present the first prediction of LIER-like emission from post-AGB stars that is based on fully self-consistent stellar evolution and photoionization models. We show that models where post-AGB stars are the dominant source of ionizing photons reproduce the nebular emission signatures observed in ETGs, including LIER-like emission line ratios in standard optical diagnostic diagrams and H$\alpha$ equivalent widths of order 0.1-3 angstroms. We test the sensitivity of LIER-like emission to the details of post-AGB models, including the mass loss efficiency and convective mixing efficiency, and show that line strengths are relatively insensitive to post-AGB timescale variations. Finally, we examine the UV-optical colors of the models and the stellar populations responsible for the UV-excess observed in some ETGs. We find that allowing as little as 3% of the HB population to be uniformly distributed to very hot temperatures (30,000 K) produces realistic UV colors for old, quiescent ETGs.

Stellar population synthesis (SPS) models have long struggled to reproduce observed optical through near-IR (NIR) spectral energy distributions (SED) of massive quiescent galaxies. We revisit this issue using a novel approach that combines the diagnostic power of full-spectrum fitting with recently updated stellar spectral libraries. First, we perform full-spectrum fitting of continuum-normalized stacked SDSS spectra in bins of velocity dispersion to infer their stellar population properties, such as the elemental abundances and age. Next, we use the resulting best-fit parameters to compute $ugriz$ colors, which are then compared to observed colors of the same galaxies. With this approach we are able to predict the $ugriz$ SEDs of low and high mass galaxies at the $\lesssim 0.03$ mag level in nearly all cases. We find that the full optical through NIR SEDs of quiescent galaxies can be reproduced only when the spectrum is fit with a flexibility that is able to capture the behavior of the entire optical absorption line spectrum. The models include variations in individual elemental abundances, nebular emission lines, and the presence of young stellar components. The successful prediction of the SED shape from continuum-normalized spectra implies that the continuum information is largely contained in the narrow absorption features. These results also imply that attempts to model broadband photometry of quiescent systems will suffer from potentially significant biases if the detailed abundance patterns are not taken into account.

The semi-empirical initial-final mass relation (IFMR) connects spectroscopically analyzed white dwarfs in star clusters to the initial masses of the stars that formed them. Most current stellar evolution models, however, predict that stars will evolve to white dwarfs $\sim$0.1 M$_\odot$ less massive than that found in the IFMR. We first look at how varying theoretical mass-loss rates, third dredge-up efficiencies, and convective-core overshoot may help explain the differences between models and observations. These parameters play an important role at the lowest masses (M$_{\rm initial}$ $<$ 3 M$_\odot$). At higher masses, only convective-core overshoot meaningfully affects white dwarf mass, but alone it likely cannot explain the observed white dwarf masses nor why the IFMR scatter is larger than observational errors predict. These higher masses, however, are also where rotational mixing in main sequence stars begins to create more massive cores, and hence more massive white dwarfs. This rotational mixing also extends a star's lifetime, making faster rotating progenitors appear like less massive stars in their semi-empirical age analysis. Applying the observed range of young B-dwarf rotations to the MIST or SYCLIST rotational models demonstrates a marked improvement in reproducing both the observed IFMR data and its scatter. The incorporation of both rotation and efficient convective-core overshoot significantly improves the match with observations. This work shows that the IFMR provides a valuable observational constraint on how rotation and convective-core overshoot affect the core evolution of a star.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

The near-Earth asteroid 3200 Phaethon (1983 TB) is an attractive object not only from a scientific viewpoint but also because of JAXA's DESTINY+ target. The rotational lightcurve and spin properties were investigated based on the data obtained in the ground-based observation campaign of Phaethon. We aim to refine the lightcurves and shape model of Phaethon using all available lightcurve datasets obtained via optical observation, as well as our time-series observation data from the 2017 apparition. Using eight 1-2-m telescopes and an optical imager, we acquired the optical lightcurves and derived the spin parameters of Phaethon. We applied the lightcurve inversion method and SAGE algorithm to deduce the convex and non-convex shape model and pole orientations. We analysed the optical lightcurve of Phaethon and derived a synodic and a sidereal rotational period of 3.6039 h, with an axis ratio of a/b = 1.07. The ecliptic longitude (lambda) and latitude (beta) of the pole orientation were determined as (308, -52) and (322, -40) via two independent methods. A non-convex model from the SAGE method, which exhibits a concavity feature, is also presented.

We present the initial-final mass relation (IFMR) based on the self-consistent analysis of Sirius B and 79 white dwarfs from 13 star clusters. We have also acquired additional signal on eight white dwarfs previously analyzed in the NGC 2099 cluster field, four of which are consistent with membership. These reobserved white dwarfs have masses ranging from 0.72 to 0.97 M$_\odot$, with initial masses from 3.0 to 3.65 M$_\odot$, where the IFMR has an important change in slope that these new data help to observationally confirm. In total, this directly measured IFMR has small scatter ($\sigma$ = 0.06 M$_\odot$) and spans from progenitors of 0.85 to 7.5 M$_\odot$. Applying two different stellar evolutionary models to infer two different sets of white dwarf progenitor masses shows that when the same model is also used to derive the cluster ages, the resulting IFMR has weak sensitivity to the adopted model at all but the highest initial masses ($>$5.5 M$_\odot$). The non-linearity of the IFMR is also clearly observed with moderate slopes at lower masses (0.08 M$_{\rm final}$/M$_{\rm initial}$) and higher masses (0.11 M$_{\rm final}$/M$_{\rm initial}$) that are broken up by a steep slope (0.19 M$_{\rm final}$/M$_{\rm initial}$) between progenitors from 2.85 to 3.6 M$_\odot$. This IFMR shows total stellar mass loss ranges from 33\% of M$_{\rm initial}$ at 0.83 M$_\odot$ to 83\% of M$_{\rm initial}$ at 7.5 M$_\odot$. Testing this total mass loss for dependence on progenitor metallicity, however, finds no detectable sensitivity across the moderate range of -0.15 $<$ [Fe/H] $<$ +0.15.

Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical simulation of reionization and galaxy formation in the Local Universe, in a volume large enough to model reionization globally but with enough resolving power to follow all atomic-cooling galactic halos in that volume. For every halo identified at a given time, we find the redshift at which the surrounding IGM reionized, along with its instantaneous star formation rate (`SFR') and baryonic gas-to-dark matter ratio ($M_\text{gas}/M_\text{DM}$). The average SFR per halo with $M < 10^9 \text{ M}_\odot$ was steady in regions not yet reionized, but declined sharply following local reionization. For $M > 10^{10} \text{ M}_\odot$, this SFR continued through local reionization, increasing with time, instead. For $10^9 < M < 10^{10} \text{ M}_\odot$, the SFR generally increased modestly through reionization, followed by a modest decline. In general, halo SFRs were higher for regions that reionized earlier. A similar pattern was found for $M_\text{gas}/M_\text{DM}$, which declined sharply following local reionization for $M < 10^9 \text{ M}_\odot$. Local reionization time correlates with local matter overdensity, which determines the local rates of structure formation and ionizing photon consumption. The earliest patches to develop structure and reionize ultimately produced more stars than they needed to finish and maintain their own reionization, exporting their `surplus' starlight to help reionize regions that developed structure later.

The Hyades, Praesepe, and Pleiades are well studied stellar clusters that anchor important secondary stellar age indicators. Recent studies have shown that main sequence turn off-based ages for these clusters may depend on the degree of rotation in the underlying stellar models. Rotation induces structural instabilities that can enhance the chemical mixing of a star, extending its fuel supply. In addition, rotation introduces a modulation of the star's observed magnitude and color due to the effects of gravity darkening. We aim to investigate the extent to which stellar rotation affects the age determination of star clusters. We utilize the MESA stellar evolution code to create models that cover a range of rotation rates corresponding to $\Omega/\Omega_c=0.0$ to $0.6$ in $0.1$ dex steps, allowing the assessment of variations in this dimension. The statistical analysis package, MATCH, is employed to derive ages and metallicities by fitting our MESA models to Tycho $B_T$, $V_T$ and 2MASS $J$, $K_s$ color-magnitude diagrams. We find that the derived ages are relatively insensitive to the effects of rotation. For the Hyades, Praesepe, and Pleiades, we derive ages based on synthetic populations that model a distribution of rotation rates or a fixed rate. Across each case, derived ages tend to agree roughly within errors, near $680$, $590$, and $110-160$ Myr for the Hyades, Praesepe, and Pleiades, respectively. These ages are in agreement with Li depletion boundary-based ages and previous analyses that used non-rotating isochrones. Our methods do not provide a strong constraint on the metallicities of these clusters.

We use the framework developed as part of the MESA Isochrones and Stellar Tracks (MIST) project to assess the utility of several types of observables in jointly measuring the age and 1D stellar model parameters in star clusters. We begin with a pedagogical overview summarizing the effects of stellar model parameters, such as the helium abundance, mass-loss efficiency, and the mixing length parameter, on observational diagnostics including the color-magnitude diagram, mass-radius relation, and surface abundances, amongst others. We find that these parameters and the stellar age influence observables in qualitatively distinctive, degeneracy-breaking ways. To assess the current state of affairs, we use the recent Gaia Data Release 2 (DR2) along with data from the literature to investigate three well-studied old open clusters---NGC6819, M67, NGC6791---as case studies. Although there is no obvious tension between the existing observations and the MIST models for NGC6819, there are interesting discrepancies in the cases of M67 and NGC6791. At this time, parallax zero point uncertainties in Gaia DR2 remain one of the limiting factors in the analysis of these clusters. With a combination of exquisite photometry, parallax distances, and cluster memberships from Gaia at the end of its mission, we anticipate precise and accurate ages for these and other star clusters in the Galaxy.

Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magnetic flux in higher-order multipolar components of the magnetic field. In this work, we present a new model for stellar spin-down that, for the first time, accounts for the stellar surface magnetic field configuration. We show how a magnetic complexity that evolves from complex toward simple configurations as a star spins down can explain the salient features of stellar rotation evolution, including the bimodal distribution of both slow and fast rotators seen in young open clusters.

The effective temperature (Teff) distribution of stellar evolution models along the red giant branch (RGB) is sensitive to a number of parameters including the overall metallicity, elemental abundance patterns, the efficiency of convection, and the treatment of the surface boundary condition. Recently there has been interest in using observational estimates of the RGB Teff to place constraints on the mixing length parameter, a_MLT, and possible variation with metallicity. Here we use 1D MESA stellar evolution models to explore the sensitivity of the RGB Teff to the treatment of the surface boundary condition. We find that different surface boundary conditions can lead to +/- 100 K metallicity-dependent offsets on the RGB relative to one another in spite of the fact that all models can reproduce the properties of the Sun. Moreover, for a given atmosphere T-tau relation, we find that the RGB Teff is also sensitive to the optical depth at which the surface boundary condition is applied in the stellar model. Nearly all models adopt the photosphere as the location of the surface boundary condition but this choice is somewhat arbitrary. We compare our models to stellar parameters derived from the APOGEE-Kepler sample of first ascent red giants and find that systematic uncertainties in the models due to treatment of the surface boundary condition place a limit of ~100 K below which it is not possible to make firm conclusions regarding the fidelity of the current generation of stellar models.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We report on the RXTE detection of a sudden increase in the absorption column density, $N_\mathrm{H}$, during the 2011 May outburst of GX 304-1. The $N_\mathrm{H}$ increased up to ${\sim}16\times 10^{22}$ atoms cm$^{-2}$, which is a factor of 3-4 larger than what is usually measured during the outbursts of GX 304-1 as covered by RXTE. Additionally, an increase in the variability of the hardness ratio as calculated from the energy resolved RXTE-PCA light curves is measured during this time range. We interpret these facts as an occultation event of the neutron star by material in the line of sight. Using a simple 3D model of an inclined and precessing Be disk around the Be type companion, we are able to qualitatively explain the $N_\mathrm{H}$ evolution over time. We are able to constrain the Be-disk density to be on the order of $10^{-11}$ g cm$^{-3}$. Our model strengthens the idea of inclined Be disks as origin of double-peaked outbursts as the derived geometry allows accretion twice per orbit under certain conditions.

We investigate the integrated properties of massive (>10 Msun), rotating, single-star stellar populations for a variety of initial rotation rates (v/vcrit=0.0, 0.2, 0.4, 0.5, and 0.6). We couple the new MESA Isochrone and Stellar Tracks (MIST) models to the Flexible Stellar Population Synthesis (FSPS) package, extending the stellar population synthesis models to include the contributions from very massive stars (>100 Msun), which can be significant in the first ~4 Myr after a starburst. These models predict ionizing luminosities that are consistent with recent observations of young nuclear star clusters. We also construct composite stellar populations assuming a distribution of initial rotation rates. Even in low-metallicity environments where rotation has a significant effect on the evolution of massive stars, we find that stellar population models require a significant contribution from fast-rotating (v/vcrit>0.4) stars in order to sustain the production of ionizing photons beyond a few Myr following a starburst. These results have potentially important implications for cosmic reionization by massive stars and the interpretation of nebular emission lines in high-redshift star-forming galaxies.

We perform smoothed particle hydrodynamics (SPH) simulations of an isolated galaxy with a new treatment for dust formation and destruction. To this aim, we treat dust and metal production self-consistently with star formation and supernova feedback. For dust, we consider a simplified model of grain size distribution by representing the entire range of grain sizes with large and small grains. We include dust production in stellar ejecta, dust destruction by supernova (SN) shocks, grain growth by accretion and coagulation, and grain disruption by shattering. We find that the assumption of fixed dust-to-metal mass ratio becomes no longer valid when the galaxy is older than 0.2 Gyr, at which point the grain growth by accretion starts to contribute to the nonlinear rise of dust-to-gas ratio. As expected in our previous one-zone model, shattering triggers grain growth by accretion since it increases the total surface area of grains. Coagulation becomes significant when the galaxy age is greater than $\sim$ 1 Gyr: at this epoch the abundance of small grains becomes high enough to raise the coagulation rate of small grains. We further compare the radial profiles of dust-to-gas ratio $(\mathcal{D})$ and dust-to-metal ratio $(\mathcal{D}/Z)$ (i.e., depletion) at various ages with observational data. We find that our simulations broadly reproduce the radial gradients of dust-to-gas ratio and depletion. In the early epoch ($\lesssim 0.3$ Gyr), the radial gradient of $\mathcal{D}$ follows the metallicity gradient with $\mathcal{D}/Z$ determined by the dust condensation efficiency in stellar ejecta, while the $\mathcal{D}$ gradient is steeper than the $Z$ gradient at the later epochs because of grain growth by accretion. The framework developed in this paper is applicable to any SPH-based galaxy evolution simulations including cosmological ones.

Unlike spiral galaxies such as the Milky Way, the majority of the stars in massive elliptical galaxies were formed in a short period early in the history of the Universe. The duration of this formation period can be measured using the ratio of magnesium to iron abundance ([Mg/Fe]), which reflects the relative enrichment by core-collapse and type Ia supernovae. For local galaxies, [Mg/Fe] probes the combined formation history of all stars currently in the galaxy, including younger and metal-poor stars that were added during late-time mergers. Therefore, to directly constrain the initial star-formation period, we must study galaxies at earlier epochs. The most distant galaxy for which [Mg/Fe] had previously been measured is at z~1.4, with [Mg/Fe]=0.45(+0.05,-0.19). A slightly earlier epoch (z~1.6) was probed by stacking the spectra of 24 massive quiescent galaxies, yielding an average [Mg/Fe] of 0.31+/-0.12. However, the relatively low S/N of the data and the use of index analysis techniques for both studies resulted in measurement errors that are too large to allow us to form strong conclusions. Deeper spectra at even earlier epochs in combination with analysis techniques based on full spectral fitting are required to precisely measure the abundance pattern shortly after the major star-forming phase (z>2). Here we report a measurement of [Mg/Fe] for a massive quiescent galaxy at z=2.1. With [Mg/Fe]=0.59+/-0.11, this galaxy is the most Mg-enhanced massive galaxy found so far, having twice the Mg enhancement of similar-mass galaxies today. The abundance pattern of the galaxy is consistent with enrichment exclusively by core-collapse supernovae and with a star-formation timescale of 0.1-0.5 Gyr - characteristics that are similar to population II stars in the Milky Way. With an average past SFR of 600-3000 Msol/yr, this galaxy was among the most vigorous star-forming galaxies in the Universe.

Density inhomogeneity in the intergalactic medium (IGM) can boost the recombination rate of ionized gas substantially, affecting the growth of HII regions during reionization. Previous attempts to quantify this effect typically failed to resolve down to the Jeans scale in the pre-ionization IGM, which is important in establishing this effect, along with the hydrodynamical back-reaction of reionization on it. Towards that end, we perform a set of fully-coupled, radiation-hydrodynamics simulations from cosmological initial conditions, extending the mass resolution of previous work to the scale of minihalos. Pre-reionization structure is evolved until a redshift $z_i$ at which the ionizing radiation from external sources arrives to sweep an R-type ionization front supersonically across the volume in a few Myr, until it is trapped on the surfaces of minihalos and converted to D-type, after which the minihalo gas is removed by photoevaporative winds. Small-scale density structures during this time lead to a high ($>$10) clumping factor for ionized gas, which hugely boosts the recombination rate until the structures are disrupted by the hydrodynamic feedback after $\sim 10-100~\rm{Myr}$. For incoming stellar radiation with intensity $J_{21}$ in a $200~h^{-1}~\rm{kpc}$ box with the mean density contrast $\bar\delta$, the number of extra recombinations per H atom, on top of what is expected from homogeneously distributed gas, is given by $0.32[J_{21}]^{0.12}[(1+z_i)/11]^{-1.7}[1+\bar\delta]^{2.5}$. In models in which most of the volume is ionized toward the end of reionization, this can add more than one recombination per H atom to the ionizing photon budget to achieve reionization.

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

We employ high-resolution cosmological zoom-in simulations focusing on a high-sigma peak and an average cosmological field at $z\sim 6-12$, in order to investigate the influence of environment and baryonic feedback on galaxy evolution in the reionization epoch. Strong feedback, e.g., galactic winds, caused by elevated star formation rates (SFRs) is expected to play an important role in this evolution. We compare different outflow prescriptions: (i) constant wind velocity (CW), (ii) variable wind scaling with galaxy properties (VW), and (iii) no outflows (NW). The overdensity leads to accelerated evolution of dark matter and baryonic structures, absent in the "normal" region, and to shallow galaxy stellar mass functions at the low-mass end. Although CW shows little dependence on both environments, the more physically motivated VW model does exhibit this effect. In addition, VW can reproduce the observed specific SFR (sSFR) and the sSFR-stellar mass relation, which CW and NW fail to satisfy simultaneously. Winds also differ substantially in affecting the state of the intergalactic medium (IGM). The difference lies in volume-filling factor of hot, high-metallicity gas which is near unity for CW, while it remains confined in massive filaments for VW, and locked up in galaxies for NW. Such gas is nearly absent in the normal region. Although all wind models suffer from deficiencies, the VW model seems to be promising in correlating the outflow properties to those of host galaxies. Further constraints on the state of the IGM at high-$z$ are needed to separate different wind models.

We report the characterization of a massive (m_p=3.9 +- 1.4 M_jup) microlensing planet (OGLE-2015-BLG-0954Lb) orbiting an M dwarf host (M=0.33 +- 0.12 M_sun) at a distance toward the Galactic bulge of 0.6 (+0.4,-0.2) kpc, which is extremely nearby by microlensing standards. The planet-host projected separation is a_perp ~ 1.2 AU. The characterization was made possible by the wide-field (4 sq. deg.) high cadence (Gamma = 6/hr) monitoring of the Korea Microlensing Telescope Network (KMTNet), which had two of its three telescopes in commissioning operations at the time of the planetary anomaly. The source crossing time t_* = 16 min is among the shortest ever published. The high-cadence, wide-field observations that are the hallmark of KMTNet are the only way to routinely capture such short crossings. High-cadence resolution of short caustic crossings will preferentially lead to mass and distance measurements for the lens. This is because the short crossing time typically implies a nearby lens, which enables the measurement of additional effects (bright lens and/or microlens parallax). When combined with the measured crossing time, these effects can yield planet/host masses and distance.

The results of a search for solar axions from the Korea Invisible Mass Search (KIMS) experiment at the Yangyang Underground Laboratory are presented. Low-energy electron-recoil events would be produced by conversion of solar axions into electrons via the axio-electric effect in CsI(Tl) crystals. Using data from an exposure of 34,596 $\rm kg \cdot days$, we set a 90 \% confidence level upper limit on the axion-electron coupling, $g_{ae}$, of $1.39 \times 10^{-11}$ for an axion mass less than 1 keV/$\rm c^2$. This limit is lower than the indirect solar neutrino bound, and fully excludes QCD axions heavier than 0.48 eV/$\rm c^2$ and 140.9 eV/$\rm c^2$ for the DFSZ and KSVZ models respectively.

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

We report on the mass and distance measurements of two single-lens events from the 2015 \emph{Spitzer} microlensing campaign. With both finite-source effect and microlens parallax measurements, we find that the lens of OGLE-2015-BLG-1268 is very likely a brown dwarf. Assuming that the source star lies behind the same amount of dust as the Bulge red clump, we find the lens is a $45\pm7$ $M_{\rm J}$ brown dwarf at $5.9\pm1.0$ kpc. The lens of of the second event, OGLE-2015-BLG-0763, is a $0.50\pm0.04$ $M_\odot$ star at $6.9\pm1.0$ kpc. We show that the probability to definitively measure the mass of isolated microlenses is dramatically increased once simultaneous ground- and space-based observations are conducted.

Spitzer microlensing parallax observations of OGLE-2015-BLG-1212 decisively breaks a degeneracy between planetary and binary solutions that is somewhat ambiguous when only ground-based data are considered. Only eight viable models survive out of an initial set of 32 local minima in the parameter space. These models clearly indicate that the lens is a stellar binary system possibly located within the bulge of our Galaxy, ruling out the planetary alternative. We argue that several types of discrete degeneracies can be broken via such space-based parallax observations.

We report the discovery of a microlensing exoplanet OGLE-2012-BLG-0563Lb with the planet-star mass ratio ~1 x 10^{-3}. Intensive photometric observations of a high-magnification microlensing event allow us to detect a clear signal of the planet. Although no parallax signal is detected in the light curve, we instead succeed at detecting the flux from the host star in high-resolution JHK'-band images obtained by the Subaru/AO188 and IRCS instruments, allowing us to constrain the absolute physical parameters of the planetary system. With the help of a spectroscopic information about the source star obtained during the high-magnification state by Bensby et al., we find that the lens system is located at 1.3^{+0.6}_{-0.8} kpc from us, and consists of an M dwarf (0.34^{+0.12}_{-0.20} M_sun) orbited by a Saturn-mass planet (0.39^{+0.14}_{-0.23} M_Jup) at the projected separation of 0.74^{+0.26}_{-0.42} AU (close model) or 4.3^{+1.5}_{-2.5} AU (wide model). The probability of contamination in the host star's flux, which would reduce the masses by a factor of up to three, is estimated to be 17%. This possibility can be tested by future high-resolution imaging. We also estimate the (J-Ks) and (H-Ks) colors of the host star, which are marginally consistent with a low metallicity mid-to-early M dwarf, although further observations are required for the metallicity to be conclusive. This is the fifth sub-Jupiter-mass (0.2<m_p/M_Jup<1) microlensing planet around an M dwarf with the mass well constrained. The relatively rich harvest of sub-Jupiters around M dwarfs is contrasted with a possible paucity of ~1--2 Jupiter-mass planets around the same type of star, which can be explained by the planetary formation process in the core-accretion scheme.

We report the detection and mass measurement of a binary lens OGLE-2015-BLG-1285La,b, with the more massive component having $M_1>1.35\,M_\odot$ (80% probability). A main-sequence star in this mass range is ruled out by limits on blue light, meaning that a primary in this mass range must be a neutron star or black hole. The system has a projected separation $r_\perp= 6.1\pm 0.4\,{\rm AU}$ and lies in the Galactic bulge. These measurements are based on the "microlens parallax" effect, i.e., comparing the microlensing light curve as seen from $Spitzer$, which lay at $1.25\,{\rm AU}$ projected from Earth, to the light curves from four ground-based surveys, three in the optical and one in the near infrared. Future adaptive optics imaging of the companion by 30m class telescopes will yield a much more accurate measurement of the primary mass. This discovery both opens the path and defines the challenges to detecting and characterizing black holes and neutron stars in wide binaries, with either dark or luminous companions. In particular, we discuss lessons that can be applied to future $Spitzer$ and $Kepler$ K2 microlensing parallax observations.

We use cosmological adaptive mesh refinement (AMR) code Enzo zoom-in simulations to study the long term evolution of the collapsing gas within dark matter (DM) halos at high redshifts. This direct collapse process is a leading candidate for rapid formation of supermassive black hole (SMBH) seeds at high z. To circumvent the Courant condition at small radii, we have used the sink particle method, and focus on the evolution on scales ~0.01-10 pc. The collapse proceeds in two stages, with the secondary runaway happening within the central 10 pc, and with no detected fragmentation. The sink particles form when the collapsing gas requires additional refinement of the grid size at the highest refinement level. Their mass never exceeds ~10^3 Mo, with the sole exception of the central seed which grows dramatically to ~ 2 x 10^6 Mo in ~2 Myr, confirming the feasibility of this path to the SMBH. The time variability of angular momentum axis in the accreted gas results in the formation of two misaligned disks --- a small inner disk, and a more massive, outer disk which is inclined by ~45^o to the inner disk. The self-gravity of these disks is heavily diluted --- both disks lie within the Roche limit of the central seed. While the inner disk is geometrically thin and weakly asymmetric, the outer disk flares due to turbulent motions as a result of the massive inflow along a pair of penetrating filaments. The geometry of inflow via filaments determines the dominant and secondary Fourier modes in this disk --- these modes have a non-self-gravitational origin. We do not confirm that m=1 is a principal mode that drives the inflow in the presence of a central massive object. While the positions of the disks depend on the scale chosen to break the self-similar collapse, the overall configuration appears to be generic, and is expected to form when the central seed becomes sufficiently massive.

Long period variable stars arise in the final stages of the asymptotic giant branch phase of stellar evolution. They have periods of up to ~1000d and amplitudes that can exceed a factor of three in the I-band flux. These stars pulsate predominantly in their fundamental mode, which is a function of mass and radius, and so the pulsation periods are sensitive to the age of the underlying stellar population. The overall number of long period variables in a population is directly related to their lifetime, which is difficult to predict from first principles because of uncertainties associated with stellar mass-loss and convective mixing. The time variability of these stars has not been previously taken into account when modeling the spectral energy distributions of galaxies. Here we construct time-dependent stellar population models that include the effects of long period variable stars, and report the ubiquitous detection of this expected `pixel shimmer' in the massive metal-rich galaxy M87. The pixel light curves display a variety of behaviors, including linearly rising and falling curves, semi-periodic curves, and sudden increases or decreases in the flux level. The observed variation of 0.1-1% is very well matched to the predictions of our models. The data provide a strong and novel constraint on the properties of variable stars in an old and metal-rich stellar population, and we infer that the lifetime of long period variables in M87 is shorter by approximately 30% compared to predictions from the latest stellar evolution models.

For all exoplanet candidates, the reliability of a claimed detection needs to be assessed through a careful study of systematic errors in the data to minimize the false positives rate. We present a method to investigate such systematics in microlensing datasets using the microlensing event OGLE-2013-BLG-0446 as a case study. The event was observed from multiple sites around the world and its high magnification (A_{max} \sim 3000) allowed us to investigate the effects of terrestrial and annual parallax. Real-time modeling of the event while it was still ongoing suggested the presence of an extremely low-mass companion (\sim 3M_\oplus ) to the lensing star, leading to substantial follow-up coverage of the light curve. We test and compare different models for the light curve and conclude that the data do not favour the planetary interpretation when systematic errors are taken into account.

We report the detection of a Cold Neptune m_planet=21+/-2MEarth orbiting a 0.38MSol M dwarf lying 2.5-3.3 kpc toward the Galactic center as part of a campaign combining ground-based and Spitzer observations to measure the Galactic distribution of planets. This is the first time that the complex real-time protocols described by Yee et al. (2015), which aim to maximize planet sensitivity while maintaining sample integrity, have been carried out in practice. Multiple survey and follow-up teams successfully combined their efforts within the framework of these protocols to detect this planet. This is the second planet in the Spitzer Galactic distribution sample. Both are in the near-to-mid disk and clearly not in the Galactic bulge.

Motivated by the recent, serendipitous discovery of the densest known galaxy, M60-UCD1, we present two initial findings from a follow-up search, using the Sloan Digital Sky Survey, Subaru/Suprime-Cam and Hubble Space Telescope imaging, and SOuthern Astrophysical Research (SOAR)/Goodman spectroscopy. The first object discovered, M59-UCD3, has a similar size to M60-UCD1 (half-light radius of r_h ~ 20 pc) but is 40% more luminous (M_V ~ -14.6), making it the new densest-known galaxy. The second, M85-HCC1, has a size like a typical globular cluster (GC; r_h ~ 1.8 pc) but is much more luminous (M_V ~ -12.5). This hypercompact cluster is by far the densest confirmed free-floating stellar system, and is equivalent to the densest known nuclear star clusters. From spectroscopy, we find that both objects are relatively young (~9 Gyr and ~3 Gyr, respectively), with metal-abundances that resemble those of galaxy centers. Their host galaxies show clear signs of large-scale disturbances, and we conclude that these dense objects are the remnant nuclei of recently accreted galaxies. M59-UCD3 is an ideal target for follow-up with high-resolution imaging and spectroscopy to search for an overweight central supermassive black hole as was discovered in M60-UCD1. These findings also emphasize the potential value of ultra-compact dwarfs and massive GCs as tracers of the assembly histories of galaxies.

We report the discovery of an extrasolar planet, KMT-2015-1b, that was detected using the microlensing technique. The planetary lensing event was observed by KMTNet survey that has commenced in 2015. With dense coverage by using network of globally distributed telescopes equipped with very wide-field cameras, the short planetary signal is clearly detected and precisely characterized. We find that KMT-2015-1b is a giant planet orbiting a low-mass M-dwarf host star. The planet has a mass about twice that of Jupiter and it is located beyond the snow line of the host star. With the improvement of existing surveys and the advent of new surveys, future microlensing planet samples will include planets not only in greatly increased number but also in a wide spectrum of hosts and planets, helping us to have a better and comprehensive understanding about the formation and evolution of planets.

The co-planarity of solar-system planets led Kant to suggest that they formed from an accretion disk, and the discovery of hundreds of such disks around young stars as well as hundreds of co-planar planetary systems by the Kepler satellite demonstrate that this formation mechanism is extremely widespread. Many moons in the solar system, such as the Galilean moons of Jupiter, also formed out of the accretion disks that coalesced into the giant planets. We report here the discovery of an intermediate system OGLE-2013-BLG-0723LB/Bb composed of a Venus-mass planet orbiting a brown dwarf, which may be viewed either as a scaled down version of a planet plus star or as a scaled up version of a moon plus planet orbiting a star. The latter analogy can be further extended since they orbit in the potential of a larger, stellar body. For ice-rock companions formed in the outer parts of accretion disks, like Uranus and Callisto, the scaled masses and separations of the three types of systems are similar, leading us to suggest that formation processes of companions within accretion disks around stars, brown dwarfs, and planets are similar.

We study the early stage of the formation of seed supermassive black holes via direct collapse in dark matter (DM) halos, in the cosmological context. We perform high-resolution zoom-in simulations of such collapse at high-$z$. Using the adaptive mesh refinement code ENZO, we resolve the formation and growth of a DM halo, until its virial temperature reaches $\sim 10^4$K, atomic cooling turns on, and collapse ensues. We demonstrate that direct collapse proceeds in two stages, although they are not well separated. The first stage is triggered by the onset of atomic cooling, and leads to rapidly increasing accretion rate with radius, from $\dot M\sim 0.1\,M_\odot {\rm yr^{-1}}$ at the halo virial radius to few $M_\odot \,{\rm yr^{-1}}$, around the scale radius $R_{\rm s}\sim 30$pc of the NFW DM density profile. The second stage of the collapse commences when the gas density takes precedence over the DM density. This is associated with the gas decoupling from the DM gravitational potential. The ensuing collapse approximates that of an isothermal sphere with $\dot M ( r )\sim $const. We confirm that the gas loses its angular momentum through non-axisymmetric perturbations and gravitational torques, to overcome the centrifugal barrier. During the course of the collapse, this angular momentum transfer process happens on nearly all spatial scales, and the angular momentum vector of the gas varies with position and time. Collapsing gas also exhibits supersonic turbulent motions which suppress gas fragmentation, and are characterized by density PDF consisting of a lognormal part and a high-density power law tail.

We reanalyze microlensing events in the published list of anomalous events that were observed from the OGLE lensing survey conducted during 2004-2008 period. In order to check the existence of possible degenerate solutions and extract extra information, we conduct analyses based on combined data from other survey and follow-up observation and consider higher-order effects. Among the analyzed events, we present analyses of 8 events for which either new solutions are identified or additional information is obtained. We find that the previous binary-source interpretations of 5 events are better interpreted by binary-lens models. These events include OGLE-2006-BLG-238, OGLE-2007-BLG-159, OGLE-2007-BLG-491, OGLE-2008-BLG-143, and OGLE-2008-BLG-210. With additional data covering caustic crossings, we detect finite-source effects for 6 events including OGLE-2006-BLG-215, OGLE-2006-BLG-238, OGLE-2006-BLG-450, OGLE-2008-BLG-143, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. Among them, we are able to measure the Einstein radii of 3 events for which multi-band data are available. These events are OGLE-2006-BLG-238, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. For OGLE-2008-BLG-143, we detect higher-order effect induced by the changes of the observer's position caused by the orbital motion of the Earth around the Sun. In addition, we present degenerate solutions resulting from the known close/wide or ecliptic degeneracy. Finally, we note that the masses of the binary companions of the lenses of OGLE-2006-BLG-450 and OGLE-2008-BLG-210 are in the brown-dwarf regime.

We report the discovery of a Jupiter-mass planet orbiting an M-dwarf star that gave rise to the microlensing event OGLE-2011-BLG-0265. Such a system is very rare among known planetary systems and thus the discovery is important for theoretical studies of planetary formation and evolution. High-cadence temporal coverage of the planetary signal combined with extended observations throughout the event allows us to accurately model the observed light curve. The final microlensing solution remains, however, degenerate yielding two possible configurations of the planet and the host star. In the case of the preferred solution, the mass of the planet is $M_{\rm p} = 0.9\pm 0.3\ M_{\rm J}$, and the planet is orbiting a star with a mass $M = 0.22\pm 0.06\ M_\odot$. The second possible configuration (2$\sigma$ away) consists of a planet with $M_{\rm p}=0.6\pm 0.3\ M_{\rm J}$ and host star with $M=0.14\pm 0.06\ M_\odot$. The system is located in the Galactic disk 3 -- 4 kpc towards the Galactic bulge. In both cases, with an orbit size of 1.5 -- 2.0 AU, the planet is a "cold Jupiter" -- located well beyond the "snow line" of the host star. Currently available data make the secure selection of the correct solution difficult, but there are prospects for lifting the degeneracy with additional follow-up observations in the future, when the lens and source star separate.

We have studied channeling effects in a Cesium Iodide (CsI) crystal that is similar in composition to the ones being used in a search for Weakly Interacting Massive Particles (WIMPs) dark matter candidates, and measured its energy-dependent quenching factor, the relative scintillation yield for electron and nuclear recoils. The experimental results are reproduced with a GEANT4 simulation that includes a model of the scintillation efficiency as a function of electronic stopping power. We present the measured and simulated quenching factors and the estimated effects of channeling.

High-magnification gravitational microlensing events provide an important channel of detecting planetary systems with multiple giants located at their birth places. In order to investigate the potential existence of additional planets, we reanalyze the light curves of the eight high-magnification microlensing events for each of which a single planet was previously detected. The analyzed events include OGLE-2005-BLG-071, OGLE-2005-BLG-169, MOA-2007-BLG-400, MOA-2008-BLG-310, MOA-2009-BLG-319, MOA-2009-BLG-387, MOA-2010-BLG-477, and MOA-2011-BLG-293. We find that including an additional planet improves fits with $\Delta\chi^2 < 80$ for seven out of eight analyzed events. For MOA-2009-BLG-319, the improvement is relatively big with $\Delta\chi^2 \sim 143$. From inspection of the fits, we find that the improvement of the fits is attributed to systematics in data. Although no clear evidence of additional planets is found, it is still possible to constrain the existence of additional planets in the parameter space. For this purpose, we construct exclusion diagrams showing the confidence levels excluding the existence of an additional planet as a function of its separation and mass ratio. We also present the exclusion ranges of additional planets with 90\% confidence level for Jupiter, Saturn, and Uranus-mass planets.

We report the first mass and distance measurement of a caustic-crossing binary system OGLE-2014-BLG-1050L using the space-based microlens parallax method. \emph{Spitzer} captured the second caustic-crossing of the event, which occurred $\sim$10 days before that seen from Earth. Due to the coincidence that the source-lens relative motion was almost parallel to the direction of the binary-lens axis, the four-fold degeneracy, which was known before only to occur in single-lens events, persists in this case, leading to either a lower-mass (0.2 $M_\odot$ and 0.07 $M_\odot$) binary at $\sim$1.1 kpc or a higher-mass (0.9 $M_\odot$ and 0.35 $M_\odot$) binary at $\sim$3.5 kpc. However, the latter solution is strongly preferred for reasons including blending and lensing probability. OGLE-2014-BLG-1050L demonstrates the power of microlens parallax in probing stellar and substellar binaries.

We constructed a neutron calibration facility based on a 300-mCi Am-Be source in conjunction with a search for weakly interacting massive particle candidates for dark matter. The facility is used to study the response of CsI(Tl) crystals to nuclear recoils induced by neutrons from the Am-Be source and comparing them with the response to electron recoils produced by Compton scattering of 662-keV $\gamma$-rays from a $^{137}$Cs source. The measured results on pulse shape discrimination (PSD) between nuclear- and electron-recoil events are quantified in terms of quality factors. A comparison with similar result from a neutron reactor demonstrate the feasibility of performing calibrations of PSD measurements using neutrons from a Am-Be source.

We present the analysis of the gravitational microlensing event OGLE-2013-BLG-0102. The light curve of the event is characterized by a strong short-term anomaly superposed on a smoothly varying lensing curve with a moderate magnification $A_{\rm max}\sim 1.5$. It is found that the event was produced by a binary lens with a mass ratio between the components of $q = 0.13$ and the anomaly was caused by the passage of the source trajectory over a caustic located away from the barycenter of the binary. From the analysis of the effects on the light curve due to the finite size of the source and the parallactic motion of the Earth, the physical parameters of the lens system are determined. The measured masses of the lens components are $M_{1} = 0.096 \pm 0.013~M_{\odot}$ and $M_{2} = 0.012 \pm 0.002~M_{\odot}$, which correspond to near the hydrogen-burning and deuterium-burning mass limits, respectively. The distance to the lens is $3.04 \pm 0.31~{\rm kpc}$ and the projected separation between the lens components is $0.80 \pm 0.08~{\rm AU}$.

We present a search for low-mass ($\leq 20 GeV/c^{2}$) weakly interacting massive particles(WIMPs), strong candidates of dark matter particles,using the low-background CsI(Tl) detector array of the Korea Invisible Mass Search (KIMS) experiment. With a total data exposure of 24,324.3kg$\cdot$days,we search for WIMP interaction signals produced by nuclei recoiling from WIMP-nuclear elastic scattering with visible energies between 2 and 4keV. The observed energy distribution of candidate events is consistent with null signals, and upper limits of the WIMP-proton spin-independent interaction are set with a 90% confidence level. The observed limit rejects most of the low mass region of parameter space favored by the DAMA annual modulation signal.

We present results from modeling the optical spectra of a large sample of quiescent galaxies between 0.1 < z < 0.7 from the Sloan Digital Sky Survey (SDSS) and the AGN and Galaxy Evolution Survey (AGES). We examine how the stellar ages and abundance patterns of galaxies evolve over time as a function of stellar mass from 10^{9.6}-10^{11.8} Msun. Galaxy spectra are stacked in bins of mass and redshift, and modeled over a wavelength range from 4000 A to 5500 A. Full spectrum stellar population synthesis modeling provides estimates of the age and the abundances of the elements Fe, Mg, C, N, and Ca. We find negligible evolution in elemental abundances at fixed stellar mass over roughly 7 Gyr of cosmic time. In addition, the increase in stellar ages with time for massive galaxies is consistent with passive evolution since z = 0.7. Taken together, these results favor a scenario in which the inner ~ 0.3-3 R_e of massive quiescent galaxies have been passively evolving over the last half of cosmic time. Interestingly, the derived stellar ages are considerably younger than the age of the universe at all epochs, consistent with an equivalent single-burst star formation epoch of z < 1.5. These young stellar population ages coupled with the existence of massive quiescent galaxies at z > 1 indicate the inhomogeneous nature of the z < 0.7 quiescent population. The data also permit the addition of newly-quenched galaxies at masses below ~10^{10.5} Msun at z < 0.7. Additionally, we analyze very deep Keck DEIMOS spectra of the two brightest quiescent galaxies in a cluster at z = 0.83. There is tentative evidence that these galaxies are older than their counterparts in low-density environments. In the Appendix, we demonstrate that our full spectrum modeling technique allows for accurate and reliable modeling of galaxy spectra to low S/N (~20 A^{-1}) and/or low spectral resolution (R ~ 500).

We developed a cryogenic system on a rotating table that achieves sub-Kelvin conditions. The cryogenic system consists of a helium sorption cooler and a pulse tube cooler in a cryostat mounted on a rotating table. Two rotary-joint connectors for electricity and helium gas circulation enable the coolers to be operated and maintained with ease. We performed cool-down tests under a condition of continuous rotation at 20 rpm. We obtained a temperature of 0.23 K with a holding time of more than 24 hours, thus complying with catalog specifications. We monitored the system's performance for four weeks; two weeks with and without rotation. A few-percent difference in conditions was observed between these two states. Most applications can tolerate such a slight difference. The technology developed is useful for various scientific applications requiring sub-Kelvin conditions on rotating platforms.

We have explored prevailing modes of galaxy growth for redshifts z ~ 6-14, comparing substantially overdense and normal regions of the universe, using high-resolution zoom-in cosmological simulations. Such rare overdense regions have been projected to host high-z quasars. We demonstrate that galaxies in such environments grow predominantly by a smooth accretion from cosmological filaments which dominates the mass input from major, intermediate and minor mergers. We find that by z ~6, the accumulated galaxy mass fraction from mergers falls short by a factor of 10 of the cumulative accretion mass for galaxies in the overdense regions, and by a factor of 5 in the normal environments. Moreover, the rate of the stellar mass input from mergers also lies below that of an in-situ star formation (SF) rate. The fraction of stellar masses in galaxies contributed by mergers in overdense regions is ~12%, and ~33% in the normal regions, at these redshifts. Our median SF rates for ~few X 10^9 Mo galaxies agrees well with the recently estimated rates for z ~ 7 galaxies from Spitzer's SURF-UP survey. Finally, we find that the main difference between the normal and overdense regions lies in the amplified growth of massive galaxies in massive dark matter halos. This leads to the formation of >= 10^{10} Mo galaxies due to the ~100-fold increase in mass during the above time period. Such galaxies are basically absent in the normal regions at these redshifts.

We detect a cold, terrestrial planet in a binary-star system using gravitational microlensing. The planet has low mass (2 Earth masses) and lies projected at $a_{\perp,ph}$ ~ 0.8 astronomical units (AU) from its host star, similar to the Earth-Sun distance. However, the planet temperature is much lower, T<60 Kelvin, because the host star is only 0.10--0.15 solar masses and therefore more than 400 times less luminous than the Sun. The host is itself orbiting a slightly more massive companion with projected separation $a_{\perp,ch}=$10--15 AU. Straightforward modification of current microlensing search strategies could increase their sensitivity to planets in binary systems. With more detections, such binary-star/planetary systems could place constraints on models of planet formation and evolution. This detection is consistent with such systems being very common.

We report the discovery of MOA-2013-BLG-220Lb, which has a super-Jupiter mass ratio $q=3.01\pm 0.02\times 10^{-3}$ relative to its host. The proper motion, $\mu=12.5\pm 1\, {\rm mas}\,{\rm yr}^{-1}$, is one of the highest for microlensing planets yet discovered, implying that it will be possible to separately resolve the host within $\sim 7$ years. Two separate lines of evidence imply that the planet and host are in the Galactic disk. The planet could have been detected and characterized purely with follow-up data, which has important implications for microlensing surveys, both current and into the LSST era.

Characterizing a microlensing planet is done from modeling an observed lensing light curve. In this process, it is often confronted that solutions of different lensing parameters result in similar light curves, causing difficulties in uniquely interpreting the lens system, and thus understanding the causes of different types of degeneracy is important. In this work, we show that incomplete coverage of a planetary perturbation can result in degenerate solutions even for events where the planetary signal is detected with a high level of statistical significance. We demonstrate the degeneracy for an actually observed event OGLE-2012-BLG-0455/MOA-2012-BLG-206. The peak of this high-magnification event $(A_{\rm max}\sim400)$ exhibits very strong deviation from a point-lens model with $\Delta\chi^{2}\gtrsim4000$ for data sets with a total number of measurement 6963. From detailed modeling of the light curve, we find that the deviation can be explained by four distinct solutions, i.e., two very different sets of solutions, each with a two-fold degeneracy. While the two-fold (so-called "close/wide") degeneracy is well-understood, the degeneracy between the radically different solutions is not previously known. The model light curves of this degeneracy differ substantially in the parts that were not covered by observation, indicating that the degeneracy is caused by the incomplete coverage of the perturbation. It is expected that the frequency of the degeneracy introduced in this work will be greatly reduced with the improvement of the current lensing survey and follow-up experiments and the advent of new surveys.

The mass of the lenses giving rise to Galactic microlensing events can be constrained by measuring the relative lens-source proper motion and lens flux. The flux of the lens can be separated from that of the source, companions to the source, and unrelated nearby stars with high-resolution images taken when the lens and source are spatially resolved. For typical ground-based adaptive optics (AO) or space-based observations, this requires either inordinately long time baselines or high relative proper motions. We provide a list of microlensing events toward the Galactic Bulge with high relative lens-source proper motion that are therefore good candidates for constraining the lens mass with future high-resolution imaging. We investigate all events from 2004 -- 2013 that display detectable finite-source effects, a feature that allows us to measure the proper motion. In total, we present 20 events with mu >~ 8 mas/yr. Of these, 14 were culled from previous analyses while 6 are new, including OGLE-2004-BLG-368, MOA-2005-BLG-36, OGLE-2012-BLG-0211, OGLE-2012-BLG-0456, MOA-2012-BLG-532, and MOA-2013-BLG-029. In <~12 years the lens and source of each event will be sufficiently separated for ground-based telescopes with AO systems or space telescopes to resolve each component and further characterize the lens system. Furthermore, for the most recent events, comparison of the lens flux estimates from images taken immediately to those estimated from images taken when the lens and source are resolved can be used to empirically check the robustness of the single-epoch method currently being used to estimate lens masses for many events.

One important strength of the microlensing method in detecting extrasolar planets is its high sensitivity to low-mass planets. However, it is often believed that microlensing detections of Earth-mass planets from ground-based observation would be difficult due to the limit set by finite-source effects. This view comes from the previous estimation of the planet detection probability based on the fractional deviation of planetary signals, but proper probability estimation requires to additionally consider the source brightness, which is directly related to the photometric precision. In this paper, we reevaluate the feasibility of low-mass planet detections considering photometric precision for different populations of source stars. From this, it is found that contribution of the improved photometric precision to the planetary signal of a giant-source event is big enough to compensate the decrease of the magnification excess caused by finite-source effects. As a result, we find that giant-source events are suitable targets for Earth-mass planet detections with significantly higher detection probability than events involved with source stars of smaller radii and predict that Earth-mass planets would be detected by prospective high-cadence surveys.

Recent observations have indicated the existence of dust in high-redshift galaxies, however, the dust properties in them are still unknown. Here we present theoretical constraints on dust properties in Lyman break galaxies (LBGs) at z=3 by post-processing a cosmological smoothed particle hydrodynamics simulation with radiative transfer calculations. We calculate the dust extinction in 2800 dark matter halos using the metallicity information of individual gas particles in our simulation. We use only bright galaxies with rest-frame UV magnitude M_1700 < -20 mag, and study the dust size, dust-to-metal mass ratio, and dust composition. From the comparison of calculated color excess between B and V-band (i.e., E(B-V)) and the observations, we constrain the typical dust size, and show that the best-fitting dust grain size is ~ 0.05 micron, which is consistent with the results of theoretical dust models for Type-II supernova. Our simulation with the dust extinction effect can naturally reproduce the observed rest-frame UV luminosity function of LBGs at z=3 without assuming an ad hoc constant extinction value. In addition, in order to reproduce the observed mean E(B-V), we find that the dust-to-metal mass ratio needs to be similar to that of the local galaxies, and that the graphite dust is dominant or at least occupy half of dust mass.

It has been shown observationally that star formation (SF) correlates tightly with the presence of molecular hydrogen (H$_2$). Therefore it would be important to investigate its implication on galaxy formation in a cosmological context. In the present work, we track the H$_2$ mass fraction within our cosmological smoothed particle hydrodynamics (SPH) code GADGET-3 using an equilibrium analytic model by Krumholz et al. This model allows us to regulate the star formation in our simulation by the local abundance of H$_2$ rather than the total cold gas density, and naturally introduce the dependence of star formation on metallicity. We investigate implications of the equilibrium H$_2$-based SF model on galaxy population properties, such as the stellar-to-halo mass ratio (SHMR), baryon fraction, cosmic star formation rate density (SFRD), galaxy specific SFR, galaxy stellar mass functions (GSMF), and Kennicutt-Schmidt (KS) relationship. The advantage of our work over the previous ones is having a large sample of simulated galaxies in a cosmological volume from high-redshift to $z=0$. We find that low-mass halos with M$_{DM}<10^{10.5}$ M$_\odot$ are less efficient in producing stars in the H$_2$-based SF model at $z\geq6$, which brings the simulations to a better agreement with observational estimates of SHMR and GSMF at the low-mass end. This is particularly evident by a reduction in the number of low-mass galaxies at M$_\star\leq10^{8}$ M$_\odot$ in the GSMF. The overall SFRD is also reduced at high-$z$ in the H$_2$ run, which results in slightly higher SFRD at low-redshift due to more abundant gas available for star formation at later times. This new H$_2$ model is able to reproduce the empirical KS relationship at $z=0$ naturally without the need for setting its normalization by hand, and overall it seems to have more advantages than the previous pressure-based SF model.

We present a detailed analysis of survey and follow-up observations of microlensing event OGLE-2012-BLG-0406 based on data obtained from 10 different observatories. Intensive coverage of the lightcurve, especially the perturbation part, allowed us to accurately measure the parallax effect and lens orbital motion. Combining our measurement of the lens parallax with the angular Einstein radius determined from finite-source effects, we estimate the physical parameters of the lens system. We find that the event was caused by a $2.73\pm 0.43\ M_{\rm J}$ planet orbiting a $0.44\pm 0.07\ M_{\odot}$ early M-type star. The distance to the lens is $4.97\pm 0.29$\ kpc and the projected separation between the host star and its planet at the time of the event is $3.45\pm 0.26$ AU. We find that the additional coverage provided by follow-up observations, especially during the planetary perturbation, leads to a more accurate determination of the physical parameters of the lens.

In the field of radiowave detection, enlarging the receiver aperture to enhance the amount of light detected is essential for greater scientific achievements. One challenge in using radio transmittable apertures is keeping the detectors cool. This is because transparency to thermal radiation above the radio frequency range increases the thermal load. A technology that maintains cold conditions while allowing larger apertures has been long-awaited. We propose radio-transparent multi-layer insulation (RT-MLI), composed from a set of stacked insulating layers. The insulator is transparent to radio frequencies, but not transparent to infrared radiation. The basic idea for cooling is similar to conventional multi-layer insulation. It leads to a reduction in thermal radiation while maintaining a uniform surface temperature. The advantage of this technique over other filter types is that no thermal links are required. As insulator material, we used foamed polystyrene; its low index of refraction makes an anti-reflection coating unnecessary. We measured the basic performance of RT-MLI to confirm that thermal loads are lowered with more layers. We also confirmed that our RT-MLI has high transmittance to radiowaves, but blocks infrared radiation. For example, RT-MLI with 12 layers has a transmittance greater than 95% (lower than 1%) below 200 GHz (above 4 THz). We demonstrated its effects in a system with absorptive-type filters, where aperture diameters were 200 mm. Low temperatures were successfully maintained for the filters. We conclude that this technology significantly enhances the cooling of radiowave receivers, and is particularly suitable for large-aperture systems. This technology is expected to be applicable to various fields, including radio astronomy, geo-environmental assessment, and radar systems.

We analyze the planetary microlensing event MOA-2010-BLG-328. The best fit yields host and planetary masses of Mh = 0.11+/-0.01 M_{sun} and Mp = 9.2+/-2.2M_Earth, corresponding to a very late M dwarf and sub-Neptune-mass planet, respectively. The system lies at DL = 0.81 +/- 0.10 kpc with projected separation r = 0.92 +/- 0.16 AU. Because of the host's a-priori-unlikely close distance, as well as the unusual nature of the system, we consider the possibility that the microlens parallax signal, which determines the host mass and distance, is actually due to xallarap (source orbital motion) that is being misinterpreted as parallax. We show a result that favors the parallax solution, even given its close host distance. We show that future high-resolution astrometric measurements could decisively resolve the remaining ambiguity of these solutions.

A planetary microlensing signal is generally characterized by a short-term perturbation to the standard single lensing light curve. A subset of binary-source events can produce perturbations that mimic planetary signals, thereby introducing an ambiguity between the planetary and binary-source interpretations. In this paper, we present analysis of the microlensing event MOA-2012-BLG-486, for which the light curve exhibits a short-lived perturbation. Routine modeling not considering data taken in different passbands yields a best-fit planetary model that is slightly preferred over the best-fit binary-source model. However, when allowed for a change in the color during the perturbation, we find that the binary-source model yields a significantly better fit and thus the degeneracy is clearly resolved. This event not only signifies the importance of considering various interpretations of short-term anomalies, but also demonstrates the importance of multi-band data for checking the possibility of false-positive planetary signals.

Observations of accretion disks around young brown dwarfs have led to the speculation that they may form planetary systems similar to normal stars. While there have been several detections of planetary-mass objects around brown dwarfs (2MASS 1207-3932 and 2MASS 0441-2301), these companions have relatively large mass ratios and projected separations, suggesting that they formed in a manner analogous to stellar binaries. We present the discovery of a planetary-mass object orbiting a field brown dwarf via gravitational microlensing, OGLE-2012-BLG-0358Lb. The system is a low secondary/primary mass ratio (0.080 +- 0.001), relatively tightly-separated (~0.87 AU) binary composed of a planetary-mass object with 1.9 +- 0.2 Jupiter masses orbiting a brown dwarf with a mass 0.022 M_Sun. The relatively small mass ratio and separation suggest that the companion may have formed in a protoplanetary disk around the brown dwarf host, in a manner analogous to planets.

We use numerical simulations to explore whether direct collapse can lead to the formation of SMBH seeds at high-z. We follow the evolution of gas within DM halos of 2 x 10^8 Mo and 1 kpc. We adopt cosmological density profiles and j-distributions. Our goal is to understand how the collapsing flow overcomes the centrifugal barrier and whether it is subject to fragmentation. We find that the collapse leads either to a central runaway or to off-center fragmentation. A disk-like configuration is formed inside the centrifugal barrier. For more cuspy DM distribution, the gas collapses more and experiences a bar-like perturbation and a central runaway. We have followed this inflow down to ~10^{-4} pc. The flow remains isothermal and the specific angular momentum is efficiently transferred by gravitational torques in a cascade of nested bars. This cascade supports a self-similar, disk-like collapse. In the collapsing phase, virial supersonic turbulence develops and fragmentation is damped. For larger initial DM cores the timescales become longer. In models with more organized initial rotation, a torus forms and appears to be supported by turbulent motions. The evolution depends on the competition between two timescales, corresponding to the onset of the central runaway and off-center fragmentation. For less organized rotation, the torus is greatly weakened, the central accretion timescale is shortened, and off-center fragmentation is suppressed --- triggering the central runaway even in previously `stable' models. The resulting SMBH masses lie in the range 2 x 10^4 Mo - 2 x 10^6 Mo, much higher than for Population III remnants. We argue that the above upper limit appears to be more realistic mass. Corollaries of this model include a possible correlation between SMBH and DM halo masses, and similarity between the SMBH and halo mass functions, at time of formation.

Gravitational microlensing events produced by lenses composed of binary masses are important because they provide a major channel to determine physical parameters of lenses. In this work, we analyze the light curves of two binary-lens events OGLE-2006-BLG-277 and OGLE-2012-BLG-0031 for which the light curves exhibit strong deviations from standard models. From modeling considering various second-order effects, we find that the deviations are mostly explained by the effect of the lens orbital motion. We also find that lens parallax effects can mimic orbital effects to some extent. This implies that modeling light curves of binary-lens events not considering orbital effects can result in lens parallaxes that are substantially different from actual values and thus wrong determinations of physical lens parameters. This demonstrates the importance of routine consideration of orbital effects in interpreting light curves of binary-lens events. It is found that the lens of OGLE-2006-BLG-277 is a binary composed of a low-mass star and a brown dwarf companion.

We developed a system that continuously maintains a cryocooler for long periods on a rotating table. A cryostat that holds the cryocooler is set on the table. A compressor is located on the ground and supplies high-purity (> 99.999%) and high-pressure (1.7 MPa) helium gas and electricity to the cryocooler. The operation of the cryocooler and other instruments requires the development of interface components between the ground and rotating table. A combination of access holes at the center of the table and two rotary joints allows simultaneous circulation of electricity and helium gas. The developed system provides two innovative functions under the rotating condition; cooling from room temperature and the maintenance of a cold condition for long periods. We have confirmed these abilities as well as temperature stability under a condition of continuous rotation at 20 revolutions per minute. The developed system can be applied in various fields; e.g., in tests of Lorentz invariance, searches for axion, radio astronomy and cosmology, and application of radar systems. In particular, there is a plan to use this system for a radio telescope observing cosmic microwave background radiation.

Although many models have been proposed, the physical mechanisms responsible for the formation of low-mass brown dwarfs are poorly understood. The multiplicity properties and minimum mass of the brown-dwarf mass function provide critical empirical diagnostics of these mechanisms. We present the discovery via gravitational microlensing of two very low-mass, very tight binary systems. These binaries have directly and precisely measured total system masses of 0.025 Msun and 0.034 Msun, and projected separations of 0.31 AU and 0.19 AU, making them the lowest-mass and tightest field brown-dwarf binaries known. The discovery of a population of such binaries indicates that brown dwarf binaries can robustly form at least down to masses of ~0.02 Msun. Future microlensing surveys will measure a mass-selected sample of brown-dwarf binary systems, which can then be directly compared to similar samples of stellar binaries.

We present the analysis of the gravitational microlensing event OGLE-2011-BLG-0251. This anomalous event was observed by several survey and follow-up collaborations conducting microlensing observations towards the Galactic Bulge. Based on detailed modelling of the observed light curve, we find that the lens is composed of two masses with a mass ratio q=1.9 x 10^-3. Thanks to our detection of higher-order effects on the light curve due to the Earth's orbital motion and the finite size of source, we are able to measure the mass and distance to the lens unambiguously. We find that the lens is made up of a planet of mass 0.53 +- 0.21,M_Jup orbiting an M dwarf host star with a mass of 0.26 +- 0.11 M_Sun. The planetary system is located at a distance of 2.57 +- 0.61 kpc towards the Galactic Centre. The projected separation of the planet from its host star is d=1.408 +- 0.019, in units of the Einstein radius, which corresponds to 2.72 +- 0.75 AU in physical units. We also identified a competitive model with similar planet and host star masses, but with a smaller orbital radius of 1.50 +- 0.50 AU. The planet is therefore located beyond the snow line of its host star, which we estimate to be around 1-1.5 AU.

We present 248 precise Doppler measurements of Barnard's Star (Gl 699), the second nearest star system to Earth, obtained from Lick and Keck Observatories during 25 years between 1987 and 2012. The early precision was 20 \ms{} but was 2 \ms{} during the last 8 years, constituting the most extensive and sensitive search for Doppler signatures of planets around this stellar neighbor. We carefully analyze the 136 Keck radial velocities spanning 8 years by first applying a periodogram analysis to search for nearly circular orbits. We find no significant periodic Doppler signals with amplitudes above $\sim$2 \ms{}, setting firm upper limits on the minimum mass (\msini) of any planets with orbital periods from 0.1 to 1000 days. Using a Monte Carlo analysis for circular orbits, we determine that planetary companions to Barnard's Star with masses above 2 \mearth{} and periods below 10 days would have been detected. Planets with periods up to 2 years and masses above 10 \mearth{} (0.03 \mjup) are also ruled out. A similar analysis allowing for eccentric orbits yields comparable mass limits. The habitable zone of Barnard's Star appears to be devoid of roughly Earth-mass planets or larger, save for face-on orbits. Previous claims of planets around the star by van de Kamp are strongly refuted. The radial velocity of Barnard's Star increases with time at $4.515\pm0.002$ \msy{}, consistent with the predicted geometrical effect, secular acceleration, that exchanges transverse for radial components of velocity.

We construct several Milky Way-like galaxy models containing a gas halo (as well as gaseous and stellar disks, a dark matter halo, and a stellar bulge) following either an isothermal or an NFW density profile with varying mass and initial spin. In addition, galactic winds associated with star formation are tested in some of the simulations. We evolve these isolated galaxy models using the GADGET-3 $N$-body/hydrodynamic simulation code, paying particular attention to the effects of the gas halo on the evolution. We find that the evolution of the models is strongly affected by the adopted gas halo component. The model without a gas halo shows an increasing star formation rate (SFR) at the beginning of the simulation for some hundreds of millions of years and then a continuously decreasing rate to the end of the run at 3 Gyr. On the other hand, the SFRs in the models with a gas halo emerge to be either relatively flat throughout the simulations or increasing over a gigayear and then decreasing to the end. The models with the more centrally concentrated NFW gas halo show overall higher SFRs than those with the isothermal gas halo of the equal mass. The gas accretion from the halo onto the disk also occurs more in the models with the NFW gas halo, however, this is shown to take place mostly in the inner part of the disk and not to contribute significantly to the star formation unless the gas halo has very high density at the central part. The rotation of a gas halo is found to make SFR lower in the model. The SFRs in the runs including galactic winds are found to be lower than the same runs but without winds. We conclude that the effects of a hot gaseous halo on the evolution of galaxies are generally too significant to be simply ignored, and expect that more hydrodynamical processes in galaxies could be understood through numerical simulations employing both gas disk and gas halo components.

Microlensing can provide an important tool to study binaries, especially those composed of faint or dark objects. However, accurate analysis of binary-lens light curves is often hampered by the well-known degeneracy between close (s<1) and wide (s>1) binaries, which can be very severe due to an intrinsic symmetry in the lens equation. Here s is the normalized projected binary separation. In this paper, we propose a method that can resolve the close/wide degeneracy using the effect of a lens orbital motion on lensing light curves. The method is based on the fact that the orbital effect tends to be important for close binaries while it is negligible for wide binaries. We demonstrate the usefulness of the method by applying it to an actually observed binary-lens event MOA-2011-BLG-040/OGLE-2011-BLG-0001, which suffers from severe close/wide degeneracy. From this, we are able to uniquely specify that the lens is composed of K and M-type dwarfs located at ~3.5 kpc from the Earth.

We present an analysis of the anomalous microlensing event, MOA-2010-BLG-073, announced by the Microlensing Observations in Astrophysics survey on 2010-03-18. This event was remarkable because the source was previously known to be photometrically variable. Analyzing the pre-event source lightcurve, we demonstrate that it is an irregular variable over time scales >200d. Its dereddened color, $(V-I)_{S,0}$, is 1.221$\pm$0.051mag and from our lens model we derive a source radius of 14.7$\pm$1.3 $R_{\odot}$, suggesting that it is a red giant star. We initially explored a number of purely microlensing models for the event but found a residual gradient in the data taken prior to and after the event. This is likely to be due to the variability of the source rather than part of the lensing event, so we incorporated a slope parameter in our model in order to derive the true parameters of the lensing system. We find that the lensing system has a mass ratio of q=0.0654$\pm$0.0006. The Einstein crossing time of the event, $T_{\rm{E}}=44.3$\pm$0.1d, was sufficiently long that the lightcurve exhibited parallax effects. In addition, the source trajectory relative to the large caustic structure allowed the orbital motion of the lens system to be detected. Combining the parallax with the Einstein radius, we were able to derive the distance to the lens, $D_L$=2.8$\pm$0.4kpc, and the masses of the lensing objects. The primary of the lens is an M-dwarf with $M_{L,p}$=0.16$\pm0.03M_{\odot}$ while the companion has $M_{L,s}$=11.0$\pm2.0M_{\rm{J}}$ putting it in the boundary zone between planets and brown dwarfs.

We report the discovery of a planetary system from observation of the high-magnification microlensing event OGLE-2012-BLG-0026. The lensing light curve exhibits a complex central perturbation with multiple features. We find that the perturbation was produced by two planets located near the Einstein ring of the planet host star. We identify 4 possible solutions resulting from the well-known close/wide degeneracy. By measuring both the lens parallax and the Einstein radius, we estimate the physical parameters of the planetary system. According to the best-fit model, the two planet masses are ~0.11 M_Jupiter and 0.68 M_Jupiter and they are orbiting a G-type main sequence star with a mass ~0.82 M_Sun. The projected separations of the individual planets are beyond the snow line in all four solutions, being ~3.8 AU and 4.6 AU in the best-fit solution. The deprojected separations are both individually larger and possibly reversed in order. This is the second multi-planet system with both planets beyond the snow line discovered by microlensing. This is the only such a system (other than the Solar System) with measured planet masses without sin(i) degeneracy. The planetary system is located at a distance 4.1 kpc from the Earth toward the Galactic center. It is very likely that extra light from stars other than the lensed star comes from the lens itself. If this is correct, it will be possible to obtain detailed information about the planet-host star from follow-up observation.

Brown dwarfs are important objects because they may provide a missing link between stars and planets, two populations that have dramatically different formation history. In this paper, we present the candidate binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered during 2004 - 2011 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-2011-BLG-149, and MOA-201-BLG-278/OGLE-2011-BLG-012N. Among them, we are able to confirm that the companions of the lenses of MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149 are brown dwarfs by determining the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The measured mass of the brown dwarf companions are (0.02 +/- 0.01) M_Sun and (0.019 +/- 0.002) M_Sun for MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149, respectively, and both companions are orbiting low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of lensing events with well covered light curves increases with new generation searches.

We study the effect of local stellar radiation and UVB on the physical properties of DLAs and LLSs at z=3 using cosmological SPH simulations. We post-process our simulations with the ART code for radiative transfer of local stellar radiation and UVB. We find that the DLA and LLS cross sections are significantly reduced by the UVB, whereas the local stellar radiation does not affect them very much except in the low-mass halos. This is because clumpy high-density clouds near young star clusters effectively absorb most of the ionizing photons from young stars. We also find that the UVB model with a simple density threshold for self-shielding effect can reproduce the observed column density distribution function of DLAs and LLSs very well, and we validate this model by direct radiative transfer calculations of stellar radiation and UVB with high angular resolution. We show that, with a self-shielding treatment, the DLAs have an extended distribution around star-forming regions typically on ~ 10-30 kpc scales, and LLSs are surrounding DLAs on ~ 30-60 kpc scales. Our simulations suggest that the median properties of DLA host haloes are: Mh = 2.4*10^10 Msun, SFR = 0.3 Msun/yr, M* = 2.4*10^8 Msun, and Z/Zsun = 0.1. About 30 per cent of DLAs are hosted by haloes having SFR = 1 - 20 Msun/yr, which is the typical SFR range for LBGs. More than half of DLAs are hosted by the LBGs that are fainter than the current observational limit. Our results suggest that fractional contribution to LLSs from lower mass haloes is greater than for DLAs. Therefore the median values of LLS host haloes are somewhat lower with Mh = 9.6*10^9 Msun, SFR = 0.06 Msun/yr, M* = 6.5*10^7 Msun and Z/Zsun = 0.08. About 80 per cent of total LLS cross section are hosted by haloes with SFR < 1 Msun/yr, hence most LLSs are associated with low-mass halos with faint LBGs below the current detection limit.

In searches for rare phenomena such as elastic scattering of dark matter particles or neutrinoless double beta decay, alpha decays of Rn222 progenies attached to the surfaces of the detection material have been identified as a serious source of background. In measurements with CsI(Tl) scintillator crystals, we demonstrate that alpha decays of surface contaminants produce fast signals with a characteristic mean-time distribution that is distinct from those of neutron- and gamma-induced events.

We examine the duty cycle and the history of star formation (SFH) for high-redshift galaxies at z>=6 using cosmological hydrodynamic simulations. We find that, even though individual galaxies have bursty SFH, the averaged SFH between z~15 to z=6 can be characterized well by either an exponentially increasing functional form with characteristic time-scales of 70 Myr to 200 Myr for galaxies with stellar masses Ms~10^6 Msun to >10^10 Msun respectively, or by a simple power-law form which exhibits a similar mass dependent time-scales. Using the SFH of individual galaxies, we measure the duty cycle of star formation (DC_SFH); i.e., the fraction of time a galaxy of a particular mass spends above a star formation rate (SFR) threshold which would make it observable to the Hubble Space Telescope (HST) during a given epoch. We also examine the fraction of galaxies at a given redshift that are brighter than a rest-frame UV magnitude (Muv ~ -18), which is sufficient enough to make them observable (DC_Muv). We find that both DC_SFH and DC_Muv make a sharp transition from zero (for galaxies with Ms <= 10^7 Msun) to unity (for Ms > 10^9 Msun). The measured duty cycle is also manifested in the intrinsic scatter in the Ms-SFR relationship (~ 1 dex) and Ms-Muv relationship (\Delta Muv ~ +-1 mag). We provide analytic fits to the DC as a function of Ms using a sigmoid function, which can be used to correct for catalogue incompleteness. We consider the effects of duty cycle to the observational estimate of galaxy stellar mass functions (GSMF) and the star formation rate density (SFRD), and find that it results in a much shallower low-mass end slopes of the GSMF and a reduction of >~ 70% of our intrinsic SFRD, making our simulation results more compatible with observational estimates.

Despite astrophysical importance of binary star systems, detections are limited to those located in small ranges of separations, distances, and masses and thus it is necessary to use a variety of observational techniques for a complete view of stellar multiplicity across a broad range of physical parameters. In this paper, we report the detections and measurements of 2 binaries discovered from observations of microlensing events MOA-2011-BLG-090 and OGLE-2011-BLG-0417. Determinations of the binary masses are possible by simultaneously measuring the Einstein radius and the lens parallax. The measured masses of the binary components are 0.43 $M_{\odot}$ and 0.39 $M_{\odot}$ for MOA-2011-BLG-090 and 0.57 $M_{\odot}$ and 0.17 $M_{\odot}$ for OGLE-2011-BLG-0417 and thus both lens components of MOA-2011-BLG-090 and one component of OGLE-2011-BLG-0417 are M dwarfs, demonstrating the usefulness of microlensing in detecting binaries composed of low-mass components. From modeling of the light curves considering full Keplerian motion of the lens, we also measure the orbital parameters of the binaries. The blended light of OGLE-2011-BLG-0417 comes very likely from the lens itself, making it possible to check the microlensing orbital solution by follow-up radial-velocity observation. For both events, the caustic-crossing parts of the light curves, which are critical for determining the physical lens parameters, were resolved by high-cadence survey observations and thus it is expected that the number of microlensing binaries with measured physical parameters will increase in the future.

New limits are presented on the cross section for Weakly Interacting Massive Particle (WIMP) nucleon scattering in the KIMS CsI(T) detector array at the Yangyang Underground Laboratory. The exposure used for these results is 24524.3 kg\cdotdays. Nuclei recoiling from WIMP interactions are identified by a pulse shape discrimination method. A low energy background due to alpha emitters on the crystal surfaces is identified and taken into account in the analysis. The detected numbers of nuclear recoils are consistent with zero and 90% confidence level upper limits on the WIMP interaction rates are set for electron equivalent energies from 3 keV to 11 keV. The 90% upper limit of NR event rate for 3.6-5.8 keV corresponding to 2-4 keV in NaI(T) is 0.0098 counts/kg/keV/day which is below the annual modulation amplitude reported by DAMA. This is incompatible with interpretations that enhance the modulation amplitude such as inelastic dark matter models. We establish the most stringent cross section limits on spin-dependent WIMP-proton elastic scattering for the WIMP masses greater than 20 GeV/c2.

High-magnification microlensing events provide an important channel to detect planets. Perturbations near the peak of a high-magnification event can be produced either by a planet or a binary companion. It is known that central perturbations induced by both types of companions can be generally distinguished due to the basically different magnification pattern around caustics. In this paper, we present a case of central perturbations for which it is difficult to distinguish the planetary and binary interpretations. The peak of a lensing light curve affected by this perturbation appears to be blunt and flat. For a planetary case, this perturbation occurs when the source trajectory passes the negative perturbation region behind the back end of an arrowhead-shaped central caustic. For a binary case, a similar perturbation occurs for a source trajectory passing through the negative perturbation region between two cusps of an astroid-shaped caustic. We demonstrate the degeneracy for 2 high-magnification events of OGLE-2011-BLG-0526 and OGLE-2011-BLG-0950/MOA-2011-BLG-336. For OGLE-2011-BLG-0526, the $\chi^2$ difference between the planetary and binary model is $\sim$ 3, implying that the degeneracy is very severe. For OGLE-2011-BLG-0950/MOA-2011-BLG-336, the stellar binary model is formally excluded with $\Delta \chi^2 \sim$ 105 and the planetary model is preferred. However, it is difficult to claim a planet discovery because systematic residuals of data from the planetary model are larger than the difference between the planetary and binary models. Considering that 2 events observed during a single season suffer from such a degeneracy, it is expected that central perturbations experiencing this type of degeneracy is common.

We present the analysis of the light curves of 9 high-magnification single-lens gravitational microlensing events with lenses passing over source stars, including OGLE-2004-BLG-254, MOA-2007-BLG-176, MOA-2007-BLG-233/OGLE-2007-BLG-302, MOA-2009-BLG-174, MOA-2010-BLG-436, MOA-2011-BLG-093, MOA-2011-BLG-274, OGLE-2011-BLG-0990/MOA-2011-BLG-300, and OGLE-2011-BLG-1101/MOA-2011-BLG-325. For all events, we measure the linear limb-darkening coefficients of the surface brightness profile of source stars by measuring the deviation of the light curves near the peak affected by the finite-source effect. For 7 events, we measure the Einstein radii and the lens-source relative proper motions. Among them, 5 events are found to have Einstein radii less than 0.2 mas, making the lenses candidates of very low-mass stars or brown dwarfs. For MOA-2011-BLG-274, especially, the small Einstein radius of $\theta_{\rm E}\sim 0.08$ mas combined with the short time scale of $t_{\rm E}\sim 2.7$ days suggests the possibility that the lens is a free-floating planet. For MOA-2009-BLG-174, we measure the lens parallax and thus uniquely determine the physical parameters of the lens. We also find that the measured lens mass of $\sim 0.84\ M_\odot$ is consistent with that of a star blended with the source, suggesting that the blend is likely to be the lens. Although we find planetary signals for none of events, we provide exclusion diagrams showing the confidence levels excluding the existence of a planet as a function of the separation and mass ratio.

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness of the light curve appears to be compatible both with binary-lens and binary-source models, including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated real-time detection of weak microlensing anomalies with immediate feedback is feasible, efficient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that are not easily resolved from photometric light curves, 3) a modelling approach that finds all features of parameter space rather than just the `favourite model' is required, and 4) the data quality is most crucial, where systematics can be confused with real features, in particular small higher-order effects such as orbital motion signatures. It moreover becomes apparent that events with weak signatures are a silver mine for statistical studies, although not easy to exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source microlensing events might hide here.

Microlensing can provide a useful tool to probe binary distributions down to low-mass limits of binary companions. In this paper, we analyze the light curves of 8 binary lensing events detected through the channel of high-magnification events during the seasons from 2007 to 2010. The perturbations, which are confined near the peak of the light curves, can be easily distinguished from the central perturbations caused by planets. However, the degeneracy between close and wide binary solutions cannot be resolved with a $3\sigma$ confidence level for 3 events, implying that the degeneracy would be an important obstacle in studying binary distributions. The dependence of the degeneracy on the lensing parameters is consistent with a theoretic prediction that the degeneracy becomes severe as the binary separation and the mass ratio deviate from the values of resonant caustics. The measured mass ratio of the event OGLE-2008-BLG-510/MOA-2008-BLG-369 is $q\sim 0.1$, making the companion of the lens a strong brown-dwarf candidate.

We present the results of a numerical study comparing photometric and physical properties of simulated z=6-9 galaxies to the observations taken by the WFC3 instrument aboard the Hubble Space Telescope. Using cosmological hydrodynamical simulations we find good agreement with observations in color-color space at all studied redshifts. We also find good agreement between observations and our Schechter luminosity function fit in the observable range, Muv<= -18, provided that a moderate dust extinction effect exists for massive galaxies. However beyond what currently can be observed, simulations predict a very large number of low-mass galaxies and evolving steep faint-end slopes from alpha_L = -2.15 at z=6 to alpha_L = -2.64 at z=9, with a dependence of |alpha_L| \propto (1+z)^0.59. During the same epoch, the normalization phi* increases and the characteristic magnitude Muv* becomes moderately brighter with decreasing redshift. We find similar trends for galaxy stellar mass function with evolving low-mass end slope from alpha_M = - 2.26 at z=6 to alpha_M = -2.87 at z=9, with a dependence of |alpha_M| \propto (1+z)^0.65. Together with our recent result on the high escape fraction of ionizing photons for low-mass galaxies, our results suggest that the low-mass galaxies are important contributor of ionizing photons for the reionisation of the Universe at z>=6.

There are mainly two different approaches to measure the cosmic star formation history: direct star formation rate density (SFRD) and stellar mass density rhostar as functions of redshift. Compilations of current observations seem to show a disparity in the two quantities, in the sense that the integral of SFRD is higher than the observed rhostar (after considering gas recycling). Using cosmological smoothed particle hydrodynamics simulations based on the concordance Lambda cold dark matter model, we show that the two quantities become more consistent with each other when we consider the observed galaxy mass limit. The comparison between simulations and (dust corrected) observed cosmic SFRD shows a good agreement, while the observed rhostar is significantly lower than the simulation results. This can be reconciled if the current high-$z$ galaxy surveys are missing faint low-mass galaxies due to their flux limit. Our simulated GSMFs have steep low-mass end slopes of alpha < -2 at z>3, and when these numerous low-mass galaxies are included, the total rhostar matches with the integral of SFRD.

Using a high-resolution cosmological numerical simulation, we have analyzed the evolution of galaxies at z~10 in a highly overdense region of the universe. These objects could represent the high redshift galaxies recently observed by the Hubble's WFC3, and be as well possible precursors of QSOs at z~6-7. To overcome the sampling and resolution problems in cosmological simulations, we have used the Constrained Realizations method. Our main result for z~10 shows the region of 3.5h^{-1}Mpc radius in comoving coordinates completely dominated by disk galaxies in the total mass range of >=10^9h^{-1}Mo. We have verified that the gaseous and stellar disks we identify are robust morphological features, capable of surviving the ongoing merger process at these redshifts. Below this mass range, we find a sharp decline in the disk fraction to negligible numbers. At this redshift, the disks appear to be gas-rich and the dark matter halos baryon-rich, by a factor of ~2-3 above the average fraction of baryons in the universe. The prevalence of disk galaxies in the high density peaks during the epoch of reionization is contrary to the morphology-density trend observed at low redshifts.

Prompted by various analyses of long (Type II) GRB rates and their relationship to the cosmic star-formation history, metallicity and luminosity function evolution, we systematically analyze these effects with a Monte Carlo code. We test various cosmic star-formation history models including analytical and empirical models as well as those derived from cosmological simulations. We also explore expressions for metallicity enhancement of the GRB rate with redshift, as presented in the literature, and discuss improvements to these analytic expressions from the point of view of galactic evolution. These are also compared to cosmological simulations on metal enrichment. Additionally we explore possible evolutionary effects of the GRB rate and luminosity function with redshift. The simulated results are tested with the observed Swift sample including the L, z, and peak flux (log N-log P) distributions. The observational data imply that an increase in the GRB rate is necessary to account for the observations at high redshift, although the form of this enhancement is unclear. A rate increase due to lower metallicity at higher redshift may not be the singular cause and is subject to a variety of uncertainties. Alternatively, evolution of the GRB luminosity function break with redshift shows promise as a possible alternative.

The Seyfert 1 galaxy Zw 229-015 is among the brightest active galaxies being monitored by the Kepler mission. In order to determine the black hole mass in Zw 229-015 from H-beta reverberation mapping, we have carried out nightly observations with the Kast Spectrograph at the Lick 3m telescope during the dark runs from June through December 2010, obtaining 54 spectroscopic observations in total. We have also obtained nightly V-band imaging with the Katzman Automatic Imaging Telescope at Lick Observatory and with the 0.9m telescope at the Brigham Young University West Mountain Observatory over the same period. We detect strong variability in the source, which exhibited more than a factor of 2 change in broad H-beta flux. From cross-correlation measurements, we find that the H-beta light curve has a rest-frame lag of 3.86(+0.69,-0.90) days with respect to the V-band continuum variations. We also measure reverberation lags for H-alpha and H-gamma and find an upper limit to the H-delta lag. Combining the H-beta lag measurement with a broad H-beta width of sigma = 1590+/-47 km/s measured from the root-mean-square variability spectrum, we obtain a virial estimate of M_BH = 1.00(-0.24,+0.19)*10^7 solar masses for the black hole in Zw 229-015. As a Kepler target, Zw 229-015 will eventually have one of the highest-quality optical light curves ever measured for any active galaxy, and the black hole mass determined from reverberation mapping will serve as a benchmark for testing relationships between black hole mass and continuum variability characteristics in active galactic nuclei.

We study the impact of ultraviolet background (UVB) radiation field and the local stellar radiation on the H_I column density distribution f(N_HI) of damped Ly-alpha systems (DLAs) and sub-DLAs at z=3 using cosmological smoothed particle hydrodynamics simulations. We find that, in the previous simulations with an optically thin approximation, the UVB was sinking into the H_I cloud too deeply, and therefore we underestimated the f(N_HI) at 19 < log(N_HI) < 21.2 compared to the observations. If the UVB is shut off in the high-density regions with n_gas > 6 x 10^{-3} cm^{-3}, then we reproduce the observed f(N_HI) at z=3 very well. We also investigate the effect of local stellar radiation by post-processing our simulation with a radiative transfer code, and find that the local stellar radiation does not change the f(N_HI) very much. Our results show that the shape of f(N_HI) is determined primarily by the UVB with a much weaker effect by the local stellar radiation and that the optically thin approximation often used in cosmological simulation is inadequate to properly treat the ionization structure of neutral gas in and out of DLAs. Our result also indicates that the DLA gas is closely related to the transition region from optically-thick neutral gas to optically-thin ionized gas within dark matter halos.

We study the relationship between the metallicity of gamma-ray burst (GRB) progenitors and the probability distribution function (PDF) of GRB host galaxies as a function of luminosity using cosmological hydrodynamic simulations of galaxy formation. We impose a maximum limit to the gas metallicity in which GRBs can occur, and examine how the predicted luminosity PDF of GRB host galaxies changes in the simulation. We perform the Kolmogorov-Smirnov test, and show that the result from our simulation agrees with the observed luminosity PDF of core-collapse supernovae (SNe) host galaxies when we assume that the core-collapse SNe trace star formation. When we assume that GRBs occur only in a low-metallicity environment with $Z\lesssim 0.1 \Zsun$, GRBs occur in lower luminosity galaxies, and the simulated luminosity PDF becomes quantitatively consistent with the observed luminosity PDF. The observational bias against the host galaxies of optically dark GRBs owing to dust extinction may be another reason for the lower luminosities of GRB host galaxies, but the observed luminosity PDF of GRB host galaxies cannot be reproduced solely by the dust bias in our simulation.

Combing the three-dimensional radiative transfer (RT) calculation and cosmological SPH simulations, we study the escape fraction of ionizing photons (f_esc) of high-redshift galaxies at z=3-6. Our simulations cover the halo mass range of M_h = 10^9 - 10^12 M_sun. We postprocess several hundred simulated galaxies with the Authentic Radiative Transfer (ART) code to study the halo mass dependence of f_esc. In this paper, we restrict ourselves to the transfer of stellar radiation from local stellar population in each dark matter halo. We find that the average f_esc steeply decreases as the halo mass increases, with a large scatter for the lower mass haloes. The low mass haloes with M_h ~ 10^9 M_sun have large values of f_esc (with an average of ~ 0.4), whereas the massive haloes with M_h ~ 10^11 M_sun show small values of f_esc (with an average of ~ 0.07). This is because in our simulations, the massive haloes show more clumpy structure in gas distribution, and star-forming regions are embedded inside these clumps, making it more difficult for the ionizing photons to escape. On the other hand, in low mass haloes, there are often conical regions of highly ionized gas due to the shifted location of young star clusters from the center of dark matter halo, which allows the ionizing photons to escape more easily than in the high-mass haloes. By counting the number of escaped ionizing photons, we show that the star-forming galaxies can ionize the intergalactic medium at z=3-6. The main contributor to the ionizing photons is the haloes with M_h < 10^10 M_sun owing to their high f_esc. The large dispersion in f_esc suggests that there may be various sizes of H{\sc ii} bubbles around the haloes even with the same mass in the early stages of reionization. We also examine the effect of UV background radiation field on f_esc using simple, four different treatment of UV background.

We develop a new ``Multicomponent and Variable Velocity'' (MVV) galactic outflow model for cosmological smoothed particle hydrodynamic (SPH) simulations. The MVV wind model reflects the fact that the wind material can arise from different phases in the interstellar medium (ISM), and the mass-loading factor in the MVV model is a function of galaxy stellar mass. We find that the simulation with the MVV outflow has the following characteristics: (i) the intergalactic medium (IGM) is hardly heated up, and the mean IGM temperature is almost the same as in the no-wind run; (ii) it has lower cosmic star formation rates (SFRs) compared to the no-wind run, but higher SFRs than the constant velocity wind run; (iii) it roughly agrees with the observed IGM metallicity, and roughly follows the observed evolution of Omega(Civ); (iv) the lower mass galaxies have larger mass-loading factors, and the low-mass end of galaxy stellar mass function is flatter than in the previous simulations. Therefore, the MVV outflow model mildly alleviates the problem of too steep galaxy stellar mass function seen in the previous SPH simulations. In summary, the new MVV outflow model shows reasonable agreement with observations, and gives better results than the constant velocity wind model.

We compute the infrared (IR) emission from high-redshift galaxies in cosmological smoothed particle hydrodynamics simulations by coupling the output of the simulation with the population synthesis code `GRASIL' by Silva et al. Based on the stellar mass, metallicity and formation time of each star particle, we estimate the full spectral energy distribution of each star particle from ultraviolet to IR, and compute the luminosity function of simulated galaxies in the Spitzer broadband filters for direct comparison with the available Spitzer observations.

We investigate the effects of the change of cosmological parameters and star formation (SF) models on the cosmic SF history using cosmological smoothed particle hydrodynamics (SPH) simulations based on the cold dark matter (CDM) model. We vary the cosmological parameters within 1-sigma error from the WMAP best-fit parameters, and find that such changes in cosmological parameters mostly affect the amplitude of the cosmic SF history. At high redshift (hereafter high-z), the star formation rate (SFR) is sensitive to the cosmological parameters that control the small-scale power of the primordial power spectrum, while the cosmic matter content becomes important at lower redshifts. We also test two new SF models: 1) the `Pressure' model based on the work by Schaye & Dalla Vecchia (2008), and 2) the `Blitz' model that takes the effect of molecular hydrogen formation into account, based on the work by Blitz & Rosolowsky (2006). Compared to the previous conventional SF model, the Pressure model reduces the SFR in low-density regions and shows better agreement with the observations of the Kennicutt-Schmidt law. This model also suppresses the early star formation and shifts the peak of the cosmic SF history toward lower redshift, more consistently with the recent observational estimates of cosmic SFR density. The simulations with the new SF model also predict lower global stellar mass densities at high-z, larger populations of low-mass galaxies and a higher gas fraction in high-z galaxies. Our results suggest that there is room left in the model uncertainties to reconcile the discrepancy that was found between the theory and observations of cosmic SF history and stellar mass density. Nevertheless, our simulations still predict higher stellar mass densities than most of the observational estimates.

We investigate the physical mechanisms of tidal heating and satellite disruption in cold dark matter host haloes using N-body simulations based on cosmological initial conditions. We show the importance of resonant shocks and resonant torques with the host halo to satellite heating. A resonant shock (torque) couples the radial (tangential) motion of a satellite in its orbit to its phase space. For a satellite on a circular orbit, an ILR-like resonance dominates the heating and this heating results in continuous satellite mass loss. We estimate the requirements for simulations to achieve these dynamics using perturbation theory. Both resonant shocks and resonant torques affect satellites on eccentric orbits. We demonstrate that satellite mass loss is an outside-in process in energy space; a satellite's stars and gas are thus protected by their own halo against tidal stripping. We simulate the evolution of a halo similar to the Large Magellanic Cloud (LMC) in our Galactic dark matter halo and conclude that the LMC stars have not yet been stripped. Finally, we present a simple algorithm for estimating the evolution of satellite mass that includes both shock heating and resonant torques.

We present the results of a numerical study on the effects of metal enrichment and metal cooling on galaxy formation and cosmic star formation (SF) history using cosmological hydrodynamic simulations. We find following differences in the simulation with metal cooling when compared to the run without it: (1) the cosmic star formation rate (SFR) is enhanced by about 50 & 20% at z=1 & 3, respectively; (2) the gas mass fraction in galaxies is lower; (3) the total baryonic mass function (gas + star) at z=3 does not differ significantly, but shows an increase in the number of relatively massive galaxies at z=1; (4) the baryonic mass fraction of intergalactic medium (IGM) is reduced at z<3 due to more efficient cooling and gas accretion onto galaxies. Our results suggest that the metal cooling enhances the galaxy growth by two different mechanisms: (1) increase of SF efficiency in the local interstellar medium (ISM), and (2) increase of IGM accretion onto galaxies. The former process is effective throughout most of the cosmic history, while the latter is effective only at z<3 when the IGM is sufficiently enriched by metals owing to feedback.

The Korea Invisible Mass Search(KIMS) experiment presents new limits on WIMP-nucleon cross section using the data from an exposure of 3409 kgd taken with low background CsI(Tl) crystals at Yangyang underground laboratory. The most stringent limit on the spin dependent interaction for pure proton case is obtained. The DAMA signal region for both spin independent and spin dependent interactions for the WIMP mass higher than 20 GeV/c^2are excluded by the single experiment with crystal scintillators.

(Abbreviated) We investigate the dynamical mechanisms responsible for producing tidal tails from dwarf satellites using N-body simulations. We identify two important dynamical co-conspirators: 1) the points where the attractive force of the host halo and satellite are balanced do not occur at equal distances from the satellite centre or at the same equipotential value for massive satellites, breaking the morphological symmetry of the leading and trailing tails; and 2) the escaped ejecta in the leading (trailing) tail continues to be decelerated (accelerated) by the satellite's gravity leading to large offsets of the ejecta orbits from the satellite orbit. The effect of the satellite's self gravity decreases only weakly with a decreasing ratio of satellite mass to host halo mass, demonstrating the importance of these effects over a wide range of subhalo masses. Not only will the morphology of the leading and trailing tails for massive satellites be different, but the observed radial velocities of the tails will be displaced from that of the satellite orbit; both the displacement and the peak radial velocity is proportional to satellite mass. If the tails are assumed to follow the progenitor satellite orbits, the tails from satellites with masses greater than 0.0001 of the host halo virial mass in a spherical halo will appear to indicate a flattened halo. Therefore, a constraint on the Milky Way halo shape using tidal streams requires mass-dependent modelling. Similarly, we compute the the distribution of tail orbits both in E_{r}-r^{-2} space and in E-L_{z} space, advocated for identifying satellite stream relics. The acceleration of ejecta by a massive satellite during escape spreads the velocity distribution and obscures the signature of a well-defined ``moving group'' in phase space.

This is the first of two papers reporting observations and analysis of the unusually bright (m_b=14.4), luminous (M_B=-25.5), nearby (z=0.192) narrow-line quasar PHL 1811, focusing on the X-ray properties and the spectral energy distribution. Two Chandra observations reveal a weak X-ray source with a steep spectrum. Variability by a factor of 4 between the two observations separated by 12 days suggest that the X-rays are not scattered emission. The XMM-Newton spectra are modelled in the 0.3--5 keV band by a steep power law with \Gamma = 2.3\pm 0.1, and the upper limit on intrinsic absorption is 8.7 x 10^{20} cm^{-2}. The spectral slopes are consistent with power law indices commonly observed in NLS1s, and it appears that we observe the central engine X-rays directly. Including two recent Swift ToO snapshots, a factor of ~5 variability was observed among the five X-ray observations reported here. In contrast, the UV photometry obtained by the XMM-Newton OM and Swift UVOT, and the HST spectrum reveal no significant UV variability. The \alpha_{ox} inferred from the Chandra and contemporaneous HST spectrum is -2.3 \pm 0.1, significantly steeper than observed from other quasars of the same optical luminosity. The steep, canonical X-ray spectra, lack of absorption, and significant X-ray variability lead us to conclude that PHL 1811 is intrinsically X-ray weak. We also discuss an accretion disk model, and the host galaxy of PHL 1811.

We consider the sensitivity of the circular-orbit adiabatic contraction approximation to the baryon condensation rate and the orbital structure of dark matter halos in the $\Lambda$CDM paradigm. Using one-dimensional hydrodynamic simulations including the dark matter halo mass accretion history and gas cooling, we demonstrate that the adiabatic approximation is approximately valid even though halos and disks may assemble simultaneously. We further demonstrate the validity of the simple approximation for $\Lambda$CDM halos with isotropic velocity distributions using three-dimensional N-body simulations. This result is easily understood: an isotropic velocity distribution in a cuspy halo requires more circular orbits than radial orbits. Conversely, the approximation is poor in the extreme case of a radial orbit halo. It overestimates the response a core dark matter halo, where radial orbit fraction is larger. Because no astronomically relevant models are dominated by low-angular momentum orbits in the vicinity of the disk and the growth time scale is never shorter than a dynamical time, we conclude that the adiabatic contraction approximation is useful in modeling the response of dark matter halos to the growth of a disk.

The Korea Invisible Mass Search (KIMS) collaboration has been carrying out WIMP search experiment with CsI(T$\ell$)crystal detectors at the YanYang Underground Laboratory. A successful reduction of the internal background of the crystal is done and a good pulse shape discrimination is achieved. We report the first result on WIMP search obtained with 237 kg$\cdot$days data using one full-size CsI(T$\ell$)crystal of 6.6 kg mass.