We installed two sets of Astronomical Site Monitoring System(ASMS) at Lijiang Observatory(GMG), for the running of the 2.4-meter Lijiang optical telescope(LJT) and the 1.6-meter Multi-channel Photometric Survey Telescope (Mephisto). The Mephistro is under construction. ASMS has been running on robotic mode since 2017. The core instruments: Cloud Sensor, All-Sky Camera and Autonomous-DIMM that are developed by our group, together with the commercial Meteorological Station and Sky Quality Meter, are combined into the astronomical optical site monitoring system. The new Cloud Sensor's Cloud-Clear Relationship is presented for the first time, which is used to calculate the All-Sky cloud cover. We designed the Autonomous-DIMM located on a tower, with the same height as LJT. The seeing data have been observed for a full year. ASMS's data for the year 2019 are also analysed in detail, which are valuable to observers.

We continuously monitored the blazar S5 0716+714 in the optical $g$, $r$ and $i$ bands from Nov. 10, 2017 to Jun. 06, 2019. The total number of observations is 201 nights including 26973 data points. This is a very large quasi-simultaneous multicolor sample for the blazar. The average time spans and time resolutions are 3.4 hours and 2.9 minutes per night, respectively. During the period of observations, the target source in the $r$ band brightens from $14^{\rm m}.16$ to $12^{\rm m}.29$ together with five prominent sub-flares, and then first becomes fainter to $14^{\rm m}.76$ and again brightens to $12^{\rm m}.94$ with seven prominent sub-flares. For the long-term variations, we find a strong flatter when brighter (FWB) trend at a low flux state and then a weak FWB trend at a higher flux state. A weak FWB trend at a low flux state and then a strong FWB trend at a higher flux state are also reported. Most of sub-flares show the strong FWB trends, except for two flares with a weak FWB trend. The particle acceleration and cooling mechanisms together with the superposition of different FWB-slopes from sub-flares are likely to explain the optical color behaviours. A scenario of bent jet is discussed.

The origin of the multi-band activities (outbursts/flares) of blazars is still a heavily debated topic. Shock and magnetic reconnection have long been considered as possible triggers for the multi-band activities. In this paper, we present an exploration of the origin of multi-band activities for a high-redshift (z =1.8385) FSRQ PKS 1502+106. Utilizing multi-band data from radio to $\gamma$-ray and optical polarization observations, we investigate two dramatic activities in detail: a $\gamma$-ray dominated outburst in 2015 and an optical dominated outburst in 2017. Our main results are as follows. (I) A fast $\gamma$-ray flare with a flux-doubling time-scale as short as 1-hr in 2015 is discovered. Based on the variability time-scale, the physical parameters of the flaring region (e.g, minimum Doppler factor, emission region size, etc.) are constrained. At the peak of the flare, the $\gamma$-ray spectrum hardens to $\Gamma_{\gamma} = 1.82\pm0.04$ and exhibits an obvious curvature/break characteristic that is caused by the typical "cooling break". Modelings of multi-band SEDs reveal a very hard electronic energy spectrum with the electronic spectral index of $1.07\pm0.53$. This result suggests that this fast $\gamma$-ray flare may be triggered by magnetic reconnection. (II) During the outburst in 2017, the optical polarization degree and optical fluxes show a very tight correlation. By analyzing Stokes parameters of polarization observations, our results show that this outburst could be triggered by a transverse shock with a compression ratio of $\eta> 2.2$, and the magnetic field intensity of the shock emission region is about $0.032$ G.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

A central compact object (CCO, e.g. a black hole) with an accretion disk has been suggested as the common central engine of various astrophysical phenomena, such as gamma-ray bursts (GRBs), tidal disruption events (TDEs) and active galactic nuclei (AGNs). A jet powered by such a system might precess due to the misalignment of the angular momenta of the CCO and accretion disk. Some quasi-periodic behaviors observed in the light curves of these phenomena can be well interpreted within the framework of a precessing jet model. In this paper, we study the emission polarization of precessing jets in the three kinds of phenomena. The polarization angle also shows a gradual change for the synchrotron emission in both the random and toroidal magnetic field configurations with the precessing jet, while it can only change abruptly by $90^\circ$ for the non-precessing top-hat jet. Polarization properties are periodic due to the assumptions made in our model. The polarization observations are crucial to confirm the precession nature of jets in GRBs, TDEs and AGNs.

