The impact of the streaming between baryons and dark matter on the first structures has been actively explored by recent studies. We investigate how the key results are affected by two popular approximations. One is to implement the streaming by accounting for only the relative motion while assuming "baryons trace dark matter" spatially at the initialization of simulation. This neglects the smoothing on the gas density taking place before the initialization. In our simulation initialized at $z_i=200$, it overestimates the gas density power spectrum up to 40\% at $k\approx10^2~h~\mbox{Mpc}^{-1}$ at $z=20$. Halo mass ($M_h$) and baryonic fraction in halos ($f_{b,h}$) are also overestimated, but the relation between the two remains unchanged. The other approximation tested is to artificially amplify the density/velocity fluctuations in the cosmic mean density to simulate the first minihalos that form in overdense regions. This gives a head start to the halo growth while the subsequent growth is similar to that in the mean density. The growth in a true overdense region, on the other hand, is accelerated gradually in time. For example, raising $\sigma_8$ by 50\% effectively transforms $z\rightarrow\sqrt{1.5}z$ in the halo mass growth history while in 2-$\sigma$ overdensity, the growth is accelerated by a constant in redshift: $z\rightarrow{z+4.8}$. As a result, halos are more grown in the former than in the latter before $z\approx27$ and vice versa after. The $f_{b,h}$-$M_h$ relation is unchanged in those cases as well, suggesting that the Pop III formation rate for a given $M_h$ is insensitive to the tested approximations.

The Virtual Telescope for X-ray Observations (VTXO) will use lightweight Phase Frensel Lenses (PFLs) in a virtual X-ray telescope with 1 km focal length and with nearly 50 milli-arcsecond angular resolution. Laboratory characterization of PFLs have demonstrated near diffraction-limited angular resolution in the X-ray band, but they require long focal lengths to achieve this quality of imaging. VTXO is formed by using precision formation flying of two SmallSats: a smaller, 6U OpticsSat that houses the PFLs and navigation beacons while a larger, ESPA-class DetectorSat contains an X-ray camera, a charged-particle radiation monitor, a precision star tracker, and the propulsion for the formation flying. The baseline flight dynamics uses a highly-elliptical supersynchronous geostationary transfer orbit to allow the inertial formation to form and hold around the 90,000 km apogee for 10 hours of the 32.5-hour orbit with nearly a year mission lifetime. The guidance, navigation, and control (GN&C) for the formation flying uses standard CubeSat avionics packages, a precision star tracker, imaging beacons on the OpticsSat, and a radio ranging system that also serves as an inter-satellite communication link. VTXO's fine angular resolution enables measuring the environments nearly an order of magnitude closer to the central engines of bright compact X-ray sources compared to the current state of the art. This X-ray imaging capability allows for the study of the effects of dust scattering nearer to the central objects such as Cyg X-3 and GX 5-1, for the search for jet structure nearer to the compact object in X-ray novae such as Cyg X-1 and GRS 1915+105, and for the search for structure in the termination shock of in the Crab pulsar wind nebula. The VTXO development was supported as one of the selected 2018 NASA Astrophysics SmallSat Study (AS3) missions.

The tomographic AP method is so far the best method in separating the Alcock-Paczynski (AP) signal from the redshift space distortion (RSD) effects and deriving powerful constraints on cosmological parameters using the $\lesssim40h^{-1}\ \rm Mpc$ clustering region. To guarantee that the method can be easily applied to the future large scale structure (LSS) surveys, we study the possibility of estimating the systematics of the method using fast simulation method. The major contribution of the systematics comes from the non-zero redshift evolution of the RSD effects, which is quantified by $\hat\xi_{\Delta s}(\mu,z)$ in our analysis, and estimated using the BigMultidark exact N-body simulation and approximate COLA simulation samples. We find about 5\%/10\% evolution when comparing the $\hat\xi_{\Delta s}(\mu,z)$ measured as $z=0.5$/$z=1$ to the measurements at $z=0$. We checked the inaccuracy in the 2pCFs computed using COLA, and find it 5-10 times smaller than the intrinsic systematics of the tomographic AP method, indicating that using COLA to estimate the systematics is good enough. Finally, we test the effect of halo bias, and find $\lesssim$1.5\% change in $\hat\xi_{\Delta s}$ when varying the halo mass within the range of $2\times 10^{12}$ to $10^{14}$ $M_{\odot}$. We will perform more studies to achieve an accurate and efficient estimation of the systematics in redshift range of $z=0-1.5$.

The redshift-space distortion (RSD) in the observed distribution of galaxies is known as a powerful probe of cosmology. Observations of large-scale RSD have given tight constraints on the linear growth rate of the large-scale structures in the universe. On the other hand, the small-scale RSD, caused by galaxy random motions inside clusters, has not been much used in cosmology, but also has cosmological information because universes with different cosmological parameters have different halo mass functions and virialized velocities. We focus on the projected correlation function $w(r_p)$ and the multipole moments $\xi_l$ on small scales ($1.4$ to $30\ h^{-1}\rm{Mpc}$). Using simulated galaxy samples generated from a physically motivated most bound particle (MBP)-galaxy correspondence scheme in the Multiverse Simulation, we examine the dependence of the small-scale RSD on the cosmological matter density parameter $\Omega_m$, the satellite velocity bias with respect to MBPs, $b_v^s$, and the merger-time-scale parameter $\alpha$. We find that $\alpha=1.5$ gives an excellent fit to the $w(r_p)$ and $\xi_l$ measured from the SDSS-KIAS value added galaxy catalog. We also define the ``strength'' of Fingers-of-God as the ratio of the parallel and perpendicular size of the contour in the two-point correlation function set by a specific threshold value and show that the strength parameter helps constraining $(\Omega_m, b_v^s, \alpha)$ by breaking the degeneracy among them. The resulting parameter values from all measurements are $(\Omega_m,b_v^s)=(0.272\pm0.013,0.982\pm0.040)$, indicating a slight reduction of satellite galaxy velocity relative to the MBP. However, considering that the average MBP speed inside haloes is $0.94$ times the dark matter velocity dispersion, the main drivers behind the galaxy velocity bias are gravitational interactions, rather than baryonic effects.

COSINE-200 is the next phase of the ongoing COSINE-100 experiment. The main purpose of the experiment is the performance of an unambiguous verification of the annual modulation signals observed by the DAMA experiment. The success of the experiment critically depends on the production of a 200 kg array of ultra-pure NaI(Tl) crystal detectors that have lower backgrounds than the DAMA crystals. The purification of raw powder is the initial but important step toward the production of ultra-pure NaI(Tl) detectors. We have already demonstrated that fractional recrystallization from water solutions is an effective method for the removal of the problematic K and Pb elements. For the mass production of purified powder, a clean facility for the fractional recrystallization had been constructed at the Institute for Basic Science (IBS), Korea. Here, we report the design of the purification process, material recovery, and performance of the NaI powder purification facility.

The Center for Underground Physics (CUP) of the Institute for Basic Science (IBS) is searching for the neutrinoless double-beta decay (0{\nu}\b{eta}\b{eta}) of 100Mo in the molybdate crystals of the AMoRE experiment. The experiment requires pure scintillation crystals to minimize internal backgrounds that can affect the 0{\nu}\b{eta}\b{eta} signal. For the last few years, we have been growing and studying Li2MoO4 crystals in a clean-environment facility to minimize external contamination during the crystal growth. Before growing Li2100MoO4 crystal, we have studied Li2natMoO4 crystal growth by a conventional Czochralski (CZ) grower. We grew a few different kinds of Li2natMO4 crystals using different raw materials in a campaign to minimize impurities. We prepared the fused Al2O3 refractories for the growth of ingots. Purities of the grown crystals were measured with high purity germanium detectors and by inductively coupled plasma mass spectrometry. The results show that the Li2MoO4 crystal has purity levels suitable for rare-event experiments. In this study, we present the growth of Li2MoO4 crystals at CUP and their purities.

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.

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.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

We present newly discovered dwarf galaxy candidates in deep and wide-field images of NGC 1291 obtained with the Korea Microlensing Telescope Network. We identify 15 dwarf galaxy candidates by visual inspection. Using imaging simulations, we demonstrate that the completeness rate of our detection is greater than 70% for the central surface brightness value of ${\mu}_{0,R}\lesssim26$ mag arcsec$^{-2}$ and for magnitudes $M_R\lesssim-10$ mag. The structural and photometric properties of the dwarf galaxy candidates appear to be broadly consistent with those of ordinary dwarf galaxies in nearby groups and clusters, with ${\mu}_{0,R}\sim$ 22.5 to 26.5 mag arcsec$^{-2}$ and effective radii of 200 pc to 1 kpc. The dwarf galaxy candidates show a concentration towards NGC 1291 and tend to be redder the closer they are to the center, possibly indicating that they are associated with NGC 1291. The dwarf candidates presented in this paper appear to be bluer than those in denser environments, revealing that the quenching of star formation in dwarf galaxies is susceptible to the environment, while the morphology shaping is not.

The surface $^{210}$Pb is one of the main background sources for dark-matter-search experiments using NaI(Tl) crystals, its spectral features associated with the beta-decay events for energies less than 60 keV depends on the depth distribution of $^{210}$Pb in the surface of an NaI(Tl) crystal. Therefore, we must understand the profile of surface $^{210}$Pb to precisely model the background measurement in the low-energy region for the low-background experiment using NaI(Tl) crystals. We estimate the depth profile of the surface $^{210}$Pb contamination by modeling the measured spectrum of the alpha emission from the decay of $^{210}$Po at the decay sequence of the surface $^{210}$Pb contamination that is obtained using an $^{222}$Rn-contaminated crystal. Moreover, we describe the energy spectra of the surface contaminations by performing Geant4 Monte Carlo simulations that are used for fitting an exponential curve, which is a function of the surface depth, with the measured data, by using the log likelihood method. The low- and high-energy events from the beta decay of surface $^{210}$Pb are also modeled to improve the depth profile for shallow depths. We simulate the energy spectra from beta decays of $^{210}$Pb that are exponentially distributed in the surface by following two exponential functions that consider mean-depth parameters in the data fitting; we observed that the energy spectra are in good agreement with the measured data.

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.

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.

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.

We present our discovery and analysis of dwarf galaxies in the NGC 3585 galaxy group by the KMTNet Supernova Program. Using deep stack images reaching $\simeq$ 28 mag arcsec$^{-2}$ in $BVI$, we discovered 46 dwarf galaxy candidates distributed in a 7 square degree field. The dwarf galaxy candidates exhibit central surface brightness as faint as $\mu_{0,V} = 26.2$ mag arcsec$^{-2}$, with effective radii larger than 150 pc and total absolute magnitudes brighter than $M_V \approx -10$ mag, if at the distance of NGC 3585. The dwarf galaxy surface number density decreases with projected distance from NGC 3585. We estimate the background contamination to be about 20% based both on the number density profile and on diffuse galaxy counts in a control field. The dwarf galaxy colors and Sersic structural parameters are consistent with those found for other dwarf galaxies. Unusually, there is no indication of a change of color or brightness in the dwarf galaxy candidates with projected distance from the group center. Approximately 20% of them contain an unresolved nucleus. The nucleated fraction is larger for brighter (and redder) galaxies but is independent of distance from the group center. We identify four ultra-diffuse galaxy candidates, all near the group center. We interpret these spatial properties as suggesting that the NGC 3585 group might be dynamically younger than the typical group. The galaxy luminosity function of the NGC 3585 group has a faint-end slope of $\alpha\approx -1.39$, which is roughly consistent with the slopes found for other nearby groups. The possible dependence of the slope on global group properties is still unclear, and continues to motivate our homogeneous deep survey of dozens of nearby groups.

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.

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.

We present new constraints on the dark matter-induced annual modulation signal using 1.7 years, of COSINE-100 data with a total exposure of 97.7 kg$\cdot$years. The COSINE-100 experiment, consisting of 106 kg of NaI(Tl) target material, is designed to carry out a model-independent test of DAMA/LIBRA's claim of WIMP discovery by searching for the same annual modulation signal using the same NaI(Tl) target. The crystal data show a 2.7 cpd/kg/keV background rate on average in the 2--6 keV energy region of interest. Using a $\chi$-squared minimization method we observe best fit values for modulation amplitude and phase of 0.0092$\pm$0.0067 cpd/kg/keV and 127.2$\pm$45 d, respectively.

We present new constraints on the dark matter-induced annual modulation signal using 1.7 years, of COSINE-100 data with a total exposure of 97.7 kg$\cdot$years. The COSINE-100 experiment, consisting of 106 kg of NaI(Tl) target material, is designed to carry out a model-independent test of DAMA/LIBRA's claim of WIMP discovery by searching for the same annual modulation signal using the same NaI(Tl) target. The crystal data show a 2.7 cpd/kg/keV background rate on average in the 2--6 keV energy region of interest. Using a $\chi$-squared minimization method we observe best fit values for modulation amplitude and phase of 0.0092$\pm$0.0067 cpd/kg/keV and 127.2$\pm$45 d, respectively.

We present results from a search for solar axions with the COSINE-100 experiment. We find no evidence of solar axion events from a data-set of 6,303.9 kg$\cdot$days exposure and set a 90\,\% confidence level upper limit on the axion-electron coupling, $g_{ae}$, at 1.70~$\times$~$10^{-11}$ for an axion mass less than 1\,keV/c$^2$. This limit excludes QCD axions heavier than 0.59\,eV/c$^2$ in the DFSZ model and 168.1\,eV/c$^2$ in the KSVZ model.

We present results from a search for solar axions with the COSINE-100 experiment. We find no evidence of solar axion events from a data-set of 6,303.9 kg$\cdot$days exposure and set a 90\,\% confidence level upper limit on the axion-electron coupling, $g_{ae}$, at 1.70~$\times$~$10^{-11}$ for an axion mass less than 1\,keV/c$^2$. This limit excludes QCD axions heavier than 0.59\,eV/c$^2$ in the DFSZ model and 168.1\,eV/c$^2$ in the KSVZ model.

