Precise measurements of the energy spectra of cosmic rays (CRs) show various kinds of features deviating from single power-laws, which give very interesting and important implications on their origin and propagation. Previous measurements from a few balloon and space experiments indicate the existence of spectral softenings around 10 TV for protons (and probably also for Helium nuclei). Very recently, the DArk Matter Particle Explorer (DAMPE) measurement about the proton spectrum clearly reveals such a softening with a high significance. Here we study the implications of these new measurements, as well as the groundbased indirect measurements, on the origin of CRs. We find that a single component of CRs fails to fit the spectral softening and the air shower experiment data simultaneously. In the framework of multiple components, we discuss two possible scenarios, the multiple source population scenario and the background plus nearby source scenario. Both scenarios give reasonable fits to the wide-band data from TeV to 100 PeV energies. Considering the anisotropy observations, the nearby source model is favored.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

The future of exoplanet detection lies in the mid-infrared (MIR). The MIR region contains the blackbody peak of both hot and habitable zone exoplanets, making the contrast between starlight and planet light less extreme. It is also the region where prominent chemical signatures indicative of life exist, such as ozone at 9.7 microns. At a wavelength of 4 microns the difference in emission between an Earth-like planet and a star like our own is 80 dB. However a jovian planet, at the same separation exhibits 60 dB of contrast, or only 20 dB if it is hot due to its formation energy or being close to its host star. A two dimensional nulling interferometer, made with chalcogenide glass, has been measured to produce a null of 20 dB, limited by scattered light. Measures to increase the null depth to the theoretical limit of 60 dB are discussed.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

One of the main purposes of the DArk Matter Particle Explorer (DAMPE) is to measure the cosmic ray nuclei up to several tens of TeV or beyond, whose origin and propagation remains a hot topic in astrophysics. The Plastic Scintillator Detector (PSD) on top of DAMPE is designed to measure the charges of cosmic ray nuclei from H to Fe and serves as a veto detector for discriminating gamma-rays from charged particles. We propose in this paper a charge reconstruction procedure to optimize the PSD performance in charge measurement. Essentials of our approach, including track finding, alignment of PSD, light attenuation correction, quenching and equalization correction are described detailedly in this paper after a brief description of the structure and operational principle of the PSD. Our results show that the PSD works very well and almost all the elements in cosmic rays from H to Fe are clearly identified in the charge spectrum.

DArk Matter Particle Explorer (DAMPE) is a space-borne apparatus for detecting the high-energy cosmic-rays like electrons, $\gamma$-rays, protons and heavy-ions. Plastic Scintillator Detector (PSD) is the top-most sub-detector of the DAMPE. The PSD is designed to measure the charge of incident high-energy particles and it also serves as a veto detector for discriminating $\gamma$-rays from charged particles. In this paper, PSD on-orbit calibration procedure is described, which includes five steps of pedestal, dynode correlation, response to minimum-ionizing particles (MIPs), light attenuation function and energy reconstruction. A method for reconstructing the charge of incident high energy cosmic-ray particles is introduced. The detection efficiency of each PSD strip is verified to be above 99.5%, the total efficiency of the PSD for charged particles is above 99.99%.

The Plastic Scintillator Detector (PSD) of the DArk Matter Particle Explorer (DAMPE) is designed to measure cosmic ray charge (Z) and to act as a veto detector for gamma-ray identification. In order to fully exploit the charge identification potential of the PSD and to enhance its capability to identify the gamma ray events, we develop a PSD detector alignment method. The path length of a given track in the volume of a PSD bar is derived taking into account the shift and rotation alignment corrections. By examining energy spectra of corner-passing events and fully contained events, position shifts and rotations of all PSD bars are obtained, and found to be on average about 1mm and 0.0015 radian respectively. To validate the alignment method, we introduce the artificial shifts and rotations of PSD bars in the detector simulation. These shift and rotation parameters can be recovered successfully by the alignment procedure. As a result of the PSD alignment procedure, the charge resolution of the PSD is improved from $4\%$ to $8\%$ depending on the nuclei.

