Gamma-ray Luminosity Research Articles

SIGNIFICANT CONTRIBUTION OF THE CERENKOV LINE-LIKE RADIATION TO THE BROAD EMISSION LINES OF QUASARS

The Cerenkov line-like radiation in a dense gas (N H > 1013 cm–3) is potentially important in the exploration of the optical broad emission lines of quasars and Seyfert 1 galaxies. With this quasi-line emission mechanism, some long standing puzzles in the study of quasars could be resolved. In this paper, we calculate the power of the Cerenkov line-like radiation in dense gas and compare with the powers of other radiation mechanisms by a fast electron to confirm its importance. From the observed gamma-ray luminosity of 3C 279, we show that the total number of fast electrons is sufficiently high to allow effective operation of the quasi-line emission. We present a model calculation for the luminosity of the Cerenkov Lyα line of 3C 279, which is high enough to compare with observations. We therefore conclude that the broad line of quasars may be a blend of the Cerenkov emission line with the real line produced by the bound-bound transition. A new approach to the absorption of the Cerenkov line is presented with the method of escape probability, which markedly simplifies the computation in the optically thick case. The revised set of formulae for the Cerenkov line-like radiation is more convenient in applications.

DISCOVERY OF GeV EMISSION FROM THE CIRCINUS GALAXY WITH THEFERMILARGE AREA TELESCOPE

We report the discovery of gamma-ray emission from the Circinus galaxy using the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. Circinus is a nearby (~4 Mpc) starburst with a heavily obscured Seyfert-type active nucleus, bipolar radio lobes perpendicular to the spiral disk, and kpc-scale jet-like structures. Our analysis of 0.1-100 GeV events collected during 4 years of LAT observations reveals a significant (~ 7.3 sigma) excess above the background. We find no indications of variability or spatial extension beyond the LAT point-spread function. A power-law model used to describe the 0.1-100 GeV gamma-ray spectrum yields a flux of (18.8+/-5.8)x10^{-9} ph cm^{-2} s^{-1} and photon index 2.19+/-0.12, corresponding to an isotropic gamma-ray luminosity of 3 x 10^{40} erg s^{-1}. This observed gamma-ray luminosity exceeds the luminosity expected from cosmic-ray interactions in the interstellar medium and inverse Compton radiation from the radio lobes. Thus the origin of the GeV excess requires further investigation.

ON THE ORIGIN OF X-RAY SPECTRA IN LUMINOUS BLAZARS

Gamma-ray luminosities of some quasar-associated blazars imply jet powers reaching values comparable to the accretion power even if assuming very strong Doppler boosting and very high efficiency of gamma-ray production. With much lower radiative efficiencies of protons than of electrons, and the recent reports of very strong coupling of electrons with shock-heated protons indicated by Particle-in-Cell (PIC) simulations, the leptonic models seem to be strongly favored over the hadronic ones. However, the electron-proton coupling combined with the ERC (External-Radiation-Compton) models of gamma-ray production in leptonic models predict extremely hard X-ray spectra, with energy indices about 0. This is inconsistent with the observed 2-10 keV slopes of blazars, which cluster around an index value of 0.6. This problem can be resolved by assuming that electrons can be cooled down radiatively to non-relativistic energies, or that blazar spectra are entirely dominated by the SSC (Synchrotron-Self Compton) component up to at least 10 keV. Here, we show that the required cooling can be sufficiently efficient only at distances r < 0.03pc. SSC spectra, on the other hand, can be produced roughly co-spatially with the observed synchrotron and ERC components, which are most likely located roughly at a parsec scale. We show that the dominant SSC component can also be produced much further than the dominant synchrotron and ERC components, at distances larger than 10 parsecs. Hence, depending on the spatial distribution of the energy dissipation along the jet, one may expect to see gamma-ray/optical events with either correlated or uncorrelated X-rays. In all cases the number of electron-positron pairs per proton is predicted to be very low. The direct verification of the proposed SSC scenario requires sensitive observations in the hard X-ray band which is now possible with the NuSTAR satellite.

X-RAY MEASUREMENT OF THE SPIN-DOWN OF CALVERA: A RADIO- AND GAMMA-RAY-QUIET PULSAR

