Flux and spectral variation characteristics of 3C 454.3 at the GeV band

We analyze the long-term lightcurve of 3C 454.3 observed with Fermi/LAT and investigate its relation to flux in the radio,optical and X-ray bands.By fitting the 1-day binned Ge V lightcurve with multiple Gaussian functions(MGF),we propose that the typical variability timescale in the Ge V band is 1–10 d.The Ge V flux variation is accompanied by the spectral variation characterized as fluxtracking,i.e.,"harder when brighter."The Ge V flux is correlated with the optical and X-ray fluxes,and a weak correlation betweenγ-ray flux and radio flux is also observed.Theγ-ray flux is not correlated with the optical linear polarization degree for the global lightcurves,but they show a correlation for the lightcurves before MJD 56000.The power density spectrum of the global lightcurve shows an obvious turnover at~7.7 d,which may indicate a typical variability timescale of 3C 454.3 in theγ-ray band.This is also consistent with the derived timescales by fitting the global lightcurve with MGF.The spectral evolution and an increase in the optical linear polarization degree along with the increase inγ-ray flux may indicate that the radiation particles are accelerated and the magnetic field is ordered by the shock processes during the outbursts.In addition,the nature of 3C 454.3 may be consistent with a self-organized criticality system,similar to Sagittariusand thus the outbursts could be from plasmoid ejections driven by magnetic reconnection.This may further support the idea that the jet radiation regions are magnetized.     

Testing a solar coronal magnetic field extrapolation code with the Titov–Dmoulin magnetic flux rope model

In the solar corona, the magnetic flux rope is believed to be a fundamental structure that accounts for magnetic free energy storage and solar eruptions. Up to the present, the extrapolation of the magnetic field from boundary data has been the primary way to obtain fully three-dimensional magnetic information about the corona. As a result, the ability to reliably recover the coronal magnetic flux rope is important for coronal field extrapolation. In this paper, our coronal field extrapolation code is examined with an analytical magnetic flux rope model proposed by Titov D′emoulin, which consists of a bipolar magnetic configuration holding a semi-circular line-tied flux rope in force-free equilibrium. By only using the vector field at the bottom boundary as input, we test our code with the model in a representative range of parameter space and find that the model field can be reconstructed with high accuracy. In particular, the magnetic topological interfaces formed between the flux rope and the surrounding arcade, i.e., the "hyperbolic flux tube" and "bald patch separatrix surface," are also reliably reproduced. By this test, we demonstrate that our CESE–MHD–NLFFF code can be applied to recovering the magnetic flux rope in the solar corona as long as the vector magnetogram satisfies the force-free constraints.     

Evolution of isolated G-band bright points: size,intensity and velocity

We study the evolution pattern of isolated G-band bright points(GBPs)in terms of their size, intensity and velocity. Using a high resolution image sequence taken with the Hinode/Solar Optical Telescope(SOT), we detect GBPs in each image by the Laplacian and Morphological Dilation algorithm, and track their evolutions by a 26-adjacent method in a three-dimensional space-time cube. For quantifying the evolution, we propose a quantification method based on lifetime normalization which aligns the different lifetimes to common stages. The quantification results show that, on average, the diameter of isolated GBPs changes from 166 to 173 km, then down to 165 km; the maximum intensity contrast changes from 1.012 to 1.027, then down to 1.011; however, the velocity changes from 1.709 to 1.593 km s-1, then up to 1.703 km s-1. The results indicate that the evolution follows a pattern such that the GBP is small, faint and fast-moving at the birth stage, becomes big, bright and slow-moving at the middle stage, then gets small, faint and fast-moving at the decay stage until disappearance. Although the differences are very small, a two-sample t-test is used to demonstrate there are significant differences in means between the distributions of the different stages. Furthermore, we quantify the relationship between the lifetimes of GBPs and their properties. It is found that there are positive correlations between the lifetimes and their sizes and intensities with correlation coefficients of0.83 and 0.65, respectively; however, there is a negative correlation between the lifetimes and velocities with a correlation coefficient of –0.49. In summary, the longer the GBP persists, the bigger, brighter and slower it will be.     

Evolution of isolated G-band bright points： size,intensity and velocity