We develop an improved Alcock-Paczynski (AP) test method that uses the redshift-space two-point correlation function (2pCF) of galaxies. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF, which should not change apart from the expected small non-linear evolution. An incorrect choice of cosmology used to convert redshift to comoving distance will manifest itself as redshift-dependent 2pCF. Our method decomposes the redshift difference of the two-dimensional correlation function into the Legendre polynomials whose amplitudes are modeled by radial fitting functions. Our likelihood analysis with this 2-D fitting scheme tightens the constraints on $\Omega_m$ and ${w}$ by $\sim 40\%$ compared to the method of Li et al. (2016, 2017, 2018) that uses one dimensional angular dependence only. We also find that the correction for the non-linear evolution in the 2pCF has a non-negligible cosmology dependence, which has been neglected in previous similar studies by Li et al.. With an accurate accounting for the non-linear systematics and use of full two-dimensional shape information of the 2pCF down to scales as small as $5~h^{-1}{\rm Mpc}$ it is expected that the AP test with redshift-space galaxy clustering anisotropy can be a powerful method to constrain the expansion history of the universe.

Assuming a standard Maxwellian for the WIMP velocity distribution, we obtain the bounds from null WIMP search results of 59.5 days of COSINE-100 data on the DAMA/LIBRA-phase2 modulation effect within the context of the non-relativistic effective theory of WIMP-nucleus scattering. Here, we systematically assume that one of the effective operators allowed by Galilean invariance dominates in the effective Hamiltonian of a spin-1/2 dark matter (DM) particle. We find that, although DAMA/LIBRA and COSINE-100 use the same sodium-iodide target, the comparison of the two results still depends on the particle-physics model. This is mainly due to two reasons: i) the WIMP signal spectral shape; ii) the expected modulation fractions, when the upper bound on the time-averaged rate in COSINE-100 is converted into a constraint on the annual modulation component in DAMA/LIBRA. We find that the latter effect is the dominant one. For several effective operators the expected modulation fractions are larger than in the standard spin-independent or spin-dependent interaction cases. As a consequence, compatibility between the modulation effect observed in DAMA/LIBRA and the null result from COSINE-100 is still possible for several non-relativistic operators. At low WIMP masses such relatively high values of the modulation fractions arise because COSINE-100 is mainly sensitive to WIMP-sodium scattering events, due to the higher threshold compared to DAMA/LIBRA. A next COSINE analysis is expected to have a full sensitivity for the 5$\sigma$ region of DAMA/LIBRA.

Assuming a standard Maxwellian for the WIMP velocity distribution, we obtain the bounds from null WIMP search results of 59.5 days of COSINE-100 data on the DAMA/LIBRA-phase2 modulation effect within the context of the non-relativistic effective theory of WIMP-nucleus scattering. Here, we systematically assume that one of the effective operators allowed by Galilean invariance dominates in the effective Hamiltonian of a spin-1/2 dark matter (DM) particle. We find that, although DAMA/LIBRA and COSINE-100 use the same sodium-iodide target, the comparison of the two results still depends on the particle-physics model. This is mainly due to two reasons: i) the WIMP signal spectral shape; ii) the expected modulation fractions, when the upper bound on the time-averaged rate in COSINE-100 is converted into a constraint on the annual modulation component in DAMA/LIBRA. We find that the latter effect is the dominant one. For several effective operators the expected modulation fractions are larger than in the standard spin-independent or spin-dependent interaction cases. As a consequence, compatibility between the modulation effect observed in DAMA/LIBRA and the null result from COSINE-100 is still possible for several non-relativistic operators. At low WIMP masses such relatively high values of the modulation fractions arise because COSINE-100 is mainly sensitive to WIMP-sodium scattering events, due to the higher threshold compared to DAMA/LIBRA. A next COSINE analysis is expected to have a full sensitivity for the 5$\sigma$ region of DAMA/LIBRA.

We present the discovery and the photometric and spectroscopic study of H-rich Type II supernova (SN) KSP-SN-2016kf (SN2017it) observed in the KMTNet Supernova Program in the outskirts of a small irregular galaxy at $z\simeq0.043$ within a day from the explosion. Our high-cadence, multi-color ($BVI$) light curves of the SN show that it has a very long rise time ($t_\text{rise}\simeq 20$ days in $V$ band), a moderately luminous peak ($M_V\simeq -$17.6 mag), a notably luminous and flat plateau ($M_V\simeq -$17.4 mag and decay slope $s\simeq0.53$ mag per 100 days), and an exceptionally bright radioactive tail. Using the color-dependent bolometric correction to the light curves, we estimate the $^{56}$Ni mass powering the observed radioactive tail to be $0.10\pm0.01$ M$_\odot$, making it a H-rich Type II SN with one of the largest $^{56}$Ni masses observed to date. The results of our hydrodynamic simulations of the light curves constrain the mass and radius of the progenitor at the explosion to be $\sim$15 M$_\odot$ (evolved from a star with an initial mass of $\sim$ 18.8 M$_\odot$) and $\sim1040$ R$_\odot$, respectively, with the SN explosion energy of $\sim 1.3\times 10^{51}$ erg s$^{-1}$. The above-average mass of the KSP-SN-2016kf progenitor, together with its low metallicity $ Z/Z_\odot \simeq0.1-0.4$ obtained from spectroscopic analysis, is indicative of a link between the explosion of high-mass red supergiants and their low-metallicity environment. The early part of the observed light curves shows the presence of excess emission above what is predicted in model calculations, suggesting there is interaction between the ejecta and circumstellar material. We further discuss the implications of the high progenitor initial mass and low-metallicity environment of KSP-SN-2016kf on our understanding of the origin of Type II SNe.

We present a multi-color, high-cadence photometric study of a distant dwarf nova KSP-OT-201611a discovered by the Korea Microlensing Telescope Network Supernova Program. From October 2016 to May 2017, two outbursts, which comprises a super/long outburst followed by a normal/short outburst separated by $\sim$91 days, were detected in the $BVI$ bands. The shapes and amplitudes of the outbursts reveal the nature of KSP-OT-201611a to be an SU UMa- or U Gem-type dwarf nova. Color variations of periodic humps in the super/long outburst possibly indicate that KSP-OT-201611a is an SU UMa-type dwarf nova. The super and normal outbursts show distinctively different color evolutions during the outbursts due most likely to the difference of time when the cooling wave is formed in the accretion disk. The outburst peak magnitudes and the orbital period of the dwarf nova indicate that it is at a large Galactocentric distance ($\sim$13.8 kpc) and height ($\sim$1.7 kpc) from the Galactic plane. KSP-OT-201611a, therefore, may provide a rare opportunity to study the accretion disk process of Population II dwarf novae.

Observations of galaxies and primordial radiation suggest that the Universe is made mostly of non-luminous dark matter. Several types of new fundamental particles have been proposed as candidates for dark matter such as weakly interacting massive particles (WIMPs) but no definitive signal has been seen despite concerted efforts by many collaborations. One exception is the much-debated claim by the DAMA collaboration of a statistically significant annual modulation in the event rate of their experiment with a period and phase consistent with that expected from WIMP dark matter. Several groups have been working to develop experiments with the aim of reproducing DAMA's results using the same target medium. Here we report results from the initial operation of the COSINE-100 experiment. COSINE-100 uses sodium iodide as the target medium-the same medium as DAMA-and is designed to carry out a model-independent test of DAMA's claim. Initial data based on the first 59.5 days indicate that there is no excess of events over the expected background, confirming that DAMA's annual modulation signal is in severe tension with results from other experiments under the assumption of dark matter having spin independent interactions and the Standard Halo Model. COSINE-100 is now taking data to study the presence of dark matter-induced annual modulation in the event rate of the sodium iodide detectors.

Observations of galaxies and primordial radiation suggest that the Universe is made mostly of non-luminous dark matter. Several types of new fundamental particles have been proposed as candidates for dark matter such as weakly interacting massive particles (WIMPs) but no definitive signal has been seen despite concerted efforts by many collaborations. One exception is the much-debated claim by the DAMA collaboration of a statistically significant annual modulation in the event rate of their experiment with a period and phase consistent with that expected from WIMP dark matter. Several groups have been working to develop experiments with the aim of reproducing DAMA's results using the same target medium. Here we report results from the initial operation of the COSINE-100 experiment. COSINE-100 uses sodium iodide as the target medium-the same medium as DAMA-and is designed to carry out a model-independent test of DAMA's claim. Initial data based on the first 59.5 days indicate that there is no excess of events over the expected background, confirming that DAMA's annual modulation signal is in severe tension with results from other experiments under the assumption of dark matter having spin independent interactions and the Standard Halo Model. COSINE-100 is now taking data to study the presence of dark matter-induced annual modulation in the event rate of the sodium iodide detectors.

The tomographic Alcock-Paczynski (AP) method can result in tight cosmological constraints by using small and intermediate clustering scales of the large scale structure (LSS) of the galaxy distribution. By focusing on the redshift dependence, the AP distortion can be distinguished from the distortions produced by the redshift space distortions (RSD). In this work, we combine the tomographic AP method with other recent observational datasets of SNIa+BAO+CMB+$H_0$ to reconstruct the dark energy equation-of-state $w$ in a non-parametric form. The result favors a dynamical DE at $z\lesssim1$, and shows a mild deviation ($\lesssim2\sigma$) from $w=-1$ at $z=0.5-0.7$. We find the addition of the AP method improves the low redshift ($z\lesssim0.7$) constraint by $\sim50\%$.

We survey globular clusters (GCs) in M85 using $ugi$-band images of a $1^{\circ} \times 1^{\circ}$ field obtained with the MegaCam at the 3.6 m Canada-France-Hawaii Telescope. We identify 1318 GC candidates with 20.0 mag $< g_0 <$ 23.5 mag in the entire survey region. Their radial number density profile is well fit by a S{\'e}rsic profile with $n$ = 2.58$^{+0.43}_{-0.33}$ and effective radius $R_{\rm e,GCS}$ = 4$\rlap{.}{'}$14 (= 22 kpc), showing that the candidates at $R < 20'$ are mostly genuine GCs in M85. We estimate the total number of GCs, $N$(total) = $1216^{+82}_{-50}$, and the specific frequency, $S_N = 1.41^{+0.10}_{-0.06}$. The overall color distribution of the GCs in M85 is bimodal, but the GCs in the central region at $R < 2'$ do not show a bimodal distribution clearly. The radial number density profile and surface number density map of the blue GCs (BGCs) show more extended structures than those of the red GCs (RGCs). The spatial distributions of both BGCs and RGCs are elongated, similar to that of the galaxy stellar light. The number fraction of the RGCs in the central region is much smaller compared to those in other early-type galaxies of similar luminosity. The mean $(g-i)_0$ color of the RGCs in M85 is about 0.1 mag bluer than typical values for other Virgo early-type galaxies of similar luminosity, indicating that a significant fraction of the RGCs in M85 may be younger than typical GCs. These results indicate that M85 might have undergone a major wet merger recently.

A search for inelastic boosted dark matter (iBDM) using the COSINE-100 detector with 59.5 days of data is presented. This relativistic dark matter is theorized to interact with the target material through inelastic scattering with electrons, creating a heavier state that subsequently produces standard model particles, such as an electron-positron pair. In this study, we search for this electron-positron pair in coincidence with the initially scattered electron as a signature for an iBDM interaction. No excess over the predicted background event rate is observed. Therefore, we present limits on iBDM interactions under various hypotheses, one of which allows us to explore an area of the experimental search for iBDM using a terrestrial detector.

A search for inelastic boosted dark matter (iBDM) using the COSINE-100 detector with 59.5 days of data is presented. This relativistic dark matter is theorized to interact with the target material through inelastic scattering with electrons, creating a heavier state that subsequently produces standard model particles, such as an electron-positron pair. In this study, we search for this electron-positron pair in coincidence with the initially scattered electron as a signature for an iBDM interaction. No excess over the predicted background event rate is observed. Therefore, we present limits on iBDM interactions under various hypotheses, one of which allows us to explore an area of the experimental search for iBDM using a terrestrial detector.

With the advent of high power lasers, new opportunities have opened up for simulating astrophysical processes in the laboratory. We show that 2nd-order Fermi acceleration can be directly investigated at the National Ignition Facility, Livermore. This requires measuring the momentum-space diffusion of 3 MeV protons produced within a turbulent plasma generated by a laser. Treating Fermi acceleration as a biased diffusion process, we show analytically that a measurable broadening of the initial proton distribution is then expected for particles exiting the plasma.

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.

We present a sample of 54 disk galaxies which have well developed extraplanar structures. We selected them using visual inspections from the color images of the Sloan Digital Sky Survey. Since the sizes of the extraplanar structures are comparable to the disks, they are considered as prominent stellar halos rather than large bulges. A single S\'ersic profile fitted to the surface brightness along the minor-axis of the disk shows a luminosity excess in the central regions for the majority of sample galaxies. This central excess is considered to be caused by the central bulge component. The mean S\'ersic index of the single component model is $1.1\pm0.9$. A double S\'ersic profile model that employs $n=1$ for the inner region, and varying $n$ for the outer region, provides a better fit than the single S\'ersic profile model. For a small fraction of galaxies, a S\'ersic profile fitted with $n=4$ for the inner region gives similar results. There is a weak tendency of increasing $n$ with increasing luminosity and central velocity dispersion, but there is no dependence on the local background density.

COSINE-100 is a dark matter direct detection experiment designed to test the annual modulation signal observed by the DAMA/LIBRA experiment. COSINE-100 consists of 8 NaI(Tl) crystals with a total mass of 106 kg, a 2200 L liquid scintillator veto, and 37 muon detector panels. We present details of the data acquisition system of COSINE-100, including waveform storage using flash analog-to-digital converters for crystal events and integrated charge storage using charge-sensitive analog-to-digital converters for liquid scintillator and plastic scintillator muon veto events. We also discuss several trigger conditions developed in order to distinguish signal events from photomultiplier noise events. The total trigger rate observed for the crystal/liquid scintillator (plastic scintillator) detector is 15 Hz (24 Hz).

COSINE-100 is a dark matter direct detection experiment designed to test the annual modulation signal observed by the DAMA/LIBRA experiment. COSINE-100 consists of 8 NaI(Tl) crystals with a total mass of 106 kg, a 2200 L liquid scintillator veto, and 37 muon detector panels. We present details of the data acquisition system of COSINE-100, including waveform storage using flash analog-to-digital converters for crystal events and integrated charge storage using charge-sensitive analog-to-digital converters for liquid scintillator and plastic scintillator muon veto events. We also discuss several trigger conditions developed in order to distinguish signal events from photomultiplier noise events. The total trigger rate observed for the crystal/liquid scintillator (plastic scintillator) detector is 15 Hz (24 Hz).

