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 present the analysis of microlensing event OGLE-2006-BLG-284, which has a lens system that consists of two stars and a gas giant planet with a mass ratio of $q_p = (1.26\pm 0.19) \times 10^{-3}$ to the primary. The mass ratio of the two stars is $q_s = 0.289\pm 0.011$, and their projected separation is $s_s = 2.1\pm 0.7\,$AU, while the projected separation of the planet from the primary is $s_p = 2.2\pm 0.8\,$AU. For this lens system to have stable orbits, the three-dimensional separation of either the primary and secondary stars or the planet and primary star must be much larger than that these projected separations. Since we do not know which is the case, the system could include either a circumbinary or a circumstellar planet. Because there is no measurement of the microlensing parallax effect or lens system brightness, we can only make a rough Bayesian estimate of the lens system masses and brightness. We find host star and planet masses of $M_{L1} = 0.35^{+0.30}_{-0.20}\,M_\odot$, $M_{L2} = 0.10^{+0.09}_{-0.06}\,M_\odot$, and $m_p = 144^{+126}_{-82}\,M_\oplus$, and the $K$-band magnitude of the combined brightness of the host stars is $K_L = 19.7^{+0.7}_{-1.0}$. The separation between the lens and source system will be $\sim 90\,$mas in mid-2020, so it should be possible to detect the host system with follow-up adaptive optics or Hubble Space Telescope observations.

We present the analysis of the binary-lens microlensing event OGLE-2013-BLG-0911. The best-fit solutions indicate the binary mass ratio of q~0.03 which differs from that reported in Shvartzvald+2016. The event suffers from the well-known close/wide degeneracy, resulting in two groups of solutions for the projected separation normalized by the Einstein radius of s~0.15 or s~7. The finite source and the parallax observations allow us to measure the lens physical parameters. The lens system is an M-dwarf orbited by a massive Jupiter companion at very close (M_{host}=0.30^{+0.08}_{-0.06} M_{Sun}, M_{comp}=10.1^{+2.9}_{-2.2} M_{Jup}, a_{exp}=0.40^{+0.05}_{-0.04} au) or wide (M_{host}=0.28^{+0.10}_{-0.08} M_{Sun}, M_{comp}=9.9^{+3.8}_{-3.5}M_{Jup}, a_{exp}=18.0^{+3.2}_{-3.2} au) separation. Although the mass ratio is slightly above the planet-brown dwarf (BD) mass-ratio boundary of q=0.03 which is generally used, the median physical mass of the companion is slightly below the planet-BD mass boundary of 13M_{Jup}. It is likely that the formation mechanisms for BDs and planets are different and the objects near the boundaries could have been formed by either mechanism. It is important to probe the distribution of such companions with masses of ~13M_{Jup} in order to statistically constrain the formation theories for both BDs and massive planets. In particular, the microlensing method is able to probe the distribution around low-mass M-dwarfs and even BDs which is challenging for other exoplanet detection methods.

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.

We report the discovery of a gas giant planet orbiting a low-mass host star in the microlensing event MOA-bin-29 that occurred in 2006. We find five degenerate solutions with the planet/host-star mass ratio of $q \sim 10^{-2}$. The Einstein radius crossing time of all models are relatively short ($\sim 4-7$ days), which indicates that the mass of host star is likely low. The measured lens-source proper motion is $5-9$ ${\rm mas}\ {\rm yr}^{-1}$ depending on the models. Since only finite source effects are detected, we conduct a Bayesian analysis in order to obtain the posterior probability distribution of the lens physical properties. As a result, we find the lens system is likely to be a gas giant orbiting a brown dwarf or a very late M-dwarf in the Galactic bulge. The probability distributions of the physical parameters for the five degenerate models are consistent within the range of error. By combining these probability distributions, we conclude that the lens system is a gas giant with a mass of $M_{\rm p} = 0.63^{+1.13}_{-0.39}\ M_{\rm Jup}$ orbiting a brown dwarf with a mass of $M_{\rm h} = 0.06^{+0.11}_{-0.04}\ M_\odot$ at a projected star-planet separation of $r_\perp = 0.53^{+0.89}_{-0.18}\ {\rm au}$. The lens distance is $D_{\rm L} = 6.89^{+1.19}_{-1.19}\ {\rm kpc}$, i.e., likely within the Galactic bulge.

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.

