The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. We construct heavy neutron stars by introducing non-trivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

M31N 2017-09a is a classical nova and was observed for some 160 days following its initial eruption, during which time it underwent a number of bright secondary outbursts. The light-curve is characterized by continual variation with excursions of at least 0.5 magnitudes on a daily time-scale. The lower envelope of the eruption suggests that a single power-law can describe the decline rate. The eruption is relatively long with $t_2 = 111$, and $t_3 = 153$ days.

We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

In this paper we report the results of the first $\sim$four years of spectroscopic and photometric monitoring of the Type IIn supernova SN 2015da (also known as PSN J13522411+3941286, or iPTF16tu). The supernova exploded in the nearby spiral galaxy NGC 5337 in a relatively highly extinguished environment. The transient showed prominent narrow Balmer lines in emission at all times and a slow rise to maximum in all bands. In addition, early observations performed by amateur astronomers give a very well-constrained explosion epoch. The observables are consistent with continuous interaction between the supernova ejecta and a dense and extended H-rich circumstellar medium. The presence of such an extended and dense medium is difficult to reconcile with standard stellar evolution models, since the metallicity at the position of SN 2015da seems to be slightly subsolar. Interaction is likely the mechanism powering the light curve, as confirmed by the analysis of the pseudo bolometric light curve, which gives a total radiated energy $\gtrsim10^{51}\,\rm{erg}$. Modeling the light curve in the context of a supernova shock breakout through a dense circumstellar medium allowed us to infer the mass of the prexisting gas to be $\simeq8\,\rm{M_{\odot}}$, with an extreme mass-loss rate for the progenitor star $\simeq0.6\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, suggesting that most of the circumstellar gas was produced during multiple eruptive events. Near- and mid-infrared observations reveal a flux excess in these domains, similar to those observed in SN 2010jl and other interacting transients, likely due to preexisting radiatively heated dust surrounding the supernova. By modeling the infrared excess, we infer a mass $\gtrsim0.4\times10^{-3}\,\rm{M_{\odot}}$ for the dust.

We study the tidal deformations of various known black hole and wormhole solutions in a simple context of warped compactification --- Randall-Sundrum theory in which the four-dimensional spacetime geometry is that of a brane embedded in five-dimensional Anti-de Sitter spacetime. The linearized gravitational perturbation theory generically reduces to either an inhomogeneous second-order ODE or a homogeneous third-order ODE of which indicial roots associated with an expansion about asymptotic infinity can be related to Tidal Love Numbers. We describe various tidal-deformed metrics, classify their indicial roots, and find that in particular the quadrupolar TLN is generically non-vanishing. Thus it could be a signature of a braneworld by virtue of its potential appearance in gravitational waveforms emitted in binary merger events.

We study one-loop corrections to the two-point correlation function of tensor perturbations in primordial cosmology induced by massless spectator matter fields. Using the Schwinger-Keldysh formalism in cosmological perturbation theory, we employ dimensional regularization and cutoff regularization to compute the finite quantum corrections at one-loop arising from massless isocurvature fields of various spins <2 which are freely propagating on the FRW spacetime. For all cases, we find a logarithmic running of the form $\frac{C}{q^3} \frac{H^4}{M^4_p} \log \left( \frac{H}{\mu } \right)$ where $C$ is a negative constant, $H$ is the Hubble parameter at horizon exit and $\mu$ is the renormalization scale.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

We report on a multi-observatory campaign to examine asteroid 3200 Phaethon during its December 2017 close approach to Earth, in order to improve our measurements of its fundamental parameters, and to search for surface variations, cometary activity and fragmentation. The mean colors of Phaethon are B-V = 0.702 +/- 0.004, V-R = 0.309 +/- 0.003, R-I = 0.266 +/- 0.004, neutral to slightly blue, consistent with previous classifications of Phaethon as a F-type or B-type asteroid. Variations in Phaethon's B-V colors (but not V-R or R-I) with observer sub-latitude are seen and may be associated with craters observed by the Arecibo radar. High cadence photometry over phases from 20 to 100 degrees allows a fit to the values of the HG photometric parameters; H = 14.57 +/- 0.02, 13.63 +/- 0.02, 13.28 +/- 0.02, 13.07 +/- 0.02; G = 0.00 +/- 0.01, -0.09 +/- 0.01, -0.10 +/- 0.01, -0.08 +/- 0.01 in the BVRI filters respectively; the negative G values are consistent with other observations of F type asteroids. Light curve variations were seen that are also consistent with concavities reported by Arecibo, indicative of large craters on Phaethon's surface whose ejecta may be the source of the Geminid meteoroid stream. A search for gas/dust production set an upper limit of 0.06 +/- 0.02 kg/s when Phaethon was 1.449 AU from the Sun, and 0.2 +/- 0.1 kg/s at 1.067 AU. A search for meter-class fragments accompanying Phaethon did not find any whose on-sky motion was not also consistent with background main belt asteroids.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

