Any high-contrast imaging instrument in a future large space-based telescope will include an integral field spectrograph (IFS) for measuring broadband starlight residuals and characterizing the exoplanet's atmospheric spectrum. In this paper, we report the development of a high-contrast integral field spectrograph (HCIFS) at Princeton University and demonstrate its application in multi-spectral wavefront control. Moreover, we propose and experimentally validate a new reduced-dimensional system identification algorithm for an IFS imaging system, which improves the system's wavefront control speed, contrast and computational and data storage efficiency.

Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

The multi-messenger observation of compact binary coalescence promises great scientific treasure. However, a synthetic observation from both gravitational wave and electromagnetic channels remains challenging. Relying on the day-to-week long macronova emission, GW170817 remains the only event with successful electromagnetic followup. In this manuscript, we explore the possibility of using the early stage X-ray afterglow to search for the electromagnetic counterpart of gravitational wave events. Two algorithms, the sequential observation and the local optimization are considered and applied to three simulated events. We consider the proposed Einstein probe as a candidate X-ray telescope. Benefiting from the large field of view and high sensitivity, we find that the sequential observation algorithm not only is easy to implement, but also promises a good chance of actual detection.

Exoplanets are very dim. An Electron Multiplying Charged Coupled Device (EMCCD) operating in photon counting (PC) mode reduces the noise level and enables their detection. Typically, PC images are added together as a co-added image before processing. We present here a signal detection and estimation technique working directly with the PC images. The method is based on the generalized maximum likelihood ratio test (GLRT) and uses a binomial distribution between PC images. It can be applied online, so it is possible to stop observations once we have reached a threshold confidence for the existence (or absence) of planets, and thus the observation time is efficiently used. It directly uses the detector model, which more accurately represents its noise characteristics. We show that it outperforms our previous GLRT method using co-added images under Gaussian noise assumption.

Optimized sensing is important for computational imaging in low-resource environments, when images must be recovered from severely limited measurements. In this paper, we propose a physics-constrained, fully differentiable, autoencoder that learns a probabilistic sensor-sampling strategy for optimized sensor design. The proposed method learns a system's preferred sampling distribution that characterizes the correlations between different sensor selections as a binary, fully-connected Ising model. The learned probabilistic model is achieved by using a Gibbs sampling inspired network architecture, and is trained end-to-end with a reconstruction network for efficient co-design. The proposed framework is applicable to sensor selection problems in a variety of computational imaging applications. In this paper, we demonstrate the approach in the context of a very-long-baseline-interferometry (VLBI) array design task, where sensor correlations and atmospheric noise present unique challenges. We demonstrate results broadly consistent with expectation, and draw attention to particular structures preferred in the telescope array geometry that can be leveraged to plan future observations and design array expansions.

We conduct a post hoc analysis of solar flare predictions made by a Long Short Term Memory (LSTM) model employing data in the form of Space-weather HMI Active Region Patches (SHARP) parameters calculated from data in proximity to the magnetic polarity inversion line where the flares originate. We train the the LSTM model for binary classification to provide a prediction score for the probability of M/X class flares to occur in next hour. We then develop a dimension-reduction technique to reduce the dimensions of SHARP parameter (LSTM inputs) and demonstrate the different patterns of SHARP parameters corresponding to the transition from low to high prediction score. Our work shows that a subset of SHARP parameters contain the key signals that strong solar flare eruptions are imminent. The dynamics of these parameters have a highly uniform trajectory for many events whose LSTM prediction scores for M/X class flares transition from very low to very high. The results demonstrate the existence of a few threshold values of SHARP parameters that when surpassed indicate a high probability of the eruption of a strong flare. Our method has distilled the knowledge of solar flare eruption learnt by deep learning model and provides a more interpretable approximation, which provides physical insight to processes driving solar flares.

We develop a mixed Long Short Term Memory (LSTM) regression model to predict the maximum solar flare intensity within a 24-hour time window 0$\sim$24, 6$\sim$30, 12$\sim$36 and 24$\sim$48 hours ahead of time using 6, 12, 24 and 48 hours of data (predictors) for each Helioseismic and Magnetic Imager (HMI) Active Region Patch (HARP). The model makes use of (1) the Space-weather HMI Active Region Patch (SHARP) parameters as predictors and (2) the exact flare intensities instead of class labels recorded in the Geostationary Operational Environmental Satellites (GOES) data set, which serves as the source of the response variables. Compared to solar flare classification, the model offers us more detailed information about the exact maximum flux level, i.e. intensity, for each occurrence of a flare. We also consider classification models built on top of the regression model and obtain better results in solar flare classifications. Our results suggest that the most efficient time period for predicting the solar activity is within 24 hours before the prediction time using the SHARP parameters and the LSTM model.

