We analyze KMT-2019-BLG-1339, a microlensing event with an obvious but incompletely resolved brief anomaly feature around the peak of the light curve. Although the origin of the anomaly is identified to be a companion to the lens with a low mass ratio $q$, the interpretation is subject to two different degeneracy types. The first type is the ambiguity in $\rho$, representing the angular source radius scaled to the angular radius of the Einstein ring, $\theta_{\rm E}$, and the other is the $s\leftrightarrow s^{-1}$ degeneracy. The former type, `finite-source degeneracy', causes ambiguities in both $s$ and $q$, while the latter induces an ambiguity only in $s$. Here $s$ denotes the separation (in units of $\theta_{\rm E}$) in projection between the lens components. We estimate that the lens components have masses $(M_1, M_2)\sim (0.27^{+0.36}_{-0.15}~M_\odot, 11^{+16}_{-7}~M_{\rm J})$ and $\sim (0.48^{+0.40}_{-0.28}~M_\odot, 1.3^{+1.1}_{-0.7}~M_{\rm J})$ according to the two solutions subject to the finite-source degeneracy, indicating that the lens comprises an M dwarf and a companion with a mass around the planet/brown dwarf boundary or a Jovian-mass planet. It is possible to lift the finite-source degeneracy by conducting future observations utilizing a high resolution instrument because the relative lens-source proper motion predicted by the solutions are widely different.

We present the analysis of a very high-magnification ($A\sim 900$) microlensing event KMT-2019-BLG-1953. A single-lens single-source (1L1S) model appears to approximately delineate the observed light curve, but the residuals from the model exhibit small but obvious deviations in the peak region. A binary lens (2L1S) model with a mass ratio $q\sim 2\times 10^{-3}$ improves the fits by $\Delta\chi^2=181.8$, indicating that the lens possesses a planetary companion. From additional modeling by introducing an extra planetary lens component (3L1S model) and an extra source companion (2L2S model), it is found that the residuals from the 2L1S model further diminish, but claiming these interpretations is difficult due to the weak signals with $\Delta\chi^2=16.0$ and $13.5$ for the 3L1S and 2L2L models, respectively. From a Bayesian analysis, we estimate that the host of the planets has a mass of $M_{\rm host}=0.31^{+0.37}_{-0.17}~M_\odot$ and that the planetary system is located at a distance of $D_{\rm L}=7.04^{+1.10}_{-1.33}~{\rm kpc}$ toward the Galactic center. The mass of the securely detected planet is $M_{\rm p}=0.64^{+0.76}_{-0.35}~M_{\rm J}$. The signal of the potential second planet could have been confirmed if the peak of the light curve had been more densely observed by followup observations, and thus the event illustrates the need for intensive followup observations for very high-magnification events even in the current generation of high-cadence surveys.

We announce the discovery of a microlensing planetary system, in which a sub-Saturn planet is orbiting an ultracool dwarf host. We detect the planetary system by analyzing the short-timescale ($t_{\rm E}\sim 4.4$~days) lensing event KMT-2018-BLG-0748. The central part of the light curve exhibits asymmetry due to the negative deviations in the rising part and the positive deviations in the falling part. We find that the deviations are explained by a binary-lens model with a mass ratio between the lens components of $q\sim 2\times 10^{-3}$. The short event timescale together with the small angular Einstein radius, $\theta_{\rm E}\sim 0.11$~mas, indicate that the mass of the planet host is very small. The Bayesian analysis conducted under the assumption that the planet frequency is independent of the host mass indicates that the mass of the planet is $M_{\rm p}=0.18^{+0.29}_{-0.10}~M_{\rm J}$, and the mass of the host, $M_{\rm h}= 0.087^{+0.138}_{-0.047}~M_\odot$, is near the star/brown dwarf boundary, but the estimated host mass is sensitive to the assumption about the planet hosting probability. High-resolution follow-up observations would lead to revealing the nature of the planet host.

Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$\sigma$. The ten years of data post-trial significance is 4.1$\sigma$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented.

Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger.

We report the discovery of a planet in the microlensing event OGLE-2018-BLG-1269, with planet-host mass ratio $q \sim 6\times10^{-4}$, i.e., $0.6$ times smaller than the Jupiter/Sun mass ratio. Combined with the $Gaia$ parallax and proper motion, a strong one-dimensional constraint on the microlens parallax vector allows us to significantly reduce the uncertainties of lens physical parameters. A Bayesian analysis that ignores any information about light from the host yields that the planet is a cold giant $(M_{2} = 0.69_{-0.22}^{+0.44}\,M_{\rm J})$ orbiting a Sun-like star $(M_{1} = 1.13_{-0.35}^{+0.72}\,M_{\odot})$ at a distance of $D_{\rm L} = 2.56_{-0.62}^{+0.92}\,{\rm kpc}$. The projected planet-host separation is $a_{\perp} = 4.61_{-1.17}^{+1.70}\,{\rm au}$. Using {\it Gaia} astrometry, we show that the blended light lies $\lesssim 12\,$mas from the host and therefore must be either the host star or a stellar companion to the host. An isochrone analysis favors the former possibility at $>99.6\%$. The host is therefore a subgiant. For host metallicities in the range of $0.0 \leq {\rm [Fe/H]} \leq +0.3$, the host and planet masses are then in the range of $1.16 \leq M_{1}/M_{\odot} \leq 1.38$ and $0.74 \leq M_{2}/M_{\rm J} \leq 0.89$, respectively. Low host metallicities are excluded. The brightness and proximity of the lens make the event a strong candidate for spectroscopic followup both to test the microlensing solution and to further characterize the system.

Stellar systems are broadly divided into collisional and non-collisional. The latter are large-N systems with long relaxation timescales and can be simulated disregarding two-body interactions, while either computationally expensive direct N-body simulations or approximate schemes are required to properly model the former. Large globular clusters and nuclear star clusters, with relaxation timescales of the order of a Hubble time, are small enough to display some collisional behaviour and big enough to be impossible to simulate with direct $N$-body codes and current hardware. We introduce a new method to simulate collisional stellar systems, and validate it by comparison with direct N-body codes on small-N simulations. The Multi-Particle collision for Dense stellar systems Code (MPCDSS) is a new code for evolving stellar systems with the Multi-Particle Collision method. Such method amounts to a stochastic collision rule that allows to conserve exactly the energy and momentum over a cluster of particles experiencing the collision. The code complexity scales with $N \log N$ in the number of particles. Unlike Monte-Carlo codes, MPCDSS can easily model asymmetric, non-homogeneous, unrelaxed and rotating systems, while allowing us to follow the orbits of individual stars. We evolve small ($N = 3.2 \times 10^4$) star clusters with MPCDSS and with the direct-summation code NBODY6, finding a similar evolution of key indicators. We then simulate different initial conditions in the $10^4 - 10^6$ star range. MPCDSS bridges the gap between small, collisional systems that can be simulated with direct $N$-body codes and large noncollisional systems. MPCDSS in principle allows us to simulate globular clusters such as Omega Cen and M54 and even the nuclear star cluster, beyond the limits of current direct N-body codes in terms of the number of particles.

The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $\mu$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $\mu$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $\sigma^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.

We report the discovery and analysis of the planetary microlensing event OGLE-2017-BLG-0406, which was observed both from the ground and by the ${\it Spitzer}$ satellite in a solar orbit. At high magnification, the anomaly in the light curve was densely observed by ground-based-survey and follow-up groups, and it was found to be explained by a planetary lens with a planet/host mass ratio of $q=7.0 \times 10^{-4}$ from the light-curve modeling. The ground-only and ${\it Spitzer}$-"only" data each provide very strong one-dimensional (1-D) constraints on the 2-D microlens parallax vector $\bf{\pi_{\rm E}}$. When combined, these yield a precise measurement of $\bf{\pi_{\rm E}}$, and so of the masses of the host $M_{\rm host}=0.56\pm0.07\,M_\odot$ and planet $M_{\rm planet} = 0.41 \pm 0.05\,M_{\rm Jup}$. The system lies at a distance $D_{\rm L}=5.2 \pm 0.5 \ {\rm kpc}$ from the Sun toward the Galactic bulge, and the host is more likely to be a disk population star according to the kinematics of the lens. The projected separation of the planet from the host is $a_{\perp} = 3.5 \pm 0.3 \ {\rm au}$, i.e., just over twice the snow line. The Galactic-disk kinematics are established in part from a precise measurement of the source proper motion based on OGLE-IV data. By contrast, the ${\it Gaia}$ proper-motion measurement of the source suffers from a catastrophic $10\,\sigma$ error.

We performed simultaneous observations of the H2O 6(1,6) - 5(2,3) (22.235080 GHz) and SiO v= 1, 2, J = 1 - 0, SiO v = 1, J = 2 - 1, 3 - 2 (43.122080, 42.820587, 86.243442, and 129.363359 GHz) masers towards the suspected D-type symbiotic star, V627 Cas, using the Korean VLBI Network. Here, we present astrometrically registered maps of the H2O and SiO v = 1, 2, J = 1 - 0, SiO v = 1, J = 2 - 1 masers for five epochs from January 2016 to June 2018. Distributions of the SiO maser spots do not show clear ring-like structures, and those of the H2O maser are biased towards the north-north-west to west with respect to the SiO maser features according to observational epochs. These asymmetric distributions of H2O and SiO masers are discussed based on two scenarios of a bipolar outflow and the presence of the hot companion, a white dwarf, in V627 Cas. We carried out ring fitting of SiO v = 1, and v = 2 masers and estimated the expected position of the cool red giant. The ring radii of the SiO v = 1 maser are slightly larger than those of the SiO v = 2 maser, as previously known. Our assumption for the physical size of the SiO maser ring of V627 Cas to be the typical size of a SiO maser ring radius (\sim4 au) of red giants yields the distance of V627 Cas to be \sim1 kpc.

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

We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the properties of the shower front and the lateral distribution function.

High-cadence observations of the Galactic bulge by the microlensing surveys led to the discovery of a handful of extremely short-timescale microlensing events that can be attributed to free-floating or wide-orbit planets. Here, we report the discovery of another strong free-floating planet candidate, which was found from the analysis of the gravitational microlensing event OGLE-2019-BLG-0551. The light curve of the event is characterized by a very short duration (<3 d) and a very small amplitude (< 0.1 mag). From modeling of the light curve, we find that the Einstein timescale, tE = 0.381 +/- 0.017 d, is much shorter, and the angular Einstein radius, thetaE = 4.35 +/- 0.34 uas, is much smaller than those of typical lensing events produced by stellar-mass lenses (tE ~ 20 d, thetaE ~ 0.3 mas), indicating that the lens is very likely to be a planetary-mass object. We conduct an extensive search for possible signatures of a companion star in the light curve of the event, finding no significant evidence for the putative host star. For the first time, we also demonstrate that the angular Einstein radius of the lens does not depend on blending in the low-magnification events with strong finite source effects.

We report the analysis of the microlensing event OGLE-2018-BLG-0677. A small feature in the light curve of the event leads to the discovery that the lens is a star-planet system. Although there are two degenerate solutions that could not be distinguished for this event, both lead to a similar planet-host mass ratio. We perform a Bayesian analysis based on a Galactic model to obtain the properties of the system and find that the planet corresponds to a super-Earth/sub-Neptune with a mass $M_{\mathrm{planet}} = {3.96}^{+5.88}_{-2.66}\mathrm{M_\oplus}$. The host star has a mass $ M_{\mathrm{host}} = {0.12}^{+0.14}_{-0.08}\mathrm{M_\odot}$. The projected separation for the inner and outer solutions are ${0.63}^{+0.20}_{-0.17}$~AU and ${0.72}^{+0.23}_{-0.19}$~AU respectively. At $\Delta\chi^2=\chi^2({\rm 1L1S})-\chi^2({\rm 2L1S})=46$, this is by far the lowest $\Delta\chi^2$ for any securely-detected microlensing planet to date, a feature that is closely connected to the fact that it is detected primarily via a "dip" rather than a "bump".

We present a new catalog and results for the cluster system of the starburst galaxy NGC 4449 based on multi-band imaging observations taken as part of the LEGUS and Halpha-LEGUS surveys. We improve the spectral energy fitting method used to estimate cluster ages and find that the results, particularly for older clusters, are in better agreement with those from spectroscopy. The inclusion of Halpha measurements, the role of stochasticity for low mass clusters, the assumptions about reddening, and the choices of SSP model and metallicity all have important impacts on the age-dating of clusters. A comparison with ages derived from stellar color-magnitude diagrams for partially resolved clusters shows reasonable agreement, but large scatter in some cases. The fraction of light found in clusters relative to the total light (i.e., T_L) in the U, B, and V filters in 25 different ~kpc-size regions throughout NGC 4449 correlates with both the specific Region Luminosity, R_L, and the dominant age of the underlying stellar population in each region. The observed cluster age distribution is found to decline over time as dN/dt ~ t^g, with g=-0.85+/-0.15, independent of cluster mass, and is consistent with strong, early cluster disruption. The mass functions of the clusters can be described by a power law with dN/dM ~ M^b and b=-1.86+/-0.2, independent of cluster age. The mass and age distributions are quite resilient to differences in age-dating methods. There is tentative evidence for a factor of 2-3 enhancement in both the star and cluster formation rate ~100 - 300 Myr ago, indicating that cluster formation tracks star formation generally. The enhancement is probably associated with an earlier interaction event.

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

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

We present a family of new rotating black hole solutions to Einstein's equations that generalizes the Kerr-Newman spacetime to include an anisotropic matter. The geometry is obtained by employing the Newman-Janis algorithm. In addition to the mass, the charge and the angular momentum, an additional hair exists thanks to the negative radial pressure of the anisotropic matter. The properties of the black hole are analyzed in detail including thermodynamics. This black hole can be used as a better engine than the Kerr-Newman one in extracting energy.

The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper.

According to the disk accretion model theory of Shakura and Sunyayev (1973), accretion disk surrounding the spinning astronomical compact objects (BH, NS) consists of fluid rings obeying the Keplerian motions. Therefore, as we approach the center of compact objects, angular momentum drops, but angular velocity and hence the dissipative heat and accretion rate grow. As a result, the shape of the accretion disk "puffs up" (gets thicker) near the surface of the spinning-compact objects. We, therefore, illustrate the geometry of this thick inner portion of the accretion disk in the context of Mukhopadhyyay's Pseudo Newtonian potential approximation for full General Theory of Relativity.

A few dozen dwarf satellite galaxies of the Milky Way have been discovered, which are often viewed as the remaining building blocks of our Galaxy. The follow-up spectroscopy showed that dwarf galaxies have a sizeable spread in their metallicities. Several scenarios were suggested to explain the metallicity spread, which can be tested by the structural patterns of stellar subpopulations with distinct metallicities. However, such chemical plus structural examination, to which we refer to as "chemostructural study" is hindered by the lack of stars with spectroscopic metallicity. Here we propose the Ca$-$$by$ photometry as an alternative way to secure metallicities for a 2$-$3 orders of magnitude larger stellar sample than the spectroscopic sample and thus enable us to perform a chemostructural study on dwarf galaxies. In particular, we use the $hk$ index [$\equiv($Ca$-b)-(b-y)$], whose validity as a photometric metallicity indicator (and crass insensitivity to age) for red-giant-branch stars was upheld via Galactic globular clusters, and observe three dwarf spheroidal galaxies$-$Draco, Sextans, and Canes Venatici I (CVnI)$-$with Subaru/Suprime-Cam. We find that in all the galaxies the metal-rich stellar populations are more centrally concentrated than the metal-poor counterparts, suggesting that the central regions of the galaxies underwent extended star formation. Such a negative radial metallicity gradient for Sextans and CVnI opposes to the traditional spectroscopic results. We also find that their metallicity distribution functions (MDFs) can be characterized by a unimodal, skewed Gaussian shape with a metal-rich peak and a metal-poor tail. We discuss their features in the chemostructure and MDF in terms of dwarf galaxy formation theories.

Only about 19 Galactic and 25 extra-galactic bona-fide Luminous Blue Variables (LBVs) are known to date. This incomplete census prevents our understanding of this crucial phase of massive star evolution which leads to the formation of heavy binary black holes via the classical channel. With large samples of LBVs one could better determine the duration and maximum stellar luminosity which characterize this phase. We search for candidate LBVs (cLBVs) in a new galaxy, NGC 7793. For this purpose, we combine high spatial resolution images from two Hubble Space Telescope (HST) programs with optical spectroscopy from the Multi Unit Spectroscopic Explorer (MUSE). By combining PSF-fitting photometry measured on F547M, F657N, and F814W images, with restrictions on point-like appearance (at HST resolution) and H{\alpha} luminosity, we find 100 potential cLBVs, 36 of which fall in the MUSE fields. Five of the latter 36 sources are promising cLBVs which have M(V) less than or equal to -7 and a combination of: H{\alpha} with a P-Cygni profile; no [O I] 6300 emission; weak or no [O III] 5007 emission; large [N II]/H{\alpha} relative to H II regions; and [S II] 6716 / [S II] 6731 similar to 1. It is not clear if these five cLBVs are isolated from O-type stars, which would favor the binary formation scenario of LBVs. Our study, which approximately covers one fourth of the optical disc of NGC 7793, demonstrates how by combining the above HST surveys with multi-object spectroscopy from 8-m class telescopes, one can efficiently find large samples of cLBVs in nearby galaxies.

At $q=1.81\pm 0.20 \times 10^{-5}$, KMT-2018-BLG-0029Lb has the lowest planet-host mass ratio $q$ of any microlensing planet to date by more than a factor of two. Hence, it is the first planet that probes below the apparent "pile-up" at $q=5$--10 $\times 10^{-5}$. The event was observed by {\it Spitzer}, yielding a microlens-parallax $\pi_{\rm E}$ measurement. Combined with a measurement of the Einstein radius $\theta_{\rm E}$ from finite-source effects during the caustic crossings, these measurements imply masses of the host $M_{\rm host}=1.14^{+0.10}_{-0.12}\, M_\odot$ and planet $M_{\rm planet} = 7.59^{+0.75}_{-0.69}\,M_\oplus$, system distance $D_L = 3.38^{+0.22}_{-0.26}\,\,{\rm kpc}$ and projected separation $a_\perp = 4.27^{+0.21}_{-0.23}\,{\rm au}$. The blended light, which is substantially brighter than the microlensed source, is plausibly due to the lens and could be observed at high resolution immediately.

We report the discovery of a super-Jovian planet in the microlensing event KMT-2016-BLG-1836, which was found by the Korea Microlensing Telescope Network's high-cadence observations (\Gamma ~ 4~{hr}^{-1}). The planet-host mass ratio q ~ 0.004. A Bayesian analysis indicates that the planetary system is composed of a super-Jovian M_{planet} = 2.2_{-1.1}^{+1.9} M_{J} planet orbiting an M or K dwarf M_{\rm host} = 0.49_{-0.25}^{+0.38} M_{Sun}, at a distance of D_{L} = 7.1_{-2.4}^{+0.8} kpc. The projected planet-host separation is 3.5^{+1.1}_{-0.9} AU, implying that the planet is located beyond the snowline of the host star. Future high-resolution images can potentially strongly constrain the lens brightness and thus the mass and distance of the planetary system. Without considering detailed detection efficiency, selection or publication biases, we find a potential "mass ratio desert" at -3.7 \lesssim \log q \lesssim -3.0 for the 31 published KMTNet planets.

The relative average minimum projected separations of star clusters in the Legacy ExtraGalactic UV Survey (LEGUS) and in tidal dwarfs around the interacting galaxy NGC 5291 are determined as a function of cluster mass to look for cluster-cluster mass segregation. Class 2 and 3 LEGUS clusters, which have a more irregular internal structure than the compact and symmetric class 1 clusters, are found to be mass segregated in low mass galaxies, which means that the more massive clusters are systematically bunched together compared to the lower mass clusters. This mass segregation is not present in high-mass galaxies nor for class 1 clusters. We consider possible causes for this segregation including differences in cluster formation and scattering in the shallow gravitational potentials of low mass galaxies.

We discuss the sensitivity of the present and near-future axion dark matter experiments to a halo of axions or axion-like particles gravitationally bound to the Earth or the Sun. The existence of such halos, assuming they are formed, renders a significant gain in the sensitivity of axion searches while satisfying all the present experimental bounds. The structure and coherence properties of these halos also imply novel signals, which can depend on the latitude or orientation of the detector. We demonstrate this by analysing the sensitivity of several distinct types of axion dark matter experiments.

With $q=(4.09\pm 0.27) \times 10^{-5}$, KMT-2019-BLG-0842Lb is the second cold planet with planet-host mass ratio $q$ below that of Uranus relative to the Sun ($q=4.37\times 10^{-5}$), the first having been discovered only in 2018. Hence, such cold, sub-Uranus-ratio planets are probably not rare. We discuss possible improvements in the search for these planets. The Bayesian estimates for the host mass, planet mass, system distance, and planet-host projected separation are $M_{\rm host}=0.76\pm 0.40 M_\odot$, $M_{\rm planet}=10.3\pm 5.5 M_\oplus$, $D_L = 3.3\pm 1.3$ kpc, and $a_\perp = 3.3\pm 1.4$ AU, respectively. The blended light is plausibly due to the host star, an identification that could be tested by high resolution images taken during the next (2020) bulge season.

Cassiopeia A is a nearby young supernova remnant that provides a unique laboratory for the study of core-collapse supernova explosions. Cassiopeia A is known to be a Type IIb supernova from the optical spectrum of its light echo, but the immediate progenitor of the supernova remains uncertain. Here we report results of near-infrared, high-resolution spectroscopic observations of Cassiopeia A where we detected the pristine circumstellar material of the supernova progenitor. Our observations revealed a strong emission line of iron (Fe) from a circumstellar clump that has not yet been processed by the supernova shock wave. A comprehensive analysis of the observed spectra, together with an HST image, indicates that the majority of Fe in this unprocessed circumstellar material is in the gas phase, not depleted onto dust grains as in the general interstellar medium. This result is consistent with a theoretical model of dust condensation in material that is heavily enriched with CNO-cycle products, supporting the idea that the clump originated near the He core of the progenitor. It has been recently found that Type IIb supernovae can result from the explosion of a blue supergiant with a thin hydrogen envelope, and our results support such a scenario for Cassiopeia A.

The cosmological relaxion can address the hierarchy problem, while its coherent oscillations can constitute dark matter in the present universe. We consider the possibility that the relaxion forms gravitationally bound objects that we denote as relaxion stars. The density of these stars would be higher than that of the local dark matter density, resulting in enhanced signals in table-top detectors, among others. Furthermore, we raise the possibility that these objects may be trapped by an external gravitational potential, such as that of the Earth or the Sun. This leads to formation of relaxion halos of even greater density. We discuss several interesting implications of relaxion halos, as well as detection strategies to probe them.

We present the analysis of the microlensing event OGLE-2016-BLG-1227. The light curve of this short-duration event appears to be a single-lens event affected by severe finite-source effects. Analysis of the light curve based on single-lens single-source (1L1S) modeling yields very small values of the event timescale, $t_{\rm E}\sim 3.5$ days, and the angular Einstein radius, $\theta_{\rm E}\sim 0.009$ mas, making the lens a candidate of a free-floating planet. Close inspection reveals that the 1L1S solution leaves small residuals with amplitude $\Delta I\lesssim 0.03$ mag. We find that the residuals are explained by the existence of an additional widely-separated heavier lens component, indicating that the lens is a wide-separation planetary system rather than a free-floating planet. From Bayesian analysis, it is estimated that the planet has a mass of $M_{\rm p} = 0.79^{+1.30}_{-0.39} M_{\rm J}$ and it is orbiting a low-mass host star with a mass of $M_{\rm host}=0.10^{+0.17}_{-0.05} M_\odot$ located with a projected separation of $a_\perp=3.4^{+2.1}_{-1.0}$ au. The planetary system is located in the Galactic bulge with a line-of-sight separation from the source star of $D_{\rm LS}=1.21^{+0.96}_{-0.63}$ kpc. The event shows that there are a range of deviations in the signatures of host stars for apparently isolated planetary lensing events and that it is possible to identify a host even when a deviation is subtle.

Exothermic scatterings of dark matter (DM) produce DM particles with significant kick velocities inside DM halos. In collaboration with DM self-interaction, the excess kinetic energy of the produced DM particles is distributed to the others, which self-heats the DM particles as a whole. The DM self-heating is efficient towards the halo center, and the heat injection is used to enhance the formation of a uniform density core inside halos. The effect of DM self-heating is expected to be more significant in smaller halos for two reasons: 1) the exothermic cross section times the relative velocity, $\left\langle\sigma_{\rm exo}v_{\rm rel}\right\rangle$, is constant; 2) and the ratio of the injected heat to the velocity dispersion squared gets larger towards smaller halos. For the first time, we quantitatively investigate the core formation from DM self-heating for halos in a wide mass range ($10^{9}$-$10^{15}\,{\rm M}_\odot$) using the gravothermal fluid formalism. Notably, we demonstrate that the core formation is sharply escalating towards smaller halos by taking the self-heating DM (i.e., DM that semi-annihilates and self-interacts) as an example. We show that the sharp escalation of core formation may cause a tension in simultaneously explaining the observed central mass deficit of Milky Way satellites, and field dwarf/low surface brightness spiral galaxies. While the details of the self-heating effect may differ among models, we expect that the sharp halo-mass dependence of the core formation is a general feature of exothermic DM.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

We report the discovery and the analysis of the short (tE < 5 days) planetary microlensing event, OGLE-2015-BLG-1771. The event was discovered by the Optical Gravitational Lensing Experiment (OGLE), and the planetary anomaly (at I ~ 19) was captured by The Korea Microlensing Telescope Network (KMTNet). The event has three surviving planetary models that explain the observed light curves, with planet-host mass ratio q \~ 5.4 * 10^{-3}, 4.5 * 10^{-3} and 4.5 * 10^{-2}, respectively. The first model is the best-fit model, while the second model is disfavored by \Delta\chi^2 ~ 3. The last model is strongly disfavored by \Delta\chi^2 ~ 15 but not ruled out. A Bayesian analysis using a Galactic model indicates that the first two models are probably composed of a Saturn-mass planet orbiting a late M dwarf, while the third one could consist of a super-Jovian planet and a mid-mass brown dwarf. The source-lens relative proper motion is mu_rel ~ 9 mas/yr, so the source and lens could be resolved by current adaptive-optics (AO) instruments in 2021 if the lens is luminous.

Short-timescale microlensing events are likely to be produced by substellar brown dwarfs (BDs), but it is difficult to securely identify BD lenses based on only event timescales $t_{\rm E}$ because short-timescale events can also be produced by stellar lenses with high relative lens-source proper motions. In this paper, we report three strong candidate BD-lens events found from the search for lensing events not only with short timescales ($t_{\rm E} \lesssim 6~{\rm days}$) but also with very small angular Einstein radii ($\theta_{\rm E}\lesssim 0.05~{\rm mas}$) among the events that have been found in the 2016--2019 observing seasons. These events include MOA-2017-BLG-147, MOA-2017-BLG-241, and MOA-2019-BLG-256, in which the first two events are produced by single lenses and the last event is produced by a binary lens. From the Bayesian analysis conducted with the combined $t_{\rm E}$ and $\theta_{\rm E}$ constraint, it is estimated that the lens masses of the individual events are $0.051^{+0.100}_{-0.027}~M_\odot$, $0.044^{+0.090}_{-0.023}~M_\odot$, and $0.046^{+0.067}_{-0.023}~M_\odot/0.038^{+0.056}_{-0.019}~M_\odot$ and the probability of the lens mass smaller than the lower limit of stars is $\sim 80\%$ for all events. We point out that routine lens mass measurements of short time-scale lensing events require survey-mode space-based observations.

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

We present a detailed near-infrared chemical abundance analysis of 10 red giant members of the Galactic open cluster NGC 752. High-resolution (R$\simeq$45000) near-infrared spectral data were gathered with the Immersion Grating Infrared Spectrograph (IGRINS), providing simultaneous coverage of the complete H and K bands. We derived the abundances of H-burning (C, N, O), $\alpha$ (Mg, Si, S, Ca), light odd-Z (Na, Al, P, K), Fe-group (Sc, Ti, Cr, Fe, Co, Ni) and neutron-capture (Ce, Nd, Yb) elements. We report the abundances of S, P, K, Ce, and Yb in NGC 752 for the first time. Our analysis yields solar metallicity and solar abundance ratios for almost all of the elements heavier than the CNO group in NGC 752. O and N abundances were measured from a number of OH and CN features in the $H$ band, and C abundances were determined mainly from CO molecular lines in the K band. High excitation \ion{C}{i} lines present in both near-infrared and optical spectra were also included in the C abundance determinations. Carbon isotopic ratios were derived from the R-branch band heads of first overtone (2-0) and (3$-$1) $^{12}$CO and (2-0) $^{13}$CO lines near 23440 \AA and (3-1) $^{13}$CO lines at about 23730 \AA. The CNO abundances and $^{12}$C/$^{13}$C ratios are all consistent with our giants having completed "first dredge-up" envelope mixing of CN-cyle products. We independently assessed NGC 752 stellar membership from Gaia astrometry, leading to a new color-magnitude diagram for this cluster. Applications of Victoria isochrones and MESA models to these data yield an updated NGC 752 cluster age (1.52 Gyr) and evolutionary stage indications for the program stars. The photometric evidence and spectroscopic light element abundances all suggest that the most, perhaps all of the program stars are members of the helium-burning red clump in this cluster.

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

The morphology of HII regions around young star clusters provides insight into the timescales and physical processes that clear a cluster's natal gas. We study ~700 young clusters (<10Myr) in three nearby spiral galaxies (NGC 7793, NGC 4395, and NGC 1313) using Hubble Space Telescope (HST) imaging from LEGUS (Legacy ExtraGalactic Ultraviolet Survey). Clusters are classified by their H$\alpha$ morphology (concentrated, partially exposed, no-emission) and whether they have neighboring clusters (which could affect the clearing timescales). Through visual inspection of the HST images, and analysis of ages, reddenings, and stellar masses from spectral energy distributions fitting, together with the (U-B), (V-I) colors, we find: 1) the median ages indicate a progression from concentrated (~3 Myr), to partially exposed (~4 Myr), to no H$\alpha$ emission (>5Myr), consistent with the expected temporal evolution of HII regions and previous results. However, 2) similarities in the age distributions for clusters with concentrated and partially exposed H$\alpha$ morphologies imply a short timescale for gas clearing (<1Myr). 3) our cluster sample's median mass is ~1000 M, and a significant fraction (~20%) contain one or more bright red sources (presumably supergiants), which can mimic reddening effects. Finally, 4) the median E(B-V) values for clusters with concentrated H$\alpha$ and those without H$\alpha$ emission appear to be more similar than expected (~0.18 vs. ~0.14, respectively), but when accounting for stochastic effects, clusters without H$\alpha$ emission are less reddened. To mitigate stochastic effects, we experiment with synthesizing more massive clusters by stacking fluxes of clusters within each H$\alpha$ morphological class. Composite isolated clusters also reveal a color and age progression for H$\alpha$ morphological classes, consistent with analysis of the individual clusters.

In this work, BV RI light curves of 55 Type II supernovae (SNe II) from the Lick Observatory Supernova Search program obtained with the Katzman Automatic Imaging Telescope and the 1 m Nickel telescope from 2006 to 2018 are presented. Additionally, more than 150 spectra gathered with the 3 m Shane telescope are published. We conduct an analyse of the peak absolute magnitudes, decline rates, and time durations of different phases of the light and colour curves. Typically, our light curves are sampled with a median cadence of 5.5 days for a total of 5093 photometric points. In average V-band plateau declines with a rate of 1.29 mag (100 days)-1, which is consistent with previously published samples. For each band, the plateau slope correlates with the plateau length and the absolute peak magnitude: SNe II with steeper decline have shorter plateau duration and are brighter. A time-evolution analysis of spectral lines in term of velocities and pseudoequivalent widths is also presented in this paper. Our spectroscopic sample ranges between 1 and 200 days post-explosion and has a median ejecta expansion velocity at 50 days post-explosion of 6500 km/s (Halpha line) and a standard dispersion of 2000 km/s. Nebular spectra are in good agreement with theoretical models using a progenitor star having a mass <16 Msol. All the data are available to the community and will help to understand SN II diversity better, and therefore to improve their utility as cosmological distance indicators.

We present BVRI and unfiltered light curves of 93 Type Ia supernovae (SNe Ia) from the Lick Observatory Supernova Search (LOSS) follow-up program conducted between 2005 and 2018. Our sample consists of 78 spectroscopically normal SNe Ia, with the remainder divided between distinct subclasses (three SN 1991bg-like, three SN 1991T-like, four SNe Iax, two peculiar, and three super-Chandrasekhar events), and has a median redshift of 0.0192. The SNe in our sample have a median coverage of 16 photometric epochs at a cadence of 5.4 days, and the median first observed epoch is ~4.6 days before maximum B-band light. We describe how the SNe in our sample are discovered, observed, and processed, and we compare the results from our newly developed automated photometry pipeline to those from the previous processing pipeline used by LOSS. After investigating potential biases, we derive a final systematic uncertainty of 0.03 mag in BVRI for our dataset. We perform an analysis of our light curves with particular focus on using template fitting to measure the parameters that are useful in standardising SNe Ia as distance indicators. All of the data are available to the community, and we encourage future studies to incorporate our light curves in their analyses.

We report two microlensing planet candidates discovered by the KMTNet survey in $2017$. However, both events have the 2L1S/1L2S degeneracy, which is an obstacle to claiming the discovery of the planets with certainty unless the degeneracy can be resolved. For KMT-2017-BLG-0962, the degeneracy cannot be resolved. If the 2L1S solution is correct, KMT-2017-BLG-0962 might be produced by a super Jupiter-mass planet orbiting a mid-M dwarf host star. For KMT-2017-BLG-1119, the light curve modeling favors the 2L1S solution but higher-resolution observations of the baseline object tend to support the 1L2S interpretation rather than the planetary interpretation. This degeneracy might be resolved by a future measurement of the lens-source relative proper motion. This study shows the problem of resolving 2L1S/1L2S degeneracy exists over a much wider range of conditions than those considered by the theoretical study of Gaudi (1998).

We present the Strategic Scientific Plan (SSP) for the direction and activities of the Gemini Observatory in the 2020s. The overarching goal is to ensure that Gemini best uses the available resources to serve the needs of its international user community throughout the coming decade. The actionable items fall into three general categories: (1) preserving Gemini's current facilities and strengths; (2) developing instrumentation and software systems, including data pipelines, to enable new scientific capabilities that build on those strengths; (3) strategizing how visiting instruments can deliver additional valuable capabilities. We provide a high-level timeline (schematically illustrated in one figure) for the main developments discussed in this SSP. The schedule is ambitious, but in light of the recent Gemini in the Era of Multi-Messenger Astronomy (GEMMA) award from the NSF, the plan becomes achievable. Lists of milestones are given for gauging progress. As these milestones are reached and new instruments become available, some current instruments will need to be retired; we make recommendations in this regard. The final section concludes by reemphasizing the importance of a strong partnership committed to the needs of all members.

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

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

We report the discovery of a planet in a binary that was discovered from the analysis of the microlensing event OGLE-2018-BLG-1700. We identify the triple nature of the lens from the fact that the complex anomaly pattern can be decomposed into two parts produced by two binary-lens events, in which one binary pair has a very low mass ratio of $\sim 0.01$ between the lens components and the other pair has a mass ratio of $\sim 0.3$. We find two sets of degenerate solutions, in which one solution has a projected separation between the primary and its stellar companion less than the angular Einstein radius $\thetae$ (close solution), while the other solution has a separation greater than $\thetae$ (wide solution). From the Bayesian analysis with the constraints of the event time scale and angular Einstein radius together with the location of the source lying in the far disk behind the bulge, we find that the planet is a super-Jupiter with a mass of $4.4^{+3.0}_{-2.0}~M_{\rm J}$ and the stellar binary components are early and late M-type dwarfs with masses $0.42^{+0.29}_{-0.19}~M_\odot$ and $0.12^{+0.08}_{-0.05}~M_\odot$, respectively, and the planetary system is located at a distance of $D_{\rm L}=7.6^{+1.2}_{-0.9}~{\rm kpc}$. The planet is a circumstellar planet according to the wide solution, while it is a circumbinary planet according to the close solution. The projected primary-planet separation is $2.8^{+3.2}_{-2.5}~{\rm au}$ commonly for the close and wide solutions, but the primary-secondary binary separation of the close solution, $0.75^{+0.87}_{-0.66}~{\rm au}$, is widely different from the separation, $10.5^{+12.1}_{-9.2}~{\rm au}$, of the wide solution.

The widespread presence of ribonucleic acid (RNA) catalysts and cofactors in Earth's biosphere today suggests that RNA was the first biopolymer to support Darwinian evolution. However, most "path-hypotheses" to generate building blocks for RNA require reduced nitrogen-containing compounds not made in useful amounts in the CO2-N2-H2O atmospheres of the Hadean. We review models for Earth's impact history that invoke a single ~10^23 kg impactor (Moneta) to account for measured amounts of platinum, gold, and other siderophilic ("iron-loving") elements on the Earth and Moon. If it were the last sterilizing impactor, Moneta would have reduced the atmosphere but not its mantle, opening a "window of opportunity" for RNA synthesis, a period when RNA precursors rained from the atmosphere to land holding oxidized minerals that stabilize advanced RNA precursors and RNA. Surprisingly, this combination of physics, geology, and chemistry suggests a time when RNA formation was most probable, ~120 +/- 100 million years after Moneta's impact, or ~4.36 +/- 0.1 billion years ago. Uncertainties in this time are driven by uncertainties in rates of productive atmosphere loss and amounts of sub-aerial land.

We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3\sigma$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago.

Scintillating NaI(Tl) crystals are used for various rare decay experiments, such as dark matter searches. The hygroscopicity of NaI(Tl) crystal makes the construction of crystal detectors in these experiments challenging and requires a tight encapsulation to prevent from air contact. More importantly, in a low radioactivity measurement, identification of external radiations and surface contamination is crucial to characterize the origin of total crystal radioactivities. Studies for NaI(Tl) crystal encapsulation with active organic scintillator vetoes have been performed to mitigate the above-mentioned issues simultaneously. A bare crystal is directly coupled with liquid and plastic scintillators to tag external radiations that penetrate from the outer part of the crystal. We report the pulse shape discrimination for organic scintillator pulses from those of the crystal scintillator in a single detector setup which makes external gammas identifiable and long-term stability tests of the detector setup.

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.

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

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

We present the analysis of the event OGLE-2017-BLG-1186 from the 2017 Spitzer microlensing campaign. This is a remarkable microlensing event because its source is photometrically bright and variable, which makes it possible to perform an asteroseismic analysis using ground-based data. We find that the source star is an oscillating red giant with average timescale of $\sim 9$ d. The asteroseismic analysis also provides us source properties including the source angular size ($\sim 27~\mu{\rm as}$) and distance ($\sim 11.5$ kpc), which are essential for inferring the properties of the lens. When fitting the light curve, we test the feasibility of Gaussian Processes (GPs) in handling the correlated noise caused by the variable source. We find that the parameters from the GP model are generally more loosely constrained than those from the traditional $\chi^2$ minimization method. We note that this event is the first microlensing system for which asteroseismology and GPs have been used to account for the variable source. With both finite-source effect and microlens parallax measured, we find that the lens is likely a $\sim 0.045~M_{\odot}$ brown dwarf at distance $\sim 9.0$ kpc, or a $\sim 0.073~M_{\odot}$ ultracool dwarf at distance $\sim 9.8$ kpc. Combining the estimated lens properties with a Bayesian analysis using a Galactic model, we find a $\sim 35$ per cent probability for the lens to be a bulge object and $\sim 65$ per cent to be a background disc object.

As part of the National Science Foundation funded "Gemini in the Era of MultiMessenger Astronomy" (GEMMA) program, Gemini Observatory is developing GNAO, a widefield adaptive optics (AO) facility for Gemini-North on Maunakea, the only 8m-class open-access telescope available to the US astronomers in the northern hemisphere. GNAO will provide the user community with a queue-operated Multi-Conjugate AO (MCAO) system, enabling a wide range of innovative solar system, Galactic, and extragalactic science with a particular focus on synergies with JWST in the area of time-domain astronomy. The GNAO effort builds on institutional investment and experience with the more limited block-scheduled Gemini Multi-Conjugate System (GeMS), commissioned at Gemini South in 2013. The project involves close partnerships with the community through the recently established Gemini AO Working Group and the GNAO Science Team, as well as external instrument teams. The modular design of GNAO will enable a planned upgrade to a Ground Layer AO (GLAO) mode when combined with an Adaptive Secondary Mirror (ASM). By enhancing the natural seeing by an expected factor of two, GLAO will vastly improve Gemini North's observing efficiency for seeing-limited instruments and strengthen its survey capabilities for multi-messenger astronomy.

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

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

We report the discovery of the microlensing planet OGLE-2018-BLG-0740Lb. The planet is detected with a very strong signal of $\Delta\chi^2\sim 4630$, but the interpretation of the signal suffers from two types of degeneracies. One type is caused by the previously known close/wide degeneracy, and the other is caused by an ambiguity between two solutions, in which one solution requires to incorporate finite-source effects, while the other solution is consistent with a point-source interpretation. Although difficult to be firmly resolved based on only the photometric data, the degeneracy is resolved in strong favor of the point-source solution with the additional external information obtained from astrometric and spectroscopic observations. The small astrometric offset between the source and baseline object supports that the blend is the lens and this interpretation is further secured by the consistency of the spectroscopic distance estimate of the blend with the lensing parameters of the point-source solution. The estimated mass of the host is $1.0\pm 0.1~M_\odot$ and the mass of the planet is $4.5\pm 0.6~M_{\rm J}$ (close solution) or $4.8\pm 0.6~M_{\rm J}$ (wide solution) and the lens is located at a distance of $3.2\pm 0.5$~kpc. The bright nature of the lens, with $I\sim 17.1$ ($V\sim 18.2$), combined with its dominance of the observed flux suggest that radial-velocity (RV) follow-up observations of the lens can be done using high-resolution spectrometers mounted on large telescopes, e.g., VLT/ESPRESSO, and this can potentially not only measure the period and eccentricity of the planet but also probe for close-in planets. We estimate that the expected RV amplitude would be $\sim 60\sin i ~{\rm m~s}^{-1}$.

We report a multiplanetary system found from the analysis of microlensing event OGLE-2018-BLG-1011, for which the light curve exhibits a double-bump anomaly around the peak. We find that the anomaly cannot be fully explained by the binary-lens or binary-source interpretations and its description requires the introduction of an additional lens component. The 3L1S (3 lens components and a single source) modeling yields three sets of solutions, in which one set of solutions indicates that the lens is a planetary system in a binary, while the other two sets imply that the lens is a multiplanetary system. By investigating the fits of the individual models to the detailed light curve structure, we find that the multiple-planet solution with planet-to-host mass ratios $\sim 9.5\times 10^{-3}$ and $\sim 15\times 10^{-3}$ are favored over the other solutions. From the Bayesian analysis, we find that the lens is composed of two planets with masses $1.8^{+3..4}_{-1.1}~M_{\rm J}$ and $2.8^{+5.1}_{-1.7}~M_{\rm J}$ around a host with a mass $0.18^{+0.33}_{-0.10}~M_\odot$ and located at a distance $7.1^{+1.1}_{-1.5}~{\rm kpc}$. The estimated distance indicates that the lens is the farthest system among the known multiplanetary systems. The projected planet-host separations are $a_{\perp,2}=1.8^{+2.1}_{-1.5}~{\rm au}$ ($0.8^{+0.9}_{-0.6}~{\rm au}$) and $a_{\perp,3}=0.8^{+0.9}_{-0.6}~{\rm au}$, where the values of $a_{\perp,2}$ in and out the parenthesis are the separations corresponding to the two degenerate solutions, indicating that both planets are located beyond the snow line of the host, as with the other four multiplanetary systems previously found by microlensing.

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

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

We present the analysis of the microlensing event OGLE-2015-BLG-1670, detected in a high-extinction field, very close to the Galactic plane. Due to the dust extinction along the line of sight, this event was too faint to be detected before it reached the peak of magnification. The microlensing light-curve models indicate a high-magnification event with a maximum of $A_\mathrm{max}\gtrsim200$, very sensitive to planetary deviations. An anomaly in the light curve has been densely observed by the microlensing surveys MOA, KMTNet, and OGLE. From the light-curve modeling, we find a planetary anomaly characterized by a planet-to-host mass ratio, $q=\left(1.00^{+0.18}_{-0.16}\right)\times 10^{-4}$, at the peak recently identified in the mass-ratio function of microlensing planets. Thus, this event is interesting to include in future statistical studies about planet demography. We have explored the possible degeneracies and find two competing planetary models resulting from the $s\leftrightarrow1/s$ degeneracy. However, because the projected separation is very close to $s=1$, the physical implications for the planet for the two solutions are quite similar, except for the value of $s$. By combining the light-curve parameters with a Galactic model, we have estimated the planet mass $M_2=17.9^{+9.6}_{-8.8}\,\mathrm{M}_\oplus$ and the lens distance $D_\mathrm{L}=6.7^{+1.0}_{-1.3}\,\mathrm{kpc}$, corresponding to a Neptune-mass planet close to the Galactic bulge. Such events with a low absolute latitude ($|b|\approx 1.1\,\mathrm{deg}$) are subject to both high extinction and more uncertain source distances, two factors that may affect the mass measurements in the provisional Wide Field Infrared Survey Telescope fields. More events are needed to investigate the potential trade-off between the higher lensing rate and the difficulty in measuring masses in these low-latitude fields.

We propose a mechanism of elementary thermal dark matter with mass up to $10^{14}$ GeV, within a standard cosmological history, whose relic abundance is determined solely by its interactions with the Standard Model, without violating the perturbative unitarity bound. The dark matter consists of many nearly degenerate particles which scatter with the Standard Model bath in a nearest-neighbor chain, and maintain chemical equilibrium with the Standard Model bath by in-equilibrium decays and inverse decays. The phenomenology includes super heavy elementary dark matter and heavy relics that decay at various epochs in the cosmological history, with implications for CMB, structure formation and cosmic ray experiments.

The morphological properties of the outflowing circumstellar envelopes surrounding mass-losing stars in eccentric binary systems are presented from a set of three-dimensional hydrodynamical model simulations. Based on four template models of the envelope viewed for a range of inclination angles of the systems, we implement interpretative tools for observations at high spectral/angular resolutions (as illustrated via velocity channel maps as well as position-velocity, radius-velocity, and angle-radius diagrams). Within this framework, the image and kinematical structures can be used to place constraints on the orbital parameters of the system. Specifically, three unique characteristic patterns in the envelopes are found that distinguish these systems from those in binary systems in circular orbits. Bifurcation of the spiral pattern, asymmetry in the interarm density depression, and a concurrent spiral/ring appearance all point to a binary system with an eccentric orbit. The methodology presented in this paper is illustrated in an analysis of recent radio observations of several asymptotic giant branch stars.

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

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

We show that ultra-light scalar dark matter (fuzzy dark matter) in galaxies has a quantum mechanical typical acceleration scale about $10^{-10}\,\mbox{ms}^{-2}$, which leads to the baryonic Tully-Fisher relation. Baryonic matter at central parts of galaxies acts as a boundary condition for dark matter wave equation and influences stellar rotation velocities in halos. Without any modification of gravity or mechanics this model also explains the radial acceleration relation and MOND-like behavior of gravitational acceleration found in galaxies having flat rotation curves. This analysis can be extended to the Faber-Jackson relation.

We report the discovery of a $Spitzer$ microlensing planet OGLE-2018-BLG-0596Lb, with preferred planet-host mass ratio $q \sim 2\times10^{-4}$. The planetary signal, which is characterized by a short $(\sim 1~{\rm day})$ "bump" on the rising side of the lensing light curve, was densely covered by ground-based surveys. We find that the signal can be explained by a bright source that fully envelops the planetary caustic, i.e., a "Hollywood" geometry. Combined with the source proper motion measured from $Gaia$, the $Spitzer$ satellite parallax measurement makes it possible to precisely constrain the lens physical parameters. The preferred solution, in which the planet perturbs the minor image due to lensing by the host, yields a Uranus-mass planet with a mass of $M_{\rm p} = 13.9\pm1.6~M_{\oplus}$ orbiting a mid M-dwarf with a mass of $M_{\rm h} = 0.23\pm0.03~M_{\odot}$. There is also a second possible solution that is substantially disfavored but cannot be ruled out, for which the planet perturbs the major image. The latter solution yields $M_{\rm p} = 1.2\pm0.2~M_{\oplus}$ and $M_{\rm h} = 0.15\pm0.02~M_{\odot}$. By combining the microlensing and $Gaia$ data together with a Galactic model, we find in either case that the lens lies on the near side of the Galactic bulge at a distance $D_{\rm L} \sim 6\pm1~{\rm kpc}$. Future adaptive optics observations may decisively resolve the major image/minor image degeneracy.

We report the discovery and characterization of KMT-2018-BLG-1990Lb, a Jovian planet $(m_p=0.57_{-0.25}^{+0.79}\,M_J)$ orbiting a late M dwarf $(M=0.14_{-0.06}^{+0.20}\,M_\odot)$, at a distance $(D_L=1.23_{-0.43}^{+1.06}\,\kpc)$, and projected at $2.6\pm 0.6$ times the snow line distance, i.e., $a_{\rm snow}\equiv 2.7\,\au (M/M_\odot)$, This is the second Jovian planet discovered by KMTNet in its low cadence ($0.4\,{\rm hr}^{-1}$) fields, demonstrating that this population will be well characterized based on survey-only microlensing data.

We report the discovery of the planet OGLE-2018-BLG-0532Lb, with very obvious signatures in the light curve that lead to an estimate of the planet-host mass ratio $q=M_{\rm planet}/M_{\rm host}\simeq 1\times10^{-4}$. Although there are no obvious systematic residuals to this double-lens/single-source (2L1S) fit, we find that $\chi^2$ can be significantly improved by adding either a third lens (3L1S, $\Delta\chi^2=81$) or second source (2L2S, $\Delta\chi^2=65$) to the lens-source geometry. After thorough investigation, we conclude that we cannot decisively distinguish between these two scenarios and therefore focus on the robustly-detected planet. However, given the possible presence of a second planet, we investigate to what degree and with what probability such additional planets may affect seemingly single-planet light curves. Our best estimates for the properties of the lens star and the secure planet are: a host mass $M\sim 0.25\,M_\odot$, system distance $D_L\sim 1\,$kpc and planet mass $m_{p,1}= 8\,M_\oplus$ with projected separation $a_{1,\perp}=1.4\,$au. However, there is a relatively bright $I=18.6$ (and also relatively blue) star projected within $<50\,$mas of the lens, and if future high-resolution images show that this is coincident with the lens, then it is possible that it is the lens, in which case, the lens would be both more massive and more distant than the best-estimated values above.

We report the discovery of KMT-2018-BLG-1292Lb, a super-Jovian $M_{\rm planet} = 4.5\pm 1.3\,M_J$ planet orbiting an F or G dwarf $M_{\rm host} = 1.5\pm 0.4\,M_\odot$, which lies physically within ${\cal O}(10\,\pc)$ of the Galactic plane. The source star is a heavily extincted $A_I\sim 5.2$ luminous giant that has the lowest Galactic latitude, $b=-0.28^\circ$, of any planetary microlensing event. The relatively blue blended light is almost certainly either the host or its binary companion, with the first explanation being substantially more likely. This blend dominates the light at $I$ band and completely dominates at $R$ and $V$ bands. Hence, the lens system can be probed by follow-up observations immediately, i.e., long before the lens system and the source separate due to their relative proper motion. The system is well characterized despite the low cadence $\Gamma=0.15$--$0.20\,{\rm hr^{-1}}$ of observations and short viewing windows near the end of the bulge season. This suggests that optical microlensing planet searches can be extended to the Galactic plane at relatively modest cost.

We report the mass and distance measurements of two single-lens events from the 2017 Spitzer microlensing campaign. The ground-based observations yield the detection of finite-source effects, and the microlens parallaxes are derived from the joint analysis of ground-based observations and Spitzer observations. We find that the lens of OGLE-2017-BLG-1254 is a $0.60 \pm 0.03 M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.11~\text{kpc}$, where $D_{\rm LS}$ is the distance between the lens and the source. The second event, OGLE-2017-BLG-1161, is subject to the known satellite parallax degeneracy, and thus is either a $0.51^{+0.12}_{-0.10} M_{\odot}$ star with $D_{\rm LS} = 0.40 \pm 0.12~\text{kpc}$ or a $0.38^{+0.13}_{-0.12} M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.19~\text{kpc}$. Both of the lenses are therefore isolated stars in the Galactic bulge. By comparing the mass and distance distributions of the eight published Spitzer finite-source events with the expectations from a Galactic model, we find that the Spitzer sample is in agreement with the probability of finite-source effects occurrence in single lens events.

In the present work, we would like to elaborate on Shakura Sunyaev Alpha disk model. To be more concrete, inclusion of the effect of the radiation pressure on the accretion phenomena have been carefully studied. To summarize, the present modification of Shakura Sunyaev Alpha disk model of disk accretion can be regarded as a new candidate theory of ADAF. Reduction of accreting matter's angular velocity leads to the suppression of the disk luminosity and the surface temperature which essentially implies the transition of the standard accretion disk model from CDAF to ADAF.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

In this work, we present the analysis of the binary microlensing event OGLE-2018-BLG-0022 that is detected toward the Galactic bulge field. The dense and continuous coverage with the high-quality photometry data from ground-based observations combined with the space-based {\it Spitzer} observations of this long time-scale event enables us to uniquely determine the masses $M_1=0.40 \pm 0.05~M_\odot$ and $M_2=0.13\pm 0.01~M_\odot$ of the individual lens components. Because the lens-source relative parallax and the vector lens-source relative proper motion are unambiguously determined, we can likewise unambiguously predict the astrometric offset between the light centroid of the magnified images (as observed by the {\it Gaia} satellite) and the true position of the source. This prediction can be tested when the individual-epoch {\it Gaia} astrometric measurements are released.

We present the detection of RR Lyrae variable stars in the field of the Sagittarius II (Sgr II) ultra-faint dwarf (UFD) galaxy. Using B, V time-series photometry obtained with the Korea Microlensing Telescope Network (KMTNet) 1.6 m telescope at CTIO and G-band data from Gaia Data Release 2 (DR2), we identified and characterized two ab-type and four c-type RR Lyrae variables. Five out of the six stars are clustered within three half-light radii (~4.8') of the galaxy indicating their association with Sgr II, while the RRab star V4 is located ~22' from the galaxy center. By excluding V4, the high c-type fraction (0.8) and the period of the only RRab star V3 (P_V3 = 0.666 days) suggest an Oosterhoff II (Oo II) classification for Sgr II. Located close to the locus of Oo II clusters in the period-amplitude diagram, V3 is similar to RRab stars in other UFDs having Oosterhoff-intermediate and Oo II properties. Sgr II is, however, more compact than usual UFDs, placed in between star clusters and dwarf galaxies in the size-luminosity plane, and therefore spectroscopic studies are eventually required to ascertain the true nature of this stellar system. We derive the metallicity ([Fe/H]_RRab = -2.1 +- 0.3) and heliocentric distance (~64 +- 3 kpc) of Sgr II from the RR Lyrae stars, and estimate its age (~12 Gyr) based on our stellar population models. The Oosterhoff properties of UFDs can be explained with the evolution effect of RR Lyrae stars in the instability strip.

This paper reports a design study for a space-based decametric wavelength telescope. While not a new concept, this design study focused on many of the operational aspects that would be required for an actual mission. This design optimized the number of spacecraft to insure good visibility of approx. 80% of the radio galaxies -- the primary science target for the mission. A 5,000 km lunar orbit was selected to guarantee minimal gravitational perturbations from Earth and lower radio interference. Optimal schemes for data downlink, spacecraft ranging, and power consumption were identified. An optimal mission duration of 1 year was chosen based on science goals, payload complexity, and other factors. Finally, preliminary simulations showing image reconstruction were conducted to confirm viability of the mission. This work is intended to show the viability and science benefits of conducting multi-spacecraft networked radio astronomy missions in the next few years.

We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the waveforms sensitive to the shower muon content. A total number of two photon candidates found in the search is fully compatible with the expected background. The $95\%\,$CL limits on the diffuse flux of the photons with energies greater than $10^{18.0}$, $10^{18.5}$, $10^{19.0}$, $10^{19.5}$ and $10^{20.0}$ eV are set at the level of $0.067$, $0.012$, $0.0036$, $0.0013$, $0.0013~\mbox{km}^{-2}\mbox{yr}^{-1}\mbox{sr}^{-1}$ correspondingly.

The aim of this work is to create a new catalog of reliable AGN candidates selected from the AKARI NEP-Deep field. Selection of the AGN candidates was done by applying a fuzzy SVM algorithm, which allows to incorporate measurement uncertainties into the classification process. The training dataset was based on the spectroscopic data available for selected objects in the NEP-Deep and NEP-Wide fields. The generalization sample was based on the AKARI NEP-Deep field data including objects without optical counterparts and making use of the infrared information only. A high quality catalog of previously unclassified 275 AGN candidates was prepared.

We present near-infrared spectroscopic analysis of 12 red giant members of the Galactic open cluster NGC 6940. High-resolution (R$\simeq$45000) and high signal-to-noise ratio (S/N > 100) near-infrared H and K band spectra were gathered with the Immersion Grating Infrared Spectrograph (IGRINS) on the 2.7m Smith Telescope at McDonald Observatory. We obtained abundances of H-burning (C, N, O), ${\alpha}$ (Mg, Si, S, Ca), light odd-Z (Na, Al, P, K), Fe-group (Sc, Ti, Cr, Fe, Co, Ni) and neutron-capture (Ce, Nd, Yb) elements. We report the abundances of S, P, K, Ce, and Yb in NGC 6940 for the first time. Many OH and CN features in the H band were used to obtain O and N abundances. C abundances were measured from four different features: CO molecular lines in the K band, high excitation C I lines present in both near-infrared and optical, CH and $C_2$ bands in the optical region. We have also determined $^{12}C/^{13}C$ ratios from the R-branch band heads of first overtone (2-0) and (3-1) $^{12}CO$ (2-0) $^{13}CO$ lines near 23440 $\overset{\lower.5em\circ}{\mathrm{A}}$ and (3-1) $^{13}CO$ lines at about 23730 $\overset{\lower.5em\circ}{\mathrm{A}}$. We have also investigated the HF feature at 23358.3 $\overset{\lower.5em\circ}{\mathrm{A}}$, finding solar fluorine abundances without ruling out a slight enhancement. For some elements (such as the ${\alpha}$ group), IGRINS data yield more internally self-consistent abundances. We also revisited the CMD of NGC 6940 by determining the most probable cluster members using Gaia DR2. Finally, we applied Victoria isochrones and MESA models in order to refine our estimates of the evolutionary stages of our targets.

In the circumnuclear region of the nearby radio galaxy NGC~1052. Our $\lambda$3mm VLBI observations with the Korean VLBI Network have spatially resolved the broad HCO$^{+}$ absorption at a velocity range of 1350--1850 km~s$^{-1}$ against a double-sided nuclear jet, and have revealed that the HCO$^{+}$ absorption is concentrated on the receding jet and the nuclear components. The distribution of the HCO$^{+}$ absorbing gas strongly supports the circumnuclear torus surrounding the supermassive black hole. From the estimations of the column density and the volume density of molecular hydrogen, the size of the molecular gas region in the torus is at least 1~pc. The broad spectral profile of HCO$^{+}$ is likely to be a blend of multiple gas clumps with various velocities. The HCO$^{+}$ absorption of NGC~1052 could trace complex kinematics in the vicinity of the active galactic nucleus, such as inflow, outflow, turbulence, and so on.

We report two microlensing events, KMT-2017-BLG-1038 and KMT-2017-BLG-1146 that are caused by planetary systems. These events were discovered by KMTNet survey observations from the $2017$ bulge season. The discovered systems consist of a planet and host star with mass ratios, $5.3_{-0.4}^{+0.2} \times 10^{-3}$ and $2.0_{-0.1}^{+0.6} \times 10^{-3}$, respectively. Based on a Bayesian analysis assuming a Galactic model without stellar remnant hosts, we find that the planet, KMT-2017-BLG-1038Lb, is a super Jupiter-mass planet ($M_{\rm p}= 2.04_{-1.15}^{+2.02}\,M_{\rm J}$) orbiting a mid-M dwarf host ($M_{\rm h}= 0.37_{-0.20}^{+0.36}\, M_{\odot}$) that is located at $6.01_{-1.72}^{+1.27}$ kpc toward the Galactic bulge. The other planet, KMT-2017-BLG-1146Lb, is a sub Jupiter-mass planet ($M_{\rm p}= 0.71_{-0.42}^{+0.80}\,M_{\rm J}$) orbiting a mid-M dwarf host ($M_{\rm h}= 0.33_{-0.20}^{+0.36}\,M_{\odot}$) at a distance toward the Galactic bulge of $6.50_{-2.00}^{+1.38}$ kpc. Both are potentially gaseous planets that are beyond their hosts' snow lines. These typical microlensing planets will be routinely discovered by second-generation microlensing surveys, rapidly increasing the number of detections.

We analyze the binary microlensing event MOA-2016-BLG-231, which was observed from the ground and from Spitzer. The lens is composed of very low-mass brown dwarfs (BDs) with $M_1 = 21^{+12}_{-5} \ M_J$ and $M_2 = 9^{+5}_{-2}\ M_J$, and it is located in the Galactic disk $D_{\rm L} = 2.85^{+0.88}_{-0.50}\ {\rm kpc}$. This is the fifth binary brown dwarf discovered by microlensing, and the BD binary is moving counter to the orbital motion of disk stars. Constraints on the lens physical properties come from late time, non-caustic-crossing features of the Spitzer light curve. Thus, MOA-2016-BLG-231 shows how Spitzer plays a crucial role in resolving the nature of BDs in binary BD events with short timescale ($\lesssim 10$ days).

Much of the star formation is obscured by dust. For the complete understanding of the cosmic star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP) field (5.4 deg$^2$). This took significant amount of satellite's lifetime, $\sim$10\% of the entire pointed observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry, we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical imaging ($\sim$25.9ABmag) in a small area of 1.0 deg$^2$. As a result, there remained thousands of AKARI's infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru telescope, we obtained deep enough optical images of the entire AKARI NEP field in 5 broad bands ($g\sim$27.5mag). These provided photometric redshift, and thereby IR luminosity for the previously undetected faint AKARI IR sources. Combined with the accurate mid-IR luminosity measurement, we constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. We have measured restframe 8$\mu$m, 12$\mu$m luminosity functions (LFs), and estimated total infrared LFs at 0.35$<$z$<$2.2. Our results are consistent with our previous work, but with much reduced statistical errors thanks to the large area coverage of the new data. We have possibly witnessed the turnover of CSFH at $z\sim$2.

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses.

Intensive Monitoring Survey of Nearby Galaxies (IMSNG) is a high cadence observation program monitoring nearby galaxies with high probabilities of hosting supernovae (SNe). IMSNG aims to constrain the SN explosion mechanism by inferring sizes of SN progenitor systems through the detection of the shock-heated emission that lasts less than a few days after the SN explosion. To catch the signal, IMSNG utilizes a network of 0.5-m to 1-m class telescopes around the world and monitors the images of 60 nearby galaxies at distances D < 50 Mpc to a cadence as short as a few hours. The target galaxies are bright in near-ultraviolet (NUV) with M_NUV < -18.4 AB mag and have high probabilities of hosting SNe (0.06 SN/yr per galaxy). With this strategy, we expect to detect the early light curves of 3.4 SNe per year to a depth of R ~ 19.5 mag, enabling us to detect the shock-heated emission from a progenitor star with a radius as small as 0.1 R_sun. The accumulated data will be also useful for studying faint features around the target galaxies and other science projects. So far, 18 SNe have occurred in our target fields (16 in IMSNG galaxies) over 5 years, confirming our SN rate estimate of 0.06 SN/yr per galaxy.

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

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

The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events induced by the primary protons and iron. An average atomic mass of UHECR is presented for energies $10^{18.0}-10^{20.0}\ \mbox{eV}$. The average atomic mass of primary particles shows no significant energy dependence and corresponds to $\langle \ln A \rangle = 2.0 \pm 0.1 (stat.) \pm 0.44 (syst.)$. The result is compared to the mass composition obtained by the Telescope Array with $\mbox{X}_{\mbox{max}}$ technique along with the results of other experiments. Possible systematic errors of the method are discussed.

Planet formation theories predict the existence of free-floating planets that have been ejected from their parent systems. Although they emit little or no light, they can be detected during gravitational microlensing events. Microlensing events caused by rogue planets are characterized by very short timescales $t_{\rm E}$ (typically below two days) and small angular Einstein radii $\theta_{\rm E}$ (up to several uas). Here we present the discovery and characterization of two ultra-short microlensing events identified in data from the Optical Gravitational Lensing Experiment (OGLE) survey, which may have been caused by free-floating or wide-orbit planets. OGLE-2012-BLG-1323 is one of the shortest events discovered thus far ($t_{\rm E}$=0.155 +/- 0.005 d, $\theta_{\rm E}$=2.37 +/- 0.10 uas) and was caused by an Earth-mass object in the Galactic disk or a Neptune-mass planet in the Galactic bulge. OGLE-2017-BLG-0560 ($t_{\rm E}$=0.905 +/- 0.005 d, $\theta_{\rm E}$=38.7 +/- 1.6 uas) was caused by a Jupiter-mass planet in the Galactic disk or a brown dwarf in the bulge. We rule out stellar companions up to a distance of 6.0 and 3.9 au, respectively. We suggest that the lensing objects, whether located on very wide orbits or free-floating, may originate from the same physical mechanism. Although the sample of ultrashort microlensing events is small, these detections are consistent with low-mass wide-orbit or unbound planets being more common than stars in the Milky Way.

We analyze the gravitational binary-lensing event OGLE-2016-BLG-0156, for which the lensing light curve displays pronounced deviations induced by microlens-parallax effects. The light curve exhibits 3 distinctive widely-separated peaks and we find that the multiple-peak feature provides a very tight constraint on the microlens-parallax effect, enabling us to precisely measure the microlens parallax $\pi_{\rm E}$. All the peaks are densely and continuously covered from high-cadence survey observations using globally located telescopes and the analysis of the peaks leads to the precise measurement of the angular Einstein radius $\theta_{\rm E}$. From the combination of the measured $\pi_{\rm E}$ and $\theta_{\rm E}$, we determine the physical parameters of the lens. It is found that the lens is a binary composed of two M dwarfs with masses $M_1=0.18\pm 0.01\ M_\odot$ and $M_2=0.16\pm 0.01\ M_\odot$ located at a distance $D_{\rm L}= 1.35\pm 0.09\ {\rm kpc}$. According to the estimated lens mass and distance, the flux from the lens comprises an important fraction, $\sim 25\%$, of the blended flux. The bright nature of the lens combined with the high relative lens-source motion, $\mu=6.94\pm 0.50\ {\rm mas}\ {\rm yr}^{-1}$, suggests that the lens can be directly observed from future high-resolution follow-up observations.

The kinematics of isolated brown dwarfs in the Galaxy, beyond the solar neighborhood, is virtually unknown. Microlensing has the potential to probe this hidden population, as it can measure both the mass and five of the six phase-space coordinates (all except the radial velocity) even of a dark isolated lens. However, the measurements of both the microlens parallax and finite-source effects are needed in order to recover the full information. Here, we combine $Spitzer$ satellite parallax measurement with the ground-based light curve, which exhibits strong finite-source effects, of event OGLE-2017-BLG-0896. We find two degenerate solutions for the lens (due to the known satellite-parallax degeneracy), which are consistent with each other except for their proper motion. The lens is an isolated brown dwarf with a mass of either $18\pm1M_J$ or $20\pm1M_J$. This is the lowest isolated-object mass measurement to date, only $\sim$45\% more massive than the theoretical deuterium-fusion boundary at solar metallicity, which is the common definition of a free-floating planet. The brown dwarf is located at either $3.9\pm0.1$ kpc or $4.1\pm0.1$ kpc toward the Galactic bulge, but with proper motion in the opposite direction of disk stars, with one solution suggesting it is moving within the Galactic plane. While it is possibly a halo brown dwarf, it might also represent a different, unknown population.

The latest Lyman-$\alpha$ forest data severely constrain the conventional warm dark matter solution to small-scale issues in the cold dark matter paradigm. It has been also reported that unconstrained astrophysical processes may address the issues. In response to this situation, we revisit the decaying dark matter solution to the issues, discussing possible signatures to discriminate decaying dark matter from astrophysical processes as a solution to small-scale issues. We consider an axinolike particle (ALPino) decaying into an axionlike particle (ALP) and gravitino with the lifetime around the age of the Universe. The ALPino mass is sub-PeV and slightly ($\Delta m/m\sim 10^{-4}$) larger than the gravitino mass, and thus the dark matter abundance does not alter virtually after the ALPino decays. On the other hand, the gravitino produced from the ALPino decay obtains a kick velocity of $\sim 30 \,{\rm km / s}$, which is sufficiently larger than a circular velocity of dwarf galaxies to impact their dark matter distributions. The Lyman-$\alpha$ forest constraints are relieved since only a small fraction ($\sim10$%) of dark matter experiences the decay at that time. Decaying dark matter is thus promoted to a viable solution to small-scale issues. The ALPino relic abundance is determined predominantly by the decay of the lightest ordinary supersymmetric particle. The monochromatic ALP emission from the ALPino decay is converted to $\sim 50 \,{\rm GeV}$ photon under the Galactic magnetic field. The morphology of the gamma-ray flux shows a distinctive feature of the model when compared to decaying dark matter that directly decays into photons. Once detected, such distinctive signals discriminate the decaying dark matter solution to small-scale issues from unconstrained astrophysical processes.

We report a possible discovery of a relic star cluster in the Sextans dwarf spheroidal galaxy. Using the \textit{hk} index ($\equiv$($Ca-b$)$-$($b-y$)) as a photometric metallicity indicator, we discriminate the metal-poor and metal-rich stars in the galaxy and find unexpected number density excess of metal-poor stars located 7'.7 ($\sim$190 pc in projected distance) away from the known galactic center. The $V-I$ color$-$magnitude diagram (CMD) for stars around the density excess reveals that both the main sequence and the giant branch are considerably narrower and redder than the bulk of field stars in Sextans. Our stellar population models show (a) that the narrow CMD is best reproduced by a simple stellar population with an age of $\sim$13 Gyr and [Fe/H] of $\sim$$-$2.3 dex, and (b) that the redder $V-I$ color of the $hk$-weak population is explained $only$ if it is $\sim$2 Gyr older than the field stars. The results lead us to conclude that the off-centered density peak is likely associated with an old, metal-poor globular cluster. The larger spatial extent ($>$ 80 pc in radius) and the smaller number of stars ($\sim$1000) than typical globular clusters point to a star cluster that is in the process of dissolution. The finding serves as the first detection of a surviving star cluster in Sextans, supporting previous suggestions of the presence of star cluster remnants in the galaxy. If confirmed, the survival of a relic star cluster until now implies a $cored$ dark matter halo profile for this dwarf galaxy.

We report the discovery of a low mass-ratio planet $(q = 1.3\times10^{-4})$, i.e., 2.5 times higher than the Neptune/Sun ratio. The planetary system was discovered from the analysis of the KMT-2017-BLG-0165 microlensing event, which has an obvious short-term deviation from the underlying light curve produced by the host of the planet. Although the fit improvement with the microlens parallax effect is relatively low, one component of the parallax vector is strongly constrained from the light curve, making it possible to narrow down the uncertainties of the lens physical properties. A Bayesian analysis yields that the planet has a super-Neptune mass $(M_{2}=34_{-12}^{+15}~M_{\oplus})$ orbiting a Sun-like star $(M_{1}=0.76_{-0.27}^{+0.34}~M_{\odot})$ located at $4.5~{\rm kpc}$. The blended light is consistent with these host properties. The projected planet-host separation is $a_{\bot}={3.45_{-0.95}^{+0.98}}~{\rm AU}$, implying that the planet is located outside the snowline of the host, i.e., $a_{sl}\sim2.1~{\rm AU}$. KMT-2017-BLG-0165Lb is the sixteenth microlensing planet with mass ratio $q<3\times10^{-4}$. Using the fifteen of these planets with unambiguous mass-ratio measurements, we apply a likelihood analysis to investigate the form of the mass-ratio function in this regime. If we adopt a broken power law for the form of this function, then the break is at $q_{\rm br}\simeq0.55\times10^{-4}$, which is much lower than previously estimated. Moreover, the change of the power law slope, $\zeta>3.3$ is quite severe. Alternatively, the distribution is also suggestive of a "pile-up" of planets at Neptune-like mass ratios, below which there is a dramatic drop in frequency.

We present a study correlating the spatial locations of young star clusters with those of molecular clouds in NGC~5194, in order to investigate the timescale over which clusters separate from their birth clouds. The star cluster catalogues are from the Legacy ExtraGalactic UV Survey (LEGUS) and the molecular clouds from the Plateau de Bure Interefrometer Arcsecond Whirpool Survey (PAWS). We find that younger star clusters are spatially closer to molecular clouds than older star clusters. The median ages for clusters associated with clouds is 4~Myr whereas it is 50~Myr for clusters that are sufficiently separated from a molecular cloud to be considered unassociated. After $\sim$6~Myr, the majority of the star clusters lose association with their molecular gas. Younger star clusters are also preferentially located in stellar spiral arms where they are hierarchically distributed in kpc-size regions for 50-100~Myr before dispersing. The youngest star clusters are more strongly clustered, yielding a two-point correlation function with $\alpha=-0.28\pm0.04$, than the GMCs ($\alpha=-0.09\pm0.03$) within the same PAWS field. However, the clustering strength of the most massive GMCs, supposedly the progenitors of the young clusters for a star formation efficiency of a few percent, is comparable ($\alpha=-0.35\pm0.05$) to that of the clusters. We find a galactocentric-dependence for the coherence of star formation, in which clusters located in the inner region of the galaxy reside in smaller star-forming complexes and display more homogeneous distributions than clusters further from the centre. This result suggests a correlation between the survival of a cluster complex and its environment.

We show that microlensing event KMT-2016-BLG-1107 displays a new type of degeneracy between wide-binary and close-binary Hollywood events in which a giant-star source envelops the planetary caustic. The planetary anomaly takes the form of a smooth, two-day "bump" far out on the falling wing of the light curve, which can be interpreted either as the source completely enveloping a minor-image caustic due to a close companion with mass ratio $q=0.036$, or partially enveloping a major-image caustic due to a wide companion with $q=0.004$. The best estimates of the companion masses are both in the planetary regime ($3.3^{+3.5}_{-1.8}\,M_{\rm jup}$ and $0.090^{+0.096}_{-0.037}\,M_{\rm jup}$) but differ by an even larger factor than the mass ratios due to different inferred host masses. We show that the two solutions can be distinguished by high-resolution imaging at first light on next-generation ("30m") telescopes. We provide analytic guidance to understand the conditions under which this new type of degeneracy can appear.

We report the discovery of a giant planet in the KMT-2016-BLG-1397 microlensing event, which was found by The Korea Microlensing Telescope Network (KMTNet) alone. The time scale of this event is t_E = 40.0 +- 0.5 days and the mass ratio between the lens star and its companion is q = 0.016 +- 0.002. The planetary perturbation in the light curve is a smooth bump, resulting in the classical binary-lens/binary-source (2L1S/1L2S) degeneracy. We measure the V - I color of the (putative) two sources in the 1L2S model, and then effectively rule out the binary source solution. The finite-source effect is marginally detected. Combined with the limits on the blend flux and the probability distribution of the source size normalized by the Einstein radius rho, a Bayesian analysis yields the lens mass M_L = 0.45+0.33-0.28 M_sun, at distance of D_L = 6.60+1.10-1.30 kpc. Thus the companion is a super-Jupiter of a mass m_p = 7.0+5.2-4.3 M_J , at a projected separation r = 5.1+1.5-1.7 AU, indicating that the planet is well beyond the snow line of the host star.

The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did not optimize the parameter values but kept them fixed to the best-fit values found by Auger, namely 9.7% for the anisotropic fraction of cosmic rays assumed to originate from the SBGs in the list and 12.9{\deg} for the angular scale of the correlations. The energy threshold we adopted is 43 EeV, corresponding to 39 EeV in Auger when taking into account the energy-scale difference between two experiments. We find that the TA data is compatible with isotropy to within 1.1{\sigma} and with the Auger result to within 1.4{\sigma}, meaning that it is not capable to discriminate between these two hypotheses.

We present the analysis of the caustic-crossing binary microlensing event OGLE-2017-BLG-0039. Thanks to the very long duration of the event, with an event time scale $t_{\rm E}\sim 130$ days, the microlens parallax is precisely measured despite its small value of $\pie\sim 0.06$. The analysis of the well-resolved caustic crossings during both the source star's entrance and exit of the caustic yields the angular Einstein radius $\thetae\sim 0.6$~mas. The measured $\pie$ and $\thetae$ indicate that the lens is a binary composed of two stars with masses $\sim 1.0~M_\odot$ and $\sim 0.15~M_\odot$, and it is located at a distance of $\sim 6$ kpc. From the color and brightness of the lens estimated from the determined lens mass and distance, it is expected that $\sim 2/3$ of the $I$-band blended flux comes from the lens. Therefore, the event is a rare case of a bright lens event for which high-resolution follow-up observations can confirm the nature of the lens.

We present the analysis of the planetary microlensing event MOA-2016-BLG-319. The event light curve is characterized by a brief ($\sim 3$ days) anomaly near the peak produced by minor-image perturbations. From modeling, we find two distinct solutions that describe the observed light curve almost equally well. From the investigation of the lens-system configurations, we find that the confusion in the lensing solution is caused by the degeneracy between the two solutions resulting from the source passages on different sides of the planetary caustic. These degeneracies can be severe for major-image perturbations but it is known that they are considerably less severe for minor-image perturbations. From the comparison of the lens-system configuration with those of two previously discovered planetary events, for which similar degeneracies were reported, we find that the degeneracies are caused by the special source trajectories that passed the star-planet axes at approximately right angles. By conducting a Bayesian analysis, it is estimated that the lens is a planetary system in which a giant planet with a mass $M_{\rm p}=0.62^{+1.16}_{-0.33}~M_{\rm J}$ ($0.65^{+1.21}_{-0.35}~M_{\rm J}$) is orbiting a low-mass M-dwarf host with a mass $M_{\rm h}=0.15^{+0.28}_{-0.08}~M_\odot$. Here the planet masses in and out of the parentheses represent the masses for the individual degenerate solutions. The projected host-planet separations are $a_\perp\sim 0.95$ au and $\sim 1.05$ au for the two solutions. The identified degeneracy indicates the need to check similar degeneracies in future analyses of planetary lensing events with minor-image perturbations.

We present a long-exposure (~10 hr) image of the supernova (SN) remnant Cassiopeia A (Cas A) obtained with the UKIRT 3.8-m telescope using a narrow band filter centered at 1.644 um emission. The passband contains [Fe II] 1.644 um and [Si I] 1.645 um lines, and our `deep [Fe II]+[Si I] image' provides an unprecedented panoramic view of Cas A, showing both shocked and unshocked SN ejecta together with shocked circumstellar medium at subarcsec (~0.7 arcsec or 0.012 pc) resolution. The diffuse emission from the unshocked SN ejecta has a form of clumps, filaments, and arcs, and their spatial distribution correlates well with that of the Spitzer [Si II] infrared emission, suggesting that the emission is likely due to [Si I] line not [Fe II] line as in shocked material. The structure of the optically-invisible western area of Cas A is clearly seen for the first time. The area is filled with many Quasi-Stationary Flocculi (QSFs) and fragments of the disrupted ejecta shell. We suggest that the anomalous radio properties in this area could be due to the increased number of such dense clumps. We identified 309 knots in the deep [Fe II]+[Si I] image and classified them into QSFs and fast-moving knots (FMKs). The total H+He mass of QSFs is ~0.23 Msun, implying that the mass fraction of dense clumps in the progenitor's red-supergiant wind is 4--13%. The spatial distribution of QSFs suggests that there had been a highly asymmetric mass loss $10^4$--$10^5$ yr before the SN explosion. The mass of the [Fe II] line-emitting, shocked dense Fe ejecta is ~3x$10^{-5}$ Msun. The comparison with the ionic S-line dominated Hubble Space Telescope WFC3/IR image suggests that the outermost FMKs in the southeastern area are Fe-rich.

