Showing 1–50 of 127 results for author: Choi, J

The daily modulations and broadband strategy in axion searches. An application with CAST-CAPP detector

Authors: C. M. Adair, K. Altenmüller, V. Anastassopoulos, S. Arguedas Cuendis, J. Baier, K. Barth, A. Belov, D. Bozicevic, H. Bräuninger, G. Cantatore, F. Caspers, J. F. Castel, S. A. Çetin, W. Chung, H. Choi, J. Choi, T. Dafni, M. Davenport, A. Dermenev, K. Desch, B. Döbrich, H. Fischer, W. Funk, J. Galan, A. Gardikiotis , et al. (38 additional authors not shown)

Abstract: It has been previously advocated that the presence of the daily and annual modulations of the axion flux on the Earth's surface may dramatically change the strategy of the axion searches. The arguments were based on the so-called Axion Quark Nugget (AQN) dark matter model which was originally put forward to explain the similarity of the dark and visible cosmological matter densities… ▽ More It has been previously advocated that the presence of the daily and annual modulations of the axion flux on the Earth's surface may dramatically change the strategy of the axion searches. The arguments were based on the so-called Axion Quark Nugget (AQN) dark matter model which was originally put forward to explain the similarity of the dark and visible cosmological matter densities $Ω_{\rm dark}\sim Ω_{\rm visible}$. In this framework, the population of galactic axions with mass $ 10^{-6} {\rm eV}\lesssim m_a\lesssim 10^{-3}{\rm eV}$ and velocity $\langle v_a\rangle\sim 10^{-3} c$ will be accompanied by axions with typical velocities $\langle v_a\rangle\sim 0.6 c$ emitted by AQNs. Furthermore, in this framework, it has also been argued that the AQN-induced axion daily modulation (in contrast with the conventional WIMP paradigm) could be as large as $(10-20)\%$, which represents the main motivation for the present investigation. We argue that the daily modulations along with the broadband detection strategy can be very useful tools for the discovery of such relativistic axions. The data from the CAST-CAPP detector have been used following such arguments. Unfortunately, due to the dependence of the amplifier chain on temperature-dependent gain drifts and other factors, we could not conclusively show the presence or absence of a dark sector-originated daily modulation. However, this proof of principle analysis procedure can serve as a reference for future studies. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 18 pages, 8 figures

Effects of pair freeze-out on photon distributions in BBN epoch

Authors: Jeongyoon Choi, Dukjae Jang, Youngshin Kwon, Gwangeon Seong, Myeong Hwan Mun, Young-Min Kim, Kyujin Kwak, Myung-Ki Cheoun

Abstract: We investigate the evolution of non-extensivity in the photon distribution during the Big Bang Nucleosynthesis (BBN) epoch using Tsallis statistics. Assuming a minimal deviation from the Planck distribution, we construct the perturbed Boltzmann equation for photons, including the collision terms for pair creation and annihilation processes. We analyze the possibility that these collisions could ca… ▽ More We investigate the evolution of non-extensivity in the photon distribution during the Big Bang Nucleosynthesis (BBN) epoch using Tsallis statistics. Assuming a minimal deviation from the Planck distribution, we construct the perturbed Boltzmann equation for photons, including the collision terms for pair creation and annihilation processes. We analyze the possibility that these collisions could cause a slight increase in the number of high-frequency photons within the BBN era, and consequently, the primordial plasma might be temporarily placed in a state of chemical non-equilibrium. We also discuss the restoration of the photon distribution to an equilibrium state as the Universe enters the matter-dominated era. These findings, which suggest possible changes in the photon distribution during the epoch between the BBN and the recombination, offer insights that support the previously proposed ansatz solution to the primordial lithium problem in arXiv:1812.09472. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 9 pages, 2 figures

Alpha backgrounds in NaI(Tl) crystals of COSINE-100

Authors: G. Adhikari, N. Carlin, D. F. F. S. Cavalcante, J. Y. Cho, J. J. Choi, S. Choi, A. C. Ezeribe, L. E. Franca, C. Ha, I. S. Hahn, S. J. Hollick, E. J. Jeon, H. W. Joo, W. G. Kang, M. Kauer, B. H. Kim, H. J. Kim, J. Kim, K. W. Kim, S. H. Kim, S. K. Kim, S. W. Kim, W. K. Kim, Y. D. Kim, Y. H. Kim , et al. (38 additional authors not shown)

Abstract: COSINE-100 is a dark matter direct detection experiment with 106 kg NaI(Tl) as the target material. 210Pb and daughter isotopes are a dominant background in the WIMP region of interest and are detected via beta decay and alpha decay. Analysis of the alpha channel complements the background model as observed in the beta/gamma channel. We present the measurement of the quenching factors and Monte Ca… ▽ More COSINE-100 is a dark matter direct detection experiment with 106 kg NaI(Tl) as the target material. 210Pb and daughter isotopes are a dominant background in the WIMP region of interest and are detected via beta decay and alpha decay. Analysis of the alpha channel complements the background model as observed in the beta/gamma channel. We present the measurement of the quenching factors and Monte Carlo simulation results and activity quantification of the alpha decay components of the COSINE-100 NaI(Tl) crystals. The data strongly indicate that the alpha decays probabilistically undergo two possible quenching factors but require further investigation. The fitted results are consistent with independent measurements and improve the overall understanding of the COSINE-100 backgrounds. Furthermore, the half-life of 216Po has been measured to be 143.4 +/- 1.2 ms, which is consistent with and more precise than recent measurements. △ Less

Submitted 30 January, 2024; v1 submitted 8 November, 2023; originally announced November 2023.

Performance of an ultra-pure NaI(Tl) detector produced by an indigenously-developed purification method and crystal growth for the COSINE-200 experiment

Authors: Hyun Seok Lee, Byung Ju Park, Jae Jin Choi, Olga Gileva, Chang Hyon Ha, Alain Iltis, Eun Ju Jeon, Dae Yeon Kim, Kyung Won Kim, Sung Hyun Kim, Sun Kee Kim, Yeong Duk Kim, Young Ju Ko, Cheol Ho Lee, Hyun Su Lee, In Soo Lee, Moo Hyun Lee, Se Jin Ra, Ju Kyung Son, Keon Ah Shin

