The uncertainty of cosmological data generated from complex processes, such as observational Hubble parameter data (OHD) and the Type Ia supernova (SN Ia) data, cannot be accurately modeled by simple analytical probability distributions, e.g. Gaussian distribution. This necessitates the use of likelihood-free inference, which bypasses the direct calculation of likelihood. In this paper, we propose a new procedure to perform likelihood-free cosmological inference using two artificial neural networks (ANN), the Masked Autoregressive Flow (MAF) and denoising autoencoder (DAE). Our procedure is the first to use DAE to extract features from data in order to simplify the structure of MAF needed to estimate the posterior. Tested on simulated Hubble parameter data with a simple Gaussian likelihood, the procedure shows the capability of extracting feature from data and estimating posteriors without the need for tractable likelihood. We demonstrate that it can accurately approximate the real posterior, and achieves a performance comparable to the traditional MCMC method. We also discuss the application of the proposed procedure on OHD and Pantheon SN Ia data, and use them to constrain cosmological parameters from the non-flat $\Lambda$CDM model. With little prior knowledge added, we found constraints on $H_0, \Omega_m, \Omega_\Lambda$ similar to relevant work. In addition, this work is also the first to use Gaussian process in the procedure of OHD simulation.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The concept of the Solar Ring mission was gradually formed from L5/L4 mission concept, and the proposal of its pre-phase study was funded by the National Natural Science Foundation of China in November 2018 and then by the Strategic Priority Program of Chinese Academy of Sciences in space sciences in May 2019. Solar Ring mission will be the first attempt to routinely monitor and study the Sun and inner heliosphere from a full 360-degree perspective in the ecliptic plane. The current preliminary design of the Solar Ring mission is to deploy six spacecraft, grouped in three pairs, on a sub-AU orbit around the Sun. The two spacecraft in each group are separated by about 30 degrees and every two groups by about 120 degrees. This configuration with necessary science payloads will allow us to establish three unprecedented capabilities: (1) determine the photospheric vector magnetic field with unambiguity, (2) provide 360-degree maps of the Sun and the inner heliosphere routinely, and (3) resolve the solar wind structures at multiple scales and multiple longitudes. With these capabilities, the Solar Ring mission aims to address the origin of solar cycle, the origin of solar eruptions, the origin of solar wind structures and the origin of severe space weather events. The successful accomplishment of the mission will advance our understanding of the star and the space environment that hold our life and enhance our capability of expanding the next new territory of human.

Although the physical origin of prompt emission in gamma-ray bursts (GRBs) remains inconclusive, previous studies have considered the synchrotron radiation of relativistic electrons as a promising mechanism. These works usually adopted a invariable injection rate of electrons ($Q$) which may be discordant with that in a Poynting-flux dominated jet. In a Poynting-flux dominated jet (e.g., ICMART model, Zhang & Yan 2011), the number of magnetic reconnections occurred simultaneously may grow rapidly with time and results in an increase of $Q$ with time. This paper is dedicated to study the synchrotron radiation spectrum in this scenario. It is found that the radiation spectrum would obviously get harder if an increasing $Q$ is adopted and a Band-like radiation spectrum can be obtained if the increase of $Q$ is fast enough. The latter is related to the fact that a bump-shape rather than a power-law spectrum appears in the low-energy regime of the obtained electron spectrum. This effect can strongly harden the low-energy radiation spectrum. It indicates that an increasing $Q$ can help to alleviate the "fast-cooling problem" of synchrotron radiation for GRBs. Our studies also reveal that a Poynting-flux dominated jet with a large emission radius, a small length of the magnetic reconnection region, or a low-minimum energy of injected electron would prefer to form a Band-like radiation spectrum. We suggest that the Band spectrum found in GRBs may be the synchrotron emission of the electrons with a bump-shape distribution in its low-energy regime.

We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

We study the timing stability of three black widow pulsars, both in terms of their long-term spin evolution and their shorter-term orbital stability. The erratic timing behaviour and radio eclipses of the first two black widow pulsar systems discovered (PSRs B1957+20 and J2051$-$0827) was assumed to be representative for this class of pulsars. With several new black widow systems added to this population in the last decade, there are now several systems known that do not show these typical orbital variations or radio eclipses. We present timing solutions using 7$-$8 yrs of observations from four of the European Pulsar Timing Array telescopes for PSRs\,J0023+0923, J2214+3000 and J2234+0944, and confirm that two of these systems do not show any significant orbital variability over our observing time span, both in terms of secular or orbital parameters. The third pulsar PSR J0023+0923 shows orbital variability and we discuss the implications for the timing solution. Our results from the long-term timing of these pulsars provide several new or improved parameters compared to earlier works. We discuss our results regarding the stability of these pulsars, and the stability of the class of black widow pulsars in general, in the context of the binary parameters, and discuss the potential of the Roche-lobe filling factor of the companion star being an indicator for stability of these systems.

We report the discovery of an eclipsing binary millisecond pulsar in the globular cluster M92 (NGC6341) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J1717+4308A, or M92A, has a pulse frequency of 316.5~Hz (3.16~ms) and a dispersion measure of 35.45 pc cm$^{-3}$. The pulsar is a member of a binary system with an orbital period of 0.20~days around a low-mass companion which has a median mass of $\sim$0.18~\Ms. From observations so far, at least two eclipsing events have been observed in each orbit. The longer one lasted for ~5000~s in the orbital phase range 0.1--0.5. The other lasted for ~500~s and occurred between 1000--2000~s before or after the longer eclipsing event. The lengths of these two eclipsing events also change. These properties suggest that J1717+4308A is a ``red-back'' system with a low-mass main sequence or sub-giant companion. Timing observations of the pulsar and further searches of the data for additional pulsars are ongoing.

Monitoring chromospheric and photospheric indexes of magnetic activity can provide valuable information, especially the interaction between different parts of the atmosphere and their response to magnetic fields. We extract chromospheric indexes, S and Rhk+, for 59,816 stars from LAMOST spectra in the LAMOST-Kepler program, and photospheric index, Reff, for 5575 stars from Kepler light curves. The log Reff shows positive correlation with log Rhk+. We estimate the power-law indexes between Reff and Rhk+ for F-, G-, and K-type stars, respectively. We also confirm the dependence of both chromospheric and photospheric activity on stellar rotation. Ca II H and K emissions and photospheric variations generally decrease with increasing rotation periods for stars with rotation periods exceeding a few days. The power-law indexes in exponential decay regimes show different characteristics in the two activity-rotation relations. The updated largest sample including the activity proxies and reported rotation periods provides more information to understand the magnetic activity for cool stars.

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

A second peak in extreme ultraviolet sometimes appears during the gradual phase of solar flares, which is known as EUV late phase (ELP). Stereotypically ELP is associated with two separated sets of flaring loops with distinct sizes, and it has been debated whether ELP is caused by additional heating or extended plasma cooling in the longer loop system. Here we carry out a survey of 55 M-and-above GOES-class flares with ELP during 2010--2014. Based on the flare-ribbon morphology, these flares are categorized as circular-ribbon (19 events), two-ribbon (23 events), and complex-ribbon (13 events) flares. Among them, 22 events (40%) are associated with coronal mass ejections, while the rest are confined. An extreme ELP, with the late-phase peak exceeding the main-phase peak, is found in 48% of two-ribbon flares, 37% of circular-ribbon flares, and 31% of complex-ribbon flares, suggesting that additional heating is more likely present during ELP in two-ribbon than in circular-ribbon flares. Overall, cooling may be the dominant factor causing the delay of the ELP peak relative to the main-phase peak, because the loop system responsible for the ELP emission is generally larger than, and well separated from, that responsible for the main-phase emission. All but one of the circular-ribbon flares can be well explained by a composite "dome-plate" quasi-separatrix layer (QSL). Only half of these show a magnetic null point, with its fan and spine embedded in the dome and plate, respectively. The dome-plate QSL, therefore, is a general and robust structure characterizing circular-ribbon flares.

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

Rotating radio transients (RRATs) are peculiar astronomical objects whose emission mechanism remains under investigation. In this paper, we present observations of three RRATs, J1538+2345, J1854+0306 and J1913+1330, observed with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Specifically, we analyze the mean pulse profiles and temporal flux density evolutions of the RRATs. Owing to the high sensitivity of FAST, the derived burst rates of the three RRATs are higher than those in previous reports. RRAT J1854+0306 exhibited a time-dynamic mean pulse profile, whereas RRAT J1913+1330 showed distinct radiation and nulling segments on its pulse intensity trains. The mean pulse profile variation with frequency is also studied for RRAT J1538+2345 and RRAT J1913+1330, and the profiles at different frequencies could be well fitted with a cone-core model and a conal-beam model, respectively.

