A study of active galactic nuclei in low surface brightness galaxies with Sloan Digital Sky Survey spectroscopy

Active galactic nuclei （AGNs） in low surface brightness galaxies （LSBGs） have received little attention in previous studies. We present a detailed spectral analysis of 194 LSBGs from the Impey et al. （1996） APM LSBG sample which has been observed spectroscopically by the Sloan Digital Sky Survey Data Release 5 （SDSS DR5）. Our elaborate spectral analysis enables us to carry out, for the first time, reliable spectral classification of nuclear processes in LSBGs based on the standard emission line diagnostic diagrams in a rigorous way. Star-forming galaxies are common, as found in about 52% of LSBGs. We find that, contrary to some previous claims, the fraction of galaxies that contain AGNs is significantly lower than that found in nearby normal galaxies of high surface brightness. This is qualitatively in line with the finding of Impey et al. This result holds true even within each morphological type from Sa to Sc. LSBGs that have larger central stellar velocity dispersions or larger physical sizes tend to have a higher chance of harboring an AGN. For three AGNs with broad emission lines, the black hole masses estimated from the emission lines are broadly consistent with the well known M-σ. relation established for normal galaxies and AGNs.

Do Radio-loud Active Galactic Nuclei really follow the same M_(BH)-σ_* Relation as Normal Galaxies?

In an examination of the relationship between the black hole mass MBH and stellar velocity dispersion σ＊ in radio-loud active galactic nuclei （AGNs）, we studied two effects which may cause uncertainties in the black hole mass estimates of radio-loud AGNs： the relativistic beaming effect on the observed optical continuum radiation and the orientation effect on the broad emission line width. After correcting these two effects, we re-examined the MBH- σ[OIII] relation for a sample of radio-loud and radio-quiet AGNs, and found the relation for radio-loud AGNs still deviated from that for nearby normal galaxies and radio-quiet AGNs. We also found there is no significant correlation between radio jet power and narrow [OIII] line width, indicating absence of strong interaction between radio jet and narrow line region. It may be that the deviation of the MBH-σ＊ relation of radio-loud AGNs is intrinsic, or that the [OIII] line width is not a good indicator of σ＊ for radio-loud AGNs.

Feedback of Active Galactic Nuclei in Seyfert 2 Galaxies

It is well accepted that feedback from active galactic nuclei （AGNs） plays an important role in the coevolution of the supermassive black hole （SMBH） and its host galaxy, but the concrete mechanism of feedback remains unclear. A considerable body of evidence suggests that AGN feedback suppresses star formation in the host galaxy. We assemble a sample of Seyfert 2 galaxies with recent observational data of compact nuclear starbursts and estimate the gas surface density as a function of column density to illuminate the relation between feedback and AGN properties. Although there are some uncertainties, our data still imply the deviation from the star formation law （Kennicutt-Schmidt law）. Further, they indicate that： （1） Feedback correlates with the Eddington ratio, rather than with the mass of SMBH, as a result of decreasing star formation efficiency. （2） The SMBH and the torus are probably undergoing coevolution. Conclusions presented here can be refined through future high resolution CO or HCN observations.

Four Late-type Galaxies with Double Radio Lobes and Properties of Such Galaxies

The triggering mechanism for radio lobes from late-type galaxies is not fully understood.More samples are desired for a thorough investigation and statistics.By utilizing the optical data from the newly released Dark Energy Spectroscopic Instrument imaging surveys and the radio sources from the NRAO VLA Sky Survey and the Faint Images of the Radio Sky at Twenty-centimeter,we identify four Late-type Galaxies with double Radio Lobes(La GRLs):J0217-3645,J0947+6220,J1412+3723 and J1736+5108.Including previously known La GRLs,we confirm the correlation between radio power P_(1.4GHz)and stellar mass M_(*)of host galaxies.Most(25/35)La GRLs belong to the blue cloud galaxies,while the newly identified cases in this work are located within the region of the red sequence.We find a clear correlation between the differential radio power,i.e.,the offset from the P_(1.4GHz)-M_(*)relation,and the galaxy color,indicating that bluer galaxies at a fixed M_(*)tend to host more powerful radio lobes.Furthermore,the majority(31/36)of La GRLs are either located in a galaxy group or displaying a disturbed morphology.We suggest that all of the galaxy mass,color and surrounding environment could play important roles in triggering radio lobes in late-type galaxies.