Abstract: The COSINE-100 experiment has been operating with 106 kg of low-background NaI(Tl) detectors to test the results from the DAMA/LIBRA experiment, which claims to have observed dark matter. However, since the background of the NaI(Tl) crystals used in the COSINE-100 experiment is 2-3 times higher than that in the DAMA detectors, no conclusion regarding the claimed observation from the DAMA/LIBRA exp… ▽ More The COSINE-100 experiment has been operating with 106 kg of low-background NaI(Tl) detectors to test the results from the DAMA/LIBRA experiment, which claims to have observed dark matter. However, since the background of the NaI(Tl) crystals used in the COSINE-100 experiment is 2-3 times higher than that in the DAMA detectors, no conclusion regarding the claimed observation from the DAMA/LIBRA experiment could be reached. Therefore, we plan to upgrade the current COSINE-100 experiment to the next phase, COSINE-200, by using ultra-low background NaI(Tl) detectors. The basic principle was already proved with the commercially available Astro-grade NaI powder from Sigma-Aldrich company. However, we have developed a mass production process of ultra-pure NaI powder at the Center for Underground Physics (CUP) of the Institute for Basic Science (IBS), Korea, using the direct purification of the raw NaI powder. We plan to produce more than 1,000 kg of ultra-pure powder for the COSINE200 experiment. With our crystal grower installed at CUP, we have successfully grown a low-background crystal using our purification technique for the NaI powder. We have assembled a low-background NaI(Tl) detector. In this article, we report the performance of this ultra-pure NaI(Tl) crystal detector produced at IBS, Korea. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Journal ref: Frontiers in Physics 11 2023

An induced annual modulation signature in COSINE-100 data by DAMA/LIBRA's analysis method

Authors: G. Adhikari, N. Carlin, J. J. Choi, S. Choi, A. C. Ezeribe, L. E. Franca, C. Ha, I. S. Hahn, S. J. Hollick, E. J. Jeon, J. H. Jo, H. W. Joo, W. G. Kang, M. Kauer, B. H. Kim, H. J. Kim, J. Kim, K. W. Kim, S. H. Kim, S. K. Kim, W. K. Kim, Y. D. Kim, Y. H. Kim, Y. J. Ko, D. H. Lee , et al. (32 additional authors not shown)

Abstract: The DAMA/LIBRA collaboration has reported the observation of an annual modulation in the event rate that has been attributed to dark matter interactions over the last two decades. However, even though tremendous efforts to detect similar dark matter interactions were pursued, no definitive evidence has been observed to corroborate the DAMA/LIBRA signal. Many studies assuming various dark matter mo… ▽ More The DAMA/LIBRA collaboration has reported the observation of an annual modulation in the event rate that has been attributed to dark matter interactions over the last two decades. However, even though tremendous efforts to detect similar dark matter interactions were pursued, no definitive evidence has been observed to corroborate the DAMA/LIBRA signal. Many studies assuming various dark matter models have attempted to reconcile DAMA/LIBRA's modulation signals and null results from other experiments, however no clear conclusion can be drawn. Apart from the dark matter hypothesis, several studies have examined the possibility that the modulation is induced by variations in their detector's environment or their specific analysis methods. In particular, a recent study presents a possible cause of the annual modulation from an analysis method adopted by the DAMA/LIBRA experiment in which the observed annual modulation could be reproduced by a slowly varying time-dependent background. Here, we study the COSINE-100 data using an analysis method similar to the one adopted by the DAMA/LIBRA experiment and observe a significant annual modulation, although the modulation phase is almost opposite to that of the DAMA/LIBRA data. Assuming the same background composition for COSINE-100 and DAMA/LIBRA, simulated experiments for the DAMA/LIBRA without dark matter signals also provide significant annual modulation with an amplitude similar to DAMA/LIBRA with opposite phase. Even though this observation does not explain the DAMA/LIBRA's results directly, this interesting phenomenon motivates deeper studies of the time-dependent DAMA/LIBRA background data. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Journal ref: Sci. Rep. 13, 4676 (2023)

Three-year annual modulation search with COSINE-100

Authors: COSINE-100 Collaboration, :, G. Adhikari, E. Barbosa de Souza, N. Carlin, J. J. Choi, S. Choi, A. C. Ezeribe, L. E. França, C. Ha, I. S. Hahn, S. J. Hollick, E. J. Jeon, J. H. Jo, H. W. Joo, W. G. Kang, M. Kauer, H. Kim, H. J. Kim, J. Kim, K. W. Kim, S. H. Kim, S. K. Kim, W. K. Kim, Y. D. Kim , et al. (34 additional authors not shown)

Abstract: COSINE-100 is a direct detection dark matter experiment that aims to test DAMA/LIBRA's claim of dark matter discovery by searching for a dark matter-induced annual modulation signal with NaI(Tl) detectors. We present new constraints on the annual modulation signal from a dataset with a 2.82 yr livetime utilizing an active mass of 61.3 kg, for a total exposure of 173 kg$\cdot$yr. This new result fe… ▽ More COSINE-100 is a direct detection dark matter experiment that aims to test DAMA/LIBRA's claim of dark matter discovery by searching for a dark matter-induced annual modulation signal with NaI(Tl) detectors. We present new constraints on the annual modulation signal from a dataset with a 2.82 yr livetime utilizing an active mass of 61.3 kg, for a total exposure of 173 kg$\cdot$yr. This new result features an improved event selection that allows for both lowering the energy threshold to 1 keV and a more precise time-dependent background model. In the 1-6 keV and 2-6 keV energy intervals, we observe best-fit values for the modulation amplitude of 0.0067$\pm$0.0042 and 0.0051$\pm$0.0047 counts/(day$\cdot$kg$\cdot$keV), respectively, with a phase fixed at 152.5 days. △ Less

Submitted 28 October, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 11 pages, 8 figures

Journal ref: Phys. Rev. D 106, 052005 (2022)

Rotation acceleration of asteroids (10115) 1992 SK, (1685) Toro, and (1620) Geographos due to the YORP effect

Authors: J. Durech, D. Vokrouhlicky, P. Pravec, Yu. N. Krugly, M. J. Kim, D. Polishook, V. V. Ayvazian, T. Bonev, Y. J. Choi, D. G. Datashvili, Z. Donchev, S. A. Ehgamberdiev, K. Hornoch, R. Ya. Inasaridze, G. V. Kapanadze, D. H. Kim, H. Kucakova, A. V. Kusakin, P. Kusnirak, H. J. Lee, I. E. Molotov, H. K. Moon, S. S. Mykhailova, I. V. Nikolenko, A. Novichonok , et al. (6 additional authors not shown)