PSR B0919+06 is known for its abnormal emission phenomenon, where the pulse emission window occasionally shifts progressively in longitude and returns afterwards. The physical mechanism behind this phenomenon is still under investigation. In this paper, we present our ultra-wideband observation of this pulsar using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), with simultaneous measurements in the frequency ranges 280-780 MHz and 1250-1550 MHz. We have identified three abnormal events, each of which becomes less apparent as the frequency decreases. At 1400 MHz, the averaged profile slightly shifted after the first and third abnormal events, implying a relationship between abnormal event and profile variation. We also found a linear trend in the left-edge position of the averaged profiles between the first and third event as well as after the third event, suggesting the existence of a slow-drifting mode between the two major events. The second event has a comparatively small shift in phase and is thus categorized as a `small flare state'. During the third event, a sequence of approximately nine pulses was seen to significantly weaken in all frequency bands, likely associated with the pseudo-nulling observed at 150 MHz. A three-component de-composition analysis of the normal averaged profiles shows that the trailing component is dominant at our observing frequencies, while the centre component has a comparatively steeper spectrum. We found the overall flux density in an abnormal event to slightly differ from that in an ordinary state, and the difference shows a frequency dependence. A comparison of the normal, abnormal and dimmed averaged profile indicates that the leading component is likely to be stable in all states.

We present a catalog of stellar age and mass estimates for a sample of 640\,986 red giant branch (RGB) stars of the Galactic disk from the LAMOST Galactic Spectroscopic Survey (DR4). The RGB stars are distinguished from the red clump stars utilizing period spacing derived from the spectra with a machine learning method based on kernel principal component analysis (KPCA). Cross-validation suggests our method is capable of distinguishing RC from RGB stars with only 2 per cent contamination rate for stars with signal-to-noise ratio (SNR) higher than 50. The age and mass of these RGB stars are determined from their LAMOST spectra with KPCA method by taking the LAMOST - $Kepler$ giant stars having asteroseismic parameters and the LAMOST-TGAS sub-giant stars based on isochrones as training sets. Examinations suggest that the age and mass estimates of our RGB sample stars with SNR $>$ 30 have a median error of 30 per cent and 10 per cent, respectively. Stellar ages are found to exhibit positive vertical and negative radial gradients across the disk, and the age structure of the disk is strongly flared across the whole disk of $6<R<13$\,kpc. The data set demonstrates good correlations among stellar age, [Fe/H] and [$\alpha$/Fe]. There are two separate sequences in the [Fe/H] -- [$\alpha$/Fe] plane: a high--$\alpha$ sequence with stars older than $\sim$\,8\,Gyr and a low--$\alpha$ sequence composed of stars with ages covering the whole range of possible ages of stars. We also examine relations between age and kinematic parameters derived from the Gaia DR2 parallax and proper motions. Both the median value and dispersion of the orbital eccentricity are found to increase with age. The vertical angular momentum is found to fairly smoothly decrease with age from 2 to 12\,Gyr, with a rate of about $-$50\,kpc\,km\,s$^{-1}$\,Gyr$^{-1}$. A full table of the catalog is public available online.

We have studied 4265 giant pulses (GPs) from the millisecond pulsar B1937+21; the largest-ever sample gathered for this pulsar, in observations made with the Large European Array for Pulsars. The pulse energy distribution of GPs associated with the interpulse are well-described by a power law, with index $\alpha = -3.99 \pm 0.04$, while those associated with the main pulse are best-described by a broken power law, with the break occurring at $\sim7$ Jy $\mu$s, with power law indices $\alpha_{\text{low}} = -3.48 \pm 0.04$ and $\alpha_{\text{high}} = -2.10 \pm 0.09$. The modulation indices of the GP emission are measured, which are found to vary by $\sim0.5$ at pulse phases close to the centre of the GP phase distributions. We find the frequency-resolved structure of GPs to vary significantly, and in a manner that cannot be attributed to the interstellar medium influence on the observed pulses. We examine the distribution of polarisation fractions of the GPs and find no correlation between GP emission phase and fractional polarisation. We use the GPs to time PSR B1937+21 and although the achievable time of arrival precision of the GPs is approximately a factor of two greater than that of the average pulse profile, there is a negligible difference in the precision of the overall timing solution when using the GPs.

Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.

The Greenland Telescope project has recently participated in an experiment to image the supermassive black hole shadow at the center of M87 using Very Long Baseline Interferometry technique in April of 2018. The antenna consists of the 12-m ALMA North American prototype antenna that was modified to support two auxiliary side containers and to withstand an extremely cold environment. The telescope is currently at Thule Air Base in Greenland with the long-term goal to move the telescope over the Greenland ice sheet to Summit Station. The GLT currently has a single cryostat which houses three dual polarization receivers that cover 84-96 GHz, 213-243 GHz and 271-377 GHz bands. A hydrogen maser frequency source in conjunction with high frequency synthesizers are used to generate the local oscillator references for the receivers. The intermediate frequency outputs of each receiver cover 4-8 GHz and are heterodyned to baseband for digitization within a set of ROACH-2 units then formatted for recording onto Mark-6 data recorders. A separate set of ROACH-2 units operating in parallel provides the function of auto-correlation for real-time spectral analysis. Due to the stringent instrumental stability requirements for interferometry a diagnostic test system was incorporated into the design. Tying all of the above equipment together is the fiber optic system designed to operate in a low temperature environment and scalable to accommodate a larger distance between the control module and telescope for Summit Station. A report on the progress of the above electronics instrumentation system will be provided.

We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $\mu$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $\mu$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz.

Stellar systems composed of single, double, triple or high-order systems are rightfully regarded as the fundamental building blocks of the Milky Way. Binary stars play an important role in formation and evolution of the Galaxy. Through comparing the radial velocity variations from multi-epoch observations, we analyze the binary fraction of dwarf stars observed with the LAMOST. Effects of different model assumptions such as orbital period distributions on the estimate of binary fractions, are investigated. The results based on log-normal distribution of orbital periods reproduce the previous complete analyses better than the power-law distribution. We find that the binary fraction increases with $T_{\rm eff}$ and decreases with [Fe/H]. We first investigate the relation between $\alpha$-elements and binary fraction in such a large sample as the LAMOST. The old stars with high [$\alpha$/Fe] dominate higher binary fraction than young stars with low [$\alpha$/Fe]. At the same mass, former forming stars possess a higher binary fraction than newly forming ones, which may be related with the evolution of the Galaxy.

Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.

In this paper, we present the detailed analysis of recurrent homologous jets originating from an emerging negative magnetic flux at the edge of an Active Region. The observed jets show multi-thermal features. Their evolution shows high consistence with the characteristic parameters of the emerging flux, suggesting that with more free magnetic energy, the eruptions tend to be more violent, frequent and blowout-like. The average temperature, average electron number density and axial speed are found to be similar for different jets, indicating that they should have been formed by plasmas from similar origins. Statistical analysis of the jets and their footpoint region conditions reveals a strong positive relationship between the footpoint-region total 131 {\AA} intensity enhancement and jets' length/width. Stronger linearly positive relationships also exist between the total intensity enhancement/thermal energy of the footpoint regions and jets' mass/kinetic/thermal energy, with higher cross-correlation coefficients. All the above results, together, confirm the direct relationship between the magnetic reconnection and the jets, and validate the important role of magnetic reconnection in transporting large amount of free magnetic energy into jets. It is also suggested that there should be more free energy released during the magnetic reconnection of blowout than of standard jet events.

Obtaining accurate and precise masses and ages for large numbers of giant stars is of great importance for unraveling the assemblage history of the Galaxy. In this paper, we estimate masses and ages of 6940 red giant branch (RGB) stars with asteroseismic parameters deduced from \emph{Kepler} photometry and stellar atmospheric parameters derived from LAMOST spectra. The typical uncertainties of mass is a few per cent, and that of age is $\sim$\,20 per cent. The sample stars reveal two separate sequences in the age -- [$\alpha$/Fe] relation -- a high--$\alpha$ sequence with stars older than $\sim$\,8\,Gyr and a low--$\alpha$ sequence composed of stars with ages ranging from younger than 1\,Gyr to older than 11\,Gyr. We further investigate the feasibility of deducing ages and masses directly from LAMOST spectra with a machine learning method based on kernel based principal component analysis, taking a sub-sample of these RGB stars as a training data set. We demonstrate that ages thus derived achieve an accuracy of $\sim$\,24 per cent. We also explored the feasibility of estimating ages and masses based on the spectroscopically measured carbon and nitrogen abundances. The results are quite satisfactory and significantly improved compared to the previous studies.

Main sequence turn-off (MSTO) stars have advantages as indicators of Galactic evolution since their ages could be robustly estimated from atmospheric parameters. Hundreds of thousands of MSTO stars have been selected from the LAMOST Galactic sur- vey to study the evolution of the Galaxy, and it is vital to derive accurate stellar parameters. In this work, we select 150 MSTO star candidates from the MSTO stars sample of Xiang that have asteroseismic parameters and determine accurate stellar parameters for these stars combing the asteroseismic parameters deduced from the Kepler photometry and atmospheric parameters deduced from the LAMOST spectra.With this sample, we examine the age deter- mination as well as the contamination rate of the MSTO stars sample. A comparison of age between this work and Xiang shows a mean difference of 0.53 Gyr (7%) and a dispersion of 2.71 Gyr (28%). The results show that 79 of the candidates are MSTO stars, while the others are contaminations from either main sequence or sub-giant stars. The contamination rate for the oldest stars is much higher than that for the younger stars. The main cause for the high contamination rate is found to be the relatively large systematic bias in the LAMOST surface gravity estimates.

The Large European Array for Pulsars (LEAP) is a European Pulsar Timing Array project that combines the Lovell, Effelsberg, Nan\c{c}ay, Sardinia, and Westerbork radio telescopes into a single tied-array, and makes monthly observations of a set of millisecond pulsars (MSPs). The overview of our experiment is presented in Bassa et al. (2016). Baseband data are recorded at a central frequency of 1396 MHz and a bandwidth of 128 MHz at each telescope, and are correlated offline on a cluster at Jodrell Bank Observatory using a purpose-built correlator, detailed in Smits et al. (2017). LEAP offers a substantial increase in sensitivity over that of the individual telescopes, and can operate in timing and imaging modes (notably in observations of the galactic centre radio magnetar; Wucknitz 2015). To date, 4 years of observations have been reduced. Here, we report on the scientific projects which have made use of LEAP data.

Polarised radio emission from PSR J1745-2900 has already been used to investigate the strength of the magnetic field in the Galactic Centre, close to Sagittarius A*. Here we report how persistent radio emission from this magnetar, for over four years since its discovery, has revealed large changes in the observed Faraday rotation measure, by up to 3500 rad m$^{-2}$ (a five per cent fractional change). From simultaneous analysis of the dispersion measure, we determine that these fluctuations are dominated by variations in the projected magnetic field, rather than the integrated free electron density, along the changing line of sight to the rapidly moving magnetar. From a structure function analysis of rotation measure variations, and a recent epoch of rapid change of rotation measure, we determine a minimum scale of magnetic fluctuations of size ~ 2 au at the Galactic Centre distance, inferring PSR J1745-2900 is just ~ 0.1 pc behind an additional scattering screen.

Coronal magnetic flux ropes are closely related to large-scale solar activities. Using a 2.5-dimensional time-dependent ideal magnetohydrodynamic (MHD) model in Cartesian coordinates, we carry out numerical simulations to investigate the evolution of a magnetic system consisting of a flux rope embedded in a fully-closed quadrupolar magnetic field with different photospheric flux distributions. It is found that, when the photospheric flux is not concentrated too much towards the polarity inversion line (PIL) and the constraint exerted by the background field is not too weak, the equilibrium states of the system are divided into two branches: the rope sticks to the photosphere for the lower branch and levitates in the corona for the upper branch. These two branches are connected by an upward catastrophe (from the lower branch to the upper) and a downward catastrophe (from the upper branch to the lower). Our simulations reveal that there exist both upward and downward catastrophes in quadrupolar fields, which may be either force-free or non-force-free. The existence and the properties of these two catastrophes are influenced by the photospheric flux distribution, and a downward catastrophe is always paired with an upward catastrophe. Comparing the decay indices in catastrophic and non-catastrophic cases, we infer that torus unstable might be a necessary but not sufficient condition for a catastrophic system.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

In this paper, we identified the magnetic source locations of 142 quasi- homologous (QH) coronal mass ejections (CMEs), of which 121 are from solar cycle (SC) 23 and 21 from SC 24. Among those CMEs, 63% originated from the same source location as their predecessor (defined as S-type), while 37% origi- nated from a different location within the same active region as their predecessor (defined as D-type). Their distinctly different waiting time distribution, peaking around 7.5 and 1.5 hours for S- and D-type CMEs, suggests that they might involve different physical mechanisms with different characteristic time scales. Through detailed analysis based on non-linear force free (NLFF) coronal mag- netic field modeling of two exemplary cases, we propose that the S-type QH CMES might involve a recurring energy release process from the same source location (by magnetic free energy replenishment), whereas the D-type QH CMEs can happen when a flux tube system disturbed by a nearby CME.

We studied the background field for 60 two-ribbon flares of M-and-above classes during 2011--2015. These flares are categorized into two groups, i.e., \emph{eruptive} and \emph{confined} flares, based on whether a flare is associated with a coronal mass ejection or not. The background field of source active regions is approximated by a potential field extrapolated from the $B_z$ component of vector magnetograms provided by the Helioseismic and Magnetic Imager. We calculated the decay index $n$ of the background field above the flaring polarity inversion line, and defined a critical height $h_\mathrm{crit}$ corresponding to the theoretical threshold ($n_\mathrm{crit}=1.5$) of the torus instability. We found that $h_\mathrm{crit}$ is approximately half of the distance between the centroids of opposite polarities in active regions, and that the distribution of $h_\mathrm{crit}$ is bimodal: it is significantly higher for confined flares than for eruptive ones. The decay index increases monotonously with increasing height for 86\% (84\%) of the eruptive (confined) flares but displays a saddle-like profile for the rest 14\% (16\%), which are found exclusively in active regions of multipolar field configuration. Moreover, $n$ at the saddle bottom is significantly smaller in confined flares than that in eruptive ones. These results highlight the critical role of background field in regulating the eruptive behavior of two-ribbon flares.

High precision and long-lasting \emph{Kepler} data enabled us to estimate stellar properties with asteroseismology as an accurate tool. We performed asteroseismic analysis on six solar-like stars observed by the \emph{Kepler} mission: KIC 6064910, KIC 6766513, KIC 7107778, KIC 10079226, KIC 10147635 and KIC 12069127. The extraction of seismic information includes two parts. First, we obtained two global asteroseismic parameters, mean large separation $\Delta\nu$ and frequency of maximum power $\nu_{\rm{max}}$, with autocorrelation function and collapsed autocorrelation function. Second, we extracted individual oscillation modes $\nu_{nl}$ with low-$l$ degree using a least-squares fit. Stellar grid models were built with Yale Rotating Stellar Evolution Code (YREC) to analyze stellar properties. They covered the range of $M=0.8\sim1.8$ $M_{\odot}$ with a step of 0.02 $M_{\odot}$ and $\rm{[Fe/H]}=-0.3\sim0.4$ $\rm{dex}$ with a step of 0.1 $\rm{dex}$. We used a Bayesian approach to estimate stellar fundamental parameters of the six stars, under the constraints of asteroseismic parameters ($\Delta\nu$, $\nu_{\rm{max}}$) and non-asteroseismic parameters ($T_{\rm{eff}}$, $\rm{[Fe/H]}$). We discover that the six targets include five sub-giant stars with $1.2\sim1.5$ $M_{\odot}$ and one main-sequence star with $1.08M_{\odot}$, and with ages in the range of $3\sim5$ $\rm{Gyr}$.

The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tied-array beam from the raw voltages from each of telescopes. We explain the concepts and techniques involved in the process of adding the raw voltages coherently. We further present the software processing pipeline that is specifically designed to deal with data from widely spaced, inhomogeneous radio telescopes and describe the steps involved in preparing, correlating and creating the tied-array beam. This includes polarization calibration, bandpass correction, frequency dependent phase correction, interference mitigation and pulsar gating. A link is provided where the software can be obtained.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

We intent to use stellar models with and without alpha-enhancement, as well as asteroseismic analysis, to study two alpha-enhanced stars, KIC 7976303 and KIC 8694723. For the alpha-enhanced models, we adopt [alpha/Fe] = 0.2, and 0.4. For the evolved sub-giant KIC 7976303 with mixed-modes, alpha-enhanced models fit the observations better than the models without alpha-enhancement, and point to a star with M = 1.20 $\pm$ 0.02 M_odot, t = 4.88 $\pm$ 0.08 Gyr, R = 2.04 $\pm$ 0.01 R_odot, and L = 5.07 $\pm$ 0.05 L_odot. For the post turnoff star KIC 8694723, we find that the models fit the observations well in all three cases ([alpha/Fe] = 0.0, 0.2, 0.4). According to the observed metal abundances, only alpha- enhanced models are used to estimate the stellar parameters, which are M = 1.13 $\pm$ 0.06 M_odot, t = 5.47 $\pm$ 0.43 Gyr, R = 1.55 $\pm$ 0.13 R_odot, and L = 3.18 $\pm$ 0.26 L_odot. Our alpha-enhanced models indicate KIC 7976303 has a larger mass and younger age than in previous works based on standard models. Furthermore, the differences of estimated mass and age between the three cases are $\sim$ 0.1 M_odot and $\sim$ 0.5 - 1.3 Gyr. These results suggest that we include alpha-enhancement in the modeling of alpha-enhanced stars, such as members of GCs (Globular Clusters), metal-poor stars in the disc and in the halo.

Oxygen and carbon are important elements in stellar populations. Their behavior refers to the formation history of the stellar populations. C and O abundances would also obviously influence stellar opacities and the overall metal abundance $Z$. With observed high-quality spectroscopic properties, we construct stellar models with C and O elements, to give more accurate ages for 70 metal-poor dwarfs, which have been determined to be high-$\alpha$ halo, low-$\alpha$ halo and thick-disk stars. Our results show that high-$\alpha$ halo stars are relatively older than low-$\alpha$ halo stars by around 2.0 Gyr. The thick-disk population has an age range between the two halo populations. The age distribution profiles indicate that high-$\alpha$ halo and low-$\alpha$ halo stars match the in situ accretion simulation by Zolotov et al., and the thick-disk stars might be formed in a relatively quiescent and long-lasting process. We also note that stellar ages are very sensitive to O abundance, since the ages clearly increase with increasing [O/Fe] values. Additionally, we obtain several stars with peculiar ages, including 2 young thick-disk stars and 12 stars older than the universe age.

Typical solar flares display two quasi-parallel, bright ribbons on the chromosphere. In between is the polarity inversion line (PIL) separating concentrated magnetic fluxes of opposite polarity in active regions (ARs). Intriguingly a series of flares exhibiting X-shaped ribbons occurred at the similar location on the outskirts of NOAA AR 11967, where magnetic fluxes were scattered, yet three of them were alarmingly energetic. The X shape, whose center coincided with hard X-ray emission, was similar in UV/EUV, which cannot be accommodated in the standard flare model. Mapping out magnetic connectivities in potential fields, we found that the X morphology was dictated by the intersection of two quasi-separatrix layers, i.e., a hyperbolic flux tube (HFT), within which a separator connecting a double null was embedded. This topology was not purely local but regulated by fluxes and flows over the whole AR. The nonlinear force-free field model suggested the formation of a current layer at the HFT, where the current dissipation can be mapped to the X-shaped ribbons via field-aligned heat conduction. These results highlight the critical role of HFTs in 3D magnetic reconnection and have important implications for astrophysical and laboratory plasmas.

Single pulses preserve information about the pulsar radio emission and propagation in the pulsar magnetosphere, and understanding the behaviour of their variability is essential for estimating the fundamental limit on the achievable pulsar timing precision. Here we report the findings of our analysis of single pulses from PSR J1713+0747 with data collected by the Large European Array for Pulsars (LEAP). We present statistical studies of the pulse properties that include distributions of their energy, phase and width. Two modes of systematic sub-pulse drifting have been detected, with a periodicity of 7 and 3 pulse periods. The two modes appear at different ranges of pulse longitude but overlap under the main peak of the integrated profile. No evidence for pulse micro-structure is seen with a time resolution down to 140 ns. In addition, we show that the fractional polarisation of single pulses increases with their pulse peak flux density. By mapping the probability density of linear polarisation position angle with pulse longitude, we reveal the existence of two orthogonal polarisation modes. Finally, we find that the resulting phase jitter of integrated profiles caused by single pulse variability can be described by a Gaussian probability distribution only when at least 100 pulses are used for integration. Pulses of different flux densities and widths contribute approximately equally to the phase jitter, and no improvement on timing precision is achieved by using a sub-set of pulses with a specific range of flux density or width.

This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.

A background of nanohertz gravitational waves from supermassive black hole binaries could soon be detected by pulsar timing arrays, which measure the times-of-arrival of radio pulses from millisecond pulsars with very high precision. The European Pulsar Timing Array uses five large European radio telescopes to monitor high-precision millisecond pulsars, imposing in this way strong constraints on a gravitational wave background. To achieve the necessary precision needed to detect gravitational waves, the Large European Array for Pulsars (LEAP) performs simultaneous observations of pulsars with all five telescopes, which allows us to coherently add the radio pulses, maximize the signal-to-noise of pulsar signals and increase the precision of times-of-arrival. We report on the progress made and results obtained by the LEAP collaboration, and in particular on the addition of the Sardinia Radio Telescope to the LEAP observations during its scientific validation phase. In addition, we discuss how LEAP can be used to monitor strong-gravity systems such as double neutron star systems and impose strong constraints on post-keplerian parameters.

We investigate the evolutions of two prominences (P1,P2) and two bundles of coronal loops (L1,L2), observed with SDO/AIA near the east solar limb on 2012 September 22. It is found that there were large-amplitude oscillations in P1 and L1, but no detectable motions in P2 and L2. These transverse oscillations were triggered by a large-scale coronal wave, originating from a large flare in a remote active region behind the solar limb. By carefully comparing the locations and heights of these oscillating and non-oscillating structures, we conclude that the propagating height of the wave is between 50 Mm and 130 Mm. The wave energy deposited in the oscillating prominence and coronal loops is at least of the order of $10^{28}$ erg. Furthermore, local magnetic field strength and Alfv\'{e}n speeds are derived from the oscillating periods and damping time scales, which are extracted from the time series of the oscillations. It is demonstrated that oscillations can be used in not only coronal seismology, but also revealing the properties of the wave.

The Double Pulsar system PSR J0737-3039A/B has proven to be an excellent laboratory for high precision tests of general relativity. With additional years of timing measurements and new telescopes like the Square Kilometre Array (SKA), the precision of these tests will increase and new effects like the Lense-Thirring precession of the orbit will become measurable. Here, we discuss the prospects of measuring the Lense-Thirring effect and thereby constraining the equations of state at supra-nuclear densities in neutron stars using the Double Pulsar.

The 2011 January 28 M1.4 flare exhibits two side-by-side candle-flame-shaped flare loop systems underneath a larger cusp-shaped structure during the decay phase, as observed at the northwestern solar limb by the Solar Dynamics Observatory (SDO). The northern loop system brightens following the initiation of the flare within the southern loop system, but all three cusp-shaped structures are characterized by ~ 10 MK temperatures, hotter than the arch-shaped loops underneath. The "Ahead" satellite of the Solar Terrestrial Relations Observatory (STEREO) provides a top view, in which the post-flare loops brighten sequentially, with one end fixed while the other apparently slipping eastward. By performing stereoscopic reconstruction of the post-flare loops in EUV and mapping out magnetic connectivities, we found that the footpoints of the post-flare loops are slipping along the footprint of a hyperbolic flux tube (HFT) separating the two loop systems, and that the reconstructed loops share similarity with the magnetic field lines that are traced starting from the same HFT footprint, where the field lines are relatively flexible. These results argue strongly in favor of slipping magnetic reconnection at the HFT. The slipping reconnection was likely triggered by the flare and manifested as propagative dimmings before the loop slippage is observed. It may contribute to the late-phase peak in Fe XVI 33.5 nm, which is even higher than its main-phase counterpart, and may also play a role in the density and temperature asymmetry observed in the northern loop system through heat conduction.

We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.

Since the release of LAMOST (the Large Sky Area Multi-Object Fiber Spectroscopic Telescope) catalog, we have the opportunity to use the LAMOST DR2 stellar catalog and \emph{WISE All-sky catalog} to search for 22 $\mu$m excess candidates. In this paper, we present 10 FGK candidates which show an excess in the infrared (IR) at 22 $\mu$m. The ten sources are all the newly identified 22 $\mu$m excess candidates. Of these 10 stars, 5 stars are F type and 5 stars are G type. The criterion for selecting candidates is $K_s-[22]_{\mu m}\geq0.387$. In addition, we present the spectral energy distributions (SEDs) covering wavelength from optical to mid-infrared band. Most of them show an obvious excess from 12 $\mu$m band and three candidates even show excess from 3.4 $\mu$m. To characterize the amount of dust, we also estimate the fractional luminosity of ten 22 $\mu$m excess candidates.

We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.

The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.

SDO/EVE provides rich information of the thermodynamic processes of solar activities, particularly of solar flares. Here, we develop a method to construct thermodynamic spectrum (TDS) charts based on the EVE spectral lines. This tool could be potentially useful to the EUV astronomy to learn the eruptive activities on the distant astronomical objects. Through several cases, we illustrate what we can learn from the TDS charts. Furthermore, we apply the TDS method to 74 flares equal to or greater than M5.0-class, and reach the following statistical results. First, EUV peaks are always behind the soft X-ray (SXR) peaks and stronger flares tend to have a faster cooling rate. There is a power law correlation between the peak delay times and the cooling rates, suggesting a coherent cooling process of flares from SXR to EUV emissions. Second, there are two distinct temperature drift patterns, called Type I and Type II. For Type I flares, the enhanced emission drifts from high to low temperature like a quadrilateral, whereas for Type II flares, the drift pattern looks like a triangle. Statistical analysis suggests that Type II flares are more impulsive than Type I flares. Third, for late-phase flares, the peak intensity ratio of the late phase to the main phase is roughly correlated with the flare class, and the flares with a strong late phase are all confined. We believe that the re-deposition of the energy carried by a flux rope, that unsuccessfully erupts out, into thermal emissions is responsible for the strong late phase found in a confined flare. Besides, with some cases we illustrate the signatures of the flare thermodynamic process in the chromosphere and transition region in TDS charts. These results provide new clues to advance our understanding of the thermodynamic processes of solar flares and associated solar eruptions, e.g., coronal mass ejections.

The Large European Array for Pulsars (LEAP) is an experiment that harvests the collective power of Europe's largest radio telescopes in order to increase the sensitivity of high-precision pulsar timing. As part of the ongoing effort of the European Pulsar Timing Array (EPTA), LEAP aims to go beyond the sensitivity threshold needed to deliver the first direct detection of gravitational waves. The five telescopes presently included in LEAP are: the Effelsberg telescope, the Lovell telescope at Jodrell Bank, the Nan\c cay radio telescope, the Sardinia Radio Telescope and the Westerbork Synthesis Radio Telescope. Dual polarization, Nyquist-sampled time-series of the incoming radio waves are recorded and processed offline to form the coherent sum, resulting in a tied-array telescope with an effective aperture equivalent to a 195-m diameter circular dish. All observations are performed using a bandwidth of 128 MHz centered at a frequency of 1396 MHz. In this paper, we present the design of the LEAP experiment, the instrumentation, the storage and transfer of data, and the processing hardware and software. In particular, we present the software pipeline that was designed to process the Nyquist-sampled time-series, measure the phase and time delays between each individual telescope and a reference telescope and apply these delays to form the tied-array coherent addition. The pipeline includes polarization calibration and interference mitigation. We also present the first results from LEAP and demonstrate the resulting increase in sensitivity, which leads to an improvement in the pulse arrival times.

We present the multi-point and multi-wavelength observation and analysis on a solar coronal jet and coronal mass ejection (CME) event in this paper. Employing the GCS model, we obtained the real (three-dimensional) heliocentric distance and direction of the CME and found it propagate in a high speed over 1000 km/s . The jet erupted before and shared the same source region with the CME. The temporal and spacial relation- ship between them guide us the possibility that the jet triggered the CME and became its core. This scenario could promisingly enrich our understanding on the triggering mechanism of coronal mass ejections and their relations with coronal large-scale jets. On the other hand, the magnetic field configuration of the source region observed by the SDO/HMI instrument and the off- limb inverse Y-shaped configuration observed by SDO/AIA 171 A passband, together provide the first detailed observation on the three-dimensional reconnection process of large-scale jets as simulated in Pariat et al. 2009. The erupting process of the jet highlights that filament-like materials are important during the eruption not only of small-scale X-ray jets (Sterling et al. 2015) but also probably of large-scale EUV jets. Based on our observation and analysis, we propose a most possible mechanism for the whole event with a blob structure overlaying the three-dimensional structure of the jet to describe the interaction between the jet and the CME.

We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.

We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $\Omega_\mathrm{gw}(f)h^2 < 1.1\times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $G\mu/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $G\mu/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $\Omega^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.

Two-ribbon brightenings are one of the most remarkable characteristics of an eruptive solar flare and are often used for predicting the occurrence of coronal mass ejections (CMEs). Nevertheless, it was called in question recently whether all two-ribbon flares are eruptive. In this paper, we investigate a two ribbon-like white-light (WL) flare that is associated with a failed magnetic flux rope (MFR) eruption on 2015 January 13, which has no accompanying CME in the WL coronagraph. Observations by \textit{Optical and Near-infrared Solar Eruption Tracer} and \textit{Solar Dynamics Observatory} reveal that, with the increase of the flare emission and the acceleration of the unsuccessfully erupting MFR, two isolated kernels appear at the WL 3600 {\AA} passband and quickly develop into two elongated ribbon-like structures. The evolution of the WL continuum enhancement is completely coincident in time with the variation of \textit{Fermi} hard X-ray 26--50 keV flux. Increase of continuum emission is also clearly visible at the whole FUV and NUV passbands observed by \textit{Interface Region Imaging Spectrograph}. Moreover, in one WL kernel, the \ion{Si}{4}, \ion{C}{2}, and \ion{Mg}{2} h/k lines display significant enhancement and non-thermal broadening. However, their Doppler velocity pattern is location-dependent. At the strongly bright pixels, these lines exhibit a blueshift; while at moderately bright ones, the lines are generally redshifted. These results show that the failed MFR eruption is also able to produce a two-ribbon flare and high-energy electrons that heat the lower atmosphere, causing the enhancement of the WL and FUV/NUV continuum emissions and chromospheric evaporation.

The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.

Solar flares typically have an impulsive phase that followed by a gradual phase as best seen in soft X-ray emissions. A recent discovery based on the EUV Variability Experiment (EVE) observations onboard the Solar Dynamics Observatory (SDO) reveals that some flares exhibit a second large peak separated from the first main phase peak by tens of minutes to hours, which is coined as the flare's EUV late phase. In this paper, we address the origin of the EUV late phase by analyzing in detail two late phase flares, an M2.9 flare on 2010 October 16 and an M1.4 flare on 2011 February 18, using multi-passband imaging observations from the Atmospheric Imaing Assembly (AIA) onboard SDO. We find that: (1) the late phase emission originates from a different magnetic loop system, which is much larger and higher than the main phase loop system. (2) The two loop systems have different thermal evolution. While the late phase loop arcade reaches its peak brightness progressively at a later time spanning for more than one hour from high to low temperatures, the main phase loop arcade reaches its peak brightness at almost the same time (within several minutes) in all temperatures. (3) Nevertheless, the two loop systems seem to be connected magnetically, forming an asymmetric magnetic quadruple configuration. (4) Further, the footpoint brightenings in UV wavelengths show a systematic delay of about one minute from the main flare region to the remote footpoint of the late phase arcade system. We argue that the EUV late phase is the result of a long-lasting cooling process in the larger magnetic arcade system.

The second peak in the Fe XVI 33.5 nm line irradiance observed during solar flares by Extreme ultraviolet Variability Experiment (EVE) is known as Extreme UltraViolet (EUV) late phase. Our previous paper (Liu et al. 2013) found that the main emissions in the late phase are originated from large-scale loop arcades that are closely connected to but different from the post flare loops (PFLs), and we also proposed that a long cooling process without additional heating could explain the late phase. In this paper, we define the extremely large late phase because it not only has a bigger peak in the warm 33.5 irradiance profile, but also releases more EUV radiative energy than the main phase. Through detailedly inspecting the EUV images from three point-of-view, it is found that, besides the later phase loop arcades, the more contribution of the extremely large late phase is from a hot structure that fails to erupt. This hot structure is identified as a flux rope, which is quickly energized by the flare reconnection and later on continuously produces the thermal energy during the gradual phase. Together with the late-phase loop arcades, the fail to erupt flux rope with the additional heating create the extremely large EUV late phase.

Millisecond pulsars (MSPs) are known as highly stable celestial clocks. Nevertheless, recent studies have revealed the unstable nature of their integrated pulse profiles, which may limit the achievable pulsar timing precision. In this paper, we present a case study on the pulse profile variability of PSR J1022+1001. We have detected approximately 14,000 sub-pulses (components of single pulses) in 35-hr long observations, mostly located at the trailing component of the integrated profile. Their flux densities and fractional polarisation suggest that they represent the bright end of the energy distribution in ordinary emission mode and are not giant pulses. The occurrence of sub-pulses from the leading and trailing components of the integrated profile is shown to be correlated. For sub-pulses from the latter, a preferred pulse width of approximately 0.25 ms has been found. Using simultaneous observations from the Effelsberg 100-m telescope and the Westerbork Synthesis Radio Telescope, we have found that the integrated profile varies on a timescale of a few tens of minutes. We show that improper polarisation calibration and diffractive scintillation cannot be the sole reason for the observed instability. In addition, we demonstrate that timing residuals generated from averages of the detected sub-pulses are dominated by phase jitter, and place an upper limit of ~700 ns for jitter noise based on continuous 1-min integrations.

The discovery and timing of radio pulsars within the Galactic centre is a fundamental aspect of the SKA Science Case, responding to the topic of "Strong Field Tests of Gravity with Pulsars and Black Holes" (Kramer et al. 2004; Cordes et al. 2004). Pulsars have in many ways proven to be excellent tools for testing the General theory of Relativity and alternative gravity theories (see Wex (2014) for a recent review). Timing a pulsar in orbit around a companion, provides a unique way of probing the relativistic dynamics and spacetime of such a system. The strictest tests of gravity, in strong field conditions, are expected to come from a pulsar orbiting a black hole. In this sense, a pulsar in a close orbit ($P_{\rm orb}$ < 1 yr) around our nearest supermassive black hole candidate, Sagittarius A* - at a distance of ~8.3 kpc in the Galactic centre (Gillessen et al. 2009a) - would be the ideal tool. Given the size of the orbit and the relativistic effects associated with it, even a slowly spinning pulsar would allow the black hole spacetime to be explored in great detail (Liu et al. 2012). For example, measurement of the frame dragging caused by the rotation of the supermassive black hole, would allow a test of the "cosmic censorship conjecture." The "no-hair theorem" can be tested by measuring the quadrupole moment of the black hole. These are two of the prime examples for the fundamental studies of gravity one could do with a pulsar around Sagittarius A*. As will be shown here, SKA1-MID and ultimately the SKA will provide the opportunity to begin to find and time the pulsars in this extreme environment.

The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.

The anticipated discovery of a pulsar in orbit with a black hole is expected to provide a unique laboratory for black hole physics and gravity. In this context, the next generation of radio telescopes, like the Five-hundred-metre Aperture Spherical radio Telescope (FAST) and the Square Kilometre Array (SKA), with their unprecedented sensitivity, will play a key role. In this paper, we investigate the capability of future radio telescopes to probe the spacetime of a black hole and test gravity theories, by timing a pulsar orbiting a stellar-mass-black-hole (SBH). Based on mock data simulations, we show that a few years of timing observations of a sufficiently compact pulsar-SBH (PSR-SBH) system with future radio telescopes would allow precise measurements of the black hole mass and spin. A measurement precision of one per cent can be expected for the spin. Measuring the quadrupole moment of the black hole, needed to test GR's no-hair theorem, requires extreme system configurations with compact orbits and a large SBH mass. Additionally, we show that a PSR-SBH system can lead to greatly improved constraints on alternative gravity theories even if they predict black holes (practically) identical to GR's. This is demonstrated for a specific class of scalar-tensor theories. Finally, we investigate the requirements for searching for PSR-SBH systems. It is shown that the high sensitivity of the next generation of radio telescopes is key for discovering compact PSR-SBH systems, as it will allow for sufficiently short survey integration times.

In recent years, instrumentation enabling pulsar observations with unprecedentedly high fractional bandwidth has been under development which can be used to substantially improve the precision of pulsar timing experiments. The traditional template-matching method used to calculate pulse times-of-arrival (ToAs), may not function effectively on these broadband data due to a variety of effects such as diffractive scintillation in the interstellar medium, profile variation as a function of frequency, dispersion measure (DM) evolution and so forth. In this paper, we describe the channelised Discrete Fourier Transform method that can greatly mitigate the influence of the aforementioned effects when measuring ToAs from broadband timing data. The method is tested on simulated data, and its potential in improving timing precision is shown. We further apply the method to PSR J1909$-$3744 data collected at the Nan\c{c}ay Radio Telescope with the Nan\c{c}ay Ultimate Pulsar Processing Instrument. We demonstrate a removal of systematics due to the scintillation effect as well as improvement on ToA measurement uncertainties. Our method also determines temporal variations in dispersion measure, which are consistent with multi-channel timing approaches used earlier.

The aim of this paper is to determinate the fundamental parameters of six exoplanet host (EH) stars and their planets. While techniques for detecting exoplanets yield properties of the planet only as a function of the properties of the host star, hence, we must accurately determine parameters of EH stars at first. For this reason, we constructed a grid of stellar models including diffusion and rotation-induced extra-mixing with given ranges of input parameters (i.e. mass, metallicity, and initial rotation rate). In addition to the commonly used observational constraints such as the effective temperature T_{eff}, luminosity L and metallicity [Fe/H], we added two observational constraints, the lithium abundance log N (Li) and the rotational period P_{rot}. These two additional observed parameters can make further constrains on the model due to their correlations with mass, age and other stellar properties. Hence, our estimations of fundamental parameters for these EH stars and their planets are with higher precision than previous works. Therefore, the combination of rotational period and lithium help us to obtain more accurate parameters for stars, leading to an improvement of the knowledge of the physical state about the EH stars and their planets.

We present multi-wavelength observations of a prominence eruption originating from a quadrupolar field configuration, in which the prominence was embedded in a side-arcade. Within the two-day period prior to its eruption on 2012 October 22, the prominence was perturbed three times by chromospheric fibrils underneath, which rose upward, became brightened, and merged into the prominence, resulting in horizontal flows along the prominence axis, suggesting that the fluxes carried by the fibrils were incorporated into the magnetic field of the prominence. These perturbations caused the prominence to oscillate and to rise faster than before. The absence of intense heating within the first two hours after the onset of the prominence eruption, which followed an exponential increase in height, indicates that ideal instability played a crucial role. The eruption involved interactions with the other side-arcade, leading up to a twin coronal mass ejection, which was accompanied by transient surface brightenings in the central arcade, followed by transient dimmings and brightenings in the two side-arcades. We suggest that flux feeding from chromospheric fibrils might be an important mechanism to trigger coronal eruptions.

We report the observation of an X-class long-duration flare which is clearly confined. It appears as a compact-loop flare in the traditional EUV passbands (171 and 195 {\AA}), but in the passbands sensitive to flare plasmas (94 and 131 {\AA}), it exhibits a cusp-shaped structure above an arcade of loops like other long-duration events. Inspecting images in a running difference approach, we find that the seemingly diffuse, quasi-static cusp-shaped structure consists of multiple nested loops that repeatedly rise upward and disappear approaching the cusp edge. Over the gradual phase, we detect numerous episodes of loop rising, each lasting minutes. A differential emission measure analysis reveals that the temperature is highest at the top of the arcade and becomes cooler at higher altitudes within the cusp-shaped structure, contrary to typical long-duration flares. With a nonlinear force-free model, our analysis shows that the event mainly involves two adjacent sheared arcades separated by a T-type hyperbolic flux tube (HFT). One of the arcades harbors a magnetic flux rope, which is identified with a filament that survives the flare owing to the strong confining field. We conclude that a new emergence of magnetic flux in the other arcade triggers the flare, while the preexisting HFT and flux rope dictate the structure and dynamics of the flare loops and ribbons during the long-lasting decay phase, and that a quasi-separatrix layer high above the HFT could account for the cusp-shaped structure.

In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists of two major steps. First, the spectral index and amplitude of dispersion measure variations are measured via a time-domain spectral analysis. Second, the linear optimal filter is constructed based on the model parameters found in the first step, and is used to extract the dispersion measure variation waveforms. Compared to current existing methods, this method has better time resolution for the study of short timescale dispersion variations, and generally produces smaller errors in waveform estimations. This method can process irregularly sampled data without any interpolation because of its time-domain nature. Furthermore, it offers the possibility to interpolate or extrapolate the waveform estimation to regions where no data is available. Examples using simulated data sets are included for demonstration.

The oscillations of the solar-like star HD 49933 have been observed thoroughly by CoRot. Two dozens of frequency shifts, which are closely related with the change in magnetic activity, have been measured. To explore the effects of the magnetic activity on the frequency shifts, we calculate frequency shifts for the radial and $l = 1$ p-modes of HD 49933 with the general variational method, which evaluates the shifts using a spatial integral of the product of a kernel and some sources. The theoretical frequency shifts well reproduce the observation. The magnitudes and positions of the sources are determined according to the $\chi^2$ criterion. We predict the source that contributes to both $l = 0$ and $l = 1$ modes is located at $0.48 - 0.62$Mm below the stellar surface. In addition, based on the assumption that $A_{0}$ is proportional to the change in the MgII activity index $\Delta{i}_{MgII}$, we obtained that the change of MgII index between minimum and maximum of HD 49933 cycle period is about 0.665. The magnitude of the frequency shifts compared to the Sun already told us that HD 49933 is much more active than the Sun, which is further confirmed in this paper. Furthermore, our calculation on the frequency shifts of $l = 1$ modes indicates the variation of turbulent velocity in the stellar convective zone may be an important source for the $l = 1$ shifts.

A jet is a considerable amount of plasma being ejected from the chromosphere or lower corona into the higher corona and is a common phenomenon. Usually, a jet is triggered by a brightening or a flare, which provides the first driving force to push plasma upward. In this process, magnetic reconnection is thought to be the mechanism to convert magnetic energy into thermal, nonthermal, and kinetic energies. However, most jets could reach an unusual high altitude and end much later than the end of its associated flare. This fact implies that there is another way to continuously transfer magnetic energy into kinetic energy even after the reconnection. The picture described above is well known in the community, but how and how much magnetic energy is released through a way other than reconnection is still unclear. By studying a prominence-like jet observed by SDO/AIA and STEREO-A/EUVI, we find that the continuous relaxation of the post-reconnection magnetic field structure is an important process for a jet to climb up higher than it could through only reconnection. The kinetic energy of the jet gained through the relaxation is 1.6 times that gained from the reconnection. The resultant energy flux is hundreds of times larger than the flux required for the local coronal heating, suggesting that such jets are a possible source to keep the corona hot. Furthermore, rotational motions appear all the time during the jet. Our analysis suggests that torsional Alfv\'en waves induced during reconnection could not be the only mechanism to release magnetic energy and drive jets.

The frequency ratios and of HD 49933 exhibit an increase at high frequencies. This behavior also exists in the ratios of other stars, which is considered to result from the low signal-to-noise ratio and the larger line width at the high-frequency end and could not be predicted by stellar models in previous work. Our calculations show that the behavior not only can be reproduced by stellar models, but can be predicted by asymptotic formulas of the ratios. The frequency ratios of the Sun, too, can be reproduced well by the asymptotic formulas. The increased behavior derives from the fact that the gradient of mean molecular weight at the bottom of the radiative region hinders the propagation of p-modes, while the hindrance does not exist in the convective core. This behavior should exist in the ratios of stars with a large convective core. The characteristic of the ratios at high frequencies provides a strict constraint on stellar models and aids in determining the size of the convective core and the extent of overshooting. Observational constraints point to a star with $M=1.28\pm0.01 M_{\odot}$, $R=1.458\pm0.005 R_{\odot}$, $t=1.83\pm0.1$ Gyr, $r_{cc}=0.16\pm0.02 R_{\odot}$, $\alpha=1.85\pm0.05$, and $\delta_{ov}=0.6\pm0.2$ for HD 49933.

We extend a previous statistical solar flare study of 155 GOES M- and X-class flares observed with AIA/SDO (Aschwanden 2012) to all 7 coronal wavelengths (94, 131, 171, 193, 211, 304, 335 \ang) to test the wavelength-dependence of scaling laws and statistical distributions. Except for the 171 and 193 \ang\ wavelengths, which are affected by EUV dimming caused by coronal mass ejections (CMEs), we find near-identical size distributions of geometric (lengths $L$, flare areas $A$, volumes $V$, fractal dimension $D_2$), temporal (flare durations $T$), and spatio-temporal parameters (diffusion coefficient $\kappa$, spreading exponent $\beta$, and maximum expansion velocities $v_{max}$) in different wavelengths, which are consistent with the universal predictions of the fractal-diffusive avalanche model of a slowly-driven self-organized criticality (FD-SOC) system, i.e., $N(L) \propto L^{-3}$, $N(A) \propto A^{-2}$, $N(V) \propto V^{-5/3}$, $N(T) \propto T^{-2}$, $D_2=3/2$, for a Euclidean dimension $d=3$. Empirically we find also a new strong correlation $\kappa \propto L^{0.94\pm0.01}$ and the 3-parameter scaling law $L \propto \kappa\ T^{0.1}$, which is more consistent with the logistic-growth model than with classical diffusion. The findings suggest long-range correlation lengths in the FD-SOC system that operate in the vicinity of a critical state, which could be used for predictions of individual extreme events. We find also that eruptive flares (with accompanying CMEs), have larger volumes $V$, longer flare durations $T$, higher EUV and soft X-ray fluxes, and somewhat larger diffusion coefficients $\kappa$ than confined flares (without CMEs).

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

The discovery of a pulsar (PSR) in orbit around a black hole (BH) is expected to provide a superb new probe of relativistic gravity and BH properties. Apart from a precise mass measurement for the BH, one could expect a clean verification of the dragging of space-time caused by the BH spin. In order to measure the quadrupole moment of the BH for testing the no-hair theorem of general relativity (GR), one has to hope for a sufficiently massive BH. In this respect, a PSR orbiting the super-massive BH in the center of our Galaxy would be the ultimate laboratory for gravity tests with PSRs. But even for gravity theories that predict the same properties for BHs as GR, a PSR-BH system would constitute an excellent test system, due to the high grade of asymmetry in the strong field properties of these two components. Here we highlight some of the potential gravity tests that one could expect from different PSR-BH systems, utilizing present and future radio telescopes, like FAST and SKA.

Radio pulsars in relativistic binary systems are unique tools to study the curved space-time around massive compact objects. The discovery of a pulsar closely orbiting the super-massive black hole at the centre of our Galaxy, Sgr A*, would provide a superb test-bed for gravitational physics. To date, the absence of any radio pulsar discoveries within a few arc minutes of Sgr A* has been explained by one principal factor: extreme scattering of radio waves caused by inhomogeneities in the ionized component of the interstellar medium in the central 100 pc around Sgr A*. Scattering, which causes temporal broadening of pulses, can only be mitigated by observing at higher frequencies. Here we describe recent searches of the Galactic centre region performed at a frequency of 18.95 GHz with the Effelsberg radio telescope.

We calculated the populations of core-helium-burning (CHeB) stars and found that the secondary red clump (SRC) stars can form an SRC peak in the distributions of the frequency of maximum seismic amplitude ($\nu_{max}$) and mean large-frequency separation ($\Delta\nu$) of CHeB stars when metallicity $Z \geq$ 0.02. The $\nu_{max}$ and $\Delta\nu$ of CHeB stars are dependent not only on He core mass but on H-shell burning. The SRC peak is composed of the CHeB stars with mass roughly between the critical mass M_{Hef} and M_{Hef}+0.2 while He core mass is between about 0.33 and 0.36 M_{sun}. The location of the SRC peak can be affected by the mixing-length parameter $\alpha$, metallicity $Z$, and overshooting parameter $\delta_{ov}$. A decrease in $\alpha$ or increase in $Z$ or $\delta_{ov}$ leads to a movement of the SRC peak towards a lower frequency. However, the change in $Z$ and $\alpha$ only slightly affects the value of M_{Hef} but the variation in $\delta_{ov}$ can significantly affects the value of M_{Hef}. Thus the SRC peak might aid in determining the value of M_{Hef} and calibrating $\delta_{ov}$. In addition, the effects of convective acceleration of SRC stars and the $\nu_{max}$ of `semi-degenerate' stars decreasing with mass result in the appearance of a shoulder between about 40 and 50 $\mu$hz in the \dnu{} distribution. However, the convective acceleration of stars with M < M_{Hef} leads to the deficit in the $\nu_{max}$ distribution between about 9 and 20 $\mu$hz{}. Moreover, the value of the parameter $b$ of the relation between $\nu_{max}$ and $\Delta\nu$ for the populations with M > M_{Hef} is obviously larger than that for the populations with $M <$ \dmhef{}.