The COSINE-100 dark matter search experiment is an array of NaI(Tl) crystal detectors located in the Yangyang Underground Laboratory (Y2L). To understand measured backgrounds in the NaI(Tl) crystals we have performed Monte Carlo simulations using the Geant4 toolkit and developed background models for each crystal that consider contributions from both internal and external sources, including cosmogenic nuclides. The background models are based on comparisons of measurement data with Monte Carlo simulations that are guided by a campaign of material assays and are used to evaluate backgrounds and identify their sources. The average background level for the six crystals (70 kg total mass) that are studied is 3.5 counts/day/keV/kg in the (2-6) keV energy interval. The dominant contributors in this energy region are found to be $^{210}$Pb and $^3$H.

The COSINE-100 dark matter search experiment is an array of NaI(Tl) crystal detectors located in the Yangyang Underground Laboratory (Y2L). To understand measured backgrounds in the NaI(Tl) crystals we have performed Monte Carlo simulations using the Geant4 toolkit and developed background models for each crystal that consider contributions from both internal and external sources, including cosmogenic nuclides. The background models are based on comparisons of measurement data with Monte Carlo simulations that are guided by a campaign of material assays and are used to evaluate backgrounds and identify their sources. The average background level for the six crystals (70 kg total mass) that are studied is 3.5 counts/day/keV/kg in the (2-6) keV energy interval. The dominant contributors in this energy region are found to be $^{210}$Pb and $^3$H.

It has been known that decays of daughter elements of $^{222}$Rn on the surface of a detector cause significant background at energies below 10 keV. In particular $^{210}$Pb and $^{210}$Po decays on the crystal surface result in significant background for dark matter search experiments with NaI(Tl) crystals. In this report, measurement of $^{210}$Pb and $^{210}$Po decays on surfaces are obtained by using a $^{222}$Rn contaminated crystal. Alpha decay events of $^{210}$Po on the surface are measured by coincidence requirements of two attached crystals. Due to recoiling of $^{206}$Pb, rapid nuclear recoil events are observed. A mean time characterization demonstrates that $^{206}$Pb recoil events can be statistically separated from those of sodium or iodine nuclear recoil events, as well as electron recoil events.

M85 is a peculiar S0 galaxy in Virgo and is a well-known merger remnant. In this paper, we present the first spectroscopic study of globular clusters (GCs) in M85. We obtain spectra for 21 GC candidates and the nucleus of M85 using the Gemini Multi-Object Spectrograph on the Gemini North 8.1 m telescope. From their radial velocities, 20 of the GCs are found to be members of M85. We find a strong rotation signal of the M85 GC system with a rotation amplitude of 235 km s$^{-1}$. The rotation axis of the GC system has a position angle of about 161$^{\circ}$, which is 51$\rlap{.}{^\circ}$5 larger than that of the stellar light. The rotation-corrected radial velocity dispersion of the GC system is estimated to be $\sigma_{\rm r,cor} = $ 160 km s$^{-1}$. The rotation parameter $\Omega R_{\rm icor}/\sigma_{\rm r,cor}$ of the GC system is derived to be 1.47$^{+1.05}_{-0.48}$, which is one of the largest among known early-type galaxies. The ages and metallicities of the GCs, which show the same trend as the results based on Lick indices, are derived from full spectrum fitting (ULySS). About a half of the GCs are an intermediate-age population of which the mean age is $\sim$ 3.7 $\pm$ 1.9 Gyr, having a mean [Fe/H] value of --0.26. The other half are old and metal-poor. These results suggest that M85 experienced a wet merging event about 4 Gyr ago, forming a significant population of star clusters. The strong rotational feature of the GC system can be explained by an off-center major merging.

We investigate surface brightness profiles (SBPs) of dwarf galaxies in field, group, and cluster environments. With deep $BVI$ images from the Korea Microlensing Telescope Network Supernova Program, SBPs of 38 dwarfs in the NGC 2784 group are fitted by a single exponential or double exponential model. We find that 53% of the dwarfs are fitted with single exponential profiles ("Type I"), while 47% of the dwarfs show double exponential profiles. 37% of all dwarfs have smaller sizes of the outer part than inner part ("Type II"), while 10% have a larger sized outer than inner part ("Type III"). We compare these results with those in the field and in the Virgo cluster, where SBP types of 102 field dwarfs are compiled from a previous study and SBP types of 375 cluster dwarfs are measured using SDSS $r$-band images. As a result, the distributions of SBP types are different in three environments. Common SBP types of the field, the NGC 2784 group, and the Virgo cluster are Type II, Type I and II, and Type I and III profiles, respectively. After comparing sizes of dwarfs in different environments, we suggest that since sizes of some dwarfs are changed due to the environmental effects, SBP types are able to be transformed and the distributions of SBP types in three environments are different. We discuss possible environmental mechanisms on the transformation of SBP types.

We perform an anisotropic clustering analysis of 1,133,326 galaxies from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR) 12 covering the redshift range $0.15<z<0.69$. The geometrical distortions of the galaxy positions, caused by incorrect cosmological model assumptions, are captured in the anisotropic two-point correlation function on scales 6 -- 40 $h^{-1}\rm Mpc$. The redshift evolution of this anisotropic clustering is used to place constraints on the cosmological parameters. We improve the methodology of Li et al. 2016, to enable efficient exploration of high dimensional cosmological parameter spaces, and apply it to the Chevallier-Polarski-Linder parametrization of dark energy, $w=w_0+w_a{z}/({1+z})$. In combination with the CMB, BAO, SNIa and $H_0$ from Cepheid data, we obtain $\Omega_m = 0.301 \pm 0.008,\ w_0 = -1.042 \pm 0.067,\ $ and $w_a = -0.07 \pm 0.29$ (68.3\% CL). Adding our new AP measurements to the aforementioned results reduces the error bars by $\sim$30 -- 40\% and improves the dark energy figure of merit by a factor of $\sim$2. We check the robustness of the results using realistic mock galaxy catalogues.

COSINE is a dark matter search experiment based on an array of low background NaI(Tl) crystals located at the Yangyang underground laboratory. The assembly of COSINE-100 was completed in the summer of 2016 and the detector is currently collecting physics quality data aimed at reproducing the DAMA/LIBRA experiment that reported an annual modulation signal. Stable operation has been achieved and will continue for at least two years. Here, we describe the design of COSINE-100, including the shielding arrangement, the configuration of the NaI(Tl) crystal detection elements, the veto systems, and the associated operational systems, and we show the current performance of the experiment.

COSINE is a dark matter search experiment based on an array of low background NaI(Tl) crystals located at the Yangyang underground laboratory. The assembly of COSINE-100 was completed in the summer of 2016 and the detector is currently collecting physics quality data aimed at reproducing the DAMA/LIBRA experiment that reported an annual modulation signal. Stable operation has been achieved and will continue for at least two years. Here, we describe the design of COSINE-100, including the shielding arrangement, the configuration of the NaI(Tl) crystal detection elements, the veto systems, and the associated operational systems, and we show the current performance of the experiment.

The hierarchical galaxy formation scenario in the Cold Dark Matter cosmogony with a non-vanishing cosmological constant and geometrically flat space has been very successful in explaining the large-scale distribution of galaxies. However, there have been claims that the scenario predicts too many satellite galaxies associated with massive galaxies compared to observations, called the missing satellite galaxy problem. Isolated groups of galaxies hosted by passively evolving massive early-type galaxies are ideal laboratories for finding the missing physics in the current theory. Here we report from a deep spectroscopic survey of such satellite systems that isolated massive early-type galaxies with no recent star formation through wet mergers or accretion have almost no satellite galaxies fainter than the r-band absolute magnitude of about Mr =-14. If only early-type satellites are used, the cutoff is at somewhat brighter magnitude of about Mr =-15. Such a cutoff has not been found in other nearby satellite galaxy systems hosted by late-type galaxies or those with merger features. Various physical properties of satellites depend strongly on the host-centric distance. Our observation indicates that the satellite galaxy luminosity function is largely determined by the interaction of satellites with the environment provided by their host, which sheds light on the missing satellite galaxy problem.

The COSINE-100 dark matter search experiment has started taking physics data with the goal of performing an independent measurement of the annual modulation signal observed by DAMA/LIBRA. A muon detector was constructed by using plastic scintillator panels in the outermost layer of the shield surrounding the COSINE-100 detector. It is used to detect cosmic ray muons in order to understand the impact of the muon annual modulation on dark matter analysis. Assembly and initial performance test of each module have been performed at a ground laboratory. The installation of the detector in Yangyang Underground Laboratory (Y2L) was completed in the summer of 2016. Using three months of data, the muon underground flux was measured to be 328 $\pm$ 1(stat.)$\pm$ 10(syst.) muons/m$^2$/day. In this report, the assembly of the muon detector and the results from the analysis are presented.

The COSINE-100 dark matter search experiment has started taking physics data with the goal of performing an independent measurement of the annual modulation signal observed by DAMA/LIBRA. A muon detector was constructed by using plastic scintillator panels in the outermost layer of the shield surrounding the COSINE-100 detector. It is used to detect cosmic ray muons in order to understand the impact of the muon annual modulation on dark matter analysis. Assembly and initial performance test of each module have been performed at a ground laboratory. The installation of the detector in Yangyang Underground Laboratory (Y2L) was completed in the summer of 2016. Using three months of data, the muon underground flux was measured to be 328 $\pm$ 1(stat.)$\pm$ 10(syst.) muons/m$^2$/day. In this report, the assembly of the muon detector and the results from the analysis are presented.

Poorly understood "baryonic physics" impacts our ability to predict the power spectrum of the kinetic Sunyaev-Zel'dovich (kSZ) effect. We study this in one sample high resolution simulation of galaxy formation and feedback, Illustris. The high resolution of Illustris allows us to probe the kSZ power spectrum on multipoles $\ell =10^3-3\times 10^4$. Strong AGN feedback in Illustris nearly wipes out gas fluctuations at $k\gtrsim1~h~\rm{Mpc}^{-1}$ and at late times, likely somewhat under predicting the kSZ power generated at $z\lesssim 1$. The post-reionization kSZ power spectrum for Illustris is well-fit by $\mathcal{D}^{z<6}_{\ell} = 1.38[\ell/3000]^{0.21}~\mu K^2$ over $3000\lesssim\ell\lesssim10000$, somewhat lower than most other reported values but consistent with the analysis of Shaw et al. Our analysis of the bias of free electrons reveals subtle effects associated with the multi-phase gas physics and stellar fractions that affect even linear scales. In particular there are fewer electrons in biased galaxies, due to gas cooling and star formation, and this leads to an electron bias less than one even at low wavenumbers. The combination of bias and electron fraction that determines the overall suppression is relatively constant, $f_e^2b^2_{e0} \sim 0.7$, but more simulations are needed to see if this is Illustris-specific. By separating the kSZ power into different terms, we find at least $6\, (10)\%$ of the signal at $\ell=3000\, (10000)$ comes from non-Gaussian connected four-point density and velocity correlations, $\left<\delta v \delta v\right>_{c}$, even without correcting for the Illustris simulation box size. A challenge going forward will be to accurately model long-wave velocity modes simultaneously with Illustris-like high resolution to capture the complexities of galaxy formation and its correlations with large scale flows.

We present $BVI$ surface photometry of 31 dwarf galaxy candidates discovered in a deep image stack from the KMTNet Supernova Program of $\sim$ 30 square degrees centered on the nearby NGC 2784 galaxy group. Our final images have a 3$\sigma$ surface brightness detection limit of $\mu_V\approx 28.5$ mag arcsec$^{-2}$. The faintest central surface brightness that we measure is $\mu_{0,V} = 26.1$ mag arcsec$^{-2}$. If these candidates are at the distance of NGC 2784, then they have absolute magnitudes greater than $M_V = -9.5$ mag and effective radii larger than 170 pc. Their radial number density decreases exponentially with distance from the center of NGC 2784 until it flattens beyond a radius of 0.5 Mpc. We interpret the baseline density level to represent the background contamination and so estimate that 22 of the 31 new candidates are dwarf members of the group. The candidate's average color, $\langle (B-V)_0\rangle\approx 0.7$, and Sersic structural parameters are consistent with those parameters for the dwarf populations of other groups. We find that the central population of dwarfs is redder and brighter than the rest of the population. The measured faint end slope of the luminosity function, $\alpha\approx-1.33$, is steeper than that of the Local Group but consistent with published results for other groups. Such comparisons are complicated by systematic differences among different studies, but will be simpler when the KMTNet survey, which will provide homogenous data for 15 to 20 groups, is completed.

We develop a methodology to use the redshift dependence of the galaxy 2-point correlation function (2pCF) across the line-of-sight, $\xi(r_{\bot})$, as a probe of cosmological parameters. The positions of galaxies in comoving Cartesian space varies under different cosmological parameter choices, inducing a {\it redshift-dependent scaling} in the galaxy distribution. This geometrical distortion can be observed as a redshift-dependent rescaling in the measured $\xi(r_{\bot})$. We test this methodology using a sample of 1.75 billion mock galaxies at redshifts 0, 0.5, 1, 1.5, 2, drawn from the Horizon Run 4 N-body simulation. The shape of $\xi(r_{\bot})$ can exhibit a significant redshift evolution when the galaxy sample is analyzed under a cosmology differing from the true, simulated one. Other contributions, including the gravitational growth of structure, galaxy bias, and the redshift space distortions, do not produce large redshift evolution in the shape. We show that one can make use of this geometrical distortion to constrain the values of cosmological parameters governing the expansion history of the universe. This method could be applicable to future large scale structure surveys, especially photometric surveys such as DES, LSST, to derive tight cosmological constraints. This work is a continuation of our previous works as a strategy to constrain cosmological parameters using redshift-invariant physical quantities.