We measure spin-down of the 59 ms X-ray pulsar Calvera by comparing the XMM-Newton discovery data from 2009 with new Chandra timing observations taken in 2013. Its period derivative is P-dot = (3.19+/-0.08)e-15, which corresponds to spin-down luminosity E-dot = 6.1e35 erg/s, characteristic age tau_c = P/2P-dot = 2.9e5 yr, and surface dipole magnetic field strength B_s = 4.4e11 G. These values rule out a mildly recycled pulsar, but Calvera could be an orphaned central compact object (anti-magnetar), with a magnetic field that was initially buried by supernova debris and is now reemerging and approaching normal strength. We also performed unsuccessful searches for high-energy gamma-rays from Calvera in both imaging and timing of >100 MeV Fermi photons. Even though the distance to Calvera is uncertain by an order of magnitude, an upper limit of d < 2 kpc inferred from X-ray spectra implies a gamma-ray luminosity limit of < 3.3e32 erg/s, which is less than that of any pulsar of comparable E-dot. Calvera shares some properties with PSR J1740+1000, a young radio pulsar that we show by virtue of its lack of proper motion was born outside of the Galactic disk. As an energetic, high-Galactic-latitude pulsar, Calvera is unique in being undetected in both radio and gamma-rays to faint limits, which should place interesting constraints on models for particle acceleration and beam patterns in pulsar magnetospheres.

The Correlation between the Gamma-Ray Luminosity and the Core-Dominance Parameters for a Fermi Blazar Sample

AbstractIn this work, we investigated the correlation between the γ-ray luminosity, logLγ and the core-dominance parameter, log (1+R), for a sample of 124 Fermi blazars with available core and extended radio emissions. Our analysis shows that there is no correlation between the γ-ray luminosity, log Lγ and the core-dominance parameter, log (1+R). However, there is a closely linear correlation between log Lγ − log LExt and log (1+R), log Lγ−log LExt = (0.95 ± 0.08) log (1+R) + (2.72 ± 0.11), for the whole sample. The result suggests that the γ-ray emissions are composed of two components, one is beamed, the other is unbeamed.

Flat Spectrum Radio Quasar Evolution and the Gamma-ray Background

AbstractWe discuss an analysis of the redshift evolutions and distributions of the gamma-ray luminosity and photon spectral index of flat spectrum radio quasar (FSRQ) type blazars. We utilize data from the Fermi Gamma-ray Space Telescope, with redshfits as determined from optical spectroscopy by Shaw et al. We find that FSRQs have evolved significantly in luminosity but negligibly in photon index, and contribute in toto roughly 20% of the total gamma-ray output of the Universe.

MODELING THE NON-RECYCLEDFERMIGAMMA-RAY PULSAR POPULATION

We use Fermi Gamma-ray Space Telescope detections and upper limits on non-recycled pulsars obtained from the Large Area Telescope (LAT) to constrain how the gamma-ray luminosity L depends on the period P and the period derivative \dot{P}. We use a Bayesian analysis to calculate a best-fit luminosity law, or dependence of L on P and \dot{P}, including different methods for modeling the beaming factor. An outer gap (OG) magnetosphere geometry provides the best-fit model, which is L \propto P^{-a} \dot{P}^{b} where a=1.36\pm0.03 and b=0.44\pm0.02, similar to but not identical to the commonly assumed L \propto \sqrt{\dot{E}} \propto P^{-1.5} \dot{P}^{0.5}. Given upper limits on gamma-ray fluxes of currently known radio pulsars and using the OG model, we find that about 92% of the radio-detected pulsars have gamma-ray beams that intersect our line of sight. By modeling the misalignment of radio and gamma-ray beams of these pulsars, we find an average gamma-ray beaming solid angle of about 3.7{\pi} for the OG model, assuming a uniform beam. Using LAT-measured diffuse fluxes, we place a 2{\sigma} upper limit on the average braking index and a 2{\sigma} lower limit on the average surface magnetic field strength of the pulsar population of 3.8 and 3.2 X 10^{10} G, respectively. We then predict the number of non-recycled pulsars detectable by the LAT based on our population model. Using the two-year sensitivity, we find that the LAT is capable of detecting emission from about 380 non-recycled pulsars, including 150 currently identified radio pulsars. Using the expected five-year sensitivity, about 620 non-recycled pulsars are detectable, including about 220 currently identified radio pulsars. We note that these predictions significantly depend on our model assumptions.

RADIATION MECHANISM AND JET COMPOSITION OF GAMMA-RAY BURSTS AND GeV-TeV-SELECTED RADIO-LOUD ACTIVE GALACTIC NUCLEI