The discovery of radio pulsars in compact orbits around Sgr A* would allow an unprecedented and detailed investigation of the spacetime of the supermassive black hole. This paper shows that pulsar timing, including that of a single pulsar, has the potential to provide novel tests of general relativity, in particular its cosmic censorship conjecture and no-hair theorem for rotating black holes. These experiments can be performed by timing observations with 100 micro-second precision, achievable with the Square Kilometre Array for a normal pulsar at frequency above 15 GHz. Based on the standard pulsar timing technique, we develop a method that allows the determination of the mass, spin, and quadrupole moment of Sgr A*, and provides a consistent covariance analysis of the measurement errors. Furthermore, we test this method in detailed mock data simulations. It seems likely that only for orbital periods below ~0.3 yr is there the possibility of having negligible external perturbations. For such orbits we expect a ~10^-3 test of the frame dragging and a ~10^-2 test of the no-hair theorem within 5 years, if Sgr A* is spinning rapidly. Our method is also capable of identifying perturbations caused by distributed mass around Sgr A*, thus providing high confidence in these gravity tests. Our analysis is not affected by uncertainties in our knowledge of the distance to the Galactic center, R0. A combination of pulsar timing with the astrometric results of stellar orbits would greatly improve the measurement precision of R0.

Millisecond pulsars (MSPs) have been studied in detail since their discovery in 1982. The integrated pulse profiles of MSPs appear to be stable, which enables precision monitoring of the pulse times of arrival (TOAs). However, for individual pulses the shape and arrival phase can vary dramatically, which is known as pulse jitter. In this paper, we investigate the stability of integrated pulse profiles for 5 MSPs, and estimate the amount of jitter for PSR J0437-4715. We do not detect intrinsic profile shape variation based on integration times from ~10 to ~100 s with the provided instrumental sensitivity. For PSR J0437-4715 we calculate the jitter parameter to be f_J=0.067+-0.002, and demonstrate that the result is not significantly affected by instrumental TOA uncertainties. Jitter noise is also found to be independent of observing frequency and bandwidth around 1.4 GHz on frequency scales of <100 MHz, which supports the idea that pulses within narrow frequency scale are equally jittered. In addition, we point out that pulse jitter would limit TOA calculation for the timing observations with future telescopes like the Square Kilometre Array and the Five hundred metre Aperture Spherical Telescope. A quantitative understanding of pulse profile stability and the contribution of jitter would enable improved TOA calculations, which are essential for the ongoing endeavours in pulsar timing, such as the detection of the stochastic gravitational wave background.