We report the discovery of an exoplanet in microlensing event OGLE-2015-BLG-1649. The planet/host-star mass ratio is $q =7.2 \times 10^{-3}$ and the projected separation normalized by the Einstein radius is $s = 0.9$. The upper limit of the lens flux is obtained from adaptive optics observations by IRCS/Subaru, which excludes the probability of a G-dwarf or more massive host star and helps to put a tighter constraint on the lens mass as well as commenting on the formation scenarios of giant planets orbiting low-mass stars. We conduct a Bayesian analysis including constraints on the lens flux to derive the probability distribution of the physical parameters of the lens system. We thereby find that the masses of the host star and planet are $M_{L} = 0.34 \pm 0.19 M_{\odot}$ and $M_{p} = 2.5^{+1.5}_{-1.4} M_{Jup}$, respectively. The distance to the system is $D_{L} = 4.23^{+1.51}_{-1.64}$kpc. The projected star-planet separation is $a_{\perp} = 2.07^{+0.65}_{-0.77}$AU. The lens-source relative proper motion of the event is quite high, at $\sim 7.1 \, {\rm mas/yr}$. Therefore, we may be able to determine the lens physical parameters uniquely or place much stronger constraints on them by measuring the color-dependent centroid shift and/or the image elongation with additional high resolution imaging already a few years from now.

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

We 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 the standard model of cosmology, the Universe is static in comoving coordinates; expansion occurs homogeneously and is represented by a global scale factor. The baryon acoustic oscillation (BAO) peak location is a statistical tracer that represents, in the standard model, a fixed comoving-length standard ruler. Recent gravitational collapse should modify the metric, rendering the effective scale factor, and thus the BAO standard ruler, spatially inhomogeneous. Using the Sloan Digital Sky Survey, we show to high significance (P < 0.001) that the spatial compression of the BAO peak location increases as the spatial paths' overlap with superclusters increases. Detailed observational and theoretical calibration of this BAO peak location environment dependence will be needed when interpreting the next decade's cosmological surveys.

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

We search for correlations between positions of extragalactic objects and arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) with primary energy $E \ge 40$ EeV as observed by the surface detector array of the Telescope Array (TA) experiment during the first 40 months of operation. We examined several public astronomical object catalogs, including the Veron-Cetty and Veron catalog of active galactic nuclei. We counted the number of TA events correlated with objects in each catalog as a function of three parameters: the maximum angular separation between a TA event and an object, the minimum energy of the events, and the maximum redshift of the objects. We determine combinations of these parameters which maximize the correlations, and calculate the chance probabilities of having the same levels of correlations from an isotropic distribution of UHECR arrival directions. No statistically significant correlations are found when penalties for scanning over the above parameters and for searching in several catalogs are taken into account.

Flat space models with multiply connected topology, which have compact dimensions, are tested against the distribution of high-redshift ($z \geq 4$) quasars of the Sloan Digital Sky Survey (SDSS). When the compact dimensions are smaller in size than the observed universe, topological lensing occurs, in which multiple images of single objects (ghost images) are observed. We improve on the recently introduced method to identify ghost images by means of four-point statistics. Our method is valid for any of the 17 multiply connected flat models, including nontoroial ones that are compactified by screw motions or glide reflection. Applying the method to the data revealed one possible case of topological lensing caused by sixth-turn screw motion, however, it is consistent with the simply connected model by this test alone. Moreover, simulations suggest that we cannot exclude the other space models despite the absence of their signatures. This uncertainty mainly originates from the patchy coverage of SDSS in the South Galactic cap, and the situation will be improved by future wide-field spectroscopic surveys.