Abstract: The rotation state of small asteroids is affected by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, which is a net torque caused by solar radiation directly reflected and thermally reemitted from the surface. Due to this effect, the rotation period slowly changes, which can be most easily measured in light curves because the shift in the rotation phase accumulates over time quadratically… ▽ More The rotation state of small asteroids is affected by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, which is a net torque caused by solar radiation directly reflected and thermally reemitted from the surface. Due to this effect, the rotation period slowly changes, which can be most easily measured in light curves because the shift in the rotation phase accumulates over time quadratically. We collected archived light curves and carried out new photometric observations for asteroids (10115) 1992 SK, (1620) Geographos, and (1685) Toro. We applied the method of light curve inversion to fit observations with a convex shape model. The YORP effect was modeled as a linear change of the rotation frequency $\upsilon \equiv \mathrm{d}ω/ \mathrm{d}t$ and optimized together with other spin and shape parameters. We detected the acceleration $\upsilon = (8.3 \pm 0.6) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ of the rotation for asteroid (10115) 1992 SK. This observed value agrees well with the theoretical value of YORP-induced spin-up computed for our shape and spin model. For (1685) Toro, we obtained $\upsilon = (3.3 \pm 0.3) \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$, which confirms an earlier tentative YORP detection. For (1620) Geographos, we confirmed the previously detected YORP acceleration and derived an updated value of $\upsilon$ with a smaller uncertainty. We also included the effect of solar precession into our inversion algorithm, and we show that there are hints of this effect in Geographos' data. The detected change of the spin rate of (10115) 1992 SK has increased the total number of asteroids with YORP detection to ten. In all ten cases, the $\mathrm{d}ω/ \mathrm{d}t$ value is positive, so the rotation of these asteroids is accelerated. It is unlikely to be just a statistical fluke, but it is probably a real feature that needs to be explained. △ Less

Submitted 13 October, 2021; originally announced October 2021.

Journal ref: A&A 657, A5 (2022)

Searching for low-mass dark matter via Migdal effect in COSINE-100

Authors: G. Adhikari, N. Carlin, J. J. Choi, S. Choi, A. C. Ezeribe, L. E. Franca, C. Ha, I. S. Hahn, S. J. Hollick, E. J. Jeon, J. H. Jo, H. W. Joo, W. G. Kang, M. Kauer, H. Kim, H. J. Kim, J. Kim, K. W. Kim, S. H. Kim, S. K. Kim, W. K. Kim, Y. D. Kim, Y. H. Kim, Y. J. Ko, H. J. Kwon , et al. (31 additional authors not shown)

Abstract: We report on the search for weakly interacting massive particle (WIMP) dark matter candidates in the galactic halo that interact with sodium and iodine nuclei in the COSINE-100 experiment and produce energetic electrons that accompany recoil nuclei via the the Migdal effect. The WIMP mass sensitivity of previous COSINE-100 searches that relied on the detection of ionization signals produced by tar… ▽ More We report on the search for weakly interacting massive particle (WIMP) dark matter candidates in the galactic halo that interact with sodium and iodine nuclei in the COSINE-100 experiment and produce energetic electrons that accompany recoil nuclei via the the Migdal effect. The WIMP mass sensitivity of previous COSINE-100 searches that relied on the detection of ionization signals produced by target nuclei recoiling from elastic WIMP-nucleus scattering was restricted to WIMP masses above $\sim$5 GeV/$c^2$ by the detectors' 1 keVee energy-electron-equivalent threshold. The search reported here looks for recoil signals enhanced by the Migdal electrons that are ejected during the scattering process. This is particularly effective for the detection of low-mass WIMP scattering from the crystals' sodium nuclei in which a relatively larger fraction of the WIMP's energy is transferred to the nucleus recoil energy and the excitation of its orbital electrons. In this analysis, the low-mass WIMP search window of the COSINE-100 experiment is extended to WIMP mass down to 200 MeV/$c^2$. The low-mass WIMP sensitivity will be further improved by lowering the analysis threshold based on a multivariable analysis technique. We consider the influence of these improvements and recent developments in detector performance to re-evaluate sensitivities for the future COSINE-200 experiment. With a 0.2 keVee analysis threshold and high light-yield NaI(Tl) detectors (22 photoelectrons/keVee), the COSINE-200 experiment can explore low-mass WIMPs down to 20 MeV/$c^2$ and probe previously unexplored regions of parameter space. △ Less

Submitted 10 January, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 10 pages, 8 figures

Journal ref: Phys. Rev. D 105, 042006 (2022)

Strong constraints from COSINE-100 on the DAMA dark matter results using the same sodium iodide target

Authors: G. Adhikari, E. Barbosa de Souza, N. Carlin, J. J. Choi, S. Choi, M. Djamal, A. C. Ezeribe, L. E. França, C. Ha, I. S. Hahn, E. J. Jeon, J. H. Jo, H. W. Joo, W. G. Kang, M. Kauer, H. Kim, H. J. Kim, K. W. Kim, S. H. Kim, S. K. Kim, W. K. Kim, Y. D. Kim, Y. H. Kim, Y. J. Ko, E. K. Lee , et al. (28 additional authors not shown)

Abstract: We present new constraints on dark matter interactions using 1.7 years of COSINE-100 data. The COSINE-100 experiment, consisting of 106 kg of tallium-doped sodium iodide (NaI(Tl)) target material, is aimed at testing DAMA's claim of dark matter observation using the same NaI(Tl) detectors. Improved event selection requirements, a more precise understanding of the detector background and the use of… ▽ More We present new constraints on dark matter interactions using 1.7 years of COSINE-100 data. The COSINE-100 experiment, consisting of 106 kg of tallium-doped sodium iodide (NaI(Tl)) target material, is aimed at testing DAMA's claim of dark matter observation using the same NaI(Tl) detectors. Improved event selection requirements, a more precise understanding of the detector background and the use of a larger data set considerably enhances the COSINE-100 sensitivity for dark matter detection. No signal consistent with the dark matter interaction is identified, and rules out model-dependent dark matter interpretations of the DAMA signals in the specific context of standard halo model with the same NaI(Tl) target for various interaction hypotheses. △ Less

Submitted 26 August, 2021; v1 submitted 8 April, 2021; originally announced April 2021.

