The interstellar oxygen isotopic ratio of 18O/17O can reflect the relative amount of the secular enrichment by ejecta from high-mass versus intermediate-mass stars. Previous observations found a Galactic gradient of 18O/17O, i.e., low ratios in the Galactic center and large values in the Galactic disk, which supports the insideout formation scenario of our Galaxy. However, the observed objects are not many and, in particular, not so many at large galactocentric distances. Thus we started a systematic study on Galactic interstellar 18O/17O, through observations of C18O and C17O multi-transition lines toward a large sample of 286 sources (at least one order of magnitude larger than previous ones), from the Galactic center region to the far outer Galaxy (~22 kpc). Here we present our observations of J=1-0 lines of C18O and C17O, with the ARO12m and the IRAM 30m telescope. We detected successfully both C18O and C17O 1-0 lines for 34 sources among our IRAM30m sample of 50 targets and for 166 sources among our ARO12m sample of 260 targets. The C18O optical depth effect on our ratio results, evaluated by fitting results of C17O spectra with hyperfine components and our RADEX non-LTE model calculation for the strongest source, was found to be insignificant. Beam dilution does not seem to be a problem either, which was supported by the fact of no systematic variation between the isotopic ratio and the heliocentric distance, and consistent measured ratios from two telescopes for most of those detected sources. Our results, though there are still very few detections made for sources in the outer Galaxy, confirm the apparent 18O/17O gradient of 18O/17O = (0.10+-0.03)R_GC+(2.95+-0.30), with a Pearson's rank correlation coefficient R = 0.69. This is supported by the newest Galactic chemical evolution model including the impact of massive stellar rotators and novae.

The gravity field maps of the satellite gravimetry missions GRACE (Gravity Recovery and Climate Experiment) and GRACE Follow-On are derived by means of precise orbit determination. The key observation is the biased inter-satellite range, which is measured primarily by a K-Band Ranging system (KBR) in GRACE and GRACE Follow-On. The GRACE Follow-On satellites are additionally equipped with a Laser Ranging Interferometer (LRI), which provides measurements with lower noise compared to the KBR. The biased range of KBR and LRI needs to be converted for gravity field recovery into an instantaneous range, i.e. the biased Euclidean distance between the satellites' center-of-mass at the same time. One contributor to the difference between measured and instantaneous range arises due to the non-zero travel time of electro-magnetic waves between the spacecraft. We revisit the calculation of the light time correction (LTC) from first principles considering general relativistic effects and state-of-the-art models of Earth's potential field. The novel analytical expressions for the LTC of KBR and LRI can circumvent numerical limitations of the classical approach. The dependency of the LTC on geopotential models and on the parameterization is studied, and afterwards the results are compared against the LTC provided in the official datasets of GRACE and GRACE Follow-On. It is shown that the new approach has a significantly lower noise, well below the instrument noise of current instruments, especially relevant for the LRI, and even if used with kinematic orbit products. This allows calculating the LTC accurate enough even for the next generation of gravimetric missions.

