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 examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

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.

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted.

Feedback is indispensable in galaxy formation. However, lacking resolutions, cosmological simulations often use ad hoc feedback parameters. Conversely, small-box simulations, while better resolving the feedback, cannot capture gas evolution beyond the simulation domain. We aim to bridge the gap by implementing small-box results of supernovae-driven outflows into dark matter halo-scale simulations and studying their impact on large scales. Galactic outflows are multiphase, but small-box simulations show that the hot phase (T$\approx$ 10$^{6-7}$ K) carries the majority of energy and metals. We implement hot outflows in idealized simulations of the Milky Way halo, and examine how they impact the circumgalactic medium (CGM). In this paper, we discuss the case when the star formation surface density is low and therefore the emerging hot outflows are gravitationally bound by the halo. We find that outflows form a large-scale, metal-enriched atmosphere with fountain motions. As hot gas accumulates, the inner atmosphere becomes "saturated". Cool gas condenses, with a rate balancing the injection of the hot outflows. This balance leads to a universal density profile of the hot atmosphere, independent of mass outflow rate. The atmosphere has a radially-decreasing temperature, naturally producing the observed X-ray luminosity and column densities of O VI, O VII, O VIII. The self-regulated atmosphere has a baryon and a metal mass of $(0.5-1.2)\times 10^{10}M_\odot$ and $(0.6-1.4)\times 10^8 M_\odot$, respectively, small compared to the ``missing" baryons and metals from the halo. We conjecture that the missing materials reside at even larger radii, ejected by more powerful outflows in the past.

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.

The Chern-Simons term, through which the cosmic Axion-like field couples to the electromagnetic field, has the effect to rotate CMB polarization directions and to break the CPT symmetry. This rotation will change the CMB power spectra, no matter isotropic or anisotropic the rotation angle is. In this paper we revisit this issue by further considering the correlations between the (anisotropic) rotation angle $\alpha$ and the CMB temperature and (unrotated) $E$ polarization fields. These correlations could be generated in the Axion-like models with nonzero potential under the adiabatic initial condition. We first investigate how these correlations contribute further modifications to the CMB power spectra, then calculate the CMB bispectra for the temperature and rotated polarization fields. These bispectra would vanish if the $T\alpha$ and $E\alpha$ correlations are absent. So, they are useful in searching for CPT violation and the $T\alpha$ and $E\alpha$ correlations arisen in the Axion-like models.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

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 use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

Partial condensation of dust from the Solar nebula is likely responsible for the diverse chemical compositions of chondrites and rocky planets/planetesimals in the inner Solar system. We present a forward physical-chemical model of a protoplanetary disc to predict the chemical compositions of planetary building blocks that may form from such a disc. Our model includes the physical evolution of the disc and the condensation, partial advection, and decoupling of the dust within it. The chemical composition of the condensate changes with time and radius. We compare the results of two dust condensation models: one where an element condenses when the midplane temperature in the disc is lower than the 50\% condensation temperature ($\rm T_{50}$) of that element and the other where the condensation of the dust is calculated by a Gibbs free energy minimization technique assuming chemical equilibrium at local disc temperature and pressure. The results of two models are generally consistent with some systematic differences of $\sim 10$\% depending upon the radial distance and an element's condensation temperature. Both models predict compositions similar to CM, CO, and CV chondrites provided that the decoupling timescale of the dust is on the order of the evolution timescale of the disc or longer. If the decoupling timescale is too short, the composition deviates significantly from the measured values. These models may contribute to our understanding of the chemical compositions of chondrites, and ultimately the terrestrial planets in the solar system, and may constrain the potential chemical compositions of rocky exoplanets.

In this paper, a new measurement system based on integration method is presented,which can provide all-weather dependability and higher precision for the measurement of FAST's feed support system. The measurement system consists of three types of measuring equipments, and a processing software with the core data fusion algorithm. The Strapdown Inertial Navigation System(SINS) can autonomously measure the position, speed and attitude of the carrier. Its own shortcoming is the measurement data diverges rapidly over time. SINS must combine the Global Positioning System(GPS) and the Total Station(TS) to obtain high-precision measurement data. Kalman filtering algorithm is adopted for the integration measurement system, which is an optimal algorithm to estimate the measurement errors. To evaluate the performance, series of tests are carried out. For the feed cabin, the maximum RMS of the position is 14.56mm, the maximum RMS of the attitude is 0.095{\deg}, these value are less than 15mm and 0.1{\deg} as the precision for measuring the feed cabin. For the Stewart manipulator, the maximum RMS of the position is 2.99mm, the maximum RMS of the at titude is 0.093{\deg}, these value are less than 3mm and 0.1{\deg} as the precision for measuring the Stewart manipulator. As a result, the new measurement meets the requirement of measurement precision for FAST's feed support system.

Asteroid impacts pose a major threat to all life on Earth. The age of the dinosaurs was abruptly ended by a 10-km-diameter asteroid. Currently, a nuclear device is the only means of deflecting large Potentially Hazardous Asteroids (PHAs) away from an Earth-impacting trajectory. The Enhanced Kinetic Impactor (EKI) concept is proposed to deflect large PHAs via maneuvering space rocks. First, an unmanned spacecraft is launched to rendezvous with an intermediate Near-Earth Asteroid (NEA). Then, more than one hundred tons of rocks are collected from the NEA as the EKI. The NEA can also be captured as the EKI if the NEA is very small. Finally, the EKI is maneuvered to impact the PHA at a high speed, resulting in a significant deflection of the PHA. For example, to deflect Apophis, as much as 200 t of rocks could be collected from a NEA as the EKI based on existing engineering capabilities. The EKI can produce a velocity increment (delta-v) of 39.81 mm/s in Apophis, thereby increasing the minimum geocentric distance during the close encounter in 2029 by 1,866.93 km. This mission can be completed in 3.96 years with a propellant cost of 2.98 t. Compared with a classic kinetic impactor, the deflection distance can be increased one order of magnitude. The EKI concept breaks through the limitation of the ground-based launch capability, which can significantly increase the mass of the impactor. We anticipate that our research will be a starting point for efficient planetary defense against large PHAs.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

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 on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

Supernovae (SNe) drive multiphase galactic outflows, impacting galaxy formation; however, cosmological simulations mostly use \textit{ad hoc} feedback models for outflows, making outflow-related predictions from first principles problematic. Recent small-box simulations resolve individual SNe remnants in the interstellar medium (ISM), naturally driving outflows and permitting a determination of the wind loading factors of energy $\eta_E$, mass $\eta_m$, and metals $\eta_Z$. In this Letter, we compile small-box results, and find consensus that the hot outflows are much more powerful than the cool outflows: (i) hot outflows generally dominate the energy flux, and (ii) their specific energy $e_{s,h}$ is 10-1000 times higher than cool outflows. Moreover, the properties of hot outflows are remarkably simple: $e_{s,h} \propto \eta_{E,h}/\eta_{m,h}$ is almost invariant over four orders of magnitude of star formation surface density. Also, we find tentatively that $\eta_{E,h}/\eta_{Z,h} \sim$ 0.5. If corroborated by more simulation data, these correlations reduce the three hot phase loading factors into one. Finally, this one parameter is closely related to whether the ISM has a "breakout" condition. The narrow range of \esh\ indicates that hot outflows cannot escape dark matter halos with log $M_{\rm{halo}}\ [M_\odot] \gtrsim 12$. This mass is also where the galaxy mass-metallicity relation reaches its plateau, implying a deep connection between \textit{hot} outflows and galaxy formation. We argue that hot outflows should be included explicitly in cosmological simulations and (semi-)analytic modeling of galaxy formation.

The circumgalactic medium (CGM) is closely linked to galaxy formation and evolution, but difficult to characterize observationally and typically poorly resolved in cosmological simulations. We use spherically-symmetric, idealized, high-resolution simulations of the CGM in $10^{12}M_\odot$ and $10^{11}M_\odot$ dark matter halos to characterize the gas pressure, turbulent and radial velocities, and degree of thermal and effective dynamic pressure support in the overall CGM as well as in its high- and low-temperature phases. We find that the $10^{12}M_\odot$ halo contains a CGM mostly formed of a hot gas halo in hydrostatic equilibrium out of which cold gas condenses and falls onto the central galaxy, while the $10^{11}M_\odot$ halo's CGM is not in hydrostatic equilibrium, has a wider spread of properties at a given galactocentric radius, does not have a clear separation of hot and cold phases, and is dominated by bulk motions. We also find that the degree of pressure support in the $10^{11}M_\odot$ halo is strongly dependent on the parameters of the galactic winds of the central galaxy. These results promote the idea that there is no "average" CGM and care must be taken when setting the initial conditions for a small-box simulation of a patch of the CGM.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

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.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

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.