At which height will a prominence inclined to be unstable, or where is the most probable critical height for the prominence destabilization? This question is statistically studied based on 362 solar limb prominences well-recognized by SLIPCAT from 2007 April to the end of 2009. We found that there are about 71% disrupted prominences (DPs), among which about 42% of them did not erupt successfully and about 89% of them experienced a sudden destabilization (SD) process. After a comprehensive analysis of the DPs, the following findings are discovered. (1) Most DPs become unstable at the height of 0.06 -- 0.14 R$_\odot$ from the solar surface, and there are two most probable critical heights, at which a prominence is much likely to get unstable; the primary one is 0.13 R$_\odot$ and the secondary one is 0.19 R$_\odot$. (2) There exists upper limit for the erupting velocity of eruptive prominences (EPs), which decreases following a power law with increasing height and mass; the kinetic energy of EPs accordingly has an upper limit too, which decreases as critical height increases. (3) Stable prominences (SPs) are generally longer and heavier than DPs, and not higher than 0.4 R$_\odot$. (4) About 62% of EPs were associated with CMEs; but there is no difference in apparent properties between EPs associated with and without CMEs.

Timing pulses of pulsars has proved to be a most powerful technique useful to a host of research areas in astronomy and physics. Importantly, the precision of this timing is not only affected by radiometer noise, but also by intrinsic pulse shape changes, interstellar medium (ISM) evolution, instrumental distortions, etc. In this paper we review the known causes of pulse shape variations and assess their effect on the precision and accuracy of a single measurement of pulse arrival time with current instrumentation. Throughout this analysis we use the brightest and most precisely timed millisecond pulsar (MSP), PSR J0437-4715, as a case study, and develop a set of diagnostic tools to evaluate profile stability in timing observations. We conclude that most causes of distortion can be either corrected by state-of-the-art techniques or taken into account in the estimation of time-of-arrival (TOA) uncertainties. The advent of a new generation of radio telescopes (e.g. the Square Kilometre Array, SKA), and their increase in collecting area has sparked speculation about the timing precision achievable through increases in gain. Based on our analysis of current data, we predict that for normal-brightness MSPs a TOA precision of between 80 and 230 ns can be achieved at 1.4 GHz with 10-minute integrations by the SKA. The actual rms timing residuals for each pulsar will be approximately at the same level only if all the other influences on timing precision (e.g. ISM, spin noise) are either corrected, modelled, or negligible.