We study the evolution pattern of isolated G-band bright points(GBPs)in terms of their size, intensity and velocity. Using a high resolution image sequence taken with the Hinode/Solar Optical Telescope(SOT), we detect GBPs in each image by the Laplacian and Morphological Dilation algorithm, and track their evolutions by a 26-adjacent method in a three-dimensional space-time cube. For quantifying the evolution, we propose a quantification method based on lifetime normalization which aligns the different lifetimes to common stages. The quantification results show that, on average, the diameter of isolated GBPs changes from 166 to 173 km, then down to 165 km; the maximum intensity contrast changes from 1.012 to 1.027, then down to 1.011; however, the velocity changes from 1.709 to 1.593 km s-1, then up to 1.703 km s-1. The results indicate that the evolution follows a pattern such that the GBP is small, faint and fast-moving at the birth stage, becomes big, bright and slow-moving at the middle stage, then gets small, faint and fast-moving at the decay stage until disappearance. Although the differences are very small, a two-sample t-test is used to demonstrate there are significant differences in means between the distributions of the different stages. Furthermore, we quantify the relationship between the lifetimes of GBPs and their properties. It is found that there are positive correlations between the lifetimes and their sizes and intensities with correlation coefficients of0.83 and 0.65, respectively; however, there is a negative correlation between the lifetimes and velocities with a correlation coefficient of –0.49. In summary, the longer the GBP persists, the bigger, brighter and slower it will be.     

Local sunspot oscillations and umbral dots

Data analysis of sunspot oscillations based on a 6-hr SDO run of an observation showed that low frequency(0.2 ω 1 m Hz) oscillations are locally similar to three and five minute oscillations.The oscillations in the sunspot are concentrated in cells of a few arcsec, each of which has its own oscillation spectrum. The analysis of two scenarios for sunspot oscillations leads to a conclusion that local sunspot oscillations occur due to a subphotospheric resonator for slow MHD waves. Empirical models of a sunspot atmosphere and the theory of slow waves in thin magnetic flux tubes are applied to modeling the subphotospheric resonator. The spectrum of local oscillations consists of a great number of lines. This kind of spectrum can occur only if the subphotospheric resonator is a magnetic tube with a rather weak magnetic field. Magnetic tubes of this sort are umbral dots that appear due to the convective tongues in monolithic sunspots. The interrelation of local oscillations with umbral dots and wavefronts of traveling waves in sunspots is discussed.     

Statistical improvement in detection level of gravitational microlensing events from their light curves

In astronomy,the brightness of a source is typically expressed in terms of magnitude.Conventionally,the magnitude is defined by the logarithm of received flux.This relationship is known as the Pogson formula.For received flux with a small signal to noise ratio(S/N),however,the formula gives a large magnitude error.We investigate whether the use of Inverse Hyperbolic Sine function(hereafter referred to as the Asinh magnitude)in the modified formulae could allow for an alternative calculation of magnitudes for small S/N flux,and whether the new approach is better for representing the brightness of that region.We study the possibility of increasing the detection level of gravitational microlensing using 40 selected microlensing light curves from the 2013 and 2014 seasons and by using the Asinh magnitude.Photometric data of the selected events are obtained from the Optical Gravitational Lensing Experiment(OGLE).We found that utilization of the Asinh magnitude makes the events brighter compared to using the logarithmic magnitude,with an average of about 3.42×10~(-2)magnitude and an average in the difference of error between the logarithmic and the Asinh magnitude of about 2.21×10(-2)magnitude.The microlensing events OB140847 and OB140885 are found to have the largest difference values among the selected events.Using a Gaussian fit to find the peak for OB140847 and OB140885,we conclude statistically that the Asinh magnitude gives better mean squared values of the regression and narrower residual histograms than the Pogson magnitude.Based on these results,we also attempt to propose a limit in magnitude value for which use of the Asinh magnitude is optimal with small S/N data.     

On the intrinsic shape of the gamma-ray spectrum for Fermi blazars

The curvature of the γ-ray spectrum in blazars may reflect the intrinsic distribution of emitting electrons, which will further give some information on the possible acceleration and cooling processes in the emitting region. The γ-ray spectra of Fermi blazars are normally fitted either by a single power-law(PL) or a log-normal(call Logarithmic Parabola, LP) form. The possible reason for this difference is not clear. We statistically explore this issue based on the different observational properties of 1419 Fermi blazars in the 3 LAC Clean Sample. We find that the γ-ray flux(100 Me V–100 Ge V) and variability index follow bimodal distributions for PL and LP blazars, where the γ-ray flux and variability index show a positive correlation. However, the distributions of γ-ray luminosity and redshift follow a unimodal distribution. Our results suggest that the bimodal distribution of γ-ray fluxes for LP and PL blazars may not be intrinsic and all blazars may have an intrinsically curved γ-ray spectrum, and the PL spectrum is just caused by the fitting effect due to less photons.     