Journal ref: Sci. Adv. 7, eabk2699 (2021)

Background modeling for dark matter search with 1.7 years of COSINE-100 data

Authors: G. Adhikari, P. Adhikari, E. Barbosa de Souza, N. Carlin, J. J. Choi, S. Choi, M. Djamal, A. C. Ezeribe, L. E. Franca, C. Ha, I. S. Hahn, E. J. Jeon, J. H. Jo, W. G. Kang, M. Kauer, G. S. Kim, H. Kim, H. J. Kim, K. W. Kim, N. Y. Kim, S. K. Kim, Y. D. Kim, Y. H. Kim, Y. J. Ko, V. A. Kudryavtsev , et al. (25 additional authors not shown)

Abstract: We present a background model for dark matter searches using an array of NaI(Tl) crystals in the COSINE-100 experiment that is located in the Yangyang underground laboratory. The model includes background contributions from both internal and external sources, including cosmogenic radionuclides and surface $^{210}$Pb contamination. To build the model in the low energy region, with a threshold of 1… ▽ More We present a background model for dark matter searches using an array of NaI(Tl) crystals in the COSINE-100 experiment that is located in the Yangyang underground laboratory. The model includes background contributions from both internal and external sources, including cosmogenic radionuclides and surface $^{210}$Pb contamination. To build the model in the low energy region, with a threshold of 1 keV, we used a depth profile of $^{210}$Pb contamination in the surface of the NaI(Tl) crystals determined in a comparison between measured and simulated spectra. We also considered the effect of the energy scale errors propagated from the statistical uncertainties and the nonlinear detector response at low energies. The 1.7 years COSINE-100 data taken between October 21, 2016 and July 18, 2018 were used for this analysis. Our Monte Carlo simulation provides a non-Gaussian peak around 50 keV originating from beta decays of bulk $^{210}$Pb in a good agreement with the measured background. This model estimates that the activities of bulk $^{210}$Pb and $^{3}$H are dominating the background rate that amounts to an average level of 2.85$\pm$0.15 counts/day/keV/kg in the energy region of (1-6) keV, using COSINE-100 data with a total exposure of 97.7 kg$\cdot$years. △ Less

Submitted 6 September, 2022; v1 submitted 27 January, 2021; originally announced January 2021.

Supernova Model Discrimination with Hyper-Kamiokande

Authors: Hyper-Kamiokande Collaboration, :, K. Abe, P. Adrich, H. Aihara, R. Akutsu, I. Alekseev, A. Ali, F. Ameli, I. Anghel, L. H. V. Anthony, M. Antonova, A. Araya, Y. Asaoka, Y. Ashida, V. Aushev, F. Ballester, I. Bandac, M. Barbi, G. J. Barker, G. Barr, M. Batkiewicz-Kwasniak, M. Bellato, V. Berardi, M. Bergevin , et al. (478 additional authors not shown)

Abstract: Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-colla… ▽ More Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-collapse supernovae is not yet well understood. Hyper-Kamiokande is a next-generation neutrino detector that will be able to observe the neutrino flux from the next galactic core-collapse supernova in unprecedented detail. We focus on the first 500 ms of the neutrino burst, corresponding to the accretion phase, and use a newly-developed, high-precision supernova event generator to simulate Hyper-Kamiokande's response to five different supernova models. We show that Hyper-Kamiokande will be able to distinguish between these models with high accuracy for a supernova at a distance of up to 100 kpc. Once the next galactic supernova happens, this ability will be a powerful tool for guiding simulations towards a precise reproduction of the explosion mechanism observed in nature. △ Less

Submitted 20 July, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 21 pages, 7 figures. Article based on thesis published as arXiv:2002.01649. v2: added references and some explanations in response to reviewer comments

Journal ref: Astrophys.J. 916 (2021) 15

GroundBIRD : A CMB polarization experiment with MKID arrays

Authors: Kyungmin Lee, Jihoon Choi, Ricardo Tanausú Génova-Santos, Makoto Hattori, Masashi Hazumi, Shunsuke Honda, Takuji Ikemitsu, Hidesato Ishida, Hikaru Ishitsuka, Yonggil Jo, Kenichi Karatsu, Kenji Kiuchi, Junta Komine, Ryo Koyano, Hiroki Kutsuma, Satoru Mima, Makoto Minowa, Joonhyeok Moon, Makoto Nagai, Taketo Nagasaki, Masato Naruse, Shugo Oguri, Chiko Otani, Michael Peel, Rafael Rebolo , et al. (9 additional authors not shown)

Abstract: GroundBIRD is a ground-based experiment for the precise observation of the polarization of the cosmic microwave background (CMB). To achieve high sensitivity at large angular scale, we adopt three features in this experiment: fast rotation scanning, microwave kinetic inductance detector (MKID) and cold optics. The rotation scanning strategy has the advantage to suppress $1/f$ noise. It also provid… ▽ More GroundBIRD is a ground-based experiment for the precise observation of the polarization of the cosmic microwave background (CMB). To achieve high sensitivity at large angular scale, we adopt three features in this experiment: fast rotation scanning, microwave kinetic inductance detector (MKID) and cold optics. The rotation scanning strategy has the advantage to suppress $1/f$ noise. It also provides a large sky coverage of 40\%, which corresponds to the large angular scales of $l \sim 6$. This allows us to constrain the tensor-to-scalar ratio by using low $l$ B-mode spectrum. The focal plane consists of 7 MKID arrays for two target frequencies, 145 GHz and 220 GHz band. There are 161 pixels in total, of which 138 are for 144 GHz and 23 are for 220 GHz. This array is currently under development and the prototype will soon be evaluated in telescope. The GroundBIRD telescope will observe the CMB at the Teide observatory. The telescope was moved from Japan to Tenerife and is now under test. We present the status and plan of the GroundBIRD experiment. △ Less

Submitted 15 November, 2020; originally announced November 2020.

Comments: 7 pages, 5 figures, LTD18 proceeding, Published in JLTP

Journal ref: J Low Temp Phys 200, 384 (2020)

Precision cosmology with voids in the final BOSS data

Authors: Nico Hamaus, Alice Pisani, Jin-Ah Choi, Guilhem Lavaux, Benjamin D. Wandelt, Jochen Weller