A deep learning network, Long-Short Term Memory (LSTM) network, is used in this work to predict whether the maximum flare class an active region (AR) will produce in the next 24 hours is class $\Gamma$. We considered $\Gamma$ are $\ge M$, $\ge C$ and any flare class. The essence of using LSTM, which is a recurrent neural network, is its capability to capture temporal information of the data samples. The input features are time sequences of 20 magnetic parameters from SHARPs - Space-weather HMI Active Region Patches. We analyzed active regions from June 2010 to Dec 2018, using the Geostationary Operational Environmental Satellite (GOES) X-ray flare catalogs and label the data samples with identified ARs in the GOES X-ray flare catalogs. Our results (i) shows consistent skill scores with recently published results using LSTMs and better than the previous work using single time input (eg. DeFN) (ii) The skill scores from the model show essential differences when different years of data was chosen for training and testing.

Future large space telescopes will be equipped with adaptive optics (AO) to overcome wavefront aberrations and achieve high contrast for imaging faint astronomical objects, such as earth-like exoplanets and debris disks. In contrast to AO that is widely used in ground telescopes, space-based AO systems will use focal plane wavefront sensing to measure the wavefront aberrations. Focal plane wavefront sensing is a class of techniques that reconstruct the light field based on multiple focal plane images distorted by deformable mirror (DM) probing perturbations. In this paper, we report an efficient focal plane wavefront sensing approach for space-based AO which optimizes the DM probing perturbation and thus also the integration time for each image. Simulation of the AO system equipped with a vortex coronagraph has demonstrated that our new approach enables efficient information acquisition and significantly reduces the time needed for achieving high contrast in space.

Two bright X-ray transients were reported from the Chandra Deep Field South archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is not identified, the latter was suggested as an excellent candidate for a rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we propose a unified model to interpret both transients within the framework of the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed from the "free zone" where the magnetar spindown powered X-ray emission escapes freely, whereas CDF-S XT1 originates from the "trapped zone" where the X-ray emission is initially blocked by the dynamical ejecta and becomes transparent after the ejecta is pushed to a distance where Thomson optical depth drops below unity. We fit the magnetar model to the light curves of both transients and derived consistent parameters for the two events, with magnetic field, initial spin period and X-ray emission efficiency being ($B_p=10^{16}\,G$, $P=1.2\,\rm ms$, $\eta = 0.001$) and ($B_p=10^{15.8}\,G$, $P=4.4\, \rm ms$, $\eta = 0.001$) for XT1 and XT2, respectively. The "isotropic equivalent" ejecta mass of XT1 is $M_{\rm ej} \sim 10^{-3}$ $M_{\odot}$, while it is not constrained for XT2. Our results suggest that more extreme magnetar parameters are required to have XT1 detected from the trapped zone. The model parameters for both events are generally consistent with those derived from SGRB X-ray plateau observations. The host galaxy properties of both transients are also consistent with those of SGRBs. The event rate densities of both XT1 and XT2 are consistent with that of BNS mergers.

Simulation of plasmas in the electromagnetic fields requires to solve numerically a kinetic equation, describing the time evolution of the particle distribution function. Here, we propose a finite volume scheme based on the integral relation for the Poisson bracket to solve the most fundamental kinetic equation, namely, the Liouville equation. The proposed scheme conserves the number of particles, maintains the total-variation-diminishing (TVD) property, and provides high-quality numerical results. Some other types of kinetic equations may be also formulated in terms of the Poisson brackets and solved with the proposed method. Among them is the focused transport equation describing the acceleration and propagation of the Solar Energetic Particles (SEPs), which is of practical importance, since the high energy SEPs produce radiation hazards. The newly proposed scheme is demonstrated to be accurate and efficient, which makes it applicable to global simulation systems analysing the space weather. We also discuss a role of focused transport and the accuracy of the diffusive approximation, in application to the SEPs

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817.

All coronagraphic instruments for exoplanet high-contrast imaging need wavefront correction systems to reject optical aberrations and create sufficiently dark holes. Since the most efficient wavefront correction algorithms (controllers and estimators) are usually model-based, the modeling accuracy of the system influences the ultimate wavefront correction performance. Currently, wavefront correction systems are typically approximated as linear systems using Fourier optics. However, the Fourier optics model is usually biased due to inaccuracies in the layout measurements, the imperfect diagnoses of inherent optical aberrations, and a lack of knowledge of the deformable mirrors (actuator gains and influence functions). Moreover, the telescope optical system varies over time because of instrument instabilities and environmental effects. In this paper, we present an expectation-maximization (E-M) approach for identifying and real-time adapting the linear telescope model from data. By iterating between the E-step (a Kalman filter and a Rauch smoother) and the M-step (analytical or gradient-based optimization), the algorithm is able to recover the system even if the model depends on the electric fields, which are unmeasurable hidden variables. Simulations and experiments in Princeton's High Contrast Imaging Lab demonstrate that this algorithm improves the model accuracy and increases the efficiency and speed of the wavefront correction.