We present the analyses of two short-timescale $(t_{\rm E} \sim 5~{\rm days})$ microlensing events, KMT-2016-BLG-1820 and KMT-2016-BLG-2142. In both light curves, the brief anomalies were clearly captured and densely covered by the Korea Microlensing Telescope Network survey. From these analyses, we find that both events have small Einstein radii of $\theta_{\rm E} = 0.12~{\rm mas}$, suggesting that the binary-lens systems are composed of very low-mass components and/or are located much closer to the lensed stars than to Earth. From Bayesian analyses, we find that these binaries have total system masses of $0.043_{-0.018}^{+0.043}~M_{\odot}$ and $0.088_{-0.041}^{+0.120}~M_{\odot}$, implying that they are well within the very-low-mass regime. The estimated lens-component masses indicate that the binary lenses consist of a giant-planet/brown-dwarf pair (KMT-2016-BLG-1820), and a dark/faint object pair (KMT-2016-BLG-2140) that are located near the deuterium-burning and hydrogen-burning mass limits, respectively. Both lens systems are likely to be in the Galactic disk with estimated distances of about $6$ kpc and $7$ kpc. The projected lens-components separations are $1.1$ AU and $0.8$ AU, and the mass ratios are $0.11$ and $0.20$. These prove that the microlensing method is effective to identify these closely-separated very-low-mass binaries having low mass-ratios.