Galactic outflows produced by stellar feedback are known to be multiphase in nature. Both observations and simulations indicate that the material within several kpc of galactic disk mid-planes consists of warm clouds embedded within a hot wind. A theoretical understanding of the outflow phenomenon, including both winds and fountain flows, requires study of the interactions among thermal phases. We develop a method to quantify these interactions via measurements of mass, momentum, and energy flux exchanges using temporally and spatially averaged quantities and conservation laws. We apply this method to a star-forming ISM MHD simulation based on the TIGRESS framework, for Solar neighbourhood conditions. To evaluate the extent of interactions among the phases, we first examine the validity of the ``ballistic model,'' which predicts trajectories of the warm phase ($5050\,\rm{K}<T<2\times10^4\,\rm{K}$) treated as non-interacting clouds. This model is successful at intermediate vertical velocities ($ 50$ km s$^{-1}$ $\lesssim |v_z| \lesssim 100 $ km s$^{-1}$), but at higher velocities we observe an excess in simulated warm outflow compared to the ballistic model. This discrepancy cannot be fully accounted for by cooling of high-velocity intermediate-temperature ($2\times10^4\,\rm{K}<T<5\times10^5\,\rm{K}$) gas. By examining the fluxes of mass, momentum and energy, we conclude that warm phase gains mass via cooling of the intermediate phase, while momentum transfer occurs from the hot ($T>5\times10^5\,\rm{K}$) to the warm phase. The large energy flux from the hot outflow that is transferred to the warm and intermediate phases is quickly radiated away. A simple interaction model implies an effective warm cloud size in the fountain flow of a few 100~pc, showing that warm-hot flux exchange mainly involves a few large clouds rather than many small ones.

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

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.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The low energy (LE) X-ray telescope is one of the three main instruments of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT). It is equipped with Swept Charge Device (SCD) sensor arrays with a total geometrical area of 384 cm2 and an energy band from 0.7 keV to 13 keV. In order to evaluate the particle induced X-ray background and the cosmic X-ray background simultaneously, LE adopts collimators to define four types of Field Of Views (FOVs). LE is constituted of three detector boxes (LEDs) and an electric control box (LEB) and achieves a good energy resolution of 140 eV at 5.9 keV, an excellent time resolution of 0.98 ms, as well as an extremely low pileup (<1% at 18000 cts/s). Detailed performance tests and calibration on the ground have been performed, including energy-channel relation, energy response, detection efficiency and time response.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

Atmospheric muons are one of the main backgrounds for current Water- and Ice-Cherenkov neutrino telescopes designed to detect astrophysical neutrinos. The inclusive fluxes of atmospheric muons and neutrinos from hadronic interactions of cosmic rays have been extensively studied with Monte Carlo and cascade equation methods, for example, CORSIKA and MCEq. However, the muons that are pair produced in electromagnetic interaction of high energy photons are quantitatively not well understood. We present new simulation results and assess the model dependencies of the high-energy atmospheric muon flux including those from electromagnetic interactions, using a new numerical electromagnetic cascade equation solver EmCa that can be easily coupled with the hadronic solver MCEq. Both codes are in active development with the particular aim to become part of the next generation CORSIKA 8 air shower simulation package. The combination of EmCa and MCEq accounts for material effects that have not been previously included in most of the available codes. Hence, the influence of these effects on the air showers will also be briefly discussed.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

The Medium Energy X-ray telescope (ME) is one of the three main telescopes on board the Insight Hard X-ray Modulation Telescope (Insight-HXMT) astronomy satellite. ME contains 1728 pixels of Si-PIN detectors sensitive in 5-30 keV with a total geometrical area of 952 cm2. Application Specific Integrated Circuit (ASIC) chips, VA32TA6, is used to achieve low power consumption and low readout noise. The collimators define three kinds of field of views (FOVs) for the telescope, 1{\deg}{\times}4{\deg}, 4{\deg}{\times}4{\deg}, and blocked ones. Combination of such FOVs can be used to estimate the in-orbit X-ray and particle background components. The energy resolution of ME is ~3 keV at 17.8 keV (FWHM) and the time resolution is 255 {\mu}s. In this paper, we introduce the design and performance of ME.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

Cosmological CPT violation will rotate the polarized direction of CMB photons, convert partial CMB E mode into B mode and vice versa. It will generate non-zero EB, TB spectra and change the EE, BB, TE spectra. This phenomenon gives us a way to detect the CPT-violation signature from CMB observations, and also provides a new mechanism to produce B mode polarization. In this paper, we perform a global analysis on tensor-to-scalar ratio $r$ and polarization rotation angles based on current CMB datasets with both low $\ell$ (Planck, BICEP2/Keck Array) and high $\ell$ (POLARBEAR, SPTpol, ACTPol). Benefited from the high precision of CMB data, we obtain the isotropic rotation angle $\bar{\alpha} = -0.01^\circ \pm 0.37^\circ $ at 68% C.L., the variance of the anisotropic rotation angles $C^{\alpha}(0)<0.0032\,\mathrm{rad}^2$, the scale invariant power spectrum $D^{\alpha\alpha}_{\ell \in [2, 350]}<4.71\times 10^{-5} \,\mathrm{rad}^2$ and $r<0.057$ at 95% C.L.. Our result shows that with the polarization rotation effect, the 95% upper limit on $r$ gets tightened by 17%.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

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 on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

In this paper, three ensemble methods: Random Forest, XGBoost, and a Hybrid Ensemble method were implemented to classify imbalanced pulsar candidates. To assist these methods, tree models were used to select features among 30 features of pulsar candidates from references. The skewness of the integrated pulse profile, chi-squared value for sine-squared fit to amended profile and best S/N value play important roles in Random Forest, while the skewness of the integrated pulse profile is one of the most significant features in XGBoost. More than 20 features were selected by their relative scores and then applied in three ensemble methods. In the Hybrid Ensemble method, we combined Random Forest and XGBoost with EasyEnsemble. By changing thresholds, we tried to make a trade-off between Recall and Precision to make them approximately equal and as high as possible. Experiments on HTRU 1 and HTRU 2 datasets show that the Hybrid Ensemble method achieves higher Recall than the other two algorithms. In HTRU 1 dataset, Recall, Precision, and F-Score of the Hybrid Ensemble method are $0.967$, $0.971$, and $0.969$, respectively. In HTRU 2 dataset, the three values of that are $0.920$, $0.917$, and $0.918$, respectively.

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

As performance of dedicated facilities continually improved, massive pulsar candidates are being received, which makes selecting valuable pulsar signals from candidates challenging. In this paper, we designed a deep convolutional neural network (CNN) with 11 layers for classifying pulsar candidates. Compared to artificial designed features, CNN chose sub-integrations plot and sub-bands plot in each candidate as inputs without carrying biases. To address the imbalanced problem, data augmentation method based on synthetic minority samples is proposed according to characteristics of pulsars. The maximum pulses of pulsar candidates were first translated to the same position, then new samples were generated by adding up multiple subplots of pulsars. The data augmentation method is simple and effective for obtaining varied and representative samples which keep pulsar characteristics. In the experiments on HTRU 1 dataset, it shows that this model can achieve recall as 0.962 while precision as 0.963.

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