Abstract: We report novel cosmological constraints obtained from cosmic voids in the final BOSS DR12 dataset. They arise from the joint analysis of geometric and dynamic distortions of average void shapes (i.e., the stacked void-galaxy cross-correlation function) in redshift space. Our model uses tomographic deprojection to infer real-space void profiles and self-consistently accounts for the Alcock-Paczyns… ▽ More We report novel cosmological constraints obtained from cosmic voids in the final BOSS DR12 dataset. They arise from the joint analysis of geometric and dynamic distortions of average void shapes (i.e., the stacked void-galaxy cross-correlation function) in redshift space. Our model uses tomographic deprojection to infer real-space void profiles and self-consistently accounts for the Alcock-Paczynski (AP) effect and redshift-space distortions (RSD) without any prior assumptions on cosmology or structure formation. It is derived from first physical principles and provides an extremely good description of the data at linear perturbation order. We validate this model with the help of mock catalogs and apply it to the final BOSS data to constrain the RSD and AP parameters $f/b$ and $D_AH/c$, where $f$ is the linear growth rate, $b$ the linear galaxy bias, $D_A$ the comoving angular diameter distance, $H$ the Hubble rate, and $c$ the speed of light. In addition, we include two nuisance parameters in our analysis to marginalize over potential systematics. We obtain $f/b=0.540\pm0.091$ and $D_AH/c=0.588\pm0.004$ from the full void sample at a mean redshift of $z=0.51$. In a flat $Λ$CDM cosmology, this implies $Ω_\mathrm{m}=0.312\pm0.020$ for the present-day matter density parameter. When we use additional information from the survey mocks to calibrate our model, these constraints improve to $f/b=0.347\pm0.023$, $D_AH/c=0.588\pm0.003$, and $Ω_\mathrm{m}=0.310\pm0.017$. However, we emphasize that the calibration depends on the specific model of cosmology and structure formation assumed in the mocks, so the calibrated results should be considered less robust. Nevertheless, our calibration-independent constraints are among the tightest of their kind to date, demonstrating the immense potential of using cosmic voids for cosmology in current and future data. △ Less

Submitted 19 November, 2020; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: 46 pages, 18 figures, 1 table. JCAP accepted version after peer review, typos corrected

Journal ref: Journal of Cosmology and Astroparticle Physics, Issue 12, article id. 023 (2020)

Improving the light collection using a new NaI(Tl)crystal encapsulation

Authors: J. J. Choi, B. J. Park, C. Ha, K. W. Kim, S. K. Kim, Y. D. Kim, Y. J. Ko, H. S. Lee, S. H. Lee, S. L. Olsen

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

Submitted 12 August, 2020; v1 submitted 3 June, 2020; originally announced June 2020.

Comments: 15 pages, 7 figures (minor edits)

Improved axion haloscope search analysis

Authors: S. Ahn, S. Lee, J. Choi, B. R. Ko, Y. K. Semertzidis

Abstract: One of the most significant and practical figures of merit in axion haloscope searches is the scanning rate, because of the unknown axion mass. Under the best experimental parameters, the only way to improve the figure of merit is to increase the experimentally designed signal to noise ratio in the axion haloscope search analysis procedure. In this paper, we report an improved axion haloscope sear… ▽ More One of the most significant and practical figures of merit in axion haloscope searches is the scanning rate, because of the unknown axion mass. Under the best experimental parameters, the only way to improve the figure of merit is to increase the experimentally designed signal to noise ratio in the axion haloscope search analysis procedure. In this paper, we report an improved axion haloscope search analysis using the data taken by the CAPP-8TB haloscope. By correcting for the background biased by the background parametrizations in the presence of axion signals, we realized a signal to noise ratio efficiency of about 100\%. Given the axion haloscope search analyses to date, the scanning rate can be improved by 21\%, with about a 10\% improvement in the signal to noise ratio. This improvement is another low cost innovation in axion haloscope searches, where all the experimental parameters are currently at their best. △ Less

Submitted 7 April, 2021; v1 submitted 16 April, 2020; originally announced April 2020.

Comments: 15 pages, 11 figures, JHEP accepted

Development of ultra-pure NaI(Tl) detectors for the COSINE-200 experiment

Authors: B. J. Park, J. J. Choe, J. S. Choi, O. Gileva, C. Ha, A. Iltis, E. J. Jeon, D. Y. Kim, K. W. Kim, S. K. Kim, Y. D. Kim, Y. J. Ko, C. H. Lee, H. S. Lee, I. S. Lee, M. H. Lee, S. H. Lee, S. J. Ra, J. K. Son, K. A. Shin

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

Submitted 31 August, 2020; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: 12 pages, 12 figures

Journal ref: Eur. Phys. J. C 80, 814 (2020)

The COSINE-100 Liquid Scintillator Veto System

Authors: G. Adhikari, E. Barbosa de Souza, N. Carlin, J. J. Choi, S. Choi, M. Djamal, A. C. Ezeribe, L. E. Franca, C. Ha, I. S. Hahn, E. J. Jeon, J. H. Jo, W. G. Kang, M. Kauer, H. Kim, H. J. Kim, K. W. Kim, S. K. Kim, Y. D. Kim, Y. H. Kim, Y. J. Ko, E. K. Lee, H. S. Lee, J. Lee, J. Y. Lee , et al. (21 additional authors not shown)

Abstract: This paper describes the liquid scintillator veto system for the COSINE-100 dark matter experiment and its performance. The COSINE-100 detector consists of eight NaI(Tl) crystals immersed in 2200~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 mode… ▽ More This paper describes the liquid scintillator veto system for the COSINE-100 dark matter experiment and its performance. The COSINE-100 detector consists of eight NaI(Tl) crystals immersed in 2200~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. △ Less

Submitted 14 May, 2021; v1 submitted 5 April, 2020; originally announced April 2020.

Journal ref: Nucl. Instrum. Methods A 1006, 165431 (2021)

The AGORA high-resolution galaxy simulations comparison project: Public data release

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Joel R. Primack, Michael J. Butler, Daniel Ceverino, Jun-Hwan Choi, Robert Feldmann, Ben W. Keller, Alessandro Lupi, Kentaro Nagamine, Thomas R. Quinn, Yves Revaz, Romain Teyssier, Spencer C. Wallace

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

Submitted 16 January, 2020; v1 submitted 13 January, 2020; originally announced January 2020.

Comments: AGORA links: http://www.AGORAsimulations.org/ and http://sites.google.com/site/santacruzcomparisonproject/blogs/quicklinks/