Gamma-ray bursts (GRBs) and GeV–TeV-selected radio-loud active galactic nuclei (AGNs) are compared based on our systematic modeling of the observed spectral energy distributions of a sample of AGNs with a single-zone leptonic model. We show that the correlation between the jet power (Pjet) and the prompt gamma-ray luminosity (Ljet) of GRBs is consistent, within the uncertainties, with the correlation between jet power and the synchrotron peak luminosity (Ls, jet) of flat spectrum radio quasars (FSRQs). Their radiation efficiencies (ε) are also comparable (>10% for most sources), which increase with the bolometric jet luminosity (Lbol, jet) for FSRQs and with the Ljet for GRBs with similar power-law indices. BL Lac objects (BL Lacs) do not follow the Pjet–Ls, jet relation of FSRQs. They have lower ε and Lbol, jet values than FSRQs, and a tentative Lbol, jet–ε relation is also found, with a power-law index different from that of the FSRQs. The magnetization parameters (σ) of FSRQs are on average larger than that of BL Lacs. They are anti-correlated with ε for the FSRQs, but positively correlated with ε for the BL Lacs. GeV narrow-line Seyfert 1 galaxies potentially share similar properties with FSRQs. Based on the analogy between GRBs and FSRQs, we suggest that the prompt gamma-ray emission of GRBs is likely produced by the synchrotron process in a magnetized jet with high radiation efficiency, similar to FSRQs. The jets of BL Lacs, on the other hand, are less efficient and are likely more matter-dominated.

Detection prospects for GeV neutrinos from collisionally heated gamma-ray bursts with IceCube/DeepCore.

Jet reheating via nuclear collisions has recently been proposed as the main mechanism for gamma-ray burst (GRB) emission. In addition to producing the observed gamma rays, collisional heating must generate 10-100 GeV neutrinos, implying a close relation between the neutrino and gamma-ray luminosities. We exploit this theoretical relation to make predictions for possible GRB detections by IceCube + DeepCore. To estimate the expected neutrino signal, we use the largest sample of bursts observed by the Burst and Transient Source Experiment in 1991-2000. GRB neutrinos could have been detected if IceCube + DeepCore operated at that time. Detection of 10-100 GeV neutrinos would have significant implications, shedding light on the composition of GRB jets and their Lorentz factors. This could be an important target in designing future upgrades of the IceCube + DeepCore observatory.

SEARCH FOR GAMMA-RAY EMISSION FROM FOUR ACCRETING MILLISECOND PULSARS WITHFERMI/LAT

We report our search for \gamma-ray emission in the energy range from 100 MeV to 300 GeV from four Accreting Millisecond Pulsars (AMPs), SAX J1808.4-3658, IGR J00291+5934, XTE J1814-338, and XTE J0929-314. The data are from four-year observations carried out by Large Area Telescope (LAT) onboard the Fermi \gamma-ray Space Telescope. The AMPs were not detected, and their \gamma-ray luminosity upper limits we obtain are 5.1*10^33 ergs/s for SAX J1808.4-3658, 2.1*10^33 ergs/s for IGR J00291+5934, 1.2*10^34 ergs/s for XTE J1814-338, and 2.2*10^33 ergs/s for XTE J0929-314. We compare our results with \gamma-ray irradiation luminosities required for producing optical modulations seen from the companions in the AMPs, which has been suggested by Takata et al. (2012), and our upper limits have excluded \gamma-ray emission as the heating source in these systems except XTE J0929-314, the upper limit of which is not deep enough. Our results also do not support the model proposed by Takata et al. (2012) that relatively strong \gamma-ray emission could arise from the outer gap of a high-mass neutron star controlled by the photon-photon pair-creation for the AMPs. Two AMPs, SAX J1808.4-3658 and IGR J00291+5934, have the measurements of their spin-down rates, and we derive the upper limits of their \gamma-ray conversion efficiencies, which are 57% and 3%, respectively. We discuss the implications to the AMP systems by comparing the efficiency upper limit values with that of 20 \gamma-ray millisecond pulsars (MSP) detected by Fermi and the newly discovered transitional MSP binary J1023+0038.

Gamma rays from cosmic rays in supernova remnants