Starting from hydrodynamic equations, we have established a set of hydrodynamic equations for average flow and a set of dynamic equations of auto- and cross-correlations of turbulent velocity and temperature fluctuations, following the classic Reynold's treatment of turbulence. The combination of the two sets of equations leads to a complete and self-consistent mathematical expressions ready for the calculations of stellar structure and oscillations. In this paper, non-locality and anisotropy of turbulent convection are concisely presented, together with defining and calibrating of the three convection parameters ($c_1$, $c_2$ and $c_3$) included in the algorithm. With the non-local theory of convection, the structure of the convective envelope and the major characteristics of non-adiabatic linear oscillations are demonstrated by numerical solutions. Great effort has been exercised to the choice of convection parameters and pulsation instabilities of the models, the results of which show that within large ranges of all three parameters ($c_1$, $c_2$ and $c_3$) the main properties of pulsation stability keep unchanged.

By using a non-local and time-dependent convection theory, we have calculated radial and low-degree non-radial oscillations for stellar evolutionary models with $M=1.4$--3.0\,$\mathrm{M}_\odot$. The results of our study predict theoretical instability strips for $\delta$ Scuti and $\gamma$ Doradus stars, which overlap with each other. The strip of $\gamma$ Doradus is slightly redder in colour than that of $\delta$ Scuti. We have paid great attention to the excitation and stabilization mechanisms for these two types of oscillations, and we conclude that radiative $\kappa$ mechanism plays a major role in the excitation of warm $\delta$ Scuti and $\gamma$ Doradus stars, while the coupling between convection and oscillations is responsible for excitation and stabilization in cool stars. Generally speaking, turbulent pressure is an excitation of oscillations, especially in cool $\delta$ Scuti and $\gamma$ Doradus stars and all cool Cepheid- and Mira-like stars. Turbulent thermal convection, on the other hand, is a damping mechanism against oscillations that actually plays the major role in giving rise to the red edge of the instability strip. Our study shows that oscillations of $\delta$ Scuti and $\gamma$ Doradus stars are both due to the combination of $\kappa$ mechanism and the coupling between convection and oscillations, and they belong to the same class of variables at the low-luminosity part of the Cepheid instability strip. Within the $\delta$ Scuti--$\gamma$ Doradus instability strip, most of the pulsating variables are very likely hybrids that are excited in both p and g modes.

We have computed linear non-adiabatic oscillations of luminous red giants using a non-local and anisotropic time-dependent theory of convection. The results show that low-order radial modes can be self-excited. Their excitation is the result of radiation and the coupling between convection and oscillations. Turbulent pressure has important effects on the excitation of oscillations in red variables.

Within the framework of non-local time-dependent stellar convection theory, we study in detail the effect of turbulent anisotropy on stellar pulsation stability. The results show that anisotropy has no substantial influence on pulsation stability of g modes and low-order (radial order $n_\mathrm{r}<5$) p modes. The effect of turbulent anisotropy increases as the radial order increases. When turbulent anisotropy is neglected, most of high-order ($n_\mathrm{r}>5$) p modes of all low-temperature stars become unstable. Fortunately, within a wide range of the anisotropic parameter $c_3$, stellar pulsation stability is not sensitive to the specific value of $c_3$. Therefore it is safe to say that calibration errors of the convective parameter $c_3$ do not cause any uncertainty in the calculation of stellar pulsation stability.

