In this brief report, we investigate the existence of 4-dimensional static spherically symmetric black holes (BHs) in the Einstein-complex-scalar-Gauss-Bonnet (EcsGB) gravity with an arbitrary potential $V(\phi)$ and a coupling $f(\phi)$ between the scalar field $\phi$ and the Gauss-Bonnet (GB) term. We find that static regular BH solutions with complex scalar hairs do not exist. This conclusion does not depend on the coupling between the GB term and the scalar field, nor on the scalar potential $V(\phi)$ and the presence of a cosmological constant $\Lambda$ (which can be either positive or negative), as longer as the scalar field remains complex and is regular across the horizon.

As an important source for large geomagnetic storms, an "ICME-in-sheath" is a completely shocked interplanetary coronal mass ejection (ICME) stuck in the sheath between a shock and host ejecta. Typical characteristics are identified from coordinated multi-sets of observations: (1) it is usually short in duration and lasts a few hours at 1 AU; (2) its solar wind parameters, in particular the magnetic field, seem to keep enhanced for a large range of distances; and (3) common ICME signatures are often lost. The host ejecta could be a single ICME or a complex ejecta, being fast enough to drive a shock. These results clarify previous misinterpretations of this phenomenon as a normal part of a sheath region. The "ICME-in-sheath" phenomenon, together with a preconditioning effect, produced an extreme set of the magnetic field, speed and density near 1 AU in the 2012 July 23 case, all around their upper limits at the same time. This is probably the most extreme solar wind driving at 1 AU and enables us to estimate the plausible upper limit for geomagnetic storm activity. With an appropriate modification in the southward field, we suggest that a geomagnetic storm with a minimum $D_{\rm st}$ of about $-2000$ nT could occur in principle. The magnetopause would be compressed to about 3.3 Earth radii from the Earth's center, well inside the geosynchronous orbit.

This article presents an analysis and the resulting limits on light dark matter inelastically scattering off of electrons, and on dark photon and axion-like particle absorption, using a second-generation SuperCDMS high-voltage eV-resolution detector. The 0.93 gram Si detector achieved a 3 eV phonon energy resolution; for a detector bias of 100 V, this corresponds to a charge resolution of 3% of a single electron-hole pair. The energy spectrum is reported from a blind analysis with 1.2 gram-days of exposure acquired in an above-ground laboratory. With charge carrier trapping and impact ionization effects incorporated into the dark matter signal models, the dark matter-electron cross section $\bar{\sigma}_{e}$ is constrained for dark matter masses from 0.5--$10^{4} $MeV$/c^{2}$; in the mass range from 1.2--50 eV$/c^{2}$ the dark photon kinetic mixing parameter $\varepsilon$ and the axioelectric coupling constant $g_{ae}$ are constrained. The minimum 90% confidence-level upper limits within the above mentioned mass ranges are $\bar{\sigma}_{e}\,=\,8.7\times10^{-34}$ cm$^{2}$, $\varepsilon\,=\,3.3\times10^{-14}$, and $g_{ae}\,=\,1.0\times10^{-9}$.

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.

In big data era, the special data with rare characteristics may be of great significations. However, it is very difficult to automatically search these samples from the massive and high-dimensional datasets and systematically evaluate them. The DoPS, our previous work [2], provided a search method of rare spectra with double-peaked profiles from massive and high-dimensional data of LAMOST survey. The identification of the results is mainly depended on visually inspection by astronomers. In this paper, as a follow-up study, a new lattice structure named SVM-Lattice is designed based on SVM(Support Vector Machine) and FCL(Formal Concept Lattice) and particularly applied in the recognition and evaluation of rare spectra with double-peaked profiles. First, each node in the SVM-Lattice structure contains two components: the intents are defined by the support vectors trained by the spectral samples with the specific characteristics, and the relevant extents are all the positive samples classified by the support vectors. The hyperplanes can be extracted from every lattice node and used as classifiers to search targets by categories. A generalization and specialization relationship is expressed between the layers, and higher layers indicate higher confidence of targets. Then, including a SVM-Lattice building algorithm, a pruning algorithm based on association rules, and an evaluation algorithm, the supporting algorithms are provided and analysed. Finally, for the recognition and evaluation of spectra with double-peaked profiles, several data sets from LAMOST survey are used as experimental dataset. The results exhibit good consistency with traditional methods, more detailed and accurate evaluations of classification results, and higher searching efficiency than other similar methods.

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.

In this paper, we systematically study spherically symmetric static spacetimes in the framework of Einstein-aether theory, and pay particular attention to the existence of black holes (BHs). In the theory, one scalar and one vector two gravitational modes additionally appear. To avoid the vacuum gravi-\v{C}erenkov radiation, they must all propagate with speeds greater than or at least equal to the speed of light. In the spherical case, only the scalar mode is relevant, so BH horizons are defined by this mode. In the present studies we first clarify several subtle issues. In particular, we find that, out of the five non-trivial field equations, only three are independent, so the problem is well-posed, as now there are only three unknown functions, two, $F(r)$ and $B(r)$, describe the spacetime and one, $A(r)$, describes the aether field. In addition, the two second-order differential equations for $A$ and $F$ are independent of $B$, and once they are found, $B$ is given simply by an algebraic expression of $F,\; A$ and their derivatives. To simplify the problem further, we explore the symmetry of field redefinitions, and work first with the redefined metric and aether field, and then obtain the physical ones by the inverse transformations. These clarifications significantly simplify the computational labor, which is important, as the problem is highly involved mathematically. In fact, it is exactly because of these, we find various numerical BH solutions with an accuracy that is at least two orders higher than previous ones. More important, these BH solutions are the only ones that satisfy the self-consistent conditions and meantime are consistent with all the observational constraints obtained so far. The locations of universal horizons are also identified, together with several other observationally interesting quantities.

Rotation and mass loss are crucially interlinked properties of massive stars, strongly affecting their evolution and ultimate fate. Massive stars rotating near their breakup limit shed mass centrifugally, creating Be stars with circumstellar disks and possibly driving outbursts. Using the MESA stellar evolution code, we examine the effects of efficient angular momentum transport on the main sequence and post-main sequence rotational evolution of massive stars. In rapid rotators, angular momentum transported from the contracting core to the expanding envelope can spin up the surface layers past the breakup rate, particularly for stars near (or beyond) the end of the main sequence and in low-metallicity environments. We also demonstrate that centrifugal instabilities could arise in rapidly rotating massive stars, potentially triggering the S Doradus outbursts observed in luminous blue variable stars. Prior mass accretion from a binary companion increases both the likelihood and the intensity of centrifugal mass loss. We discuss implications for massive stellar evolution, Be stars, and luminous blue variables.

In modeling the X-ray spectra of active galactic nuclei (AGNs), the inclination angle is a parameter that can play an important role in analyzing the X-ray spectra of AGN, but it has never been studied in detail. We present a broadband X-ray spectral analysis of the joint NuSTAR-XMM-Newton observations of 13 sources with [OIII] measured inclinations determined by Fischer et al. 2013. By freezing the inclination angles at the [OIII] measured values when modeling the observations, the spectra are well fitted and the geometrical properties of the obscuring structure of the AGNs are slightly better constrained than those fitted when the inclination angles are left free to vary. We also test if one could freeze the inclinations at other specific angles in fitting the AGN X-ray spectra as commonly did in the literatures. We find that one should always let the inclination angle free to vary in modeling the X-ray spectra of AGNs, while fixing the inclination angle at [OIII] measured values and fixing the inclination angle at 60$^\circ$ also present correct fits of the sources in our sample.Correlations between the covering factor and the average column density of the obscuring torus with respect to the Eddington ratio are also measured, suggesting that the distribution of the material in the obscuring torus is regulated by the Eddington ratio, which is in agreement with previous studies. In addition, no geometrical correlation is found between the narrow line region of the AGN and the obscuring torus, suggesting that the geometry might be more complex than what is assumed in the simplistic unified model.

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.

Testing gravity by binary pulsars nowadays becomes a key issue. Screened modified gravity is a kind of scalar-tensor theory with screening mechanism in order to satisfy the tight Solar System tests. In this paper, we investigate how the screening mechanism affects the orbital dynamics of binary pulsars, and calculate in detail the five post-Keplerian (PK) parameters in this theory. These parameters differ from those of general relativity (GR), and the differences are quantified by the scalar charges, which lead to the dipole radiation in this theory. We combine the observables of PK parameters for the ten binary pulsars, respectively, to place the constraints on the scalar charges and possible deviations from GR. The dipole radiation in the neutron star (NS) - white dwarf (WD) binaries leads to more stringent constraints on deviations from GR. The most constraining systems for the scalar charges of NS and WD are PSR~B1913$+$16 and PSR~J1738$+$0333, respectively. The results of all tests exclude significant strong-field deviations and show good agreement with GR.

Using the multiconfiguration Dirac-Hartree-Fock and the relativistic configuration interaction methods, a consistent set of transition energies and radiative transition data for the lowest 546 (623, 701, 745) states of the $3p^4 3d$, $3s 3p^2 3d^2$, $3s 3p^3 4p$, $3s 3p^4$, $3s^2 3d^3$, $3s^2 3p^2 3d$, $3s^2 3p^2 4d$, $3s^2 3p^2 4s$, $3p^3 3d^2$, $3p^5$, $3s 3p 3d^3$, $3s 3p^3 3d$, $3s 3p^3 4s$, $3s^2 3p 3d^2$, %$3s^2 3p^2 4f$, $3s^2 3p^2 4p$, $3s^2 3p^3$ configurations in Mn~XI (Fe~XII, Co~XIII, Ni~XIV) is provided. The comparison between calculated excitation energies for the $n=4$ states and available experimental values for Fe XII indicate that the calculations are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Lines from these states are prominent in the soft X-rays. With the present calculations, several recent new identifications are confirmed. Other identifications involving $3p^2 4d$ levels in Fe~XII that were found questionable are discussed and a few new assignments are recommended. As some $n=4$ states of the other ions also show large discrepancies between experimental and calculated energies, we reassess their identification. The present study provides highly accurate atomic data for the $n=4$ states of P-like ions of astrophysical interest, for which experimental data are scarce.

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.

We present a detailed spectral analysis of the black hole candidate MAXI J1836-194. The source was caught in the intermediate state during its 2011 outburst by Suzaku and RXTE. We jointly fit the X-ray data from these two missions using the relxill model to study the reflection component, and a steep inner emissivity profile indicating a compact corona as the primary source is required in order to achieve a good fit. In addition, a reflection model with a lamp-post configuration (relxilllp), which is normally invoked to explain the steep emissivity profile, gives a worse fit and is excluded at 99% confidence level compared to relxill. We also explore the effect of the ionization gradient on the emissivity profile by fitting the data with two relativistic reflection components, and it is found that the inner emissivity flattens. These results may indicate that the ionization state of the disc is not constant. All the models above require a supersolar iron abundance higher than 4.5. However, we find that the high-density version of reflionx can describe the same spectra even with solar iron abundance well. A moderate rotating black hole (a* = 0.84-0.94) is consistently obtained by our models, which is in agreement with previously reported values.

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.