We present a library of high-resolution (R $\equiv$ $\lambda$/$\Delta$$\lambda$ $\sim$ 45,000) and high signal-to-noise ratio (S/N $\geq$ 200) near-infrared spectra for stars of a wide range of spectral types and luminosity classes. The spectra were obtained with the Immersion GRating INfrared Spectrograph (IGRINS) covering the full range of the H (1.496-1.780 $\mu$m) and K (2.080-2.460 $\mu$m) atmospheric windows. The targets were primarily selected for being MK standard stars covering a wide range of effective temperatures and surface gravities with metallicities close to the Solar value. Currently, the library includes flux-calibrated and telluric-absorption-corrected spectra of 84 stars, with prospects for expansion to provide denser coverage of the parametric space. Throughout the H and K atmospheric windows, we identified spectral lines that are sensitive to $T_\mathrm{eff}$ or $\log g$ and defined corresponding spectral indices. We also provide their equivalent widths. For those indices, we derive empirical relations between the measured equivalent widths and the stellar atmospheric parameters. Therefore, the derived empirical equations can be used to calculate $T_\mathrm{eff}$ and $\log g$ of a star without requiring stellar atmospheric models.

We have derived elemental abundances of three field red horizontal branch stars using high-resolution (R$\simeq$ 45,000), high signal-to-noise ratio (S/N $\gtrsim$ 200) $H$ and $K$ band spectra obtained with the Immersion Grating Infrared Spectrograph (IGRINS). We have determined the abundances of 21 elements including $\alpha$ (Mg, Si, Ca, S), odd-Z (Na, Al, P, K), Fe-group (Sc, Ti, Cr, Co, Ni), neutron-capture (Ce, Nd, Yb), and CNO group elements. S, P and K are determined for the first time in these stars. $H$ and $K$ band spectra provide a substantial number of S I lines, which potentially can lead to a more robust exploration of the role of sulfur in the cosmochemical evolution of the Galaxy. We have also derived $^{12}$C/$^{13}$C ratios from synthetic spectra of the first overtone (2$-$0) and (3$-$1) $^{12}$CO and (2$-$0) $^{13}$CO lines near 23440 \AA\ and $^{13}$CO (3$-$1) lines at about 23730 \AA. Comparison of our results with the ones obtained from the optical region suggests that the IGRINS high-resolution $H$ and $K$ band spectra offer more internally self-consistent atomic lines of the same species for several elements, especially the $\alpha$ elements. This in turn provides more reliable abundances for the elements with analytical difficulties in the optical spectral range.