The tomographic Alcock-Paczynski (AP) method utilizes the redshift evolution of the AP distortion to place constraints on cosmological parameters. It has proved to be a robust method that can separate the AP signature from the redshift space distortion (RSD) effect, and deliver powerful cosmological constraints using the $\lesssim 40h^{-1}\ \rm Mpc$ clustering region. In previous works, the tomographic AP method was performed via the anisotropic 2-point correlation function statistic. In this work we consider the feasibility of conducting the analysis in the Fourier domain and examine the pros and cons of this approach. We use the integrated galaxy power spectrum (PS) as a function of direction, $\hat P_{\Delta k}(\mu)$, to quantify the magnitude of anisotropy in the large-scale structure clustering, and use its redshift variation to do the AP test. The method is tested on the large, high resolution Big-MultiDark Planck (BigMD) simulation at redshifts $z=0-1$, using the underlying true cosmology $\Omega_m=0.3071,\ w=-1$. Testing the redshift evolution of $\hat P_{\Delta k}(\mu)$ in the true cosmology and cosmologies deviating from the truth with $\delta \Omega_m=0.1,\ \delta w=0.3$, we find that the redshift evolution of the AP distortion overwhelms the effects created by the RSD by a factor of $\sim1.7-3.6$. We test the method in the range of $k\in(0.2,1.8)\ h\ \rm Mpc^{-1}$, and find that it works well throughout the entire regime. We tune the halo mass within the range $2\times 10^{13}$ to $10^{14}\ M_{\odot}$, and find that the change of halo bias results in $\lesssim 5 \%$ change in $\hat P_{\Delta k}(\mu)$, which is less significant compared with the cosmological effect. Our work shows that it is feasible to conduct the tomographic AP analysis in the Fourier space.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), 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 selected 90 massive star-forming clumps with strong N2H+, HCO+, HCN, and HNC emission from the Millimetre Astronomy Legacy Team 90 GHz survey. We obtained Herschel data for all 90 sources and NRAO VLA Sky Survey data for 51 of them. We convolved and regridded all images to the same resolution and pixel size and derived the temperature, H2 column density, molecules' abundances and abundance, and ratios of each pixel. Our analysis yields three main conclusions. First, the abundances of N2H+, HCO+, HCN, and HNC increase when the column density decreases and the temperature increases, with spatial variations in their abundances dominated by changes in the H2 column density. Second, the abundance ratios between N2H+, HCO+, HCN, and HNC also display systemic variations as a function of the column density due to the chemical properties of these molecules. Third, the sources associated with the 20 cm continuum emission can be classified into four types based on the behavior of the abundances of the four molecules considered here as a function of this emission. The variations of the first three types could also be attributed to the variation of the H2 column density.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

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.

The analysis of spatial observations on a sphere is important in areas such as geosciences, physics and embryo research, just to name a few. The purpose of the package rcosmo is to conduct efficient information processing, visualisation, manipulation and spatial statistical analysis of Cosmic Microwave Background (CMB) radiation and other spherical data. The package was developed for spherical data stored in the Hierarchical Equal Area isoLatitude Pixelation (Healpix) representation. rcosmo has more than 100 different functions. Most of them initially were developed for CMB, but also can be used for other spherical data as rcosmo contains tools for transforming spherical data in cartesian and geographic coordinates into the HEALPix representation. We give a general description of the package and illustrate some important functionalities and benchmarks.

Neutrino annihilation process on a hyperaccreting disk is one of the leading models to explain the generation of relativistic jets of gamma-ray bursts (GRBs). However, there is still uncertainty regarding the NAE of neutron star-accretion disc (NS-disc) system because of complicated microphysics processes and effects of strong magnetic field. In this paper, we investigate the NAE by assuming that the prompt jet of GRB 070110 is driven by neutrino pair annihilation in the NS-disc system. Our calculation shows $\eta_{\nu\bar{\nu}}> 1.2 \times 10^{-3}$ under the estimated accretion rate $\dot{M}\simeq 0.04\rm M_{\odot}\cdot s^{-1}$. Independent of the detailed accretion disc models, our result shows that the magnetic field may play an important role in the neutrino annihilation process on the hyperaccreting magnetized accretion disc. Compared with the theoretical value of $\eta_{\rm \nu \bar{\nu}}$ of the nonmagnetized BH-disc system, the NAE should increase significantly in the case of the NS-disc system if the GRB is powered by magnetar-disc system.

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.

We find an interesting fact that fictitious retrograde co-orbitals of Saturn, or small bodies inside the retrograde 1:1 resonance with Saturn, are highly unstable in our numerical simulations. It is shown that in the presence of Jupiter, the retrograde co-orbitals will get ejected from Saturn's co-orbital space within a timescale of 10 Myr. This scenario reminds us of the instability of Saturn Trojans caused by both the Great Inequality and the secular resonances. Therefore, we carry out in-depth inspections on both mechanisms and prove that the retrograde resonance overlap, raised by Great Inequality, cannot serve as an explanation for the instability of retrograde co-orbitals, due to the weakness of the retrograde 2:5 resonance with Jupiter at a low eccentricity. However, we discover that both $\nu_5$ and $\nu_6$ secular resonances contribute to the slow growth of the eccentricity, therefore, are possibly the primary causes of the instability inside Saturn's retrograde co-orbital space.

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 analysis of stellar binary microlensing event OGLE-2015-BLG-0060 based on observations obtained from 13 different telescopes. Intensive coverage of the anomalous parts of the light curve was achieved by automated follow-up observations from the robotic telescopes of the Las Cumbres Observatory. We show that, for the first time, all main features of an anomalous microlensing event are well covered by follow-up data, allowing us to estimate the physical parameters of the lens. The strong detection of second-order effects in the event light curve necessitates the inclusion of longer-baseline survey data in order to constrain the parallax vector. We find that the event was most likely caused by a stellar binary-lens with masses $M_{\star1} = 0.87 \pm 0.12 M_{\odot}$ and $M_{\star2} = 0.77 \pm 0.11 M_{\odot}$. The distance to the lensing system is 6.41 $\pm 0.14$ kpc and the projected separation between the two components is 13.85 $\pm 0.16$ AU. Alternative interpretations are also considered.

We use the Richardson-Lucy deconvolution algorithm to extract one dimensional (1D) spectra from LAMOST spectrum images. Compared with other deconvolution algorithms, this algorithm is much more fast. The practice on a real LAMOST image illustrates that the 1D resulting spectrum of this method has a higher SNR and resolution than those extracted by the LAMOST pipeline. Furthermore, our algorithm can effectively depress the ringings that are often shown in the 1D resulting spectra of other deconvolution methods.

Gravitational softening length is one of the key parameters to properly set up a cosmological $N$-body simulation. In this paper, we perform a large suit of high-resolution $N$-body simulations to revise the optimal softening scheme proposed by Power et al. (P03). Our finding is that P03 optimal scheme works well but is over conservative. Using smaller softening lengths than that of P03 can achieve higher spatial resolution and numerically convergent results on both circular velocity and density profiles. However using an over small softening length overpredicts matter density at the inner most region of dark matter haloes. We empirically explore a better optimal softening scheme based on P03 form and find that a small modification works well. This work will be useful for setting up cosmological simulations.

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.

The measurement of the speed of gravitational waves (GWs) is useful to distinguish general relativity from massive gravity. We propose a new model-independent strategy to measure the speed of GWs with the distorted pulsars. Theoretically, when the the GW frequencies from a distorted pulsar are twice the frequencies of EM pulses, they should be emitted at the same time. By measuring the arrival times of these two signals emitted at the same time, the speed of GWs can be calculated with the time difference. Specifically, when the glitches of pulsars are taken into consideration, some pulsars in our Galaxy and nearby galaxies are potential to test our new strategy at high accuracy in the foreseeable future. On the other hand, the new method is meaningful as a motivation for future design of instruments.

The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets' masses, obstructing detection of Earth twins. Since 2015 July we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-velocity measurements and cross-correlation function (CCF) parameters obtained in the first 3 years of observation, interpreting them in the context of spatially-resolved solar observations. We describe a Bayesian mixture-model approach to automated data-quality monitoring. We provide dynamical and daily differential-extinction corrections to place the radial velocities in the heliocentric reference frame, and the CCF shape parameters in the sidereal frame. We achieve a photon-noise limited radial-velocity precision better than 0.43 m s$^{-1}$ per 5-minute observation. The day-to-day precision is limited by zero-point calibration uncertainty with an RMS scatter of about 0.4 m s$^{-1}$. We find significant signals from granulation and solar activity. Within a day, granulation noise dominates, with an amplitude of about 0.4 m s$^{-1}$ and an autocorrelation half-life of 15 minutes. On longer timescales, activity dominates. Sunspot groups broaden the CCF as they cross the solar disc. Facular regions temporarily reduce the intrinsic asymmetry of the CCF. The radial-velocity increase that accompanies an active-region passage has a typical amplitude of 5 m s$^{-1}$ and is correlated with the line asymmetry, but leads it by 3 days. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity.

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

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