The gravitational wave data of GW170817 favor the equation of state (EoS) models that predict compact neutron stars (NSs), consistent with the radius constraints from X-ray observations. Motivated by such a remarkable progress, we examine the fate of the remnants formed in NS mergers and focus on the roles of the angular momentum and the mass distribution of the binary NSs. In the mass shedding limit (for which the dimensionless angular momentum equals to the Keplerian value, i.e., $j=j_{\rm Kep}$), the adopted { seven EoS models, except H4 and ALF2,} yield supramassive NSs in more than half of the mergers. However, for $j\lesssim 0.7j_{\rm Kep}$, the presence or absence of a non-negligible fraction of supramassive NSs formed in the mergers depends sensitively on both the EoS and the mass distribution of the binary systems. The NS mergers with a total gravitational mass $\leq 2.6M_\odot$ are found to be able to shed valuable light on both the EoS model and the angular momentum of the remnants if supramassive NSs are still absent. We have also discussed the uncertainty on estimating the maximum gravitational mass of non-rotating NSs ($M_{\rm max}$) due to the unknown $j$ of the pre-collapse remnants. With the data of GW170817 and the assumption of the mass loss of $0.03M_\odot$, we have $M_{\rm max}<(2.19,~2.32)M_\odot$ (90\% confidence level) for $j=(1.0,~0.8)j_{\rm Kep}$, respectively.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

N-body simulations predict that galaxies at the Milky Way scale host a large number of dark matter (DM) subhalos. Some of these subhalos, if they are massive enough or close enough to the Earth, might be detectable in $\gamma$ rays due to the DM annihilation. 3FGL J2212.5+0703, an unidentified gamma-ray source, has been suggested to be the counterpart candidate of a DM subhalo by Bertoni et al. (2015, 2016). In this work we analyze the Fermi-LAT Pass 8 data of 3FGL J2212.5+0703 to independently test the DM subhalo hypothesis of this source. In order to suppress the possible contamination from two nearby very-bright blazars, we just take into account the front-converting gamma-rays which have better angular resolutions than that of the back-converting photons. In addition to the spatial distribution analysis, we have extended the spectrum analysis down to the energies of $\sim 100$ MeV, and thoroughly examined the variability of the emission during the past 8 years. We confirm that 3FGL J2212.5+0703 is a steady and spatially-extended gamma-ray emitter at a high confidence level. The spectrum is well consistent with that expected from DM annihilation into $b\bar{b}$. The introduction of a phenomenological LogParabola spectrum just improves the fit slightly. All these results suggest that 3FGL J2212.5+0703 could be indicative of a DM subhalo.

Photonic integrated circuits are established as the technique of choice for a number of astronomical processing functions due to their compactness, high level of integration, low losses, and stability. Temperature control, mechanical vibration and acoustic noise become controllable for such a device enabling much more complex processing than can realistically be considered with bulk optics. To date the benefits have mainly been at wavelengths around 1550 nm but in the important Mid-Infrared region, standard photonic chips absorb light strongly. Chalcogenide glasses are well known for their transparency to beyond 10000 nm, and the first results from coupler devices intended for use in an interferometric nuller for exoplanetary observation in the Mid-Infrared L band (3800-4200 nm) are presented here showing that suitable performance can be obtained both theoretically and experimentally for the first fabricated devices operating at 4000 nm.

The $^{14}$O($\alpha$,$p$)$^{17}$F reaction is one of the key reactions involved in the breakout from the hot-CNO cycle to the rp-process in type I x-ray bursts (XRBs). The resonant properties in the compound nucleus $^{18}$Ne have been investigated through resonant elastic scattering of $^{17}$F+$p$. The radioactive $^{17}$F beam was separated by the CNS Radioactive Ion Beam separator (CRIB) and bombarded a thick H$_2$ gas target at 3.6 MeV/nucleon. The recoiling light particles were measured by three ${\Delta}$E-E silicon telescopes at laboratory angles of $\theta$$_{lab}$$\approx$3$^\circ$, 10$^\circ$ and 18$^\circ$, respectively. Five resonances at $E_{x}$=6.15, 6.28, 6.35, 6.85, and 7.05 MeV were observed in the excitation functions, and their spin-parities have been determined based on an $R$-matrix analysis. In particular, $J^{\pi}$=1$^-$ was firmly assigned to the 6.15-MeV state which dominates the thermonuclear $^{14}$O($\alpha$,$p$)$^{17}$F rate below 2 GK. As well, a possible new excited state in $^{18}$Ne was observed at $E_{x}$=6.85$\pm$0.11 MeV with tentative $J$=0 assignment. This state could be the analog state of the 6.880 MeV (0$^{-}$) level in the mirror nucleus $^{18}$O, or a bandhead state (0$^+$) of the six-particle four-hole (6$p$-4$h$) band. A new thermonuclear $^{14}$O($\alpha$,$p$)$^{17}$F rate has been determined, and the astrophysical impact of multiple recent rates has been examined using an XRB model. Contrary to previous expectations, we find only modest impact on predicted nuclear energy generation rates from using reaction rates differing by up to several orders of magnitude.