Double or extended main-sequence turn-offs (DMSTOs) and dual red clump (RC) were observed in intermediate-age clusters, such as in NGC 1846 and 419. the DMSTOs are interpreted as that the cluster has two distinct stellar popula- tions with differences in age of about 200-300 Myr but with the same metallicity. The dual RC is interpreted as a result of a prolonged star formation. Using a stellar population-synthesis method, we calculated the evolutions of binary-star stellar population (BSP). We found that binary interactions and merging can reproduce the dual RC in the color-magnitude diagrams of an intermediate-age cluster, whereas in actuality only a single population exists. Moreover, the binary interactions can lead to an extended MSTO rather than DMSTOs. However, the rest of main sequence, subgiant branch and first giant branch are hardly spread by the binary interactions. Part of the observed dual RC and extended MSTO may be the results of binary interactions and merger.

We search for persistent and quasi-periodic release events of streamer blobs during 2007 with the Large Angle Spectrometric Coronagraph on the \textit{Solar and Heliospheric Observatory} and assess the velocity of the slow solar wind along the plasma sheet above the corresponding streamer by measuring the dynamic parameters of blobs. We find 10 quasi-periodic release events of streamer blobs lasting for three to four days. In each day of these events, we observe three-five blobs. The results are in line with previous studies using data observed near the last solar minimum. Using the measured blob velocity as a proxy for that of the mean flow, we suggest that the velocity of the background slow solar wind near the Sun can vary significantly within a few hours. This provides an observational manifestation of the large velocity variability of the slow solar wind near the Sun.

The pulsar PSR B1828-11 has long-term, highly periodic and correlated variations in both pulse shape and the rate of slow-down. This phenomenon may provide evidence for precession of the pulsar as suggested previously within the framework of free precession as well as forced one. On a presumption of forced precession, we propose a quark planet model to this precession henomenon instead, in which the pulsar is torqued by a quark planet. We construct this model by constraining mass of the pulsar ($M_{\rm psr}$), mass of the planet ($M_{\rm pl}$) and orbital radius of the planet ($r_{\rm pl}$). Five aspects are considered: derived relation between $M_{\rm psr}$ and $r_{\rm pl}$, movement of the pulsar around the center of mass, ratio of $M_{\rm psr}$ and $M_{\rm pl}$, gravitational wave radiation timescale of the planetary system, and death-line criterion. We also calculate the range of precession period derivative and gravitational wave strength (at earth) permitted by the model. Under reasonable parameters, the observed phenomenon can be understood by a pulsar ($10^{-4}\sim10^{-1}M_{\odot}$) with a quark planet ($10^{-8}\sim10^{-3}M_{\odot}$) orbiting it. According to the calculations presented, the pulsar would be a quark star because of its low mass, which might eject a lump of quark matter (to become a planet around) during its birth.