In order to search for the evidence of quasi-periodic oscillation (QPO) in blazar, multicolor optical observations of the BL Lacertae object S5 0716+714 were performed from 2005 to 2012. For $I$ band on March 19 2010 with low variability amplitude and low flux level, the same quasi-periodic oscillation $\simeq50$ minutes with 99% significance levels is confirmed by ZDCF method, Lomb-Scargle method, REDFIT and fitting sinusoidal curves. The observed QPO is likely to be explained by accretion disk variability. If the observed QPO indicates an innermost stable orbital period from the accretion disk, the QPO $\simeq50$ min corresponds to a black hole mass of $5.03\times10^6 M_\odot$ for a non-rotating Schwarzschild black hole and $3.2\times10^7 M_\odot$ for a maximally rotating Kerr black hole.

We present our observations of the optical intra-day variability (IDV) in $\gamma$-ray BL Lac object Mrk 501. The observations were run with the 1.02 m and 2.4 m optical telescopes at Yunnan Observatories from 2005 April to 2012 May. The light curve at the $R$ band on 2010 May 15 passes both variability tests (the $F$ test and the ANOVA test). A flare within the light curve on 2010 May 15 has a magnitude change $\Delta m = 0.03 \pm 0.005_{\rm{stat}} \pm 0.007_{\rm{sys}}$ mag, \textbf{a darkening timescale of $\tau_{\rm{d}}=$ 26.7 minutes}, and an amplitude of IDV $Amp=2.9\% \pm0.7\%$. A decline \textbf{described by 11 consecutive flux measurements} within the flare can be fitted linearly with a Pearson's correlation coefficient $r = 0.945$ at the confidence level of $> 99.99\%$. Under the assumptions that the IDV is tightly connected to the mass of the black hole, \textbf{and that the flare duration, being two times $\tau_{\rm{d}}$, is representative of the minimum characteristic timescale, we can derive upper bounds to the mass of the black hole}. In the case of the Kerr black hole, the timescale of $\Delta t_{\rm{min}}^{\rm{ob}}=$ 0.89 hours gives $M_{\bullet}\la 10^{9.20} M_{\odot}$, which is consistent with measurements reported in the literature. This agreement indicates that the hypothesis about $M_{\bullet}$ and $\Delta t_{\rm{min}}^{\rm{ob}}$ is consistent with the measurements/data.

Bolometric luminosity is a basic physical parameter that is widely used in the study of blazars. Due to the lack of simultaneous full wavelength data, several estimators of the bolometric luminosity are being used in practice. In this paper, we study and evaluate the reliability and significance of six estimators, the $5GHz$ luminosity, the 1keV luminosity, the $\gamma$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity+the $\gamma$-ray luminosity and the 1keV luminosity+the $\gamma$-ray luminosity, by analyzing the linear correlations between the integrated bolometric luminosity and them. Our main results are as follows. (i) All the six estimators are reliable in the sense that they are all significant correlated with the bolometric luminosity. (ii) Ranking from the higher significance of the reliability to lower one the six estimators are the the $5GHz$ luminosity+the $\gamma$-ray luminosity, the 1keV luminosity+the $\gamma$-ray luminosity, $\gamma$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity and the 1keV luminosity. (iii) We suggest that the bolometric luminosity can be well estimated by the $\gamma$-ray luminosity using the best linear equation that given in this paper for Fermi FSRQs. (iv) According to the linear regressions obtained in the analysis, we provide calibration for each estimator.

We have monitored the BL Lacertae object S5 0716+714 in the optical bands during 2012 January and February with long time spans on intraday timescales ($>$5 hr) and high time resolutions. During this monitoring period, the object shows violent flaring activity both in short and intraday timescales. The object has high value of duty cycle. The light curves detected as intraday variability (IDV) show variability of various shapes. The variability amplitude is from 12.81 per cent to 33.22 per cent, and the average value $19.92\pm5.87$ per cent. The overall magnitude variabilities are $\bigtriangleup B=1^{\rm m}.24$, $\bigtriangleup V=1^{\rm m}.42$, $\bigtriangleup R=1^{\rm m}.3$, $\bigtriangleup I=1^{\rm m}.23$ respectively. During the observations, the average change rate is $<CR>=0.035\pm0.009$ Mag/h during the ascent and $<CR>=0.035\pm0.014$ Mag/h during the descent. However, different cases are found on certain nights. There are good inter-bands correlations but not significant time lags for intraday and short timescales. The results of the autocorrelation function show that the variability timescales range from 0.054 day to 0.134 day. Most of nights show bluer when brighter (BWB) chromatic trend; a weak redder with brighter (RWB) trend is found; a few nights show no correlations between magnitude and color index. The BWB trend appears in the short timescales. During the flare the spectral index exhibits a clockwise loop for inter-nights. A shock-in-jet model and the shock wave propagating along a helical path are likely to explain the variability and color index variability.