Comparison between DAMA/LIBRA and COSINE-100 in the light of Quenching Factors

Authors: Y. J. Ko, K. W. Kim, G. Adhikari, P. Adhikari, E. Barbosa de Souza, N. Carlin, J. J. Choi, S. Choi, M. Djamal, A. C. Ezeribe, C. Ha, I. S. Hahn, E. J. Jeon, J. H. Jo, W. G. Kang, M. Kauer, G. S. Kim, H. Kim, H. J. Kim, N. Y. Kim, S. K. Kim, Y. D. Kim, Y. H. Kim, E. K. Lee, H. S. Lee , et al. (24 additional authors not shown)

Abstract: 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… ▽ More 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$σ$ regions from the DAMA/LIBRA data are still generally excluded by the COSINE-100 data. △ Less

Submitted 23 October, 2019; v1 submitted 10 July, 2019; originally announced July 2019.

Comments: Accepted in JCAP

Journal ref: JCAP 11 (2019) 008

Study of cosmogenic radionuclides in the COSINE-100 NaI(Tl) detectors

Authors: E. Barbosa de Souza, B. J. Park, G. Adhikari, P. Adhikari, N. Carlin, J. J. Choi, S. Choi, M. Djamal, A. C. Ezeribe, C. Ha, I. S. Hahn, E. J. Jeon, J. H. Jo, W. G. Kang, M. Kauer, G. S. Kim, H. Kim, H. J. Kim, K. W. Kim, N. Y. Kim, S. K. Kim, Y. D. Kim, Y. H. Kim, Y. J. Ko, V. A. Kudryavtsev , et al. (26 additional authors not shown)

Abstract: 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… ▽ More 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 $γ$/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. △ Less

Submitted 15 September, 2019; v1 submitted 30 May, 2019; originally announced May 2019.

Self-consistent predictions for LIER-like emission lines from post-AGB stars

Authors: Nell Byler, Julianne J. Dalcanton, Charlie Conroy, Benjamin D. Johnson, Jieun Choi, Aaron Dotter, Philip Rosenfield

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

Submitted 24 April, 2019; originally announced April 2019.

Comments: ApJ accepted. 20 pages, 8 figures

The Imprint of Element Abundance Patterns on Quiescent Galaxy SEDs

Authors: Jieun Choi, Charlie Conroy, Benjamin D. Johnson

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

Submitted 18 January, 2019; originally announced January 2019.

Comments: 9 pages, 5 figures. Accepted to ApJ

A Novel Approach to Constrain Rotational Mixing & Convective-Core Overshoot in Stars Using the Initial-Final Mass Relation

Authors: Jeffrey D. Cummings, Jason S. Kalirai, Jieun Choi, Cyril Georgy, Pier-Emmanuel Tremblay, Enrico Ramirez-Ruiz

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

Submitted 9 January, 2019; originally announced January 2019.

Comments: 7 pages, 5 figures, accepted in ApJ Letters

Optical observations of NEA 3200 Phaethon (1983 TB) during the 2017 apparition

Authors: Myung-Jin Kim, Hee-Jae Lee, Sang-Min Lee, Dong-Heun Kim, Fumi Yoshida, Przemyslaw Bartczak, Grzegorz Dudzinski, Jintae Park, Young-Jun Choi, Hong-Kyu Moon, Hong-Suh Yim, Jin Choi, Eun-Jung Choi, Joh-Na Yoon, Alexander Serebryanskiy, Maxim Krugov, Inna Reva, Kamoliddin E. Ergashev, Otabek Burkhonov, Shuhrat A. Ehgamberdiev, Yunus Turayev, Zhong-Yi Lin, Tomoko Arai, Katsuhito Ohtsuka, Takashi Ito , et al. (2 additional authors not shown)

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

Submitted 16 September, 2018; originally announced September 2018.

Comments: 14 pages, 4 figures, 1 figure in Appendix A. Accepted for publication in Astronomy & Astrophysics (A&A)

Journal ref: A&A 619, A123 (2018)

The White Dwarf Initial-Final Mass Relation for Progenitor Stars From 0.85 to 7.5 M$_\odot$

Authors: Jeffrey D. Cummings, Jason S. Kalirai, P. -E. Tremblay, Enrico Ramirez-Ruiz, Jieun Choi

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

Submitted 5 September, 2018; originally announced September 2018.

Comments: 14 pages, 6 figures, accepted for publication in ApJ

Star Cluster Ages in the Gaia Era

Authors: Jieun Choi, Charlie Conroy, Yuan-Sen Ting, Phillip A. Cargile, Aaron Dotter, Benjamin D. Johnson

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

Submitted 10 July, 2018; originally announced July 2018.

Comments: 25 pages, 18 figures. Accepted to ApJ

Suppression of Star Formation in Low-Mass Galaxies Caused by the Reionization of their Local Neighborhood

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Dominique Aubert, Nicolas Gillet, Jun-Hwan Choi, Ilian T. Iliev, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Anson D'Aloisio, Hyunbae Park, Yehuda Hoffman

Abstract: 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 s… ▽ More 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. △ Less

Submitted 25 July, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 15 pages, 15 figures, accepted by MNRAS 07/17/18, comments welcome

Journal ref: MNRAS, 480, 1740-1753 (2018)

Hyper-Kamiokande Design Report

Authors: Hyper-Kamiokande Proto-Collaboration, :, K. Abe, Ke. Abe, H. Aihara, A. Aimi, R. Akutsu, C. Andreopoulos, I. Anghel, L. H. V. Anthony, M. Antonova, Y. Ashida, V. Aushev, M. Barbi, G. J. Barker, G. Barr, P. Beltrame, V. Berardi, M. Bergevin, S. Berkman, L. Berns, T. Berry, S. Bhadra, D. Bravo-Berguño, F. d. M. Blaszczyk , et al. (291 additional authors not shown)

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

Submitted 28 November, 2018; v1 submitted 9 May, 2018; originally announced May 2018.

Comments: 325 pages

On the Red Giant Branch: Ambiguity in the Surface Boundary Condition Leads to ~100 K Uncertainty in Model Effective Temperatures

Authors: Jieun Choi, Aaron Dotter, Charlie Conroy, Yuan-Sen Ting

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

Submitted 10 May, 2018; originally announced May 2018.

Comments: 9 pages, 5 figures, accepted to ApJ. Figure 5 shows the main result