We present the analysis of a binary microlensing event KMT-2016-BLG-2052, for which the lensing-induced brightening of the source star lasted for 2 seasons. We determine the lens mass from the combined measurements of the microlens parallax $\pie$ and angular Einstein radius $\thetae$. The measured mass indicates that the lens is a binary composed of M dwarfs with masses of $M_1\sim 0.34~M_\odot$ and $M_2\sim 0.17~M_\odot$. The measured relative lens-source proper motion of $\mu\sim 3.9~{\rm mas}~{\rm yr}^{-1}$ is smaller than $\sim 5~{\rm mas}~{\rm yr}^{-1}$ of typical Galactic lensing events, while the estimated angular Einstein radius of $\thetae\sim 1.2~{\rm mas}$ is substantially greater than the typical value of $\sim 0.5~{\rm mas}$. Therefore, it turns out that the long time scale of the event is caused by the combination of the slow $\mu$ and large $\thetae$ rather than the heavy mass of the lens. From the simulation of Galactic lensing events with very long time scales ($t_{\rm E}\gtrsim 100$ days), we find that the probabilities that long time-scale events are produced by lenses with masses $\geq 1.0~M_\odot$ and $\geq 3.0~M_\odot$ are $\sim 19\%$ and 2.6\%, respectively, indicating that events produced by heavy lenses comprise a minor fraction of long time-scale events. The results indicate that it is essential to determine lens masses by measuring both $\pie$ and $\thetae$ in order to firmly identify heavy stellar remnants such as neutron stars and black holes.

