The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 \mu\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

We observed a nearby radio galaxy, the Centaurus A (Cen A), three times with Suzaku in 2009, and measured the wide-band X-ray spectral variability more accurately than the previous measurements. The Cen A was in the active phase in 2009, and the flux became higher by a factor of 1.5--2.0 and the spectrum became harder than that in 2005. The Fe-K line intensity increased by 20--30% from 2005 to 2009. The correlation of the count rate between the XIS 3--8 keV and PIN 15--40 keV band showed a complex behavior with a deviation from a linear relation. The wide-band X-ray continuum in 2--200 keV can be fitted with an absorbed powerlaw model plus a reflection component, or a powerlaw with a partial covering Compton-thick absorption. The difference spectra between high and low flux periods in each observation were reproduced by a powerlaw with a partial covering Compton-thick absorption. Such a Compton-thick partial covering absorber was for the first time observed for the Cen A. The powerlaw photon index of the difference spectra in 2009 is almost the same as that of the time-averaged spectra in 2005, but steeper by $\sim0.2$ than that of the time-averaged spectra in 2009. This suggests an additional hard powerlaw component with a photon index of $<1.6$ in 2009. This hard component could be a lower part of the inverse-Compton-scattered component from the jet, whose gamma-ray emission has recently been detected with the Fermi/LAT.

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

Improvements of in-orbit calibration of GSO scintillators in the Hard X-ray Detector on board Suzaku are reported. To resolve an apparent change of the energy scale of GSO which appeared across the launch for unknown reasons, consistent and thorough re-analyses of both pre-launch and in-orbit data have been performed. With laboratory experiments using spare hardware, the pulse height offset, corresponding to zero energy input, was found to change by ~0.5 of the full analog voltage scale, depending on the power supply. Furthermore, by carefully calculating all the light outputs of secondaries from activation lines used in the in-orbit gain determination, their energy deposits in GSO were found to be effectively lower, by several percent, than their nominal energies. Taking both these effects into account, the in-orbit data agrees with the on-ground measurements within ~5%, without employing the artificial correction introduced in the previous work (Kokubun et al. 2007). With this knowledge, we updated the data processing, the response, and the auxiliary files of GSO, and reproduced the HXD-PIN and HXD-GSO spectra of the Crab Nebula over 12-300 keV by a broken powerlaw with a break energy of ~110 keV.

We systematically analyzed the high-quality Suzaku data of 88 Seyfert galaxies. We obtained a clear relation between the absorption column density and the equivalent width of the 6.4 keV line above 10$^{23}$ cm$^{-2}$, suggesting a wide-ranging column density of $10^{23-24.5}$ cm$^{-2}$ with a similar solid and a Fe abundance of 0.7--1.3 solar for Seyfert 2 galaxies. The EW of the 6.4 keV line for Seyfert 1 galaxies are typically 40--120 eV, suggesting the existence of Compton-thick matter like the torus with a column density of $>10^{23}$ cm$^{-2}$ and a solid angle of $(0.15-0.4)*4pi$, and no difference of neutral matter is visible between Seyfert 1 and 2 galaxies. An absorber with a lower column density of $10^{21-23}$ cm$^{-2}$ for Compton-thin Seyfert 2 galaxies is suggested to be not a torus but an interstellar medium. These constraints can be understood by the fact that the 6.4 keV line intensity ratio against the 10--50 keV flux is almost identical within a range of 2--3 in many Seyfert galaxies. Interestingly, objects exist with a low EW, 10--30 eV, of the 6.4 keV line, suggesting that those torus subtends only a small solid angle of $<0.2*4pi$. Ionized Fe-K$\alpha$ emission or absorption lines are detected from several percents of AGNs. Considering the ionization state and equivalent width, emitters and absorbers of ionized Fe-K lines can be explained by the same origin, and highly ionized matter is located at the broad line region. The rapid increase in EW of the ionized Fe-K emission lines at $N_{H}>10^{23}$ cm$^{-2}$ is found, like that of the cold material. It is found that these features seem to change for brighter objects with more than several $10^{44}$ erg/s such that the Fe-K line features become weak. We discuss this feature, together with the torus structure.