The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction is thought to be one of the key breakout reactions from the hot CNO cycles to the rp-process in type I x-ray bursts. In this work, the resonant properties of the compound nucleus $^{22}$Mg have been investigated by measuring the resonant elastic scattering of $^{21}$Na+$p$. An 89 MeV $^{21}$Na radioactive beam delivered from the CNS Radioactive Ion Beam Separator bombarded an 8.8 mg/cm$^2$ thick polyethylene (CH$_{2}$)$_{n}$ target. The $^{21}$Na beam intensity was about 2$\times$10$^{5}$ pps, with a purity of about 70% on target. The recoiled protons were measured at the center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175.2${^\circ}$, 152.2${^\circ}$, and 150.5${^\circ}$ by three sets of $\Delta E$-$E$ telescopes, respectively. The excitation function was obtained with the thick-target method over energies $E_x$($^{22}$Mg)=5.5--9.2 MeV. In total, 23 states above the proton-threshold in $^{22}$Mg were observed, and their resonant parameters were determined via an $R$-matrix analysis of the excitation functions. We have made several new $J^{\pi}$ assignments and confirmed some tentative assignments made in previous work. The thermonuclear $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate has been recalculated based on our recommended spin-parity assignments. The astrophysical impact of our new rate has been investigated through one-zone postprocessing x-ray burst calculations. We find that the $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate significantly affects the peak nuclear energy generation rate, reaction fluxes, as well as the onset temperature of this breakout reaction in these astrophysical phenomena.

The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction provides a pathway for breakout from the hot CNO cycles to the $rp$-process in type I x-ray bursts. To better determine this astrophysical reaction rate, the resonance parameters of the compound nucleus $^{22}$Mg have been investigated by measuring the resonant elastic scattering of $^{21}$Na+$p$. An 89 MeV $^{21}$Na radioactive ion beam was produced at the CNS Radioactive Ion Beam Separator and bombarded an 8.8 mg/cm$^2$ thick polyethylene target. The recoiled protons were measured at scattering angles of $\theta_{c.m.}$$\approx 175 {^\circ}$ and 152${^\circ}$ by three $\Delta E$-$E$ silicon telescopes. The excitation function was obtained with a thick-target method over energies $E_x$($^{22}$Mg)=5.5--9.2 MeV. The resonance parameters have been determined through an $R$-matrix analysis. For the first time, the $J^{\pi}$ values for ten states above the alpha threshold in $^{22}$Mg have been experimentally determined in a single consistent measurement. We have made three new $J^{\pi}$ assignments and confirmed seven of the ten tentative assignments in the previous work. The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction rate has been recalculated, and the astrophysical impact of our new rate has been investigated through one-zone postprocessing x-ray burst calculations. We find that the $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate significantly affects the peak nuclear energy generation rate and the onset temperature of this breakout reaction in these phenomena.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

We consider the estimation of the regularization parameter for the simultaneous deblurring of multiple noisy images via Tikhonov regularization. We approach the problem in three ways. We first reduce the problem to a single-image deblurring for which the regularization parameter can be estimated through a classic generalized cross-validation (GCV) method. A modification of this function is used for correcting the undersmoothing typical of the original technique. With a second method, we minimize an average least-squares fit to the images and define a new GCV function. In the last approach, we use the classical $GCV$ on a single higher-dimensional image obtained by concatanating all the images into a single vector. With a reliable estimator of the regularization parameter, one can fully exploit the excellent computational characteristics typical of direct deblurring methods, which, especially for large images, makes them competitive with the more flexible but much slower iterative algorithms. The performance of the techniques is analyzed through numerical experiments. We find that under the independent homoscedastic and Gaussian assumptions made on the noise, the three approaches provide almost identical results with the first single image providing the practical advantage that no new software is required and the same image can be used with other deblurring algorithms.