We present a BVI optical photometric study of the old open cluster Ruprecht 6 using the data obtained with the SMARTS 1.0 m telescope at the CTIO, Chile. Its color-magnitude diagrams show the clear existence of the main-sequence stars, of which turn-off point is located around V ~ 18.45 mag and B-V ~ 0.85 mag. Three red clump (RC) stars are identified at V = 16.00 mag, I = 14.41 mag and B-V = 1.35 mag. From the mean Ks-band magnitude of RC stars (Ks=12.39 +- 0.21 mag) in Ruprecht 6 from 2MASS photometry and the known absolute magnitudes of the RC stars (M_Ks = -1.595 +- 0.025 mag), we obtain the distance modulus to Ruprecht 6 (m-M)_0 = 13.84 +- 0.21 mag (d=5.86 +- 0.60 kpc). From the (J-K_s) and (B-V) colors of the RC stars, comparison of the (B-V) and (V-I) colors of the bright stars in Ruprecht 6 with those of the intrinsic colors of dwarf and giant stars, and the PARSEC isochrone fittings, we derive the reddening values of E(B-V) = 0.42 mag and E(V-I) = 0.60 mag. Using the PARSEC isochrone fittings onto the color-magnitude diagrams, we estimate the age and metallicity to be: log (t) =9.50 +- 0.10 (t =3.16 +- 0.82 Gyr) and [Fe/H] = -0.42 +- 0.04 dex. We present the Galactocentric radial metallicity gradient analysis for old (age > 1 Gyr) open clusters of Dias et al. catalog, which likely follow a single relation of [Fe/H] =(-0.034 +- 0.007) R_GC + (0.190 +- 0.080) (rms = 0.201) for the whole radial range or dual relation of [Fe/H] =(-0.077 +- 0.017) R_GC + (0.609 +- 0.161) (rms = 0.152) and constant ([Fe/H] ~ -0.3 dex) value, inside and outside of R_GC ~ 12 kpc, respectively. The metallicity and Galactocentric radius (13.28 +- 0.54 kpc) of Ruprecht 6 obtained in this study seem to be consistent with both of the relations.

The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics.

Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $\mu$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization of the Universe.

We present the results from a wide-field spectroscopic survey of globular clusters (GCs) in the Virgo Cluster. We obtain spectra for 201 GCs and 55 ultracompact dwarfs (UCDs) using the Hectospec on the Multiple Mirror Telescope, and derive their radial velocities. We identify 46 genuine intracluster GCs (IGCs), not associated with any Virgo galaxies, using the 3D GMM test on the spatial and radial velocity distribution.They are located at the projected distance 200 kpc $\lesssim$ R $\lesssim$ 500 kpc from the center of M87. The radial velocity distribution of these IGCs shows two peaks, one at $v_{\rm r}$ = 1023 km s$^{-1}$ associated with the Virgo main body, and another at $v_{\rm r}$ = 36 km s$^{-1}$ associated with the infalling structure. The velocity dispersion of the IGCs in the Virgo main body is $\sigma_{\rm{GC}} \sim$ 314 km s$^{-1}$, which is smoothly connected to the velocity dispersion profile of M87 GCs, but much lower than that of dwarf galaxies in the same survey field, $\sigma_{\rm{dwarf}} \sim$ 608 km s$^{-1}$. The UCDs are more centrally concentrated on massive galaxies, M87, M86, and M84. The radial velocity dispersion of the UCD system is much smaller than that of dwarf galaxies. Our results confirm the large-scale distribution of Virgo IGCs indicated by previous photometric surveys. The color distribution of the confirmed IGCs shows a bimodality similar to that of M87 GCs. This indicates that most IGCs are stripped off from dwarf galaxies and some from massive galaxies in the Virgo.

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.

We apply the methodology developed in \cite{Li2014,Li2015} to BOSS DR12 galaxies and derive cosmological constraints from the redshift dependence of the Alcock-Paczynski (AP) effect. The apparent anisotropy in the distribution of observed galaxies arise from two main sources, the redshift-space distortion (RSD) effect due to the galaxy peculiar velocities, and the geometric distortion when incorrect cosmological models are assumed for transforming redshift to comoving distance, known as the AP effect. Anisotropies produced by the RSD effect are, although large, maintaining a nearly uniform magnitude over a large range of redshift, while the degree of anisotropies from the AP effect varies with redshift by much larger magnitude. We split the DR12 galaxies into six redshift bins, measure the 2-point correlation function in each bin, and assess the redshift evolution of anisotropies. We obtain constraints of $\Omega_m=0.290 \pm 0.053,\ \ w = -1.07 \pm 0.15$, which are comparable with the current constraints from other cosmological probes such as type Ia supernovae, cosmic microwave background, and baryon acoustic oscillation (BAO). Combining these cosmological probes with our method yield tight constraints of $ \Omega_m = 0.301 \pm 0.006,\ w=-1.054 \pm 0.025$. Our method is complementary to the other large scale structure probes like BAO and topology. We expect this technique will play an important role in deriving cosmological constraints from large scale structure surveys.

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 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.

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.

The Korea Invisible Mass Search (KIMS) collaboration has developed low-background NaI(Tl) crystals that are suitable for the direct detection of WIMP dark matter. With experience built on the KIMS-CsI programs, the KIMS-NaI experiment will consist of a 200~kg NaI(Tl) crystal array surrounded by layers of shielding structures and will be operated at the Yangyang underground laboratory. The goal is to provide an unambiguous test of the DAMA/LIBRA's annual modulation signature. Measurements of six prototype crystals show progress in the reduction of internal contaminations of radioisotopes. Based on our understanding of these measurements, we expect to achieve a background level in the final detector configuration that is less than 1~count/day/keV/kg for recoil energies around 2~keV. The annual modulation sensitivity for the KIMS-NaI experiment shows that an unambiguous 7$\sigma$ test of the DAMA/LIBRA signature would be possible with a 600~kg$\cdot$year exposure with this system.

The Korea Invisible Mass Search (KIMS) collaboration has developed low-background NaI(Tl) crystals that are suitable for the direct detection of WIMP dark matter. With experience built on the KIMS-CsI programs, the KIMS-NaI experiment will consist of a 200~kg NaI(Tl) crystal array surrounded by layers of shielding structures and will be operated at the Yangyang underground laboratory. The goal is to provide an unambiguous test of the DAMA/LIBRA's annual modulation signature. Measurements of six prototype crystals show progress in the reduction of internal contaminations of radioisotopes. Based on our understanding of these measurements, we expect to achieve a background level in the final detector configuration that is less than 1~count/day/keV/kg for recoil energies around 2~keV. The annual modulation sensitivity for the KIMS-NaI experiment shows that an unambiguous 7$\sigma$ test of the DAMA/LIBRA signature would be possible with a 600~kg$\cdot$year exposure with this system.

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.

Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langle\delta\delta v v\rangle$. In the post-reionization epoch, nonlinear effects dominate in the power spectrum. We use perturbation theory and cosmological $N$-body simulations to calculate this nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular emphasis on the connected term that has been ignored in the literature. While the contribution from the connected term on small scales ($k>0.1\,h\,\rm{Mpc}^{-1}$) is subdominant relative to the unconnected term, we find that its contribution to the kSZ power spectrum at $\ell = 3000$ at $z<6$ can be as large as ten percent of the unconnected term, which would reduce the allowed contribution from the reionization epoch ($z>6$) by twenty percent. The power spectrum of transverse momentum on large scales is expected to scale as $k^2$ as a consequence of momentum conservation. We show that both the leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of the unconnected and connected terms are necessary to reproduce $k^2$.

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.

We report on the response of a high light-output NaI(Tl) crystal to nuclear recoils induced by neutrons from an Am-Be source and compare the results with the response to electron recoils produced by Compton scattered 662 keV $\gamma$-rays from a $^{137}$Cs source. The measured pulse-shape discrimination (PSD) power of the NaI(Tl) crystal is found to be significantly improved because of the high light output of the NaI(Tl) detector. We quantify the PSD power with a quality factor and estimate the sensitivity to the interaction rate for weakly interacting massive particles (WIMPs) with nucleons, and the result is compared with the annual modulation amplitude observed by the DAMA/LIBRA experiment. The sensitivity to spin-independent WIMP-nucleon interactions based on 100 kg$\cdot$year of data from NaI detectors is estimated with simulated experiments, using the standard halo model.

We report on the response of a high light-output NaI(Tl) crystal to nuclear recoils induced by neutrons from an Am-Be source and compare the results with the response to electron recoils produced by Compton scattered 662 keV $\gamma$-rays from a $^{137}$Cs source. The measured pulse-shape discrimination (PSD) power of the NaI(Tl) crystal is found to be significantly improved because of the high light output of the NaI(Tl) detector. We quantify the PSD power with a quality factor and estimate the sensitivity to the interaction rate for weakly interacting massive particles (WIMPs) with nucleons, and the result is compared with the annual modulation amplitude observed by the DAMA/LIBRA experiment. The sensitivity to spin-independent WIMP-nucleon interactions based on 100 kg$\cdot$year of data from NaI detectors is estimated with simulated experiments, using the standard halo model.

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.

Collisionless shocks can be produced as a result of strong magnetic fields in a plasma flow, and therefore are common in many astrophysical systems. The Weibel instability is one candidate mechanism for the generation of sufficiently strong fields to create a collisionless shock. Despite their crucial role in astrophysical systems, observation of the magnetic fields produced by Weibel instabilities in experiments has been challenging. Using a proton probe to directly image electromagnetic fields, we present evidence of Weibel-generated magnetic fields that grow in opposing, initially unmagnetized plasma flows from laser-driven laboratory experiments. Three-dimensional particle-in-cell simulations reveal that the instability efficiently extracts energy from the plasma flows, and that the self-generated magnetic energy reaches a few percent of the total energy in the system. This result demonstrates an experimental platform suitable for the investigation of a wide range of astrophysical phenomena, including collisionless shock formation in supernova remnants, large-scale magnetic field amplification, and the radiation signature from gamma-ray bursts.

Determining physical parameters of binary microlenses is hampered by the lack of information about the angular Einstein radius due to the difficulty of resolving caustic crossings. In this paper, we present the analysis of the binary microlensing event OGLE-2013-BLG-0578, for which the caustic exit was precisely predicted in advance from real-time analysis, enabling to densely resolve the caustic crossing and to measure the Einstein radius. From the mass measurement of the lens system based on the Einstein radius combined with the additional information about the lens parallax, we identify that the lens is a binary that is composed of a late-type M-dwarf primary and a substellar brown-dwarf companion. The event demonstrates the capability of current real-time microlensing modeling and the usefulness of microlensing in detecting and characterizing faint or dark objects in the Galaxy.

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.

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.

Among the direct search experiments for WIMP dark matter, the DAMA experiment observed an annual modulation signal interpreted as WIMP interactions with 9.2$\sigma$ significance. However, this result is contradictory with other direct search experiments reporting null signals in the same parameter space allowed by the DAMA observation, necessitating clarification of the origin of the modulation signal observed using the NaI(Tl) crystals of the DAMA experiment independently. Here, we report the first results of NaI(Tl) crystal measurement at the Yangyang Underground Laboratory to grow ultra-low-background NaI(Tl) crystal detectors.

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 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.

We present a spectroscopic study of the four extended star clusters (ESCs) in NGC 6822 based on the data obtained with Gemini Multi-Object Spectrograph (GMOS) on the Gemini-South 8.1m telescope. Their radial velocities derived from the spectra range from $-61.2 \pm 20.4$ km s$^{-1}$ (for C1) to $-115.34 \pm 57.9$ km s$^{-1}$ (for C4) and, unlike the intermediate age carbon stars, they do not display any sign of systematic rotation around NGC 6822. The ages and metallicities derived using the Lick indices show that the ESCs are old ($\geq 8$ Gyr) and metal poor ([Fe/H]$\lesssim -1.5$). NGC 6822 is found to have both metal poor ([Fe/H]$\approx -2.0$) and metal rich ([Fe/H]$\approx -0.9$) star clusters within 15 arcmin (2 kpc) from the center, while only metal poor clusters are observed in the outer halo with $r \geq 20$ arcmin (2.6 kpc). The kinematics, old ages, and low metallicities of ESCs suggest that ESCs may have accreted into the halo of NGC 6822. Based on the velocity distribution of ESCs, we have determined the total mass and the mass-to-light ratio of NGC 6822: $M_{N6822} = 7.5^{+4.5}_{-0.1} \times 10^{9} M_{\odot}$ and $(M/L)_{N6822} = 75^{+45}_{-1} (M/L)_{\odot}$. It shows that NGC 6822 is one of the most dark matter dominated dwarf galaxies in the Local Group.

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.

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 present a determination of the two-dimensional shape parameters of the blue and red globular cluster systems (GCSs) in a large number of elliptical galaxies and lenticular galaxies (early-type galaxies, called ETGs). We use a homogeneous data set of the globular clusters in 23 ETGs obtained from the HST/ACS Virgo Cluster Survey. The position angles of both blue and red GCSs show a correlation with those of the stellar light distribution, showing that the major axes of the GCSs are well aligned with those of their host galaxies. However, the shapes of the red GCSs show a tight correlation with the stellar light distribution as well as with the rotation property of their host galaxies, while the shapes of the blue GCSs do much less. These provide clear geometric evidence that the origins of the blue and red globular clusters are distinct and that ETGs may have dual halos: a blue (metal-poor) halo and a red (metal-rich) halo. These two halos show significant differences in metallicity, structure, and kinematics, indicating that they are formed in two distinguishable ways. The red halos might have formed via dissipational processes with rotation, while the blue halos are through accretion.

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.

Supersonic plasma outflows driven by multi-beam, high-energy lasers, such as Omega and NIF, have been and will be used as platforms for a variety of laboratory astrophysics experiments. Here we propose a new way of launching high density and high velocity, plasma jets using multiple intense laser beams in a hollow ring formation. We show that such jets provide a more flexible and versatile platform for future laboratory astrophysics experiments. Using high resolution hydrodynamic simulations, we demonstrate that the collimated jets can achieve much higher density, temperature and velocity when multiple laser beams are focused to form a hollow ring pattern at the target, instead of focused onto a single spot. We carried out simulations with different ring radii and studied their effects on the jet properties. Implications for laboratory collisionless shock experiments are discussed.

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 study chromospheric oscillations including umbral flashes and running penumbral waves in a sunspot of active region (AR) using scanning spectroscopy in H-alpha and Ca II 8542A, with the Fast Imaging Solar Spectrograph (FISS) at the 1.6 meter New Solar Telescope at Big Bear Solar Observatory. A bisector method is applied to spectral observations to construct chromospheric Doppler velocity maps. Temporal sequence analysis of these shows enhanced high-frequency oscillations inside the sunspot umbra in both lines. Their peak frequency gradually decreases outward from the umbra. The oscillation power is found to be associated with magnetic-field strength and inclination, with different relationships in different frequency bands.