We report on a detailed investigation of the high-energy gamma-ray emission from NGC\,1275, a well-known radio galaxy hosted by a giant elliptical located at the center of the nearby Perseus cluster. With the increased photon statistics, the center of the gamma-ray emitting region is now measured to be separated by only 0.46' from the nucleus of NGC1275, well within the 95% confidence error circle with radius ~1.5'. Early Fermi-LAT observations revealed a significant decade-timescale brightening of NGC1275 at GeV photon energies, with a flux about seven times higher than the one implied by the upper limit from previous EGRET observations. With the accumulation of one-year of Fermi-LAT all-sky-survey exposure, we now detect flux and spectral variations of this source on month timescales, as reported in this paper. The average >100 MeV gamma-ray spectrum of NGC1275 shows a possible deviation from a simple power-law shape, indicating a spectral cut-off around an observed photon energy of E = 42.2+-19.6 GeV, with an average flux of F = (2.31+-0.13) X 10^{-7} ph/cm^2/s and a power-law photon index, Gamma = 2.13+-0.02. The largest gamma-ray flaring event was observed in April--May 2009 and was accompanied by significant spectral variability above E > 1-2 GeV. The gamma-ray activity of NGC1275 during this flare can be described by a hysteresis behavior in the flux versus photon index plane. The highest energy photon associated with the gamma-ray source was detected at the very end of the observation, with the observed energy of E = 67.4GeV and an angular separation of about 2.4' from the nucleus. In this paper we present the details of the Fermi-LAT data analysis, and briefly discuss the implications of the observed gamma-ray spectral evolution of NGC1275 in the context of gamma-ray blazar sources in general.

We report on a study of the large-scale temperature structure of the Perseus cluster with Suzaku, using the observational data of four pointings of 30' offset regions, together with the data from the central region. Thanks to the Hard X-ray Detector (HXD-PIN: 10 - 60 keV), Suzaku can determine the temperature of hot galaxy clusters. We performed the spectral analysis, by considering the temperature structure and the collimator response of the PIN correctly. As a result, we found that the upper limit of the temperature in the outer region is $\sim$ 14 keV, and an extremely hot gas, which was reported for RXJ 1347.5-1145 and A 3667, was not found in the Perseus cluster. This indicates that the Perseus cluster has not recently experienced a major merger.

NGC 4636, an X-ray bright elliptical galaxy, was observed for 70 ks with Suzaku. The low background and good energy resolution of the XIS enable us to estimate the foreground Galactic emission accurately and hence measure, for the first time, the O, Mg, Si and Fe abundances out to a radius of ~28 arcmin ($\simeq$ 140 kpc). These metal abundances are as high as $>$1 solar within the central 4' and decrease by ~50% towards the outer regions. Further, the O to Fe abundance ratio is about 0.60--1.0 solar in all regions analyzed, indicating that the products of both SNe II and SNe Ia have mixed and diffused to the outer regions of the galaxy. The O and Fe metal mass-to-light-ratios (MLR) of NGC 4636 are 2--3 times larger than those of NGC 1399 implying that metal distributions in NGC 4636 are less extended than those in NGC 1399, possibly due to environmental factors, such as frequency of galaxy interaction. We also found that the MLRs of NGC 4636 at 0.1 $r_{180}$ are $\sim$5 times smaller than those of clusters of galaxies, possibly consistent with the correlation between temperature and MLR of other spherically symmetric groups of galaxies. We also confirmed a resonant scattering signature in the Fe$_{XV II}}$ line in the central region, as previously reported using the XMM-Newton RGS.

Suzaku Hard X-ray Detector (HXD) achieved the lowest background level than any other previously or currently operational missions sensitive in the energy range of 10--600 keV, by utilizing PIN photodiodes and GSO scintillators mounted in the BGO active shields to reject particle background and Compton-scattered events as much as possible. Because it does not have imaging capability nor rocking mode for the background monitor, the sensitivity is limited by the reproducibility of the non X-ray background (NXB) model. We modeled the HXD NXB, which varies with time as well as other satellites with a low-earth orbit, by utilizing several parameters, including particle monitor counts and satellite orbital/attitude information. The model background is supplied as an event file in which the background events are generated by random numbers, and can be analyzed in the same way as the real data. The reproducibility of the NXB model depends on the event selection criteria (such as cut-off rigidity and energy band) and the integration time, and the 1sigma systematic error is estimated to be less than 3% (PIN 15--40 keV) and 1% (GSO 50--100 keV) for more than 10 ksec exposure.

Clusters of galaxies are among the best candidates for particle acceleration sources in the universe, a signature of which is non-thermal hard X-ray emission from the accelerated relativistic particles. We present early results on Suzaku observations of non-thermal emission from Abell 3376, which is a nearby on-going merger cluster. Suzaku observed the cluster twice, focusing on the cluster center containing the diffuse radio emission to the east, and cluster peripheral region to the west. For both observations, we detect no excess hard X-ray emission above the thermal cluster emission. An upper limit on the non-thermal X-ray flux of $2.1\times10^{-11}$ erg cm$^{-2}$ s$^{-1}$ (15--50 keV) at the 3$\sigma$ level from a $34\times34$ arcmin$^2$ region, derived with the Hard X-ray Detector (HXD), is similar to that obtained with the BeppoSAX/PDS. Using the X-ray Imaging Spectrometer (XIS) data, the upper limit on the non-thermal emission from the West Relic is independently constrained to be $<1.1\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$ (4$-$8 keV) at the 3$\sigma$ level from a 122 arcmin$^2$ region. Assuming Compton scattering between relativistic particles and the cosmic microwave background (CMB) photons, the intracluster magnetic field $B$ is limited to be $>0.03\mu$G (HXD) and $>0.10\mu$G (XIS).