The synchrotron radiation from secondary electrons and positrons (SEPs) generated by hadronic interactions in the shock of supernova remnant (SNR) could be a distinct evidence of cosmic ray (CR) production in SNR shocks. Here we provide a method where the observed gamma-ray flux from SNRs, created by pion decays, is directly used to derive the SEP distribution and hence the synchrotron spectrum. We apply the method to three gamma-ray bright SNRs. In the young SNR RX J1713.7-3946, if the observed GeV-TeV gamma-rays are of hadronic origin and the magnetic field in the SNR shock is $B\gtrsim 0.5$mG, the SEPs may produce a spectral bump at $10^{-5}-10^{-2}$eV, exceeding the predicted synchrotron component of the leptonic model, and a soft spectral tail at $\gtrsim 100$keV, distinct from the hard spectral slope in the leptonic model. In the middle-aged SNRs IC443 and W44, if the observed gamma-rays are of hadronic origin, the SEP synchrotron radiation with $B\sim 400 - 500 \mu$G can well account for the observed radio flux and spectral slopes, supporting the hadronic origin of gamma-rays. Future microwave to far-infrared and hard X-ray (>100keV) observations are encouraged to constraining the SEP radiation and the gamma-ray origin in SNRs.

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.

We review and discuss insights on ideal magnetohydrodynamic (MHD) instabilities that can play a role in destabilizing solar coronal flux rope structures. For single flux ropes, failed or actual eruptions may result from internal or external kink evolutions, or from torus unstable configurations. We highlight recent findings from 3D magnetic field reconstructions and simulations where kink and torus instabilities play a prominent role. For interacting current systems, we critically discuss different routes to coronal dynamics and global eruptions, due to current channel coalescence or to tilt-kink scenarios. These scenarios involve the presence of two nearby current channels and are clearly distinct from the popular kink or torus instability. Since the solar corona is pervaded with myriads of magnetic loops -- creating interacting flux ropes typified by parallel or antiparallel current channels as exemplified in various recent observational studies -- coalescence or tilt-kink evolutions must be very common for destabilizing adjacent flux rope systems. Since they also evolve on ideal MHD timescales, they may well drive many sympathetic eruptions witnessed in the solar corona. Moreover, they necessarily lead to thin current sheets that are liable to reconnection. We review findings from 2D and 3D MHD simulations for tilt and coalescence evolutions, as well as on particle acceleration aspects derived from computed charged particle motions in tilt-kink disruptions and coalescing flux ropes. The latter were recently studied in two-way coupled kinetic-fluid models, where the complete phase-space information of reconnection is incorporated.

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.

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.

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

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.

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.

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.

In this paper, we study the gravitational waveforms, polarizations and radiation powers of the first relativistic triple systems PSR J0337 + 1715, observed in 2014, by using the post-Newtonian approximations to their lowest order. Although they cannot be observed either by current or next generation of the detectors, they do provide useful information to test different theories of gravity. In particular, we carry out the studies in three different theories, general relativity (GR), Brans-Dicke (BD) gravity, and Einstein-aether ($\ae$) theory. The tensor modes $h_{+}$ and $h_{\times}$ exist in all three theories and have almost equal amplitudes. Their frequencies are all peaked at two locations, $ f^{(+, \times)}_1 \simeq 0.0686656 \mu$Hz and $ f^{(+, \times)}_2 \simeq 14.2138 \mu$Hz, which are about twice of the outer and inner orbital frequencies of the triple system, as predicted in GR. In $\ae$-theory, all the six polarization modes are different from zero, but the breathing $h_b$ and longitudinal $h_L$ modes are not independent and also peaked at two frequencies, but at the frequencies, $ f^{(b, L, \ae)}_{1} \simeq 0.0457771 \mu$Hz and $f^{(b, L, \ae)}_{ 2} \simeq 7.09545 \mu$Hz. A similar phenomenon is also observed in BD gravity, in which only the three modes $h_{+},\; h_{\times}$ and $h_{b}$ exit. We also study the radiation powers, and find that the quadrupole emission in each of the three theories has almost the same amplitude, but the dipole emission can be as big as the quadrupole emission in $\ae$-theory. This can provide a very promising window to obtain severe constraints on $\ae$-theory by the multi-band gravitational wave astronomy.

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.

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.

An unexpected strong geomagnetic storm occurred on 2018 August 26, which was caused by a slow coronal mass ejection (CME) from a gradual eruption of a large quiet-region filament. We investigate the eruption and propagation characteristics of this CME in relation to the strong geomagnetic storm with remote sensing and in situ observations. Coronal magnetic fields around the filament are extrapolated and compared with EUV observations. We determine the propagation direction and tilt angle of the CME flux rope near the Sun using a graduated cylindrical shell (GCS) model and the Sun-to-Earth kinematics of the CME with wide-angle imaging observations from STEREO A. We reconstruct the flux-rope structure using a Grad-Shafranov technique based on the in situ measurements at the Earth and compare it with those from solar observations and the GCS results. Our conclusions are as follows: (1) the eruption of the filament was unusually slow and occurred in the regions with relatively low critical heights of the coronal field decay index; (2) the axis of the CME flux rope rotated in the corona as well as in interplanetary space, which tended to be aligned with the local heliospheric current sheet; (3) the CME was bracketed between slow and fast solar winds, which enhanced the magnetic field inside the CME at 1 AU; (4) the geomagnetic storm was caused by the enhanced magnetic field and a southward orientation of the flux rope at 1 AU from the rotation of the flux rope.

Future direct searches for low-mass dark matter particles with germanium detectors, such as SuperCDMS SNOLAB, are expected to be limited by backgrounds from radioactive isotopes activated by cosmogenic radiation inside the germanium. There are limited experimental data available to constrain production rates and a large spread of theoretical predictions. We examine the calculation of expected production rates, and analyze data from the second run of the CDMS low ionization threshold experiment (CDMSlite) to estimate the rates for several isotopes. We model the measured CDMSlite spectrum and fit for contributions from tritium and other isotopes. Using the knowledge of the detector history, these results are converted to cosmogenic production rates at sea level. The production rates in atoms/(kg$\cdot$day) are 74$\pm$9 for $^3$H, 1.5$\pm$0.7 for $^{55}$Fe, 17$\pm$5 for $^{65}$Zn, and 30$\pm$18 for $^{68}$Ge.

In order to have a comprehensive view of the propagation and evolution of coronal mass ejections (CMEs) from the Sun to deep interplanetary space beyond 1 au, we carry out a kinematic analysis of 7 CMEs in solar cycle 23. The events are required to have coordinated coronagraph observations, interplanetary type II radio bursts, and multi-point in-situ measurements at the Earth and Ulysses. A graduated cylindrical shell model, an analytical model without free parameters and a magnetohydrodynamic model are used to derive CME kinematics near the Sun, to quantify the CME/shock propagation in the Sun-Earth space, and to connect in-situ signatures at the Earth and Ulysses, respectively. We find that each of the 7 CME-driven shocks experienced a major deceleration before reaching 1 au and thereafter propagated with a gradual deceleration from the Earth to larger distances. The resulting CME/shock propagation profile for each case is roughly consistent with all the data, which verifies the usefulness of the simple analytical model for CME/shock propagation in the heliosphere. The statistical analysis of CME kinematics indicates a tendency that the faster the CME, the larger the deceleration, and the shorter the deceleration time period within 1 au. For several of these events, the associated geomagnetic storms were mainly caused by the southward magnetic fields in the sheath region. In particular, the interaction between a CME-driven shock and a preceding ejecta significantly enhanced the preexisting southward magnetic fields and gave rise to a severe complex geomagnetic storm.

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 joint NuSTAR and XMM-Newton spectral analysis in the 0.6-70 keV band of three candidate Compton thick (CT-) AGN selected in the 100-month Swift-BAT catalog. These objects were previously classified as CT-AGNs based on low quality Swift-XRT and Swift-BAT data, and had soft photon indices (Gamma>2.2) that suggested a potential overestimation of the line of sight column density. Thanks to the high-quality NuSTAR and XMM-Newton data we were able to determine that in all three objects the photon index was significantly overestimated, and two out of three sources are reclassified from CT to Compton thin, confirming a previously observed trend, i.e., that a significant fraction of BAT-selected, candidate CT-AGNs with poor soft X-ray data are reclassified as Compton thin when the NuSTAR data are added to the fit. Finally, thanks to both the good XMM-Newton spatial resolution and the high NuSTAR and XMM-Newton spectral quality, we found that the third object in our sample was associated to the wrong counterpart: the correct one, 2MASX J10331570+5252182, has redshift z=0.14036, which makes it one of the very few candidate CT-AGNs in the 100-month BAT catalog detected at z>0.1, and a rare CT quasar.

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.

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.

Type Ia supernovae (SNe Ia) play key roles in revealing the accelerating expansion of the universe, but our knowledge about their progenitors is still very limited. Here we report the discovery of a rigid dichotomy in circumstellar (CS) environments around two subclasses of type Ia supernovae (SNe Ia) as defined by their distinct photospheric velocities. For the SNe Ia with high photospheric velocities (HV), we found a significant excess flux in blue light during 60-100 days past maximum, while this phenomenon is absent for SNe with normal photospheric velocity (Normal). This blue excess can be attributed to light echoes by circumstellar dust located at a distance of about 1-3x10^{17} cm from the HV subclass. Moreover, we also found that the HV SNe Ia show systematically evolving Na I absorption line by performing a systematic search of variable Na I absorption lines in spectra of all SNe Ia, whereas this evolution is rarely seen in Normal ones. The evolving Na I absorption can be modeled in terms of photoionization model, with the location of the gas clouds at a distance of about 2x10^{17} cm, in striking agreement with the location of CS dust inferred from B-band light curve excess. These observations show clearly that the progenitors of HV and Normal subclasses are systematically different, suggesting that they are likely from single and double degenerate progenitor systems, respectively.

We present the first limits on inelastic electron-scattering dark matter and dark photon absorption using a prototype SuperCDMS detector having a charge resolution of 0.1 electron-hole pairs (CDMS HVeV, a 0.93 gram CDMS HV device). These electron-recoil limits significantly improve experimental constraints on dark matter particles with masses as low as 1 MeV/$\mathrm{c^2}$. We demonstrate a sensitivity to dark photons competitive with other leading approaches but using substantially less exposure (0.49 gram days). These results demonstrate the scientific potential of phonon-mediated semiconductor detectors that are sensitive to single electronic excitations.

The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) searches for interactions between dark matter particles and germanium nuclei in cryogenic detectors. The experiment has achieved a low energy threshold with improved sensitivity to low-mass (<10 GeV/c$^2$) dark matter particles. We present an analysis of the final CDMSlite data set, taken with a different detector than was used for the two previous CDMSlite data sets. This analysis includes a data "salting" method to protect against bias, improved noise discrimination, background modeling, and the use of profile likelihood methods to search for a dark matter signal in the presence of backgrounds. We achieve an energy threshold of 70 eV and significantly improve the sensitivity for dark matter particles with masses between 2.5 and 10 GeV/c$^2$ compared to previous analyses. We set an upper limit on the dark matter-nucleon scattering cross section in germanium of 5.4$\times$10$^{-42}$ cm$^2$ at 5 GeV/c$^2$, a factor of $\sim$2.5 improvement over the previous CDMSlite result.

We conduct forward-modeling analysis based on our 2.5 dimensional magnetohydrodynamics (MHD) simulation of magnetic flux rope (MFR) formation and eruption driven by photospheric converging motion. The current sheet (CS) evolution during the MFR formation and eruption process in our MHD simulation can be divided into four stages. The current sheet (CS) evolution during the MFR formation and eruption process in our MHD simulation can be divided into four stages. The first stage shows the CS forming and gradually lengthening. Resistive instabilities that disrupt the CS mark the beginning of the second stage. Magnetic islands disappear in the third stage and reappear in the fourth stage. Synthetic images and light curves of the seven Solar Dynamics Observatory/Atmospheric Imaging Assembly (AIA) channels, i.e., 9.4 nm, 13.1 nm, 17.1 nm, 19.3 nm, 21.1 nm, 30.4 nm, and 33.5 nm, and the 3-25 keV thermal X-ray are obtained with forward-modeling analysis. The loop-top source and the coronal sources of the soft X-ray are reproduced in forward modeling. The light curves of the seven SDO/AIA channels start to rise once resistive instabilities develop. The light curve of the 3-25 keV thermal X-ray starts to go up when the reconnection rate reaches one of its peaks. Quasiperiodic pulsations (QPPs) appear twice in the SDO/AIA 17.1 nm, 21.1 nm, and 30.4 nm channels, corresponding to the period of chaotic (re)appearance and CS-guided displacements of the magnetic islands. QPPs appear once in the SDO/AIA 9.4 nm and 33.5 nm channels after the disruption of the CS by resistive instabilities and in the 19.3 nm channel when the chaotic motion of the magnetic islands reappears.

