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

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

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.

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the MPIA 2.2m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measure time delays between all pairs of multiple images with only one or two seasons of monitoring at the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate is obtained for the delay between image A and image B of DES 0407-5006, with $\tau_{AB} = -128.4^{+3.5}_{-3.8}$ days (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combine our high-cadence data with measurements from decade-long light-curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 day on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily.

The massive star-forming region W75N~(B) is thought to host a cluster of massive protostars (VLA~1, VLA~2, and VLA~3) undergoing different evolutionary stages. In this work, we present radio continuum data with the highest sensitivity and angular resolution obtained to date in this region, using the VLA-A and covering a wide range of frequencies (4-48~GHz), which allowed us to study the morphology and the nature of the emission of the different radio continuum sources. We also performed complementary studies with multi-epoch VLA data and ALMA archive data at 1.3 mm wavelength. We find that VLA~1 is driving a thermal radio jet at scales of $\approx$0.1 arcsec ($\approx$130 au), but also shows signs of an incipient hyper-compact HII region at scales of $\lesssim$ 1 arcsec ($\lesssim$ 1300~au). VLA~3 is also driving a thermal radio jet at scales of a few tenths of arcsec (few hundred of au). We conclude that this jet is shock-exciting the radio continuum sources Bc and VLA~4 (obscured HH objects), which show proper motions moving outward from VLA~3 at velocities of $\approx$112--118~km/s. We have also detected three new weak radio continuum sources, two of them associated with millimeter continuum cores observed with ALMA, suggesting that these two sources are also embedded YSOs in this massive star-forming region.

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.

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.

The recent discovery of a millisecond radio burst from the Galactic magnetar SGR 1935+2154 has initiated an intense discussion about the connection between magnetars and fast radio bursts (FRBs). Although some properties of the radio burst from SGR 1935+2154 are not the same as cosmological FRBs, there are theoretical models which propose a connection between magnetars and FRBs (see review by Margalit et al. 2020). In particular, the role of a magnetar wind nebula is included in some models, and therefore it is worthwhile to investigate the X-ray environment of SGR 1935+2154 in more detail. Here, we report on the discovery of an X-ray transient feature near SGR 1935+2154 using archival Chandra data and discuss its possible origin.

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.

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

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

LDN 1157, is one of the several clouds situated in the cloud complex, LDN 1147/1158, represents a coma-shaped morphology with a well-collimated bipolar outflow emanating from a Class 0 protostar, LDN 1157-mm. The main goals of this work are (a) to map the inter-cloud magnetic field (ICMF) geometry of the region surrounding LDN 1157 to investigate its relationship with the cloud morphology, with the outflow direction and with the core magnetic field (CMF) geometry inferred from the mm- and sub-mm polarization results from the literature, and (b) to investigate the kinematic structure of the cloud. We carried out R-band polarization observations of the stars projected on the cloud to map the pc-scale magnetic field geometry and made spectroscopic observations of the entire cloud in 12CO, C18O and N2H+ (J=1-0) lines to investigate its kinematic structure. We obtained a distance of 340$\pm$3 pc to the LDN 1147/1158, complex based on the Gaia DR2 parallaxes and proper motion values of the three YSOs associated with the complex. A single filament of $\sim1.2$ pc in length and $\sim0.09$ pc in width is found to run all along the coma-shaped cloud. Based on the relationships between the ICMF, CMF, filament orientations, outflow direction, and the presence of an hour-glass morphology of the magnetic field, it is likely that the magnetic field had played an important role in the star formation process in LDN 1157. Combining the proper motions of the YSOs and the radial velocity of LDN 1147/1158 and another complex LDN 1172/1174 which is situated $\sim2$\dgr~east of it, we found that both the complexes are moving collectively toward the Galactic plane. The filamentary morphology of the east-west segment of LDN 1157 may have formed as a result of mass lost by ablation due to the interaction of the moving cloud with the ambient interstellar medium.

In this work, we study the X-ray bow-shock nebula powered by the mature pulsar PSR B1929+10 using data from XMM-Newton, with an effective exposure of $\sim$ 300 ks, offering the deepest investigation of this system thus far. We found the X-ray axial outflow extends as long as $\sim$ 8 arc minute behind the proper motion direction, which is a factor of two longer than the result reported in the previous study. Furthermore, we found evidence of two faint lateral outflows extending laterally with respect to the proper motion. We also found indications of spectral hardening along the axial outflow, suggesting that certain acceleration processes might occur along this feature.

We report the results of searching pulsar-like candidates from the unidentified objects in the $3^{\rm rd}$ Catalog of Hard Fermi-LAT sources (3FHL). Using a machine-learning based classification scheme with a nominal accuracy of $\sim98\%$, we have selected 27 pulsar-like objects from 200 unidentified 3FHL sources for an identification campaign. Using archival data, X-ray sources are found within the $\gamma-$ray error ellipses of 10 3FHL pulsar-like candidates. Within the error circles of the much better constrained X-ray positions, we have also searched for the optical/infrared counterparts and examined their spectral energy distributions. Among our short-listed candidates, the most secure identification is the association of 3FHL J1823.3-1339 and its X-ray counterpart with the globular cluster Mercer 5. The $\gamma-$rays from the source can be contributed by a population of millisecond pulsars residing in the cluster. This makes Mercer 5 as one of the slowly growing hard $\gamma-$ray population of globular clusters with emission $>10$ GeV. Very recently, another candidate picked by our classification scheme, 3FHL J1405.1-6118, has been identified as a new $\gamma-$ray binary with an orbital period of $13.7$ days. Our X-ray analysis with a short Chandra observation has found a possible periodic signal candidate of $\sim1.4$ hrs and a putative extended X-ray tail of $\sim20$ arcsec long. Spectral energy distribution of its optical/infrared counterpart conforms with a blackbody of $T_{\rm bb}\sim40000$ K and $R_{\rm bb}\sim12R_{\odot}$ at a distance of 7.7 kpc. This is consistent with its identification as an early O star as found by infrared spectroscopy.

We aim to investigate variations in the arrival time of coherent stellar pulsations due to the light-travel time effect to test for the presence of sub-stellar companions. Those companions are the key to one possible formation scenario of apparently single sub-dwarf B stars. We made use of an extensive set of ground-based observations of the four large amplitude p-mode pulsators DW Lyn, V1636 Ori, QQ Vir, and V541 Hya. Observations of the TESS space telescope are available on two of the targets. The timing method compares the phase of sinusoidal fits to the full multi-epoch light curves with phases from the fit of a number of subsets of the original time series. Observations of the TESS mission do not sample the pulsations well enough to be useful due to the (currently) fixed two-minute cadence. From the ground-based observations, we infer evolutionary parameters from the arrival times. The residual signals show many statistically significant periodic signals, but no clear evidence for changes in arrival time induced by sub-stellar companions. The signals can be explained partly by mode beating effects. We derive upper limits on companion masses set by the observational campaign.

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.

Fast radio bursts (FRBs) are mysterious radio bursts with a time scale of approximately milliseconds. Two populations of FRB, namely repeating and non-repeating FRBs, are observationally identified. However, the differences between these two and their origins are still cloaked in mystery. Here we show the time-integrated luminosity-duration ($L_{\nu}$-$w_{\rm int,rest}$) relations and luminosity functions (LFs) of repeating and non-repeating FRBs in the FRB Catalogue project. These two populations are obviously separated in the $L_{\nu}$-$w_{\rm int,rest}$ plane with distinct LFs, i.e., repeating FRBs have relatively fainter $L_{\nu}$ and longer $w_{\rm int,rest}$ with a much lower LF. In contrast with non-repeating FRBs, repeating FRBs do not show any clear correlation between $L_{\nu}$ and $w_{\rm int,rest}$. These results suggest essentially different physical origins of the two. The faint ends of the LFs of repeating and non-repeating FRBs are higher than volumetric occurrence rates of neutron-star mergers and accretion-induced collapse (AIC) of white dwarfs, and are consistent with those of soft gamma-ray repeaters (SGRs), type Ia supernovae, magnetars, and white-dwarf mergers. This indicates two possibilities: either (i) faint non-repeating FRBs originate in neutron-star mergers or AIC and are actually repeating during the lifetime of the progenitor, or (ii) faint non-repeating FRBs originate in any of SGRs, type Ia supernovae, magnetars, and white-dwarf mergers. The bright ends of LFs of repeating and non-repeating FRBs are lower than any candidates of progenitors, suggesting that bright FRBs are produced from a very small fraction of the progenitors regardless of the repetition. Otherwise, they might originate in unknown progenitors.

Investigating the Gunn-Peterson trough of high redshift quasars (QSOs) is a powerful way to reveal the cosmic reionization. As one of such attempts, we perform a series of analyses to examine the absorption lines observed with one of the highest redshift QSOs, PSO J006.1240+39.2219, which we previously discovered at z = 6.62. Using the Subaru telescope, we obtained medium-resolution spectrum with a total exposure time of 7.5 hours. We calculate the Ly$\alpha$ transmission in different redshift bins to determine the near zone radius and the optical depth at 5.6$<$z$<$6.5. We find a sudden change in the Ly$\alpha$ transmission at 5.75$<$z$<$5.86, which is consistent with the result from the literature. The near zone radius of the QSO is 5.79$\pm$0.09 $p$Mpc, within the scatter of the near zone radii of other QSOs measured in previous studies. We also analyze the dark gap distribution to probe the neutral hydrogen fractions beyond the saturation limit of the Gunn-Peterson trough. We extend the measurement of the dark gaps to 5.7$<$z$<$6.3. We find that the gap widths increase with increasing redshifts, suggesting more neutral Universe at higher redshifts. However, these measurements strongly depend on the continuum modeling. As a continuum model-free attempt, we also perform the dark-pixel counting analysis, to find the upper limit of $\langle x_{\rm H I} \rangle \sim$0.6 (0.8) at $z<$5.8 ($z>$5.8). All three analyses based on this QSO show increasingly neutral hydrogen towards higher redshifts, adding precious measurements up to z$\sim$6.5.

Within the so-called cosmological collider program, imprints of new particles on primordial non-Gaussianities have been studied intensively. In particular, their non-analytic features in the soft limit provide a smoking gun for new particles at the inflation scale. While this approach is very powerful to probe particles of the mass near the Hubble scale, the signal is exponentially suppressed for heavy particles. In this paper, to enlarge the scope of the cosmological collider, we explore a new approach to probing spins of heavy particles from signs of Wilson coefficients of the inflaton effective action and the corresponding primordial non-Gaussianities. As a first step, we focus on the regime where the de Sitter conformal symmetry is weakly broken. It is well known that the leading order effective operator $(\partial_\mu\phi\partial^\mu\phi)^2$ is universally positive as a consequence of unitarity. In contrast, we find that the sign of the six derivative operator $(\nabla_\mu\partial_\nu\phi)^2(\partial_\rho\phi)^2$ is positive for intermediate heavy scalars, whereas it is negative for intermediate heavy spinning states. Therefore, under the assumption of tree-level UV completion, the sign can be used to probe spins of heavy particles generating the effective interaction. We also study phenomenology of primordial non-Gaussianities thereof.

We show that dynamical evolution in a strong (Galactic Centre-like) tidal field can create clusters that would appear to have very top-heavy IMFs. The tidal disruption of single star forming events can leave several bound `clusters' spread along 20 pc of the orbit within 1-2 Myr. These surviving (sub)clusters tend to contain an over-abundance of massive stars, with low-mass stars tending to be spread along the whole `tidal arm'. Therefore observing a cluster in a strong tidal field with a top-heavy IMF might well not mean the stars formed with a top-heavy IMF.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

The intergalactic medium (IGM) at $z\sim$ 5 to 6 is largely ionized, and yet the main source for the IGM ionization in the early universe is uncertain. Of the possible contributors are faint quasars with $-26 \lesssim M_{\rm 1450} \lesssim -23$, but their number density is poorly constrained at $z\sim5$. In this paper, we present our survey of faint quasars at $z\sim5$ in the European Large-Area {\it ISO} Survey-North 1 (ELAIS-N1) field over a survey area of 6.51 deg$^2$ and examine if such quasars can be the dominant source of the IGM ionization. We use the deep optical/near-infrared data of the ELAIS-N1 field as well as the additional medium-band observations to find $z \sim 5$ quasars through a two-step approach using the broadband color selection, and SED fitting with the medium-band information included. Adopting Bayesian information criterion, we identify ten promising quasar candidates. Spectra of three of the candidates are obtained, confirming all of them to be quasars at $z\sim5$ and supporting the reliability of the quasar selection. Using the promising candidates, we derive the $z\sim5$ quasar luminosity function at $-26 \lesssim M_{\rm 1450} \lesssim -23$. The number density of faint $z\sim5$ quasars in the ELAIS-N1 field is consistent with several previous results that quasars are not the main contributors to the IGM-ionizing photons at $z\sim5$

The CS molecule is known to be absorbed onto dust in the cold and dense conditions, causing it to get significantly depleted in the central region of cores. This study is aimed to investigate the depletion of the CS molecule using the optically thin C$^{34}$S molecular line observations. We mapped five prestellar cores, L1544, L1552, L1689B, L694-2, and L1197 using two molecular lines, C$^{34}$S $(J=2-1)$ and N$_2$H$^+$ $(J=1-0)$ with the NRO 45-m telescope, doubling the number of cores where the CS depletion was probed using C$^{34}$S. In most of our targets, the distribution of C$^{34}$S emission shows features that suggest that the CS molecule is generally depleted in the center of the prestellar cores. The radial profile of the CS abundance with respect to H$_2$ directly measured from the CS emission and the Herschel dust emission indicates that the CS molecule is depleted by a factor of $\sim$3 toward the central regions of the cores with respect to their outer regions. The degree of the depletion is found to be even more enhanced by an order of magnitude when the contaminating effect introduced by the presence of CS molecules in the surrounding envelope that lie along the line-of-sight is removed. Except for L1197 which is classified as relatively the least evolved core in our targets based on its observed physical parameters, we found that the remaining four prestellar cores are suffering from significant CS depletion at their central region regardless of the relative difference in their evolutionary status.

We first present new Subaru narrow-band observations of the Ly$\alpha$ halo around the quasi-stellar object (QSO) CFHQ J232908$-$030158 at $z=6.42$, which appears the most luminous and extended halo at $z>5$ ($L_{Ly\alpha}=9.8\times10^{43}$ erg s$^{-1}$ within $37$ pkpc diameter). Then, combining these measurements with available data in the literature, we find two different evolutions of QSOs' Ly$\alpha$ halos. First is a possible short-term evolution with QSO age seen in four $z>6$ QSOs. We find the anti-correlation between the Ly$\alpha$ halo scales with QSOs' IR luminosity, with J2329-0301's halo being the brightest and largest. It indicates that ionizing photons escape more easily out to circum-galactic regions when host galaxies are less dusty. We also find a positive correlation between IR luminosity and black hole mass ($M_\text{BH}$). Given $M_\text{BH}$ as an indicator of QSO age, we propose a hypothesis that a large Ly$\alpha$ halo mainly exists around QSOs in the young phase of their activity due to a small amount of dust. The second is an evolution with cosmic time seen over $z\sim2-5$. We find the increase of surface brightness toward lower-redshift with a similar growth rate to that of dark matter halos (DHs) which evolve to $M_\text{DH}=10^{12}-10^{13}$ M$_\odot$ at $z=2$. The extent of Ly$\alpha$ halos is also found to increase at a rate scaling with the virial radius of growing DHs, $r_\text{vir} \propto M_\text{DH}^{1/3}$($1+z$)$^{-1}$. These increases are consistent with a scenario that the CGM around QSOs evolves in mass and size keeping pace with hosting DHs.

We investigate broadband emission properties of the pulsar wind nebula (PWN) 3C 58 using a spectral energy distribution (SED) model. We attempt to match simultaneously the broadband SED and spatial variations of X-ray emission in the PWN. We further the model to explain a possible far-IR feature of which a hint is recently suggested in 3C 58: a small bump at $\sim$$10^{11}$ GHz in the PLANCK and Herschel band. While external dust emission may easily explain the observed bump, it may be internal emission of the source implying an additional population of particles. Although significance for the bump is not high, here we explore possible origins of the IR bump using the emission model and find that a population of electrons with GeV energies can explain the bump. If it is produced in the PWN, it may provide new insights into particle acceleration and flows in PWNe.

We present the development of Linear Astigmatism Free - Three Mirror System (LAF-TMS). This is a prototype of an off-axis telescope that enables very wide field of view (FoV) infrared satellites that can observe Paschen-$\alpha$ emission, zodiacal light, integrated star light, and other infrared sources. It has the entrance pupil diameter of 150 mm, the focal length of 500 mm, and the FoV of 5.5$^\circ$ $\times$ 4.1$^\circ$. LAF-TMS is an obscuration-free off-axis system with minimal out-of-field baffling and no optical support structure diffraction. This optical design is analytically optimized to remove linear astigmatism and to reduce high-order aberrations. Sensitivity analysis and Monte-Carlo simulation reveal that tilt errors are the most sensitive alignment parameters that allow $\sim$1$^\prime$. Optomechanical structure accurately mounts aluminum mirrors, and withstands satellite-level vibration environments. LAF-TMS shows optical performance with 37 $\mu$m FWHM of the point source image satisfying Nyquist sampling requirements for typical 18 $\mu$m pixel Infrared array detectors. The surface figure errors of mirrors and scattered light from the tertiary mirror with 4.9 nm surface micro roughness may affect the measured point spread function (PSF). Optical tests successfully demonstrate constant optical performance over wide FoV, indicating that LAF-TMS suppresses linear astigmatism and high-order aberrations.

We consider gravitational radiation and electromagnetic radiation from point mass binary with electric charges in a Keplerian orbit, and calculate the merger rate distribution of primordial black hole binaries with charges and a general mass function by taking into account gravitational torque and electromagnetic torque by the nearest primordial black hole. We apply the formalism to the extremal charged case and find that $\alpha=-(m_i+m_j)^2\partial^2 \ln {\cal R}(m_i,m_j)/\partial m_i \partial m_j=12/11$, which is independent of the mass function.

Large liquid argon detectors offer one of the best avenues for the detection of galactic weakly interacting massive particles (WIMPs) via their scattering on atomic nuclei. The liquid argon target allows exquisite discrimination between nuclear and electron recoil signals via pulse-shape discrimination of the scintillation signals. Atmospheric argon (AAr), however, has a naturally occurring radioactive isotope, $^{39}$Ar, a $\beta$ emitter of cosmogenic origin. For large detectors, the atmospheric $^{39}$Ar activity poses pile-up concerns. The use of argon extracted from underground wells, deprived of $^{39}$Ar, is key to the physics potential of these experiments. The DarkSide-20k dark matter search experiment will operate a dual-phase time projection chamber with 50 tonnes of radio-pure underground argon (UAr), that was shown to be depleted of $^{39}$Ar with respect to AAr by a factor larger than 1400. Assessing the $^{39}$Ar content of the UAr during extraction is crucial for the success of DarkSide-20k, as well as for future experiments of the Global Argon Dark Matter Collaboration (GADMC). This will be carried out by the DArT in ArDM experiment, a small chamber made with extremely radio-pure materials that will be placed at the centre of the ArDM detector, in the Canfranc Underground Laboratory (LSC) in Spain. The ArDM LAr volume acts as an active veto for background radioactivity, mostly $\gamma$-rays from the ArDM detector materials and the surrounding rock. This article describes the DArT in ArDM project, including the chamber design and construction, and reviews the background required to achieve the expected performance of the detector.

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

Young planets with masses approaching Jupiter's have tides strong enough to clear gaps around their orbits in the protostellar disk. Gas flow through the gaps regulates the planets' further growth and governs the disks' evolution. Magnetic forces may drive that flow if the gas is sufficiently ionized to couple to the fields. We compute the ionizing effects of the X-rays from the central young star, using Monte Carlo radiative transfer calculations to find the spectrum of Compton-scattered photons reaching the planet's vicinity. The scattered X-rays ionize the gas at rates similar to or greater than the interstellar cosmic ray rate near planets the mass of Saturn and of Jupiter, located at 5 au and at 10 au, in disks with the interstellar mass fraction of sub-micron dust and with the dust depleted a factor 100. Solving a gas-grain recombination reaction network yields charged particle populations whose ability to carry currents is sufficient to partly couple the magnetic fields to the gas around the planet. Most cases can undergo Hall shear instability, and some can launch magnetocentrifugal winds. However the material on the planet's orbit has diffusivities so large in all the cases we examine, that magneto-rotational turbulence is prevented and the non-ideal terms govern the magnetic field's evolution. Thus the flow of gas in the gaps opened by the young giant planets depends crucially on the finite conductivity.

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.

The preacceleration of electrons through reflection and shock drift acceleration (SDA) is essential for the diffusive shock acceleration (DSA) of nonthermal electrons in collisionless shocks. Previous studies suggested that, in weak quasi-perpendicular ($Q_\perp$) shocks in the high-$\beta$ ($\beta=P_{\rm gas}/P_{\rm B}$) intracluster medium (ICM), the temperature anisotropy due to SDA-reflected electrons can drive the electron firehose instability, which excites oblique nonpropagating waves in the shock foot. In this paper, we investigate, through a linear analysis and particle-in-cell (PIC) simulations, the firehose instabilities driven by an electron temperature anisotropy (ETAFI) and also by a drifting electron beam (EBFI) in $\beta\sim100$ ICM plasmas. The EBFI should be more relevant in describing the self-excitation of upstream waves in $Q_\perp$-shocks, since backstreaming electrons in the shock foot behave more like an electron beam rather than an anisotropic bi-Maxwellian population. We find that the basic properties of the two instabilities, such as the growth rate, $\gamma$, and the wavenumber of fast-growing oblique modes are similar in the ICM environment, with one exception; while the waves excited by the ETAFI are nonpropagating ($\omega_r=0$), those excited by the EBFI have a non-zero frequency ($\omega_r\neq0$). However, the frequency is small with $\omega_r<\gamma$. Thus, we conclude that the interpretation of previous studies for the nature of upstream waves based on the ETAFI remains valid in $Q_\perp$-shocks in the ICM.

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.

Compilation of papers presented by the JEM-EUSO Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

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

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.

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works.

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

We examine whether galaxy environments directly affect triggering nuclear activity in Sloan Digital Sky Survey (SDSS) local spiral galaxies using a volume-limited sample with the $r$-band absolute magnitude $M_{r} < -19.0$ and $0.02 < z < 0.055$ selected from the SDSS Data Release 7. To avoid incompleteness of the central velocity dispersion $\sigma$ of the volume-limited sample and to fix the black hole mass affecting AGN activity, we limit the sample to a narrow $\sigma$ range of $130$ km s$^{-1}<\sigma<200$ km s$^{-1}$. We define a variety of environments as a combination of neighbour interactions and local density on a galaxy. After the central star formation rate (which is closely related to AGN activity level) is additionally restricted, the direct impact of the environment is unveiled. In the outskirts of rich clusters, red spiral galaxies show a significant excess of the AGN fraction despite the lack of central gas. We argue that they have been pre-processed before entering the rich clusters, and due to mergers or strong encounters in the in-fall region, their remaining gases efficiently lose angular momentum. We investigate an environment in which many star-forming galaxies coexist with a few starburst-AGN composite hosts having the highest [OIII] luminosity. We claim that they are a gas-rich merger product in groups or are group galaxies in-falling into clusters, indicating that many AGN signatures may be obscured following the merger events.

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.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

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 conducted a survey of open clusters within 1 kpc from the Sun using the astrometric and photometric data of the Gaia Data Release 2. We found 655 cluster candidates by visual inspection of the stellar distributions in proper motion space and spatial distributions in l-b space. All of the 655 cluster candidates have a well defined main-sequence except for two candidates if we consider that the main sequence of very young clusters is somewhat broad due to differential extinction. Cross-matching of our 653 open clusters with known open clusters in various catalogs resulted in 207 new open clusters. We present the physical properties of the newly discovered open clusters. The majority of the newly discovered open clusters are of young to intermediate age and have less than ~50 member stars.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

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

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

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.

In our recent report, observational evidence supports that the rotational direction of a galaxy tends to be coherent with the average motion of its nearby neighbors within 1 Mpc. We extend the investigation to neighbors at farther distances, in order to examine if such dynamical coherence is found even in large scales. The Calar Alto Legacy Integral Field Area (CALIFA) survey data and the NASA-Sloan Atlas (NSA) catalog are used. From the composite map of velocity distribution of 'neighbor' galaxies within 15 Mpc from the CALIFA galaxies, the composite radial profiles of the luminosity-weighted mean velocity of neighbors are derived. These profiles show unexpectedly strong evidence of the dynamical coherence between the rotation of the CALIFA galaxies and the average line-of-sight motion of their neighbors within several Mpc distances. Such a signal is particularly strong when the neighbors are limited to red ones: the luminosity-weighted mean velocity at 1 < D <= 6 Mpc is as large as 30.6+/-10.9 km/s (2.8-sigma significance to random spin-axis uncertainty) for central rotation (R <= Re). In the comparison of several subsamples, the dynamical coherence tends to be marginally stronger for the diffuse or kinematically-well-aligned CALIFA galaxies. For this mysterious coherence in large scales, we cautiously suggest a scenario that it results from a possible relationship between the long-term motion of a large-scale structure and the rotations of galaxies in it.

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.

Event rate of long Gamma-Ray Bursts (GRBs) is expected to be an useful tracer of the cosmic star-formation history. For this purpose, it is necessary to understand what kind of star formations/galaxies are traced by GRBs. Here we report rest-frame far-infrared (FIR) continuum detections of GRB 070521 and 080207 host galaxies at $z\sim2$ with ALMA band 8 and 9. The FIR photometries provide the reliable star-formation rates (SFRs), because FIR emission is free from dust extinction and possible radio contamination from long-lived afterglows of GRBs. The spectral energy distribution fittings indicate 49.85$^{+72.33}_{-2.86}$ and 123.4$^{+25.19}_{-21.78}$ $M_{\odot}$ yr$^{-1}$ for the 070521 and 080207 hosts, respectively. The derived SFRs place them on the \lq \lq main sequence\rq \rq\ of normal star-forming galaxies at $z\sim2$. The derived SFRs are significantly lower than that of radio observations. It is inferred that the observed radio fluxes in a previous study are contaminated by the afterglows. ALMA marginally detected [C~{\sc ii}]\,158\,$\mu$m emission line from the GRB 080207 host galaxy with S/N $\sim$ 4. This is the first detection of [C~{\sc ii}]\,158\,$\mu$m of a GRB host at $z>2$, and the second detection among known GRBs. The luminosity ratio of [C~{\sc~ii}]\,158$\mu$m to FIR is 7.5$\times 10^{-4}$, which is one of the smallest values among galaxies at $z\sim 1-2$ with the same FIR luminosity. The \lq \lq [C~{\sc ii}] deficit\rq \rq\ could be a new physical property to characterise GRB hosts at $z\sim1-2$. Possible parameters controlling the deficit include the metallicity, initial mass function, and gas density.

Galaxies evolve from a blue star-forming phase into a red quiescent one by quenching their star formation activity. In high density environments, this galaxy evolution proceeds earlier and more efficiently. Therefore, local galaxy clusters are dominated by well-evolved red, elliptical galaxies. The fraction of blue galaxies in clusters monotonically declines with decreasing redshift, i.e., the Butcher-Oemler effect. In the local Universe, observed blue fractions of massive clusters are as small as $\lesssim$ 0.2. Here we report a discovery of a \lq \lq blue cluster\rq \rq, that is a local galaxy cluster with an unprecedentedly high fraction of blue star-forming galaxies yet hosted by a massive dark matter halo. The blue fraction is 0.57, which is 4.0 $\sigma$ higher than those of the other comparison clusters under the same selection and identification criteria. The velocity dispersion of the member galaxies is 510 km s$^{-1}$, which corresponds to a dark matter halo mass of 2.0$^{+1.9}_{-1.0}\times 10^{14}$ M$_{\odot}$. The blue fraction of the cluster is more than 4.7 $\sigma$ beyond the standard theoretical predictions including semi-analytic models of galaxy formation. The probability to find such a high blue fraction in an individual cluster is only 0.003\%, which challenges the current standard frameworks of the galaxy formation and evolution in the $\Lambda$CDM Universe. The spatial distribution of galaxies around the blue cluster suggests that filamentary cold gas streams can exist in massive halos even in the local Universe. However these cold streams have already disappeared in the theoretically simulated local universes.

Nature of dark energy remains unknown. Especially, to constrain the time variability of the dark-energy, a new, standardisable candle that can reach more distant Universe has been awaited. Here we propose a new distance measure using fast radio bursts (FRBs), which are a new emerging population of $\sim$ ms time scale radio bursts that can reach high-$z$ in quantity. We show an empirical positive correlation between the time-integrated luminosity (L$_{\nu}$) and rest-frame intrinsic duration ($w_{\rm int,rest}$) of FRBs. The L$_{\nu}-w_{\rm int,rest}$ correlation is with a weak strength but statistically very significant, i.e., Pearson coefficient is $\sim$ 0.5 with p-value of $\sim$0.038, despite the smallness of the current sample. This correlation can be used to measure intrinsic luminosity of FRBs from the observed $w_{\rm int,rest}$. By comparing the luminosity with observed flux, we measure luminosity distances to FRBs, and thereby construct the Hubble diagram. This FRB cosmology with the L$_{\nu}-w_{\rm int,rest}$ relation has several advantages over SNe Ia, Gamma-Ray Burst (GRB), and well-known FRB dispersion measure (DM)-$z$ cosmology; (i) access to higher redshift Universe beyond the SNe Ia, (ii) high event rate that is $\sim$ 3 order of magnitude more frequent than GRBs, and (iii) it is free from the uncertainty from intergalactic electron density models, i.e., we can remove the largest uncertainty in the well-debated DM-$z$ cosmology of FRB. Our simulation suggests that the L$_{\nu}-w_{\rm int,rest}$ relation provides us with useful constraints on the time variability of the dark energy when the next generation radio telescopes start to find FRBs in quantity.

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.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

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.

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

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

We investigate the stellar population properties of passive spiral galaxies in the CALIFA survey. Nine spiral galaxies that have NUV-r > 5 and no/weak nebular emission lines in their spectra are selected as passive spirals. Our passive spirals lie in the redshift range of 0.001 < z < 0.021 and have stellar mass range of 10.2 < log(M_{\star}/M_{\odot}) < 10.8. They clearly lie in the domain of early-type galaxies in the WISE IR color-color diagram. We analyze the stellar populations out to two effective radius, using the best-fitting model to the measured absorption line-strength indices in the Lick/IDS system. We find that stellar populations of the passive spirals span a wide range, even in their centers, and hardly show any common trend amongst themselves either. We compare the passive spirals with S0s selected in the same mass range. S0s cover a wide range in age, metallicity, and [a/Fe], and stellar populations of the passive spirals are encompassed in the spread of the S0 properties. However, the distribution of passive spirals are skewed toward higher values of metallicity, lower [a/Fe], and younger ages at all radii. These results show that passive spirals are possibly related to S0s in their stellar populations. We infer that the diversity in the stellar populations of S0s may result from different evolutionary pathways of S0 formation, and passive spirals may be one of the possible channels.

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

The presence of blue cores in some dwarf early-type galaxies (dEs) in galaxy clusters suggests the scenario of late-type galaxy infall and subsequent transformation into red, quiescent dEs. We present Gemini Multi-Object Spectrographs long-slit spectroscopy of two dEs with blue cores (dE(bc)s), EVCC 591 and EVCC 516, located at the core and outskirt of the Virgo cluster, respectively. We obtained their internal kinematics along the major axis out to, at least, ~ 1 effective radius. EVCC 591 shows evidence of a kinematically decoupled core (KDC) with a size of 2" (160 pc), exhibiting an inverted pattern for velocity with respect to the main body of its host galaxy. The rotation curve of the stellar component in the inner region of EVCC 591 is steeper than that in the rest of the galaxy. On the other hand, overall velocity profiles of stellar and ionized gas components of EVCC 516 show no signature of significant rotation. The occurrence of a KDC and zero rotation in the internal kinematics along with the central star formation support the scenario of gas-rich dwarf-dwarf mergers in the formation of these two dE(bc)s. Furthermore, evolution of dE(bc)s in a cluster environment into ordinary dEs with KDCs is possible based on their structural properties. We suggest that at least some of the dE(bc)s in the Virgo cluster were formed through dwarf-dwarf mergers in lower density environments before they subsequently fell into the cluster; they were then quenched by subsequent effects within the cluster environment.

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.

Recent observations of gravitational waves from binary mergers of black holes or neutron stars and the rapid development of ultra-intense lasers lead strong field physics to a frontier of new physics in the 21st century. Strong gravity phenomena are most precisely described by general relativity, and lasers that are described by another most precisely tested quantum electrodynamics (QED) can be focused into a tiny area in a short period through the chirped pulse amplification and generate extremely high intensity electromagnetic (EM) fields beyond the conventional methods. It is physically interesting to study QED phenomena in curved spacetimes, in which both strong gravitational and electromagnetic fields play important roles. There are many sources for strong gravitational and electromagnetic fields in the sky or universe, such highly magnetized neutron stars, magnetized black holes, and the early universe. We review quantum field theoretical frameworks for QED both in the Minkowski spacetime and curved spacetimes, in particular, charged black holes and the early universe, and discuss the QED physics in strong EM fields, such as the vacuum polarization and Schwinger pair production and their implications to astrophysics and cosmology.

We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723.

The flat spectrum radio quasar 4C 38.41 showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional monitoring observations at various bands across the electromagnetic spectrum. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components and the evolution of the physical properties seem to be in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1 mG, do not seem to dramatically change during the flares, and are much smaller, by a factor 10,000, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting region.

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

"TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C$^{18}$O (1-0) mainly tracing low density cloud and about 460 square arcminute area using N$_{2}$H$^{+}$ (1-0) mainly tracing dense cores. CS (2-1) and SO $(3_{2}-2_{1})$ were also used simultaneously to map $\sim$440 square arcminute area of this region. We identified 10 filaments by applying the dendrogram technique to the C$^{18}$O data-cube and 8 dense N$_{2}$H$^{+}$ cores by using {\sc FellWalker}. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that L1478 consists of several filaments with slightly different velocities. Especially the filaments which are supercritical are found to contain dense cores detected in N$_{2}$H$^{+}$. Comparison of non-thermal velocity dispersions derived from C$^{18}$O and N$_{2}$H$^{+}$ for the filaments and dense cores indicates that some of dense cores share similar kinematics with those of the surrounding filaments while several dense cores have different kinematics with those of their filaments. This suggests that the formation mechanism of dense cores and filaments can be different in individual filaments depending on their morphologies and environments.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

An increasing number of white dwarf stars show atmospheric chemical composition polluted by heavy elements accreted from debris disk material. The existence of such debris disks strongly suggests the presence of one or more planet(s) whose gravitational interaction with rocky planetesimals is responsible for their disruption by tidal effect. The ZZ Ceti pulsator and polluted DAZ white dwarf GD 133 is a good candidate for searching for such a potential planet. We started in 2011 a photometric follow-up of its pulsations. As a result of this work in progress, we used the data gathered from 2011 to 2015 to make an asteroseismological analysis of GD 133, providing the star parameters from a best fit model with $M$/$M_{\odot}$ = 0.630 $\pm$ 0.002, $T_{\rm eff}$ = 12400 K $\pm$ 70 K, log($M_{\rm He}/M$) = -2.00 $\pm$ 0.02, log($M_{\rm H}/M$) = -4.50 $\pm$ 0.02 and determining a rotation period of $\approx$ 7 days.

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

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 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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 study dressed inflation with a cuscuton and find a novel exact power-law solution. It is well known that the conventional power-law inflation is inconsistent with the Planck data. In contrast to this standard lore, we find that power-law inflation with a cuscuton can be reconciled with the Planck data. Moreover, we argue that the cuscuton generally ameliorates inflation models so that predictions are consistent with observations.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

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

A linear correlation has been proposed between the CO luminosity ($\rm{L}^{\prime}_{\rm{CO}}$) and full-width at half maximum (FWHM) for high-redshift (z > 1) submillimeter galaxies. However, the controversy concerning the $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlation seems to have been caused by the use of heterogeneous samples (e.g., different transition lines) and/or data with large measurement uncertainties. In order to avoid the uncertainty caused by using different rotational transitions, in this work we make an extensive effort to select only CO($J = 1-0$) data from the literature. We separate these wide-ranging redshift data into two samples : the low-redshift (z < 1) and high-redshift (z > 1) samples. The samples are corrected for lensing magnification factors if gravitational-lensing effects appeared in the observations. The correlation analysis shows that there exists significant $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlations for both the low-redshift and high-redshift samples. A comparison of the low- and high-redshift $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlations does not show strong evolution with redshift. Assuming that there is no evolution, we can use this relation to determine the model independent distances of high-redshift galaxies. We then constrain cosmological models with the calibrated high-redshift CO data and the sample of Type Ia supernovae in the Union 2.1 compilation. In the constraint for wCDM with our samples, the derived values are w_{0} = -1.02 {\pm} 0.17, {\Omega}_{m0} = 0.30{\pm}0.02, and H_{0} = 70.00 {\pm}0.60 km\,s^{-1}\,Mpc^{-1}.

Revealing what fraction of galaxies harbor AGN is central in understanding black hole accretion history of the Universe. However, optical and soft X-ray surveys miss the most highly obscured AGNs. Infrared (IR), instead, is more robust against absorption. Previous IR photometric surveys, however, only had 4 or 5 filters in mid-IR. Our AKARI North Ecliptic Pole (NEP) wide field sample has 18 filters in mid-IR (9 from AKARI, 4 from WISE, and 5 from Spitzer), for the first time, allowing a sophisticated mid-IR SED fitting diagnosis for a statistical number of sources (89178 over 5.4 deg$^2$). By using a SED fitting technique, we investigate the evolution of AGN fraction as a function of redshift and IR (8-1000 $\mu$m) luminosity in an extinction-free way. We found that the AGN fraction (F$_{\rm AGN}$) shows no sign of strong redshift evolution. Instead, F$_{\rm AGN}$ increases with increasing IR luminosity in all redshifts bins ($0<z<2$).

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.

Evolutionary properties of infrared (IR) luminous galaxies are important keys to understand dust-obscured star formation history and galaxy evolution. Based on the near- to mid-IR imaging with 9 continuous filters of AKARI space telescope, we present the characteristics of dusty star-forming (SF) galalxies showing polycyclic aromatic hydrocarbon (PAH) features observed by the North Ecliptic Pole (NEP) wide field survey of AKARI and Herschel. All the sample galaxies from the AKARI/NEP-Wide data are selected based both on the Herschel/SPIRE 250 {\mu}m detection and optical spectroscopic redshift data. The physical modelling of spectral energy distribution (SED) using all available data points from u to sub-mm 500 {\mu}m band, including WISE and PACS data where available, takes unique advantages of the continuous near to mid-IR coverage, reliable constraint on far-IR peak, spectroscopically determined accurate redshifts, as well as energy balance principle by MAGPHYS. This enables us to derive physically meaningful and accurate total infrared luminosity and 8 {\mu}m (or PAH) luminosity consistently. Our sample galaxies are in the redshift range z <1, and majority of them appear to be normal SF/spiral populations showing PAH features near the 8 {\mu}m. These SF galaxies showing PAHs in the mid-IR include various types from quiescent to starbursts. Some of our sample show shortage of 8 {\mu}m luminosity compared to the total IR luminosity and this PAH deficit gets severe in more luminous IR galaxies, suggesting PAH molecules in these galaxies destroyed by strong radiation field from SF region or a large amount of cold dust in ISM. The specific SFR of our sample shows mass dependent time evolution which is consistent with downsizing evolutionary pattern.

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.

Analyzing multiband observations of dust continuum emission is one of the useful tools to constrain dust properties which help us to understand the physical properties of the disks. We perform the synthetic ALMA multiband analysis to find the best ALMA band set for constraining the dust properties of the TW Hya protoplanetary disk. We find that the Band [10,6,3] set is the best set among the possible combinations of ALMA Band [3,4,5,6,7,8,9,10]. We also find two conditions for the good ALMA band sets providing narrow constraint ranges on dust properties; (1) Band 9 or 10 is included in the band set and (2) Enough frequency intervals between the bands. These are related with the conditions which give good constraints on dust properties: the combination of optically thick and thin bands are required, and large beta (beta is the power-law index of dust opacity), and low dust temperature are preferable. To examine our synthetic analysis results, we apply the multiband analysis to ALMA archival data of the TW Hya disk at Band 4, 6, 7, and 9. Band [9,6,4] set provides the dust properties close to the model profile, while Band [7,6,4] set gives the dust properties deviating from the model at all radii with too broad constraint range to specify the accurate values of dust temperature, optical depth, and \b{eta}. Since these features are expected by the synthetic multiband analysis, we confirm that the synthetic multiband analysis is well consistent with the results derived from real data.

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.

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.

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.

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.

We present the analysis of a new near-infrared (NIR) spectrum of a recently discovered $z=6.621$ quasar PSO J006+39 in an attempt to explore the early growth of supermassive black holes (SMBHs). This NIR (rest-frame ultraviolet, UV) spectrum shows blue continuum slope and rich metal emission lines in addition to Ly$\alpha$ line. We utilize the MgII line width and the rest frame luminosity $L_\text{3000\AA}$ to find the mass of SMBH ($M_\text{BH}$) to be $\sim 10^8 M_\odot$, making this one of the lowest mass quasars at high redshift. The power-law slope index ($\alpha_\lambda$) of the continuum emission is $-2.94\pm0.03$, significantly bluer than the slope of $\alpha_\lambda=-7/3$ predicted from standard thin disc models. We fit the spectral energy distribution (SED) using a model which can fit local SMBHs, which includes warm and hot Comptonisation powered by the accretion flow as well as an outer standard disc. The result shows that the very blue slope is probably produced by a small radial ($\sim230$ gravitational radius, $R_\text{g}$) extent of the standard accretion disc. All plausible SED models require that the source is super-Eddington ($L_\text{bol}/L_\text{Edd} \gtrsim 9$), so the apparently small disc may simply be the inner funnel of a puffed up flow, and clearly the SMBH in this quasar is in a rapid growth phase. We also utilize the rest-frame UV emission lines to probe the chemical abundance in the broad line region (BLR) of this quasar. We find that this quasar has super solar metallicity through photoionization model calculations.

We present the results of a systematic search for molecular outflows in 68 Very Low Luminosity Objects (VeLLOs) from single-dish observations in CO isotopologues which find 16 VeLLOs showing clear outflow signatures in the CO maps. With additional three VeLLOs from the literature, we analyzed the outflow properties for 19 VeLLOs, identifying 15 VeLLOs as proto-Brown Dwarf (BD) candidates and four VeLLOs as likely faint protostar candidates. The proto-BD candidates are found to have a mass accretion rate ($\sim 10^{-8} - 10^{-7}$ $\rm M_{\odot}$ yr$^{-1}$) lower than that of the protostar candidates ($\gtrsim 10^{-6}$ $\rm M_{\odot}$ yr$^{-1}$). Their accretion luminosities are similar to or smaller than their internal luminosities, implying that many proto-BD candidates might have had either small accretion activity in a quiescent manner throughout their lifetime, or be currently exhibiting a relatively higher (or episodic) mass accretion than the past. There are strong trends that outflows of many proto-BDs are less active if they are fainter or have less massive envelopes. The outflow forces and internal luminosities for more than half of the proto-BD candidates seem to follow an evolutionary track of a protostar with its initial envelope mass of $\sim$0.08 $\rm M_{\odot}$, indicating that some BDs may form in less massive dense cores in a way similar to normal stars. But, because there also exists a significant fraction (about 40%) of proto-BDs with much weaker outflow force than expected by the relations for protostars, we should not rule out the possibility of other formation mechanism for the BDs.

Ultra-high-energy cosmic rays (UHECRs) are known to come from outside of our Galaxy, but their origin still remains unknown. The Telescope Array (TA) experiment recently identified a high concentration in the arrival directions of UHECRs with energies above $5.7 \times 10^{19} eV$, called hotspot. We here report the presence of filaments of galaxies, connected to the Virgo Cluster, on the sky around the hotspot, and a statistically significant correlation between hotspot events and the filaments. With 5-year TA data, the maximum significance of binomial statistics for the correlation is estimated to be 6.1 $\sigma$ at correlation angle 3.4 degree. The probability that the above significance appears by chance is $\sim 2.0 \times 10^{-8}$ (5.6 $\sigma$). Based on this finding, we suggest a model for the origin of TA hotspot UHECRs; they are produced at sources in the Virgo Cluster, and escape to and propagate along filaments, before they are scattered toward us. This picture requires the filament magnetic fields of strength $\gtrsim 20$ nG, which need to be confirmed in future observations.

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 study time evolution of Schrodinger-Newton system using the self-consistent Crank-Nicolson method to understand the dynamical characteristics of nonlinear systems. Compactifying the radial coordinate by a new one, which brings the spatial infinity to a finite value, we are able to impose the boundary condition at infinity allowing for a numerically exact treatment of the Schrodinger-Newton equation. We study patterns of gravitational cooling starting from exponentially localized initial states. When the gravitational attraction is strong enough, we find that a small-sized oscillatory solitonic core is forming quickly, which is surrounded by a growing number of temporary halo states. In addition a significant fraction of particles escape to asymptotic regions. The system eventually settles down to a stable solitonic core state while all the excess kinetic energy is carried away by the escaping particles, which is a phenomenon of gravitational cooling.

The short gamma-ray burst (GRB) 170817A was the first GRB associated with a gravitational-wave event. Due to the exceptionally low luminosity of the prompt $\gamma$-ray and the afterglow emission, the origin of both radiation components is highly debated. The most discussed models for the burst and the afterglow include a regular GRB jet seen off-axis and the emission from the cocoon encompassing a "choked" jet. Here, we report low radio-frequency observations at 610 and 1390~MHz obtained with the Giant Metrewave Radio Telescope (GMRT). Our observations span a range of $\sim7$ to $\sim152$ days after the burst. The afterglow started to emerge at these low frequencies about 60~days after the burst. The $1390$~MHz light curve barely evolved between 60 and 150 days, but its evolution is also marginally consistent with a $F_\nu\propto t^{0.8}$ rise seen in higher frequencies. We model the radio data and archival X-ray, optical and high-frequency radio data with models of top-hat and Gaussian structured GRB jets. We performed a Markov Chain Monte Carlo analysis of the structured-jet parameter space. Though highly degenerate, useful bounds on the posterior probability distributions can be obtained. Our bounds of the viewing angle are consistent with that inferred from the gravitational wave signal. We estimate the energy budget in prompt emission to be an order of magnitude lower than that in the afterglow blast-wave.

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

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

The faint quasars with $M_{1450}>-24$ mag are known to hold the key to the determination of the ultraviolet emissivity for the cosmic re-ionization. But only a few have been identified so far because of the limitations on the survey data. Here, we present the first results of the $z\sim5$ faint quasar survey with the Infrared Medium-deep Survey (IMS), which covers $\sim100$ deg$^{2}$ areas in $J$-band to the depths of $J_{\rm AB}\sim23$ mag. To improve selection methods, the medium-band follow-up imaging has been carried out using the SED camera for QUasars in Early uNiverse (SQUEAN) on the Otto Struve 2.1 m Telescope. The optical spectra of the candidates were obtained with 8-m class telescopes. We newly discovered 10 quasars with $-25<M_{1450}<-23$ at $z\sim5$, among which three have been missed in a previous survey using the same optical data over the same area, implying the necessity for improvements in high redshift faint quasars selection. We derived photometric redshifts from the medium-band data, and find that they have high accuracies of $\langle|\Delta z|/(1+z)\rangle=0.016$. The medium-band-based approach allows us to rule out many of the interlopers that contaminate $\gtrsim20~\%$ of the broad-band-selected quasar candidates. These results suggest that the medium-band-based approach is a powerful way to identify $z\sim5$ quasars and measure their redshifts at high accuracy (1-2 %). It is also a cost-effective way to understand the contribution of quasars to the cosmic re-ionization history.

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 on the analysis of 34 years of photometric observations of the pulsating helium atmosphere white dwarf GD358. The complete data set includes archival data from 1982-2006, and 1195.2 hours of new observations from 2007- 2016. From this data set, we extract 15 frequencies representing g-mode pulsation modes, adding 4 modes to the 11 modes known previously. We present evidence that these 15 modes are ell = 1 modes, 13 of which belong to a consecutive sequence in radial overtone k. We perform a detailed asteroseismic analysis using models that include parameterized, complex carbon and oxygen core composition profiles to fit the periods. Recent spectroscopic analyses place GD358 near the red edge of the DBV instability strip, at 24,000 plus or minus 500 K and a log g of 7.8 plus or minus 0.08 dex. The surface gravity translates to a mass range of 0.455 to 0.540 solar masses. Our best fit model has a temperature of 23,650 K and a mass of 0.5706 solar masses. That is slightly more massive than suggested by most the recent spectroscopy. We find a pure helium layer mass of 10^-5.50, consistent with the result of previous studies and the outward diffusion of helium over time.

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.

In hierarchical galaxy formation models, galaxies evolve through mergers and accretions. Tidally-disrupted debris from these processes can remain as diffuse, faint structures, which can provide useful insight into the assembly history of galaxies. To investigate the properties of the faint structures in outskirts of nearby galaxies, we conduct deep and wide-field imaging survey with KMTNet. We present our observing strategy and optimal data reduction process to recover the faint extended features in the imaging data of NGC 1291 taken with KMTNet. Through the dark sky flat-fielding and optimal sky subtraction, we can effectively remove inhomogeneous patterns. In the combined images, the peak-to-peak global sky gradients were reduced to less than $\sim0.5$% and $\sim0.3$% of the original $B$- and $R$-band sky levels, respectively. However, we find local spatial fluctuations in the background sky which can affect the precise measurement of the sky value. Consequently, we can reach the surface brightness of $\mu_{B,1\sigma} \sim$ 29.5 and $\mu_{R,1\sigma}\sim$ 28.5 mag arcsec$^{-2}$ in azimuthally averaged one-dimensional surface brightness profiles, that is mainly limited by the uncertainty in the sky determination. These results suggest that the deep imaging data produced by KMTNet are suitable to study the faint features of nearby galaxies such as outer disks and dwarf companions, but unideal (not impossible) to detect stellar halos. The one-dimensional profile revealed that NGC 1291 appeared to have Type I disk out to $R$ $\sim$ 30 kpc with no obvious color gradient and excess light due to a stellar halo was undetected.

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.