Stellar Populations of AGN-host Dwarf Galaxies Selected with Different Methods

In this paper we investigate the stellar populations and star formation histories of 235 active galactic nucleus(AGN)-host dwarf galaxies,consisting of four samples identified separately with different methods(i.e.,radio,X-ray,mid-IR and variability),utilizing the synthesis code STARLIGHT and spectra from the Sloan Digital Sky Survey Data Release 8.Our results show that the variability sample is the oldest,while the mid-IR sample is the youngest,for which the luminosity at 4020?is dominated(>50%)by the young population(t<10~8yr).The light-weighted mean stellar age of the whole sample is in general about 0.7 dex younger than the optical sample studied in Cai et al.We compare the population results between fitting models with and without a power-law(PL)component and find that the neglect of a PL component would lead to an under-and over-estimation by 0.2 and0.1 dex for the light-and mass-weighted mean stellar age,respectively,for our sample of dwarf galaxies,which has a mean fractional contribution of~16%from the AGN.In addition,we obtain further evidence for a possible suppression of star formation in the host galaxy by the central AGN.We also find that there exists an anticorrelation between the extinction-corrected[O III]luminosity and light-weighted mean stellar age,confirming our previous finding that there is a physical connection between AGN and star-forming activities in AGN-host dwarfs.

On the M-σ Relationship and SMBH Mass Estimates of Selected Nearby Galaxies

Super massive black holes are believed to influence galactic evolution and dynamics. A histogram of SMBH masses for different redshift regimes may reveal clues on how the SMBH evolve in time. A prominent method for SMBH mass estimation is based on the linear correlation between the bulge velocity dispersion and the SMBH mass. Known as M-σ relationship, this method is known to provide reasonable but not very accurate mass estimates due to considerable scatter in data. In order to increase the precision, we surveyed the literature and gathered SMBH and velocity dispersion data for low redshift (z 0.02) spiral galaxies. We report the M-σ relationship for low redshift spiral galaxies as, By using this refined M-σ relationship we measured 32 SMBH masses and determined upper and lower mass boundaries and the mass histogram for spiral galaxies in a narrow redshift regime (0.016 z < 0.017). The spectroscopic data are obtained from The SLOAN Digital Survey and The National Observatory of Turkey (TUG). The targets are selected within a low redshift range for discernible [OIII] lines. TUG observations are carried out on the RTT150 1.5 m telescope using TUG Faint Object Spectrographic Camera and the SLOAN data are obtained from the 7th data release of the survey. We measured the bandwidths of narrow [OIII] lines, which are shown to be indicative in estimating stellar bulge velocity dispersion and estimated the central black hole masses from the refined version of the empirical M-σ relationship. The estimated masses vary between 9.51 × 106 - 2.36 × 108 solar masses.

Comparisons of Jet Properties between GeV Radio Galaxies and Blazars

We compile a sample of spectral energy distributions （SEDs） of 12 GeV radio galaxies （RGs）, including eight FR I RGs and four FR II RGs. These SEDs can be represented with the one-zone leptonic model. No significant unification, as expected in the unification model, is found for the derived jet parameters between FR I RGs and BL Lacertae objects （BL Lacs） and between FR II RGs and flat spectrum radio quasars （FSRQs）. However, on average FR I RGs have a larger ＇Tb （break Lorentz factor of electrons） and lower B （magnetic field strength） than FR II RGs, analogous to the differences be- tween BL Lacs and FSRQs. The derived Doppler factors （~） of RGs are on average smaller than those of blazars, which is consistent with the unification model such that RGs are the misaligned parent pop- ulations of blazars with smaller tS. On the basis of jet parameters from SED fits, we calculate their jet powers and the powers carded by each component, and compare their jet compositions and radiation efficiencies with blazars. Most of the RG jets may be dominated by particles, like BL Lacs, not FSRQs. However, the jets of RGs with higher radiation efficiencies tend to have higher jet magnetization. A strong anticorrelation between synchrotron peak frequency and jet power is observed for GeV RGs and blazars in both the observer and co-moving frames, indicating that the ＂sequence＂ behavior among blazars, together with the GeV RGs, may be intrinsically dominated by jet power.

Are Seyfert 2 Galaxies without Polarized Broad Emission Lines More Obscured?

