In this work, we collect quasi-simultaneous infrared, optical, X-ray and $\gamma$-ray data of 60 $Fermi$-4LAC flat spectrum radio quasars (FSRQs). In the framework of the conventional one-zone leptonic model, we investigate the physical properties of $Fermi$-4LAC FSRQs' jets by modeling their quasi-simultaneous spectral energy distributions (SEDs). Our main results are summarized as follows. (1) There is a linear correlation between synchrotron peak frequency and curvature of the electron energy distribution. As suggested by previous works, the slope of the best linear fitting equation of this correlation is consistent with statistic acceleration which needs a fluctuation of fractional acceleration gain. (2) The gamma-ray dissipation regions are located at the range from 0.1 to 10 pc away from the super-massive black hole, and located outside the broad-line region (BLR) and within the dusty torus (DT). (3) A size relation $P_{\rm e}$ (the kinetic power carried in relativistic electrons) $\sim$ $P_{\rm B}$ (Poynting flux) $\leq$ $P_{\rm r}$ (the radiative power ) $<$ $P_{\rm p}$ (the kinetic power in cold protons) is found in our modeling. Among them, $P_{\rm e}\sim P_{\rm B}$ suggests that SEDs of almost all FSRQs with parameters are close to equipartition between the magnetic field and the relativistic electrons. The $P_{\rm e} < P_{\rm r}$ suggest that the most energy of the relativistic electrons are dissipated by EC radiation for FSRQs. (4) There is an anti-correlation between the peak energy of SEDs ($\gamma_{\rm peak}$) and the jet power ($P_{\rm jet}$), which is consistent with the blazar sequence.

We collect a sample of 78 Fermi detected blazars with broad line region (BLR) data and the quasi-simultaneous multi-wave band data. By analyzing the sample, we find: (1) For whole blazar sample, there exist significant correlations between radio, $\gamma$-ray luminosity and BLR luminosity. The slope of the best-fit linear regression equation is close to the theoretical value 1. Our results provided another support of jet-disk symbiosis. (2) For FSRQ sample, we find a significant correlation between radio luminosity and optical luminosity, and weaker correlations between X-ray luminosity and radio, optical luminosity. It indicates that the X-ray band radiation mechanism is slightly different from radio and optical bands. The X-ray emission is diluted by inverse Compton scattering (IC). There is a significant correlation between BLR luminosity and $\gamma$-ray luminosity, but no correlations between X-ray luminosity and BLR, $\gamma$-ray luminosity, which suggests that the seed soft photons are different for IC process. From these results we suggest that the X-ray band emission is diluted by synchrotron self-Compton (SSC) process, $\gamma$-ray band emission is produced by external Compton (EC) process.

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