A magnetic confinement nuclear fusion mechanism for solar flares

We propose a magnetic confinement nuclear fusion mechanism for the evolution of a solar flare in the solar atmosphere.The mechanism agrees with two observed characteristics of explosive flares and coronal mass ejections(CMEs) that have proved to be very difficult to explain with previous mechanisms:the huge enrichments of3 He and the high energy gamma ray radiation.The twisted magnetic flux rope is a typical structure during the solar flares,which is closely related to the solar active region that magnetic fields have almost complete control over the plasma.Consequently,the plasma inside the flux rope is heated to more than 1.0×107 K by an adiabatic compression process,and then the thermonuclear fusion can take place in the flux rope accompanied with high energy gamma rays.We utilize the time-dependent ideal 2.5-dimensional magnetohydrodynamic(MHD) simulation to demonstrate the physical mechanism for producing flares,which reveals three stages of flare development with the process of magnetic energy conversion and intense release during the solar flares and CMEs in the solar atmosphere.Furthermore,we discuss the relationship between magnetic reconnection and solar eruptions.     

TeV cosmic-ray proton and helium spectra in the myriad model Ⅱ

Recent observations show that the spectra of cosmic ray nuclei start to harden above ～ 102 GeV, which contradicts the conventional steady-state cosmic ray model. We had suggested that this anomaly is due to the propagation effect of cosmic rays released from local young cosmic ray sources; the total flux of cosmic rays should be computed with the Myriad Model, where a contribution from sources in the local catalog is added to the background. However, although the hardening could be elegantly explained in this model, the model parameters obtained from the fit are skewed toward a region with fast diffusion and a low supernova rate in the Galaxy, in disagreement with other observations. We further explore this model in order to set up a concordant picture. Two possible improvements related to cosmic ray sources are considered. First, instead of the usual axisymmetric disk model, we examine a spiral model for the source distribution. Second, for nearby and young sources which are necessary to explain the hardening, we allow for an energy-dependent escape. We find that a major improvement comes from incorporating an energy-dependent escape time for local sources, and with both modifications not only are the cosmic ray proton and helium anomalies resolved, but also the parameters attain values in a reasonable range compatible with other analyses.     

Size-Flux Relation in Solar Active Regions

We present a study of the relationship between integral area and corresponding total magnetic flux for solar active regions. It is shown that some of these relationships are satisfied to simple power laws. Fractal examination showed that some of these power laws can not be justified inside the simple models of stationary magnetic flux tube aggregation. All magnetic fluxes and corresponding areas were calculated using the data measured with the Solar Magnetic Field Telescope of the Huairou Solar Observing Station in Beijing.     

Studies of grand minima in sunspot cycles by using a flux transport solar dynamo model

We propose that grand minima in solar activity are caused by simultane- ous fluctuations in the meridional circulation and the Babcock-Leighton mechanism for the poloidal field generation in the flux transport dynamo model. We present the following results： （a） fluctuations in the meridional circulation are more effective in producing grand minima; （b） both sudden and gradual initiations of grand minima are possible; （c） distributions of durations and waiting times between grand minima seem to be exponential; （d） the coherence time of the meridional circulation has an effect on the number and the average duration of grand minima, with a coherence time of about 30 yr being consistent with observational data. We also study the occurrence of grand maxima and find that the distributions of durations and waiting times between grand maxima are also exponential, like the grand minima. Finally we address the question of whether the Babcock-Leighton mechanism can be operative during grand minima when there are no sunspots. We show that an a-effect restricted to the upper portions of the convection zone can pull the dynamo out of the grand minima and can match various observational requirements if the amplitude of this a-effect is suitably fine-tuned.     

Similarity of jet radiation between flat spectrum radio quasars and GeV narrow-line Seyfert 1 galaxies: a universal δ-L_c correlation