We calculate the angular power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations induced by the kinetic Sunyaev-Zel'dovich (kSZ) effect from the epoch of reionization (EOR). We use detailed N-body+radiative transfer simulations to follow inhomogeneous reionization of the intergalactic medium (IGM). For the first time we take into account the "self-regulation" of reionization: star formation in low-mass dwarf galaxies (10^8 M_\sun \lesssim M \lesssim 10^9 M_\sun) or minihalos (10^5 M_\sun \lesssim M \lesssim 10^8 M_\sun) is suppressed if these halos form in the regions that were already ionized or Lyman-Werner dissociated. Some previous work suggested that the amplitude of the kSZ power spectrum from the EOR can be described by a two-parameter family: the epoch of half ionization and the duration of reionization. However, we argue that this picture applies only to simple forms of the reionization history which are roughly symmetric about the half-ionization epoch. In self-regulated reionization, the universe begins to be ionized early, maintains a low level of ionization for an extended period, and then finishes reionization as soon as high-mass atomically-cooling halos dominate. While inclusion of self-regulation affects the amplitude of the kSZ power spectrum only modestly (\sim 10 %), it can change the duration of reionization by a factor of more than two. We conclude that the simple two-parameter family does not capture the effect of a physical, yet complex, reionization history caused by self-regulation. When added to the post-reionization kSZ contribution, our prediction for the total kSZ power spectrum is below the current upper bound from the South Pole Telescope. Therefore, the current upper bound on the kSZ effect from the EOR is consistent with our understanding of the physics of reionization.

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.

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

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 present a spectroscopic analysis of the metallicities, ages, and alpha-elements of the globular clusters (GCs) in the giant elliptical galaxy (gE) NGC 4636 in the Virgo cluster. Line indices of the GCs are measured from the integrated spectra obtained with Faint Object Camera and Spectrograph (FOCAS) on the Subaru 8.2 m telescope. We derive [Fe/H] values of 59 GCs based on the Brodie & Huchra method, and [Z/H], age, and [a/Fe] values of 33 GCs from the comparison of the Lick line indices with single stellar population models. The metallicity distribution of NGC 4636 GCs shows a hint of a bimodality with two peaks at [Fe/H]=-1.23 (sigma=0.32) and -0.35 (sigma=0.19). The age spread is large from 2 Gyr to 15 Gyr and the fraction of young GCs with age < 5 Gyr is about 27%. The [a/Fe] of the GCs shows a broad distribution with a mean value [a/Fe]~0.14 dex. The dependence of these chemical properties on the galactocentric radius is weak. We also derive the metallicities, ages, and [a/Fe] values for the GCs in other nearby gEs (M87, M49, M60, NGC 5128, NGC 1399, and NGC 1407) from the line index data in the literature using the same methods as used for NGC 4636 GCs. The metallicity distribution of GCs in the combined sample of seven gEs including NGC 4636 is found to be bimodal, supported by the KMM test with a significance level of >99.9%. All these gEs harbor some young GCs with ages less than 5 Gyr. The mean age of the metal-rich GCs ([Fe/H] > -0.9) is about 3 Gyr younger than that of the metal-poor GCs. The mean value of [a/Fe] of the gE GCs is smaller than that of the Milky Way GCs. We discuss these results in the context of GC formation in gEs.

We present the first spectroscopic study of the globular clusters (GCs) in the giant elliptical galaxy (gE) M86 in the Virgo cluster. Using spectra obtained in the Multi-Object Spectroscopy (MOS) mode of the Faint Object Camera and Spectrograph (FOCAS) on the Subaru telescope, we measure the radial velocities for 25 GCs in M86. The mean velocity of the GCs is derived to be vp = -354+81 -79 km/s, which is different from the velocity of the M86 nucleus (vgal = -234\pm41 km/s). We estimate the velocity dispersion of the GCs, {\sigma}p = 292+32-32 km/s, and find a hint of rotation of the M86 GC system. A comparison of the observed velocity dispersion profiles of the GCs and stars with a prediction based on the stellar mass profile strongly suggests the existence of an extended dark matter halo in M86. We also estimate the metallicities and ages for 16 and 8 GCs, respectively. The metallicities of M86 GCs are in the range -2.0 < [Fe/H] < -0.2 with a mean value of -1.13\pm0.47. These GCs show a wide age distribution from 4 to 15 Gyr.

We present a photometric study of the globular clusters (GCs) in the Virgo giant elliptical galaxy M86 based on Washington CT1 images. The colors of the GCs in M86 show a bimodal distribution with a blue peak at (C -T1) = 1.30 and a red peak at (C -T1) = 1.72. The spatial distribution of the red GCs is elongated similarly to that of the stellar halo, while that of the blue GCs is roughly circular. The radial number density profile of the blue GCs is more extended than that of the red GCs. The radial number density profile of the red GCs is consistent with the surface brightness profile of the M86 stellar halo. The GC system has a negative radial color gradient, which is mainly due to the number ratio of the blue GCs to the red GCs increasing as galactocentric radius increase. The bright blue GCs in the outer region of M86 show a blue tilt: the brighter they are, the redder their mean colors get. These results are discussed in comparison with other Virgo giant elliptical galaxies.

We report the discovery of two new globular clusters in the remote halos of M81 and M82 in the M81 Group based on Hubble Space Telescope archive images. They are brighter than typical globular clusters (MV = -9.34 mag for GC-1 and M_V = -10.51 mag for GC-2), and much larger than known globular clusters with similar luminosity in the MilkyWay Galaxy and M81. Radial surface brightness profiles for GC-1 and GC-2 do not show any features of tidal truncation in the outer part. They are located much farther from both M81 and M82 in the sky, compared with previously known star clusters in these galaxies. Color-magnitude diagrams of resolved stars in each cluster show a well-defined red giant branch (RGB), indicating that they are metal-poor and old. We derive a low metallicity with [Fe/H] $\simeq -2.3$ and an old age ~14 Gyr for GC-2 from the analysis of the absorption lines in its spectrum in the Sloan Digital Sky Survey in comparison with the simple stellar population models. The I-band magnitude of the tip of the RGB for GC-2 is 0.26 mag fainter than that for the halo stars in the same field, showing that GC-2 is ~400 kpc behind the M81 halo along our line of sight. The deprojected distance to GC-2 from M81 is much larger than any other known globular clusters in the local universe. This shows that GC-2 is the most isolated globular cluster in the local universe.

We present a photometric estimation of the distance and reddening values to the dwarf irregular galaxy NGC 1156, which is one of the best targets to study the isolated dwarf galaxies in the nearby universe. We have used the imaging data sets of the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) High Resolution Channel (HRC) of the central region of NGC 1156 (26" X 29") available in the HST archive for this study. From the (U-B, B-V) color-color diagram, we first estimate the total (foreground + internal) reddening toward NGC 1156 of E(B-V) =0.35 +/- 0.05 mag, whereas only the foreground reddening was previously known to be E(B-V)=0.16 mag (Burstein & Heiles) or 0.24 mag (Schlegel, Finkbeiner, & Davis). Based on the brightest stars method, selecting the three brightest blue supergiant (BSG) stars with mean B magnitude of <B(3B)> = 21.94 mag and the three brightest red supergiant (RSG) stars with mean V magnitude of <V(3R)> = 22.76 mag, we derive the distance modulus to NGC 1156 to be (m-M)_{0,BSG} = 29.55 mag and (m-M)_{0,RSG} = 29.16 mag. By using weights of 1 and 1.5 for the distance moduli from using the BSGs and the RSGs, respectively, we finally obtain the weighted mean distance modulus to NGC 1156 (m-M)_0 = 29.39 +/- 0.20 mag (d = 7.6 +/- 0.7 Mpc), which is in very good agreement with the previous estimates. Combining the photometry data of this study with those of Karachentsev et al. gives smaller distance to NGC 1156, which is discussed together with the limits of the data.

A pixel analysis is carried out on the interacting face-on spiral galaxy NGC 5194 (M51A), using the HST/ACS images in the F435W, F555W and F814W (BVI) bands. After 4x4 binning of the HST/ACS images to secure a sufficient signal-to-noise ratio for each pixel, we derive several quantities describing the pixel color-magnitude diagram (pCMD) of NGC 5194: blue/red color cut, red pixel sequence parameters, blue pixel sequence parameters and blue-to-red pixel ratio. The red sequence pixels are mostly older than 1 Gyr, while the blue sequence pixels are mostly younger than 1 Gyr, in their luminosity-weighted mean stellar ages. The color variation in the red pixel sequence from V = 20 mag arcsec^(-2) to V = 17 mag arcsec^(-2) corresponds to a metallicity variation of \Delta[Fe/H] ~ 2 or an optical depth variation of \Delta\tau_V ~ 4 by dust, but the actual sequence is thought to originate from the combination of those two effects. At V < 20 mag arcsec^(-2), the color variation in the blue pixel sequence corresponds to an age variation from 5 Myr to 300 Myr under the assumption of solar metallicity and \tau_V = 1. To investigate the spatial distributions of stellar populations, we divide pixel stellar populations using the pixel color-color diagram and population synthesis models. As a result, we find that the pixel population distributions across the spiral arms agree with a compressing process by spiral density waves: dense dust \rightarrow newly-formed stars. The tidal interaction between NGC 5194 and NGC 5195 appears to enhance the star formation at the tidal bridge connecting the two galaxies. We find that the pixels corresponding to the central active galactic nucleus (AGN) area of NGC 5194 show a tight sequence at the bright-end of the pCMD, which are in the region of R ~ 100 pc and may be a photometric indicator of AGN properties.

We present a study on four new star clusters discovered in the halo of the intriguing dwarf irregular galaxy NGC 6822 from a wide field survey covering 3 deg x 3 deg area carried out with MegaCam at Canada-France-Hawaii Telescope (CFHT). The star clusters have extended structures with half-light radii R_h ~ 7.5 -- 14.0 pc, larger than typical Galactic globular clusters and other known globular clusters in NGC 6822. The integrated colors and color magnitude diagrams (CMD) of resolved stars suggest that the new star clusters are 2 -- 10 Gyr old and relatively metal poor with Z=0.0001--0.004 based on the comparison with theoretical models. The projected distance of each star cluster from the galaxy center ranges from 10.7 arcmin (~ 1.5 kpc) to 77 arcmin (~ 11 kpc), far beyond the optical body of the galaxy. Interestingly, the new star clusters are aligned along the elongated old stellar halo of NGC 6822, which is almost perpendicular to the HI gas distribution where young stellar populations exist. We also find that the colors and half-light radii of the new clusters are correlated with the galactocentric distance: clusters farther from the galaxy center are larger and bluer than those closer to the galaxy center. We discuss the stellar structure and evolution of NGC 6822 implied by these new extended star clusters in the halo. We also discuss the current status of observational and theoretical understandings regarding the origin of extended star clusters in NGC 6822 and other galaxies.

In the last decade we witness an advent of new types of dwarf stellar systems in cluding ultra-compact dwarfs, ultra-faint dwarf spheroidals, and exotic globular clusters, breaking the old simple paradigm for dwarf galaxies and globular clusters. These objects become more intriguing, and understanding of these new findings be comes more challenging. Recently we discovered a new type of large scale structure in the Virgo cluster of galaxies: it is composed of globular clusters. Globular clusters in Virgo are found wandering between galaxies (intracluster globular clusters) as well as in galaxies. These intracluster globular clusters fill a significant fraction in the area of the Virgo cluster and they are dominated by blue globular clusters. These intracluster globular clusters may be closely related with the first dwarf galaxies in the universe.

We extend the concept of galaxy environment from the local galaxy number density to the gravitational potential and its functions like the shear tensor. For this purpose we examine whether or not one can make an accurate estimation of the gravitational potential from an observational sample which is finite in volume, biased due to galaxy biasing, and subject to redshift space distortion. Dark halos in a $\Lambda$CDM simulation are used in this test. We find that one needs to stay away from the sample boundaries by more than 30$h^{-1}$Mpc to reduce the error within 20% of the root mean square values of the potential or the shear tensor. The error due to the galaxy biasing can be significantly reduced by using the galaxy mass density field instead of the galaxy number density field. The error caused by the redshift space distortion can be effectively removed by correcting galaxy positions for the peculiar velocity effects. We inspect the dependence of dark matter halo properties on four environmental parameters; local density, gravitational potential, and the ellipticity and prolateness of the shear tensor. We find the local density has the strongest correlation with halo properties. This is evidence that the internal physical properties of dark halos are mainly controlled by small-scale physics. In high density regions dark halos are on average more massive and spherical, and have higher spin parameter and velocity dispersion. In high density regions dark halos are on average more massive and spherical, and have higher spin parameter and velocity dispersion. We also study the relation between the environmental parameters and the subtypes of dark halos. The spin parameter of satellite halos depends only weakly on the local density for all mass ranges studied while that of isloated or central halos depends more sensitively on the local density. (abridged)

The Cosmic Ray Energetics And Mass (CREAM) balloon experiment had two successful flights in 2004/05 and 2005/06. It was designed to perform energy measurements from a few GeV up to 1000 TeV, taking advantage of different detection techniques. The first instrument, CREAM-1, combined a transition radiation detector with a calorimeter to provide independent energy measurements of cosmicraynuclei. Each detector was calibrated with particle beams in a limited range of energies. In order to assess the absolute energy scale of the instrument and to investigate the systematic effects of each technique, a cross-calibration was performed by comparing the two independent energy estimates on selected samples of oxygen and carbon nuclei.

The CREAM calorimeter, designed to measure the spectra of cosmic-ray nuclei from under 1 TeV to 1000 TeV, is a 20 radiation length (X0) deep sampling calorimeter. The calorimeter is comprised of 20 layers of tungsten interleaved with 20 layers of scintillating fiber ribbons, and is preceded by a pair of graphite interaction targets providing about 0.42 proton interaction lengths (\lambda int). The calorimeter was placed in one of CERN's SPS accelerator beams for calibration and testing. Beams of 150 GeV electrons were used for calibration, and a variety of electron, proton, and nuclear fragment beams were used to test the simulation model of the detector. In this paper we discuss the performance of the calorimeter in the electron beam and compare electron beam data with simulation results.