Context: Cosmic rays are thought to be accelerated at supernova remnant (SNR) shocks, but conclusive evidence is lacking. Aims: New data from ground-based gamma-ray telescopes and the Large Area Telescope on the Fermi Gamma-ray Space Telescope are used to test this hypothesis. A simple model for gamma-ray production efficiency is compared with measured gamma-ray luminosities of SNRs, and the GeV to TeV fluxes ratios of SNRs are examined for correlations with SNR ages. Methods: The supernova explosion is modeled as an expanding spherical shell of material that sweeps up matter from the surrounding interstellar medium (ISM). The accumulated kinetic energy of the shell, which provides the energy available for nonthermal particle acceleration, changes when matter is swept up from the ISM and the SNR shell decelerates. A fraction of this energy is assumed to be converted into the energy of cosmic-ray electrons or protons. Three different particle radiation processes---nuclear pion-production interactions, nonthermal electron bremsstrahlung, and Compton scattering---are considered. Results: The efficiencies for gamma-ray production by these three processes are compared with gamma-ray luminosities of SNRs. Our results suggest that SNRs become less gamma-ray luminous at >~ 10^4 yr, and are consistent with the hypothesis that supernova remnants accelerate cosmic rays with an efficiency of ~10% for the dissipation of kinetic energy into nonthermal cosmic rays. Weak evidence for an increasing GeV to TeV flux ratio with SNR age is found.

THE NATURE OF γ-RAY LOUD NARROW-LINE SEYFERT I GALAXIES PKS 1502+036 AND PKS 2004–447

Variable gamma-ray emission has been discovered in five Radio-loud Narrow Line Seyfert 1 (NLSy1) galaxies by the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope. This has clearly demonstrated that these NLSy1 galaxies do have relativistic jets similar to two other cases of gamma-ray emitting Active Galactic Nuclei (AGN), namely blazars and radio galaxies. We present here our results on the multi-band analysis of two gamma-ray emitting NLSy1 galaxies namely PKS 1502+036 (z = 0.409) and PKS 2004-447 (z = 0.240) using archival data. We generate multi-band long term light curves of these sources, build their spectral energy distribution (SED) and model them using an one zone leptonic model. They resemble more to the SEDs of the flat spectrum radio quasar (FSRQ) class of AGN. We then compare the SEDs of these two sources with two other Fermi detected AGN along the traditional blazar sequence, namely the BL Lac Mrk 421 (z = 0.03) and the FSRQ 3C 454.3 (z = 0.86). The SEDs of both PKS 1502+036 and PKS 2004-447 are found to be intermediate to the SEDs of Mrk 421 and 3C 454.3. In the gamma-ray spectral index v/s gamma-ray luminosity plane, both these NLSy1 galaxies occupy a distinct position, wherein, they have luminosity between Mrk 421 and 3C 454.3, however steep gamma-ray spectra similar to 3C 454.3. Their Compton dominance as well as their X-ray spectral slope also lie between Mrk 421 and 3C 454.3. We argue that the physical properties of both PKS 1502+036 and PKS 2004$-$447 are in general similar to blazars and intermediate between FSRQs and BL Lac objects and these sources thus could fit into the traditional blazar sequence.

STELLAR ENCOUNTER RATE IN GALACTIC GLOBULAR CLUSTERS

The high stellar densities in the cores of globular clusters cause significant stellar interactions. These stellar interactions can produce close binary mass-transferring systems involving compact objects and their progeny, such as X-ray binaries and radio millisecond pulsars. Comparing the numbers of these systems and interaction rates in different clusters drives our understanding of how cluster parameters affect the production of close binaries. In this paper we estimate stellar encounter rates (Gamma) for 124 Galactic globular clusters based on observational data as opposed to the methods previously employed, which assumed "King-model" profiles for all clusters. By deprojecting cluster surface brightness profiles to estimate luminosity density profiles, we treat "King-model" and "core-collapsed" clusters in the same way. In addition, we use Monte-Carlo simulations to investigate the effects of uncertainties in various observational parameters (distance, reddening, surface brightness) on Gamma, producing the first catalog of GC stellar encounter rates with estimated errors. Comparing our results with published observations of likely products of stellar interactions (numbers of X-ray binaries, numbers of radio millisecond pulsars, and gamma-ray luminosity) we find both clear correlations and some differences with published results.

LUMINOSITY EVOLUTION OF GAMMA-RAY PULSARS