The Super Cryogenic Dark Matter Search experiment (SuperCDMS) at the Soudan Underground Laboratory studied energy loss associated with Frenkel defect formation in germanium crystals at mK temperatures using in situ $^{210}$Pb sources. We examine the spectrum of $^{206}$Pb nuclear recoils near its expected 103 keV endpoint energy and determine an energy loss of $\left(6.08\pm0.18\right)$ %, which we attribute to defect formation. From this result and using TRIM simulations, we extract the first experimentally determined average displacement threshold energy of $\left(19.7^{+0.6}_{-0.5}\right)$ eV for germanium. This has implications for the analysis thresholds of future germanium-based dark matter searches.

We investigate the prompt emission and the afterglow properties of short duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multi-wavelength facilities including the GTC 10.4m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL/SPI/ACS, Swift/BAT and Fermi/GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB130603B, our new afterglow photometric data constraints the pre jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new as well as previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for 8 other sGRBs we are able to; place tight constraints on the non-detection of the afterglow, host galaxy or any underlying kilonova emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs.

It is paramount from both scientific and societal perspectives to understand the generation of extreme space weather. We discuss the formation of solar superstorms based on a comparative study of the 2012 July 23 and 2017 July 23 eruptions. The first one is Carrington-class, and the second could rival the 1989 March event that caused the most intense geomagnetic storm of the space age. Observations of these events in the historically weak solar cycle 24 indicate that a solar superstorm can occur in any solar cycle and at any phase of the cycle. Recurrent patterns are identified in both cases, including the long-lived eruptive nature of the active region, a complex event composed of successive eruptions from the same active region, and in-transit interaction between the successive eruptions resulting in exceptionally strong ejecta magnetic fields at 1 AU. Each case also shows unique characteristics. Preconditioning of the upstream solar wind leading to unusually high solar wind speeds at 1 AU is observed in the first case whereas not in the latter. This may suggest that the concept of "preconditioning" appears to be necessary for making a Carrington-class storm. We find a considerable deflection by nearby coronal holes in the second case but not in the first. On the basis of these results, we propose a hypothesis for further investigation that superstorms are "perfect storms" in nature, i.e., a combination of circumstances that results in an event of unusual magnitude. Historical records of some extreme events seem to support our hypothesis.

Gravitationally bound hierarchies containing three or more components are very common in our Universe. In this paper we study {\em periodic} gravitational wave (GW) form, their polarizations, response function, its Fourier transform, and energy loss rate of a triple system through three different channels of radiation, the scalar, vector and tensor modes, in Einstein-aether theory of gravity. In the weak-field approximations and with the recently obtained constraints of the theory, we first analyze the energy loss rate of a binary system, and find that the dipole contributions from the scalar and vector modes could be of the order of ${\cal{O}}\left(c_{14}\right){\cal{O}}\left(G_Nm/d\right)^2$, where $c_{14} \; (\equiv c_{1} + c_{4})$ is constrained to $c_{14} \lesssim {\cal{O}}\left(10^{-5}\right)$ by current observations, where $c_i$'s are the four coupling constants of the theory. On the other hand, the "strong-field" effects for a binary system of neutron stars are about six orders lower than that of GR. So, in this paper we ignore these "strong-field" effects and first develop the general formulas to the lowest post-Newtonian order, by taking the coupling of the aether field with matter into account. Within this approximation, we find that the scalar breather mode and the scalar longitudinal mode are all suppressed by a factor of ${\cal{O}}\left(c_{14}\right)$ with respect to the transverse-traceless modes ($h_{+}$ and $h_{\times}$), while the vectorial modes $(h_{X}$ and $h_{Y}$) are suppressed by a factor of $c_{13} \lesssim {\cal{O}}\left(10^{-15}\right)$. Applying the general formulas to a triple system with periodic orbits, we find that the corresponding GW form, response function, and its Fourier transform depend sensitively on both the configuration of the triple system and their orientations with respect to the detectors.

In this paper, we present extensive optical and ultraviolet (UV) observations of the type Ia supernova (SN Ia) 2013gs discovered during the Tsinghua-NAOC Transient Survey. The photometric observations in the optical show that the light curves of SN 2013gs is similar to that of normal SNe Ia, with an absolute peak magnitude of $M_{B}$ = $-$19.25 $\pm$ 0.15 mag and a post-maximum decline rate $\Delta$m$_{15}$(B) = 1.00 $ \pm $ 0.05 mag. \emph{Gehrels Swift} UVOT observations indicate that SN 2013gs shows unusually strong UV emission (especially in the $uvw1$ band) at around the maximum light (M$_{uvw1}$ $\sim$ $-$18.9 mag). The SN is characterized by relatively weak Fe~{\sc ii} {\sc iii} absorptions at $\sim$ 5000{\AA} in the early spectra and a larger expansion velocity ($v_{Si}$ $\sim$ 13,000 km s$^{-1}$ around the maximum light) than the normal-velocity SNe Ia. We discuss the relation between the $uvw1-v$ color and some observables, including Si~{\sc ii} velocity, line strength of Si~{\sc ii} $\lambda$6355, Fe~{\sc ii}/{\sc iii} lines and $\Delta m_{15}$(B). Compared to other fast-expanding SNe Ia, SN 2013gs exhibits Si and Fe absorption lines with similar strength and bluer $uvw1-v$ color. We briefly discussed the origin of the observed UV dispersion of SNe Ia.

We present the analysis of a sample of 35 candidate Compton thick (CT-) active galactic nuclei (AGNs) selected in the nearby Universe (average redshift <z>~0.03) with the Swift-BAT 100-month survey. All sources have available NuSTAR data, thus allowing us to constrain with unprecedented quality important spectral parameters such as the obscuring torus line-of-sight column density (N_{H, z}), the average torus column density (N_{H, tor}) and the torus covering factor (f_c). We compare the best-fit results obtained with the widely used MyTorus (Murphy et al. 2009) model with those of the recently published borus02 model (Balokovic et al. 2018) used in the same geometrical configuration of MyTorus (i.e., with f_c=0.5). We find a remarkable agreement between the two, although with increasing dispersion in N_{H, z} moving towards higher column densities. We then use borus02 to measure f_c. High-f_c sources have, on average, smaller offset between N_{H, z} and N_{H, tor} than low-f_c ones. Therefore, low f_c values can be linked to a "patchy torus" scenario, where the AGN is seen through an over-dense region in the torus, while high-f_c objects are more likely to be obscured by a more uniform gas distribution. Finally, we find potential evidence of an inverse trend between f_c and the AGN 2-10 keV luminosity, i.e., sources with higher f_c values have on average lower luminosities.

We present the 2--78 keV spectral analysis of the deep NuSTAR and XMM-Newton observation of a nearby Seyfert 2 galaxy, ESO 116-G018, which is selected as a candidate Compton-thick (CT-) active galactic nucleus (AGN) based on a previous Chandra-Swift/BAT study. Through our analysis, the source is for the first time confirmed to be a CT-AGN at a $>$3 $\sigma$ confidence level, with the 'line-of-sight' column density N$_{\rm H,Z}$ = [2.46-2.76] $\times$ $10^{24}$ cm$^{-2}$. The 'global average' column density of the obscuring torus is N$_{\rm H,S}$ = [0.46-0.62] $\times$ $10^{24}$ cm$^{-2}$, which suggests a clumpy, rather than uniform, distribution of the obscuring material surrounding the accreting supermassive black hole. The excellent-quality data given by the combined NuSTAR and XMM-Newton observations enable us to produce a strong constraint on the covering factor of the torus of ESO 116-G018, which is found to be $f_c$ = [0.13-0.15]. We also estimate the bolometric luminosity from the broad-band X-ray spectrum, being L$\rm _{bol}$ = [2.57-3.41] $\times$ 10$^{44}$ erg s$^{-1}$.

A powerful coronal mass ejection (CME) occurred on 2017 September 10 near the end of the declining phase of the historically weak solar cycle 24. We obtain new insights concerning the geometry and kinematics of CME-driven shocks in relation to their heliospheric impacts from the optimal, multi-spacecraft observations of the eruption. The shock, which together with the CME driver can be tracked from the early stage to the outer corona, shows a large oblate structure produced by the vast expansion of the ejecta. The expansion speeds of the shock along the radial and lateral directions are much larger than the translational speed of the shock center, all of which increase during the flare rise phase, peak slightly after the flare maximum and then decrease. The near simultaneous arrival of the CME-driven shock at the Earth and Mars, which are separated by 156.6$^{\circ}$ in longitude, is consistent with the dominance of expansion over translation observed near the Sun. The shock decayed and failed to reach STEREO A around the backward direction. Comparison between ENLIL MHD simulations and the multi-point in situ measurements indicates that the shock expansion near the Sun is crucial for determining the arrival or non-arrival and space weather impact at certain heliospheric locations. The large shock geometry and kinematics have to be taken into account and properly treated for accurate predictions of the arrival time and space weather impact of CMEs.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.

We present the combined NuSTATR and XMM-Newton 0.6-79 keV spectral analysis of a Seyfert 2 galaxy, NGC 1358, which we selected as a candidate Compton thick (CT-) active galactic nucleus (AGN) on the basis of previous Swift/BAT and Chandra studies. According to our analysis, NGC 1358 is confirmed to be a CT-AGN using physical motivated models, at >3 $\sigma$ confidence level. Our best-fit shows that the column density along the 'line-of-sight' of the obscuring material surrounding the accreting super-massive black hole is N$\rm _H$ = [1.96--2.80] $\times$ 10$^{24}$ cm$^{-2}$. The high-quality data from NuSTAR gives the best constraints on the spectral shape above $\sim$10 keV to date on NGC 1358. Moreover, by combining NuSTAR and XMM-Newton data, we find that the obscuring torus has a low covering factor ($f_c$ <0.17), and the obscuring material is distributed in clumps, rather than uniformly. We also derive an estimate of NGC 1358's Eddington ratio, finding it to be $\lambda_{\rm Edd}$ $\sim$$4.7_{-0.3}^{+0.3}$ $\times$ 10$^{-2}$, which is in acceptable agreement with previous measurements. Finally, we find no evidence of short-term variability, over a $\sim$100 ks time-span, in terms of both 'line-of-sight' column density and flux.