We have monitored the quasar 3C 273 in optical $V$, $R$ and $I$ bands from 2005 to 2016. Intraday variability (IDV) is detected on seven nights. The variability amplitudes for most of nights are less than 10\% and four nights more than 20\%. When considering the nights with time spans $>4$ hours, the value of duty cycle (DC) is 14.17 per cent. Over the twelve years, the overall magnitude and color index variabilities are $\bigtriangleup I=0^{\rm m}.67$, $\bigtriangleup R=0^{\rm m}.72$, $\bigtriangleup V=0^{\rm m}.68$, and $\bigtriangleup (V-R)=0^{\rm m}.25$ respectively. The largest clear IDV has an amplitude of 42% over just 5.8 minutes and the weakest detected IDV is 5.4% over 175 minutes. The BWB (bluer when brighter) chromatic trend is dominant for 3C 273 and appears at different flux levels on intraday timescales. The BWB trend exists for short-term timescales and intermediate-term timescales but different timescales have different correlations. There is no BWB trend for our whole time-series data sets. A significant anti-correlation between BWB trend and length of timescales is found. Combining with $V$-band data from previous works, we find a possible quasi-periodicity of $P=3918\pm1112$ days. The possible explanations for the observed variability, BWB chromatic trend and periodicity are discussed.

We fit the spectral energy distributions (SEDs) of members of a large sample of Fermi 2LAC blazars to synchrotron and inverse Compton (IC) models. Our main results are as follows. (i) As suggested by previous works, the correlation between peak frequency and curvature can be explained by statistical or stochastic particle acceleration mechanisms. For BL Lacs, we find a linear correlation between synchrotron peak frequency and its curvature. The slope of the correlation is consistent with the stochastic acceleration mechanisms and confirm previous studies. For FSRQs, we also find a linear correlation, but its slope cannot be explained by previous theoretical models. (ii) We find a significant correlation between IC luminosity and synchrotron luminosity. The slope of the correlation of FSRQs is consistent with the EC process. And the slope of the correlation of BL Lac is consistent with the SSC process. (iii) We find several significant correlations between IC curvature and several basic parameters of blazars (black hole mass, broad line luminosity, the Lorentz factor of jet). We also find significant correlations between bolometric luminosity and these basic parameters of blazars which suggest that the origin of jet is a mixture of the mechanisms proposed by Blandford $\&$ Znajek and by Blandford $\&$ Payne.

We investigate the physical properties of \textit{Fermi} TeV BL Lac objects jets by modeling the quasi-simultaneous spectral energy distribution of 29 \textit{Fermi} TeV BL Lacs in the frame of a one-zone leptonic synchrotron self-Compton model. Our main results are the following: (i) There is a negative correlation between $B$ and $\delta$ in our sample, which suggests that $B$ and $\delta$ are dependent on each other mainly in Thomson regime. (ii) There are negative correlations between $\nu_{\text{sy}}$ and $r$, the $\nu_{\text{IC}}$ and $r$, which is a signature of the energy-dependence statistical acceleration or the stochastic acceleration. There is a significant correlation between $r$ and $s$, which suggests that the curvature of the electron energy distribution is attributed to the energy-dependence statistical acceleration mechanism. (iii) By assuming one proton per relativistic electron, we estimate the jet power and radiative power. A size relation $P_{\text{e}} \sim P_{\text{p}} > P_{\text{r}} \gtrsim P_{\text{B}}$ is found in our sample. The $P_{\text{e}}>P_{\text{B}}$ suggests that the jets are particle dominated, and the $P_{\text{e}}\sim P_{\text{p}}$ means that the mean energy of relativistic electrons approaches $m_{\text{p}}/m_{\text{e}}$. There are not significant correlations between $P_{\text{jet}}$ and black hole mass in high or low state with a sub-sample of 18 sources, which suggests that the jet power weakly depends on the black hole mass. (iv) There is a correlation between the changes in the flux density at 1 TeV and the changes in the $\gamma_{\text{peak}}$, which suggests the change/evolution of electron energy distribution may be mainly responsible for the flux variation.