New XMM-Newton data of seven Seyfert 2 galaxies with optical spectropolarimetric observations are presented. An analysis of the 0.5 - 10 keV spectra shows that four Seyfert 2 galaxies with polarized broad lines （PBLs） are absorbed by NH 〈 1024 cm^-2, while two of three Seyfert 2 galaxies without PBLs show evidence of Compton-thick obscuration, supporting the conclusion that Seyfert 2 galaxies without PBLs are more obscured than those with PBLs, Adding the measured obscuration indicators （NH, T ratio, and Fe Kα line EW） of six luminous AGNs to our previous sample improves the significance level of the difference in absorption from 92.3% to 96.3% for NH, 99.1% to 99.4% for T ratio, and 95.3% to 97.4% for Fe Kα line EW. The present results support and enhance the suggestions that the absence of PBLs in Seyfert 2 galaxies can be explained by larger viewing angles of the line of sight to the putative dusty toms, which leads to the obscuration of the broad-line scattering screen, as expected in the unification model.

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

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

Optical identification of radio-loud active galactic nuclei in the ROSAT-Green-Bank sample with SDSS spectroscopy

Results of extended and refined optical identification of 181 radio/X-ray sources in the RASS-Green Bank （RGB） catalog are presented （Brinkmann et al. 1997） which have been spectroscopically observed in the Sloan Digital Sky Survey （SDSS） DR5. The SDSS spectra of the optical counterparts are modeled in a careful and selfconsistent way by incorporating the host galaxy＇s starlight. Optical emission line parameters are presented, which are derived accurately and reliably, along with the radio 1.4- 5 GHz spectral indices estimated using （non-simultaneous） archival data. For 72 sources, the identifications are presented for the first time. It is confirmed that the majority of strong radio/X-ray emitters are radio-loud active galactic nuclei （AGNs）, particularly blazars. Taking advantage of the high spectral quality and resolution and our refined spectral modeling, we are able to disentangle narrow line radio galaxies （NLRGs）, as vaguely termed in most previous identification work, into Seyfert II galaxies and LINERs （low-ionization nuclear emission regions）, based on the standard emission line diagnostics. The NLRGs in the RGB sample, mostly belonging to ＇weak line radio galaxies＇, are found to have optical spectra consistent predominantly with LINERs, and only a small fraction with Seyfert II galaxies. A small number of LINERs have radio power as high as 10623 - 10^26 W Hz^- 1 at 5 GHz, being among the strongest radio emitting LINERs known so far. Two sources are identified with radio-loud narrow line Seyfert 1 galaxies （NLS 1 s）, a class of rare objects. The presence is also confirmed of flat-spectrum radio quasars whose radio-optical-X-ray effective spectral indices are similar to those of High-energy peaked BL Lacs （HBLs）, as suggested by Padovani et al., although it is still a debate as to whether this is the case for their actual spectral energy distributions.

Gaseous versus Stellar Velocity Dispersion in Emission-Line Galaxies

We compare the ionized gas velocity dispersion σgas with the stellar velocity dispersion σ＊ in star-forming galaxies, composite galaxies, Low Ionization Nuclear Emissionline Regions （LINERs） and Seyfert 2s, compiled from a cross-identification of Sloan Digital Sky Survey Fourth Data Release （SDSS DR4） and Point Source Catalogue Redshift Survey （PSCz）. We measure σgas from the FWHMs of emission lines （Hα, [NⅡ]λλ6548, 6583 and [SⅡ]λλ6716, 6731）. A significant correlation between the gas and stellar velocity dispersion exists, despite substantial scatter. The mean value of the gas to stellar velocity dispersion ratio is close to unity. This suggests that gas velocity dispersion can substitute for the stellar velocity dispersion as a tracer of the gravitational potential well for all the four types of galaxies, but the involved uncertainties are different from type to type. We also studied -↑σgas, as a function of the redshift and the axial ratio to test the effects of aperture and galaxy inclination, and found that both effects are weak. Finally we checked the trend of -↑σgas/σ＊ with the infrared luminosity and found no significant correlation.

The pc-scale radio structure of MIR-observed radio galaxies