We present an early-phase $g$-band light curve and visual-wavelength spectra of the normal Type Ia supernova (SN) 2013gy. The light curve is constructed by determining the appropriate S-corrections to transform KAIT natural-system $B$- and $V$-band photometry and Carnegie Supernova Project natural-system $g$-band photometry to the Pan-STARRS1 $g$-band natural photometric system. A Markov Chain Monte Carlo calculation provides a best-fit single power-law function to the first ten epochs of photometry described by an exponent of $2.16^{+0.06}_{-0.06}$ and a time of first light of MJD 56629.4$^{+0.1}_{-0.1}$, which is $1.93^{+0.12}_{-0.13}$ days (i.e., $<48$~hr) before the discovery date (2013 December 4.84 UT) and $-19.10^{+0.12}_{-0.13}$ days before the time of $B$-band maximum (MJD 56648.5$\pm0.1$). The estimate of the time of first light is consistent with the explosion time inferred from the evolution of the Si II $\lambda$6355 Doppler velocity. Furthermore, discovery photometry and previous nondetection limits enable us to constrain the companion radius down to $R_c \leq 4\,R_{\odot}$. In addition to our early-time constraints, we use a deep +235 day nebular-phase spectrum from Magellan/IMACS to place a stripped H-mass limit of $< 0.018\,M_{\odot}$. Combined, these limits effectively rule out H-rich nondegenerate companions.

We present extensive spectroscopic observations for one of the closest type Ia supernovae (SNe Ia), SN 2014J discovered in M82, ranging from 10.4 days before to 473.2 days after B-band maximum light. The diffuse interstellar band (DIB) features detected in a high-resolution spectrum allow an estimate of line-of-sight extinction as Av=1.9+/-0.6 mag. Spectroscopically, SN 2014J can be put into the high-velocity (HV) subgroup in Wang's classification with a velocity of Si~II 6355 at maximum light as about 12200 km/s, but has a low velocity gradient (LVG, following Benetti's classification) as 41+/-2 km/s/day, which is inconsistent with the trend that HV SNe Ia generally have larger velocity gradients. We find that the HV SNe Ia with LVGs tend to have relatively stronger Si III (at ~4400 Angstrom) absorptions in early spectra, larger ratios of S II 5468 to S II 5640, and weaker Si II 5972 absorptions compared to their counterparts with similar velocities but high velocity gradients. This shows that the HV+LVG subgroup of SNe Ia may have intrinsically higher photospheric temperature, which indicates that their progenitors may experience more complete burning in the explosions relative to the typical HV SNe Ia.

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 Cryogenic Dark Matter Search (CDMS II) experiment aims to detect dark matter particles that elastically scatter from nuclei in semiconductor detectors. The resulting nuclear-recoil energy depositions are detected by ionization and phonon sensors. Neutrons produce a similar spectrum of low-energy nuclear recoils in such detectors, while most other backgrounds produce electron recoils. The absolute energy scale for nuclear recoils is necessary to interpret results correctly. The energy scale can be determined in CDMS II silicon detectors using neutrons incident from a broad-spectrum $^{252}$Cf source, taking advantage of a prominent resonance in the neutron elastic scattering cross section of silicon at a recoil (neutron) energy near 20 (182) keV. Results indicate that the phonon collection efficiency for nuclear recoils is $4.8^{+0.7}_{-0.9}$% lower than for electron recoils of the same energy. Comparisons of the ionization signals for nuclear recoils to those measured previously by other groups at higher electric fields indicate that the ionization collection efficiency for CDMS II silicon detectors operated at $\sim$4 V/cm is consistent with 100% for nuclear recoils below 20 keV and gradually decreases for larger energies to $\sim$75% at 100 keV. The impact of these measurements on previously published CDMS II silicon results is small.

We report the observations of a moderate but relatively intense geo-effective solar eruption on 2015 November 4 from the peripheral diffusive polarities of active region 12443. We use space-borne Solar Dynamics Observatory and ACE observations. EUV images identified helical pattern along a filament channel and we regard this channel as flux-rope structure. Flow velocity derived from tracked magnetograms infers converging motion along the polarity inversion line beneath the filament channel. An associated magnetic cancellation process was detected in the converging region. Further, the pre-eruptive EUV brightening was observed in the converging region, the most intense part of which appeared in the magnetic cancellation region. These observations imply that the converging and cancelling flux probably contributed to the formation of the helical magnetic fields associated with the flux rope. A filament-height estimation method suggests that the middle part of the filament probably lies at a low altitude and was consistent with the initial place of the eruption. A thick current channel associated with the flux rope is also determined. For a expanding thick current channel, the critical height of the decay index for torus instability lies in the range of 37 - 47 Mm. Southward magnetic fields in the sheath and the ejecta induced a geomagnetic storm with a Dst global minimum of ~-90 nT.

We investigate how a weak coronal mass ejection (CME) launched on 2016 October 8 without obvious signatures in the low corona produced a relatively intense geomagnetic storm. Remote sensing observations from SDO, STEREO and SOHO and in situ measurements from WIND are employed to track the CME from the Sun to the Earth. Using a graduated cylindrical shell (GCS) model, we estimate the propagation direction and the morphology of the CME near the Sun. CME kinematics are determined from the wide-angle imaging observations of STEREO A and are used to predict the CME arrival time and speed at the Earth. We compare ENLIL MHD simulation results with in situ measurements to illustrate the background solar wind where the CME was propagating. We also apply a Grad--Shafranov technique to reconstruct the flux rope structure from in situ measurements in order to understand the geo-effectiveness associated with the CME magnetic field structure. Key results are obtained concerning how a weak CME can generate a relatively intense geomagnetic storm: (1) there were coronal holes at low latitudes, which could produce high speed streams (HSSs) to interact with the CME in interplanetary space; (2) the CME was bracketed between a slow wind ahead and a HSS behind, which enhanced the southward magnetic field inside the CME and gave rise to the unexpected geomagnetic storm.

We present an investigation of the rotation and non-radial motion of a coronal mass ejection (CME) from AR 12468 on 2015 December 16 using observations from SDO, SOHO, STEREO A and Wind. The EUV and HMI observations of the source region show that the associated magnetic flux rope (MFR) axis pointed to the east before the eruption. We use a nonlinear fore-free field (NLFFF) extrapolation to determine the configuration of the coronal magnetic field and calculate the magnetic energy density distributions at different heights. The distribution of the magnetic energy density shows a strong gradient toward the northeast. The propagation direction of the CME from a Graduated Cylindrical Shell (GCS) modeling deviates from the radial direction of the source region by about 45 degr in longitude and about 30 degr in latitude, which is consistent with the gradient of the magnetic energy distribution around the AR.The MFR axis determined by the GCS modeling points southward, which has rotated counterclockwise by about 95 degr compared with the orientation of the MFR in the low corona.The MFR reconstructed by a Grad-Shafranov (GS) method at 1 AU has almost the same orientation as the MFR from the GCS modeling, which indicates that the MFR rotation occurred in the low corona. It is the rotation of the MFR that caused the intense geomagnetic storm with the minimum Dst of -155 nT. These results suggest that the coronal magnetic field surrounding the MFR plays a crucial role in the MFR rotation and propagation direction.

The SuperCDMS experiment is designed to directly detect weakly interacting massive particles (WIMPs) that may constitute the dark matter in our Galaxy. During its operation at the Soudan Underground Laboratory, germanium detectors were run in the CDMSlite mode to gather data sets with sensitivity specifically for WIMPs with masses ${<}$10 GeV/$c^2$. In this mode, a higher detector-bias voltage is applied to amplify the phonon signals produced by drifting charges. This paper presents studies of the experimental noise and its effect on the achievable energy threshold, which is demonstrated to be as low as 56 eV$_{\text{ee}}$ (electron equivalent energy). The detector-biasing configuration is described in detail, with analysis corrections for voltage variations to the level of a few percent. Detailed studies of the electric-field geometry, and the resulting successful development of a fiducial parameter, eliminate poorly measured events, yielding an energy resolution ranging from ${\sim}$9 eV$_{\text{ee}}$ at 0 keV to 101 eV$_{\text{ee}}$ at ${\sim}$10 eV$_{\text{ee}}$. New results are derived for astrophysical uncertainties relevant to the WIMP-search limits, specifically examining how they are affected by variations in the most probable WIMP velocity and the Galactic escape velocity. These variations become more important for WIMP masses below 10 GeV/$c^2$. Finally, new limits on spin-dependent low-mass WIMP-nucleon interactions are derived, with new parameter space excluded for WIMP masses $\lesssim$3 GeV/$c^2$

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.

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.

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.

In the post-Newtonian (PN) regime, the timescale on which the spins of binary black holes precess is much shorter than the radiation-reaction timescale on which the black holes inspiral to smaller separations. On the precession timescale, the angle between the total and orbital angular momenta oscillates with nutation period $\tau$, during which the orbital angular momentum precesses about the total angular momentum by an angle $\alpha$. This defines two distinct frequencies that vary on the radiation-reaction timescale: the nutation frequency $\omega \equiv 2\pi/\tau$ and the precession frequency $\Omega \equiv \alpha/\tau$. We use analytic solutions for generic spin precession at 2PN order to derive Fourier series for the total and orbital angular momenta in which each term is a sinusoid with frequency $\Omega - n\omega$ for integer $n$. As black holes inspiral, they can pass through nutational resonances ($\Omega = n\omega$) at which the total angular momentum tilts. We derive an approximate expression for this tilt angle and show that it is usually less than $10^{-3}$ radians for nutational resonances at binary separations $r > 10M$. The large tilts occurring during transitional precession (near zero total angular momentum) are a consequence of such states being approximate $n=0$ nutational resonances. Our new Fourier series for the total and orbital angular momenta converge rapidly with $n$ providing an intuitive and computationally efficient approach to understanding generic precession that may facilitate future calculations of gravitational waveforms in the PN regime.

We discuss the synchrotron emission of fast cooling electrons in shocks. The fast cooling electrons behind the shocks can generate a position-dependent inhomogeneous electron distribution if they do not have enough time to mix homogeneously. This would lead to a very different synchrotron spectrum in low-frequency bands to that in the homogeneous case due to the synchrotron absorption. In this paper, we calculate the synchrotron spectrum in the inhomogeneous case in a gamma-ray burst (GRB). Both the forward shock and the reverse shock are considered. We find for the reverse shock dominated case, we would expect a "reverse shock bump" in the low-frequency spectrum. The spectral bump is due to the combination synchrotron absorption in both the forward and reverse shock regions. In the forward shock spectrum in the low frequencies has two unconventional segments with spectral slopes of $\lesssim1$ and $11/8$. The slope of $11/8$ has been found by some authors, while the slope of $\lesssim1$ is new, which is due to the approximately constant electron temperature in the optically thick region. In the future, simultaneous observations in multiple bands (especially in the low frequency bands) in the GRB early afterglow or prompt emission phases will possibly reveal these spectral characteristics and enable us to identify the reverse shock component and distinguish between the forward and reverse shock emissions. This also may be a method with which to diagnose the electron distribution status (homogeneous or inhomogeneous) after fast cooling in the relativistic shock region.

We investigate the coronal and interplanetary evolution of a coronal mass ejection (CME) launched on 2010 September 4 from a source region linking two active regions (ARs) 11101 and 11103, using extreme ultraviolet imaging, magnetogram, white-light and in situ observations from SDO, STEREO, SOHO, VEX and Wind. A potential-field source-surface model is employed to examine the configuration of the coronal magnetic field surrounding the source region. The graduated cylindrical shell model and a triangulation method are applied to determine the kinematics of the CME in the corona and interplanetary space. From the remote sensing and in situ observations we obtain some key results: (1) the CME was deflected in both the eastward and southward directions in the low corona by the magnetic pressure from the two ARs and possibly interacted with another ejection, which caused that the CME arrived at VEX that was longitudinally distant from the source region; (2) although VEX was closer to the Sun, the observed and derived CME arrival times at VEX are not earlier than those at Wind, which suggests the importance of determining both the frontal shape and propagation direction of the CME in interplanetary space; (3) the ICME was compressed in the radial direction while the longitudinal transverse size was extended.