We investigate the electrodynamic structure of a pulsar outer-magnetospheric particle accelerator and the resultant gamma-ray emission. By considering the condition for the accelerator to be self-sustained, we derive how the trans-magnetic-field thickness of the accelerator evolves with the pulsar age. It is found that the thickness is small but increases steadily if the neutron-star envelope is contaminated by sufficient light elements. For such a light element envelope, the gamma-ray luminosity of the accelerator is kept approximately constant as a function of age in the initial ten thousand years, forming the lower bound of the observed distribution of the gamma-ray luminosity of rotation-powered pulsars. If the envelope consists of only heavy elements, on the other hand, the thickness is greater but increases less rapidly than what a light element envelope has. For such a heavy element envelope, the gamma-ray luminosity decreases relatively rapidly, forming the upper bound of the observed distribution. The gamma-ray luminosity of a general pulsar resides between these two extreme cases, reflecting the envelope composition and the magnetic inclination angle with respect to the rotation axis. The cutoff energy of the primary curvature emission is regulated below several GeV even for young pulsars, because the gap thickness, and hence the acceleration electric field is suppressed by the polarization of the produced pairs.

PHASE-AVERAGED SPECTRA AND LUMINOSITIES OF GAMMA-RAY EMISSIONS FROM YOUNG ISOLATED PULSARS

We study the phase-averaged spectra and luminosities of gamma-ray emissions from young, isolated pulsars within a revised outer gap model. In the revised version of the outer gap, there are two possible cases for the outer gaps: the fractional size of the outer gap is estimated through the photon-photon pair process in the first case (Case I), and is limited by the critical field lines in the second case (Case II). The fractional size is described by Case I if the fractional size at the null charge surface in Case I is smaller than that in Case II, and vice versa. Such an outer gap can extend from the inner boundary, whose radial distance to the neutron star is less than that of the null charge surface to the light cylinder for a gamma-ray pulsar with a given magnetic inclination. When the shape of the outer gap is determined, assuming that high-energy emission at an averaged radius of the field line in the center of the outer gap, with a Gaussian distribution of the parallel electric field along the gap height, represents typical emission, the phase-averaged gamma-ray spectrum for a given pulsar can be estimated in the revised model with three model parameters. We apply the model to explain the phase-averaged spectra of the Vela (Case I) and Geminga (Case II) pulsars. We also use the model to fit the phase-averaged spectra of 54 young, isolated gamma-ray pulsars, and then calculate the gamma-ray luminosities and compare them with the observed data from Fermi-LAT.

Pulsed high-energyγ-rays from thermal populations in the current sheets of pulsar winds

Context. More than one hundred GeV pulsars have been detected up to now by the LAT telescope on the Fermi gamma-ray observatory, showing peak energies around a few GeV. Current modelling proposes that the high energy emission comes from outer magnetospheric gaps, however radiation from the equatorial current sheet which separates the two magnetic hemispheres outside the light cylinder has also been investigated. Aims. In this paper we discuss the region right outside the light cylinder, or "near wind" zone. We investigate the possibility that synchrotron radiation emitted by thermal populations in the equatorial current sheet of the pulsar wind in this region can explain the lightcurves and spectra observed by Fermi/LAT. Methods. We use analytical estimates as well as detailed numerical computation to calculate the gamma-ray luminosities, lightcurves and spectra of gamma-ray pulsars. Results. Many of the characteristics of the gamma-ray pulsars observed by Fermi/LAT can be reproduced by our model, most notably the position of these objects in the P - Pdot diagram, and the range of gamma-ray luminosities. A testable result is a sub-exponential cutoff with an index b = 0.35. We also predict the existence of a population of pulsars with cutoff energies in the MeV range. These have systematically lower spindown luminosities than the Fermi/LAT detected pulsars. Conclusions. It is possible for relativistic populations of electrons and positrons in the current sheet of a pulsar's wind right outside the light cylinder to emit synchrotron radiation that peaks in the sub-GeV to GeV regime, with gamma-ray efficiencies similar to those observed for the Fermi/LAT pulsars.

Thermal and non-thermal traces of AGN feedback: results from cosmological AMR simulations

We investigate the observable effects of feedback from Active Galactic Nuclei (AGN) on non-thermal components of the intracluster medium (ICM). We have modelled feedback from AGN in cosmological simulations with the adaptive mesh refinement code ENZO, investigating three types of feedback that are sometimes called quasar, jet and radio mode. Using a small set of galaxy clusters simulated at high resolution, we model the injection and evolution of Cosmic Rays, as well as their effects on the thermal plasma. By comparing, both, the profiles of thermal gas to observed profiles from the ACCEPT sample, and the secondary gamma-ray emission to the available upper limits from FERMI, we discuss how the combined analysis of these two observables can constrain the energetics and mechanisms of feedback models in clusters. Those modes of AGN feedback that provide a good match to X-ray observations, yield a gamma-ray luminosity resulting from secondary cosmic rays that is about below the available upper limits from FERMI. Moreover, we investigate the injection of turbulent motions into the ICM from AGN, and the detectability of these motions via the analysis of line broadening of the Fe XXIII line. In the near future, deeper observations/upper-limits of non-thermal emissions from galaxy clusters will yield stringent constraints on the energetics and modes of AGN feedback, even at early cosmic epochs.