Double neutron star (DNS) merger events are promosing candidates of short Gamma-ray Burst (sGRB) progenitors as well as high-frequecy gravitational wave (GW) emitters. On August 17, 2017, such a coinciding event was detected by both the LIGO-Virgo gravitational wave detector network as GW170817 and Gamma-Ray Monitor on board NASA's {\it Fermi} Space Telescope as GRB 170817A. Here we show that the fluence and spectral peak energy of this sGRB fall into the lower portion of the distributions of known sGRBs. Its peak isotropic luminosity is abnormally low. The estimated event rate density above this luminosity is at least $190^{+440}_{-160} {\rm Gpc^{-3} \ yr^{-1}}$, which is close to but still below the DNS merger event rate density. This event likely originates from a structured jet viewed from a large viewing angle. There are similar faint soft GRBs in the {\it Fermi} archival data, a small fraction of which might belong to this new population of nearby, low-luminosity sGRBs.

Zhang (2013) proposed a type of GRB-less X-ray transient associated with double neutron star (NS-NS) mergers under the conjecture of a rapidly-spinning magnetar merger product with the line of sight off the short GRB jet. We investigate possible light curves of these transients by considering different observer's viewing angles. We perform Monte Carlo simulations to calculate the peak luminosity function (LF) and event rate density of these X-ray transients. By considering that a fraction of massive neutron stars may be supra-massive and later collapse into black holes after spinning down, we investigate how the predicted LF depends on the equation of state (EoS) of the central object and the geometry of the system. In general, the LF can be fit by two log-normal distributions peaking around $10^{46.4}$ and $10^{49.6}$ $\rm erg\,s^{-1}$, corresponding to the trapped and free zones, respectively. For the majority of the EoS models, the current non-detection is consistent with having a free zone solid angle at most a few times of the solid angle of the short GRB jet. The event rate density of these X-ray transients is around a few tens of $\rm Gpc^{-3}yr^{-1}$ for luminosity above $10^{45}$ $\rm erg\,s^{-1}$. We predict that future X-ray telescopes (such as Einstein Probe) with sensitivity $\sim 10^{-11}$ $\rm erg\,s^{-1}\,cm^{-2}$ would detect as many as several tens of such transients per year per steradian. Within 200 Mpc, the aLIGO average range for NS-NS mergers, the estimated event rate of these transients is about 1 transient per year all sky.

GRB 160821B is a short gamma-ray burst (GRB) at redshift $z=0.16$, with a duration less than 1 second and without detection of any "extended emission" up to more than 100 seconds in both {\em Swift}/BAT and {\em Fermi}/GBM bands. An X-ray plateau with a sharp drop 180 seconds after the BAT trigger was observed with {\em Swift}/XRT. No supernova or kilo-nova signature was detected. Assuming the central engine of this SGRB is a recently born supra-massive magnetar, we can explain the SGRB as jet radiation and its X-ray plateau as the internal energy dissipation of the pulsar wind as it spins down. We constrain its surface magnetic field as $B_{\rm p}<3.12\times 10^{16}$ G and initial spin period as $P_0< 8.5\times 10^{-3}$ seconds. Its equation of state is consistent with the GM1 model with $M_{\rm TOV} \sim 2.37 M_\odot$ and ellipticity $\epsilon<0.07$. Its gravitational wave (GW) radiation may be detectable with the future Einstein Telescope, but is much weaker than the current detectability limit of advanced-LIGO. The GW radiation of such an event would be detectable by advanced-LIGO if it occurred at a distance of 100 Mpc ($z=0.023$).