Age Determinations of the Hyades, Praesepe, and Pleiades via MESA Models with Rotation

Authors: Seth Gossage, Charlie Conroy, Aaron Dotter, Jieun Choi, Philip Rosenfield, Philip Cargile, Andrew Dolphin

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

Submitted 13 July, 2018; v1 submitted 17 April, 2018; originally announced April 2018.

Comments: 20 pages, 12 figures; ApJ Accepted

The Revolution Revolution: magnetic morphology driven spin-down

Authors: C. Garraffo, J. J. Drake, A. Dotter, J. Choi, D. J. Burke, S. P. Moschou, J. D. Alvarado-Gomez, V. L. Kashyap, O. Cohen

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

Submitted 15 June, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 8 pages, 4 figures, to appear in ApJ

A precessing Be disk as a possible model for occultation events in GX 304-1

Authors: Matthias Kühnel, Richard E. Rothschild, Atsuo T. Okazaki, Sebastian Müller, Katja Pottschmidt, Ralf Ballhausen, Jieun Choi, Ingo Kreykenbohm, Felix Fürst, Diana M. Marcu-Cheatham, Paul Hemphill, Macarena Sagredo, Peter Kretschmar, Silvia Martínez-Núñez, José Miguel Torrejón, Rüdiger Staubert, Jörn Wilms

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

Submitted 29 June, 2017; originally announced June 2017.

Comments: 12 pages, 8 figures, 3 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

The Evolution and Properties of Rotating Massive Star Populations

Authors: Jieun Choi, Charlie Conroy, Nell Byler

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

Submitted 15 February, 2017; originally announced February 2017.

Comments: 16 pages, 17 figures. Resubmitted to ApJ following a referee report

A massive, quiescent, population II galaxy at a redshift of 2.1

Authors: Mariska Kriek, Charlie Conroy, Pieter G. van Dokkum, Alice E. Shapley, Jieun Choi, Naveen A. Reddy, Brian Siana, Freeke van de Voort, Alison L. Coil, Bahram Mobasher

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

Submitted 6 December, 2016; originally announced December 2016.

Comments: Published in Nature on 8 December 2016 (5 pages, 3 figures)

The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test

Authors: Ji-hoon Kim, Oscar Agertz, Romain Teyssier, Michael J. Butler, Daniel Ceverino, Jun-Hwan Choi, Robert Feldmann, Ben W. Keller, Alessandro Lupi, Thomas Quinn, Yves Revaz, Spencer Wallace, Nickolay Y. Gnedin, Samuel N. Leitner, Sijing Shen, Britton D. Smith, Robert Thompson, Matthew J. Turk, Tom Abel, Kenza S. Arraki, Samantha M. Benincasa, Sukanya Chakrabarti, Colin DeGraf, Avishai Dekel, Nathan J. Goldbaum , et al. (18 additional authors not shown)

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

Submitted 8 March, 2018; v1 submitted 10 October, 2016; originally announced October 2016.

Comments: 28 pages, 35 figures, Accepted for publication in the Astrophysical Journal, Image resolution greatly reduced, High-resolution version of this article is available at http://www.jihoonkim.org/agora/AGORA_Paper4_draft.pdf, The first paper of the AGORA Initiative is at http://adsabs.harvard.edu/abs/2014ApJS..210...14K, More information on AGORA is at http://www.AGORAsimulations.org/

Journal ref: Astrophysical Journal, 833, 202, 2016

Galaxy Simulation with Dust Formation and Destruction

Authors: Shohei Aoyama, Kuan-Chou Hou, Ikkoh Shimizu, Hiroyuki Hirashita, Keita Todoroki, Jun-Hwan Choi, Kentaro Nagamine

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

Submitted 18 January, 2017; v1 submitted 23 September, 2016; originally announced September 2016.

Comments: 18 pages, 16 figures, published on MNRAS

Journal ref: MNRAS (2017) 466 (1): 105-121

MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models

Authors: Jieun Choi, Aaron Dotter, Charlie Conroy, Matteo Cantiello, Bill Paxton, Benjamin D. Johnson

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

Submitted 28 April, 2016; originally announced April 2016.

Comments: Accepted to ApJ

Search for solar axions with CsI(Tl) crystal detectors

Authors: KIMS Collaboration, Y. S. Yoon, H. K. Park, H. Bhang, J. H. Choi, S. Choi, I. S. Hahn, E. J. Jeon, H. W. Joo, W. G. Kang, B. H. Kim, G. B. Kim, H. J. Kim, K. W. Kim, S. C. Kim, S. K. Kim, Y. D. Kim, Y. H. Kim, H. S. Lee, J. H. Lee, J. K. Lee, S. J. Lee, D. S. Leonard, J. Li, S. S. Myung , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 23 May, 2016; v1 submitted 6 April, 2016; originally announced April 2016.

Comments: 9 pages, 6 figures

A Super-Jupiter Microlens Planet Characterized by High-Cadence KMTNet Microlensing Survey Observations of OGLE-2015-BLG-0954

Authors: I. -G. Shin, Y. -H. Ryu, A. Udalski, M. Albrow, S. -M. Cha, J. -Y. Choi, S. -J. Chung, C. Han, K. -H. Hwang, Y. K. Jung, D. -J. Kim, S. -L. Kim, C. -U. Lee, Y. Lee, B. -G. Park, H. Park, R. W. Pogge, J. C. Yee, P. Pietrukowicz, P. Mróz, S. Kozłowski, R. Poleski, J. Skowron, I. Soszyński, M. K. Szymański , et al. (4 additional authors not shown)

Abstract: 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.… ▽ More 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. △ Less

Submitted 23 June, 2016; v1 submitted 29 February, 2016; originally announced March 2016.

Comments: JKAS, in press, 9 pages, 2 figures

The Hydrodynamic Feedback of Cosmic Reionization on Small-Scale Structures and Its Impact on Photon Consumption during the Epoch of Reionization

Authors: Hyunbae Park, Paul R. Shapiro, Jun-hwan Choi, Naoki Yoshida, Shingo Hirano, Kyungjin Ahn