We have monitored the BL Lac object Mrk 501 in optical $V$, $R$ and $I$ bands from 2010 to 2015. For Mrk 501, the presence of strong host galaxy component can affect the results of photometry. After subtracting the host galaxy contributions, the source shows intraday and long-term variabilities for optical flux and color indices. The average variability amplitudes of $V$, $R$ and $I$ bands are $22.05\%, 22.25\%, 23.82\%$ respectively, and the value of duty cycle 14.87 per cent. A minimal variability timescale of 106 minutes is detected. No significant time lag between $V$ and $I$ bands is found on one night. The bluer-when-brighter (BWB) trend is dominant for Mrk 501 on intermediate, short and intraday timescales which supports the shock-in-jet model. For the long timescale, Mrk 501 in different state can have different BWB trend. The corresponding results of non-correcting host galaxy contributions are also presented.

We study the $\rm{\gamma}$-ray luminosity and beaming effect for Fermi blazars. Our results are as follows. (i) There are significant correlations between $\rm{\gamma}$-ray luminosity and radio core luminosity, and between $\rm{\gamma}$-ray luminosity and $\rm{R_{v}}$, which suggests that the $\rm{\gamma}$-ray emissions have strong beaming effect. (ii) Using the $\rm{L_{ext}/M_{abs}}$ as an indicator of environment effects, we find that there have no significant correlation between $\rm{\gamma}$-ray luminosity and $\rm{L_{ext}/M_{abs}}$ for all sources when remove the effect of redshift. FSRQs considered alone also do not show a significant correlation, while BL Lacs still show a significant correlation when remove the effect of redshift. These results suggest that the $\rm{\gamma}$-ray emission may be affected by environment on the kiloparsec-scale for BL Lacs.