Propagation of coronal mass ejections (CMEs) from the Sun far into interplanetary space is not well understood due to limited observations. In this study we examine the propagation characteristics of two geo-effective CMEs, which occurred on 2005 May 6 and 13, respectively. Significant heliospheric consequences associated with the two CMEs are observed, including interplanetary CMEs (ICMEs) at the Earth and Ulysses, interplanetary shocks, a long-duration type II radio burst, and intense geomagnetic storms. We use coronagraph observations from SOHO/LASCO, frequency drift of the long-duration type II burst, in situ measurements at the Earth and Ulysses, and magnetohydrodynamic (MHD) propagation of the observed solar wind disturbances at 1 AU to track the CMEs from the Sun far into interplanetary space. We find that both of the two CMEs underwent a major deceleration within 1 AU and thereafter a gradual deceleration when they propagated from the Earth to deep interplanetary space due to interactions with the ambient solar wind. The results also reveal that the two CMEs interacted with each other in the distant interplanetary space even though their launch times on the Sun were well separated. The intense geomagnetic storm for each case was caused by the southward magnetic fields ahead of the CME, stressing the critical role of the sheath region in geomagnetic storm generation, although for the first case there is a corotating interaction region involved.

We present a detailed study of an earth-directed coronal mass ejection (Full halo CME) event happened on 2011 February 15 making use of white light observations by three coronagraphs and radio observations by Wind/WAVES. We applied three different methods to reconstruct the propagation direction and traveling distance of the CME and its driven shock. We measured the kinematics of the CME leading edge from white light images observed by STEREO A and B, tracked the CME-driven shock using the frequency drift observed by Wind/WAVES together with an interplanetary density model, and obtained the equivalent scattering centers of the CME by Polarization Ratio(PR) method. For the first time, we applied PR method to different features distinguished from LASCO/C2 polarimetric observations and calculated their projections onto white light images observed by STEREO A and B. By combining the graduated cylindrical shell (GCS) forward modeling with the PR method, we proposed a new GCS-PR method to derive 3D parameters of a CME observed from a single perspective at Earth. Comparisons between different methods show a good degree of consistence in the derived 3D results.

SuperCDMS SNOLAB will be a next-generation experiment aimed at directly detecting low-mass (< 10 GeV/c$^2$) particles that may constitute dark matter by using cryogenic detectors of two types (HV and iZIP) and two target materials (germanium and silicon). The experiment is being designed with an initial sensitivity to nuclear recoil cross sections ~ 1 x 10$^{-43}$ cm$^2$ for a dark matter particle mass of 1 GeV/c$^2$, and with capacity to continue exploration to both smaller masses and better sensitivities. The phonon sensitivity of the HV detectors will be sufficient to detect nuclear recoils from sub-GeV dark matter. A detailed calibration of the detector response to low energy recoils will be needed to optimize running conditions of the HV detectors and to interpret their data for dark matter searches. Low-activity shielding, and the depth of SNOLAB, will reduce most backgrounds, but cosmogenically produced $^{3}$H and naturally occurring $^{32}$Si will be present in the detectors at some level. Even if these backgrounds are x10 higher than expected, the science reach of the HV detectors would be over three orders of magnitude beyond current results for a dark matter mass of 1 GeV/c$^2$. The iZIP detectors are relatively insensitive to variations in detector response and backgrounds, and will provide better sensitivity for dark matter particle masses (> 5 GeV/c$^2$). The mix of detector types (HV and iZIP), and targets (germanium and silicon), planned for the experiment, as well as flexibility in how the detectors are operated, will allow us to maximize the low-mass reach, and understand the backgrounds that the experiment will encounter. Upgrades to the experiment, perhaps with a variety of ultra-low-background cryogenic detectors, will extend dark matter sensitivity down to the "neutrino floor", where coherent scatters of solar neutrinos become a limiting background.

The prediction of the arrival time for fast coronal mass ejections (CMEs) and their associated shocks is highly desirable in space weather studies. In this paper, we use two shock propagation models, i.e. Data Guided Shock Time Of Arrival (DGSTOA) and Data Guided Shock Propagation Model (DGSPM), to predict the kinematical evolution of interplanetary shocks associated with fast CMEs. DGSTOA is based on the similarity theory of shock waves in the solar wind reference frame, and DGSPM on the non-similarity theory in the stationary reference frame. The inputs are the kinematics of the CME front at the maximum speed moment obtained from the geometric triangulation method applied to STEREO imaging observations together with the Harmonic Mean approximation. The outputs provide the subsequent propagation of the associated shock. We apply these models to the CMEs on 2012 January 19, January 23, and March 7. We find that the shock models predict reasonably well the shock's propagation after the impulsive acceleration. The shock's arrival time and local propagation speed at Earth predicted by these models are consistent with in situ measurements of WIND. We also employ the Drag-Based Model (DBM) as a comparison, and find that it predicts a steeper deceleration than the shock models after the rapid deceleration phase. The predictions of DBM at 1 AU agree with the following ICME or sheath structure, not the preceding shock. These results demonstrate the applicability of the shock models used here for future arrival time prediction of interplanetary shocks associated with fast CMEs.

The absorption feature O I 7773 is an important spectral indicator for type Ia supernovae (SNe~Ia) that can be used to trace the unburned material at outer layers of the exploding white dwarf. In this work, we use a large sample of SNe~Ia to examine this absorption at early phases (i.e., -13 days <t <-7 days), and make comparisons with the absorption features of Si~II 6355 and Ca~II near-infrared (NIR) triplet. We show that for a subgroup of spectroscopically normal SNe with normal photospheric velocities (i.e., v_si < 12,500 km s^{-1} at optical maximum), the line strength of high velocity feature (HVF) of O~I is inversely correlated with that of Si~II (or Ca~II), and this feature also shows a negative correlation with the luminosity of SNe Ia. This finding, together with other features we find for the O~I HVF, reveal that for this subgroup of SNe~Ia explosive oxygen burning occurs at the outermost layer of supernova and difference in burning there could lead to the observed diversity, which are in remarkable agreement with the popular delayed-detonation model of Chandrasekhar mass WD.

Extensive and independent observations of Type Ia supernova (SN Ia) SN 2013dy are presented, including a larger set of $UBVRI$ photometry and optical spectra from a few days before the peak brightness to $\sim$ 200 days after explosion, and ultraviolet (UV) photometry spanning from $t \approx -10$ days to $t \approx\ +15$ days referring to the $B$ band maximum. The peak brightness (i.e., $M_{\rm B} = -19.65 \pm 0.40$ mag, $L_{\rm max} = [1.95 \pm 0.55] \times 10^{43}$ erg s$^{-1}$) and the mass of synthesised $^{56}$Ni (i.e., $M$($^{56}$Ni) = 0.90 $\pm$ 0.26 M$_{\odot}$) are calculated, and they conform to the expectation for a SN Ia with a slow decline rate (i.e., $\Delta m_{15}(B)$ = 0.90 $\pm$ 0.03 mag, Phillips 1993). However, the near infrared (NIR) brightness of this SN (i.e., $M_{\rm H} = -17.33 \pm 0.30$ mag) is at least 1.0 mag fainter than usual. Besides, spectroscopy classification reveals that SN 2013dy resides on the border of "core normal" and "shallow silicon" subclasses in the Branch et al. (2009) classification scheme, or on the border of the "normal velocity" SNe Ia and 91T/99aa-like events in the Wang et al. (2009a) system. These suggest that SN 2013dy is a slow-declining SN Ia located on the transitional region of nominal spectroscopic subclasses and might not be a typical normal sample of SNe Ia.

Two extrapolation models of the solar coronal magnetic field are compared using magnetogram data from the SDO/HMI instrument. The two models, a horizontal current-current sheet-source surface (HCCSSS) model and a potential field-source surface (PFSS) model differ in their treatment of coronal currents. Each model has its own critical variable, respectively the radius of a cusp surface and a source surface, and it is found that adjusting these heights over the period studied allows better fit between the models and the solar open flux at 1 AU as calculated from the Interplanetary Magnetic Field (IMF). The HCCSSS model provides the better fit for the overall period from 2010 November to 2015 May as well as for two subsets of the period - the minimum/rising part of the solar cycle, and the recently-identified peak in the IMF from mid-2014 to mid-2015 just after solar maximum. It is found that a HCCSSS cusp surface height of 1.7 Rsun provides the best fit to the IMF for the overall period, while 1.7 & 1.9 Rsun give the best fits for the two subsets. The corresponding values for the PFSS source surface height are 2.1, 2.2 and 2.0 Rsun respectively. This means that the HCCSSS cusp surface rises as the solar cycle progresses while the PFSS source surface falls.

Black holes (BHs) hide themselves behind various astronomical phenomena, and their properties, i.e., mass and spin, are usually difficult to constrain. One leading candidate for the central engine model of gamma-ray bursts (GRBs) invokes a stellar mass BH and a neutrino-dominated accretion flow (NDAF), with the relativistic jet launched due to neutrino-anti-neutrino annihilations. Such a model gives rise to a matter-dominated fireball, and is suitable to interpret GRBs with a dominant thermal component with a photospheric origin. We propose a method to constrain BH mass and spin within the framework of this model, and apply the method to a thermally-dominant GRB 101219B whose initial jet launching radius $r_0$ is constrained from the data. Using our numerical model of NDAF jets, we estimate the following constraints on the central BH: mass $M_{\rm BH} \sim 5-9~M_\odot$, spin parameter $a_* \gtrsim 0.6$, and disk mass $3~M_\odot \lesssim M_{\rm disk} \lesssim 4~M_\odot$. Our results also suggest that the NDAF model is a competitive candidate for the central engine of GRBs with a strong thermal component.

We present well-sampled optical observations of the bright Type Ia supernova (SN~Ia) SN 2011fe in M101. Our data, starting from $\sim16$ days before maximum light and extending to $\sim463$ days after maximum, provide an unprecedented time series of spectra and photometry for a normal SN~Ia. Fitting the early-time rising light curve, we find that the luminosity evolution of SN 2011fe follows a $t^n$ law, with the index $n$ being close to 2.0 in the $VRI$ bands but slightly larger in the $U$ and $B$ bands. Combining the published ultraviolet (UV) and near-infrared (NIR) photometry, we derive the contribution of UV/NIR emission relative to the optical. SN 2011fe is found to have stronger UV emission and reaches its UV peak a few days earlier than other SNe~Ia with similar $\Delta m_{15}(B)$, suggestive of less trapping of high-energy photons in the ejecta. Moreover, the $U$-band light curve shows a notably faster decline at late phases ($t\approx 100$--300 days), which also suggests that the ejecta may be relatively transparent to UV photons. These results favor the notion that SN 2011fe might have a progenitor system with relatively lower metallicity. On the other hand, the early-phase spectra exhibit prominent high-velocity features (HVFs) of O~I $\lambda$7773 and the Ca~II~NIR triplet, but only barely detectable in Si~II~6355. This difference can be caused either by an ionization/temperature effect or an abundance enhancement scenario for the formation of HVFs; it suggests that the photospheric temperature of SN 2011fe is intrinsically low, perhaps owing to incomplete burning during the explosion of the white dwarf.