The tomographic Alcock-Paczynski (AP) test is a robust large-scale structure (LSS) measurement that receives little contamination from the redshift space distortion (RSD). It has placed tight cosmological constraints by using small and intermediate clustering scales of the LSS data. However, previous works have neglected the cross-correlation among different redshift bins, which could cause the statistical uncertainty being underestimated by $\sim$20\%. In this work, we further improve this method by including this multi-redshifts full correlation. We apply it to the SDSS DR12 galaxies sample and find out that, for $\Lambda$CDM, the combination of AP with the Planck+BAO dataset slightly reduces (within 1-$\sigma$) $\Omega_m$ to $0.304\pm0.007$ (68.3\% CL). This then leads to a larger $H_0$ and also mildly affects $\Omega_b h^2$, $n_s$ and the derived parameters $z_*$, $r_*$, $z_{re}$ but not $\tau$, $A_s$ and $\sigma_8$. For the flat $w$CDM model, our measurement gives $\Omega_m=0.301\pm 0.010$ and $w=-1.090\pm 0.047$, where the additional AP measurement reduces the error budget by $\sim 25\%$. When including more parameters into the analysis, the AP method also improves the constraints on $\Omega_k$, $\sum m_\mu$, $N_{\rm eff}$ by $20-30\%$. Early universe parameters such as $dn_s/d{\rm ln}k$ and $r$, however, are unaffected. Assuming the dark energy equation of state $w=w_0+w_a \frac{z}{1+z}$, the Planck+BAO+SNIa+$H_0$+AP datasets prefer a dynamical dark energy at $\approx1.5 \sigma$ CL. Finally, we forecast the cosmological constraints expected from the DESI galaxy survey and find that combining AP with CMB+BAO method would improve the $w_0$-$w_a$ constraint by a factor of $\sim 10$.

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

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

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

The interaction between inflowing gas clouds and galactic outflows at the interface where the galactic disk transitions into the circumgalactic medium is an important process in galaxy fueling, yet remains poorly understood. Using a series of tall-box hydrodynamic ENZO simulations, we have studied the interaction between smooth gas inflow and supernovae-driven outflow at the disk-halo interface with pc-scale resolution. A realistic wind of outflowing material is generated by supernovae explosions in the disk, while inflowing gas is injected at the top boundary of the simulation box with an injection velocity ranging from $10-100 \rm \ km \ s^{-1}$. We find that cooling and hydrodynamic instabilities drive the injected gas to fragment into cold ($\sim 10^{3}$ K) cloud clumps with typical densities of $\sim 1 \rm \ cm^{-3}$. These clumps initially accelerate before interacting and partially mixing with the outflow and decelerating to velocities in the 50-100 $\rm km \ s^{-1}$ range. When the gas clumps hit the disk, $10\%-50 \%$ of the injected material is able to accrete (depending on the injection velocity). Clumps originating from gas injected with a higher initial velocity approach the disk with greater ram pressure, allowing them to penetrate through the disk in low density regions. We use (equilibrium) CLOUDY photoionization models to generate absorption and emission signatures of gas accretion, finding that our mock HI and H$\alpha$ observables are prominent and generally consistent with measurements in the Milky Way. We do not predict enhanced emission/absorption for higher ionization states such as OVI.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

The concordance $\Lambda$CDM cosmology predicts tens of satellite galaxies distributed in the inner region ($<40\ \lkpc$) of the Milky Way (MW), yet at most $12$ were discovered at present day, including 3 discovered very recently by Dark Energy Survey (DES) and 5 from other surveys (e.g. PanSTARRS, MagLiteS). We use $5$ ultra-high resolution simulations of MW sized dark matter haloes from the $\aqua$ project, combined with $\galf$ semi-analytical galaxy formation model, to investigate properties of the model satellite galaxy population inside $40\ \lkpc$ of MW sized haloes. On average, in each halo this model predicts about $20$ inner satellite galaxies, among them $5$ are comparable to the classic satellites in the luminosity, these are in stark contrast to the corresponding numbers in observations. We further investigate the survivability of these model inner satellites in the presence of a central stellar disk with a set of ideal simulations. These are done by re-evolving a quarter (30) of the whole $\aqua$ inner satellite galaxies (121) by including a static disk potential in addition to the MW halo. Our finding is that the additional disk completely disrupt 40 percents of these satellites and results in $14$ satellite galaxies within the $40\ \lkpc$ of each $\aqua$ at the end, in reasonably well agreement with observations.

The gravitational instabilities are important to the evolution of the disks and the planet formation in the disks. We calculate the evolution of the disks which form from the collapse of the molecular cloud cores. By changing the properties of the cloud cores and the hydrodynamical viscosity parameters, we explore their effects on the properties of the gravitational instabilities. We find that the disk is unstable when the angular velocity of the molecular cloud core is larger than a critical value. The time duration of the instability increases as the angular velocity of the core increases. The increase of the hydrodynamical viscosity parameter hardly affects the stability of the disk, but decreases the time duration of the critical state of the gravitational instability in the disk. The instability of the disks can happen at very early time of evolution of the disk, which is consistent with the observations.

We present an analysis of the microlensing event OGLE-2015-BLG-0232. This event is challenging to characterize for two reasons. First, the light curve is not well sampled during the caustic crossing due to the proximity of the full Moon impacting the photometry quality. Moreover, the source brightness is difficult to estimate because this event is blended with a nearby K dwarf star. We found that the light curve deviations are likely due to a close brown dwarf companion (i.e., s = 0.55 and q = 0.06), but the exact nature of the lens is still unknown. We finally discuss the potential of follow-up observations to estimate the lens mass and distance in the future.

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.

For the first time in human history, technologies have matured sufficiently to enable a mission capable of discovering and characterizing habitable planets like Earth orbiting sunlike stars other than the Sun. At the same time, such a platform would enable unique science not possible from ground-based facilities. This science is broad and exciting, ranging from new investigations of our own solar system to a full range of astrophysics disciplines. The Habitable Exoplanet Observatory, or HabEx, is one of four studies currently being undertaken by NASA in preparation for the 2020 Astrophysics Decadal Survey. HabEx has been designed to be the Great Observatory of the 2030s, with community involvement through a competed and funded Guest Observer (GO) program. This interim report describes the HabEx baseline concept, which is a space-based 4-meter diameter telescope mission concept with ultraviolet (UV), optical, and near-infrared (near-IR) imaging and spectroscopy capabilities. More information on HabEx can be found at https://www.jpl.nasa.gov/habex

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

We propose the formation of massive pristine dark-matter (DM) halos with masses of $\sim 10^8~M_\odot$, due to the dynamical effects of frequent mergers in rare regions of the Universe with high baryonic streaming velocity relative to DM. Since the streaming motion prevents gas collapse into DM halos and delays prior star formation episodes, the gas remains metal-free until the halo virial temperatures $\gtrsim 2\times 10^4~{\rm K}$. The minimum cooling mass of DM halos is boosted by a factor of $\sim 10-30$ because frequent major mergers of halos further inhibit gas collapse. We use Monte Carlo merger trees to simulate the DM assembly history under a streaming velocity of twice the root-mean-square value, and estimate the number density of massive DM halos containing pristine gas as $\simeq 10^{-4}~{\rm cMpc}^{-3}$. When the gas infall begins, efficient Ly$\alpha$ cooling drives cold streams penetrating inside the halo and feeding a central galactic disk. When one stream collides with the disk, strong shock forms a dense and hot gas cloud, where the gas never forms H$_2$ molecules due to effective collisional dissociation. As a result, a massive gas cloud forms by gravitational instability and collapses directly into a massive black hole (BH) with $M_\bullet \sim 10^5~M_\odot$. Almost simultaneously, a galaxy with $M_{\star, \rm tot}\sim 10^6~M_\odot$ composed of Population III stars forms in the nuclear region. If the typical stellar mass is as high as $\sim 100~M_\odot$, the galaxy could be detected with the James Webb Space Telescope even at $z\gtrsim 15$. These massive seed BHs would be fed by continuous gas accretion from the host galaxy, and grow to be bright quasars observed at $z\gtrsim 6$.