Globular clusters are found usually in galaxies and they are an excellent tracer of dark matter. Long ago it was suggested that there may exist intracluster globular clusters (IGCs) bound to a galaxy cluster rather than to any single galaxy. Here we present a map showing the large scale distribution of globular clusters over the entire Virgo cluster. It shows that IGCs are found out to 5 million light years from the Virgo center, and that they are concentrated in several substructures much larger than galaxies. These objects might have been mostly stripped off from low-mass dwarf galaxies.

We present a kinematic analysis of the globular cluster(GC) system in the giant elliptical galaxy (gE) NGC 4636 in the Virgo cluster. Using the photometric and spectroscopic database of 238 GCs, we have investigated the kinematics of the GC system. The NGC 4636 GC system shows weak overall rotation, which is dominated by the red GCs. However, both the blue GCs and red GCs show some rotation in the inner region at R<4.3'. The velocity dispersion for all the GCs is derived to be sigma_p = 225{+12-9} km/s. The velocity dispersion for the blue GCs (sig=251 km/s) is slightly larger than that for the red GCs (sig=205 km/s). The velocity dispersions for the blue GCs about the mean velocity and about the best fit rotation curve have a significant variation depending on the galactocentric radius. Comparison of observed stellar and GC velocity dispersion profiles with the velocity dispersion profiles calculated from the stellar mass profile shows that the mass-to-light ratio should increase as the galactocentric distance increases, indicating the existence of an extended dark matter halo. From the comparison of the observed GC velocity dispersion profiles and the velocity dispersion profiles calculated for the X-ray mass profiles in the literature, we find that the orbit of the GC system is tangential, and that the orbit of the red GCs is slightly more tangential than that of the blue GCs. We compare the GC kinematics of NGC 4636 with those of other six gEs, finding that the kinematic properties of the GCs are diverse among gEs. We find several correlations between the kinematics of the GCs and the global parameters of their host galaxies. We discuss the implication of the results for the formation models of the GC system in gEs, and suggest a mixture scenario for the origin of the GCs in gEs.

We present a spectroscopic study of the globular clusters (GCs) in the giant elliptical galaxy NGC 4636 in the Virgo cluster. We selected target GC candidates using the Washington photometry derived from the deep CCD images taken at the KPNO 4m. Then we obtained the spectra of 164 target objects in the field of NGC 4636 using the Multi-Object Spectroscopy (MOS) mode of Faint Object Camera and Spectrograph (FOCAS) on the SUBARU 8.2m Telescope. We have measured the velocities for 122 objects: 105 GCs in NGC 4636, the nucleus of NGC 4636, 11 foreground stars, 2 background galaxies, and 3 probable intracluster GCs in the Virgo cluster. The GCs in NGC 4636 are located in the projected galactocentric radius within 10arcmin (corresponding to 43 kpc). The measured velocities for the GCs range from 300km/s to 1600km/s, with a mean value of 932_{-22}^{+25} km/s, which is in good agreement with the velocity for the nucleus of NGC 4636, 928\pm 45 km/s. The velocity dispersion of the GCs in NGC 4636 is derived to be 231_{-17}^{+15} km/s and the velocity dispersion of the blue GCs is slightly larger than that of the red GCs. Combining our results with data in the literature, we produce a master catalog of radial velocities for 238 GCs in NGC 4636. The velocity dispersion of the GCs in the master catalog is found to be 225_{-9}^{+12} km/s for the entire sample, 251_{-12}^{+18} km/s for 108 blue GCs, and 205_{-13}^{+11} km/s for 130 red GCs.

We present the star formation history and chemical evolution of the Sextans dSph dwarf galaxy as a function of galactocentric distance. We derive these from the $VI$ photometry of stars in the $42' \times 28'$ field using the SMART model developed by Yuk & Lee (2007, ApJ, 668, 876) and adopting a closed-box model for chemical evolution. For the adopted age of Sextans 15 Gyr, we find that $>$84% of the stars formed prior to 11 Gyr ago, significant star formation extends from 15 to 11 Gyr ago ($\sim$ 65% of the stars formed 13 to 15 Gyr ago while $\sim$ 25% formed 11 to 13 Gyr ago), detectable star formation continued to at least 8 Gyr ago, the star formation history is more extended in the central regions than the outskirts, and the difference in star formation rates between the central and outer regions is most marked 11 to 13 Gyr ago. Whether blue straggler stars are interpreted as intermediate age main sequence stars affects conclusions regarding the star formation history for times 4 to 8 Gyr ago, but this is at most only a trace population. We find that the metallicity of the stars increased rapidly up to [Fe/H]=--1.6 in the central region and to [Fe/H]=--1.8 in the outer region within the first Gyr, and has varied slowly since then. The abundance ratios of several elements derived in this study are in good agreement with the observational data based on the high resolution spectroscopy in the literature. We conclude that the primary driver for the radial gradient of the stellar population in this galaxy is the star formation history, which self-consistently drives the chemical enrichment history.

We present a result of age estimation for star clusters in M33. We obtain color-magnitude diagrams (CMDs) of resolved stars in 242 star clusters from the HST/WFPC2 images. We estimate ages of 100 star clusters among these, by fitting the Padova theoretical isochrones to the observational CMDs. Age distribution of the star clusters shows a dominant peak at log(t) ~ 7.8. Majority of star clusters are younger than log(t) = 9.0, while ten star clusters are older than log(t) ~ 9.0. There are few clusters younger than log(t) = 7 in this study, which is in contrast with the results based on the integrated photometry of star clusters in the previous studies. Radial distribution of the cluster ages shows that young to intermediate-age clusters are found from the center to the outer region, while old clusters are distributed farther from M33 center. We discuss briefly the implication of the results with regard to the formation of M33 cluster system.

We provide an overview of the robotic Super-LOTIS (Livermore Optical Transient Imaging System) telescope and present results from gamma-ray burst (GRB) afterglow observations using Super-LOTIS and other Steward Observatory telescopes. The 0.6-m Super-LOTIS telescope is a fully robotic system dedicated to the measurement of prompt and early time optical emission from GRBs. The system began routine operations from its Steward Observatory site atop Kitt Peak in April 2000 and currently operates every clear night. The telescope is instrumented with an optical CCD camera and a four position filter wheel. It is capable of observing Swift Burst Alert Telescope (BAT) error boxes as early or earlier than the Swift UV/Optical Telescope (UVOT). Super-LOTIS complements the UVOT observations by providing early R- and I-band imaging. We also use the suite of Steward Observatory telescopes including the 1.6-m Kuiper, the 2.3-m Bok, the 6.5-m MMT, and the 8.4-m Large Binocular Telescope to perform follow-up optical and near infrared observations of GRB afterglows. These follow-up observations have traditionally required human intervention but we are currently working to automate the 1.6-m Kuiper telescope to minimize its response time.

We present a photometric study of the globular clusters in the giant elliptical galaxy M60 in the Virgo cluster, based on deep, relatively wide field Washington CT_1 CCD images. The color-magnitude diagram reveals a significant population of globular clusters in M60, and a large number of young luminous clusters in NGC 4647, a small companion spiral galaxy north-west of M60. The color distribution of the globular clusters in M60 is clearly bimodal, with a blue peak at (C-T_1)=1.37, and a red peak at (C-T_1)=1.87. We derive two new transformation relations between the (C-T_1)_0 color and [Fe/H] using the data for the globular clusters in our Galaxy and M49. Using these relations we derive the metallicity distribution of the globular clusters in M60, which is also bimodal: a dominant metal-poor component with center at [Fe/H]=-1.2, and a weaker metal-rich component with center at [Fe/H]=-0.2. The radial number density profile of the globular clusters is more extended than that of the stellar halo, and the radial number density profile of the blue globular clusters is more extended than that of the red globular clusters. The number density maps of the globular clusters show that the spatial distribution of the blue globular clusters is roughly circular, while that of the red globular cluster is elongated similarly to that of the stellar halo. We estimate the total number of the globular clusters in M60 to be 3600+/-500$,and the specific frequency to be S_N=3.8+/-0.4. The mean color of the bright blue globular clusters gets redder as they get brighter in both the inner and outer region of M60. This blue tilt is seen also in the outer region of M49, the brightest Virgo galaxy. Implications of these results are discussed.

We present a kinematic analysis of the globular cluster(GC) system in M31. Using the photometric and spectroscopic database of 504 GCs, we have investigated the kinematics of the M31 GC system. We find that the all GC system shows strong rotation, with rotation amplitude of v_rot~190km/s, and that a weak rotation persists even for the outermost samples at |Y|>5kpc. The rotation-corrected velocity dispersion for the GC system is estimated to be sigma_{p,r}~130km/s, and it increases from sigma_{p,r}~120km/s at |Y|<1kpc to sigma_{p,r}~150km/s at |Y|>5kpc. These results are very similar to those for the metal-poor GCs. This shows that there is a dynamically hot halo in M31 that is rotating but primarily pressure-supported. We have identified 50 "friendless" GCs, and they appear to rotate around the major axis of M31. For the subsamples of metal-poor and metal-rich GCs, we have found that the metal-rich GCs are more centrally concentrated than the metal-poor GCs, and both subsamples show strong rotation. For the subsamples of bright and faint GCs, it is found that the rotation for the faint GCs is stronger than that for the bright GCs. We have identified 56 GCs and GC candidates with X-ray detection. It is found that the majority of X-ray emitting GCs follow the disk rotation, and that the redder, more metal-rich, and brighter GCs are more likely to be detected as X-ray emitting GCs. We have derived a rotation curve of M31 using the GCs at |Y|<0.6kpc. We have estimated the dynamical mass of M31 using `Projected Mass Estimator(PME)' and `Tracer Mass Estimator(TME)'. We finally discuss the implication of these results and compare the kinematics of GCs with that of planetary nebulae in M31.

We present the measurement of radial velocities for globular clusters in M60, giant elliptical galaxy in the Virgo cluster. Target globular cluster candidates were selected using the Washington photometry based on the deep $16\arcmin \times 16\arcmin$ images taken at the KPNO 4m and using the $VI$ photometry derived from the HST/WFPC2 archive images. The spectra of the target objects were obtained using the Multi-Object Spectrograph (MOS) at the Canada-France-Hawaii Telescope (CFHT). We have measured the radial velocity for 111 objects in the field of M60: 93 globular clusters (72 blue globular clusters with $1.0\le(C-T_1)<1.7$ and 21 red globular clusters with $1.7\le(C-T_1)<2.4$), 11 foreground stars, 6 small galaxies, and the nucleus of M60. The measured velocities of the 93 globular clusters range from $\sim 500$ km s$^{-1}$ to $\sim 1600$ km s$^{-1}$, with a mean value of $1070_{-25}^{+27}$ km s$^{-1}$, which is in good agreement with the velocity of the nucleus of M60 ($v_{\rm gal}=1056$ km s$^{-1}$). Combining our results with data in the literature, we present a master catalog of radial velocities for 121 globular clusters in M60. The velocity dispersion of the globular clusters in the master catalog is found to be $234_{-14}^{+13}$ km s$^{-1}$ for the entire sample, $223_{-16}^{+13}$ km s$^{-1}$ for 83 blue globular clusters, and $258_{-31}^{+21}$ km s$^{-1}$ for 38 red globular clusters.

We present a kinematic analysis of the globular cluster (GC) system in the giant elliptical galaxy (gE) M60 in the Virgo cluster. Using the photometric and spectroscopic database of 121 GCs (83 blue GCs and 38 red GCs), we have investigated the kinematics of the GC system. We have found that the M60 GC system shows a significant overall rotation. The rotation amplitude of the blue GCs is slightly smaller than or similar to that of the red GCs, and their angles of rotation axes are similar. The velocity dispersions about the mean velocity and about the best fit rotation curve for the red GCs are marginally larger than those for the blue GCs. Comparison of observed stellar and GC velocity dispersion profiles with those calculated from the stellar mass profile shows that the mass-to-light ratio should be increased as the galactocentric distance increases, indicating the existence of an extended dark matter halo. The entire sample of GCs in M60 is found to have a tangentially biased velocity ellipsoid unlike the GC systems in other gEs. Two subsamples appear to have different velocity ellipsoids. The blue GC system has a modest tangentially biased velocity ellipsoid, while the red GC system has a modest radially biased or an isotropic velocity ellipsoid. From the comparison of the kinematic properties of the M60 GC system to those of other gEs (M87, M49, NGC 1399, NGC 5128, and NGC 4636), it is found that the velocity dispersion of the blue GC system is similar to or larger than that of the red GC system except for M60, and the rotation of the GC system is not negligible. The entire sample of each GC system shows an isotropic velocity ellipsoid except for M60, while the subsamples show diverse velocity ellipsoids. We discuss the implication of these results for the formation models of the GC system in gEs.

We work out the effects of hypercritical accretion, which transfers mass from the secondary to the primary (first-born) neutron star (NS) in a binary, showing that the mass of the primary would end up 0.6 M_sun greater than the secondary NS in contradiction with the very nearly equal masses of the two NS's in binaries measured to date. We discuss that the primary NS evolves into a low-mass black-hole (LMBH) due to the hypercritical accretion. Using a flat distribution in the mass ratio q (q = M_secondary/M_primary), we calculated a ratio of LMBH-NS and NS-NS systems to be 5, in rough agreement with Pinsonneault and Stanek (2006). These authors emphasize the importance of "twins", which we discuss. The two NS's in the twins would be close in mass and further increase the number of mergings.

We present the result of a wide-field survey of globular clusters (GCs) in M31 covering a 3deg x 3deg field c. We have searched for GCs on CCD images taken with Washington CMT1 filters at the KPNO 0.9 m telescope using steps: (1) inspection of morphological parameters given by the SExtractor package such as stellarity, full maximum, and ellipticity; (2) consulting the spectral types and radial velocities obtained from spectra takena spectrograph at the WIYN 3.5 m telescope; and (3) visual inspection of the images of each object. We have and GC candidates, of which 605 are newly found GCs and GC candidates and 559 are previously known GCs. Amoects there are 113 genuine GCs, 258 probable GCs, and 234 possible GCs, according to our classification critee known objects there are 383 genuine GCs, 109 probable GCs, and 67 possible GCs. In total there are 496 genprobable GCs and 301 possible GCs. Most of these newly found GCs have T1 magnitudes of 17.5 - 19.5 mag, [17.9 < V < 19.9 mag assuming (C-T1) ~ 1.5], and (C-T1) colors in the range 1 - 2.