We measure the spectrum of cosmic rays with energies greater than $10^{18.2}$ eV with the Fluorescence Detectors (FDs) and the Surface Detectors (SDs) of the Telescope Array Experiment using the data taken in our first 2.3-year observation from May 27 2008 to September 7 2010. A hybrid air shower reconstruction technique is employed to improve accuracies in determination of arrival directions and primary energies of cosmic rays using both FD and SD data. The energy spectrum presented here is in agreement with our previously published spectra and the HiRes results.

We present a measurement of the energy spectrum of ultra-high-energy cosmic rays performed by the Telescope Array experiment using monocular observations from its two new FADC-based fluorescence detectors. After a short description of the experiment, we describe the data analysis and event reconstruction procedures. Since the aperture of the experiment must be calculated by Monte Carlo simulation, we describe this calculation and the comparisons of simulated and real data used to verify the validity of the aperture calculation. Finally, we present the energy spectrum calculated from the merged monocular data sets of the two FADC-based detectors, and also the combination of this merged spectrum with an independent, previously published monocular spectrum measurement performed by Telescope Array's third fluorescence detector (Abu-Zayyad {\it et al.}, {Astropart. Phys.} 39 (2012), 109). This combined spectrum corroborates the recently published Telescope Array surface detector spectrum (Abu-Zayyad {\it et al.}, {Astrophys. Journ.} 768 (2013), L1) with independent systematic uncertainties.

The Telescope Array (TA) collaboration has measured the energy spectrum of ultra-high energy cosmic rays with primary energies above 1.6 x 10^(18) eV. This measurement is based upon four years of observation by the surface detector component of TA. The spectrum shows a dip at an energy of 4.6 x 10^(18) eV and a steepening at 5.4 x 10^(19) eV which is consistent with the expectation from the GZK cutoff. We present the results of a technique, new to the analysis of ultra-high energy cosmic ray surface detector data, that involves generating a complete simulation of ultra-high energy cosmic rays striking the TA surface detector. The procedure starts with shower simulations using the CORSIKA Monte Carlo program where we have solved the problems caused by use of the "thinning" approximation. This simulation method allows us to make an accurate calculation of the acceptance of the detector for the energies concerned.

We study the anisotropy of Ultra-High Energy Cosmic Ray (UHECR) events collected by the Telescope Array (TA) detector in the first 40 months of operation. Following earlier studies, we examine event sets with energy thresholds of 10 EeV, 40 EeV, and 57 EeV. We find that the distributions of the events in right ascension and declination are compatible with an isotropic distribution in all three sets. We then compare with previously reported clustering of the UHECR events at small angular scales. No significant clustering is found in the TA data. We then check the events with E>57 EeV for correlations with nearby active galactic nuclei. No significant correlation is found. Finally, we examine all three sets for correlations with the large-scale structure of the Universe. We find that the two higher-energy sets are compatible with both an isotropic distribution and the hypothesis that UHECR sources follow the matter distribution of the Universe (the LSS hypothesis), while the event set with E>10 EeV is compatible with isotropy and is not compatible with the LSS hypothesis at 95% CL unless large deflection angles are also assumed. We show that accounting for UHECR deflections in a realistic model of the Galactic magnetic field can make this set compatible with the LSS hypothesis.

The Telescope Array (TA) experiment, located in the western desert of Utah,USA, is designed for observation of extensive air showers from extremely high energy cosmic rays. The experiment has a surface detector array surrounded by three fluorescence detectors to enable simultaneous detection of shower particles at ground level and fluorescence photons along the shower track. The TA surface detectors and fluorescence detectors started full hybrid observation in March, 2008. In this article we describe the design and technical features of the TA surface detector.

The Telescope Array's Middle Drum fluorescence detector was instrumented with telescopes refurbished from the High Resolution Fly's Eye's HiRes-1 site. The data observed by Middle Drum in monocular mode was analyzed via the HiRes-1 profile-constrained geometry reconstruction technique and utilized the same calibration techniques enabling a direct comparison of the energy spectra and energy scales between the two experiments. The spectrum measured using the Middle Drum telescopes is based on a three-year exposure collected between December 16, 2007 and December 16, 2010. The calculated difference between the spectrum of the Middle Drum observations and the published spectrum obtained by the data collected by the HiRes-1 site allows the HiRes-1 energy scale to be transferred to Middle Drum. The HiRes energy scale is applied to the entire Telescope Array by making a comparison between Middle Drum monocular events and hybrid events that triggered both Middle Drum and the Telescope Array's scintillator Ground Array.