Photometric and spectroscopic observations of a slowly declining, luminous Type Ia supernova (SN Ia) SN 2011hr in the starburst galaxy NGC 2691 are presented. SN 2011hr is found to peak at $M_{B}=-19.84 \pm 0.40\,\rm{mag}$, with a post-maximum decline rate $\Delta$m$_{15}$(B) = 0.92 $\pm$ 0.03\,$\rm{mag}$. From the maximum-light bolometric luminosity, $L=(2.30 \pm 0.90) \times 10^{43}\,\rm{erg\,s^{-1}}$, we estimate the mass of synthesized \Nifs\ in SN 2011hr to be $M(\rm{^{56}Ni})=1.11 \pm 0.43\,M_{\sun}$. SN 2011hr appears more luminous than SN 1991T at around maximum light, and the absorption features from its intermediate-mass elements (IMEs) are noticeably weaker than the latter at similar phases. Spectral modeling suggests that SN 2011hr has the IMEs of $\sim$\,0.07 M$_{\sun}$ in the outer ejecta, which is much lower than the typical value of normal SNe Ia (i.e., 0.3 -- 0.4 M$_{\sun}$) and is also lower than the value of SN 1991T (i.e., $\sim$\,0.18 M$_{\sun}$). These results indicate that SN 2011hr may arise from a Chandrasekhar-mass white dwarf progenitor that experienced a more efficient burning process in the explosion. Nevertheless, it is still possible that SN 2011hr may serve as a transitional object connecting the SN 1991T-like SNe Ia with the superluminous subclass like SN 2007if given that the latter also shows very weak IMEs at all phases.

Observations of very early multi-wavelength afterglows are critical to reveal the properties of the radiating fireball and its environment as well as the central engine of gamma-ray bursts (GRBs). We report our optical observations of GRB 111228A from 95 sec to about 50 hours after the burst trigger and investigate its properties of the prompt gamma-rays and the ambient medium using our data and the data observed with {\em Swift} and {\em Fermi} missions. Our joint optical and X-ray spectral fits to the afterglow data show that the ambient medium features as low dust-to-gas ratio. Incorporating the energy injection effect, our best fit to the afterglow lightcurves with the standard afterglow model via the Markov Chain Monte Carlo (MCMC) technique shows that $\epsilon_e=(6.9\pm 0.3)\times 10^{-2}$, $\epsilon_B=(7.73\pm 0.62)\times 10^{-6}$, $E_{\rm K}=(6.32\pm 0.86)\times 10^{53}\rm erg$, $n=0.100\pm 0.014$ cm$^{-3}$. The low medium density likely implies that the afterglow jet may be in a halo or in a hot ISM. Achromatic shallow decay segment observed in the optical and X-ray bands is well explained with the long-lasting energy injection from the central engine, which would be a magnetar with a period of about 1.92 ms inferred from the data. The $E_p$ of its time-integrated prompt gamma-ray spectrum is $\sim 26$ KeV. Using the initial Lorentz factor ($\Gamma_0=476^{+225}_{-237}$) derived from our afterglow model fit, it is found that GRB 111228A satisfies the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ relation and bridges the typical GRBs and low luminosity GRBs in this relation.

The high-velocity features (HVFs) in optical spectra of type Ia supernovae (SNe Ia) are examined with a large sample including very early-time spectra (e.g., t < -7 days). Multiple Gaussian fits are applied to examine the HVFs and their evolutions, using constraints on expansion velocities for the same species (i.e., SiII 5972 and SiII 6355). We find that strong HVFs tend to appear in SNe Ia with smaller decline rates (e.g., dm15(B)<1.4 mag), clarifying that the finding by Childress et al. (2014) for the Ca-HVFs in near-maximum-light spectra applies both to the Si-HVFs and Ca-HVFs in the earlier phase. The Si-HVFs seem to be more common in fast-expanding SNe Ia, which is different from the earlier result that the Ca-HVFs are associated with SNe Ia having slower SiII 6355 velocities at maximum light (i.e., Vsi). This difference can be due to that the HVFs in fast-expanding SNe Ia usually disappear more rapidly and are easily blended with the photospheric components when approaching the maximum light. Moreover, SNe Ia with both stronger HVFs at early phases and larger Vsi are found to have noticeably redder B-V colors and occur preferentially in the inner regions of their host galaxies, while those with stronger HVFs but smaller Vsi show opposite tendencies, suggesting that these two subclasses have different explosion environments and their HVFs may have different origins. We further examine the relationships between the absorption features of SiII 6355 and CaII IR lines, and find that their photospheric components are well correlated in velocity and strength but the corresponding HVFs show larger scatter. These results cannot be explained with ionization and/or thermal processes alone, and different mechanisms are required for the creation of HVF-forming region in SNe Ia.

We present optical observations of a SN 2002cx-like supernova SN 2013en in UGC 11369, spanning from a phase near maximum light (t= +1 d) to t= +60 d with respect to the R-band maximum. Adopting a distance modulus of mu=34.11 +/- 0.15 mag and a total extinction (host galaxy+Milky Way) of $A_V \sim1.5$ mag, we found that SN 2013en peaked at $M(R)\sim -18.6$ mag, which is underluminous compared to the normal SNe Ia. The near maximum spectra show lines of Si II, Fe II, Fe III, Cr II, Ca II and other intermediate-mass and iron group elements which all have lower expansion velocities (i.e., ~ 6000 km/s). The photometric and spectroscopic evolution of SN 2013en is remarkably similar to those of SN 2002cx and SN 2005hk, suggesting that they are likely to be generated from a similar progenitor scenario or explosion mechanism.

The long-term optical and infrared color variability of blazars has been investigated with the SMARTS monitoring data. The sample in this study consists of 49 flat spectrum radio quasars and 22 BL Lac objects. The fractional variability amplitudes of each source have been calculated in both optical R band and infrared J band. Overall, the variability amplitudes of FSRQs are larger than those of BL Lac objects. The results also suggest that the variability amplitude of most FSRQs is larger at lower energy band (J band) than at higher one (R band), while the variability amplitude of BL Lacs are larger at higher energy band. Two types of blazars both display color variation along with the variability in brightness. However, they show different variation behaviors in general. With the whole data set, 35 FSRQs exhibit redder-when-brighter trends, and 11 FSRQs exhibit opposite trends; 11 BL Lacs follow bluer-when-brighter trends, and 7 BL Lacs follow opposite trends. The examination in detail shows that there are 10 blazars showing redder-when-brighter trend in their low state, and bluer-when-brighter or stable-when-brighter trends in their high state. Some more complicated color behaviors have also been detected in several blazars. The non-thermal jet emission and the thermal emission from the accretion disc are employed to explain the observed color behaviors.

We present extensive ultraviolet, optical, and near-infrared observations of the type IIP supernova (SN IIP) 2013ej in the nearby spiral galaxy M74. The multicolor light curves, spanning from $\sim$ 8--185 days after explosion, show that it has a higher peak luminosity (i.e., M$_{V}$ $\sim$$-$17.83 mag at maximum light), a faster post-peak decline, and a shorter plateau phase (i.e., $\sim$ 50 days) compared to the normal type IIP SN 1999em. The mass of $^{56}$Ni is estimated as 0.02$\pm$0.01 M$_{\odot}$ from the radioactive tail of the bolometric light curve. The spectral evolution of SN 2013ej is similar to that of SN 2004et and SN 2007od, but shows a larger expansion velocity (i.e., $v_{Fe II} \sim$ 4600 km s$^{-1}$ at t $\sim$ 50 days) and broader line profiles. In the nebular phase, the emission of H$\alpha$ line displays a double-peak structure, perhaps due to the asymmetric distribution of $^{56}$Ni produced in the explosion. With the constraints from the main observables such as bolometric light curve, expansion velocity and photospheric temperature of SN 2013ej, we performed hydrodynamical simulations of the explosion parameters, yielding the total explosion energy as $\sim$0.7$\times$ 10$^{51}$ erg, the radius of the progenitor as $\sim$600 R$_{\odot}$, and the ejected mass as $\sim$10.6 M$_{\odot}$. These results suggest that SN 2013ej likely arose from a red supergiant with a mass of 12--13 M$_{\odot}$ immediately before the explosion.

The Tsinghua University-National Astronomical Observatories of China (NAOC) Transient Survey (TNTS) is an automatic survey for a systematic exploration of optical transients (OTs), conducted with a 60/90 cm Schmidt telescope at Xinglong station of NAOC. This survey repeatedly covers ~ 1000 square degrees of the north sky with a cadence of 3-4 days. With an exposure of 60 s, the survey reaches a limited unfiltered magnitude of about 19.5 mag. This enables us to discover supernovae at their relatively young stages. In this paper, we describe the overall performance of our survey during the first year and present some preliminary results.

The IceCube (IC) collaboration recently reported the detection of TeV-PeV extraterrestrial neutrinos whose origin is yet unknown. By the photon-neutrino connection in $pp$ and $p\gamma$ interactions, we use the \fermi-LAT observations to constrain the origin of the IC detected neutrinos. We find that Galactic origins, i.e., the diffuse Galactic neutrinos due to cosmic ray (CR) propagation in the Milky Way, and the neutrinos from the Galactic point sources, may not produce the IC neutrino flux, thus these neutrinos should be of extragalactic origin. Moreover, the extragalactic gamma-ray bursts (GRBs) may not account for the IC neutrino flux, the jets of active galactic nuclei may not produce the IC neutrino spectrum, but the starburst galaxies (SBGs) may be promising sources. As suggested by the consistency between the IC detected neutrino flux and the Waxman-Bahcall bound, GRBs in SBGs may be the sources of both the ultrahigh energy, $>10^{19}$eV, CRs and the $1-100$~PeV CRs that produce the IC detected TeV-PeV neutrinos.

The water Cherenkov detector array (WCDA) for the large high altitude air shower observatory(LHAASO) will employ more than 3600 hemisphere 8 inch photomultiplier tubes (PMT). The good time performance of PMT, especially the transit time spread (TTS), is required for WCDA. TTS is usually defined as the TTS of single photoelectron, and usually determined by using single photoelectron counting technique. A method using the photoelectron spectrum is researched for the measurement of TTS. The method is appropriate for multi-photoelectrons and makes it possible to measure the TTS of different photoelectrons at the same time. The TTS of different photoelectrons is measured for Hamamatsu R5912 with the divider circuit designed in specifically. The TTS of single photoelectron is determined to 3.3 ns and satisfies the requirement of WCDA.

PKS 2155-304 is a well studied BL Lac object in the southern sky. The historical optical data during different period have been collected and compiled. Light curves with a time span of 35 years have been constructed. The R-band light curve has been analyzed by means of three methods: epoch folding method, Jurkevich method and discrete correlation function (DCF) method. It is derived that there is an evident periodic component of 317 days (i.e. 0.87 yr) superposed on a long-term trend with large-amplitude variation in the light curve. The variability of this source is accompanied by a slight color variation, and the brightness and color index are correlated with each other. On the long time-scale, PKS 2155-304 exhibits a tendency of bluer-when-brighter, which means the spectrum becomes flatter when the source brightens.

The Milky Way is the product of a complex evolution of generations of mergers, collapse, star formation, supernova and collisional heating, radiative and collisional cooling, and ejected nucleosynthesis. Moreover, all of this occurs in the context of the cosmic expansion, the formation of cosmic filaments, dark-matter halos, spiral density waves, and emerging dark energy. In this review we summarize observational evidence and discuss recent calculations concerning the formation, evolution nucleosynthesis in the galaxies of the Local-Group. In particular, we will briefly summarize observations and simulations for the dwarf galaxies and the two large spirals of the Local Group. We discuss how galactic halos form within the dark matter filaments that define a super-galactic plane. Gravitational interaction along this structure leads to streaming flows toward the two dominant galaxies in the cluster. These simulations and observations also suggest that a significant fraction of the Galactic halo formed as at large distances and then arrived later along these streaming flows. We also consider the insight provided by observations and simulations of nucleosynthesis both within the galactic halo and dwarf galaxies in the Local Group.