We present detailed optical, near-infrared, and radio observations of the X-ray flash 050416a obtained with Palomar and Siding Springs Observatories as well as HST and the VLA, placing this event among the best-studied X-ray flashes to date. In addition, we present an optical spectrum from Keck LRIS from which we measure the redshift of the burst, z=0.6528. At this redshift the isotropic-equivalent prompt energy release was about 10^51 erg, and using a standard afterglow synchrotron model we find that the blastwave kinetic energy is a factor of 10 larger, E_K,iso ~ 10^52 erg. The lack of an observed jet break to t ~ 20 days indicates that the opening angle is larger than 7 deg and the total beaming-corrected relativistic energy is larger than 10^50 erg. We further show that the burst produced a strong radio flare at t ~ 40 days accompanied by an observed flattening in the X-ray band which we attribute to an abrupt circumburst density jump or an episode of energy injection (either from a refreshed shock or off-axis ejecta). Late-time observations with HST show evidence for an associated supernova with peak optical luminosity roughly comparable to that of SN 1998bw. Next, we show that the host galaxy of XRF 050416a is actively forming stars at a rate of at least 2 M_solar per year with a luminosity of L_B ~ 0.5L* and metallicity of Z ~ 0.2-0.8 Z_solar. Finally, we discuss the nature of XRF 050416a in the context of short-hard gamma-ray bursts and under the framework of off-axis and dirty fireball models for X-ray flashes.

We present Follow-Up Network for Gamma-Ray Bursts (FUN GRB) Collaboration observations of the optical afterglow of GRB 021211 made between 143 seconds and 102 days after the burst. Our unique data set includes the earliest filtered detections and color information for an afterglow in the pre-Swift era. We find that the afterglow is best described by (1) propagation through a wind-swept medium, (2) a cooling break that is blueward of the observed optical frequencies, and (3) a hard electron energy distribution. However, superimposed on this "standard model" behavior we find significant chromatic variations within the first few hours after the burst. We consider possible reasons for these variations, including the possibility that they are due to a dust echo. Finally, we constrain physical parameters that describe the afterglow and surrounding medium for a variety of scenarios and find that GRB 021211's afterglow is faint for a combination of 3-4 reasons: (1) a low fraction of energy in relativistic electrons, (2) a low density for the wind-swept medium, implying either a low mass-loss rate and/or a high wind velocity for the progenitor, (3) a wide opening/viewing angle for the jet, and possibly (4) moderate source frame extinction. The jet appears to be significantly far from equipartition and magnetically dominated. More extreme versions of this might explain the darkness of many afterglows in the Swift era.

We present the results from a {\it Hubble Space Telescope/ACS} study of the supernovae associated with gamma-ray bursts 040924 ($z=0.86$) and 041006 ($z=0.71$). We find evidence that both GRBs were associated with a SN 1998bw-like supernova dimmed by $\sim 1.5$ and $\sim 0.3$ magnitudes, respectively, making GRB 040924 the faintest GRB-associated SN ever detected. We study the luminosity dispersion in GRB/XRF-associated SNe and compare to local Type Ibc supernovae from the literature. We find significant overlap between the two samples, suggesting that GRB/XRF-associated SNe are not necessarily more luminous nor produce more $^{56}$Ni than local SNe. Based on the current (limited) datasets, we find that the two samples may share a similar $^{56}$Ni production mechanism.

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.

We present BV CCD photometry for the open clusters Czernik 24 and Czernik 27. These clusters have never been studied before, and we provide, for the first time, the cluster parameters; reddening, distance, metallicity and age. Czernik 24 is an old open cluster with age 1.8 +/- 0.2 Gyr, metallicity [Fe/H]=-0.41 +/- 0.15 dex, distance modulus (m-M)_0 = 13.1 +/- 0.3 mag (d=4.1 +/- 0.5 kpc), and reddening E(B-V) = 0.54 +/- 0.12 mag. The parameters for Czernik 27 are estimated to be age = 0.63 +/- 0.07 Gyr, [Fe/H]= -0.02 +/- 0.10 dex, (m-M)_0 = 13.8 +/- 0.2 mag (d=5.8 +/- 0.5 kpc), and E(B-V) = 0.15 +/- 0.05 mag. The metallicity and distance values for Czernik 24 are consistent with the relation between the metallicity and the Galactocentric distance of other old open clusters. We find the metallicity gradient of 51 old open clusters including Czernik 24 to be Delta [Fe/H]/Delta R_gc= -0.064 +/- 0.009 dex/kpc.

We present a systematic study of the low mass X-ray binary (LMXB) populations of 6 elliptical galaxies. We utilize Chandra archival data to identify 665 X-ray point sources and HST archival data supplemented by ground observations to identify 6173 GCs. We associate 209 LMXBs with red GC (RGC) and 76 LMXBs with blue GCs (BGC), while we find no optical GC counterpart for 258 LMXBs. This is the largest GC-LMXB sample studied so far. We confirm previous reports suggesting that the fraction of GCs associated with LMXBs is 3 times larger in RGCs than in BGCs, indicating that metallicity is a primary factor in the GC-LMXB formation. While as already known, the brighter (and bigger) GCs have a higher probability to host LMXBs, we find that this optical luminosity (or mass) dependency is stronger in RGCs than in BGCs. We also find that GCs located near the galaxy center have a higher probability to harbor LMXBs compared to those in the outskirts. The radial distributions of GC-LMXBs (for both RGC and BGC) are steeper than those of the whole optical GC sample, but consistent with those of the optical halo light, suggesting that there must be another parameter (in addition to metallicity) governing LMXB formation in GCs. This second parameter must depend on the galacto-centric distance. We find no statistically significant difference in the X-ray properties among RGC-LMXBs, BGC-LMXBs and field-LMXBs. The similarity of the X-ray spectra of BGC-LMXBs and RGC-LMXBs is inconsistent with the irradiation-induced stellar wind model prediction. The similarity of the X-ray luminosity functions (XLFs) of GC-LMXBs and field-LMXBs indicates that there is no significant difference in the fraction of BH binaries present in these two populations, in contrast to Galactic LMXBs where BH binaries are not found in GCs.

We report the discovery of three new star clusters in the halo of the Local Group dwarf irregular galaxy NGC 6822. These clusters were found in the deep images taken with the MegaPrime at the CFHT, covering a total field of 2 deg x 2 deg. The most remote cluster is found to be located as far as 79 arcmin away from the center of NGC 6822. This distance is several times larger than the size of the region in NGC 6822 where star clusters were previously found. Morphological structures of the clusters and color-magnitude diagrams of the resolved stars in the clusters show that at least two of these clusters are proabably old globular clusters.

The final chapter in the long-standing mystery of the gamma-ray bursts (GRBs) centres on the origin of the short-hard class, suspected on theoretical grounds to result from the coalescence of neutron star or black hole binary systems. Numerous searches for the afterglows of short-hard bursts have been made, galvanized by the revolution in our understanding of long-duration GRBs that followed the discovery in 1997 of their broadband (X-ray, optical, and radio) afterglow emission. Here we present the discovery of the X-ray afterglow of a short-hard burst whose accurate position allows us to unambiguously associate it with a star-forming galaxy at redshift z=0.160, and whose optical lightcurve definitively excludes a supernova association. Together with results from three other recent short-hard bursts, this suggests that short-hard bursts release much less energy than the long-duration GRBs. Models requiring young stellar populations, such as magnetars and collapsars, are ruled out, while coalescing degenerate binaries remain the most promising progenitor candidates.

We present the results from an Hubble Space Telescope/ACS search for supernovae associated with X-ray flashes 020903, 040701, 040812 and 040916. We find strong evidence that XRF 020903 (z=0.25) was associated with a SN 1998bw-like supernova and confirm this using optical spectroscopy at t ~ 25 days. We find no evidence, however, for SN 1998bw-like supernovae associated with the other three events. In the case of XRF 040701 (z=0.21), we rule out even a faint supernova similar to SN 2002ap, using template light-curves for several local Type Ic supernovae. For the two cases in which the redshift is not known, XRFs 040812 and 040916, we derive robust redshift limits assuming they were accompanied by supernovae similar to SN 1998bw and compare these limits with photometric redshift constraints provided by their host galaxies. We supplement this analysis with results for three additional events (XRFs 011030, 020427 and 030723) and discuss the observed diversity of supernovae associated with X-ray flashes and gamma-ray bursts. We conclude that XRF-SNe exist, but can be significantly fainter than SN 1998bw, possibly consistent with the observed spread in local Type Ibc supernovae.

We present optical follow up observations of the long GRB 001007 between 6.14 hours and ~468 days after the event. An unusually bright optical afterglow (OA) was seen to decline following a steep power law decay with index alpha = -2.03 +/- 0.11, possibly indicating a break in the light curve at t - to < 3.5 days, as found in other bursts. Upper limits imposed by the LOTIS alerting system 6.14 hours after the gamma ray event provide tentative (1.2 sigma) evidence for a break in the optical light curve. The spectral index beta of the OA yields -1.24 +/- 0.57. These values may be explained both by several fireball jet models and by the cannonball model. Fireball spherical expansion models are not favoured. Late epoch deep imaging revealed the presence of a complex host galaxy system, composed of at least two objects located 1.2" (1.7 sigma) and 1.9" (2.7 sigma) from the afterglow position.

We present UBVI surface photometry for 20.'5 X 20.'5 area of a late-type spiral galaxy NGC 300. In order to understand the morphological properties and luminosity distribution characteristics of NGC 300, we have derived isophotal maps, surface brightness profiles, ellipticity profiles, position angle profiles, and color profiles. By merging the I-band data of our surface brightness measurements with those of Boeker et al. (2002) based on Hubble Space Telescope observations, we have made combined I-band surface brightness profiles for the region of 0."02 < r < 500" and decomposed the profiles into three components: a nucleus, a bulge, and an exponential disk.

We report 31 polarimetric observations of the afterglow of GRB 030329 with high signal-to-noise and high sampling frequency. The data imply that the afterglow magnetic field has small coherence length and is mostly random, probably generated by turbulence.

We report 31 polarimetric observations of the afterglow of GRB 030329 with high signal-to-noise and high sampling frequency. We establish the polarization light curve, detect sustained polarization at the percent level, and find significant variability of polarization degree and angle. The data imply that the afterglow magnetic field has small coherence length and is mostly random, probably generated by turbulence.

We present deep wide field $VI$ CCD photometry of the Sextans dwarf spheroidal galaxy (dSph) in the Local Group, covering a field of 42' x 28' located at the center of the galaxy. Color-magnitude diagrams of the Sextans dSph show well-defined red giant branch (RGB), blue horizontal branch (BHB), prominent red horizontal branch (RHB), asymptotic giant branch (AGB), about 120 variable star candidates including RR Lyraes and anomalous Cepheids, about 230 blue stragglers (BSs), and main sequence (MS) stars. The distance to the galaxy is derived using the $I$-band magnitude of the tip of the RGB at I(TRGB)=15.95+/-0.04: (m-M)_0=19.90+/-0.06 for an adopted reddening of E(B-V)=0.01. The mean metallicity of the RGB is estimated from the (V-I) color: [Fe/H]=-2.1+/-0.1(statistical error) +/-0.2(standard calibration error) dex, with a dispersion of sigma[Fe/H]=0.2 dex. The age of the MSTO of the main old population is estimated to be similar to that of M92, and there are seen some stellar populations with younger age. There is found to be one RGB bump at V=19.95+/-0.05 mag (M_V=0.03 mag), and a weak brighter bump at V=19.35+/-0.05 mag(Mv=-0.58 mag) which is probably an AGB bump. The V-band luminosity function of the RGB and MS stars is in general similar to that of M92, with a slight excess of stars in the magnitude range brighter than the MSTO with respect to that of M92. The bright BSs are more centrally concentrated than the faint BSs. Significant radial gradients are seen for several kinds of populations: the RHB, the red RGB, the red subgiant branch (SGB), and the bright BSs are more centrally concentrated toward the center of the galaxy, compared with the BHB, the blue RGB, the blue SGB, and the faint BSs, respectively.

We have observed a sample of Delta Scuti stars discovered by the ROTSE collaboration in 1999 with Super-LOTIS in order to characterize changes in their pulsation periods over a time baseline of roughly three years. Achieving these goals required the creation of an automated astrometric and photometric data reduction pipeline for the Super-LOTIS camera. Applying this pipeline to data from a June 2002 observing campaign, we detect pulsations in 18 objects, and find that in two cases the periods have changed significantly over the three years between the ROTSE and Super-LOTIS observations. Since theory predicts that evolutionary period changes should be quite small, sources of non-evolutionary period changes due to the interactions of pulsations modes are discussed.