By modeling the broadband spectral energy distributions(SEDs) of a typical flat spectrum radio quasar(FSRQ, 3C 279) and two GeV narrow-line Seyfert 1 galaxies(NLS1s, PMN J0948+0022 and 1H0323+342) in different flux stages with one-zone leptonic models, we find a universal correlation between their Doppler factors(δ) and peak luminosities(L_c) of external Compton scattering bumps. Compiling a combined sample of FSRQs and GeV NLS1s, it is found that both FSRQs and GeV NLS1s in different stages and in different sources follow the same δ-L_c correlation well. This indicates that the variations of observed luminosities may be essentially due to the Doppler boosting effect. The universal δ-L_c relation between FSRQs and GeV NLS1s in different stages may be further evidence that the particle acceleration and radiation mechanisms for the two kinds of sources are similar. In addition, by replacing L_c with the observed luminosity in the Fermi/LAT band(L_(LAT)), this correlation holds and it may serve as an empirical indicator of δ. We estimate the δ values with L_(LAT) for 484 FSRQs in the Fermi/LAT Catalog and they range from 3 to 41, with a median of 16, which are statistically consistent with the values derived by other methods.     

The core dominance parameter and Fermi detection of extragalactic radio sources

By cross-correlating an archive sample of 542 extragalactic radio sources with the Fermi-LAT Third Source Catalog(3FGL), we have compiled a sample of 80 γ-ray sources and 462 non-Fermi sources with available core dominance parameter(RCD), and core and extended radio luminosity; all the parameters are directly measured or derived from available data in the literature. We found that RCD has significant correlations with radio core luminosity, γ-ray luminosity and γ-ray flux; the Fermi sources have on average higher RCD than non-Fermi sources. These results indicate that the Fermi sources should be more compact,and the beaming effect should play a crucial role in the detection of γ-ray emission. Moreover, our results also show Fermi sources have systematically larger radio flux than non-Fermi sources at fixed RCD, indicating larger intrinsic radio flux in Fermi sources. These results show a strong connection between radio andγ-ray flux for the present sample and indicate that the non-Fermi sources are likely due to the low beaming effect, and/or the low intrinsic γ-ray flux. This supports a scenario that has been published in the literature:a co-spatial origin of the activity for the radio and γ-ray emission, suggesting that the origin of the seed photons for the high-energy γ-ray emission is within the jet.     

Visible and Near-Infrared Contrast of Faculae in Active Region NOAA 8518

We compare the contrast of faculae, in visible light and in the near infrared (NIR), that were associated with the active region NOAA 8518 which crossed the solar disk from April 19 to 27, 1999. We obtained NIR continuum images at 1.6 μm at the Big Bear Solar Observatory (BBSO) with an Indium Gallium Arsenide (In Ga As) NIR digital camera. We also obtained high-resolution longitudinal magnetograms and visible light filtergrams at 610.3 nm with the newly developed Digital Vector Magnetograph (DVMG). Our data show that the contrast of faculae has the same sign in both the visible and the NIR. We did not find any so-called “dark faculae”, faculae that are bright in the visible and simultaneously dark in the NIR.We determined a threshold magnetic flux density that separates pores from faculae.     

Time-resolved spectral analysis of prompt emission from long gamma-ray bursts with GeV emission

We performed detailed time-resolved spectroscopy of bright long gammaray bursts(GRBs)which show significant GeV emissions(GRB 080916C,GRB090902B and GRB 090926A).In addition to the standard Band model,we also use a model consisting of a black body and a power law to fit the spectra.We find that for the latter model there are indications of an additional soft component in the spectra.While previous studies have shown that such models are required for GRB 090902B,here we find that a composite spectral model consisting of two blackbodies and a power law adequately fits the data of all the three bright GRBs.We investigate the evolution of the spectral parameters and find several interesting features that appear in all three GRBs,like(a)temperatures of the blackbodies are strongly correlated with each other,(b)fluxes in the black body components are strongly correlated with each other,(c)the temperatures of the black body trace the profile of the individual pulses of the GRBs,and(d)the characteristics of power law components like the spectral index and the delayed onset bear a close similarity to the emission characteristics in the GeV regions.We discuss the implications of these results and the possibility of identifying the radiation mechanisms during the prompt emission of GRBs.     

A Parametric Survey of the CME Triggering Process by Numerical Simulations

Observations indicate that solar coronal mass ejections (CMEs) are closely associated with reconnection-favored flux emergence, which was explained in the emerging flux trigger mechanism for CMEs by Chen & Shibata based on numerical simulations. We present a parametric survey of the triggering agent: its polarity orientation, position, and the amount of the unsigned flux. The results suggest that whether a CME can be triggered depends on both the amount and location of the emerging flux, in addition to its polarity orientation. A diagram is presented to show the eruption and non-eruption regimes in the parameter space. The work is aimed at providing useful information for the space weather forecast.     