Extensive optical and ultraviolet (UV) observations of the type Ia supernova (SN Ia) 2012fr are presented in this paper. It has a relatively high luminosity, with an absolute $B$-band peak magnitude of about $-19.5$ mag and a smaller post-maximum decline rate than normal SNe Ia [e.g., $\Delta m _{15}$($B$) $= 0.85 \pm 0.05$ mag]. Based on the UV and optical light curves, we derived that a $^{56}$Ni mass of about 0.88 solar masses was synthesized in the explosion. The earlier spectra are characterized by noticeable high-velocity features of \ion{Si}{2} $\lambda$6355 and \ion{Ca}{2} with velocities in the range of $\sim22,000$--$25,000$ km s$^{-1}$. At around the maximum light, these spectral features are dominated by the photospheric components which are noticeably narrower than normal SNe Ia. The post-maximum velocity of the photosphere remains almost constant at $\sim$12,000 km s$^{-1}$ for about one month, reminiscent of the behavior of some luminous SNe Ia like SN 1991T. We propose that SN 2012fr may represent a subset of the SN 1991T-like SNe Ia viewed in a direction with a clumpy or shell-like structure of ejecta, in terms of a significant level of polarization reported in Maund et al. (2013).

We present multiwavelength observations of the afterglow of GRB 130427A, the brightest (in total fluence) gamma-ray burst of the past 29 years. Optical spectroscopy from Gemini-North reveals the redshift of the GRB to be z=0.340, indicating that its unprecedented brightness is primarily the result of its relatively close proximity to Earth; the intrinsic luminosities of both the GRB and its afterglow are not extreme in comparison to other bright GRBs. We present a large suite of multiwavelength observations spanning from 300 s to 130 d after the burst and demonstrate that the afterglow shows relatively simple, smooth evolution at all frequencies with no significant late-time flaring or rebrightening activity. The entire dataset from 1 GHz to 10 GeV can be modeled as synchrotron emission from a combination of reverse and forward shocks in good agreement with the standard afterglow model, providing strong support to the applicability of the underlying theory and clarifying the nature of the GeV emission observed to last for minutes to hours following other very bright GRBs. A tenuous, wind-stratified circumburst density profile is required by the observations, suggesting a massive-star progenitor with a low mass-loss rate, perhaps due to low metallicity. GRBs similar in nature to GRB 130427A, inhabiting low-density media and exhibiting strong reverse shocks, are probably not uncommon but may have been difficult to recognize in the past due to their relatively faint late-time radio emission; more such events should be found in abundance by the new generation of sensitive radio and millimeter instruments.

In the internal shock model for gamma-ray bursts (GRBs), the synchrotron spectrum from the fast cooling electrons in a homogeneous downstream magnetic field (MF) is too soft to produce the low-energy slope of GRB spectra. However the magnetic field may decay downstream with distance from the shock front. Here we show that the synchrotron spectrum becomes harder if electrons undergo synchrotron and inverse-Compton cooling in a decaying MF. To reconcile this with the typical GRB spectrum with low energy slope $\nu F_\nu\propto\nu$, it is required that the postshock MF decay time is comparable to the cooling time of the bulk electrons (corresponding to a MF decaying length typically of $\sim10^5$ skin depths); that the inverse-Compton cooling should dominate synchrotron cooling after the MF decay time; and/or that the MF decays with comoving time roughly as $B\propto t^{-1.5}$. An internal shock synchrotron model with a decaying MF can account for the majority of GRBs with low energy slopes not harder than $\nu^{4/3}$.

Thanks to the rapid progress in precision cosmology in the last few years, we now have access to physical observables that may constrain the theory of inflation through the non-Gaussianity (NG) signatures in the cosmic microwave background radiation and the distribution of large-scale structure. Numerical modeling of the NG signals from different inflation models is essential to correctly interpret current and near future data from large-scale structure surveys. In this study, we use high-resolution cosmological hydrodynamical simulations to investigate the effects of primordial NG on the formation and evolution of galaxies from the cosmic dawn to the present day. Focusing on the local type primordial NG, we find that it may affect the formation history of stars and black holes in galaxies, and their distribution. Compared to the Gaussian case, large non-Gaussian potential with $f_{NL} \gtrsim 10^3$ leads to earlier collapse of the first structures, more massive galaxies especially at high redshifts, stronger clustering of galaxies, and higher halo bias. However, for smaller NG with $f_{NL} \lesssim 10^2$, the effect is significantly weaker. Observations of the distribution and properties of high-redshift, rare objects such as the first galaxies and quasars may provide further constraints on the primordial NG.

Type Ia supernovae (SNe Ia) have been used as excellent standardizable candles for measuring cosmic expansion, but their progenitors are still elusive. Here we report that the spectral diversity of SNe Ia is tied to their birthplace environments. We find that those with high-velocity ejecta are substantially more concentrated in the inner and brighter regions of their host galaxies than are normal-velocity SNe Ia. Furthermore, the former tend to inhabit larger and more-luminous hosts. These results suggest that high-velocity SNe Ia likely originate from relatively younger and more metal-rich progenitors than normal-velocity SNe Ia, and are restricted to galaxies with substantial chemical evolution.

We present our successful Chandra program designed to identify, with sub-arcsecond accuracy, the X-ray afterglow of the short GRB 111117A, which was discovered by Swift and Fermi. Thanks to our rapid target of opportunity request, Chandra clearly detected the X-ray afterglow, though no optical afterglow was found in deep optical observations. The host galaxy was clearly detected in the optical and near-infrared band, with the best photometric redshift of z=1.31_{-0.23}^{+0.46} (90% confidence), making it one of the highest known short GRB redshifts. Furthermore, we see an offset of 1.0 +- 0.2 arcseconds, which corresponds to 8.4 +- 1.7 kpc, between the host and the afterglow position. We discuss the importance of using Chandra for obtaining sub-arcsecond X-ray localizations of short GRB afterglows to study GRB environments.

According to the second law of thermodynamics, the arrow of time points to an ever increasing entropy of the Universe. However, exactly how the entropy evolves with time and what drives the growth remain largely unknown. Here, for the first time, we quantify the evolving entropy of cosmic structures using a large-scale cosmological hydrodynamical simulation. Our simulation starts from initial conditions predicted by the leading LambdaCDM cosmology, self-consistently evolves the dynamics of both dark and baryonic matter, star formation, black hole growth and feedback processes, from the cosmic dawn to the present day. Tracing the entropy contributions of these distinct components in the simulation, we find a strong link between entropy growth and structure formation. The entropy is dominated by that of the black holes in all epochs, and its evolution follows the same path as that of galaxies: it increases rapidly from a low-entropy state at high redshift until z~2, then transits to a slower growth. Our results suggest that cosmic entropy may co-evolve with cosmic structure, and that its growth may be driven mainly by the formation of black holes in galaxies. We predict that the entropy will continue to increase in the near future, but likely at a constant rate.

We present extensive optical observations of a Type IIn supernova (SN) 2010jl for the first 1.5 years after the discovery. The UBVRI light curves demonstrated an interesting two-stage evolution during the nebular phase, which almost flatten out after about 90 days from the optical maximum. SN 2010jl has one of the highest intrinsic H_alpha luminosity ever recorded for a SN IIn, especially at late phase, suggesting a strong interaction of SN ejecta with the dense circumstellar material (CSM) ejected by the progenitor. This is also indicated by the remarkably strong Balmer lines persisting in the optical spectra. One interesting spectral evolution about SN 2010jl is the appearance of asymmetry of the Balmer lines. These lines can be well decomposed into a narrow component and an intermediate-width component. The intermediate-width component showed a steady increase in both strength and blueshift with time until t ~ 400 days after maximum, but it became less blueshifted at t ~ 500 days when the line profile appeared relatively symmetric again. Owing to that a pure reddening effect will lead to a sudden decline of the light curves and a progressive blueshift of the spectral lines, we therefore propose that the asymmetric profiles of H lines seen in SN 2010jl is unlikely due to the extinction by newly formed dust inside the ejecta, contrary to the explanation by some early studies. Based on a simple CSM-interaction model, we speculate that the progenitor of SN 2010jl may suffer a gigantic mass loss (~ 30-50 M_sun) in a few decades before explosion. Considering a slow moving stellar wind (e.g., ~ 28 km/s) inferred for the preexisting, dense CSM shell and the extremely high mass-loss rate (1-2 M_sun per yr), we suggest that the progenitor of SN 2010jl might have experienced a red supergiant stage and explode finally as a post-red supergiant star with an initial mass above 30-40 M_sun.

Previous studies have found that the width of gamma-ray burst (GRB) pulse is energy dependent and that it decreases as a power-law function with increasing photon energy. In this work we have investigated the relation between the energy dependence of pulse and the so-called Band spectrum by using a sample including 51 well-separated fast rise and exponential decay long-duration GRB pulses observed by BATSE (Burst and Transient Source Experiment on the Compton Gamma Ray Observatory). We first decompose these pulses into rise, and decay phases and find the rise widths, and the decay widths also behavior as a power-law function with photon energy. Then we investigate statistically the relations between the three power-law indices of the rise, decay and total width of pulse (denoted as $\delta_r$, $\delta_d$ and $\delta_w$, respectively) and the three Band spectral parameters, high-energy index ($\alpha$), low-energy index ($\beta$) and peak energy ($E_p$). It is found that (1)$\alpha$ is strongly correlated with $\delta_w$ and $\delta_d$ but seems uncorrelated with $\delta_r$; (2)$\beta$ is weakly correlated with the three power-law indices and (3)$E_p$ does not show evident correlations with the three power-law indices. We further investigate the origin of $\delta_d-\alpha$ and $\delta_w-\alpha$. We show that the curvature effect and the intrinsic Band spectrum could naturally lead to the energy dependence of GRB pulse width and also the $\delta_d-\alpha$ and $\delta_w-\alpha$ correlations. Our results would hold so long as the shell emitting gamma rays has a curve surface and the intrinsic spectrum is a Band spectrum or broken power law. The strong $\delta_d-\alpha$ correlation and inapparent correlations between $\delta_r$ and three Band spectral parameters also suggest that the rise and decay phases of GRB pulses have different origins.

In this paper, ten CME events viewed by the STEREO twin spacecraft are analyzed to study the deflections of CMEs during their propagation in the corona. Based on the three-dimensional information of the CMEs derived by the graduated cylindrical shell (GCS) model [Thernisien et al., 2006], it is found that the propagation directions of eight CMEs had changed. By applying the theoretical method proposed by Shen et al. [2011] to all the CMEs, we found that the deflections are consistent, in strength and direction, with the gradient of the magnetic energy density. There is a positive correlation between the deflection rate and the strength of the magnetic energy density gradient and a weak anti-correlation between the deflection rate and the CME speed. Our results suggest that the deflections of CMEs are mainly controlled by the background magnetic field and can be quantitatively described by the magnetic energy density gradient (MEDG) model.

Long-lived >100 MeV emission has been a common feature of most Fermi-LAT detected gamma-ray bursts (GRBs), e.g., detected up to ~10^3s in long GRBs 080916C and 090902B and ~10^2s in short GRB 090510. This emission is consistent with being produced by synchrotron emission of electrons accelerated to high energy by the relativistic collisionless shock propagating into the weakly magnetized medium. Here we show that this high-energy afterglow emission constrains the preshock magnetic field to satisfy 1(n/1cc)^{9/8} mG<B<10^2(n/1cc)^{3/8}mG, where n is the preshock density, more stringent than the previous constraint by X-ray afterglow observations on day scale. This suggests that the preshock magnetic field is strongly amplified, most likely by the streaming of high energy shock accelerated particles.