As the discoveries of more minor bodies in retrograde resonances with giant planets, such as 2015 BZ509 and 2006 RJ2, our curiosity about the Kozai-Lidov dynamics inside the retrograde resonance has been sparked. In this study, we focus on the 3D retrograde resonance problem and investigate how the resonant dynamics of a minor body impacts on its own Kozai-Lidov cycle. Firstly we deduce the action-angle variables and canonical transformations that deal with the retrograde orbit specifically. After obtaining the dominant Hamiltonian of this problem, we then carry out the numerical averaging process in closed form to generate phase-space portraits on a $e-\omega$ space. The retrograde 1:1 resonance is particularly scrutinized in detail, and numerical results from a CRTBP model shows a great agreement with the our semi-analytical portraits. On this basis, we inspect two real minor bodies currently trapped in retrograde 1:1 mean motion resonance. It is shown that they have different Kozai-Lidov states, which can be used to analyze the stability of their unique resonances. In the end, we further inspect the Kozai-Lidov dynamics inside the 2:1 and 2:5 retrograde resonance, and find distinct dynamical bifurcations of equilibrium points on phase-space portraits.

In this paper, we use two hydrodynamic simulations to study the barycentre offsets $r_{{\rm off}}$ between the entire halo and the gas or dark matter subhalo with largest overdensity. We restrict our analysis to those halos with large rescaled offsets, i.e. $r_{{\rm off}}/r_{200}\gtrsim 0.05$. The halos which are less massive and at higher redshifts are more likely to have larger $r_{{\rm off}}/r_{200}$. The 3D velocity dispersion of the subhalos with largest overdensity is found to have a similar mass and redshift dependence as $r_{{\rm off}}/r_{200}$. We also find that the gas subhalos with maximum overdensity and larger gas velocity dispersion $\sigma_{v}^{\rm ICM,sub}/\sigma_{200}$ tend to have smaller barycentre offsets $r_{{\rm off}}/r_{200}$. Similar result is found for the highest overdensity dark matter subhalos but with less significance. This anti-correlation is more significant for the snapshot $z=0$ than $z=0.2$ and $0.5$. The underlying physical mechanism for this anti-correlation remains to be explored.

A sample of 542 eclipsing binaries (EBs) with periods shorter than 2 d were selected from the Microlensing Observations in Astrophysics (MOA) EB catalogue (Li et al. 2017) for eclipse-time variation analysis. For this sample we were able to obtain the time series from MOA-II that span 9.5yr. We discovered 91 EBs, out of the 542 EBs, with detected light-travel-time effect signals suggesting the presence of tertiary companions of orbiting periods from 250 d-28 yr. The frequency of EBs with tertiary companions in our sample increases as the period decreases and reaches a value of 0.65 for contact binaries with periods shorter than 0.3 d. If only the contact binaries of periods < 0.26d are considered, the frequency even goes to the unit. Our results suggest that contact binaries with periods close to the 0.22-d contact binary limit are commonly accompanied by relatively close tertiary companions.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

The cosmic distance duality relation (DDR), which connects the angular diameter distance and luminosity distance through a simple formula $D_A(z)(1+z)^2/D_L(z)\equiv1$, is an important relation in cosmology. Therefore, testing the validity of DDR is of great importance. In this paper, we test the possible violation of DDR using the available local data including type Ia supernovae (SNe Ia), galaxy clusters and baryon acoustic oscillations (BAO). We write the modified DDR as $D_A(z)(1+z)^2/D_L(z)=\eta(z)$, and consider two different parameterizations of $\eta(z)$, namely $\eta_1(z)=1+\eta_0 z$ and $\eta_2(z)=1+\eta_0 z/(1+z)$. The luminosity distance from SNe Ia are compared with the angular diameter distance from galaxy clusters and BAO at the same redshift. Two different cluster data are used here, i.e. elliptical clusters and spherical clusters. The parameter $\eta_0$ is obtained using the Markov chain Monte Carlo methods. It is found that $\eta_0$ can be strictly constrained by the elliptical clusters + BAO data, with the best-fitting values $\eta_0=-0.04\pm 0.12$ and $\eta_0=-0.05\pm 0.22$ for the first and second parametrizations, respectively. However, the spherical clusters + BAO data couldn't strictly constrain $\eta_0$ due to the large intrinsic scatter. In any case studied here, no evidence for the violation of DDR is found.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

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

Aims. 2015 BZ509 is the first asteroid confirmed to be in retrograde co-orbit resonance (or 1/-1 resonance) with the giant planets in the solar system. While Saturn is the only giant planet whose Trojans are not discovered until now, we identify some small bodies among Centaurs and Damocloids that are potentially in 1/-1 resonance with Saturn in the present study. Methods. We integrate numerically the motion of the 1000 clones (include the nominal orbit) of each Centaur whose orbit has a semi-major axis between 9.3 au and 9.8 au and an inclination i > 90 deg. To confirm and evaluate the 1/-1 resonant configurations mentioned above, we introduce a useful one-degree integrable approximation for planar 1/-1 resonance. Results. We identify four candidates potentially in 1/-1 resonance with Saturn. The capture in this particular resonant state during the 40000 yr integration timespan is very common for 2006 RJ2 (906/1000 clones), 2006 BZ8 (878/1000 clones), and 2017 SV13 (998/1000 clones), and it is less likely for 2012 YE8 (426/1000 clones). According to our statistical results, 2006 RJ2 is the best candidate to be currently in a 1/-1 mean motion resonance with Saturn, and 2017 SV13 is another important potential candidate. Moreover, 2012 YE8 and 2006 BZ8 are also Centaurs of interest but their current and long-term 1/-1 resonant state with Saturn is less likely. The proportions of the clones captured in the relative long-term stable co-orbit resonance (over 10000 yr) are also given. Conclusions. Small bodies in retrograde co-orbit resonance with giant planets may be more common than previously expected. Identification of these potential mysterious minor bodies encourages the search for such objects on a larger scale in our solar system. The findings of this paper are also useful for understanding the origin and dynamical evolution of the Centaurs and Damocloids on retrograde orbits.

We study the gravitational wave (GW) production induced by the asymmetric jets of gamma-ray bursts (GRBs). The asymmetric jets result in a recoil force acted on the central compact object, whose motion leads to emission of GW. Under reasonable assumptions and simplifications, we derive the analytic form of the produce GWs. The amplitude of emitted GWs is estimated to be relatively low, but possibility exists that they can be detected by future experiments such as the Einstein Telescope. We find the dynamical properties of the central object, which is difficult to be studied via the electromagnetic (EW) channel, can be inferred by measuring the emitted GWs. Moreover, we find the emitted GWs can be used determine whether the relativistic jets is launched by the neutrino annihilation process or the Blandford-Znajek process, which cannot be clearly distinguished by the current GRB observations. Our work manifests the importance of the GW channel in multi-messenger astronomy. The physical information encoded in the GW and EW emissions of an astrophysical object is complementary to each other; in case some physics can not be effectively investigated using the EW channel alone, including the GW channel can be very helpful.

We present the analysis of planetary microlensing event MOA-2011-BLG-291, which has a mass ratio of $q=(3.8\pm0.7)\times10^{-4}$ and a source star that is redder (or brighter) than the bulge main sequence. This event is located at a low Galactic latitude in the survey area that is currently planned for NASA's WFIRST exoplanet microlensing survey. This unusual color for a microlensed source star implies that we cannot assume that the source star is in the Galactic bulge. The favored interpretation is that the source star is a lower main sequence star at a distance of $D_S=4.9\pm1.3\,$kpc in the Galactic disk. However, the source could also be a turn-off star on the far side of the bulge or a sub-giant in the far side of the Galactic disk if it experiences significantly more reddening than the bulge red clump stars. However, these possibilities have only a small effect on our mass estimates for the host star and planet. We find host star and planet masses of $M_{\rm host} =0.15^{+0.27}_{-0.10}M_\odot$ and $m_p=18^{+34}_{-12}M_\oplus$ from a Bayesian analysis with a standard Galactic model under the assumption that the planet hosting probability does not depend on the host mass or distance. However, if we attempt to measure the host and planet masses with host star brightness measurements from high angular resolution follow-up imaging, the implied masses will be sensitive to the host star distance. The WFIRST exoplanet microlensing survey is expected to use this method to determine the masses for many of the planetary systems that it discovers, so this issue has important design implications for the WFIRST exoplanet microlensing survey.