Since 2007, the Telescope Array (TA) experiment, based in Utah, USA, has been observing ultra high energy cosmic rays to understand their origins. The experiment involves a surface detector (SD) array and three fluorescence detector (FD) stations. FD stations, installed surrounding the SD array, measure the air fluorescence light emitted from extensive air showers (EASs) for precise determination of their energies and species. The detectors employed at one of the three FD stations were relocated from the High Resolution Fly's Eye experiment. At the other two stations, newly designed detectors were constructed for the TA experiment. An FD consists of a primary mirror and a camera equipped with photomultiplier tubes. To obtain the EAS parameters with high accuracies, understanding the FD optical characteristics is important. In this paper, we report the characteristics and installation of new FDs and the performances of the FD components. The results of the monitored mirror reflectance during the observation time are also described in this report.

A flat Universe model supported by recent observations has 18 possible choices for its overall topology. To detect or exclude these possibilities is one of the most important tasks in modern cosmology, but it has been very difficult for globally inhomogeneous ones because of a long-time calculation. In this brief paper we provide an object-based 3D method to overcome the problem, as an extension of Fujii & Yoshii (2011a). Though the test depends on the observer's location in the universe, this method drastically reduces calculation times to constrain inhomogeneous topologies, and will be useful in exhaustively constraining the size of the Universe.

We present a method to constrain cosmic topology from the distribution of astronomical objects projected on the celestial sphere. This is an extension of the 3D method introduced in Fujii & Yoshii (2011) that is to search for a pair of pairs of observed objects (quadruplet) linked by a holonomy, i.e., the method we present here is to search for a pair of celestial sphere $n$-tuplets for $n \geq 3$. We find, however, that this method is impractical to apply in realistic situations due to the small signal to noise ratio. We conclude therefore that it is unrealistic to constrain the topology of the Universe from the celestial distribution, and the 3D catalogs are necessary for the purpose.

A new, improved version of a cosmic crystallography method for constraining cosmic topology is introduced. Like the circles-in-the-sky method using CMB data, we work in a thin, shell-like region containing plenty of objects. Two pairs of objects (quadruplet) linked by a holonomy show a specific distribution pattern, and three filters of \emph{separation, vectorial condition}, and \emph{lifetime of objects} extract these quadruplets. Each object $P_i$ is assigned an integer $s_i$, which is the number of candidate quadruplets including $P_i$ as their members. Then an additional device of $s_i$-histogram is used to extract topological ghosts, which tend to have high values of $s_i$. In this paper we consider flat spaces with Euclidean geometry, and the filters are designed to constrain their holonomies. As the second filter, we prepared five types that are specialized for constraining specific holonomies: one for translation, one for half-turn corkscrew motion and glide reflection, and three for $n$-th turn corkscrew motion for $n=4, 3,$ and 6. {Every multiconnected space has holonomies that are detected by at least one of these five filters.} Our method is applied to the catalogs of toy quasars in flat $\Lambda$-CDM universes whose typical sizes correspond to $z\sim 5$. With these simulations our method is found to work quite well. {These are the situations in which type-II pair crystallography methods are insensitive because of the tiny number of ghosts. Moreover, in the flat cases, our method should be more sensitive than the type-I pair (or, in general, $n$-tuplet) methods because of its multifilter construction and its independence from $n$.}

A search for ultra-high energy gamma-rays emitted by the young, nearby supernova remnant that was discovered recently by the COMPTEL and ROSAT satellites was made using the JANZOS database for the period 1987-1993. A 95% confidence upper limit on the flux above 100 TeV of 3x10^(-13)cm^(-2)sec^(-1) was obtained. This is an order of magnitude below the expected flux based on the standard model of cosmic ray acceleration in supernova shocks. An optical survey of the region that has been commenced is also reported. This uses UK and ESO Schmidt plates, and CCD images by a NZ/Japan microlensing group.