We have explored the role environmental factors play in determining characteristics of young stellar objects in nearby dwarf irregular and Blue Compact Dwarf galaxies. Star clusters are characterized by concentrations, masses, and formation rates, OB associations by mass and mass surface density, O stars by their numbers and near-ultraviolet absolute magnitudes, and HII regions by Halpha surface brightnesses. These characteristics are compared to surrounding galactic pressure, stellar mass density, HI surface density, and star formation rate surface density. We find no trend of cluster characteristics with environmental properties, implying that larger scale effects are more important in determining cluster characteristics or that rapid dynamical evolution erases memory of the initial conditions. On the other hand, the most massive OB associations are found at higher pressure and HI surface density, and there is a trend of higher HII region Halpha surface brightness with higher pressure, suggesting that a higher concentration of massive stars and gas are found preferentially in regions of higher pressure. At low pressures we find massive stars but not bound clusters and OB associations. We do not find evidence for an increase of cluster formation efficiency as a function of star formation rate density. However, there is an increase in the ratio of the number of clusters to number of O stars with pressure, perhaps reflecting an increase in clustering properties with star formation rate.

We report the discovery and the analysis of the short timescale binary-lens microlensing event, MOA-2015-BLG-337. The lens system could be a planetary system with a very low mass host, around the brown dwarf/planetary mass boundary, or a brown dwarf binary. We found two competing models that explain the observed light curves with companion/host mass ratios of q~0.01 and ~0.17, respectively. From the measurement of finite source effects in the best-fit planetary model, we find a relatively small angular Einstein radius of theta_E ~ 0.03 mas which favors a low mass lens. We conduct a Bayesian analysis to obtain the probability distribution of the lens properties. The results for the planetary models strongly depend on the minimum mass, M_min, in the assumed mass function. In summary, there are two solutions of the lens system: (1) a brown dwarf/planetary mass boundary object orbited by a super-Neptune (the planetary model with M_min=0.001 M_sun) and (2) a brown dwarf binary (the binary model). If the planetary models is correct, this system can be one of a new class of planetary system, having a low host mass and also a planetary mass ratio (q <0.03) between the host and its companion. The discovery of the event is important for the study of planetary formation in very low mass objects. In addition, it is important to consider all viable solutions in these kinds of ambiguous events in order for the future comprehensive statistical analyses of planetary/binary microlensing events.