We investigated the relationship between the accretion process and jet properties by utilizing very long baseline array （VLBA） and mid-infrared （MIR） data for a sample of 45 3CRR radio galaxies selected with a flux density at 178MHz 〉 16.4Jy, 5 GHz very large array （VLA） core flux density ≥7 mJy and MIR observations. The pc-scale radio structure at 5 GHz is presented by using our VLBA observations for 21 sources acquired in February, 2016, the analysis of archival data for 16 objects and directly obtaining measurements for eight radio galaxies available from literatures. The accretion mode is constrained from the Eddington ratio with a dividing value of 0.01, which is estimated from the MIR-based bolometric luminosity and the black hole masses. While most Fanaroff-Riley type II radio galaxies （FRIIs） have higher Eddington ratio than Fanaroff-Riley type I radio galaxies （FRIs）, we found that there is indeed no single correspondence between the FR morphology and accretion mode with eight FRIIs at low accretion rate and two FRIs at high accretion rate. There is a significant correlation between the VLBA core luminosity at 5 GHz and the Eddington ratio. Various morphologies are identi- fied in our sample, including core only, single-sided core-jet and two-sided core-jet structures. We found that the higher accretion rate may be more likely related with the core-jet structure, thus generating a more extended jet. These results imply that the higher accretion rates are likely able to produce more powerful jets. There is a strong correlation between the MIR luminosity at 15 ktm and VLBA 5 GHz core luminosity, in favor of the tight relation between the accretion disk and jets. In our sample, the core brightness temperature ranges from 109 to 101338 K with a median value of 10^11.09K, indicating that systematically the beaming effect may not be significant. The exceptional cases, FRIs at high accretion rates and FRIIs at low accretion rates, are exclusively at the high and low ends, respectively, of the distribution of the flux ratio for VLBA core to 178 MHz flux density. It is not impossible that the locations of these sources are due to the recent shining or weakening of their central engines （i.e., both accretion and jet）.

Geometry of Broad Line Regions of Active Galactic Nuclei

It has long remained an open question as to the geometry of the broad line region （BLR） in active galactic nuclei （AGNs）. The reverberation mapping technique which measures the response of the broad emission lines to the ionizing continuum, when combined with multiwavelength continuum fitted by sophisticated accretion disks, provides a way of probing the BLR geometry. We analyze a sample of 35 AGNs, which have been monitored by the reverberation mapping campaign. In view of energy budget, the reverberation-based BH masses are found to be in agreement with those obtained by accretion disk models in two thirds of the present sample while the reverberation mapping methods underestimate the BH masses in about one third of objects, as also suggested by CoUin et al. in a recent work. We point out that there are obviously two kinds of BLR geometry, which are strongly dependent on the Eddington ratio, and separated by the value LBol/LEdd - 0.1. These results prefer a scenario of the disk and wind configuration of the BLR and identify the Eddington ratio as the physical driver regulating the wind in the BLR.

Search for H I emission from superdisk candidates associated with radio galaxies

Giant gaseous layers(termed “superdisks”) have been hypothesized in the past to account for the strip-like radio emission gap(or straight-edged central brightness depression) observed between twin radio lobes, in over a dozen relatively nearby powerful Fanaroff-Riley Class II radio galaxies. They could also provide a plausible alternative explanation for a range of observations. Although a number of explanations have been proposed for the origin of the superdisks, little is known about their material content. Some X-ray observations of superdisk candidates indicate the presence of hot gas, but a cool dusty medium also seems to be common. If they are entirely or partly composed of neutral gas, then it may be directly detectable and we report here a first attempt to detect/image any neutral hydrogen gas present in the superdisks that are inferred to be present in four nearby radio galaxies. We have not found a positive H I signal in any of the four sources, resulting in tight upper limits on the H I number density in the postulated superdisks,estimated directly from the central rms noise values of the final radio continuum subtracted image. The estimated ranges of the upper limit on neutral hydrogen number density and column density are 10^-4-10^-3 atoms per cm3 and 10^19-10^20 atoms per cm^2, respectively. No positive H I signal is detected even after combining all the four available H I images(with inverse variance weighting). This clearly rules out an H I dominated superdisk as a viable model to explain these structures, however, the possibility of a superdisk being composed of warm/hot gas still remains open.

Jet Formation Mechanism of the Gamma-Ray-emitting Narrow-line Seyfert 1 Galaxies