Several types of extra-galactic high-energy transients have been discovered, which include high-luminosity and low-luminosity long-duration gamma-ray bursts (GRBs), short-duration GRBs, supernova shock breakouts (SBOs), and tidal disruption events (TDEs) without or with an associated relativistic jet. In this paper, we apply a unified method to systematically study the redshift-dependent event rate densities and the global luminosity functions (ignoring redshift evolution) of these transients. We introduce some empirical formulae for the redshift-dependent event rate densities for different types of transients, and derive the local specific event rate density, which also represents its global luminosity function. Long GRBs have a large enough sample to reveal features in the global luminosity function, which is best characterized as a triple power law. All the other transients are consistent with having a single power law luminosity function. The total event rate density depends on the minimum luminosity, and we obtain the following values in units of ${\rm Gpc^{-3}~yr^{-1}}$: $0.8^{+0.1}_{-0.1}$ for high-luminosity long GRBs above $ 10^{50}~{\rm erg~s^{-1}}$, $164^{+98}_{-65}$ for low-luminosity long GRBs above $5\times 10^{46}~{\rm erg~s^{-1}}$, $1.3^{+0.4}_{-0.3}$, $1.2^{+0.4}_{-0.3}$, and $3.3^{+1.0}_{-0.8}$ above $ 10^{50}~ {\rm erg~s^{-1}}$ for short GRBs with three different merger delay models (Gaussian, log-normal, and power law), $1.9^{+2.4}_{-1.2}\times 10^4$ above $ 10^{44}~{\rm erg~s^{-1}}$ for SBOs, $ 4.8^{+3.2}_{-2.1}\times10^2 $ for normal TDEs above $10^{44}~ {\rm erg~s^{-1}}$, and $0.03^{+0.04}_{-0.02}$ above $ 10^{48}~ {\rm erg~s^{-1}}$ for TDE jets as discovered by Swift. Intriguingly, the global luminosity functions of different kinds of transients, which cover over 12 orders of magnitude, are consistent with a single power law with an index of -1.6.

Transport properties of dense liquid helium under the conditions of planet's core and cool atmosphere of white dwarfs have been investigated by using the improved centroid path-integral simulations combined with density functional theory. The self-diffusion is largely higher and the shear viscosity is notably lower predicted with the quantum mechanical description of the nuclear motion compared with the description by Newton equation. The results show that nuclear quantum effects (NQEs), which depends on the temperature and density of the matter via the thermal de Broglie wavelength and the ionization of electrons, are essential for the transport properties of dense liquid helium at certain astrophysical conditions. The Stokes-Einstein relation between diffusion and viscosity in strongly coupled regime is also examined to display the influences of NQEs.

Space debris has been posing a serious threat to human space activities and is needed to be measured and cataloged. As a new technology of space target surveillance, the measurement accuracy of DRLR (Diffuse Reflection Laser Ranging) is much higher than that of microwave radar and electro-optical measurement. Based on laser ranging data of space debris from DRLR system collected at SHAO (Shanghai Astronomical Observatory) in March-April 2013, the characteristics and precision of the laser ranging data are analyzed and its applications in OD (Orbit Determination) of space debris are discussed in this paper, which is implemented for the first time in China. The experiment indicates that the precision of laser ranging data can reach 39cm-228cm. When the data is sufficient enough (4 arcs of 3 days), the orbit accuracy of space debris can be up to 50m.

Fe exists abundantly in the universe. In particular, the dynamical structures and transport properties of warm dense Fe are crucial for understanding the evolution and structures of giant planets. In this article, we present the ionic structures, equation of states, diffusion and viscosity of Fe at two typical densities of 33.385 g/cm$^3$ and 45 g/cm$^3$ in the temperature range of 1 eV and 10 eV, giving the data by the first principles calculations using quantum Langevin molecular dynamics (QLMD). Furthermore, the validation of Stokes-Einstein (SE) relation in this regime is discussed, showing the importance of choosing the effective atomic diameter. The results remind us of the careful usage of the SE relation under extreme conditions.

A direct WIMP (Weakly Interacting Massive Particle) detector with a neutron veto system is designed to better reject neutrons. Two experimental configurations are studied in the present paper: one is for an Xe detector enclosed in a Gd-loaded scintillator and the other one is for an Xe detector placed inside a reactor neutrino detector. The Gd-doped liquid scintillator (or the neutrino detector) is used as a neutron veto device. The neutron backgrounds for the two experimental designs have been estimated using Geant4 simulations. The results show that the neutron backgrounds can decrease to O(0.1) events per year per tonne of liquid Xenon. We calculate the sensitivities to spin-independent WIMP-nucleon elastic scattering. An exposure of one tonne $\times$ year could reach a cross-section of about 6$\times$$10^{-11}$ pb.

With the development of network and the World Wide Web (WWW), the Internet has been growing and changing dramatically. More and more on-line database systems and different kinds of services are available for astronomy research. How to help users find their way through the jungle of information services becomes an important challenge. Although astronomers have been aware of the importance of interoperability and introduced the concept of Virtual Observatory as a uniform environment for future astronomical on-line resources and services, transparent access to heterogeneous on-line information is still difficult. SkyMouse is a lightweight interface for distributed astronomical on-line resources and services, which is designed and developed by us, i.e., Chinese Virtual Observatory Project. Taking advantage of screen word-capturing technology, different kinds of information systems can be queried through simple mouse actions, and results are returned in a uniform web page. SkyMouse is an easy to use application, aiming to show basic information or to create a comprehensive overview of a specific astronomical object. In this paper current status of on-line resources and services access is reviewed; system architecture, features and functions of SkyMouse are described; challenges for intelligent interface for on-line astronomical resources and services are discussed.