We describe ZStreak, a semi-real-time pipeline specialized in detecting small, fast-moving near-Earth asteroids (NEAs) that is currently operating on the data from the newly-commissioned Zwicky Transient Facility (ZTF) survey. Based on a prototype originally developed by Waszczak et al. (2017) for the Palomar Transient Factory (PTF), the predecessor of ZTF, ZStreak features an improved machine-learning model that can cope with the $10\times$ data rate increment between PTF and ZTF. Since its first discovery on 2018 February 5 (2018 CL), ZTF/ZStreak has discovered $45$ confirmed new NEAs over a total of 232 observable nights until 2018 December 31. Most of the discoveries are small NEAs, with diameters less than $\sim100$ m. By analyzing the discovery circumstances, we find that objects having the first to last detection time interval under 2 hr are at risk of being lost. We will further improve real-time follow-up capabilities, and work on suppressing false positives using deep learning.

The spectral flux density of stars can indicate their atmospheric physical properties. A detector can obtain any band flux density at the design stage. However, the band flux density is confirmed and fixed in the process of operation because of the restriction of filters. Other band flux densities cannot be obtained through the same detector. In this study, a computational model of stellar spectral flux density is established based on basic physical parameters which are effective temperature and angular parameter. The stochastic particle swarm optimization algorithm is adopted to address this issue with appropriately chosen values of the algorithm parameters. Four star catalogues are studied and consist of the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST), Wide-field Infrared Survey Explorer (WISE), Midcourse Space Experiment (MSX), and Two Micron All Sky Survey (2MASS). The given flux densities from catalogues are input parameters. Stellar effective temperatures and angular parameters are inversed using the given flux densities according to SPSO algorithm. Then the flux density is calculated according to Planck's law on the basis of stellar effective temperatures and angular parameters. The calculated flux density is compared with the given value from catalogues. It is found that the inversion results are in good agreement for all bands of the MSX and 2MASS catalogues, whereas they do not agree well in some bands of the LAMOST and WISE catalogues. Based on the results, data from the MSX and 2MASS catalogues can be used to calculate the spectral flux density at different wavelengths of given wavelength ranges. The stellar flux density is obtained and can provide data support and an effective reference for detection and recognition of stars.

Detection schemes for the quantum chromodynamics axions and other axion-like particles in light-shining-through-a-wall (LSW) experiments are based on the conversion of these particles into photons in a magnetic field. An alternative scheme may involve the detection via a resonant atomic or molecular transition induced by resonant axion absorption. The signal obtained in this process is second order in the axion-electron interaction constant but may become first order if we allow interference between the axion-induced transition amplitude and the transition amplitude induced by the electromagnetic radiation that produces the axions.

Since its discovery in 2008, the Andromeda galaxy nova M31N 2008-12a has been observed in eruption every single year. This unprecedented frequency indicates an extreme object, with a massive white dwarf and a high accretion rate, which is the most promising candidate for the single-degenerate progenitor of a type-Ia supernova known to date. The previous three eruptions of M31N 2008-12a have displayed remarkably homogeneous multi-wavelength properties: (i) From a faint peak, the optical light curve declined rapidly by two magnitudes in less than two days; (ii) Early spectra showed initial high velocities that slowed down significantly within days and displayed clear He/N lines throughout; (iii) The supersoft X-ray source (SSS) phase of the nova began extremely early, six days after eruption, and only lasted for about two weeks. In contrast, the peculiar 2016 eruption was clearly different. Here we report (i) the considerable delay in the 2016 eruption date, (ii) the significantly shorter SSS phase, and (iii) the brighter optical peak magnitude (with a hitherto unobserved cusp shape). Early theoretical models suggest that these three different effects can be consistently understood as caused by a lower quiescence mass-accretion rate. The corresponding higher ignition mass caused a brighter peak in the free-free emission model. The less-massive accretion disk experienced greater disruption, consequently delaying re-establishment of effective accretion. Without the early refueling, the SSS phase was shortened. Observing the next few eruptions will determine whether the properties of the 2016 outburst make it a genuine outlier in the evolution of M31N 2008-12a.

Considering features of stellar spectral radiation and survey explorers, we established a computational model for stellar effective temperatures, detected angular parameters, and gray rates. Using known stellar flux data in some band, we estimated stellar effective temperatures and detected angular parameters using stochastic particle swarm optimization (SPSO). We first verified the reliability of SPSO, and then determined reasonable parameters that produced highly accurate estimates under certain gray deviation levels. Finally, we calculated 177,860 stellar effective temperatures and detected angular parameters using the Midcourse Space Experiment (MSX) catalog data. These derived stellar effective temperatures were accurate when we compared them to known values from literatures. This research made full use of catalog data and presented an original technique for studying stellar characteristics. It proposed a novel method for calculating stellar effective temperatures and detected angular parameters, and provided theoretical and practical data support for finding radiation flow information for any band.