Abstract: 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 reioniza… ▽ More 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δ$, 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δ]^{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. △ Less

Submitted 17 October, 2016; v1 submitted 20 February, 2016; originally announced February 2016.

Comments: 21 pages, 19 figure, Accepted to ApJ

Spitzer Observations of OGLE-2015-BLG-1212 Reveal a New Path to Breaking Strong Microlens Degeneracies

Authors: V. Bozza, Y. Shvartzvald, A. Udalski, S. Calchi Novati, I. A. Bond, C. Han, M. Hundertmark, R. Poleski, M. Pawlak, M. K. Szymański, J. Skowron, P. Mróz, S. Kozłowski, Ł. Wyrzykowski, P. Pietrukowicz, I. Soszyński, K. Ulaczyk, C. Beichman, G. Bryden, S. Carey, M. Fausnaugh, B. S. Gaudi, A. Gould, C. B. Henderson, R. W. Pogge , et al. (76 additional authors not shown)

Abstract: 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 t… ▽ More 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. △ Less

Submitted 10 February, 2016; v1 submitted 7 January, 2016; originally announced January 2016.

Comments: 27 pages, 6 figures, accepted by ApJ

Ubiquitous Time Variability of Integrated Stellar Populations

Authors: Charlie Conroy, Pieter van Dokkum, Jieun Choi

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

Submitted 12 November, 2015; originally announced November 2015.

Comments: to appear in Nature

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Authors: Pierre Ocvirk, Nicolas Gillet, Paul R. Shapiro, Dominique Aubert, Ilian T. Iliev, Romain Teyssier, Gustavo Yepes, Jun-Hwan Choi, David Sullivan, Alexander Knebe, Stefan Gottloeber, Anson D'Aloisio, Hyunbae Park, Yehuda Hoffman, Timothy Stranex

Abstract: 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 reso… ▽ More 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) △ Less

Submitted 6 September, 2016; v1 submitted 30 October, 2015; originally announced November 2015.

Comments: 26 pages, 16 figures, accepted in MNRAS

Red noise versus planetary interpretations in the microlensing event OGLE-2013-BLG-446

Authors: E. Bachelet, D. M. Bramich, C. Han, J. Greenhill, R. A. Street, A. Gould, G. D Ago, K. AlSubai, M. Dominik, R. Figuera Jaimes, K. Horne, M. Hundertmark, N. Kains, C. Snodgrass, I. A. Steele, Y. Tsapras, M. D. Albrow, V. Batista, J. -P. Beaulieu, D. P. Bennett, S. Brillant, J. A. R. Caldwell, A. Cassan, A. Cole, C. Coutures , et al. (76 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 28 October, 2015; v1 submitted 9 October, 2015; originally announced October 2015.

Comments: accepted ApJ 2015

Mass Measurements of Isolated Objects from Space-based Microlensing

Authors: Wei Zhu, S. Calchi Novati, A. Gould, A. Udalski, C. Han, Y. Shvartzvald, C. Ranc, U. G. Jorgensen, R. Poleski, V. Bozza, C. Beichman, G. Bryden, S. Carey, B. S. Gaudi, C. B. Henderson, R. W. Pogge, I. Porritt, B. Wibking, J. C. Yee, M. Pawlak, M. K. Szymanski, J. Skowron, P. Mroz, S. Kozlowski, L. Wyrzykowski , et al. (56 additional authors not shown)

Abstract: 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$… ▽ More 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. △ Less

Submitted 21 April, 2016; v1 submitted 7 October, 2015; originally announced October 2015.

Comments: 10 papers, 4 figures, 2 tables; ApJ in press

The Baryon Cycle at High Redshifts: Effects of Galactic Winds on Galaxy Evolution in Overdense and Average Regions

Authors: Raphael Sadoun, Isaac Shlosman, Jun-Hwan Choi, Emilio Romano-Díaz

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

Submitted 21 July, 2016; v1 submitted 21 September, 2015; originally announced September 2015.

Comments: 22 pages, 15 figures, accepted for publication in the Astrophysical Journal

Spitzer Parallax of OGLE-2015-BLG-0966: A Cold Neptune in the Galactic Disk

Authors: R. A. Street, A. Udalski, S. Calchi Novati, M. P. G. Hundertmark, W. Zhu, A. Gould, J. Yee, Y. Tsapras, D. P. Bennett, The RoboNet Project, MiNDSTEp Consortium, U. G. Jorgensen, M. Dominik, M. I. Andersen, E. Bachelet, V. Bozza, D. M. Bramich, M. J. Burgdorf, A. Cassan, S. Ciceri, G. D'Ago, Subo Dong, D. F. Evans, Sheng-hong Gu, H. Harkonnen , et al. (84 additional authors not shown)

Abstract: 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 maint… ▽ More 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. △ Less

Submitted 27 August, 2015; originally announced August 2015.

Comments: 28 pages, 3 figures, 2 tables, submitted to ApJ

Spitzer Microlens Measurement of a Massive Remnant in a Well-Separated Binary

Authors: Y. Shvartzvald, A. Udalski, A. Gould, C. Han, V. Bozza, M. Friedmann, M. Hundertmark, C. Beichman, G. Bryden, S. Calchi Novati, S. Carey, M. Fausnaugh, B. S. Gaudi, C. B. Henderson, T. Kerr, R. W. Pogge, W. Varricatt, B. Wibking, J. C. Yee, W. Zhu, R. Poleski, M. Pawlak, M. K. Szymański, J. Skowron, P. Mróz , et al. (41 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 18 December, 2015; v1 submitted 26 August, 2015; originally announced August 2015.

Comments: ApJ. accepted, 34 pages, 7 figures, 2 tables

Supermassive Black Hole Seed Formation at High Redshifts: Long-Term Evolution of the Direct Collapse

Authors: Isaac Shlosman, Jun-Hwan Choi, Mitchell C. Begelman, Kentaro Nagamine

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

Submitted 1 December, 2015; v1 submitted 20 August, 2015; originally announced August 2015.

Comments: 14 pages, 11 figures, MNRAS, in press, typos corrected

KMT-2015-1b: a Giant Planet Orbiting a Low-mass Dwarf Host Star Discovered by a New High-cadence Microlensing Survey with a Global Telescope Network

Authors: K. -H. Hwang, C. Han, J. -Y. Choi, H. Park, Y. K. Jung, I. -G. Shin, M. D. Albrow, A. Gould, V. Bozza, B. -G. Park, S. -L. Kim, C. -U. Lee, S. -M. Cha, D. -J. Kim, Y. Lee

Abstract: 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… ▽ More 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. △ Less

Submitted 19 July, 2015; originally announced July 2015.

Comments: 3 figures and 1 table