Dependence of the length of solar filament threads on the magnetic configuration

High-resolution Hαobservations indicate that filaments consist of an assembly of thin threads.In quiescent filaments,the threads are generally short,whereas in active region filaments,the threads are generally long.In order to explain these observational features,we performed one-dimensional radiative hydrodynamic simulations of filament formation along a dipped magnetic flux tube in the framework of the chromospheric evaporation-coronal condensation model.The geometry of a dipped magnetic flux tube is characterized by three parameters,i.e.,the depth(D),the half-width(w)and the altitude(h)of the magnetic dip.A survey of the parameters in numerical simulations shows that when allowing the filament thread to grow in 5 days,the maximum length(Lth)of the filament thread increases linearly with w,and decreases linearly with D and h.The dependence is fitted into a linear function Lth=0.84w-0.88D-2.78h+17.31(Mm).Such a relation can qualitatively explain why quiescent filaments have shorter threads and active region filaments have longer threads.     

Quasars in the Galactic Anti-Center Area from LAMOST DR3

We present a sample of quasars discovered in an area near the Galactic Anti-Center covering150?≤ l ≤ 210?and |b| ≤ 30?,based on LAMOST Data Release 3(DR3).This sample contains 151 spectroscopically confirmed quasars.Among them 80 are newly discovered with LAMOST.All these quasars are very bright,with i magnitudes peaking around 17.5 mag.All the new quasars were discovered serendipitously from objects that were originally targeted with LAMOST as stars having bluer colors,except for a few candidates targeted as variable,young stellar objects.This bright quasar sample at low Galactic latitudes will help fill the gap in the spatial distribution of known quasars near the Galactic disk that are used to construct an astrometric reference frame for the purpose of accurate proper motion measurements that can be applied to,for example,Gaia.They are also excellent tracers to probe the kinematics and chemistry of the interstellar medium in the Milky Way disk and halo via absorption line spectroscopy.     

History of Star Formation and Chemical Enrichment in the Milky Way Disk

Based on a physical treatment of the star formation law similar to that given by Efstathiou, we have improved our two-component chemical evolution model for the Milky Way disk. Two gas infall rates are compared, one exponential, one Gaussian. It is shown that the star formation law adopted in this paper depends more strongly on the gas surface density than that in Chang et al. It has large effects on the history of star formation and gas evolution of the whole disk. In the solar neighborhood, the history of chemical evolution and star formation is not sensitive to whether the infall rate is Gaussian or exponential. For the same infall time scale, both forms predict the same behavior for the current properties of the Galactic disk. The model predictions do depend on whether or not the infall time scale varies with the radius, but current available observations cannot decide which case is the more realistic. Our results also show that it would be inadequate to describe the gradient evolution along the Gala     

Astrometric calibration of the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center （XSTPS-GAC）

We present astrometric calibration of the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center(XSTPS-GAC). XSTPS-GAC is the photometric part of the Digital Sky Survey of the Galactic Anti-center(DSS-GAC),which is a photometric and spectroscopic sky survey, in combination with LAMOST.In order to select an astrometric reference catalog, we made comparisons between the four widely used astrometric catalogs, GSC2.3, USNO-B1.0, UCAC3 and PPMXL.PPMXL shows relatively small systematic errors in positions and more homogeneous proper motion distributions toward the Galactic Anti-center(GAC), and was selected as the reference catalog. Based on the high quality and bright reference stars that were picked out from PPMXL, we performed a 4th-order polynomial fitting in image units,to construct the transformation relation between coordinates used by XSTPS-GAC and standard coordinates, and to simultaneously correct the image distortions in the CCD. Then we applied the derived relation to all sources to obtain their mean celestial coordinates based on the International Celestial Reference System. For bright point sources with r 17.0 mag, the accuracy of astrometric calibration could reach about80 mas for each of the g, r, i bands, with systematic errors being less than 10 mas. But for the faint sources at the brightness limit of the survey, which was r ～ 19.0 mag, the accuracy can still reach 200 mas. After combining all observations, the final weighted average coordinates could reach an accuracy of less than 70 mas for bright stars. For faint stars, the rms residuals of weighted coordinates decrease to ～ 110 mas. The final combined XSTPS-GAC coordinates show a good consistency with the Sloan Digital Sky Survey.