We present Near-IR photometry of the Arches cluster, a young and massive stellar cluster near the Galactic center. We have analyzed the high resolution (FWHM 0.2") H and K' band images in the Galactic Center Demonstration Science Data Set, which were obtained with the Gemini/Hokupa's adaptive optics (AO) system. We present the color-magnitude diagram, the luminosity function and the initial mass function (IMF) of the stars in the Arches cluster in comparison with the HST/NICMOS data. The IMF slope for the range of 1.0< log(M/M_sun) <2.1 is estimated to be Gamma = -0.79+/-0.16, in good agreements with the earlier result based on the HST/NICMOS data [Figer et al. 1999, ApJ, 525, 750]. These results strengthen the evidence that the IMF of the bright stars close to the Galactic center is much flatter than that for the solar neighborhood. This is also consistent with a recent finding that the IMFs of the bright stars in young clusters in M33 get flatter as the galactocentric distance decreases [Lee et al. 2001, astro-ph 0109258]. It is found that the power of the Gemini/AO system is comparable, with some limits, to that of the HST/NICMOS.

We present multi-instrument optical observations of the High Energy Transient Explorer (HETE-2)/Interplanetary Network (IPN) error box of GRB 010921. This event was the first gamma ray burst (GRB) localized by HETE-2 which has resulted in the detection of an optical afterglow. In this paper we report the earliest known observations of the GRB010921 field, taken with the 0.11-m Livermore Optical Transient Imaging System (LOTIS) telescope, and the earliest known detection of the GRB010921 optical afterglow, using the 0.5-m Sloan Digital Sky Survey Photometric Telescope (SDSS PT). Observations with the LOTIS telescope began during a routine sky patrol 52 minutes after the burst. Observations were made with the SDSS PT, the 0.6-m Super-LOTIS telescope, and the 1.34-m Tautenburg Schmidt telescope at 21.3, 21.8, and 37.5 hours after the GRB, respectively. In addition, the host galaxy was observed with the USNOFS 1.0-m telescope 56 days after the burst. We find that at later times (t > 1 day after the burst), the optical afterglow exhibited a power-law decline with a slope of $\alpha = 1.75 \pm 0.28$. However, our earliest observations show that this power-law decline can not have extended to early times (t < 0.035 day).

We present interplanetary network localization, spectral, and time history information for four short-duration, hard spectrum gamma-ray bursts, GRB000607, 001025B, 001204, and 010119. All of these events were followed up with sensitive radio and optical observations (the first and only such bursts to be followed up in the radio to date), but no detections were made, demonstrating that the short bursts do not have anomalously intense afterglows. We discuss the upper limits, and show that the lack of observable counterparts is consistent both with the hypothesis that the afterglow behavior of the short bursts is like that of the long duration bursts, many of which similarly have no detectable afterglows, as well as with the hypothesis that the short bursts have no detectable afterglows at all. Small number statistics do not allow a clear choice between these alternatives, but given the present detection rates of various missions, we show that progress can be expected in the near future.

We present a photometric study of the stars in ionizing star clusters embedded in several giant HII regions of M33 (CC93, IC 142, NGC 595, MA2, NGC 604 and NGC 588). Our photometry are based on the HST-WFPC2 images of these clusters. Color-magnitude diagrams and color-color diagrams of these clusters are obtained and are used for estimating the reddenings and ages of the clusters. The luminosity functions (LFs) and initial mass functions (IMFs) of the massive stars in these clusters are also derived. The slopes of the IMFs range from Gamma = -0.5 to -2.1. It is found interestingly that the IMFs get steeper with increasing galactocentric distance and with decreasing [O/H] abundance.

GRB 990123 established the existence of prompt optical emission from gamma-ray bursts (GRBs). The Livermore Optical Transient Imaging System (LOTIS) has been conducting a fully automated search for this kind of simultaneous low energy emission from GRBs since October 1996. Although LOTIS has obtained simultaneous, or near simultaneous, coverage of the error boxes obtained with BATSE, IPN, XTE, and BeppoSAX for several GRBs, image analysis resulted in only upper limits. The unique gamma-ray properties of GRB 990123, such as very large fluence (top 0.4%) and hard spectrum, complicate comparisons with more typical bursts. We scale and compare gamma-ray properties, and in some cases afterglow properties, from the best LOTIS events to those of GRB 990123 in an attempt to determine whether the prompt optical emission of this event is representative of all GRBs. Furthermore, using LOTIS upper limits in conjunction with the relativistic blast wave model, we weakly constrain the GRB and afterglow parameters such as density of the circumburster medium and bulk Lorentz factor of the ejecta.

The optical transient of the faint Gamma Ray Burst 990308 was detected by the QUEST camera on the Venezuelan 1-m Schmidt telescope starting 3.28 hours after the burst. Our photometry gives $V = 18.32 \pm 0.07$, $R = 18.14 \pm 0.06$, $B = 18.65 \pm 0.23$, and $R = 18.22 \pm 0.05$ for times ranging from 3.28 to 3.47 hours after the burst. The colors correspond to a spectral slope of close to $f_{\nu} \propto \nu^{1/3}$. Within the standard synchrotron fireball model, this requires that the external medium be less dense than $10^{4} cm^{-3}$, the electrons contain $> 20%$ of the shock energy, and the magnetic field energy must be less than 24% of the energy in the electrons for normal interstellar or circumstellar densities. We also report upper limits of $V > 12.0$ at 132 s (with LOTIS), $V > 13.4$ from 132-1029s (with LOTIS), $V > 15.3$ at 28.2 min (with Super-LOTIS), and a 8.5 GHz flux of $< 114 \mu Jy$ at 110 days (with the Very Large Array). WIYN 3.5-m and Keck 10-m telescopes reveal this location to be empty of any host galaxy to $R > 25.7$ and $K > 23.3$. The lack of a host galaxy likely implies that it is either substantially subluminous or more distant than a red shift of $\sim 1.2$.

We report on the very early time search for an optical afterglow from GRB 971227 with the Livermore Optical Transient Imaging System (LOTIS). LOTIS began imaging the `Original' BATSE error box of GRB 971227 approximately 14 s after the onset of gamma-ray emission. Continuous monitoring of the position throughout the evening yielded a total of 499 images (10 s integration). Analysis of these images revealed no steady optical afterglow brighter than R=12.3 +- 0.2 in any single image. Coaddition of different combinations of the LOTIS images also failed to uncover transient optical emission. In particular, assuming a constant early time flux, no optical afterglow brighter than R=14.2 +- 0.2 was present within the first 1200 s and no optical afterglow brighter than R=15.0 +- 0.2 was present in the first 6.0 h. Follow up observations by other groups revealed a likely X-ray afterglow and a possible optical afterglow. Although subsequent deeper observations could not confirm a fading source, we show that these transients are not inconsistent with our present knowledge of the characteristics of GRB afterglows. We also demonstrate that with the upgraded thermoelectrically cooled CCDs, LOTIS is capable of either detecting very early time optical afterglow or placing stringent constraints on the relationship between the gamma-ray emission and the longer wavelength afterglow in relativistic blast wave models.

We present UBVI CCD photometry of two open clusters NGC 1798 and NGC 2192 which were little studied before. Color-magnitude diagrams of these clusters show several features typical for old open clusters: a well-defined main-sequence, a red giant clump, and a small number of red giants. The main sequence of NGC 1798 shows a distinct gap at V ~ 16.2 mag. From the surface number density distribution we have measured the size of the clusters, obtaining 8'.3 (= 10.2 pc) for NGC 1798 and 7'.3 (= 7.5 pc) for NGC 2192. Then we have determined the reddening, metallicity, and distance of these clusters using the color-color diagrams and color-magnitude diagrams: E(B-V)=0. 51\pm0.04, [Fe/H] = -0.47\pm0.15 dex and (m-M)_0=13.1\pm0.2 (d=4.2\pm 0.3 kpc) for NGC 1798, and E(B-V)=0.20\pm0.03, [Fe/H] = -0.31\pm0.15 dex and (m-M)_0=12.7\pm0.2 (d=3.5\pm 0.3 kpc) for NGC 2192. The ages of these clusters have been estimated using the morphological age indicators and the isochrone fitting with the Padova isochrones: 1.4 \pm 0.2 Gyrs for NGC 1798 and 1.1 \pm 0.1 Gyrs for NGC 2192. The luminosity functions of the main sequence stars in these clusters are found to be similar to other old open clusters. The metallicity and distance of these clusters are consistent with the relation between the metallicity and galactocentric distance of other old open clusters.

We are constructing a 0.6 meter telescope system to search for early time gamma-ray burst(GRB) optical counterparts. Super-LOTIS (Super-Livermore Optical Transient Imaging System) is an automated telescope system that has a 0.8 x 0.8 deg field-of-view, is sensitive to Mv ~ 19 and responds to a burst trigger within 5 min. This telescope will record images of the gamma-ray burst coordinates that is given by the GCN (GRB Coordinate Network). A measurement of GRB light curves at early times will greatly enhance our understanding of GRB physics.

LOTIS (Livermore Optical Transient Imaging System) is a gamma-ray burst optical counterpart search experiment located near Lawrence Livermore National Laboratory in California. The system is linked to the GCN (GRB Coordinates Network) real-time coordinate distribution network and can respond to a burst trigger in 6-15 seconds. LOTIS has a total field-of-view of 17.4 deg x 17.4 deg with a completeness sensitivity of m_V ~ 11 for a 10 second integration time. Since operations began in October 1996, LOTIS has responded to over 30 GCN/BATSE GRB triggers. Seven of these triggers are considered good events subject to the criteria of clear weather conditions, < 60 s response time, and > 50% coverage of the final BATSE 3sigma error circle. We discuss results from the first year of LOTIS operations with an emphasis on the observations and analysis of GRB 971006 (BATSE trigger 6414).

LOTIS is a gamma-ray burst optical counterpart search experiment located near Lawrence Livermore National Laboratory in California. Since operations began in October 1996, LOTIS has responded to five triggers as of July 30, 1997, which occurred during good weather conditions. GRB970223 (BATSE Trigger #6100) was an exceptionally strong burst lasting $\sim30$ s with a peak at $\sim8$ s. LOTIS began imaging the error box $\sim 11$ s after the burst began, and achieved simultaneous optical coverage of 100% of the region enclosed by the BATSE $3\sigma$ error circle and the IPN annulus. No optical transients were observed brighter than the m$_V \sim 11$ completeness limit of the resulting images providing a new upper limit on the simultaneous optical to gamma-ray fluence ratio of $R_L < 1.1 \times 10^{-4}$ and on the simultaneous optical (at 700 nm) to gamma-ray (at 100 keV) flux density ratio of $R_F < 305$ for a B type spectrum and $R_F < 475$ for an M type spectrum.

The Gamma-Ray Optical Counterpart Search Experiment (GROCSE) presents new experimental upper limits on the optical flux from gamma-ray bursts (GRBs). Our experiment consisted of a fully-automated very wide-field opto-electronic detection system that imaged locations of GRBs within a few seconds of receiving trigger signals provided by BATSE's real-time burst coordinate distribution network (BACODINE). The experiment acquired ~3800 observing hours, recording 22 gamma-ray burst triggers within $\sim$30 s of the start of the burst event. Some of these bursts were imaged while gamma-ray radiation was being detected by BATSE. We identified no optical counterparts associated with gamma-ray bursts amongst these events at the m$_V$ $\sim$ 7.0 to 8.5 sensitivity level. We find the ratio of the upper limit to the V-band optical flux, F$_\nu$, to the gamma-ray fluence, $\Phi_\gamma$, from these data to be $2 \times 10^{-18} < F_\nu/\Phi_\gamma < 2 \times 10^{-16}$.

We report the detection of 45 candidate microlensing events in fields toward the Galactic bulge. These come from the analysis of 24 fields containing 12.6 million stars observed for 190 days in 1993. Many of these events are of extremely high signal to noise and are remarkable examples of gravitational microlensing. The distribution of peak magnifications is shown to be consistent with the microlensing interpretation of these events. Using a sub-sample of 1.3 million ``Clump Giant" stars whose distance and detection efficiency are well known, we find 13 events and estimate the microlensing optical depth toward the Galactic Bulge as $\tau_{\rm bulge} = 3.9 {+ 1.8 \atop - 1.2} \times 10^{-6}$ averaged over an area of $\sim 12$ square degrees centered at Galactic coordinates $\ell = 2.55^\circ$ and $b = -3.64^\circ$. This is similar to the value reported by the OGLE collaboration, and is marginally higher than current theoretical models for $\tau_{\rm bulge}$. The optical depth is also seen to increase significantly for decreasing $\vert b\vert$. These results demonstrate that obtaining large numbers of microlensing events toward the Galactic bulge is feasible, and that the study of such events will have important consequences for the structure of the Galaxy and its dark halo.

The MACHO collaboration reports on the analysis of our first year LMC data, 9.5 million light curves with an average of 235 observations each. Automated selection procedures give 3 events consistent with microlensing. We evaluate our experimental detection efficiency using a range of Monte- Carlo simulations. Using a `standard' halo density profile we find that a halo comprised entirely of Machos in the mass range 3 \ten{-4} to 0.06 \msun would predict > 15 detected events in this dataset; thus a standard spherical halo cannot be dominated by objects in this mass range. Assuming all three events are microlensing of halo objects and fitting a naive spherical halo model to our data yields a Macho halo fraction f =0.19+0.16-0.10, a total mass in Machos (inside 50 kpc) of 7.6+6-4 \ten{10} \msun, and a microlensing optical depth 8.8+7-5 \ten{-8} (68\% CL). Exploring a wide range of halo models we find that our constraints on the Macho fraction are quite model-dependent, but constraints on the total mass in Machos within 50 kpc are quite secure.

We have monitored 8.6 million stars in the Large Magellanic Cloud for 1.1 years and have found 3 events consistent with gravitational microlensing. We place strong constraints on the Galactic halo content in the form of compact lensing objects in the mass range $10^{-4} \msun$ to $10^{-1} \msun$. Three events is fewer than expected for a standard spherical halo of objects in this mass range, but appears to exceed the number expected from known Galactic populations. Fitting a naive spherical halo model to our data yields a MACHO fraction $f = 0.20^{+0.33}_{-0.14}$, which implies a total MACHO mass (inside 50 kpc) of $8.0^{+14}_{-6}\ten{10} \msun$, and a microlensing optical depth $9^{+15}_{-7} \ten{-8}$ ($\sim 68$\% CL).

There is now abundant evidence for the presence of large quantities of unseen matter surrounding normal galaxies, including our own$^{1,2}$. The nature of this `dark matter' is unknown, except that it cannot be made of normal stars, dust, or gas, as they would be easily detected. Exotic particles such as axions, massive neutrinos or other weakly interacting massive particles (collectively known as WIMPs) have been proposed)$^{3,4}$, but have yet to be detected. A less exotic alternative is normal matter in the form of bodies with masses ranging from that of a large planet to a few $ \msun$. Such objects, known collectively as massive compact halo objects$^5$ (MACHOs) might be brown dwarfs or `Jupiters' (bodies too small to produce their own energy by fusion), neutron stars, old white dwarfs, or black holes. Paczynski$^6$ suggested that MACHOs might act as gravitational microlenses, occasionally causing the apparent brightness of distant background stars temporarily to increase. We are conducting a microlensing experiment to determine whether the dark matter halo of our galaxy is made up of MACHOs. Here we report a candidate for a microlensing event, detected by monitoring the light curves of 1.8 million stars in the Large Magellanic Cloud for one year. The light curve shows no variation for most of the year of data taking, and an upward excursion lasting over 1 month, with a maximum increase of $\approx \bf 2$ mag. The most probable lens mass, inferred from the duration of the event, is $\bf \sim 0.1 \,\msun$.