General relativistic corrections to the expansion rate of the Universe arise when the Einstein equations are averaged over a spatial volume in a locally inhomogeneous cosmology. It has been suggested that they may contribute to the observed cosmic acceleration. In this paper, we propose a new scheme that utilizes numerical simulations to make a realistic estimate of the magnitude of these corrections for general inhomogeneities in (3+1) spacetime. We then quantitatively calculate the volume averaged expansion rate using N-body large-scale structure simulations and compare it with the expansion rate in a standard FRW cosmology. We find that in the weak gravitational field limit, the converged corrections are slightly larger than the previous claimed 10^{-5} level, but not large enough nor even of the correct sign to drive the current cosmic acceleration. Nevertheless, the question of whether the cumulative effect can significantly change the expansion history of the Universe needs to be further investigated with strong-field relativity.

The Lorentz factor (LF) of gamma-ray burst (GRB) ejecta may be constrained by observations of high-energy (HE) spectral attenuation. The recent Fermi-LAT observations of prompt GeV emission from several bright GRBs have leaded to conclusions of unexpectedly large LFs, $\Gamma>10^3$. Here we revisit this problem with two main concerns. (1) With one-zone assumption where all photons are assumed to be generated in the same region (radius) and time, we {\em self-consistently} calculate the $\gamma\gamma$ optical depth by adopting a target photon spectrum with HE cutoff. We find that this might be important when the GRB LF is below a few hundreds. (2) Recent Fermi-LAT observations suggest that the bulk MeV-range and HE ($\ga100$ MeV) emission may arise from different regions. We then consider a two-zone case where HE emission is generated in much larger radii than that of the MeV-range emission. We find that the HE emission may be mainly attenuated by MeV-range emission and that the attenuated HE spectrum does not show an exponential spectral cutoff but a slight steepening. This suggests that there may be no abrupt cutoff due to $\gamma\gamma$ attenuation if relaxing the one-zone assumption. By studying the spectra of three bright Fermi-LAT GRBs 080916C, 090510 and 090902B, we show that a bulk LF of\textbf{ $\Gamma\sim600$ }can be consistent with observations in the two-zone case. Even lower LFs can be obtained in the multi-zone case.

In this paper we have analyzed the temporal and spectral behavior of 52 Fast Rise and Exponential Decay (FRED) pulses in 48 long-duration gamma-ray bursts (GRBs) observed by the CGRO/BATSE, using a pulse model with two shape parameters and the Band model with three shape parameters, respectively. It is found that these FRED pulses are distinguished both temporally and spectrally from those in long-lag pulses. Different from these long-lag pulses only one parameter pair indicates an evident correlation among the five parameters, which suggests that at least $\sim$4 parameters are needed to model burst temporal and spectral behavior. In addition, our studies reveal that these FRED pulses have correlated properties: (i) long-duration pulses have harder spectra and are less luminous than short-duration pulses; (ii) the more asymmetric the pulses are the steeper the evolutionary curves of the peak energy ($E_{p}$) in the $\nu f_{\nu}$ spectrum within pulse decay phase are. Our statistical results give some constrains on the current GRB models.

The issue of the influence of coronal holes (CHs) on coronal mass ejections (CMEs) in causing solar energetic particle (SEP) events is revisited. It is a continuation and extension of our previous work (Shen et al., 2006), in which no evident effect of CHs on CMEs in generating SEPs were found by statistically investigating 56 CME events. This result is consistent with the conclusion obtained by Kahler in 2004. In this paper, we extrapolate the coronal magnetic field, define CHs as the regions consisting of only open magnetic field lines and perform a similar analysis on this issue for totally 76 events by extending the study interval to the end of 2008. Three key parameters, CH proximity, CH area and CH relative position, are involved in the analysis. The new result confirms the previous conclusion that CHs did not show any evident effect on CMEs in causing SEP events.

A model of jet precession driven by a neutrino-cooled disc around a spinning black hole is present in order to explain the temporal structure and spectral evolution of gamma-ray bursts (GRBs). The differential rotation of the outer part of a neutrino dominated accretion disc may result in precession of the inner part of the disc and the central black hole, hence drives a precessed jet via neutrino annihilation around the inner part of the disc. Both analytic and numeric results for our model are present. Our calculations show that a black hole-accretion disk system with black hole mass $M \simeq 3.66 M_\odot$, accretion rate $\dot{M} \simeq 0.54 M_\odot \rm s^{-1}$, spin parameter $a=0.9$ and viscosity parameter $\alpha=0.01$ may drive a precessed jet with period P=1 s and luminosity $L=10^{51}$ erg s$^{-1}$, corresponding to the scenario for long GRBs. A precessed jet with $P=0.1$s and $L=10^{50}$ erg s$^{-1}$ may be powered by a system with $M \simeq 5.59 M_\odot$, $\dot{M} \simeq 0.74 M_\odot \rm s^{-1}$, $a=0.1$, and $\alpha=0.01$, possibly being responsible for the short GRBs. Both the temporal and spectral evolution in GRB pulse may explained with our model. GRB central engines likely power a precessed jet driven by a neutrino-cooled disc. The global GRB lightcurves thus could be modulated by the jet precession during the accretion timescale of the GRB central engine. Both the temporal and spectral evolution in GRB pulse may be due to an viewing effect due to the jet precession.

Employing two samples containing of 56 and 59 well-separated FRED (fast rise and exponential decay) gamma-ray burst (GRB) pulses whose spectra are fitted by the Band spectrum and Compton model, respectively, we have investigated the evolutionary slope of $E_{p}$ (where $E_{p}$ is the peak energy in the $\nu F\nu$ spectrum) with time during the pulse decay phase. The bursts in the samples were observed by the Burst and Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory. We first test the $E_{p}$ evolutionary slope during the pulse decay phase predicted by Lu et al. (2007) based on the model of highly symmetric expanding fireballs in which the curvature effect of the expanding fireball surface is the key factor concerned. It is found that the evolutionary slopes are normally distributed for both samples and concentrated around the values of 0.73 and 0.76 for Band and Compton model, respectively, which is in good agreement with the theoretical expectation of Lu et al. (2007). However, the inconsistence with their results is that the intrinsic spectra of most of bursts may bear the Comptonized or thermal synchrotron spectrum, rather than the Band spectrum. The relationships between the evolutionary slope and the spectral parameters are also checked. We show the slope is correlated with $E_{p}$ of time-integrated spectra as well as the photon flux but anticorrelated with the lower energy index $\alpha$. In addition, a correlation between the slope and the intrinsic $E_{p}$ derived by using the pseudo-redshift is also identified. The mechanisms of these correlations are unclear currently and the theoretical interpretations are required.

Coronal shocks are important structures, but there are no direct observations of them in solar and space physics. The strength of shocks plays a key role in shock-related phenomena, such as radio bursts and solar energetic particle (SEP) generation. This paper presents an improved method of calculating Alfven speed and shock strength near the Sun. This method is based on using as many observations as possible, rather than one-dimensional global models. Two events, a relatively slow CME on 2001 September 15 and a very fast CME on 2000 June 15, are selected to illustrate the calculation process. The calculation results suggest that the slow CME drove a strong shock, with Mach number of 3.43 - 4.18, while the fast CME drove a relatively weak shock, with Mach number of 1.90 - 3.21. This is consistent with the radio observations, which find a stronger and longer decameter-hectometric (DH) type II radio burst during the first event, and a short DH type II radio burst during the second event. In particular, the alculation results explain the observational fact that the slow CME produced a major solar energetic particle (SEP) event, while the fast CME did not. Through a comparison of the two events, the importance of shock strength in predicting SEP events is addressed.

In this paper, standard accretion disk models of AGNs are tested using light curves of 26 objects well observed for reverberation mapping. Time scales of variations are estimated by the most common definition of the variability time scale and the zero-crossing time of the autocorrelation function of the optical light curves for each source. The measured time scales of variations by the two methods are consistent with each other. If the typical value of the viscosity parameter $\alpha \sim 0.1$ is adopted, the measured optical variability time scales are most close to the thermal time scales of the standard disks. If $\alpha$ is allowed to range from $\sim 0.03$ to $\sim 0.2$, the measured time scales are consistent with the thermal time scales of the standard disks. There is a linear relation between the measured variability time scales and black hole masses, and this linear relation is qualitatively consistent with expectation of the standard accretion disk models. The time lags measured by the ZDCF between different bands are on the order of days. The measured time lags of NGC 4151 and NGC 7469 are marginally consistent with the time lags estimated in the case of continuum thermal reprocessing for the standard accretion disk models. However, the measured time lags of NGC 5548 and Fairall 9 are unlikely to be the case of continuum thermal reprocessing. Our results are unlikely to be inconsistent with or are likely to be conditionally in favor of the standard accretion disk models of AGNs.

The high efficiency of converting kinetic energy into gamma-rays estimated with late-time afterglows in Gamma-Ray Burst (GRB) phenomenon challenges the commonly accepted internal-shock model. However, the efficiency is still highly uncertain because it is sensitive to many effects. In this Letter we study the sideways expansion effect of jets on estimating the efficiency. We find that this effect is considerable, reducing the efficiency by a factor of $\sim0.5$ for typical parameters, when the afterglow data $\sim 10$ hr after the GRB trigger are used to derive the kinetic energy. For a more dense circumburst medium, this effect is more significant. As samples, taking this effect into account, we specifically calculate the efficiency of two bursts whose parameters were well constrained. Almost the same results are derived. This suggests that the sideways expansion effect should be considered when the GRB efficiency is estimated with the late afterglow data.

White-light observations of the solar corona show that there are two characteristic types of Coronal Mass Ejections (CMEs) in terms of speed-height profiles: so-called fast CMEs that attain high speeds low in the corona and slow CMEs that gradually accelerate from low initial speeds. Low and Zhang (2002) have recently proposed that fast and slow CMEs result from initial states with magnetic configurations characterized by normal prominences (NPs) and inverse prominences (IPs), respectively. To test their theory, we employed a two-dimensional, time-dependent, resistive magnetohydrodynamic code to simulate the expulsion of CMEs in these two different prominence environments. Our numerical simulations demonstrate that (i) a CME-like expulsion is more readily produced in an NP than in an IP environment, and, (ii) a CME originating from an NP environment tends to have a higher speed early in the event than one originating from an IP environment. Magnetic reconnection plays distinct roles in the two different field topologies of these two environments to produce their characteristic CME speed-height profiles. Our numerical simulations support the proposal of Low and Zhang (2002) although the reconnection development for the NP associated CME is different from the one sketched in their theory. Observational implications of our simulations are discussed.

We present deep optical photometry of the afterglow of gamma-ray burst (GRB) 041006 and its associated hypernova obtained over 65 days after detection (55 R-band epochs on 10 different nights). Our early data (t<4 days) joined with published GCN data indicates a steepening decay, approaching F_nu ~t^{-0.6} at early times (<<1 day) and F_nu ~t^{-1.3} at late times. The break at t_b=0.16+-0.04 days is the earliest reported jet break among all GRB afterglows. During our first night, we obtained 39 exposures spanning 2.15 hours from 0.62 to 0.71 days after the burst that reveal a smooth afterglow, with an rms deviation of 0.024 mag from the local power-law fit, consistent with photometric errors. After t~4 days, the decay slows considerably, and the light curve remains approximately flat at R~24 mag for a month before decaying by another magnitude to reach R~25 mag two months after the burst. This ``bump'' is well-fitted by a k-corrected light curve of SN1998bw, but only if stretched by a factor of 1.38 in time. In comparison with the other GRB-related SNe bumps, GRB 041006 stakes out new parameter space for GRB/SNe, with a very bright and significantly stretched late-time SN light curve. Within a small sample of fairly well observed GRB/SN bumps, we see a hint of a possible correlation between their peak luminosity and their ``stretch factor'', broadly similar to the well-studied Phillips relation for the type Ia supernovae.