The orientation of galaxy spin vectors within the large scale structure has been considered an important test of our understanding of structure formation. We investigate the angular changes of galaxy spin vectors in clusters - denser environments than are normally focused upon, using hydrodynamic zoomed simulations of 17 clusters YZiCS and a set of complementary controlled simulations. The magnitude by which galaxies change their spin vector is found to be a function of their rotational support with larger cumulative angular changes of spin vectors when they have initially lower $V_{\theta}/\sigma$. We find that both mergers and tidal perturbations can significantly swing spin vectors, with larger changes in spin vector for smaller pericentre distances. Strong tidal perturbations are also correlated with the changes in stellar mass and specific angular momentum of satellite galaxies. However, changes in spin vector can often result in a canceling out of previous changes. As a result, the integrated angular change is always much larger than the angular change measured at any instant. Also, overall the majority of satellite galaxies do not undergo mergers or sufficiently strong tidal perturbation after infall into clusters, and thus they end up suffering little change to their spin vectors. Taken as a whole, these results suggest that any signatures of spin alignment from the large scale structure will be preserved in the cluster environment for many gigayears.

Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes.

We propose a mechanism generating primordial magnetic fields after the $e^+e^-$ annihilations. Our mechanism involves an ultra-light axion-like particle (ALP) which constitutes the dark matter, and a dark $U(1)_X$ gauge boson introduced to bypass the obstacle placed by the conductivity of cosmic plasma. In our scheme, a coherently oscillating ALP amplifies the dark photon field, and part of the amplified dark photon field is concurrently converted to the ordinary magnetic field through the ALP-induced magnetic mixing. For the relevant ALP mass range $10^{-21} {\rm eV}\lesssim m_\phi\lesssim 10^{-17}{\rm eV}$, our mechanism can generate $B\sim 10^{-24} \,{\rm G} \,(m_\phi/10^{-17} {\rm eV})^{5/4}$ with a coherent length $\lambda \sim (m_\phi/10^{-17} {\rm eV})^{-1/2}$ kpc, which is large enough to provide a seed of the galactic magnetic fields. The mechanism also predicts a dark $U(1)_X$ electromagnetic field $E_X \sim B_X\sim 80\,{\rm nG}\, (m_\phi/10^{-17}{\rm eV})^{-1/4}$, which can result in interesting astrophysical/cosmological phenomena by inducing the mixings between the ALP, ordinary photon, and dark photon states.