OGLE-2014-BLG-0962 (OB140962) is a stellar binary microlensing event that was well-covered by observations from the Spitzer satellite as well as ground-based surveys. Modelling yields a unique physical solution: a mid-M+M-dwarf binary with $M_{\rm prim} = 0.20 \pm 0.01 M_\odot$ and $M_{\rm sec} = 0.16 \pm 0.01 M_\odot$, with projected separation of $2.0 \pm 0.3$ AU. The lens is only $D_{\rm LS} = 0.41 \pm 0.06$ kpc in front of the source, making OB140962 a bulge lens and the most distant Spitzer binary lens to date. In contrast, because the Einstein radius ($\theta_{\rm E} = 0.143 \pm 0.007$ mas) is unusually small, a standard Bayesian analysis, conducted in the absence of parallax information, would predict a brown dwarf binary. We test the accuracy of Bayesian analysis over a set of Spitzer lenses, finding overall good agreement throughout the sample. We also illustrate the methodology for probing the Galactic distribution of planets by comparing the cumulative distance distribution of the Spitzer 2-body lenses to that of the Spitzer single lenses.

The gas is the dominant component of baryonic matter in most galaxy groups and clusters. The spatial offsets of gas centre from the halo centre could be an indicator of the dynamical state of cluster. Knowledge of such offsets is important for estimate the uncertainties when using clusters as cosmological probes. In this paper, we study the centre offsets $r_{\rm off}$ between the gas and that of all the matter within halo systems in $\Lambda$CDM cosmological hydrodynamic simulations. We focus on two kinds of centre offsets: one is the three-dimensional PB offsets between the gravitational potential minimum of the entire halo and the barycentre of the ICM, and the other is the two-dimensional PX offsets between the potential minimum of the halo and the iterative centroid of the projected synthetic X-ray emission of the halo. Halos at higher redshifts tend to have larger values of rescaled offsets $r_{\rm off}/r_{200}$ and larger gas velocity dispersion $\sigma_{v}^{\rm gas}/\sigma_{200}$. For both types of offsets, we find that the correlation between the rescaled centre offsets $r_{\rm off}/r_{200}$ and the rescaled 3D gas velocity dispersion, $\sigma_v^{\rm gas}/\sigma_{200}$ can be approximately described by a quadratic function as $r_{{\rm off}}/r_{200} \propto (\sigma_{v}^{\rm gas}/\sigma_{200} - k_2)^{2}$. A Bayesian analysis with MCMC method is employed to estimate the model parameters. Dependence of the correlation relation on redshifts and the gas mass fraction are also investigated.

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

Schwinn et al. (2017) have recently compared the abundance and distribution of massive substructures identified in a gravitational lensing analysis of Abell 2744 by Jauzac et al. (2016) and N-body simulation and found no cluster in {\Lambda}CDM simulation that is similar to Abell 2744. Schwinn et al.(2017) identified the measured projected aperture masses with the actual masses associated with subhaloes in the MXXL N-body simulation. We have used the high resolution Phoenix cluster simulations to show that such an identification is incorrect: the aperture mass is dominated by mass in the body of the cluster that happens to be projected along the line-of-sight to the subhalo. This enhancement varies from factors of a few to factors of more than 100, particularly for subhaloes projected near the centre of the cluster. We calculate aperture masses for subhaloes in our simulation and compare them to the measurements for Abell 2744. We find that the data for Abell 2744 are in excellent agreement with the matched predictions from {\Lambda}CDM. We provide further predictions for aperture mass functions of subhaloes in idealized surveys with varying mass detection thresholds.

Asteroids in mean motion resonances with giant planets are common in the solar system, but it was not until recently that several asteroids in retrograde mean motion resonances with Jupiter and Saturn were discovered. A retrograde co-orbital asteroid of Jupiter, 2015 BZ509 is confirmed to be in a long-term stable retrograde 1:1 mean motion resonance with Jupiter, which gives rise to our interests in its unique resonant dynamics. In this paper, we investigate the phase-space structure of the retrograde 1:1 resonance in detail within the framework of the circular restricted three-body problem. We construct a simple integrable approximation for the planar retrograde resonance using canonical contact transformation and numerically employ the averaging procedure in closed form. The phase portrait of the retrograde 1:1 resonance is depicted with the level curves of the averaged Hamiltonian. We thoroughly analyze all possible librations in the co-orbital region and uncover a new apocentric libration for the retrograde 1:1 resonance inside the planet's orbit. We also observe the significant jumps in orbital elements for outer and inner apocentric librations, which are caused by close encounters with the perturber.

This paper studies the binary disruption problem and asteroid capture mechanism in a sun-planet-binary four-body system. Firstly, the binary disruption condition is studied and the result shows that the binary is always disrupted at the perigee of their orbit instantaneously. Secondly, an analytic expression to describe the energy exchange between the binary is derived based on the instantaneous disruption hypothesis. The analytic result is validated through numerical integration. We obtain the energy exchange in encounters simultaneously by the analytic expression and numerical integration. The maximum deviation of this two results is always less than 25% and the mean deviation is about 8.69%. The analytic expression can give us an intuitive description of the energy exchange between the binary. It indicates that the energy change depends on the hyperbolic shape of the binary orbit with respect to the planet, the masses of planet and the primary member of the binary, the binary phase at perigee. We can illustrate the capture/escape processes and give the capture/escape region of the binary clearly by numerical simulation. We analyse the influence of some critical factors to the capture region finally.

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

We present the analysis of microlensing event MOA-2010-BLG-117, and show that the light curve can only be explained by the gravitational lensing of a binary source star system by a star with a Jupiter mass ratio planet. It was necessary to modify standard microlensing modeling methods to find the correct light curve solution for this binary-source, binary-lens event. We are able to measure a strong microlensing parallax signal, which yields the masses of the host star, $M_* = 0.58\pm 0.11 M_\odot$, and planet $m_p = 0.54\pm 0.10 M_{\rm Jup}$ at a projected star-planet separation of $a_\perp = 2.42\pm 0.26\,$AU, corresponding to a semi-major axis of $a = 2.9{+1.6\atop -0.6}\,$AU. Thus, the system resembles a half-scale model of the Sun-Jupiter system with a half-Jupiter mass planet orbiting a half-solar mass star at very roughly half of Jupiter's orbital distance from the Sun. The source stars are slightly evolved, and by requiring them to lie on the same isochrone, we can constrain the source to lie in the near side of the bulge at a distance of $D_S = 6.9 \pm 0.7\,$kpc, which implies a distance to the planetary lens system of $D_L = 3.5\pm 0.4\,$kpc. The ability to model unusual planetary microlensing events, like this one, will be necessary to extract precise statistical information from the planned large exoplanet microlensing surveys, such as the WFIRST microlensing survey.

We report on the analysis of a microlensing event OGLE-2014-BLG-1722 that showed two distinct short term anomalies. The best fit model to the observed light curves shows that the two anomalies are explained with two planetary mass ratio companions to the primary lens. Although a binary source model is also able to explain the second anomaly, it is marginally ruled out by 3.1 $\sigma$. The 2-planet model indicates that the first anomaly was caused by planet "b" with a mass ratio of $q = (4.5_{-0.6}^{+0.7}) \times 10^{-4}$ and projected separation in unit of the Einstein radius, $s = 0.753 \pm 0.004$. The second anomaly reveals planet "c" with a mass ratio of $q_{2} = (7.0_{-1.7}^{+2.3}) \times 10^{-4}$ with $\Delta \chi^{2} \sim 170$ compared to the single planet model. Its separation has a so-called close-wide degeneracy. We estimated the physical parameters of the lens system from Bayesian analysis. This gives that the masses of planet b and c are $m_{\rm b} = 56_{-33}^{+51}\,M_{\oplus}$ and $m_{\rm c} = 85_{-51}^{+86}\,M_{\oplus}$, respectively, and they orbit a late type star with a mass of $M_{\rm host} = 0.40_{-0.24}^{+0.36}\,M_{\odot}$ located at $D_{\rm L} = 6.4_{-1.8}^{+1.3}\,\rm kpc$ from us. If the 2-planet model is true, this is the third multiple planet system detected by using the microlensing method, and the first multiple planet system detected in the low magnification events, which are dominant in the microlensing survey data. The occurrence rate of multiple cold gas giant systems is estimated using the two such detections and a simple extrapolation of the survey sensitivity of 6 year MOA microlensing survey (Suzuki et al. 2016) combined with the 4 year $\mu$FUN detection efficiency (Gould et al. 2010). It is estimated that $6 \pm 2\,\%$ of stars host two cold giant planets.