We use a large sample of gamma-ray narrow-line Seyfert 1 galaxies(γ-NLS1 s) to study the jet formation mechanisms. We find that the jet power of γ-NLS1 s is lower than the maximum jet power of the Blandford–Payne(BP) mechanism. At the same time, we find that there is a significant correlation between jet power and accretion disk luminosity. Moreover, the contribution rates of the accretion to the jet power are larger than that of black hole mass to jet power. These results further suggest that the jet of γ NLS1 s is mainly produced by the BP mechanism.

East Asian VLBI Network observations of active galactic nuclei jets:imaging with KaVA+Tianma+Nanshan

The East Asian Very Long Baseline Interferometry(VLBI) Network(EAVN) is a rapidly evolving international VLBI array that is currently promoted under joint efforts among China,Japan and Korea.EAVN aims at forming a joint VLBI Network by combining a large number of radio telescopes distributed over East Asian regions.After the combination of the Korean VLBI Network(KVN) and the VLBI Exploration of Radio Astrometry(VERA) into Ka VA,further expansion with the joint array in East Asia is actively promoted.Here we report the first imaging results(at 22 and 43 GHz) of bright radio sources obtained with Ka VA connected to Tianma 65-m and Nanshan 26-m Radio Telescopes in China.To test the EAVN imaging performance for different sources,we observed four active galactic nuclei(AGN) having different brightness and morphology.As a result,we confirmed that the Tianma 65-m Radio Telescope(TMRT) significantly enhances the overall array sensitivity,a factor of 4 improvement in baseline sensitivity and 2 in image dynamic range compared to the case of Ka VA only.The addition of the Nanshan 26-m Radio Telescope(NSRT) further doubled the east-west angular resolution.With the resulting high-dynamic-range,high-resolution images with EAVN(Ka VA+TMRT+NSRT),various fine-scale structures in our targets,such as the counter-jet in M87,a kink-like morphology of the 3 C 273 jet and the weak emission in other sources are successfully detected.This demonstrates the powerful capability of EAVN to study AGN jets and to achieve other science goals in general.Ongoing expansion of EAVN will further enhance the angular resolution,detection sensitivity and frequency coverage of the network.

On the reality of broad iron L lines from the narrow line Seyfert 1 galaxies 1H0707–495 and IRAS 13224–3809

We performed time resolved spectroscopy of 1H0707-495 and IRAS 13224-3809 using long XMM-Newton observations. These are strongly variable narrow line Seyfert 1 galaxies and show broad features around 1 keV that have been interpreted as relativistically broad Fe Lc~ lines. Such features are not clearly observed in other active galactic nuclei despite sometimes having high iron abundance required by the best fitted blurred reflection models. Given the importance of these lines, we explore whether the rapid variability of spectral parameters may introduce broad bumps/dips artificially in the time averaged spectrum, which may then be mistaken as broadened lines. We tested this hypothesis by performing time resolved spectroscopy using long （〉 100 ks） XMM-Newton observations and by dividing them into segments with typical exposures of a few ks. We extracted spectra from each such segment and modeled them using a two component phenomenological model consisting of a power law to represent the hard component and a black body to represent the soft emission. As expected, both the sources showed variations in the spectral parameters. Using these variation trends, we simulated model spectra for each segment and then co-added to get a combined simulated spectrum. In the simulated spectra, we found no broad features below 1 keV and in particular no deviation near 0.9 keV as seen in the real averaged spectra. This implies that the broad Fe Lα line that is seen in the spectra of these sources is not an artifact of the variation of spectral components and, hence, provides evidence that the line is indeed genuine.

Observational signatures of close binaries of supermassive black holes in active galactic nuclei