In this article, we present the application of the weighted horizontal gradient of magnetic field (WGM) flare prediction method to 3-dimensional (3D) extrapolated magnetic configurations of 13 flaring solar active regions (ARs). The main aim is to identify an optimal height range, if any, in the interface region between the photosphere and lower corona, where the flare onset time prediction capability of WGM is best exploited. The optimal height is where flare prediction, by means of the WGM method, is achieved earlier than at the photospheric level. 3D magnetic structures, based on potential and non-linear force-free field extrapolations, are constructed to study a vertical range from the photosphere up to the low corona with a 45 km step size. The WGM method is applied as a function of height to all 13 flaring AR cases that are subject to certain selection criteria. We found that applying the WGM method between 1000 and 1800 km above the solar surface would improve the prediction of the flare onset time by around 2-8 hours.Certain caveats and an outlook for future work along these lines are also discussed.

Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $\sigma_{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm{cm}^2$ (6 $ \times 10^{-38} \mathrm{cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1~TeV annihilating into $\tau^+ \tau^-$ ($b\bar{b}$) with a local density of $\rho_{\mathrm{DM}} = 0.3~\mathrm{ GeV/cm}^3$. The constraints scale inversely with $\rho_{\mathrm{DM}}$ and are independent of the DM velocity distribution.

The structural property of the magnetic field in flare-bearing solar active regions (ARs) is one of the key aspects for understanding and forecasting solar flares. In this paper, we make a comparative analysis on the chirality and magnetic configurations associated with two X-class two-ribbon flares happening in AR 10930 and AR 11158. The photospheric magnetic fields of the two ARs were observed by space-based instruments, and the corresponding coronal magnetic fields were calculated based on the nonlinear force-free field model. The analysis shows that the electric current in the two ARs was distributed mostly around the main polarity inversion lines (PILs) where the flares happened, and the magnetic chirality (indicated by the signs of force-free factor $\alpha$) along the main PILs is opposite for the two ARs, i.e., left-handed ($\alpha<0$) for AR 10930 and right-handed ($\alpha>0$) for AR 11158. It is found that, for both the flare events, a prominent magnetic connectivity (featured by co-localized strong $\alpha$ and strong current density distributions) was formed along the main PIL before flare and was totally broken after flare eruption. The two branches of the broken magnetic connectivity, combined with the prominent magnetic connectivity before flare, compose the opposite magnetic configurations in the two ARs owing to their opposite chirality, i.e., Z-shaped configuration in AR 10930 with left-handed chirality and inverse Z-shaped configuration in AR 11158 with right-handed chirality. It is speculated that two-ribbon flares can be generally classified to these two magnetic configurations by chirality in the flare source regions of ARs.

With Aperture synthesis (AS) technique, a number of small antennas can assemble to form a large telescope which spatial resolution is determined by the distance of two farthest antennas instead of the diameter of a single-dish antenna. Different from direct imaging system, an AS telescope captures the Fourier coefficients of a spatial object, and then implement inverse Fourier transform to reconstruct the spatial image. Due to the limited number of antennas, the Fourier coefficients are extremely sparse in practice, resulting in a very blurry image. To remove/reduce blur, "CLEAN" deconvolution was widely used in the literature. However, it was initially designed for point source. For extended source, like the sun, its efficiency is unsatisfied. In this study, a deep neural network, referring to Generative Adversarial Network (GAN), is proposed for solar image deconvolution. The experimental results demonstrate that the proposed model is markedly better than traditional CLEAN on solar images.

It is well-known that there is a gradient, there will drive a flow inevitably. For example, a density-gradient may drive a diffusion flow, an electrical potential-gradient may drive an electric current in plasmas, etc. Then, what will be driven when a magnetic-gradient occurs in solar atmospheric plasmas? Considering the ubiquitous distribution of magnetic-gradient in solar plasma loops, this work demonstrates that magnetic-gradient pumping (MGP) mechanism is valid even in the partial ionized solar photosphere, chromosphere as well as in the corona. It drives energetic particle flows which carry and convey kinetic energy from the underlying atmosphere to move upwards, accumulate around the looptop and increase there temperature and pressure, and finally lead to eruptions around the looptop by triggering ballooning instabilities. This mechanism may explain the formation of the observing hot cusp-structures above flaring loops in most preflare phases, therefore, the magnetic-gradient should be a natural driver of solar eruptions. Furthermore, we may also apply to understand many other astrophysical phenomena, such as the temperature distribution above sunspots, the formation of solar plasma jets, type-II spicule, and fast solar wind above coronal holes, as well as the fast plasma jets related to white dwarfs, neutron stars and black holes.

We present a new method of DEEM, the direct energy encircling method, for characterising the performance of fibres in most astronomical spectroscopic applications. It's a versatile platform to measure focal ratio degradation (FRD), throughput, and point spread function (PSF). The principle of DEEM and the relation between the encircled energy (EE) and the spot size were derived and simulated based on the power distribution model (PDM). We analysed the errors of DEEM and pointed out the major error source for better understanding and optimisation. The validation of DEEM has been confirmed by comparing the results with conventional method which shows that DEEM has good robustness with high accuracy in both stable and complex experiment environments. Applications on the integral field unit (IFU) show that the FRD of 50$\mu$m core fibre is substandard for the requirement which requires the output focal ratio to be slower than 4.5. The homogeneity of throughput is acceptable and higher than 85 per cent. The prototype IFU of the first generation helps to find out the imperfections to optimise the new design of the next generation based on the staggered structure with 35$\mu$m core fibres of $N.A.$=0.12, which can improve the FRD performance. The FRD dependence on wavelength and core size is revealed that higher output focal ratio occurs at shorter wavelengths for large core fibres, which is in agreement with the prediction of PDM. But the dependence of the observed data is weaker than the prediction.

Focal ratio degradation (FRD) is a major contributor to throughput and light loss in a fibre spectroscopic telescope system. We combine the guided mode theory in geometric optics and a well-known model, power distribution model (PDM), to predict and explain the FRD dependence properties. We present a robust method by modifying the energy distribution method (EDM) with \emph{f-intercept} to control the input condition. This method provides a way to determine the proper position of the fibre end on the focal plane to improve energy utilization and FRD performance, which lifts the relative throughput up to 95\% with variation of output focal ratio less than 2\%. And this method can also help to optimize the arrangement of the position of focal-plane plate to enhance the coupling efficiency in a telescope. To investigate length properties, we modified PDM by introducing a new parameter, focal distance \emph{f}, into the original model to make it available for multi-position measurement system. The results show that the modified model is robust and feasible for measuring the key parameter \emph{d}$_0$ to simulate the transmission characteristics. The output focal ratio in the experiment does not follow the prediction trend but shows an interesting phenomenon that the output focal ratio increases at first to the peak, then decreases and remains stable finally with increasing fibre length longer than 15m, which provides a reference for choosing appropriate length of fibre to improve the FRD performance for the design of the fibre system in a telescope.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$.

Base on about 4,500 large tangential velocity ($V_\mathrm{tan}>0.75V_\mathrm{esc}$) with high-precision proper motions and $5\sigma$ parallaxes in Gaia DR2 5D information derived from parallax and proper motion, we identify more than 600 high velocity stars with $50\%$ unbound probability. Of these, 28 nearby (less than 6 kpc) late-type Hypervelocity stars (HVSs) with over $99\%$ possibility of unbound are discovered. In order to search for the unbound stars from the full Gaia DR2 6D phase space information derived from parallax, proper motion and radial velocity, we also identify 28 stars from the total velocity ($V_\mathrm{gc}>0.75V_\mathrm{esc}$) that have probabilities greater than $50\%$ of being unbound from the Galaxy. Of these, only three have a nearly $99\%$ probabilities of being unbound. On the whole HVSs subsample, there is 12 sources reported by other surveys. We study the spatial distribution of angular positions and angular separation of HVSs. We find the unbound HVSs are spatially anisotropic that is most significant in the Galactic longitude at more than $3\sigma$ level, and lower unbound probability HVSs are systematically more isotropic. The spatial distribution can reflect the origin of HVSs and we discuss the possible origin link with the anisotropy.

The Ultra-Long Wavelength (ULW) regime of longer than 10 m (corresponding to frequencies below 30 MHz), remains as the last virtually unexplored window in radio astronomy, and is presently attracting considerable attention as an area of potentially rewarding studies. However, the opaqueness of the Earth's ionosphere makes the ULW celestial radio emission very difficult to detect with ground-based instrumentation. The impact of the ionosphere on ULW radio emission depends on the Solar cycle activities and varies with time. In addition, the ULW spectrum region is densely populated by intensive artificial radio frequency interference (RFI). An obvious solution of these problems is to place an ULW radio telescope in space. However, this solution is expensive and poses non-negligible technological challenges. An alternative approach is triggered by recent studies showing that the period of post 2020 will be most suitable for exploratory ground-based ULW radio observations due to the expected 'calm' state of the ionosphere; the ionospheric cutoff-frequency could be well below 10 MHz, even in the day time. In anticipation of this upcoming opportunity, we propose and present in this paper a concept of an experimental ULW radio array, with the intention of setting it up in Inner Mongolia, China. This ULW facility will use the infrastructure of the currently operational Miangtu Spectra Radio Heliograph (MUSER). The proposed ULW array covers the frequency range from 1 to 72 MHz. This experimental array will be used for exploratory studies of celestial radio emission in the ULW range of the spectrum.

We determined the chemical and kinematic properties of the Galactic thin and thick disk using a sample of 307,246 A/F/G/K-type giant stars from the LAMOST spectroscopic survey and Gaia DR2 survey. Our study found that the thick disk globally exhibits no metallicity radial gradient, but the inner disk ($R \le 8$ kpc) and the outer disk ($R>8$ kpc) have different gradients when they are studied separately. The thin disk also shows two different metallicity radial gradients for the inner disk and the outer disk, and has steep metallicity vertical gradient of d[Fe/H]/d$|z|$ $=-0.12\pm0.0007$ dex kpc$^{-1}$, but it becomes flat when it is measured at increasing radial distance, while the metallicity radial gradient becomes weaker with increasing vertical distance. Adopting a galaxy potential model, we derived the orbital eccentricity of sample stars and found a downtrend of average eccentricity with increasing metallicity for the thick disk. The variation of the rotation velocity with the metallicity shows a positive gradient for the thick disk stars and a negative one for the thin disk stars. Comparisons of our observed results with models of disk formation suggest that radial migration could have influenced the chemical evolution of the thin disk. The formation of the thick disk could be affected by more than one processes: the accretion model could play an indispensable role, while other formation mechanisms, such as the radial migration or heating scenario model could also have a contribution.

In this paper, we review our current efforts to test General Relativity in the strong field regime by studying the reflection spectrum of supermassive black holes. So far we have analyzed 11 sources with observations of NuSTAR, Suzaku, Swift, and XMM-Newton. Our results are consistent with general relativity, according to which the spacetime metric around astrophysical black holes should be well approximated by the Kerr solution. We discuss the systematic uncertainties in our model and we present a preliminary study on the impact of some of them on the measurement of the spacetime metric.

We present observations of the C-band $1_{10}-1_{11}$ (4.8 GHz) and Ku-band $2_{11}-2_{12}$ (14.5 GHz) K-doublet lines of H$_2$CO and the C-band $1_{10}-1_{11}$ (4.6 GHz) line of H$_2$$^{13}$CO toward a large sample of Galactic molecular clouds, through the Shanghai Tianma 65-m radio telescope (TMRT). Our sample with 112 sources includes strong H$_2$CO sources from the TMRT molecular line survey at C-band and other known H$_2$CO sources. All three lines are detected toward 38 objects (43 radial velocity components) yielding a detection rate of 34\%. Complementary observations of their continuum emission at both C- and Ku-bands were performed. Combining spectral line parameters and continuum data, we calculate the column densities, the optical depths and the isotope ratio H$_2$$^{12}$CO/H$_2$$^{13}$CO for each source. To evaluate photon trapping caused by sometimes significant opacities in the main isotopologue's rotational mm-wave lines connecting our measured K-doublets, and to obtain $^{12}$C/$^{13}$C abundance ratios, we used the RADEX non-LTE model accounting for radiative transfer effects. This implied the use of the new collision rates from \citet{Wiesenfeld2013}. Also implementing distance values from trigonometric parallax measurements for our sources, we obtain a linear fit of $^{12}$C/$^{13}$C = (5.08$\pm$1.10)D$_{GC}$ + (11.86$\pm$6.60), with a correlation coefficient of 0.58. D$_{GC}$ refers to Galactocentric distances. Our $^{12}$C/$^{13}$C ratios agree very well with the ones deduced from CN and C$^{18}$O but are lower than those previously reported on the basis of H$_2$CO, tending to suggest that the bulk of the H$_2$CO in our sources was formed on dust grain mantles and not in the gas phase.

We present the study of a sample of seven "bare" active galactic nuclei (AGN) observed with Suzaku. We interpret the spectrum of these sources with a relativistic reflection component and we employ our model RELXILL_NK to test the Kerr nature of their supermassive black holes. We constrain the Johannsen deformation parameters $\alpha_{13}$ and $\alpha_{22}$, in which the Kerr metric is recovered when $\alpha_{13} = \alpha_{22} = 0$. All our measurements are consistent with the hypothesis that the spacetime geometry around these supermassive objects is described by the Kerr solution. For some sources, we obtain quite strong constraints on $\alpha_{13}$ and $\alpha_{22}$ when compared to those found in our previous studies. We discuss the systematic uncertainties in our tests and the implications of our results.

Final results are reported from operation of the PICO-60 C$_3$F$_8$ dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. The chamber was operated at thermodynamic thresholds as low as 1.2 keV without loss of stability. A new blind 1404-kg-day exposure at 2.45 keV threshold was acquired with approximately the same expected total background rate as the previous 1167-kg-day exposure at 3.3 keV. This increased exposure is enabled in part by a new optical tracking analysis to better identify events near detector walls, permitting a larger fiducial volume. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 2.5 $\times$ 10$^{-41}$ cm$^2$ for a 25 GeV WIMP, and improve on previous PICO results for 3-5 GeV WIMPs by an order of magnitude.

We present a study of decimetric radio activity, using the first high time cadence (0.5 s) images from the ${\textit{Giant Meterwave Radio Telescope}}$ (${\textit{GMRT}}$) at 610 MHz, associated with ${\textit{GOES}}$ C1.4 and M1.0 class solar flares, and a coronal mass ejection (CME) onset that occurred on 20 June 2015. The high spatial resolution images from ${\textit{GMRT}}$ show a strong radio source during the C1.4 flare, located ~500$''$ away from the flaring site with no corresponding bright footpoints or coronal features nearby. In contrast, however, strong radio sources are found near the flaring site during the M1.0 flare and around the CME onset time. Weak radio sources, located near the flaring site, are also found during the maximum of the C1.4 flare activity, which show a temporal association with metric type III bursts identified by the Solar Broadband Radio Spectrometer at Yunnan Astronomical Observatory. Based on a multi-wavelength analysis and magnetic potential field source surface extrapolations, we suggest that the source electrons of ${\textit{GMRT}}$ radio sources and metric type III bursts were originated from a common electron acceleration site. We also show that the strong ${\textit{GMRT}}$ radio source is generated by a coherent emission process and its apparent location far from the flaring site is possibly due to the wave-ducting effect.

The light curves of solar-type stars present both periodic fluctuation and flare spikes. The gradual periodic fluctuation is interpreted as the rotational modulation of magnetic features on the stellar surface and is used to deduce magnetic feature activity properties. The flare spikes in light curves are used to derive flare activity properties. In this paper, we analyze the light curve data of three solar-type stars (KIC 6034120, KIC 3118883, and KIC 10528093) observed with Kepler space telescope and investigate the relationship between their magnetic feature activities and flare activities. The analysis shows that: (1) both the magnetic feature activity and the flare activity exhibit long-term variations as the Sun does; (2) unlike the Sun, the long-term variations of magnetic feature activity and flare activity are not in phase with each other; (3) the analysis of star KIC 6034120 suggests that the long-term variations of magnetic feature activity and flare activity have a similar cycle length. Our analysis and results indicate that the magnetic features that dominate rotational modulation and the flares possibly have different source regions, although they may be influenced by the magnetic field generated through a same dynamo process.

In this paper, we use the equation of state based on a modification of 2+1 flavors Nambu-Jona- Lasinio (NJL) model to study the quark matter of hybrid stars. For comparison, we utilize five EOSs of the relativistic mean-field (RMF) model to describe the hadronic phase. With the 3- window crossover interpolation approach, we try to construct relatively soft hybrid EOSs but find the maximum masses of hybrid stars do not differ too much. The results are quite close to two solar mass, which is consistent with the mass constraint of PSR J0348+0432. Furthermore it is noteworthy that the heaviest stable stars have central densities higher than that of the deconfinement transition thus suggesting a pure quark core in the hybrid star.

We consider a field theory model of coupled dark energy which treats dark energy as a three-form field and dark matter as a spinor field. By assuming the effective mass of dark matter as a power-law function of the three-form field and neglecting the potential term of dark energy, we obtain three solutions of the autonomous system of evolution equations, including a de Sitter attractor, a tracking solution and an approximate solution. To understand the strength of the coupling, we confront the model with the latest Type Ia Supernova (SN Ia), Baryon Acoustic Oscillations (BAO) and Cosmic Microwave Backround (CMB) radiation observations, with the conclusion that the combination of these three databases marginalized over the present dark matter density parameter $\Omega_{m0}$ and the present three-form field $\kappa X_{0}$ gives stringent constraints on the coupling constant, $-0.017< \lambda <0.047$ ($2\sigma$ confidence level), by which we give out the model applicable parameter range.

In radio astronomy, the Ultra-Long Wavelengths (ULW) regime of longer than 10 m (frequencies below 30 MHz), remains the last virtually unexplored window of the celestial electromagnetic spectrum. The strength of the science case for extending radio astronomy into the ULW window is growing. However, the opaqueness of the Earth's ionosphere makes ULW observations by ground-based facilities practically impossible. Furthermore, the ULW spectrum is full of anthropogenic radio frequency interference (RFI). The only radical solution for both problems is in placing an ULW astronomy facility in space. We present a concept of a key element of a space-borne ULW array facility, an antenna that addresses radio astronomical specifications. A tripole-type antenna and amplifier are analysed as a solution for ULW implementation. A receiver system with a low power dissipation is discussed as well. The active antenna is optimized to operate at the noise level defined by the celestial emission in the frequency band 1 - 30 MHz. Field experiments with a prototype tripole antenna enabled estimates of the system noise temperature. They indicated that the proposed concept meets the requirements of a space-borne ULW array facility.

Solar radio type III bursts are believed to be the most sensitive signature of near-relativistic electron beam propagation in the corona. A solar radio type IIIb-III pair burst with fine frequency structures, observed by the Low Frequency Array (LOFAR) with high temporal ($\sim10$ ms) and spectral (12.5 kHz) resolutions at 30 - 80 MHz, is presented. The observations show that the type III burst consists of many striae, which have a frequency scale of about 0.1 MHz in both the fundamental (plasma) and the harmonic (double plasma) emission. We investigate the effects of background density fluctuations based on the observation of striae structure to estimate the density perturbation in solar corona. It is found that the spectral index of the density fluctuation spectrum is about $-1.7$, and the characteristic spatial scale of the density perturbation is around $700$ km. This spectral index is very close to a Kolmogorov turbulence spectral index of $-5/3$, consistent with a turbulent cascade. This fact indicates that the coronal turbulence may play the important role of modulating the time structures of solar radio type III bursts, and the fine structure of radio type III bursts could provide a useful and unique tool to diagnose the turbulence in the solar corona.

We examine the decay vacuum model with a parameter $\epsilon$ that indicates the vacuum energy decay rate. By constraining this model with cosmic microwave background radiation, baryon acoustic oscillation, type Ia supernovae and 30 H(z) cosmic chronometer data points, we find that $\epsilon=-0.0003\pm0.00024$ with the best fitted $\chi^{2}$ value slightly smaller than that in the $\Lambda$CDM model. A negative value of $\epsilon$ suggesting dark matter energy decay into vacuum energy. We also obtain the Hubble constant $H_{0}=68.05\pm0.56$ that can alleviate the current $H_{0}$ tension between the local observation by the Hubble Space Telescope and the global measurement by the Planck Satellite. Using the effective equation of state formalism, we find this model is quintessence-like.

We apply a novel pre-flare tracking of sunspot groups towards improving the estimation of flare onset time by focusing on the evolution of the 3D magnetic field construction of AR 11429. The 3D magnetic structure is based on potential field extrapolation encompassing a vertical range from the photosphere through the chromosphere and transition region into the low corona. The basis of our proxy measure of activity prediction is the so-called weighted horizontal gradient of magnetic field (WG_M) defined between spots of opposite polarities close to the polarity inversion line of an active region. The temporal variation of the distance of the barycenter of the opposite polarities is also found to possess potentially important diagnostic information about the flare onset time estimation as function of height similar to its counterpart introduced initially in an application at the photosphere only in Korsos et al. (2015). We apply the photospheric pre-flare behavioural patterns of sunspot groups to the evolution of their associated 3D-constructed AR 11429 as function of height. We found that at a certain height in the lower solar atmosphere the onset time may be estimated much earlier than at the photosphere or at any other heights. Therefore, we present a tool and recipe that may potentially identify the optimum height for flare prognostic in the solar atmosphere allowing to improve our flare prediction capability and capacity.

As a dedicated solar radio interferometer, the MingantU SpEctral RadioHeliograph (MUSER) generates massive observational data in the frequency range of 400 MHz -- 15 GHz. High-performance imaging forms a significantly important aspect of MUSER's massive data processing requirements. In this study, we implement a practical high-performance imaging pipeline for MUSER data processing. At first, the specifications of the MUSER are introduced and its imaging requirements are analyzed. Referring to the most commonly used radio astronomy software such as CASA and MIRIAD, we then implement a high-performance imaging pipeline based on the Graphics Processing Unit (GPU) technology with respect to the current operational status of the MUSER. A series of critical algorithms and their pseudo codes, i.e., detection of the solar disk and sky brightness, automatic centering of the solar disk and estimation of the number of iterations for clean algorithms, are proposed in detail. The preliminary experimental results indicate that the proposed imaging approach significantly increases the processing performance of MUSER and generates images with high-quality, which can meet the requirements of the MUSER data processing.

The reason for the occurrence of different elements of the fine structure of solar radio bursts in the decimeter and centimeter wavelength ranges has been determined based on all available data from terrestrial and satellite observations. In some phenomena, fast pulsations, a zebra structure, fiber bursts, and spikes have been observed almost simultaneously. Two phenomena have been selected to show that the pulsations of radio emission are caused by particles accelerated in the magnetic reconnection region and that the zebra structure is excited in a source, such as a magnetic trap for fast particles. The complex combination of unusual fiber bursts, zebra structure, and spikes in the phenomenon on December 1, 2004, is associated with a single source, a magnetic island formed after a coronal mass ejection.

Various solar features can be seen on maps of the Sun in the mm and sub-mm wavelength range. The recently installed Atacama Large Millimeter/submillimeter Array (ALMA) is capable of observing the Sun in that wavelength range with an unprecedented spatial, temporal and spectral resolution. To interpret solar observations with ALMA the first important step is to compare ALMA maps with simultaneous images of the Sun recorded in other spectral ranges. First we identify different structures in the solar atmosphere seen in the optical, IR and EUV parts of the spectrum (quiet Sun (QS), active regions (AR), prominences on the disc, magnetic inversion lines (IL), coronal holes (CH) and coronal bright points (CBPs)) in a full disc solar ALMA image. The second aim is to measure the intensities (brightness temperatures) of those structures and compare them with the corresponding QS level. A full disc solar image at 1.21 mm obtained on December 18, 2015 during a CSV-EOC campaign with ALMA is calibrated and compared with full disc solar images from the same day in H\alpha, in He I 1083 nm core, and with SDO images (AIA at 170 nm, 30.4 nm, 21.1 nm, 19.3 nm, and 17.1 nm and HMI magnetogram). The brightness temperatures of various structures are determined by averaging over corresponding regions of interest in the ALMA image. Positions of the QS, ARs, prominences on the disc, ILs, CHs and CBPs are identified in the ALMA image. At 1.21 mm ARs appear as bright areas (but sunspots are dark), while prominences on the disc and CHs are not discernible from the QS background, although having slightly less intensity than surrounding QS regions. ILs appear as large, elongated dark structures and CBPs correspond to ALMA bright points. These results are in general agreement with sparse earlier measurements at similar wavelengths. The identification of CBPs represents the most important new result.

New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less than 1 event per month. A blind analysis of an efficiency-corrected 1167-kg-day exposure at a 3.3-keV thermodynamic threshold reveals no single-scattering nuclear recoil candidates, consistent with the predicted background. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 3.4 $\times$ 10$^{-41}$ cm$^2$ for a 30-GeV$\thinspace$c$^{-2}$ WIMP, more than one order of magnitude improvement from previous PICO results.

The Atacama Large Millimeter-submillimeter Array (ALMA) radio telescope has commenced science observations of the Sun starting in late 2016. Since the Sun is much larger than the field of view of individual ALMA dishes, the ALMA interferometer is unable to measure the background level of solar emission when observing the solar disk. The absolute temperature scale is a critical measurement for much of ALMA solar science, including the understanding of energy transfer through the solar atmosphere, the properties of prominences, and the study of shock heating in the chromosphere. In order to provide an absolute temperature scale, ALMA solar observing will take advantage of the remarkable fast-scanning capabilities of the ALMA 12m dishes to make single-dish maps of the full Sun. This article reports on the results of an extensive commissioning effort to optimize the mapping procedure, and it describes the nature of the resulting data. Amplitude calibration is discussed in detail: a path that utilizes the two loads in the ALMA calibration system as well as sky measurements is described and applied to commissioning data. Inspection of a large number of single-dish datasets shows significant variation in the resulting temperatures, and based on the temperature distributions we derive quiet-Sun values at disk center of 7300 K at lambda=3 mm and 5900 K at lambda=1.3 mm. These values have statistical uncertainties of order 100 K, but systematic uncertainties in the temperature scale that may be significantly larger. Example images are presented from two periods with very different levels of solar activity. At a resolution of order 25 arcsec, the 1.3 mm wavelength images show temperatures on the disk that vary over about a 2000 K range.

Observations of the Sun at millimeter and submillimeter wavelengths offer a unique probe into the structure, dynamics, and heating of the chromosphere; the structure of sunspots; the formation and eruption of prominences and filaments; and energetic phenomena such as jets and flares. High-resolution observations of the Sun at millimeter and submillimeter wavelengths are challenging due to the intense, extended, low- contrast, and dynamic nature of emission from the quiet Sun, and the extremely intense and variable nature of emissions associated with energetic phenomena. The Atacama Large Millimeter/submillimeter Array (ALMA) was designed with solar observations in mind. The requirements for solar observations are significantly different from observations of sidereal sources and special measures are necessary to successfully carry out this type of observations. We describe the commissioning efforts that enable the use of two frequency bands, the 3 mm band (Band 3) and the 1.25 mm band (Band 6), for continuum interferometric-imaging observations of the Sun with ALMA. Examples of high-resolution synthesized images obtained using the newly commissioned modes during the solar commissioning campaign held in December 2015 are presented. Although only 30 of the eventual 66 ALMA antennas were used for the campaign, the solar images synthesized from the ALMA commissioning data reveal new features of the solar atmosphere that demonstrate the potential power of ALMA solar observations. The ongoing expansion of ALMA and solar-commissioning efforts will continue to enable new and unique solar observing capabilities.

We propose a new pressure-parametrization model to explain the accelerated expansion of the late-time Universe by considering the dark contents (dark matter and dark energy) as a unified dark fluid. To realize this model more physically, we reconstruct it with the quintessence and phantom scalar fields, respectively. We use the recent cosmological data to constrain this model, distinguish it from the standard cosmological model and find that the value of the Hubble constant $H_0=68.34^{+0.53}_{-0.92}$ supports the global measurement by the Planck satellite at the $1\sigma$ confidence level.

Sausage modes that are axisymmetric fast magnetoacoustic oscillations of solar coronal loops are characterized by variation of the plasma density and magnetic field, and hence cause time variations of the electron plasma frequency and cyclotron frequency. The latter parameters determine the condition for the double plasma resonance (DPR), which is responsible for the appearance of zebra-pattern (ZP) structures in time spectra of solar type IV radio bursts. We perform numerical simulations of standing and propagating sausage oscillations in a coronal loop modeled as a straight, field-aligned plasma slab, and determine the time variation of the DPR layer locations. Instant values of the plasma density and magnetic field at the DPR layers allowed us to construct skeletons of the time variation of ZP stripes in radio spectra. In the presence of a sausage oscillation, the ZP structures are shown to have characteristic wiggles with the time period prescribed by the sausage oscillation. Standing and propagating sausage oscillations are found to have different signatures in ZP patterns. We conclude that ZP wiggles can be used for the detection of short-period sausage oscillations and the exploitation of their seismological potential.

This work reports a peculiar and interesting train of microwave type III pair bursts in the impulsive rising phase of a solar flare on 2011 September 26. The observations include radio spectrometers at frequency of 0.80 - 2.00 GHz, hard X-ray (RHESSI and FERMI), EUV images of SWAP/PROBA-2 and magnetogram of HMI/SDO. By using a recently developed method (Tan et al. 2016a), we diagnosed the plasma density, temperature, plasma beta, magnetic field near the source region, the energy of energetic electrons and the distance between the acceleration region and the emission start sites of type III bursts. From the diagnostics, we find that: (1) The plasma density, temperature, magnetic field, and the distance between the acceleration region and the emission start sites almost have no obvious variations during the period of type III pair trains, while the energy of electrons has an obvious peak value which is consistent to the hard X-ray emission. (2) The plasma beta is much higher than an unity showing a highly dynamic process near the emission start site of type III bursts. (3) Although the reversed-slope type III branches drift slower at one order of magnitude than that of the normal type III branches, the related downgoing and upgoing electrons still could have same order of magnitude of energy. These facts indicate that both of the upgoing and downgoing electrons are possibly accelerated by similar mechanism and in a small source region. This diagnostics can help us to understand the microphysics in the source region of solar bursts.

The variations of the 3-D coronal magnetic fields associated with the X3.4-class flare of active region 10930 are studied in this paper. The coronal magnetic field data are reconstructed from the photospheric vector magnetograms obtained by the Hinode satellite and using the nonlinear force-free field extrapolation method developed in our previous work (He et al., 2011). The 3-D force-free factor $\alpha$, 3-D current density, and 3-D magnetic energy density are employed to analyze the coronal data. The distributions of $\alpha$ and current density reveal a prominent magnetic connectivity with strong negative $\alpha$ values and strong current density before the flare. This magnetic connectivity extends along the main polarity inversion line and is found to be totally broken after the flare. The distribution variation of magnetic energy density reveals the redistribution of magnetic energy before and after the flare. In the lower space of the modeling volume the increase of magnetic energy dominates, and in the higher space the decrease of energy dominates. The comparison with the flare onset imaging observation exhibits that the breaking site of the magnetic connectivity and site with the highest values of energy density increase coincide with the location of flare initial eruption. We conclude that a cramped positive $\alpha$ region appearing in the photosphere causes the breaking of the magnetic connectivity. A scenario for flare initial eruption is proposed in which the Lorentz force acting on the isolated electric current at the magnetic connectivity breaking site lifts the associated plasmas and causes the initial ejection.

A detailed investigation was carried out to understand a prolonged (~44 hours) weakly southward interplanetary magnetic field (IMF-Bz) condition during May 02--04, 1998. In-situ observations, during the period, showed the passage of an expanding magnetic cloud embedded in an interplanetary coronal mass ejection (ICME), followed up by a shock and an interplanetary discontinuity driven by another ICME. It is the arrival of the ICMEs and the upfront shocks that cause the prolonged southward IMF-Bz condition. The magnetic configuration of the source regions of the IMF associated with the ICME interval were also examined, which showed open magnetic field structures, emanating from a small active region on the north of the heliospheric current sheet (HCS). The structures remained constantly to the north of the HCS, both on April 29 and May 01, suggesting no change in their polarity. The draping of these outward directed radial field lines around the propagating CMEs in the shocked plasma explained the observed polarity changes of the IMF-Bz at 1 AU. In addition, multiple enhancements were also detected in the geomagnetic field variations, which showed a distinct one-to-one correspondence with the density pulses observed at 1 AU, during 0700--1700 UT on May 03. The spectral analysis of both the variations showed the same discrete frequencies of 0.48, 0.65 and 0.75 mHz, demonstrating that the solar wind density enhancements can cause detectable global geomagnetic disturbances. The observations, thus, provide a deeper insight into the possible causes and geomagnetic consequences of a prolonged weakly southward IMF-Bz condition.

After more than half a century of community support related to the science of "solar activity'', IAU's Commission 10 was formally discontinued in 2015, to be succeeded by C.E2 with the same area of responsibility. On this occasion, we look back at the growth of the scientific disciplines involved around the world over almost a full century. Solar activity and fields of research looking into the related physics of the heliosphere continue to be vibrant and growing, with currently over 2,000 refereed publications appearing per year from over 4,000 unique authors, publishing in dozens of distinct journals and meeting in dozens of workshops and conferences each year. The size of the rapidly growing community and of the observational and computational data volumes, along with the multitude of connections into other branches of astrophysics, pose significant challenges; aspects of these challenges are beginning to be addressed through, among others, the development of new systems of literature reviews, machine-searchable archives for data and publications, and virtual observatories. As customary in these reports, we highlight some of the research topics that have seen particular interest over the most recent triennium, specifically active-region magnetic fields, coronal thermal structure, coronal seismology, flares and eruptions, and the variability of solar activity on long time scales. We close with a collection of developments, discoveries, and surprises that illustrate the range and dynamics of the discipline.

Observations show that there is a proton spectral "break" with E$_{break}$ at 1-10MeV in some large CME-driven shocks. Theoretical model usually attribute this phenomenon to a diffusive shock acceleration. However, the underlying physics of the shock acceleration still remains uncertain. Although previous numerical models can hardly predict this "break" due to either high computational expense or shortcomings of current models, the present paper focuses on simulating this energy spectrum in converged two shocks by Monte Carlo numerical method. Considering the Dec 13 2006 CME-driven shock interaction with an Earth bow shock, we examine whether the energy spectral "break" could occur on an interaction between two shocks. As result, we indeed obtain the maximum proton energy up to 10MeV, which is the premise to investigate the existence of the energy spectral "break". Unexpectedly, we further find a proton spectral "break" appears distinctly at the energy $\sim$5MeV.

The measurement of positions and sizes of radio sources in the observations of solar radio spectral fine structures in an M6.5 flare on April 11, 2013 were observed simultaneously by several radio instruments at four different observatories: Chinese Solar Broadband Radio Spectrometers at Huairou (SBRS/Huairou), Ondrejov Radio spectrograph in the Czech Republic (ORSC/Ondrejov), Badary Broadband Microwave spectropolarimeter (BMS/Irkutsk), and spectrograph/IZMIRAN (Moscow, Troitsk). The fine structures include microwave zebra patterns (ZP), fast pulsations, and fibers. They were observed during the flare brightening located at the tops of a loop arcade. The dynamics of the polarization was associated with the motion of the flare exciter, which was observed in EUV images at 171A and 131A (SDO/AIA). Combining magnetograms observed by the SDO Helioseismic and Magnetic Imager (HMI) with the homologous assumption of EUV flare brightening and ZP bursts, we deduced that the observed ZPs correspond to the ordinary radio emission mode. However, future analysis needs to verify the assumption that zebra radio sources are really related to a closed magnetic loop, and are located at lower heights in the solar atmosphere than the source of pulsations.

Observations from multiple spacecraft show that there are energy spectral "breaks" at 1-10MeV in some large CME-driven shocks. However, numerical models can hardly simulate this property due to high computational expense. The present paper focuses on analyzing these energy spectral "breaks" by Monte Carlo particle simulations of an isolated CME-driven shock. Taking the Dec 14 2006 CME-driven shock as an example, we investigate the formation of this energy spectral property. For this purpose, we apply different values for the scattering time in our isolated shock model to obtain the highest energy "tails", which can potentially exceed the "break" energy range. However, we have not found the highest energy "tails" beyond the "break" energy range, but instead find that the highest energy "tails" reach saturation near the range of energy at 5MeV. So, we believe that there exists an energy spectral "cut off" in an isolated shock. If there is no interaction with another shock, there would not be formation of the energy spectral "break" property.

This work presents a systematic investigation of the influence of weather conditions on the calibration errors by using Gaussian fitness, least chi-square linear fitness and wavelet transform to analyze the calibration coefficients from observations of the Chinese Solar Broadband Radio Spectrometers (at frequency bands of 1.0-2.0 GHz, 2.6-3.8 GHz, and 5.2-7.6 GHz) during 1997-2007. We found that calibration coefficients are influenced by the local air temperature. Considering the temperature correction, the calibration error will reduce by about $10\%-20\%$ at 2800 MHz. Based on the above investigation and the calibration corrections, we further study the radio emission of the quiet-Sun by using an appropriate hybrid model of the quiet-Sun atmosphere. The results indicate that the numerical flux of the hybrid model is much closer to the observation flux than that of other ones.

Chinese Spectral Radio Heliograph (CSRH) is an advanced aperture synthesis solar radio heliograph, developed by National Astronomical Observatories, Chinese Academy of Sciences independently. It consists of 100 reflector antennas, which are grouped into two antenna arrays (CSRH-I and CSRH-II) for low and high frequency bands respectively. The frequency band of CSRH-I is 0.4-2GHz and for CSRH-II, the frequency band is 2-15GHz. In the antenna and feed system, CSRH uses an Eleven feed to receive signals coming from the Sun, the radiation pattern with lower side lobe and back lobe of the feed is well radiated. The characteristics of gain G and antenna noise temperature T effect the quality of solar radio imaging. For CSRH, measured G is larger than 60 dBi and $ T $ is less than 120K, after CSRH-I was established, we have successfully captured a solar radio burst between 1.2-1.6GHz on November 12, 2010 through this instrument and this event was confirmed through the observation of Solar Broadband Radio Spectrometer (SBRS) at 2.84GHz and Geostationary Operational Environmental Satellite (GOES). In addition, an image obtained from CSRH-I clearly reveals the profile of the solar radio burst. The other observational work is the imaging of Fengyun-2E geosynchronous satellite which is assumed to be a point source. This data processing method indicates that, the method of deleting errors about dirty image could be used for processing other surface sources.

This paper presents an overview of results obtained during the CAWSES II period on the short term variability of the Sun and how it affects the near Earth space environment. CAWSES II was planned to examine the behavior of the solar terrestrial system as the solar activity climbed to its maximum phase in solar cycle 24. After a deep minimum following cycle 23, the Sun climbed to a very weak maximum in terms of the sunspot number in cycle 24 (MiniMax24), so many of the results presented here refer to this weak activity in comparison with cycle 23. The short term variability that has immediate consequence to Earth and geospace manifests as solar eruptions from closed field regions and high speed streams from coronal holes. Both electromagnetic (flares) and mass emissions (coronal mass ejections, CMEs) are involved in solar eruptions, while coronal holes result in high speed streams that collide with slow wind forming the so called corotating interaction regions (CIRs). Fast CMEs affect Earth via leading shocks accelerating energetic particles and creating large geomagnetic storms. CIRs and their trailing high speed streams (HSSs), on the other hand, are responsible for recurrent small geomagnetic storms and extended (days) of auroral zone activity, respectively. The latter lead to the acceleration of relativistic magnetospheric killer electrons. One of the major consequences of the weak solar activity is the altered physical state of the heliosphere that has serious implications for the shock-driving and storm causing properties of CMEs. Finally, a discussion is presented on extreme space weather events prompted by the 2012 July 23 super storm event that occurred on the backside of the Sun. Many of these studies were enabled by the simultaneous availability of remote-sensing and in situ observations from multiple vantage points with respect to the Sun Earth line.

Radiobursts exibiting fine structure observed over two years during the rising phase of Cycle 24 by the SBRS are analyzed. In five events zebra structure, various fiber bursts and fast pulsations were observed. Events on 15 and 24 February 2011 are of the greatest interest. The polarization of radio emission in all three cases is related to the ordinary wave mode of radio emission.Almost all events in the microwave range contain superfine structure. It is possible that each type of fine structure is excited by the same mechanism, and the broad variety of events is related to the dynamics of flare processes.

Solar Optical Telescope onboard Hinode observed a sunspot (AR 11836) with two light bridges (LBs) on 31 Aug 2013. We analysed a 2-hour \ion{Ca}{2} H emission intensity data set and detected strong 5-min oscillation power on both LBs and in the inner penumbra. The time-distance plot reveals that 5-min oscillation phase does not vary significantly along the thin bridge, indicating that the oscillations are likely to originate from the underneath. The slit taken along the central axis of the wide light bridge exhibits a standing wave feature. However, at the centre of the wide bridge, the 5-min oscillation power is found to be stronger than at its sides. Moreover, the time-distance plot across the wide bridge exhibits a herringbone pattern that indicates a counter-stream of two running waves originated at the bridge sides. Thus, the 5-min oscillations on the wide bridge also resemble the properties of running penumbral waves. The 5-min oscillations are suppressed in the umbra, while the 3-min oscillations occupy all three cores of the sunspot's umbra, separated by the LBs. The 3-min oscillations were found to be in phase at both sides of the LBs. It may indicate that either LBs do not affect umbral oscillations, or umbral oscillations at different umbral cores share the same source. Also, it indicates that LBs are rather shallow objects situated in the upper part of the umbra. We found that umbral flashes follow the life cycles of umbral oscillations with much larger amplitudes. They cannot propagate across LBs. Umbral flashes dominate the 3-min oscillation power within each core, however, they do not disrupt the phase of umbral oscillation.

Microwave emission with spectral zebra pattern structures (ZPs) is observed frequently in solar flares and the Crab pulsar. The previous observations show that ZP is only a structure overlapped on the underlying broadband continuum with slight increments and decrements. This work reports an extremely unusual strong ZP burst occurring just at the beginning of a solar flare observed simultaneously by two radio telescopes located in China and Czech Republic and by the extreme ultraviolet (EUV) telescope on board NASA's satellite Solar Dynamics Observatory on 2013 April 11. It is a very short and super strong explosion whose intensity exceeds several times that of the underlying flaring broadband continuum emission, lasting for just 18 s. EUV images show that the flare starts from several small flare bursting points (FBPs). There is a sudden EUV flash with extra enhancement in one of these FBPs during the ZP burst. Analysis indicates that the ZP burst accompanying EUV flash is an unusual explosion revealing a strong coherent process with rapid particle acceleration, violent energy release, and fast plasma heating simultaneously in a small region with short duration just at the beginning of the flare.

We present simulations for interferometer arrays in Earth orbit and on the lunar surface to guide the design and optimization of space-based Ultra-Long Wavelength missions, such as those of China's Chang'E program. We choose parameters and present simulations using simulated data to identify inter-dependencies and constraints on science and engineering parameters. A regolith model is created for the lunar surface array simulation, the results show that the lunar regolith will have an undesirable effect on the observation. We estimate data transmission requirement, calculate sensitivities for both cases, and discuss the trade-off between brightness temperature sensitivity and angular resolution for the Earth orbit array case.

It is of great importance to track the solar wind back to its photospheric source region and identify the related current sheets; this will provide key information for investigating the origin and predictions of the solar wind. We report a statistical study relating the photospheric footpoint motion and in-situ observation of current sheets in the solar wind. We used the potential force-free source--surface (PFSS) model and the daily synoptic charts to trace the solar wind back from 1 AU, as observed by the Wind spacecraft, to the solar surface. As the footpoints move along the solar surface we obtain a time series of the jump times between different points. These jumps can be within a cell and between adjacent cells. We obtained the distribution of the jump times and the distribution for a subset of the jump times in which only jumps between adjacent cells were counted. For both cases, the distributions clearly show two populations. These distributions are compared with the distribution of in-situ current sheets reported in an earlier work of Miao Peng and Li (2011). Its implications on the origin of the current sheets are discussed.

Quasi-periodic wiggles of microwave zebra pattern structures with period range from about 0.5 s to 1.5 s are found in a X-class solar flare on 2006 December 13 at the 2.6-3.8 GHz with the Chinese Solar Broadband Radio Spectrometer (SBRS/Huairou). Periodogram and correlation analysis show that the wiggles have two-three significant periodicities and almost in phase between stripes at different frequency. The Alfven speed estimated from the zebra pattern structures is about 700 Km/s. We obtain the spatial size of the waveguiding plasma structure to be about 1 Mm with the detected period of about 1 s. It suggests the ZP wiggles can be associated with the fast mag- netoacoustic oscillations in the flaring active region. The lack of a significant phase shift between wiggles of different stripes suggests that the ZP wiggles are caused by a standing sausage oscillation.

A 3-D coronal magnetic field is reconstructed for NOAA 11158 on Feb 14, 2011. A GPU-accelerated direct boundary integral equation (DBIE) method is implemented. This is about 1000 times faster than the original DBIE used on solar NLFFF modeling. Using the SDO/HMI vector magnetogram as the bottom boundary condition, the reconstructed magnetic field lines are compared with the projected EUV loop structures from different views three-dimensionally by SDO/AIA and STEREO A/B spacecraft simultaneously for the first time. They show very good agreement so that the topological configurations of the magnetic fields can be analyzed, thus its role in the flare process of the active region can be better understood. A quantitative comparison with some stereoscopically reconstructed coronal loops shows that the present averaged misalignment angles are at the same order as the state-of-the-art results obtained with reconstructed coronal loops as prescribed conditions and better than other NLFFF methods. It is found that the observed coronal loop structures can be grouped into bundles of closed and open loops with some central bright coronal loops around the polarity inversion line (PIL). The reconstructed highly-shearing magnetic field lines agree very well with the low-lying S-shaped filament channel along PIL. They are in a pivot position to all other surrounding coronal structures, and a group of electric current lines co-aligned with the central bright EUV loops overlying the filament channel is also obtained. This central lower-lying magnetic field loop system must have played a key role in powering the flare. It should be noted that while a strand-like coronal feature along PIL may be related to the filament, one cannot simply attribute all the coronal bright features along PIL to manifestation of the filament without any stereoscopically information. It shows that DBIE is rigorous and effective.

Current sheets are ubiquitous in the solar wind.They are a major source of the solar wind MHD turbulence intermittency. They may result from non-linear interactions of the solar wind MHD turbulence or are the boundaries of flux tubes that originate from the solar surface. Some current sheets appear in pairs and are the boundaries of transient structures such as magnetic holes and reconnection exhausts, or the edges of pulsed Alfv\'{e}n waves. For an individual current sheet, discerning whether it is a flux tube boundary or due to non-linear interactions, or the boundary of a transient structure is difficult. In this work, using data from the {\sl Wind} spacecraft, we identify two three-current-sheet events. Detailed examination of these two events suggest that they are best explained by the flux tube crossing scenario. Our study provides a convincing evidence supporting the scenario that the solar wind consists of flux tubes where distinct plasmas reside.

The observation of $1.97\pm0.04$ solar-mass neutron-like star gives constraint on the equation of state (EOS) of cold, condensed matter. In this paper, the EOS for both pure quark star and hybrid star with a quark core described by quasi-particle model are considered. The parameters of quasi-particle model which affect the mass of both quark star and hybrid star can be constrained by the observation.

The morphologic properties of the magnetic networks during Carrington Rotations (CR) 1955 to 2091 (from 1999 to 2010) have been analyzed by applying the watershed algorithm to magnetograms observed by the Michelson Doppler Interferometer (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft. We find that the average area of magnetic cells on the solar surface at lower latitudes (within +-50 degree) are smaller than those at higher latitudes (beyond +-50 degree). Statistical analysis of these data indicates that the magnetic networks are of fractal in nature, and the average fractal dimension is D_f = 1.253+-0.011. We also find that both the fractal dimension and the size of the magnetic networks are anti-correlated with the sunspot area. This is perhaps because a strong magnetic field can suppress spatially modulated oscillation, compress the boundaries of network cells, leading to smoother cell boundaries. The fractal dimension of the cell deviates that predicted from an isobar of Kolmogorov homogeneous turbulence.

It is generally considered that the emission of microwave zebra pattern (ZP) structures requires high density and high temperatures, which is similar to the situation of the flaring region where primary energy releases. Therefore, the parameters analysis of ZPs may reveal the physical conditions of the flaring source region. This work investigates the variations of 74 microwave ZP structures observed by Chinese Solar Broadband Radio Spectrometer (SBRS/Huairou) at 2.6-3.8 GHz in 9 solar flares, and found that the ratio between the plasma density scale height LN and the magnetic field scale height LB in emission source displays a tendency of decrease during the flaring processes. The ratio LN/LB is about 3-5 before the maximum of flares. It decreases to about 2 after the maximum. The detailed analysis of three typical X-class flares implied that the variation of LN/LB during the flaring process is most likely due to the topological changes of the magnetic field in the flaring source region, and the step-wise decrease of LN/LB possibly reflects the magnetic field relaxation relative to the plasma density when the flaring energy released. This result may also constrain the solar flare modeling to some extent.

Observations of the interplanetary shock provide us with strong evidence of particle acceleration to multi-MeV energies, even up to GeV energy, in a solar flare or coronal mass ejection (CME). Diffusive shock acceleration is an efficient mechanism for particle acceleration. For investigating the shock structure, the energy injection and energy spectrum of a CME-driven shock, we perform dynamical Monte Carlo simulation of the 14-Dec-2006 CME-driven shock using an anisotropic scattering law. The simulated results of the shock fine structure, particle injection, and energy spectrum are presented. We find that our simulation results give a good fit to the observations from multiple spacecraft.

Based on the analysis of the microwave observations at frequency of 2.60 -- 3.80 GHz in a solar X1.3 flare event observed at Solar Broadband RadioSpectrometer in Huairou (SBRS/Huairou) on 2005 July 30, an interesting reversed drifting quasi-periodic pulsating structure (R-DPS) is confirmed. The R-DPS is mainly composed of two drifting pulsating components: one is a relatively slow very short-period pulsation (VSP) with period of about 130 -- 170 ms, the other is a relatively fast VSP with period of about 70 -- 80 ms. The R-DPS has a weak left-handed circular polarization. Based on the synthetic investigations of Reuven Ramaty High Energy Solar Spectroscopic Imaging (RHESSI) hard X-ray, Geostationary Operational Environmental Satellite (GOES) soft X-ray observation, and magnetic field extrapolation, we suggest the R-DPS possibly reflects flaring dynamic processes of the emission source regions.

The unusual radio burst of May 29, 2003 connected with the M1.5 flare in AR 10368 has been analyzed. It was observed by the Solar Broadband Radio Spectrometer (SBRS/Huairou station, Beijing) in the 5.2-7.6 GHz range. It proved to be only the third case of a neat zebra structure appearing among all observations at such high frequencies. Despite the short duration of the burst (25 s), it provided a wealth of data for studying the superfine structure with millisecond resolution (5 ms). We localize the site of emission sources in the flare region, estimate plasma parameters in the generation sites, and suggest applicable mechanisms for interpretating spikes and zebra-structure generation. Positions of radio bursts were obtained by the Siberian Solar Radio Telescope (SSRT) (5.7 GHz) and Nobeyama radioheliograph (NoRH) (17 GHz). The sources in intensity gravitated to tops of short loops at 17 GHz, and to long loops at 5.7 GHz. Short pulses at 17 GHz (with a temporal resolution of 100 ms) are registered in the R-polarized source over the N-magnetic polarity (extraordinary mode). Dynamic spectra show that all the emission comprised millisecond pulses (spikes) of 5-10 ms duration in the instantaneous band of 70 to 100 MHz, forming the superfine structure of different bursts, essentially in the form of fast or slow-drift fibers and various zebra-structure stripes. Five scales of zebra structures have been singled out. As the main mechanism for generating spikes (as the initial emission) we suggest the coalescence of plasma waves with whistlers in the pulse regime of interaction between whistlers and ion-sound waves. In this case one can explain the appearance of fibers and sporadic zebra-structure stripes exhibiting the frequency splitting.

Although diffusive shock acceleration (DSA) could be simulated by some well-established models, the assumption of the injection rate from the thermal particles to the superthermal population is still a contentious problem. But in the self-consistent Monte Carlo simulations, because of the prescribed scattering law instead of the assumption of the injected function, hence particle injection rate is intrinsically defined by the prescribed scattering law. We expect to examine the correlation of the energy injection with the prescribed multiple scattering angular distributions. According to the Rankine-Hugoniot conditions, the energy injection and the losses in the simulation system can directly decide the shock energy spectrum slope. By the simulations performed with multiple scattering law in the dynamical Monte Carlo model, the energy injection and energy loss functions are obtained. As results, the case applying anisotropic scattering law produce a small energy injection and large energy losses leading to a soft shock energy spectrum, the case applying isotropic scattering law produce a large energy injection and small energy losses leading to a hard shock energy spectrum.

According to the shock jump conditions, the total fluid's mass, momentum, and energy should be conserved in the entire simulation box. We perform the dynamical Monte Carlo simulations with the multiple scattering law for energy analysis. The various energy functions of time are obtained by monitoring the total particles' mass, momentum, and energy in the simulation box. In conclusion, the energy analysis indicates that the smaller energy losses in the prescribed scattering law are, the harder the energy spectrum produced is.

We analyse 3-min oscillations of microwave and EUV emission generated at different heights of a sunspot atmosphere, studying the amplitude and frequency modulation of the oscillations, and its relationship with the variation of the spatial structure of the oscillations. High-resolution data obtained with the Nobeyama Radioheliograph, TRACE and SDO/AIA are analysed with the use of the Pixelised Wavelet Filtering and wavelet skeleton techniques. 3-min oscillations in sunspots appear in the form of repetitive trains of the duration 8-20 min (13 min in average). The typical interval between the trains is 30-50 min. The oscillation trains are transient in frequency and power. We detected a repetitive frequency drifts of 3-min oscillations during the development of individual trains. Wavelet analysis shows three types of the frequency drift: positive, negative and fluctuations without drift. The start and end of the drifts coincide with the start time and end of the train. The comparative study of 3-min oscillations in the sequences of microwave and EUV images show the appearance of new sources of the oscillations in sunspots during the development of the trains. These structures can be interpreted as waveguides that channel upward propagating waves, responsible for 3-min oscillations. A possible explanation of the observed properties is the operation of two simultaneous factors: dispersive evolution of the upwardly-propagating wave pulses and the non-uniformity of the distribution of the oscillation power over the sunspot umbra with different wave sources corresponding to different magnetic flux tubes with different physical conditions and line-of-sight angles.

The so-called zebra structures in radio dynamic spectra, specifically their frequencies and frequency drifts of emission stripes, contain information on the plasma parameters in the coronal part of flare loops. This paper presents observations of zebra structures in a microwave range. Dynamic spectra were recorded by Chinese spectro-polarimeters in the frequency band close to the working frequencies of the Siberian Solar Radio Telescope. The emission sources are localized in the flare regions, and we are able to estimate the plasma parameters in the generation sites using X-ray data. The interpretation of the zebra structures in terms of the existing theories is discussed. The conclusion has been arrived that the preferred generation mechanism of zebra structures in the microwave range is the conversion of plasma waves to electromagnetic emission on the double plasma resonance surfaces distributed across a flare loop.

Zebra pattern structure (ZP) is the most intriguing fine structure on the dynamic spectrograph of solar microwave burst. On 15 February 2011, there erupts an X2.2 flare event on the solar disk, it is the first X-class flare since the solar Schwabe cycle 24. It is interesting that there are several microwave ZPs observed by the Chinese Solar Broadband Radiospectrometer (SBRS/Huairou) at frequency of 6.40 ~ 7.00 GHz (ZP1), 2.60 ~ 2.75 GHz (ZP2), and the Yunnan Solar Broadband Radio Spectrometer (SBRS/Yunnan) at frequency of 1.04 ~ 1.13 GHz (ZP3). The most important phenomena is the unusual high-frequency ZP structure (ZP1, up to 7.00 GHz) occurred in the early rising phase of the flare, and there are two ZP structure (ZP2, ZP3) with relative low frequencies occurred in the decay phase of the flare. By scrutinizing the current prevalent theoretical models of ZP structure generations, and comparing their estimated magnetic field strengths in the corresponding source regions, we suggest that the double plasma resonance model should be the most possible one for explaining the formation of microwave ZPs, which may derive the magnetic field strengths as about 230 - 345 G, 126 - 147 G, and 23 - 26 G in the source regions of ZP1, ZP2, and ZP3, respectively.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

We independently develop a simulation code following the previous dynamical Monte Carlo simulation of the diffusive shock acceleration under the isotropic scattering law during the scattering process, and the same results are obtained. Since the same results test the validity of the dynamical Monte Carlo method for simulating a collisionless shock, we extend the simulation toward including an anisotropic scattering law for further developing this dynamical Monte Carlo simulation. Under this extended anisotropic scattering law, a Gaussian distribution function is used to describe the variation of scattering angles in the particle's local frame. As a result, we obtain a series of different shock structures and evolutions in terms of the standard deviation values of the given Gaussian scattering angular distributions. We find that the total energy spectral index increases as the standard deviation value of the scattering angular distribution increases, but the subshock's energy spectral index decreases as the standard deviation value of the scattering angular distribution increases.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

Based the Alfven wave oscillation model (AWOM) and relativistic precession model (RPM) for twin kHz QPOs, we estimate the emission positions of most detected kHz QPOs to be at r=18+-3 km (R/15km) except Cir X-1 at r = 30\+-5 km (R/15km). For the proposed Keplerian frequency as an upper limit to kHz QPO, the spin effects in Kerr Spacetime are discussed, which have about a 5% (2%) modification for that of the Schwarzchild case for the spin frequency of 1000 (400) Hz.The application to the four typical QPO sources, Cir X-1, Sco X-1, SAX J1808.4-3658 and XTE 1807-294, is mentioned.

Based on the joint-observations of the radio broadband spectral emissions of solar eclipse on August 1, 2008 at Jiuquan (total eclipse) and Huairou (partial eclipse) at the frequencies of 2.00 -- 5.60 GHz (Jiuquan), 2.60 -- 3.80 GHZ (Chinese solar broadband radiospectrometer, SBRS/Huairou), and 5.20 -- 7.60 GHz (SBRS/Huairou), the authors assemble a successive series of broadband spectrum with a frequency of 2.60 -- 7.60 GHz to observe the solar eclipse synchronously. This is the first attempt to analyze the solar eclipse radio emission under the two telescopes located at different places with broadband frequencies in the periods of total and partial eclipse. With these analyses, the authors made a new semiempirical model of the coronal plasma density of the quiet Sun and made a comparison with the classic models.

We report the timing-analysis results obtained for RRAT J1819-1458 from regular timing observations at 1.54 GHz using the Urumqi 25 m radio telescope between 2007 April to 2008 March. RRAT J1819-1458 is a relatively young and highly magnetized neutron star discovered by its sporadic short bursts in the Parkes Multibeam Pulsar Survey data. In 94 hrs of observation data, we detected a total of 162 dispersed bursts of RRAT J1819-1458 with the signal-to-noise ratios (S/Ns) above 5-$\sigma$ threshold. Among them, 5 bursts clearly show two-component structure. The S/N of the strongest burst is 13.3. The source's DM measured through our data is 196.0 $\pm$0.4 pc cm$^{-3}$. The timing position, frequency and its first derivative were determined using standard pulsar timing techniques on the arrival times of these individual bursts. The accuracy of the solved rotating parameters are improved comparing with that in previous publication. Our timing position with 2-$\sigma$ error is consistent with the position of its X-ray counterpart CXOU J181934.1-145804. The effect of timing noise and the phase fluctuation of the individual short bursts on the timing residuals is briefly discussed. The distribution of the timing residuals is bimodal, which cannot be explained readily by timing irregularity.

An extended radio emission after a soft X-ray (SXR) maximum was detected in the active region NOAA 9077 by several observatories for the solar flare after 21:42 UT on July 10, 2000. Also some radio fine structures before the enduring radio emission were observed with the 1.0-2.0 GHz spectrometer of Beijing Astronomical Observatory (BAO) in the same time. We apply a shear-driven quadrupolar reconnection model (SQR) to analyze the fine structures and the related radio emission. We find that the footpoint shear motion of the flux loop is accompanied with the emerging up of the loop during the reconnection process. We tentatively interpret the extended radio emission as the nonthermal radiation caused by a new reconnection process between emerging flux loop and pre-existing overarching loop after the soft X-ray maximum.

Observations of the solar radio emission in the microwave frequency range show fine structures consisting of a number of the almost parallel narrow frequency bands. We interpret these bands as the cyclotron harmonics plasma emission. This emission is generated by the anisotropic electron beam, which excites longitudinal waves on the normal Doppler effect resonance. Then the longitudinal waves convert into the radio emission on the second harmonic of the longitudinal waves frequency and sometimes into the fundamental one. Estimations of the magnetic field strength made on the basis of such a model give the values of the magnetic field in the microwave burst sources as ~ 100-200 G. The estimations of density variations are also given.