In this paper, we will give a general introduction to the project of Ali CMB Polarization Telescope (AliCPT), which is a Sino-US joint project led by the Institute of High Energy Physics (IHEP) and has involved many different institutes in China. It is the first ground-based Cosmic Microwave Background (CMB) polarization experiment in China and an integral part of China's Gravitational Waves Program. The main scientific goal of AliCPT project is to probe the primordial gravitational waves (PGWs) originated from the very early Universe. The AliCPT project includes two stages. The first stage referred to as AliCPT-1, is to build a telescope in the Ali region of Tibet with an altitude of 5,250 meters. Once completed, it will be the worldwide highest ground-based CMB observatory and open a new window for probing PGWs in northern hemisphere. AliCPT-1 telescope is designed to have about 7,000 TES detectors at 90GHz and 150GHz. The second stage is to have a more sensitive telescope (AliCPT-2) with the number of detectors more than 20,000. Our simulations show that AliCPT will improve the current constraint on the tensor-to-scalar ratio $r$ by one order of magnitude with 3 years' observation. Besides the PGWs, the AliCPT will also enable a precise measurement on the CMB rotation angle and provide a precise test on the CPT symmetry. We show 3 years' observation will improve the current limit by two order of magnitude.

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $\tau \propto \lambda^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models.

In recent years there has been a lot of interest in discussing frame dependences/independences of the cosmological perturbations under the conformal transformations. This problem has previously been investigated in terms of the covariant approach for a single component universe, and it was found that the covariant approach is very powerful to pick out the perturbative variables which are both gauge and conformal invariant. In this work, we extend the covariant approach to a universe with multicomponent fluids. We find that similar results can be derived, as expected. In addition, some other interesting perturbations are also identified to be conformal invariant, such as entropy perturbation between two different components.

We report the discovery of a sub-Jupiter mass planet orbiting beyond the snow line of an M-dwarf most likely in the Galactic disk as part of the joint Spitzer and ground-based monitoring of microlensing planetary anomalies toward the Galactic bulge. The microlensing parameters are strongly constrained by the light curve modeling and in particular by the Spitzer-based measurement of the microlens parallax, $\pi_\mathrm{E}$. However, in contrast to many planetary microlensing events, there are no caustic crossings, so the angular Einstein radius, $\theta_\mathrm{E}$ has only an upper limit based on the light curve modeling alone. Additionally, the analysis leads us to identify 8 degenerate configurations: the four-fold microlensing parallax degeneracy being doubled by a degeneracy in the caustic structure present at the level of the ground-based solutions. To pinpoint the physical parameters, and at the same time to break the parallax degeneracy, we make use of a series of arguments: the $\chi^2$ hierarchy, the Rich argument, and a prior Galactic model. The preferred configuration is for a host at $D_L=3.73_{-0.67}^{+0.66}~\mathrm{kpc}$ with mass $M_\mathrm{L}=0.30_{-0.12}^{+0.15}~\mathrm{M_\odot}$, orbited by a Saturn-like planet with $M_\mathrm{planet}=0.43_{-0.17}^{+0.21}~\mathrm{M_\mathrm{Jup}}$ at projected separation $a_\perp = 1.70_{-0.39}^{+0.38}~\mathrm{au}$, about 2.1 times beyond the system snow line. Therefore, it adds to the growing population of sub-Jupiter planets orbiting near or beyond the snow line of M-dwarfs discovered by microlensing. Based on the rules of the real-time protocol for the selection of events to be followed up with Spitzer, this planet will not enter the sample for measuring the Galactic distribution of planets.

We present the analysis of the binary-microlensing event OGLE-2014-BLG-0289. The event light curve exhibits very unusual five peaks where four peaks were produced by caustic crossings and the other peak was produced by a cusp approach. It is found that the quintuple-peak features of the light curve provide tight constraints on the source trajectory, enabling us to precisely and accurately measure the microlensing parallax $\pi_{\rm E}$. Furthermore, the three resolved caustics allow us to measure the angular Einstein radius $\thetae$. From the combination of $\pi_{\rm E}$ and $\thetae$, the physical lens parameters are uniquely determined. It is found that the lens is a binary composed of two M dwarfs with masses $M_1 = 0.52 \pm 0.04\ M_\odot$ and $M_2=0.42 \pm 0.03\ M_\odot$ separated in projection by $a_\perp = 6.4 \pm 0.5$ au. The lens is located in the disk with a distance of $D_{\rm L} = 3.3 \pm 0.3$~kpc. It turns out that the reason for the absence of a lensing signal in the {\it Spitzer} data is that the time of observation corresponds to the flat region of the light curve.

We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$\beta$ lags, we measure black hole masses for all five targets. We also obtain H$\gamma$ and He{\sc ii}\,$\lambda 4686$ lags for all objects except 3C 382. The He{\sc ii}\,$\lambda 4686$ lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100--300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

We report the discovery of OGLE-2016-BLG-1190Lb, which is likely to be the first Spitzer microlensing planet in the Galactic bulge/bar, an assignation that can be confirmed by two epochs of high-resolution imaging of the combined source-lens baseline object. The planet's mass M_p= 13.4+-0.9 M_J places it right at the deuterium burning limit, i.e., the conventional boundary between "planets" and "brown dwarfs". Its existence raises the question of whether such objects are really "planets" (formed within the disks of their hosts) or "failed stars" (low mass objects formed by gas fragmentation). This question may ultimately be addressed by comparing disk and bulge/bar planets, which is a goal of the Spitzer microlens program. The host is a G dwarf M_host = 0.89+-0.07 M_sun and the planet has a semi-major axis a~2.0 AU. We use Kepler K2 Campaign 9 microlensing data to break the lens-mass degeneracy that generically impacts parallax solutions from Earth-Spitzer observations alone, which is the first successful application of this approach. The microlensing data, derived primarily from near-continuous, ultra-dense survey observations from OGLE, MOA, and three KMTNet telescopes, contain more orbital information than for any previous microlensing planet, but not quite enough to accurately specify the full orbit. However, these data do permit the first rigorous test of microlensing orbital-motion measurements, which are typically derived from data taken over <1% of an orbital period.

As one of the most popular applications of lidar systems, the atmospheric visibility is defined to be inversely proportional to the atmospheric extinction coefficient at 0.55 um. Since the laser at 1.5 um shows the highest maximum permissible exposure in the wavelength ranging from 0.3 um to 10 um, the eye-safe 1.5 um lidar can be deployed in urban areas. In such a case, the measured extinction coefficient at 1.5 um should be converted to that at 0.55 um for visibility retrieval. Although several models have been established since 1962, the accurate wavelength conversion remains a challenge. An adaptive inversion algorithm for 1.5 um visibility lidar is proposed and demonstrated by using the in situ Angstrom wavelength exponent, which is derived from either a sun photometer or an aerosol spectrometer. The impact of the particle size distribution of atmospheric aerosols and the Rayleigh backscattering of atmospheric molecules are taken into account. In comparison, the Angstrom wavelength exponent derived from the sun photometer is 7.7% higher than that derived from the aerosol spectrometer. Then, using 1.5 um visibility lidar, the visibility with a temporal resolution of 5 min is detected over 48 hours in Hefei (31.83 N, 117.25 E). The average visibility error between the new method and a visibility sensor (Vaisala, PWD52) is 5.2% with the R-square value of 0.96, while the relative error between another reference visibility lidar at 532 nm and the visibility sensor is 6.7% with the R-square value of 0.91. All results agree with each other well, demonstrating the accuracy and stability of the algorithm.

We present the discovery of a substellar companion to the primary host lens in the microlensing event MOA-2012-BLG-006. The companion-to-host mass ratio is 0.016, corresponding to a companion mass of $\approx8~M_{\rm Jup} (M_*/0.5M_\odot)$. Thus, the companion is either a high-mass giant planet or a low-mass brown dwarf, depending on the mass of the primary $M_*$. The companion signal was separated from the peak of the primary event by a time that was as much as four times longer than the event timescale. We therefore infer a relatively large projected separation of the companion from its host of $\approx10~{\rm a.u.}(M_*/0.5M_\odot)^{1/2}$ for a wide range (3-7 kpc) of host star distances from the Earth. We also challenge a previous claim of a planetary companion to the lens star in microlensing event OGLE-2002-BLG-045.