It was recently pointed out that semi-annihilating dark matter (DM) may experience a novel temperature evolution dubbed as self-heating. Exothermic semi-annihilation converts the DM mass to the kinetic energy. This yields a unique DM temperature evolution, $T_{\chi} \propto 1 / a$, in contrast to $ T_{\chi} \propto 1 / a^{2}$ for free-streaming non-relativistic particles. Self-heating continues as long as self-scattering sufficiently redistributes the energy of DM particles. In this paper, we study the evolution of cosmological perturbations in self-heating DM. We find that sub-GeV self-heating DM leaves a cutoff on the subgalactic scale of the matter power spectrum when the self-scattering cross section is $\sigma_{\rm self} / m_{\chi} \sim {\cal O} (1) \,{\rm cm}^{2} /{\rm g}$. Then we present a particle physics realization of the self-heating DM scenario. The model is based on recently proposed strongly interacting massive particles with pion-like particles in a QCD-like sector. Pion-like particles semi-annihilate into an axion-like particle, which is thermalized with dark radiation. The dark radiation temperature is smaller than the standard model temperature, evading the constraint from the effective number of neutrino degrees of freedom. It is easily realized when the dark sector is populated from the standard model sector through a small coupling.

It was recently pointed out that semi-annihilating dark matter (DM) may experience a novel temperature evolution dubbed as self-heating. Exothermic semi-annihilation converts the DM mass to the kinetic energy. This yields a unique DM temperature evolution, $T_{\chi} \propto 1 / a$, in contrast to $ T_{\chi} \propto 1 / a^{2}$ for free-streaming non-relativistic particles. Self-heating continues as long as self-scattering sufficiently redistributes the energy of DM particles. In this paper, we study the evolution of cosmological perturbations in self-heating DM. We find that sub-GeV self-heating DM leaves a cutoff on the subgalactic scale of the matter power spectrum when the self-scattering cross section is $\sigma_{\rm self} / m_{\chi} \sim {\cal O} (1) \,{\rm cm}^{2} /{\rm g}$. Then we present a particle physics realization of the self-heating DM scenario. The model is based on recently proposed strongly interacting massive particles with pion-like particles in a QCD-like sector. Pion-like particles semi-annihilate into an axion-like particle, which is thermalized with dark radiation. The dark radiation temperature is smaller than the standard model temperature, evading the constraint from the effective number of neutrino degrees of freedom. It is easily realized when the dark sector is populated from the standard model sector through a small coupling.

We present a parsec-scale molecular hydrogen (H$_2$ 1-0 S(1) at 2.12~\micron) outflow discovered from the UKIRT Widefield Infrared Survey for H$_2$. The outflow is located in the infrared dark cloud core MSXDC G053.11+00.05 MM1 at 1.7 kpc and likely associated with two young stellar objects (YSOs) at the center. The overall morphology of the outflow is bipolar along the NE-SW direction with a brighter lobe to the southwest, but the detailed structure consists of several flows and knots. With the total length of $\sim$1 pc, the outflow luminosity is fairly high with $L_{\rm H_{2}} > 6~L_{\sun}$, implying a massive outflow-driving YSO if the entire outflow is driven by a single source. The two putative driving sources, located at the outflow center, show photometric variability of $\gtrsim$1 mag in {\it H}- and {\it K}-bands. This, with their early evolutionary stage from spectral energy distribution (SED) fitting, indicates that both are capable of ejecting outflows and may be eruptive variable YSOs. The YSO masses inferred from SED fitting are $\sim$10~$M_{\sun}$ and $\sim$5~$M_{\sun}$, suggesting the association of the outflow with massive YSOs. The geometrical morphology of the outflow is well explained by the lower mass YSO by assuming a single source origin, but without kinematic information, the contribution from the higher mass YSO cannot be ruled out. Considering star formation process by fragmentation of a high-mass core into several lower mass stars, we also suggest the possible presence of another, yet-undetected driving source deeply embedded in the core.

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

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

We present the analysis of the binary-lens microlensing event OGLE-2017-BLG-0537. The light curve of the event exhibits two strong caustic-crossing spikes among which the second caustic crossing was resolved by high-cadence surveys. It is found that the lens components with a mass ratio $\sim 0.5$ are separated in projection by $\sim 1.3\thetae$, where $\thetae$ is the angular Einstein radius. Analysis of the caustic-crossing part yields $\thetae=1.77\pm 0.16$~mas and a lens-source relative proper motion of $\mu =12.4\pm 1.1~{\rm mas}~{\rm yr}^{-1}$. The measured $\mu$ is the third highest value among the events with measured proper motions and $\sim 3$ times higher than the value of typical Galactic bulge events, making the event a strong candidate for follow-up observations to directly image the lens by separating it from the source. From the angular Einstein radius combined with the microlens parallax, it is estimated that the lens is composed of two main-sequence stars with masses $M_1\sim 0.4~M_\odot$ and $M_2\sim 0.2~M_\odot$ located at a distance of $D_{\rm L}\sim 1.2$~kpc. However, the physical lens parameters are not very secure due to the weak microlens-parallax signal, and thus we cross check the parameters by conducting a Bayesian analysis based on the measured Einstein radius and event timescale combined with the blending constraint. From this, we find that the physical parameters estimated from the Bayesian analysis are consistent with those based on the measured microlens parallax. Resolving the lens from the source can be done in about 5 years from high-resolution follow-up observations and this will provide a rare opportunity to test and refine the microlensing model.

We describe a new microlensing-event alert algorithm that is tailored to the Korea Microlensing Telescope Network (KMTNet) multi-observatory system. The algorithm focuses on detecting "rising" events, i.e., events whose brightness is increasing as a function of time. The algorithm proceeds in three steps. It first identifies light curves with at least $N_{\rm high}$ points that are at least $3\sigma$ above the median and that had been taken since shortly before the previous search for new events. It then demands that there are at least $N_{\rm high}$ $\textit{consecutive}$ high points when considering any combination of light curves from one, two, or three observatories. Finally, it fits to a "rising microlensing model" consisting of a broken line, i.e., flat before some time $t_{\rm rise}$ and rising linearly afterward. If this fit is better than a flat line by $\Delta\chi^2>\Delta\chi^2_{\rm thresh}$, the light curve is sent for human review. Here, $(N_{\rm high},\Delta\chi^2_{\rm thresh})=(5,250)$ or $(10,400)$, depending on the field cadence. For 2018, KMTNet alerts will initially be restricted to a few northern bulge fields and may gradually extend to the full northern bulge. Further expansion of coverage is expected in 2019.

Microlensing is a powerful and unique technique to probe isolated objects in the Galaxy. To study the characteristics of these interesting objects based on the microlensing method, measurement of the microlens parallax is required to determine the properties of the lens. Of the various methods to measure microlens parallax, the most robust way is to make simultaneous ground- and space-based observations, i.e., by measuring the space-based microlens parallax. However, space-based campaigns usually require "expensive" resources. Gould & Yee (2012) proposed an idea called the "cheap space-based microlens parallax" that can measure the lens-parallax using only $2$ or $3$ space-based observations of high-magnification events. This cost-effective observation strategy to measure microlens parallaxes could be used by space-borne telescopes to build a complete sample for studying isolated objects. This would enable a direct measurement of the mass function including both extremely low-mass objects and high-mass stellar remnants. However, to adopt this idea requires a test to check how it would work in actual situations. Thus, we present the first practical test of the idea using the high-magnification microlensing event OGLE-2016-BLG-1045, for which a subset of Spitzer observations fortuitously duplicate the prescription of Gould & Yee (2012). From the test, we confirm that the measurement of the lens-parallax adopting this idea has sufficient accuracy to determine the physical properties of the isolated lens.

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

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

Using the kinematic position angles (PA_kin), an accurate indicator for the spin axis of a galaxy, obtained from the ATLAS3D integral-field-unit (IFU) spectroscopic data, we discovered that 57 Virgo early-type galaxies tend to prefer the specific PA_kin values of 20 degree and 100 degree, suggesting that they are kinematically aligned with each other. These kinematic alignment angles are further associated with the directions of the two distinct axes of the Virgo cluster extending east-west and north-south, strongly suggesting that the two distinct axes are the filamentary structures within the cluster as a trace of infall patterns of galaxies. Given that the spin axis of a massive early-type galaxy does not change easily even in clusters from the hydrodynamic simulations, Virgo early-type galaxies are likely to fall into the cluster along the filamentary structures while maintaining their angular momentum. This implies that many early-type galaxies in clusters are formed in filaments via major mergers before subsequently falling into the cluster. Investigating the kinematic alignment in other clusters will allow us to understand the formation of galaxy clusters and early-type galaxies.

Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

We report the detection of RR Lyrae variable stars in Crater II, a recently discovered large and diffuse satellite dwarf galaxy of the Milky Way (MW). Based on B, V time-series photometry obtained with the Korea Microlensing Telescope Network (KMTNet) 1.6-m telescope at CTIO, we identified 83 ab-type and 13 c-type pulsators by fitting template light curves. The detected RR Lyrae stars are centrally concentrated, which ensures that most of them are members of Crater II. In terms of the distribution of RRab stars in the period-amplitude diagram, Crater II is clearly different from ultra-faint dwarf (UFD) galaxies, but very similar to the two classical MW dwarf spheroidal (dSph) galaxies Draco and Carina with Oosterhoff-intermediate (Oo-int) properties. Combined with the mean period of ab-type variables (<Pab> = 0.631+-0.004 d) and the c-type fraction (~0.14) in Crater II, this suggests an Oo-int classification for Crater II and implies that its nature is more like a dSph rather than a UFD. We also estimated the mean metallicity, reddening, and distance of Crater II, from the photometric and pulsation properties of the RR Lyrae stars. The stellar population model we have constructed indicates that Crater II is dominated by an old population, but is relatively younger than the oldest globular clusters in the MW. With a lack of high-amplitude short-period RRab stars, Crater II, like most of the other less massive dSphs, is probably not a surviving counterpart of the major building blocks of the MW halo.

One of the main theories for explaining the formation of spiral arms in galaxies is the stationary density wave theory. This theory predicts the existence of an age gradient across the arms. We use the stellar cluster catalogues of the galaxies NGC 1566, M51a, and NGC 628 from the Legacy Extragalactic UV Survey (LEGUS) program. In order to test for the possible existence of an age sequence across the spiral arms, we quantified the azimuthal offset between star clusters of different ages in our target galaxies. We found that NGC 1566, a grand-design spiral galaxy with bisymmetric arms and a strong bar, shows a significant age gradient across the spiral arms that appears to be consistent with the prediction of the stationary density wave theory. In contrast, M51a with its two well-defined spiral arms and a weaker bar does not show an age gradient across the arms. In addition, a comparison with non LEGUS star cluster catalogues for M51a yields similar results. We believe that the spiral structure of M51a is not the result of a stationary density wave with a fixed pattern speed. Instead, tidal interactions could be the dominant mechanism for the formation of spiral arms. We also found no offset in the azimuthal distribution of star clusters with different ages across the weak spiral arms of NGC 628.

The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

Reconstructing the expansion history of the Universe from type Ia supernovae data, we fit the growth rate measurements and put model-independent constraints on some key cosmological parameters, namely, $\Omega_\mathrm{m},\gamma$, and $\sigma_8$. The constraints are consistent with those from the concordance model within the framework of general relativity, but the current quality of the data is not sufficient to rule out modified gravity models. Adding the condition that dark energy density should be positive at all redshifts, independently of its equation of state, further constrains the parameters and interestingly supports the concordance model.