Based on a large Fermi blazar sample, the blazar sequence (synchrotron peak frequency $\nu_{\rm peak}$ versus synchrotron peak luminosity $L_{\rm peak}$) is revisited. It is found that there is significant anti-correlation between $\nu_{\rm peak}$ and $L_{\rm peak}$ for blazars. However, after being Doppler corrected, the anti-correlation disappears. The jet cavity power ($P_{\rm jet}$) is estimated from extended radio luminosity. So it is free of beaming effect. We find that there are significant anti-correlations between $P_{\rm jet}$ and beam-corrected $\nu_{\rm peak}^{'}$ for both blazars and radio galaxies, which supports the blazar sequence and unification of blazars and radio galaxies (an alternative relationship is the correlation between jet power and $\gamma$-ray photon index).

We study the relation between accretion, black hole mass and jet power in AGN, by using a large group of blazars detected by the Fermi Large Area Telescope and radio galaxies. Our main results are as follows. (i) The jet power of FSRQs and FRII-HEG depends on the black hole mass, which suggests that the FSRQs and FRII-HEG are in Radiation-Pressure Dominated regime. The jet power of BL Lacs and FRI-LEG depends on the accretion, which suggests that the BL Lacs and FRI-LEG are in the Gas-Pressure Dominated regime. (ii) We find that most of FSRQs and BL Lacs have $\rm{P_{jet}>L_{BZ}^{max}}$, which suggests that the Blandford-Znajek mechanism is insufficient to explain the jet power of these objects. (iii) The FSRQs are roughly separated from BL Lacs by the Ledlow-Owen's dividing line in the $\rm{\log P_{jet}-\log M}$ plane, which supports the unified scheme of AGN. (iv) The FSRQs and BL Lacs have a clear division at $\rm{L_{bol}/L_{Edd}\sim0.01}$, and the distribution of Eddington ratios of BL Lacs and FSRQs exhibits a bimodal nature, which imply that the accretion mode of FSRQs may be different from that of BL Lacs. (v) We find a significant correlation between broad line luminosity and jet power, which supports a direct tight connection between jet power and accretion.

A sample of 111 Fermi blazars each with a well-established radio core luminosity, broad-line luminosity, bolometric luminosity and black hole mass has been compiled from the literatures.We present a significant correlation between radio core and broad-line emission luminosities that supports a close link between accretion processes and relativistic jets. Analysis reveals a relationship of $\rm{LogL_{BLR}\sim(0.81\pm0.06)LogL_{R}^{C}}$ which is consistant with theoretical predicted coefficient and supports that blazar jets are powered by energy extraction from a rapidly spinning Kerr black hole through the magnetic field provided by the accretion disk. Through studying the correlation between the intrinsic bolometric luminosity and the black hole mass, we find a relationship of $\rm{{Log}\frac{L_{in}}{L_{\odot}}=(0.95\pm0.26){Log}\frac{M}{M_{\odot}}+(3.53\pm2.24)}$ which supports mass-luminosity relation for Fermi blazars derived in this work is a powerlaw relation similar to that for main-sequence stars. Finally, evolutionary sequence of blazars is discussed.

We compiled a sample of 24 flat-spectrum radio quasars (FSRQs), 21 BL Lacs, 13 FRI and 12 FRII radio galaxies to study the jet power and broad-line luminosity relation, and the jet power and black hole mass relation. Furthermore, we obtain the histograms of key parameters. Our main results are as the follows:(i) We found that the FSRQs are roughly separated from BL Lacs by the Ledlow-Owen FRI/FRII dividing line in the $\rm{\log{P_{jet}}-\log{M}}$ plane. This result supports the unification model that BL Lacs unify with FRI radio galaxies and FSRQs unify with FRII radio galaxies.(ii)We found that the coefficient of the best-fit linear regression equation of $\rm{\log{L_{BLR}}-\log{P_{jet}}}$ relation is very close to 1 for our sample. The correlation between broad line luminosity and jet power is significant which supports that jet power has a close link with accretion.

By statistically analyzing a large sample which includes blazars of Fermi detection (FBs) and non-Fermi detection (NFBs), we find that there are significant differences between FBs and NFBs for redshift, black hole mass, jet kinetic power from "cavity" power, broad-line luminosity, and ratio of core luminosity to absolute V-band magnitude ($R_{\rm v}$), but not for ratio of radio core to extended flux ($R_{\rm c}$) and Eddington ratio. Compared with NFBs, FBs have larger mean jet power, $R_{\rm c}$ and $R_{\rm v}$ while smaller mean redshift, black hole mass, broad-line luminosity. These results support that the beaming effect is main reason for differences between FBs and NFBs, and that FBs are likely to have a more powerful jet. For both Fermi and non-Fermi blazars, there are significant correlations between jet power and the accretion rate (traced by the broad-emission-lines luminosity), between jet power and black hole mass; for Fermi blazars, the black hole mass does not have significant influence on jet power while for non-Fermi blazars, both accretion rate and black hole mass have contributions to the jet power. Our results support the "blazar sequence" and show that synchrotron peak frequency ($\nu_{\rm peak}$) is associated with accretion rate but not with black hole mass.

We have analyzed a large sample of clean blazars detected by Fermi Large Area Telescope (LAT). Using literature and calculation, we obtained intrinsic $\gamma$-ray luminosity excluding beaming effect, black hole mass, broad-line luminosity (used as a proxy for disk luminosity), jet kinetic power from "cavity" power and bulk Lorentz factor for parsec-scale radio emission, and studied the distributions of these parameters and relations between them. Our main results are as follows. (i) After excluding beaming effect and redshift effect, intrinsic $\gamma$-ray luminosity with broad line luminosity, black hole mass and Eddington ratio have significant correlations. Our results confirm the physical distinction between BL Lacs and FSRQs. (ii) The correlation between broad line luminosity and jet power is significant which supports that jet power has a close link with accretion. Jet power depends on both the Eddington ratio and black hole mass. We also obtain $LogL_{\rm BLR}\sim(0.98\pm0.07)Log P_{\rm jet}$ for all blazars, which is consistent with the theoretical predicted coefficient. These results support that jets are powered by energy extraction from both accretion and black hole spin (i.e., not by accretion only). (iii) For almost all BL Lacs, $P_{\rm jet}>L_{\rm disk}$; for most of FSRQs, $P_{\rm jet}<L_{\rm disk}$. The "jet-dominance" (parameterized as $\frac{P_{\rm jet}}{L_{\rm disk}}$) is mainly controlled by the bolometric luminosity. Finally, the radiative efficiency of $\gamma$-ray and properties of TeV blazars detected by Fermi LAT were discussed.

Local mixing-length theory is incapable of describing nonlocal phenomena in stellar convection, such as overshooting. Therefore standard solar models constructed with the local mixing-length theory deviate significantly from the Sun at the boundaries of the convection zone, where convection becomes less efficient and nonlocal effects are important. The differences between observed and computed frequencies come mainly from the near-surface region, while the localized sound-speed difference is just below the convective envelope. In this paper we compute a solar envelope model using Xiong's nonlocal convection theory, and carry out helioseismic analysis. The nonlocal model has a smooth transition at the base of the convection zone, as revealed by helioseismology. It reproduces solar frequencies more accurately, and reduces the localized sound-speed difference between the Sun and standard solar models.

Helioseismic inversions reveal a major discrepancy in sound speed between the Sun and the standard solar model just below the base of solar convection zone. We demonstrate that this discrepancy is caused by the inherent shortcomings of the local mixing-length theory adopted in the standard solar model. Using a self-consistent nonlocal convection theory, we construct an envelope model of the Sun for sound-speed inversion. Our solar model has a very smooth transition from the convective envelope to the radiative interior; and the convective energy flux changes sign crossing the boundaries of the convection zone. It shows evident improvement over the standard solar model, with a significant reduction in the discrepancy in sound speed between the Sun and local convection models.

This paper reports an application of gas-kinetic BGK scheme to the computation of turbulent compressible convection in the stellar interior. After incorporating the Sub-grid Scale (SGS) turbulence model into the BGK scheme, we tested the effects of numerical parameters on the quantitative relationships among the thermodynamic variables, their fluctuations and correlations in a very deep, initially gravity-stratified stellar atmosphere. Comparison indicates that the thermal properties and dynamic properties are dominated by different aspects of numerical models separately. An adjustable Deardorff constant in the SGS model $c_\mu=0.25$ and an amplitude of artificial viscosity in the gas-kinetic BGK scheme $C_2=0$ are appropriate for current study. We also calculated the density-weighted auto- and cross-correlation functions in Xiong's (\cite{xiong77}) turbulent stellar convection theories based on which the gradient type of models of the non-local transport and the anisotropy of the turbulence are preliminarily studied. No universal relations or constant parameters were found for these models.

Under nonlocal convection theory, convection extends without limit therefore no apparent boundary can be defined clearly as in the local theory. From the requirement of a similar structure for both local and non-local models having the same depth of convection zone, and taking into account the driving mechanism of turbulent convection, we argue that a proper definition of the boundary of a convective zone should be the place where the convective energy flux (i.e. the correlation of turbulent velocity and temperature) changes its sign. Therefore, it is convectively unstable region when the flux is positive, and it is convective overshooting zone when the flux becomes negative. The physical picture of the overshooting zone drawn by the usual non-local mixing-length theory is not correct. In fact, convection is already sub-adiabatic ($\nabla<\nabla_{ad}$) far before reaching the unstable boundary; while in the overshooting zone below the convective zone, convection is sub-adiabatic and super-radiative ($\nabla_{rad}<\nabla<\nabla_{ad}$). The transition between the adiabatic temperature gradient and the radiative one is continuous and smooth instead of a sudden switch. In the unstable zone the temperature gradient is approaching radiative rather than going to adiabatic. We would like to claim again that, the overshooting distance is different for different physical quantities......

Using our non-local time-dependent theory of convection, the linear non-adiabatic oscillations of 10 evolutionary model series with masses of 1--3M$_\odot$ are calculated. The results show that there is a red giant instability strip in the lower temperature side of the Hertzsprung-Russel (HR) diagram which goes along the sequences of the red giant branch (RGB) and the asymptotic giant branch (AGB). For red giants of lower luminosities, pulsation instability are found at high order overtones, the lower order modes from the fundamental to the second overtone are stable. Towards higher luminosity and lower effective temperature, instability moves to lower order modes, and the amplitude growth rate of oscillations also grows. At the high luminosity end of the strip, the fundamental and the first overtone become unstable, while all the modes above the 4th order become stable. The excitation mechanism have been sdudied in detail. It is found that turbulent pressure plays a key role for excitating of red variables. The frozen convecttin approximation is unavailaable for the low temperature stars with extended convective envelopes. In any case, this approximation can explain neither the red edge of the Cepheid instability strip, nor the blue edge of the pulsating red giant instability strip. An analytic expression of a pulsation constant as a function of stellar mass, luminosity and effective temperature is presented from this work.

We have performed a linear pulsational stability survey of 6 series of long period variable models with M=1.0 Msun, L=3000 - 8000Lsun, and (X,Z)= (0.700,0.020),(0.735,0.005). The dynamic and thermodynamic couplings between convection and oscillations are treated by using a statistical theory of nonlocal and time-dependent convection. The results show that the fundamental and all the low overtones are always pulsationally unstable for the low temperature models when the coupling between convection and oscillations is ignored. When the coupling is considered, there is indeed a "Mira" pulsationally instability region outside of the Cepheid instability strip on the H-R diagram. The coolest models near the Hayashi track are pulsationally stable. Towards high temperature the fundamental mode first becomes unstable, and then the first overtone. Some one of the 2nd -4th overtone may become unstable for the hotter models. All the modes higher than 4th (n > 4) are pulsationally stable. The position and the width of such an instability region on the H-R diagram critically depends on the mass, luminosity and metal abundance of the star.

Using a nonlocal time-dependent theory of convection, we have calculated the linear non-adiabatic oscillations of the Horizontal Branch (HB) stars, with both the dynamic and thermodynamic coupling between convection and oscillations carefully treated. Turbulent pressure and turbulent viscosity have been included consistently in our equations of non-adiabatic pulsation. When the coupling between convection and oscillations is ignored, for all models with Teff <= 7350, the fundamental through the second overtone are pulsationally unstable; while for Teff <= 6200 all the models are unstable up to (at least) the 9th overtone. When the coupling between convection and oscillations is included, the RR Lyrae instability strip is very well predicted. Within the strip the most models are pulsationally unstable only for the fundamental and the first few overtones. The turbulent viscosity is an important damping mechanism. Being exclusively distinct from the luminous red variables (long period variables), the HB stars to the right of the RR strip are pulsationally stable for the fundamental and low-order overtones, but become unstable for some oen of the high-order overtones. This may provide a valuable clue for the short period, low amplitude red variables found outside the red edge of the RR strip on the H-R diagram of globular clusters. Moreover, we present a new radiation modulated excitation mechanism functioning in a zone of radiation flux gradient. The effects of nonlocal convection and the dynamic coupling between convection and oscillations are discussed. The spatial oscillations of the thermal variables in the pulsational calculations have been effectively suppressed.