We present the discovery of two planetary systems consisting of a Saturn-mass planet orbiting an M-dwarf, which were detected in faint microlensing events OGLE-2013-BLG-0132 and OGLE-2013-BLG-1721. The planetary anomalies were covered with high cadence by OGLE and MOA photometric surveys. The light curve modeling indicates that the planet-to-host mass ratios are $(5.15 \pm 0.28)\times 10^{-4}$ and $(13.18 \pm 0.72)\times 10^{-4}$, respectively. Both events were too short and too faint to measure a reliable parallax signal and hence the lens mass. We therefore used a Bayesian analysis to estimate the masses of both planets: $0.29^{+0.16}_{-0.13}\ M_{Jup}$ (OGLE-2013-BLG-0132Lb) and $0.64^{+0.35}_{-0.31}\ M_{Jup}$ (OGLE-2013-BLG-1721Lb). Thanks to a high relative proper motion, OGLE-2013-BLG-0132 is a promising candidate for the high-resolution imaging follow-up. Both planets belong to an increasing sample of sub-Jupiter-mass planets orbiting M-dwarfs beyond the snow line.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

We present the result of microlensing event MOA-2016-BLG-290, which received observations from the two-wheel Kepler (K2), Spitzer, as well as ground-based observatories. A joint analysis of data from K2 and the ground leads to two degenerate solutions of the lens mass and distance. This degeneracy is effectively broken once the (partial) Spitzer light curve is included. Altogether, the lens is found to be an extremely low-mass star located in the Galactic bulge. MOA-2016-BLG-290 is the first microlensing event for which we have signals from three well-separated ($\sim1$ AU) locations. It demonstrates the power of two-satellite microlensing experiment in reducing the ambiguity of lens properties, as pointed out independently by S. Refsdal and A. Gould several decades ago.

Usually when applying the mimetic model to the early universe, higher derivative terms are needed to promote the mimetic field to be dynamical. However such models suffer from the ghost and/or the gradient instabilities and simple extensions cannot cure this pathology. We point out in this paper that it is possible to overcome this difficulty by considering the direct couplings of the higher derivatives of the mimetic field to the curvature of the spacetime.

We report the discovery of a planet-mass companion to the microlens OGLE-2016-BLG-0263L. Unlike most low-mass companions that were detected through perturbations to the smooth and symmetric light curves produced by the primary, the companion was discovered through the channel of a repeating event, in which the companion itself produced its own single-mass light curve after the event produced by the primary had ended. Thanks to the continuous coverage of the second peak by high-cadence surveys, the possibility of the repeating nature due to source binarity is excluded with a $96\%$ confidence level. The mass of the companion estimated by a Bayesian analysis is $M_{\rm p}=4.1_{-2.5}^{+6.5}\ M_{\rm J}$. The projected primary-companion separation is $a_\perp = 6.5^{+1.3}_{-1.9}$ au. The ratio of the separation to the snow-line distance of $a_\perp/a_{\rm sl}\sim 15.4$ corresponds to the region beyond Neptune, the outermost planet of the solar system. We discuss the importance of high-cadence surveys in expanding the range of microlensing detections of low-mass companions and future space-based microlensing surveys.

We obtained new optical observations of the X-ray source XMMU J083850.38-282756.8, the previously proposed counterpart of the gamma-ray source 3FGL J0838.8-2829. Time-series photometry in the $r'$ band reveals periodic modulation of $\approx1$ magnitude that is characteristic of the heating of the photosphere of a low-mass companion star by a compact object. The measured orbital period is 5.14817$\pm$0.00012 hr. The shape of the light curve is variable, evidently due to the effects of flaring and asymmetric heating. Spectroscopy reveals a companion of type M1 or later, having a radial velocity amplitude of 315$\pm$17 km/s, with period and phasing consistent with the heating interpretation. The mass function of the compact object is 0.69$\pm$0.11 $M_{\odot}$, which allows a neutron star in a high-inclination orbit. Variable, broad H$\alpha$ emission is seen, which is probably associated with a wind from the companion. These properties, as well as the X-ray and gamma-ray luminosities at the inferred distance of <1.7 kpc, are consistent with a redback millisecond pulsar in its non-accreting state. A search for radio pulsations is needed to confirm this interpretation and derive complete system parameters for modeling, although absorption by the ionized wind could hinder such detection.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

We present the analysis of the binary gravitational microlensing event MOA-2015-BLG-020. The event has a fairly long timescale (about 63 days) and thus the light curve deviates significantly from the lensing model that is based on the rectilinear lens-source relative motion. This enables us to measure the microlensing parallax through the annual parallax effect. The microlensing parallax parameters constrained by the ground-based data are confirmed by the Spitzer observations through the satellite parallax method. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined. It is found that the binary lens is composed of two dwarf stars with masses $M_1 = 0.606 \pm 0.028M_\odot$ and $M_2 = 0.125 \pm 0.006M_\odot$ in the Galactic disk. Assuming the source star is at the same distance as the bulge red clump stars, we find the lens is at a distance $D_L = 2.44 \pm 0.10 kpc$. In the end, we provide a summary and short discussion of all published microlensing events in which the annual parallax effect is confirmed by other independent observations.

We review the paradigm of holographic dark energy (HDE), which arises from a theoretical attempt of applying the holographic principle (HP) to the dark energy (DE) problem. Making use of the HP and the dimensional analysis, we derive the general formula of the energy density of HDE. Then, we describe the properties of HDE model, in which the future event horizon is chosen as the characteristic length scale. We also introduce the theoretical explorations and the observational constraints for this model. Next, in the framework of HDE, we discuss various topics, such as spatial curvature, neutrino, instability of perturbation, time-varying gravitational constant, inflation, black hole and big rip singularity. In addition, from both the theoretical and the observational aspects, we introduce the interacting holographic dark energy scenario, where the interaction between dark matter and HDE is taken into account. Furthermore, we discuss the HDE scenario in various modified gravity (MG) theories, such as Brans-Dicke theory, braneworld theory, scalar-tensor theory, Horava-Lifshitz theory, and so on. Besides, we introduce the attempts of reconstructing various scalar-field DE and MG models from HDE. Moreover, we introduce other DE models inspired by the HP, in which different characteristic length scales are chosen. Finally, we make comparisons among various HP-inspired DE models, by using cosmological observations and diagnostic tools.

We present the first catalogue of eclipsing binaries in two MOA fields towards the Galactic bulge, in which over 8,000 candidates, mostly contact and semi-detached binaries of periods < 1 d, were identified. In this paper, the light curves of a small number of interesting candidates including eccentric binaries, binaries with noteworthy phase modulations and eclipsing RS CVn type stars are shown as examples. In addition, we identified three triple object candidates by detecting the light-travel-time effect in their eclipse time variation curves.

We report the discovery of a super-Earth mass planet in the microlensing event MOA-2012-BLG-505. This event has the second shortest event timescale of $t_{\rm E}=10 \pm 1$ days where the observed data show evidence of planetary companion. Our 15 minute high cadence survey observation schedule revealed the short subtle planetary signature. The system shows the well known close/wide degeneracy. The planet/host-star mass ratio is $q =2.1 \times 10^{-4}$ and the projected separation normalized by the Einstein radius is s = 1.1 or 0.9 for the wide and close solutions, respectively. We estimate the physical parameters of the system by using a Bayesian analysis and find that the lens consists of a super-Earth with a mass of $6.7^{+10.7}_{-3.6}M_{\oplus}$ orbiting around a brown-dwarf or late M-dwarf host with a mass of $0.10^{+0.16}_{-0.05}M_{\odot}$ with a projected star-planet separation of $0.9^{+0.3}_{-0.2}$AU. The system is at a distance of $7.2 \pm 1.1$ kpc, i.e., it is likely to be in the Galactic bulge. The small angular Einstein radius ($\theta_{\rm E}=0.12 \pm 0.02$ mas) and short event timescale are typical for a low-mass lens in the Galactic bulge. Such low-mass planetary systems in the Bulge are rare because the detection efficiency of planets in short microlensing events is relatively low. This discovery may suggest that such low mass planetary systems are abundant in the Bulge and currently on-going high cadence survey programs will detect more such events and may reveal an abundance of such planetary systems.