We show that dense OGLE and KMTNet $I$-band survey data require four bodies (sources plus lenses) to explain the microlensing light curve of OGLE-2015-BLG-1459. However, these can equally well consist of three lenses and one source (3L1S), two lenses and two sources (2L2S) or one lens and three sources (1L3S). In the 3L1S and 2L2S interpretations, the host is a brown dwarf and the dominant companion is a Neptune-class planet, with the third body (in the 3L1S case) being a Mars-class object that could have been a moon of the planet. In the 1L3S solution, the light curve anomalies are explained by a tight (five stellar radii) low-luminosity binary source that is offset from the principal source of the event by $\sim 0.17\,\au$. These degeneracies are resolved in favor of the 1L3S solution by color effects derived from comparison to MOA data, which are taken in a slightly different ($R/I$) passband. To enable current and future ($WFIRST$) surveys to routinely characterize exomoons and distinguish among such exotic systems requires an observing strategy that includes both a cadence faster than 9 min$^{-1}$ and observations in a second band on a similar timescale.

We analyze the combined Spitzer and ground-based data for OGLE-2017-BLG-1140 and show that the event was generated by a Jupiter-class $(m_p\simeq 1.6\,M_{\rm jup})$ planet orbiting a mid-late M dwarf $(M\simeq 0.2\,M_\odot)$ that lies $D_{LS}\simeq 1.0\,\mathrm{kpc}$ in the foreground of the microlensed, Galactic-bar, source star. The planet-host projected separation is $a_\perp \simeq 1.0\,\mathrm{au}$, i.e., well-beyond the snow line. By measuring the source proper motion ${\mathbf{\mu}}_s$ from ongoing, long-term OGLE imaging, and combining this with the lens-source relative proper motion ${\mathbf{\mu}}_\mathrm{rel}$ derived from the microlensing solution, we show that the lens proper motion ${\mathbf{\mu}}_l={\mathbf{\mu}}_\mathrm{rel} + {\mathbf{\mu}}_s$ is consistent with the lens lying in the Galactic disk, although a bulge lens is not ruled out. We show that while the Spitzer and ground-based data are comparably well fitted by planetary (i.e., binary-lens, 2L1S) models and by binary-source (1L2S) models, the combination of Spitzer and ground-based data decisively favor the planetary model. This is a new channel to resolve the 2L1S/1L2S degeneracy, which can be difficult to break in some cases.

One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high $\it{muon \, purity}$. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of $\sim$65% at a lateral distance of 2000 m from the shower core is $1.72 \pm 0.10{\rm (stat.)} \pm 0.37 {\rm (syst.)}$ times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a $1.67 \pm 0.10 \pm 0.36$ excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances.

We present Korea Microlensing Telescope Network (KMTNet) light curves for microlensing-event candidates for the 2016 season, which covers an area of 97 sq.deg observed at cadences ranging from Gamma=0.2/hr to Gamma=8/hr from three southern sites in Chile, South Africa, and Australia. These 2163 light curves are comprised of 1856 "clear microlensing" and 307 "possible microlensing" events (including 265 previously released from the K2 C9 field). The data policy is very similar to the one governing the 2015 release. The changes relative to 2015 in the algorithms to find and vet microlensing events are comprehensively described.

We report the discovery of microlensing planet OGLE-2017-BLG-0373Lb. We show that while the planet-host system has an unambiguous microlens topology, there are two geometries within this topology that fit the data equally well, which leads to a factor 2.5 difference in planet-host mass ratio, i.e., $q=1.5\times 10^{-3}$ vs. $q=0.6\times 10^{-3}$. We show that this is an "accidental degeneracy" in the sense that it is due to a gap in the data. We dub it "the caustic-chirality degeneracy". We trace the mathematical origins of this degeneracy, which should enable similar degenerate solutions to be easily located in the future. A Bayesian estimate, based on a Galactic model, yields a host mass $M=0.25^{+0.30}_{-0.15} M_\odot$ at a distance $D_L=5.9^{+1.3}_{-1.95}$ kpc. The lens-source relative proper motion is relatively fast, $\mu=9$ mas/yr, which implies that the host mass and distance can be determined by high-resolution imaging after about 10 years. The same observations could in principle resolve the discrete degeneracy in $q$, but this will be more challenging.

We present a high-resolution, near-IR spectroscopic study of multiple outflows in the LkH$\alpha$ 234 star formation region using the Immersion GRating INfrared Spectrometer (IGRINS). Spectral mapping over the blueshifted emission of HH 167 allowed us to distinguish at least three separate, spatially overlapped, outflows in H${_2}$ and [Fe II] emission. We show that the H${_2}$ emission represents not a single jet, but complex multiple outflows driven by three known embedded sources: MM1, VLA 2, and VLA 3. There is a redshifted H${_2}$ outflow at a low velocity, $\VLSR$ $<$ $+$50 {\kms}, with respect to the systemic velocity of $\VLSR$ $=$ $-$11.5 {\kms}, that coincides with the H${_2}$O masers seen in earlier radio observations two arcseconds southwest of VLA 2. We found that the previously detected [Fe II] jet with $|$$\VLSR$$|$ $>$ 100 {\kms} driven by VLA 3B is also detected in H${_2}$ emission, and confirm that this jet has a position angle about 240$\degree$. Spectra of the redshifted knots at 14$\arcsec$$-$65$\arcsec$ northeast of LkH$\alpha$ 234 are presented for the first time. These spectra also provide clues to the existence of multiple outflows. We detected high-velocity (50$-$120 {\kms}) H${_2}$ gas in the multiple outflows around LkH$\alpha$ 234. Since these gases move at speeds well over the dissociation velocity ($>$ 40 {\kms}), the emission must originate from the jet itself rather than H${_2}$ gas in the ambient medium. Also, position-velocity diagrams and excitation diagram indicate that emission from knot C in HH 167 come from two different phenomena, shocks and photodissociation.

The freeze-out of dark matter (DM) depends on the evolution of the DM temperature. The DM temperature does not have to follow the standard model one, when the elastic scattering is not sufficient to maintain the kinetic equilibrium. We study the temperature evolution of the semi-annihilating DM, where a pair of the DM particles annihilate into one DM particle and another particle coupled to the standard model sector. We find that the kinetic equilibrium is maintained solely via semi-annihilation until the last stage of the freeze-out. After the freeze-out, semi-annihilation converts the mass deficit to the kinetic energy of DM, which leads to non-trivial evolution of the DM temperature. We argue that the DM temperature redshifts like radiation as long as the DM self-interaction is efficient. We dub this novel temperature evolution as self-heating. Notably, the structure formation is suppressed at subgalactic scales like keV-scale warm DM but with GeV-scale self-heating DM if the self-heating lasts roughly until the matter-radiation equality. The long duration of the self-heating requires the large self-scattering cross section, which in turn flattens the DM density profile in inner halos. Consequently, self-heating DM can be a unified solution to apparent failures of cold DM to reproduce the observed subgalactic scale structure of the Universe.

The freeze-out of dark matter (DM) depends on the evolution of the DM temperature. The DM temperature does not have to follow the standard model one, when the elastic scattering is not sufficient to maintain the kinetic equilibrium. We study the temperature evolution of the semi-annihilating DM, where a pair of the DM particles annihilate into one DM particle and another particle coupled to the standard model sector. We find that the kinetic equilibrium is maintained solely via semi-annihilation until the last stage of the freeze-out. After the freeze-out, semi-annihilation converts the mass deficit to the kinetic energy of DM, which leads to non-trivial evolution of the DM temperature. We argue that the DM temperature redshifts like radiation as long as the DM self-interaction is efficient. We dub this novel temperature evolution as self-heating. Notably, the structure formation is suppressed at subgalactic scales like keV-scale warm DM but with GeV-scale self-heating DM if the self-heating lasts roughly until the matter-radiation equality. The long duration of the self-heating requires the large self-scattering cross section, which in turn flattens the DM density profile in inner halos. Consequently, self-heating DM can be a unified solution to apparent failures of cold DM to reproduce the observed subgalactic scale structure of the Universe.

We report the discovery of a planetary system in which a super-earth orbits a late M-dwarf host. The planetary system was found from the analysis of the microlensing event OGLE-2017-BLG-0482, wherein the planet signal appears as a short-term anomaly to the smooth lensing light curve produced by the host. Despite its weak signal and short duration, the planetary signal was firmly detected from the dense and continuous coverage by three microlensing surveys. We find a planet/host mass ratio of $q\sim 1.4\times 10^{-4}$. We measure the microlens parallax $\pi_{\rm E}$ from the long-term deviation in the observed lensing light curve, but the angular Einstein radius $\theta_{\rm E}$ cannot be measured because the source trajectory did not cross the planet-induced caustic. Using the measured event timescale and the microlens parallax, we find that the masses of the planet and the host are $M_{\rm p}=9.0_{-4.5}^{+9.0}\ M_\oplus$ and $M_{\rm host}=0.20_{-0.10}^{+0.20}\ M_\odot$, respectively, and the projected separation between them is $a_\perp=1.8_{-0.7}^{+0.6}$ au. The estimated distance to the lens is $D_{\rm L}=5.8_{-2.1}^{+1.8}$ kpc. The discovery of the planetary system demonstrates that microlensing provides an important method to detect low-mass planets orbiting low-mass stars.

We present a study of the dust-to-gas ratios in five nearby galaxies NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with archival HST/ACS data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques, and correlate these maps with neutral HI and CO gas maps from literature, including The HI Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey (HERACLES). We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$ could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool (BEAST) for NGC 7793 and find systematically lower extinctions from SED-fitting as compared to isochrone matching.

We report the discovery of a giant planet in the OGLE-2017-BLG-1522 microlensing event. The planetary perturbations were clearly identified by high-cadence survey experiments despite the relatively short event timescale of $t_{\rm E} \sim 7.5$ days. The Einstein radius is unusually small, $\theta_{\rm E} = 0.065\,$mas, implying that the lens system either has very low mass or lies much closer to the microlensed source than the Sun, or both. A Bayesian analysis yields component masses $(M_{\rm host}, M_{\rm planet})=(46_{-25}^{+79}, 0.75_{-0.40}^{+1.26})~M_{\rm J}$ and source-lens distance $D_{\rm LS} = 0.99_{-0.54}^{+0.91}~{\rm kpc}$, implying that this is a brown-dwarf/Jupiter system that probably lies in the Galactic bulge, a location that is also consistent with the relatively low lens-source relative proper motion $\mu = 3.2 \pm 0.5~{\rm mas}~{\rm yr^{-1}}$. The projected companion-host separation is $0.59_{-0.11}^{+0.12}~{\rm AU}$, indicating that the planet is placed beyond the snow line of the host, i.e., $a_{sl} \sim 0.12~{\rm AU}$. Planet formation scenarios combined with the small companion-host mass ratio $q \sim 0.016$ and separation suggest that the companion could be the first discovery of a giant planet that formed in a protoplanetary disk around a brown dwarf host.

We present a general framework where the effective field theory of single field inflation is extended by the inclusion of operators with mass dimension 3 and 4 in the unitary gauge. These higher dimensional operators introduce quartic and sextic corrections to the dispersion relation. We study the regime of validity of this extended effective field theory of inflation and the effect of these higher dimensional operators on CMB observables associated with scalar perturbations, such as the speed of sound, the amplitude of the power spectrum and the tensor-to-scalar ratio. Tensor perturbations remain instead, unaltered.

We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources.

Mass measurements of gravitational microlenses require one to determine the microlens parallax $\pie$, but precise $\pie$ measurement, in many cases, is hampered due to the subtlety of the microlens-parallax signal combined with the difficulty of distinguishing the signal from those induced by other higher-order effects. In this work, we present the analysis of the binary-lens event OGLE-2017-BLG-0329, for which $\pie$ is measured with a dramatically improved precision using additional data from space-based $Spitzer$ observations. We find that while the parallax model based on the ground-based data cannot be distinguished from a zero-$\pie$ model at 2$\sigma$ level, the addition of the $Spitzer$ data enables us to identify 2 classes of solutions, each composed of a pair of solutions according to the well-known ecliptic degeneracy. It is found that the space-based data reduce the measurement uncertainties of the north and east components of the microlens-parallax vector $\pivec_{\rm E}$ by factors $\sim 18$ and $\sim 4$, respectively. With the measured microlens parallax combined with the angular Einstein radius measured from the resolved caustic crossings, we find that the lens is composed of a binary with components masses of either $(M_1,M_2)\sim (1.1,0.8)\ M_\odot$ or $\sim (0.4,0.3)\ M_\odot$ according to the two solution classes. The first solution is significantly favored but the second cannot be securely ruled out based on the microlensing data alone. However, the degeneracy can be resolved from adaptive optics observations taken $\sim 10$ years after the event.

We present the analysis of KMT-2016-BLG-0212, a low flux-variation $(I_{\rm flux-var}\sim 20$) microlensing event, which is well-covered by high-cadence data from the three Korea Microlensing Telescope Network (KMTNet) telescopes. The event shows a short anomaly that is incompletely covered due to the brief visibility intervals that characterize the early microlensing season when the anomaly occurred. We show that the data are consistent with two classes of solutions, characterized respectively by low-mass brown-dwarf $(q=0.037)$ and sub-Neptune $(q<10^{-4})$ companions, respectively. Future high-resolution imaging should easily distinguish between these solutions.

We analyze the binary gravitational microlensing event OGLE-2017-BLG-1130 (mass ratio q~0.45), the first published case in which the binary anomaly was only detected by the Spitzer Space Telescope. This event provides strong evidence that some binary signals can be missed by observations from the ground alone but detected by Spitzer. We therefore invert the normal procedure, first finding the lens parameters by fitting the space-based data and then measuring the microlensing parallax using ground-based observations. We also show that the normal four-fold space-based degeneracy in the single-lens case can become a weak eight-fold degeneracy in binary-lens events. Although this degeneracy is resolved in event OGLE-2017-BLG-1130, it might persist in other events.

We report the discovery, via the microlensing method, of a new very-low-mass binary system. By combining measurements from Earth and from the Spitzer telescope in Earth-trailing orbit, we are able to measure the microlensing parallax of the event, and find that the lens likely consists of an $(12.0 \pm 0.6) M_{\rm J}$ + $(15.7 \pm 1.5) M_{\rm J}$ super-Jupiter / brown-dwarf pair. The binary is located at a distance of $(3.08 \pm 0.18)$ kpc in the Galactic Plane, and the components have a projected separation of $(0.43 \pm 0.03)$ AU. Two alternative solutions with much lower likelihoods are also discussed, an 8- and 6-$M_{\rm J}$ model and a 90- and 70-$M_{\rm J}$ model. Although disfavored at the 3-$\sigma$ and 5-$\sigma$ levels, these alternatives cannot be rejected entirely. We show how the more-massive of these models could be tested with future direct imaging.