On proton synchrotron blazar models: the case of quasar 3C 279

In the present work we propose an innovative estimation method for the minimum Doppler factor and energy content of the gamma-ray emitting region of quasar 3C 279, using a standard proton synchrotron blazar model and the principles of automatic photon quenching. The latter becomes relevant for high enough magnetic fields and results in spontaneous annihilation of gamma-rays. The absorbed energy is then redistributed into electron-positron pairs and soft radiation. We show that as quenching sets an upper value for the source rest-frame gamma-ray luminosity, one has, by neccessity, to resort to Doppler factors that lie above a certain value in order to explain the TeV observations. The existence of this lower limit for the Doppler factor has also implications on the energetics of the emitting region. In this aspect, the proposed method can be regarded as an extension of the widely used one for estimating the equipartition magnetic field using radio observations. In our case, the leptonic synchrotron component is replaced by the proton synchrotron emission and the radio by the VHE gamma-ray observations. We show specifically that one can model the TeV observations by using parameter values that minimize both the energy density and the jet power at the cost of high-values of the Doppler factor. On the other hand, the modelling can also be done by using the minimum possible Doppler factor; this, however, leads to a particle dominated region and high jet power for a wide range of magnetic field values. Despite the fact that we have focused on the case of 3C 279, our analysis can be of relevance to all TeV blazars favoring hadronic modelling that have, moreover, simultaneous X-ray observations.

SBS 0846+513: a new γ-ray-emitting narrow-line Seyfert 1 galaxy

We report Fermi-LAT observations of the radio-loud AGN SBS 0846+513 (z=0.5835), optically classified as a Narrow-Line Seyfert 1 galaxy, together with new and archival radio-to-X-ray data. The source was not active at gamma-ray energies during the first two years of Fermi operation. A significant increase in activity was observed during 2010 October-2011 August. In particular a strong gamma-ray flare was observed in 2011 June reaching an isotropic gamma-ray luminosity (0.1-300 GeV) of 1.0x10^48 erg/s, comparable to that of the brightest flat spectrum radio quasars, and showing spectral evolution in gamma rays. An apparent superluminal velocity of (8.2+/-1.5)c in the jet was inferred from 2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. Both the power released by this object during the flaring activity and the apparent superluminal velocity are strong indications of the presence of a relativistic jet as powerful as those of blazars. In addition, variability and spectral properties in radio and gamma-ray bands indicate blazar-like behaviour, suggesting that, except for some distinct optical characteristics, SBS 0846+513 could be considered as a young blazar at the low end of the blazar's black hole mass distribution.

Constrainingγ-ray pulsar gap models with a simulated pulsar population

With the large sample of young gamma-ray pulsars discovered by the Fermi Large Area Telescope (LAT), population synthesis has become a powerful tool for comparing their collective properties with model predictions. We synthesised a pulsar population based on a radio emission model and four gamma-ray gap models (Polar Cap, Slot Gap, Outer Gap, and One Pole Caustic) normalizing to the number of detected radio pulsars in select group of surveys. The luminosity and the wide beams from the outer gaps can easily account for the number of Fermi detections in 2 years of observations. The wide slot-gap beams requires an increase by a factor of ~10 of the predicted luminosity to produce a reasonable number of gamma-ray pulsars. Such large increases in the luminosity may be accommodated by implementing offset polar caps. The narrow polar-cap beams contribute at most only a handful of LAT pulsars. Standard distributions in birth location and pulsar spin-down power (Edot) fail to reproduce the LAT findings: all models under-predict the number of LAT pulsars with high Edot, and they cannot explain the high probability of detecting both the radio and gamma-ray beams at high Edot. The beaming factor remains close to 1 over 4 decades in Edot evolution for the slot gap whereas it significantly decreases with increasing age for the outer gaps. The evolution of the slot-gap luminosity with Edot is compatible with the large dispersion of gamma-ray luminosity seen in the LAT data. The stronger evolution predicted for the outer gap, which is linked to the polar cap heating by the return current, is apparently not supported by the LAT data. The LAT sample of gamma-ray pulsars therefore provides a fresh perspective on the early evolution of the luminosity and beam width of the gamma-ray emission from young pulsars, calling for thin and more luminous gaps.