Inspired by the General Relativity for many decades,experimental physicists and astronomers have a solid dream to detect gravitational waves(GWs)from mergers of black holes,which came true until the excellent performance of the Laser Interferometer Gravitational-Wave Observatory(LIGO)at hundreds Hz.Nano-Hz GWs are expected to be radiated by close-binaries of supermassive black holes(CB-SMBHs;defined as those with separations less than^0.1 pc)formed during galaxy mergers and detected through the Pulsar Timing Array(PTA)technique.As of the writing,there remains no nano-Hz GWs detection.Searching for CB-SMBHs is also observationally elusive though there exist a number of possible candidates.In this review,we focus on observational signatures of CB-SMBHs from theoretic expectations,simulations and observations.These signatures appear in energy distributions of multiwavelength continuum,long term variations of continuum,jet morphology,reverberation delay maps and spectroastrometry of broad emission lines,AGN type transitions between type-1 and type-2(changing-look),and gaseous dynamics of circumbinary disks,etc.Unlike hundred-Hz GWs from stellar mass black hole binaries,the waveform chirping of nano-Hz GWs is too slow to detect in a reasonable human timescale.We have to resort to electromagnetic observations to measure orbital parameters of CB-SMBHs to test nano-Hz GW properties.Reverberation mapping is a powerful tool for probing kinematics and geometry of ionized gas in the gravitational well of SMBHs(single or binary)and therefore provides a potential way to determine orbital parameters of CB-SMBHs.In particular,a combination of reverberation mapping with spectroastrometry(realized at the Very Large Telescope Interferometer)will further reinforce this capability.The Atacama Large Millimeter/submillimeter Array(ALMA)and the forthcoming Square Kilometre Array(SKA)are suggested to reveal dynamics of circumbinary disks through molecular emission lines.

Radial stellar populations of AGN-host dwarf galaxies in SDSS-Ⅳ MaNGA survey

Based on MaNGA integral field unit(IFU) spectroscopy we search 60 AGN candidates,which have stellar masses M*≤5 × 10^(9) M⊙ and show AGN ionization signatures in the BPT diagram.For these AGN candidates,we derive the spatially resolved stellar population with the stellar population synthesis code STARLIGHT and measure the gradients of the mean stellar age and metallicity.We find that the gradients of mean stellar age(metallicity) of individual AGN-host dwarfs are diverse in 0-0.5 Re,0.5-1 Re and 0-1 Re.However,the overall behavior of the mean stellar age(metallicity) profiles tend to be flat,as the median values of the gradients are close to zero.We further study the overall behavior of the mean stellar age(metallicity) by plotting the co-added radial profiles for the AGN sample and compare with a control sample with similar stellar mass.We find that the median values of light-weighted mean stellar ages of the AGN sample are as old as 2-3 Gyr within 2 Re,which are about 4-7 times older than those of the control sample.Meanwhile,most of the AGN candidates are low-level AGNs,as only eight sources have L[OⅢ]>1039.5 erg s^(-1).Hence,the AGNs in dwarf galaxies might accelerate the evolution of galaxies by accelerating the consumption of the gas,resulting in an overall quenching of the dwarf galaxies,and the AGNs also become weak due to the lack of gas.The median values of mass-weighted mean stellar age of both samples within 2 Re are similar and as old as about 10 Gyr,indicating that the stellar mass is mainly contributed by old stellar populations.The gradients of co-added mean stellar metallicity for both samples tend to be negative but close to zero,and the similar mean stellar metallicity profiles for both samples indicate that the chemical evolution of the host galaxy is not strongly influenced by the AGN.

Properties of Broad Band Continuum of Narrow Line Seyfert 1 Galaxies

We have performed a statistical study of the properties of the broad band continuum of Narrow Line Seyfert 1 galaxies (NLS1s) by collecting ratio, infrared, optical and X-ray continuum data from various databases and compared the results with control samples of Broad Line Seyfert 1 galaxies (BLS1s). We find that the fraction (similar to 6%) of Radio Loud (RL) NLS1s is significantly less than that of BLS1s (similar to 13%), which is caused by the lack of radio-very-loud sources in the former. The rarity of RL NLS1s, especially radio-very-loud ones, is consistent with the scenario of small black hole and high accretion rate for NLS1s. Six new radio loud NLS1s are found and five RL NLS1 candidates are presented. In comparison with the BLS1s, the NLS1s tend to have stronger far infrared emission, cooler infrared colors and redder B - K color, which suggests that NLS1s are hosted by dust-richer nuclei. The NLS1s also show steeper soft X-ray spectrum and large soft X-ray to optical flux ratio, while a significant fraction show flat soft X-ray spectra. At least two factors can account for this, absorption and spectral variability. We also perform a correlation analysis between various broad band data. It is found that most correlations identified for NLS1s are also valid for radio quiet BLS1s: (1) the optical colors are anti-correlated with X-ray spectral index; (2) higher optical, Xray and NIR luminosity objects show bluer optical colors and red H - K color; (3) higher luminosity objects show warmer IRAS color; (4) the radio loudness correlates with B - K and X-ray to optical flux ratio. Radio loud objects behave somewhat differently in a few correlations.