\section{Introduction}
\subsection{Composite galaxies}
The Baldwin-Phillips-Terlevich (BPT) diagrams (Baldwin et al. 1981) are commonly used as a powerful
tool to classify narrow emission-line galaxies according to their dominant energy sources.
The diagrams are then refined by Veilleux \& Osterbrock (1987). In the diagnostic diagram
[\ion{O}{3}]$\lambda5007$/H$\beta$ vs. [\ion{N}{2}]$\lambda6583$/H$\alpha$,
the star-forming sequence shows a curved concentration extending from up-left
to bottom-right, while AGNs (including LINERs, Low Ionization Nuclear Emission-line Regions)
concentrate on the right part of the diagram. A theoretic demarcation curve determining the upper limit
of the theoretical stellar photoionization model is proposed by Kewley et al. (2001),
and an empirical curve dividing the pure star-forming galaxies recently by Kauffmann et al. (2003a). The
region between the two curves is occupied by the so-called composite objects (hereafter composites
for short) whose spectra are believed
to contain significant contributions from both star formation and AGN (Seyfert or LINER).

As a fundamental, tight relationship, the well-established relation between the
mass of supermassive black hole (SMBH) and velocity dispersion of
the bugle in which the SMBH resides is served as a clue of co-evolution of
AGN and star formation in the bugle
(e.g., Magorrian et al. 1998; Gebhardt et al. 2000; Granato et al. 2001;
Tremaine et al. 2002; Ferrarese \& Merritt 2001; Greene \& Ho 2006). At first,
the issue of co-evolution could be
studied in phenomenology. In fact, there are many observations of both imaging and
spectroscopy that show the evidence of
co-existence of AGN and circumnuclear star formation (e.g., Canalizo \& Stockton 2001;
Cid Fernandes et al. 2001, 2005; Hines et al. 1999; Schinerer et al. 1998;
Heckman et al. 1997; Gonzalez Delgado 2002; Zuther et al. 2007; Sarzi et al. 2007; Asari et al. 2007;
Davies et al. 2007).

We should next answer the
questions: What is the physics that links star formation and AGN phenomena? How AGN co-evolves
with star formation in details? (e.g., Ho 2005a). So far,
there are two distinct scenarios that connect AGN and
star formation phenomena: 1) starburst occurs first and triggers AGN's activity latter (e.g.,
Weedman 1983);
2) in contrast, AGN's activity occurs before the beginning of starburst
(Goncalves et al. 1999 and references therein).
If the hypothesis of co-evolution of AGN and the bulge of its host
is correct, then whatever AGN occurs before/after
starburst, the composites might represent an important state in galaxy evolution,
and are valuable for examining the issue of co-evolution in details.

\subsection{A new approach by Seyfert 1.8 and 1.9 galaxies}

In order to study the co-evolution of AGN and star formation, properties of both
AGN and its host galaxy are required to be examined carefully and simultaneously.
However, discrimination between type I AGN and its host galaxy is a difficult
task, because both the AGN's continuum and broad lines mask the stellar features in optical
wavelength. The problem of AGN-host separation could be partially solved by several
techniques listed as follows:
\begin{itemize}
\item Off-nuclear spectroscopy: By detecting the Balmer absorption lines and emission
lines of \ion{H}{2} region in the surrounding of QSOs,
Canalizo \& Stockton (2001) examined properties of circumnuclear stellar
populations in the hosts of nine local QSOs. The estimated ages range from a few Myrs
old to post-starburst. In a similar way, Nolan et al. (2001) found that the off-nuclear
($\sim$5\symbol{125}) stellar populations of optically selected QSOs are generally
dominated by old ones ($\sim$8-14Gyrs). Relatively old stellar populations ($\sim$0.1-2Gyrs)
are also frequently detected in the off-nuclear regions in a few powerful radio-loud AGNs
(Tadhunter et al. 2005).
\item AGN+post starburst: In this case the star light is easily discriminated in the 1-D spectroscopy
because of the
prominent Balmer absorptions of A-type stars. Brotherton et al. (1999) identified a 400Myrs old
instantaneous starburst in quasar \object{UN\,J1025-0040}. A $\sim100$Myrs old post-starburst is
identified in NLS1-like type I AGN \object{SDSS\,J085338.27+033246.1} (Wang \& Wei 2006).
Wang et al. (2004) identified an evident Balmer jump in type I AGN
\object{SDSS\,J022119.84+005628.4}.
\item Spectroscopy of obscured AGNs: By analyzing the spectra of narrow-line AGNs from Sloan Digital
Sky Survey (SDSS; York et al. 2000),  Heckman et al. (2004 and references therein)
and Kewley et al. (2006)
found that most of the accretion-driven growth of mass of SMBH occurs in the galaxy associated with a relatively
young stellar population (see also in Heckman \& Kauffmann 2006).
\item Infrared color: By analyzing a sample of IRAS-selected type 1.5 AGNs,
Wang et al. (2006) found that the essential Eigenvector 1 (E1) space
(Boroson \& Green 1992, hereafter BG92) is related with the age of stellar population of
the host as assessed by infrared color $\alpha(60,25)$.
\item Deconvolution of 2-D spectra: By performing the MCS deconvolution method on 2-D spectra,
Letawe et al. (2007, and references therein) recently separated the individual spectra of the
quasar and of the underlying host galaxy in a sample
of 20 low-redshift luminous quasars. Quite young stellar populations ($\sim1-2$ Gyrs) are subsequently
identified in a sample of 18 quasars by a semi-analytic decomposition (Jahnke et al. 2007).
\end{itemize}

Although the understanding of the co-evolution of AGN and star formation
has been remarkably developed over the past few years by using these approaches,
their intrinsic dis-advantages prevent the systematical and
extensive investigation on the issue of co-evolution.
We refer the readers to Letawe et al. (2007) for a comment of the drawbacks of these
approaches. In the present study, taking advantage of Seyfert 1.8 and 1.9 galaxies (i.e., partially obscured AGNs),
a new approach is adopted by us to directly obtain
properties of both central AGN and circumnuclear stellar population systematically and
simultaneously. Seyfert 1.8 galaxy is defined as one whose spectrum shows strong narrow components,
and very weak but still visible broad components of H$\alpha$ and H$\beta$.
No broad components of Balmer lines except a weak broad component at H$\alpha$ can be identified in
Seyfert 1.9 galaxy (Osterbrock 1989).
In these two types of AGN,
due to the orientation effect (e.g., Antonucci 1993), their optical continuum emitted from the central accretion
disks is blocked by a torus, and the star light from the bugle is much easier to be detected.
In addition, the broad emission lines shown in the spectrum, rather than some tracers, allow us to derive
AGN's properties directly.

As the first step, by deriving the properties of both AGN and stellar population
from the optical spectra taken by SDSS, we focus on the nature of
the composites with broad Balmer emission lines (hereafter composite
AGNs for abbreviation), and then study their implications on the issue of co-evolution
in this paper.
The issues related with other types of galaxies (e.g., Seyfert and LINER)
will be presented in subsequent paper(s).
At present, SDSS provides a large and uniform spectroscopic
sample of galaxies to make such study available. The paper is organized as follows.
The sample selection and data reduction are
described in Section 2. The physical properties of both AGN and stellar population are
derived in Section 3.
Basing upon the BPT diagrams, a further classification of the composite AGNs is shown in Section 4.
We investigate the AGN-host connection in the composite AGNs in Section 5.
The $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology
with parameter $h=0.7$, $\Omega_{\rm M}=0.3$, and $\Omega_{\Lambda}=0.7$ (Spergel et al. 2003)
is adopted throughout the paper.

\section{Sample Selection and Data Reduction}

\subsection{The Sloan Digital Sky Survey}
The SDSS is designed to eventually provide an imaging and spectroscopic survey of
one-quarter of the entire sky (e.g., York et al. 2000).
The survey is carried out by a dedicated 2.5-meter wide-field (3\symbol{23})
telescope at Apache Point Observatory. The telescope is equipped
with two fiber-fed spectrographes and a mosaic CCD camera.
Each spectrum is taken with a 3\symbol{125} diameter fiber aperture, and covers a wavelength range
from 3800\AA\ to 9200\AA\ at the observe-frame. The spectral resolution varies from 1850 to 2200,
corresponding to $\sigma_{\rm{int}}\sim65 \rm{km\ s^{-1}}$. The spectrophotometric calibrations
with accuracy of about 20\%
are performed by observing subdwarf F stars in each 3\symbol{23} field of view.
Tow pipelines, \sl spectro2d \rm and \sl specBS \rm are developed to reduce
the raw spectra (Glazebrook et al. 1998; Bromley et al. 1998).

\subsection{Selection of Seyfert 1.8 and 1.9 galaxies}

We select a sample of the composite AGNs from the SDSS-Data Release
4 (Adelman-McCarthy et al. 2006) MPA/JHU catalogs (Heckman et al. 2004; Kauffmann et al. 2003a,b,c;
Tremonti et al., 2004)\footnote{These catalogs can be
downloaded from http://www.mpa-garching.mpg.de/SDSS/.}.
The catalogs provide spectral and physical properties for 567,486 individual
narrow emission-line galaxies, including Seyferts, LINERs and star-forming galaxies.

First of all, in order to ensure that our results are not distorted by the spectra
with poor quality, we only consider
the galaxies whose spectra have median signal-to-noise ratio per pixel of the whole spectrum
$\rm{S/N}\geq20$. Secondly, the star-forming galaxies which are under the demarcation line given in
Kauffmann et al. (2003a) are excluded by simply adopting the emission
line intensities provided by the MPA/JHU catalogs. The extinction is not corrected here
since both the line ratios [\ion{O}{3}]/H$\beta$ and [\ion{N}{2}]/H$\alpha$ are not sensitive to
dust extinction. There are total 27,341 galaxies that fulfill these criterion.
An object is then classified as a composite if
\begin{equation}
\rm{\log([OIII]/H\beta)<0.61/(\log([NII]/H\alpha)-0.05)+1.3}
\end{equation}
and
\begin{equation}
\rm{\log([OIII]/H\beta)>0.61/(\log([NII]/H\alpha)-0.47)+1.19}
\end{equation}
Altogether, there are 12,343 objects that could be classified as the composites with median S/N
ratio lager than 20.

An automatic routine is then developed by us to reduce the raw data and to select
objects with broad emission lines attributed to central AGN. At the beginning,
the Galactic extinction is corrected for each spectrum by assuming an $R_V=3.1$
extinction law (Cardelli et al. 1989). The color excess $E(B-V)$ for each spectrum
is taken from the Schlegel, Finkbeiner and Davis Galactic reddening maps
(Schlegel et al. 1998). Each of the extinction-corrected spectra is then
transformed to the rest frame, along with the $k$-correction,
according to the corresponding redshift provided by the SDSS pipelines.

The contribution of star light is then separated from each
spectrum by the principle component analysis (PCA) method
(e.g., Yip et al. 2004; Li et al, 2005; Hao et al. 2005). A library of stellar absorption
spectra is built by applying the PCA technique to the standard
single stellar population (SSP) models developed by Bruszual \& Charlot
(2003, hereafter BC03). The models have a resolution of 3\AA.
The optical spectral evolution is provided for six different metallicities.
At each metallicity, there are total 221 spectra at the unequally
spaced time steps from the age of 0.1Myrs to 20 Gyrs.
The first seven eigenspectra are used to
model the starlight component simply through their linear combination.
The component of AGN's continuum is ignored in the separation in most of the cases,
since the contribution of the observed continuum of AGN is negligible in their spectra (e.g., Kauffmann et al. 2003a),
except in 7 objects. In summary, in order
to appropriately separate the starlight, our template contains the seven
eigenspectra and a Galactic extinction curve (Cardelli et al. 1989) for
most of the cases. Because only the instrumental resolution is taken into
account in the models of BC03, the template is convolved with a Gaussian
function to match the velocity dispersion in each object. The velocity dispersion
is determined by the cross-correlation method before the fitting and subtraction.
For the cases that show emission of both broad H$\beta$
line and \ion{Fe}{2} complex of AGN, both an additional
power-law continuum with form $F_\nu\propto\nu^{-0.5}$ (Richards et al. 2006; Elvis et al. 1994)
and a template of the \ion{Fe}{2} complex are required to model the spectra. We
simply adopt the empirical template of the \ion{Fe}{2} complex given in BG92. The line
width of the template is forced to be identical to that of H$\beta$ broad component (e.g., BG92).

A $\chi^2$ minimization is used to model the starlight component in each spectrum.
The fitting is performed over the rest frame wavelength
range from 3700\AA\ to 7000\AA, except for the regions around the strong emission lines:
e.g., Balmer lines, [\ion{O}{3}]$\lambda\lambda4959,5007$, [\ion{N}{2}]$\lambda\lambda6548,6583$,
[\ion{S}{2}]$\lambda\lambda6716,6731$, [\ion{O}{2}]$\lambda3727$.
The possibly broad components of Balmer lines are also considered to determine
the emission-line free regions. As an illustration, the modeling of the starlight
component is shown in the left column of Figure 1 for three typical cases.
In each panel, the upper and lower curves shows the
starlight-subtracted (emission-line) spectrum and modeled starlight spectrum, respectively.

The emission-line spectra are subsequently used to search for the composite AGNs,
according to the existence of a high velocity wing of
broad H$\alpha$ component at its blue side. The red wing is not
used in the sample selection because of the superposition of the strong [\ion{N}{2}]$\lambda$6583
emission line.
The existence of a broad wing is accepted
if the flux ratio $F_{\rm{w}}/\sigma_{\rm c}\geq3$, where $F_{\rm{w}}$ is the
flux of line wing averaged over the wavelength range from 6500\AA\ to 6530\AA\
at the rest frame, $\sigma_{\rm c}$ the standard deviation of flux within the emission-line
free region between $\lambda$5980 and $\lambda$6020. Iterative tests on
a typical small sample taken from SDSS are performed to determine
the wavelength ranges quoted above. The automatically selected composite AGNs are then inspected
by eyes. Finally, our sample totally contains 85 composite AGNs with broad H$\alpha$ emission
line.

\subsection{Emission line measurements}

The emission lines are measured once the starlight component is subtracted.
We model and measure emission lines of the 85 composite AGNs by the SPECFIT task (Kriss 1994) in the IRAF
package\footnote{IRAF is distributed by the National Optical Astronomical Observatories,
which is operated by the Association of Universities for Research in Astronomy, Inc.,
under cooperative agreement with the National Science Foundation.}.
In each spectrum, each emission line is fitted by a set of Gaussian components. The intensity
ratios of the [\ion{O}{3}] and [\ion{N}{2}] doublets are fixed to their theoretical values (e.g.,
Dimitrijevic et al. 2007).
In Figure 1, the emission-line modeling in H$\beta$ region is schemed in the middle column,
and the modeling in H$\alpha$ region in the right column.
Direct integration is used to measure the flux
of [\ion{O}{1}]$\lambda6300$ emission line in each spectrum through the SPLOT task.
The uncertainties
of the emission-line measurements are dominated by the subtraction of the starlight component, especially
for the faint H$\beta$ emission. Similar as Asari et al. (2007), a wide ($\sim200$\AA) absorption trough
in the continuum around H$\beta$ is frequently identified in our starlight-subtracted spectra, especially
in the galaxies associated with an old stellar population, which leads to a
2\%-7\% underestimate in H$\beta$ flux if the continuum is fixed to be zero.
This small mismatch is probably caused by the calibration in the STELIB library (Asari et al. 2007).
We ignore this small underestimate in this paper
because relatively young stellar populations are frequently identified in our sample
(see Figure 5 for details).

The results of the emission-line measurements are
listed in Table 1. Column (1) lists the identification of each object listed in our sample,
and Column (2) the corresponding redshift. For each object, the line ratios of [\ion{N}{2}]/H$\alpha$,
[\ion{S}{2}]$\lambda\lambda6726,6731$/H$\alpha$, [\ion{O}{1}]/H$\alpha$ and [\ion{O}{3}]/H$\beta$
are listed from Column (3) to Column (6). Column (9) and (12) shows the intrinsic luminosity
of H$\alpha$ broad component ($L_{\mathrm{H\alpha}}$) and that of [\ion{O}{3}]$\lambda$5007
($L_{\mathrm{[OIII]}}$), respectively.
Both luminosities are corrected for local extinction.
The extinction is inferred from the narrow-line ratio H$\alpha$/H$\beta$ for each object, assuming the
Balmer decrement for standard case B recombination and Galactic extinction curve with $R_V=3.1$.

\section{Deriving Physical Properties}

In this section we derive physical properties of both AGN and stellar population
basing upon the reduction and measurements described above.

\subsection{Deriving properties of AGN}

There are compelling evidence suggesting that both Eddington ratio ($L/L_{\rm E}$) and virial mass of SMBH
($M_{\rm{BH}}$) are basic parameters determining the properties of individual AGN (e.g., Boroson 2002; BG92;
Sulentic et al. 2000; Grupe 2004). In this study, we attempt to directly obtain the physical
properties of individual central AGN according to its broad emission lines.
Basing upon the calibration of the measurements of the reverberation mapping methods,
the empirical relationships that are commonly used to estimate $L/L_{\rm E}$ and $M_{\rm{BH}}$ in
typical type I AGNs have been well established in recent years
(e.g., Kaspi et al., 2000, 2005; Peterson et al. 2004;
Vestergaard et al. 2002; Wu et al. 2004). Both velocity of broad H$\beta$ emission line
and luminosity at 5100\AA\ (or luminosity of broad H$\beta$ emission) emitted from the central AGN
are required in these calibrations,
which are however not feasible in our sample both because of the
dilution of AGN's luminosity
by the host's starlight and because of the non-detectable broad H$\beta$ emission in a majority of cases.
Greene \& Ho (2005) recently extended
the estimations of both bolometric luminosity and $M_{\rm{BH}}$ to broad H$\alpha$ line
\begin{equation}
L_{5100}=2.4\times10^{43}\bigg(\frac{L_{\mathrm{H}\alpha}}{10^{42}\ \mathrm{erg\ s^{-1}}}\bigg)^{0.86}\ \mathrm{erg\ s^{-1}}
\end{equation}
\begin{equation}
M_{\mathrm{BH}}=2\times10^{6}\bigg(\frac{L_{\mathrm{H\alpha}}}{10^{42}\ \mathrm{ergs\ s^{-1}}}\bigg)^{0.55}
\bigg(\frac{\mathrm{FWHM(H\alpha)}}{1000\ \mathrm{km\ s^{-1}}}\bigg)^{2.06}M_{\odot}
\end{equation}
The H$\alpha$ emission line provides much more robust estimations of both $L/L_{\rm E}$ and $M_{\rm{BH}}$
in the cases where the contamination caused by host and/or non-thermal jet radiation is present.
The relationships allow us to estimate both $L/L_{\rm E}$ and $M_{\rm{BH}}$ in
the composite AGNs. For each object,
the estimated values of $M_{\rm{BH}}$ and $L/L_{\rm E}$ is listed in Column (10) and
(11) in Table 1, respectively.

According to the unified model of AGN (Antonucci 1993), the existence of the broad emission lines
in AGN's spectrum depends
on the orientation of the axis of tours with respect to the line-of-sight of observer. In order to assess
how much the orientation effect affects $L_{\mathrm{H\alpha}}$ (and hence $L/L_{\mathrm{Edd}}$
and $M_{\mathrm{BH}}$)
that is estimated from the broad H$\alpha$ line, the $L_{[\mathrm{OIII}]}$ vs.
$L_{\mathrm{H\alpha}}$ diagram is shown in the \sl Panel D \rm in Figure 7.
A strong correlation is found between $L_{[\mathrm{OIII}]}$
and $L_{\mathrm{H\alpha}}$. [\ion{O}{3}] emission is isotropic in AGN because the line is
emitted only from the optically thin region that is at large enough distance from the
central black hole. The
isotropy of [\ion{O}{3}] emission is also supported by the correlation between [\ion{O}{3}] and
orientation-independent [\ion{O}{2}] emission (Kuraszkiewicz et al. 2000),
although the isotropy
has been questioned in some radio-loud AGNs (Baker 1997; di Serego Alighirei et al. 1997).
The tight correlation shown in the \sl Panel D \rm in Figure 7 therefore indicates that
the variation of $L_{\mathrm{H\alpha}}$
in our sample is not entirely due to the orientation effect,
and that $L_{\mathrm{H\alpha}}$ is a good measure of AGN's power in our sample.

A simple least-square linear fit within a wide range of luminosity gives an empirical relationship
between $L_{[\mathrm{OIII}]}$ and $L_{\mathrm{H\alpha}}$
\begin{equation}
\log\left(\frac{L_{[\mathrm{OIII}]}}{10^{40}\ \mathrm{erg\ s^{-1}}}\right)=(1.518\pm0.06)+
(1.078\pm0.07)\log\left(\frac{L_{\mathrm{H\alpha}}}{10^{42}\ \mathrm{erg\ s^{-1}}}\right)
\end{equation}
Incorporating with Eq. 3 yields a relationship between continuum luminosity $L_{5100}$ and
$L_{\mathrm{[OIII]}}$:
$L_{5100}/10^{42}\mathrm{erg\ s^{-1}}\approx1.5(L_{\mathrm{[OIII]}}/10^{40}\mathrm{erg\ s^{-1}})^{0.8}$.
Taking into account of the dynamical range of $L_{\mathrm{[OIII]}}$ from
$\sim10^{40}$ to $\sim10^{42}\ \rm{erg\ s^{-1}}$,
an average ratio is $L_{5100}/L_{\mathrm{[OIII]}}\sim100$.
Although the ratio is three times lower than the mean value derived
in Kauffmann et al. (2003a) and in Heckman et al. (2004), it
is reasonable because of the large scatter in the correlation between $L_{\mathrm{[OIII]}}$
and $L_{5100}$. There are two possible
explanations for the discrepancy of the $L_{5100}/L_{\mathrm{[OIII]}}$ ratio: 1) The luminosity
of the broad H$\alpha$ line (and hence $L/L_{\mathrm{Edd}}$ and $M_{\mathrm{BH}}$)
might be underestimated due to the obscuration caused by torus in the frame of
the unified model of AGN (Antonucci 1993).  It is noted that this systematic
shift would not seriously affect our final statistic trends, both because of the
tight correlation between $L_{\mathrm{H\alpha}}$ and  $L_{\mathrm{[OIII]}}$ and
because of the weak dependencies of $L/L_{\mathrm{Edd}}$ and $M_{\mathrm{BH}}$ on $L_{\mathrm{H\alpha}}$
(i.e., $L/L_{\mathrm{Edd}}\propto L_{\mathrm{H\alpha}}^{0.31}$ and
$M_{\mathrm{BH}}\propto L_{\mathrm{H\alpha}}^{0.55}$, see Eqs. 3 and 4).
2) In additional to the AGN's contribution, star formation can also contributes to the [\ion{O}{3}] luminosity.
In fact, our sample is mainly composed of the composites in which the narrow emission lines
are likely to be mainly ionized by hot stars (e.g., Kewley et al. 2006, and see below).

We estimate the upper limits of the underestimates of both $L/L_{\mathrm{Edd}}$ and $M_{\mathrm{BH}}$
as follows. Assuming the underestimate of $L_{\mathrm{H\alpha}}$ is entirely
caused by the obscuration of torus and by the reddening in BLR, the observed broad H$\alpha$ luminosity
$L'_{\mathrm{H\alpha}}$ could be written as
$L'_{\mathrm{H\alpha}}=f_{t}f_{e}L_{\mathrm{H\alpha}}$, where $L_{\mathrm{H\alpha}}$ is the intrinsic
H$\alpha$ luminosity, $f_{t}$ and $f_{e}$ represents the decreases of the line luminosity due to the
obscuration and reddening, respectively. If the three times lower in
the $L_{5100}/L_{\mathrm{[OIII]}}$ ratio is caused by the obscuration, $f_t$ is inferred to be $\sim0.28$
according to the relation $L_{5100}\propto L_{\mathrm{H\alpha}}^{0.86}$.
The factor $f_e$ could be estimated from the 19 objects for
which the broad H$\beta$ is present. Recent studies found the mean value of the broad-line Balmer decrements
(H$\alpha$/H$\beta$) in typical type I AGNs is significantly close to the value of 3.1
predicted by the case B recombination
(e.g., Greene \& Ho 2005; Vanden Berk et al. 2001; Dong et al. 2007). However, large
Balmer decrements (7-10) are observed in Seyfert 1.9-like galaxies (e.g., Osterbrock 1981).
Given the Balmer decrement in the case B recombination,
only upper limit of $f_e$ could be obtained in our sample,
both because of the orientation effect and because of the fact that higher
density and/or higher ionization parameter is required to produce H$\beta$ than H$\alpha$.
The mean value of the broad-line Balmer decrement
$\langle$H$\alpha$/H$\beta\rangle_\mathrm{B}$ is 5.97 in current sample.
Excluding the reddening inferred from
the narrow-line Balmer decrement $\langle$H$\alpha$/H$\beta\rangle_\mathrm{N} =4.13$,
we obtain $f_e\sim0.47$.
Finally, the upper limit of the underestimate of $L/L_{\mathrm{Edd}}$ is estimated to be $\sim50$\%,
and that of $M_{\mathrm{BH}}$ to be $\sim70$\%.

\subsection{Deriving age of stellar population of host}

Both 4000\AA\ break ($\rm D_n(4000)$) and equivalent width (EW) of the H$\delta$ absorption (H$\delta_A$)
are widely used in AGN's hosts as indicators of the ages of their stellar populations (e.g., Heckman et al. 2004;
Kauffmann et al. 2003a; Kewley et al. 2006).
Incorporating with the large high-resolution stellar spectra libraries,
the latest theoretical population synthesis models of BC03 indicate that both $\rm D_n(4000)$
and H$\delta_{\rm{A}}$ are reliable age indicators until a few Gyrs after
the burst (e.g., Figure 2 in Kauffmann et al. 2003c),
although both indices are sensitive to metallicity in the very old stellar populations.

The 4000\AA\ break is a strong discontinuity in the optical spectrum of a galaxy.
The break is mainly caused by the absorption features of ionized metal lines blueward of 4000\AA.
The popular definition of $\rm D_n(4000)$ is (Balogh et al. 1999; Bruzual 1983)
\begin{equation}
D_n(4000)=\frac{\int_{4000}^{4100}f_\lambda d\lambda}{\int_{3850}^{3950}f_\lambda d\lambda}
\end{equation}
Gorgas et al. (1999) calibrated $\rm D_n(4000)$ in terms of the stellar atmospheric parameters: effective
temperature, metallicity and surface gravity. In a young stellar population,
small $\rm D_n(4000)$ arises from the
low opacity of metal since it is multiply thermal ionized in hot stars. On the contrary, large
$\rm D_n(4000)$ appears in an old stellar population in which the metal opacity is expected to be high.

The strength of the H$\delta$ absorption is an age indicator when the stellar populations are in the post-starburst phase,
i.e. 0.1-1 Gyrs after the burst when massive stars (O and B stars) terminated their evolutions
(Worthey \& Ottaviani 1997; Gonzalez Delgado et al. 1999).
The index H$\delta_{\rm{A}}$ measures the EW of absorption in A type stars, and is defined
as (Worthey \& Ottaviani 1997)
\begin{equation}
\mathrm{H}\delta_A=(4122.25-4083.50)(1-\frac{F_I}{F_C})
\end{equation}
where $F_I$ is the flux within the $\lambda\lambda4083.50-4122.25$ feature bandpass, and $F_C$ the flux of
the pseudo-continuum within two defined bandpasses: blue $\lambda\lambda4041.60-4079.75$ and red
$\lambda\lambda4128.50-4161.00$.
As the case of $\rm D_n(4000)$, similar calibration has been carried out by Worthey \& Ottaviani (1997).

We measure both indices in each starlight spectrum.
The starlight spectrum is separated from the raw data as described in Section 2.2.
For each object, the measured values of $\rm D_n(4000)$ and H$\delta_{\rm{A}}$
is listed in Column (7) and (8) in Table 1, respectively.
The comparison between the values of $D_n(4000)$ obtained in this paper and
that derived by MPA/JHU group is shown in the left panel in Figure 2.
The same comparison but for H$\delta_A$ is shown in the right panel.
The left panel shows a significant correlation (with slope $\approx1$)
between the two sets of $D_n(4000)$ values,
and the right panel a less significant, but still tight one for H$\delta_A$. The comparisons clearly
indicate that our measurements are highly consistent with those derived by MPA/JHU group, although
our H$\delta_A$ appears to be slightly smaller than those of MPA/JHU group at left-bottom corner.
This small discrepancy could be ignored in current study since
we only focus on the statistic properties of the composite AGNs.
A large discrepancy is only
found in one object SDSS\,J085338.27+033246.1. The value of H$\delta_A$ derived in this paper is
7.22 which is much larger than H$\delta_A$=2.78 obtained by MPA/JHU group.
The large discrepancy is due to the difference in the spectra components adopted to
model both continuum and absorption features. In addition to the seven eigenspectra, a
powerlaw continuum attributed to AGN is involved in the spectra modeling for the object,
which reduces its starlight continuum level, and enhances
its measured H$\delta_A$ consequently.

\section{Classification of the Composite AGNs}

In this section, we first examine the distributions of the physical properties derived in the last
section on the traditional BPT diagrams, and then perform a very fine classification
on the composite AGNs. Their host's properties are finally explored basing upon the classification.

The BPT diagrams are used as a powerful tool to diagnose dominant energy source in
emission-line galaxies. Seyferts, LINERs, and star-forming galaxies can be separated according to
their emission-line ratios. With respect to star-forming galaxies,
AGNs have larger ratios in [\ion{O}{3}]/H$\beta$,
[\ion{N}{2}]/H$\alpha$, [\ion{S}{2}]/H$\alpha$, and [\ion{O}{1}]/H$\alpha$,
because of their stronger and harder ionizing field.
AGNs therefore concentrate in the up-right corner in each BPT diagram.
The theoretical curves discriminating between AGNs and star-forming
galaxies in the BPT diagrams were recently proposed by Kewley et al. (2001).
There are many optical classification schemes to segregate LINERs from other emission-line galaxies
(e.g., Heckman et al. 1980, Veilleux \& Osterbrock 1987; Ho et al. 1997). In this study, we adopt the
new empirical classification scheme recently proposed by Kewley et al. (2006).
The classification scheme was obtained by analyzing a large sample
of emission-line galaxies selected from SDSS survey. These new classification line is very
similar to that of Heckman et al (1980) in both slope and intersection.

The BPT diagrams are displayed in Figure 3 for the composite AGNs.
The three panels correspond to the three different methods for classifying
galaxies using two pairs of line ratios. In each
panel, the theoretical classification line is shown by a solid line,
and the empirical classification line suggested by Kauffmann et al. (2003a) by a dashed line.
The classification lines used to separate Seyferts from LINERs are shown by dot-dashed lines.
The size of each point is scaled to be proportional to the corresponding value of $D_n(4000)$ simply
because of the high level of accuracy of this parameter (See Section 3.2 for details).

At first, almost
all the composite AGNs are located between the theoretical and empirical lines in the
[\ion{O}{3}]/H$\beta$ vs. [\ion{N}{2}]/H$\alpha$ diagram, which indicates that
our AGN-host separations and emission-line measurements are consistent with those of MPA/JHU groups.
In addition, it is interesting that
the objects occupy a narrow vertical stripe with nearly constant line ratio
[\ion{N}{2}]/H$\alpha$. The examination of this phenomenon is out of the scope of this paper, and
will be studied in subsequent works. Moreover, no clear trend is found along this
stripe for $D_n(4000)$, although the objects with larger $D_n(4000)$ are only found at the top end of the stripe.
In contrast, clear trends are found in the other two panels.
The objects located in the occupations of LINERs generally have larger
$D_n(4000)$ with respect to that located below the theoretical demarcation lines,
which implies that the composite AGNs at least could be divided into two sub-groups:
\ion{H}{2} region-like ionization + young stellar population and LINER-like ionization + old stellar population
\footnote{Similar phenomena and indications could be found in the
case of $L/L_{\mathrm{Edd}}$ and luminosity of [\ion{O}{3}] instead of $D_n(4000)$, with,
however, less significance.}.

Since the composite AGNs are initially selected only from the [\ion{N}{2}]/H$\alpha$ vs. [\ion{O}{3}]/H$\beta$ diagram,
a more fine classification on the objects is performed as similar as done in Kewley et al. (2006).
We separate our sample into three groups according to
the classification scheme described as follows:
a) The \ion{H}{2} composite AGNs are galaxies that lie below the theoretical demarcation curves
in all the [\ion{N}{2}]/H$\alpha$, [\ion{S}{2}]/H$\alpha$ and
[\ion{O}{1}]/H$\alpha$ diagrams. These galaxies correspond to the ``composite galaxies'' in Kewley et al. (2006);
b) In contrast, the Seyfert+LINERs are galaxies that lie above the theoretical demarcation
curves in all the three diagrams; c) The remaining are multiply-classified galaxies that are classified
differently in the three diagrams. These galaxies correspond to the case b ``ambiguous galaxies'' in
Kewley et al. (2006).
The classification is shown in the three BPT diagrams in Figure 4. The demarcation lines are
shown in the figure by the same symbols used in Figure 3.
Only 74 objects for which all the line ratios are
determined are plotted. The \ion{H}{2} composite AGNs, Seyfert+LINERs and multiply-classified
galaxies are shown by red, blue, and black points, respectively.

Our detailed classification indicates
that, among the 74 objects, six ones can be classified as Seyfert+LINERs rather than composites.
The fraction of these six objects is so small that their existence has negligible effect on
the overall properties of our sample. Totally,
there are 45 \ion{H}{2} composite AGNs and 23 multiply-classified galaxies. The comparison in
the age of stellar population between
the \ion{H}{2} composite AGNs and multiply-classified galaxies is illustrated in Figure 5.
The dashed line shows distribution of $D_n(4000)$ for the \ion{H}{2} composite AGNs, and
the solid line for the multiply-classified galaxies.
A narrow distribution with median value $D_n(4000)$=1.30 and standard deviation of 0.11
is found in the \ion{H}{2} composite AGNs.
On the contrary, a very wide and possibly
bimodal distribution is identified in the multiply-classified galaxies.

\section{Star Formation History and Evolution of AGN}

We plot our sample of the composite AGNs in the $D_n(4000)$/H$\delta_A$ plane in Figure 6.
The dashed line shows the evolution of the standard SSP model with solar metallicity, and the dot-dashed line
the evolution of the model with continuous star formation history (i.e., an exponentially declining SFR
$\psi(t)\propto e^{-t/(4\rm Gyr)}$) at solar metallicity (BC03). The $D_n(4000)$/H$\delta_A$
plane is a useful diagnosis to determine whether a galaxy has been forming stars continuously or
in a recent burst. Galaxies with continuous star formation history occupy a narrow locus
in the plane. A recent burst ended in 0.1-1 Gyrs ago results in a significant
displacement away from the locus (i.e., be larger in H$\delta_A$), because the optical spectrum
is dominated by the emission of A type stars. Even though the small discrepancy in
H$\delta_A$ is found between our measures and those in MPA/JHU groups, the composite AGNs
fall close to the model with continuous star formation history, except the
object SDSS\,J085338.27+033246.1 with significant strong H$\delta$ absorption H$\delta_A$=7.22 (see
Section 3.2 for explanation).

Our sample allows us to investigate the issue of co-evolution of AGN and star formation
in a new approach. Our study differs from the previous investigations in: 1)
it should be emphasized that each composite AGN
certainly contains a Seyfert 1 nucleus in its center, because of the existence of the broad H$\alpha$
emission. In fact, there might be two populations of Seyfert 2 galaxies: with and without a hidden
Seyfert 1 nucleus. The spectropolarimetric surveys and X-ray observations indicate that the
hidden Seyfert 1 nuclei are indeed not present in a fraction (10\%-50\%) of Seyfert 2 galaxies
(e.g., Tran 2001, 2003; Moran et al. 2000, 2001, 2007; Panessa \& Bassini 2002; Gallo et al. 2006;
Pappa et al. 2001); 2) In addition, the basic properties of AGN (i.e., $L/L_{\mathrm{Edd}}$ and $\rm M_{BH}$)
are directly obtained from the
nuclear emission rather than the tracers with large uncertainties. For instance,
the luminosity of [\ion{O}{3}] emission line is statistically used as a tracer of AGN's bolometric
luminosity in Seyfert 2 galaxies by many authors (e.g., Kauffmann et al. 2003a; Heckman et al. 2004; Kewley et al. 2006).

It was suspected that $L/L_{\mathrm{Edd}}$ might be of significance of ``age'' of AGN
(e.g., Grupe et al. 1999; Grupe 2004; Wang et al. 2006). In order to examine this point, the
$L/L_{\mathrm{Edd}}$ vs. $D_n(4000)$ diagram is presented in the \sl Panel A \rm in Figure 7.
The diagram indicates that the $L/L_{\mathrm{Edd}}$ estimated from the broad H$\alpha$
emission varies as a function of the age of the host's stellar population.
The black holes with high
specific accretion rates tend to reside in the hosts with relatively young stellar populations.
This dependence therefore strongly supports the evolutionary
scenario of AGN in which AGNs with high $L/L_{\mathrm{Edd}}$ evolve to ones with low $L/L_{\mathrm{Edd}}$
(Grupe 2004; Wang et al. 2006; Wang \& Wei 2006).
The \sl Panel C \rm in Figure 7 shows a variation of the $L/L_{\mathrm{Edd}}$ with the line ratio
[\ion{O}{1}]/H$\alpha$ which is sensitive to the hardness of the ionizing field. As shown in
the diagram, a soft ionizing spectrum tends to be associated with a relatively high
$L/L_{\mathrm{Edd}}$ on average, and a hard spectrum to be associated with an extremely low $L/L_{\mathrm{Edd}}$ only.
This trend is consistent with the strong correlation identified between the hardness of hard/soft X-ray emission
and $L/L_{\mathrm{Edd}}$ for both local and distant AGNs (Grupe 2004; Bian 2005; Shemmer et al. 2006).
The lowest value of the $L/L_{\mathrm{Edd}}$ among our sample is 0.003, which is consistent with
the prediction of the advection dominated accretion flow (ADAF) model.
The model can only produce sub-Eddington accretions ($L/L_{\mathrm{Edd}}\leq0.01-0.1$).
Because a large fraction of the energy released is advected
with hot gas into black hole rather than radiated by emission, the gas is extremely hot, and
produces a harder spectrum in a low state
(e.g., Abramowicz et al. 1995; Narayan \& Yi 1995; Esin et al. 1997; Meyer et al. 2000;
Ferreira et al. 2006).

The another candidate of ``age'' of AGN is the accretion rate onto a black hole
(i.e., luminosity released by accretion). Relatively good relations are identified between
$L_{\mathrm H\alpha}$ and $D_n(4000)$, and between $L_{\mathrm H\alpha}$ and [\ion{O}{1}]/H$\alpha$.
The two relations are presented in Figure 8.
Noted that both relations degrade at the right-bottom end. This phenomenon does not appear
when $L/L_{\mathrm{Edd}}$ is used in place of $L_{\mathrm H\alpha}$ as shown in the \sl Panel A \rm
and \sl C \rm in Figure 7, which indicates that the $L/L_{\mathrm{Edd}}$ is the major physical driver
in the evolution of AGN.

A significant correlation between $D_n(4000)$ and the line ratio [\ion{O}{1}]/H$\alpha$ is first shown in
the \sl Panel B \rm in figure 7. The tight correlation implies that AGN evolves from the high state
with a soft spectrum to the low state with a hard spectrum. The line ratio [\ion{O}{1}]/H$\alpha$ is
suggested to be used as a potential indicator of the age of stellar population.

\section{Discussion}

In addition to the AGN's contribution, other main mechanisms of the broad Balmer emission lines
are the stellar envelopes of Wolf-Rayet stars and OB associations associated multiple SN events
(Izotov et al. 2007 and references therein). Both mechanisms can produce the broad emission lines
with gaseous velocity as large as $\stackrel{>}{\scriptstyle\sim} 1000\ \mathrm{km\ s^{-1}}$.
The luminosity of the broad H$\alpha$ emission line ranges from $10^{40}$ to $10^{42}\ \mathrm{ergs\ s^{-1}}$
in the composite AGNs, which greatly exceeds the typical luminosity of stellar winds from WR or LBV stars
($10^{36}$-$10^{40}\  \mathrm{ergs\ s^{-1}}$). Moreover, the broad H$\alpha$ emission in
the composite AGNs is unlikely caused by the SN events, because of the existence of the
underlying relatively old stellar populations in the spectra.

\subsection{Aperture effect}
The redshift ranges roughly from 0.05 to 0.15 in our sample.
In this range of distance, the typical projected fiber diameter ($3$\symbol{125}) is from $\sim3$kpc to $\sim8$kpc.
When a spectrum is taken in the SDSS spectroscopic survey, the fixed fiber aperture
would result in more contamination of the light emitted from the outer part of host galaxy in
distant galaxies than in nearby galaxies.
First, the effect leads the LINER fraction decreases with redshift because the weak
emission lines are more difficult to be detected in distant galaxies than in nearby galaxies.
Such effect exists in our sample, and is shown in the left panel in Figure 9.
The fraction of the LINER-like galaxies with an old stellar population drops significantly when redshift
$z>0.1$, which is consistent with Kauffmann et al. (2003a) and Kewley et al. (2006).
Secondly, for a fixed aperture, as the contamination increases
with redshift, the value of $D_n(4000)$ would increase with redshift if a circumnuclear starburst occurs.
If so, one could see a lower boundary with positive slope in the $D_n(4000)$ vs. z diagram,
which is, however, not present in the diagram. The aperture effect can not explain why no
young stellar populations ($D_{4000}\leq1.2-1.4$) is deteced in the low-redshift galaxies.
The negative trend shown in the diagram
could be explained by the consequence of the increase of $L_{\mathrm{H\alpha}}$ with the redshift.
Moreover, the lower limit of redshift of $z\sim0.05$ ensures that the spectra
taken within $3$\symbol{125} aperture could approximately reflect the global properties (Kewley et al. 2005).
The right panel in Figure 9 shows the broad-to-narrow H$\alpha$ line ratio as a function of the redshift.
As shown by the diagram, in general, our criterion for the existence of the broad H$\alpha$ line
is free of the aperture effect.

\subsection{Star formation rate}
In section 4, we have built a sub-sample of the \ion{H}{2} composite AGNs which
are below the theoretical demarcation lines in the three BPT diagrams.
The left panel in Figure 10
shows the [\ion{O}{3}]/[\ion{O}{2}] vs. [\ion{O}{1}]/H$\alpha$ diagram for the
\ion{H}{2} composite AGNs.
The distribution shows that the \ion{H}{2} composite AGNs
lie within the occupation of the star-forming galaxies, which
additionally suggests that hot stars might be main contributor of the ionizing field for the narrow emission lines.
The luminosity of [\ion{O}{2}] emission line is
an indicator of current star formation rate (SFR) for \ion{H}{2} regions as suggested by
Kennicutt (1998) and Kewley et al. (2004). However, the origin of the [\ion{O}{2}] emission is complex in the
\ion{H}{2} composite
AGNs, since both AGN and starburst can contribute to [\ion{O}{2}] (e.g., Kim et al. 2006; Yan et al. 2006).
In order to estimate the contribution by AGN, the right panel in Figure 10 shows the luminosity of [\ion{O}{3}]
as a function of the line ration [\ion{O}{2}]/[\ion{O}{3}]. The solid line is the least-squares bisector
regression for Type I AGNs given by Kim et al. (2006). The \ion{H}{2} composite AGNs clearly show enhanced
[\ion{O}{2}]/[\ion{O}{3}] ratios. Similar as done in Kim et al. (2006),
we estimate the [\ion{O}{2}] luminosities emitted from \ion{H}{2} region by
assuming the enhanced [\ion{O}{2}]/[\ion{O}{3}] ratios are caused by star formations.
The mean value of the [\ion{O}{2}]/[\ion{O}{3}] ratio is $\sim1.0$ for the \ion{H}{2} composite AGNs.
Given the average [\ion{O}{3}] luminosity $L_{[\mathrm{OIII]}}\sim3\times10^{41} \mathrm{\ erg\ s^{-1}}$,
the [\ion{O}{2}]/[\ion{O}{3}] ratio predicted by the regression is only $\sim0.25$. Then, on average,
$\sim$75 percent of the [\ion{O}{2}] emission is estimated to come from star formation in the \ion{H}{2} composite AGNs.

We use the recent calibration (Kewley et al. 2004)
\begin{equation}
\mathrm{SFR}=7.9\times\frac{L_{\mathrm{[OII],42}}}{16.73-1.75[\log(\mathrm{O/H})+12]}M_{\odot}\rm{yr^{-1}}
\end{equation}
to estimate current SFRs, where $L_\mathrm{[OII],42}$ is the luminosity of [OII]
emission in units of $10^{42}\ \rm{ergs\ s^{-1}}$. The metallicity is fixed to be $\log (\rm{O/H})+12=9.2$,
i.e., twice of the solar, as suggested by Ho (2005b).

A correlation between the current SFRs and $L/L_{\mathrm{Edd}}$ is
shown in the \sl left \rm panel in Figure 11. The relatively large scatter could be caused by
both the underestimation of $L_{\mathrm{H\alpha}}$ and the contribution of AGN to $L_{\mathrm{[OII]}}$.
The correlation indicates that a more powerful AGN tends to be
associated with a more intense current star formation. The relation is reasonable because
the strengths of both AGN and star-forming activities depend on the consumption of gas.
Heckman et al. (2004) found a roughly constant ratio of the star formation rate to the accretion rate onto black
hole in the SDSS survey. The dependence of the SFR on the
luminosity of AGN is further supported by the correlation,
between the X-ray luminosity in 0.5-4.5 keV band and CO luminosity, found in Seyfert 1 galaxies
and in quasars (Yamada 1994).
The \sl middle \rm and \sl right \rm panels
in Figure 11 shows variations of the SFRs with $D_n(4000)$ and H$\delta_A$, respectively.
As shown by the diagrams, the current SFR decreases when stellar populations passively evolves,
because of the depletion of the gas. Both relations are consistent with Figure 6 in which a continuous star
formation history with an exponentially decaying SFR is required to match the observations to the modeled evolution locus.
In addition, the [\ion{O}{2}] inferred SFRs are typically below 10$\rm{M_{\odot}yr^{-1}}$.
The suppressed current SFRs are identified in a number of luminous AGNs by several authors (e.g., Ho 2005b;
Kim et al. 2006; Wang \& Wei 2006; Martin et al. 2007).

\subsection{Role of $L/L_{\mathrm{Edd}}$ in AGN evolution}

Similar as the scenario proposed by Weedman (1983),
Wang \& Wei (2006) recently proposed a scenario of co-evolution of AGN and star formation in its host
basing upon the differential growth of the black hole mass and bulge mass.
In the scenario, after a starburst occurs first and triggers accretion onto the central black hole latter,
AGN then shines with decreasing $L/L_{\mathrm{Edd}}$ and insignificant increasing black hole mass,
as the young stellar population ages and fades. The scenario is supported
by the prediction of the model developed by Kawakatu et al. (2003). The model predicts that the AGN-dominated
phase follows the host-dominated phase in a time scale of a few $10^8$ years. This change in phase is
also predicted in the more solid physical models developed by Granato et al. (2001, 2004).

In addition to the models, there were accumulating observations supporting
that $L/L_{\mathrm{Edd}}$
plays an important role in the evolution of galaxy in recent years. Mathur (2000) suggested
that NLS1 might be young AGN with high $L/L_{\mathrm{Edd}}$ and small black hole mass. The young stellar
populations within post-starburst phase are found to be associated with NLS1s in a few cases
(e.g., Wang et al. 2004; Wang \& Wei 2006; Zhou et al. 2005).
The evolutionary significance of $L/L_{\mathrm{Edd}}$
was suspected by Grupe (2004). In two years latter,
Wang et al. (2006) extended the well-documented E1 space into the infrared color
$\alpha(60,25)$ that is an indicator of relative strength of star formation with respect to AGN.
The infrared-color dominated E1 space
implies a relation between the $L/L_{\mathrm{Edd}}$ and star formation history. With the advent of
SDSS, Kewley et al. (2006) and Wild et al. (2007, and references therein) revealed a relationship
between the star formation history and $L_{\mathrm{[OIII]}}/\sigma_*^4$ in type II AGNs .

Taking advantage of the composite AGNs,
the evolutionary significance of $L/L_{\mathrm{Edd}}$ (or the E1 space) is obviously
demonstrated by the clear dependence of $L/L_{\mathrm{Edd}}$ on $D_n(4000)$.
In addition, the evolutionary significance of the E1 space is also implied by the tight
correlation between [\ion{O}{1}]/H$\alpha$ and $D_n(4000)$. The line ratio [\ion{O}{1}]/H$\alpha$ is known
to be sensitive to the hardness of the ionizing radiation field. The X-ray photon index is an important
component in the E1 space
(e.g., Wang et al. 1996; Sulentic et al. 2000; Xu et al. 2003; Grupe et al. 1999; Grupe 2004).
In summary, by adopting the new approach, our results are not only consistent with the previous studies, but also
strongly support the argument that AGN evolves along the E1 sequence from
high $L/L_{\mathrm{Edd}}$ end to low
$L/L_{\mathrm{Edd}}$ end as stars evolve from young to old stellar population.

\subsection{``Young'' and ``Old'' composite AGNs}

Our results indicate that $\sim$80\% of the composite AGNs are high $L/L_{\mathrm{Edd}}$
emitters associated with a post-starburst (with age less than
1\,Gyr, corresponding to a value $D_n(4000)\leq1.5$).
In the context of the evolutionary scenario, the composite AGNs associated with both
high $L/L_{\mathrm{Edd}}$ and young stellar population might be ``Young'' AGNs.
The spatial connection between the post-starburst and AGN is supported by the spatial
resolved spectroscopy in the three AGNs selected from SDSS (Goto 2006), and by recent photometric observations
with spatial resolution of 0.1-0.2\symbol{125} (Davies et al. 2007). In terms of the equivalent width
of the narrow H$\alpha$ emission, Westoby et al. (2007) found the composite galaxies peak in the valley
between the red and blue sequence (e.g., Strateva et al. 2001; Baldry et al. 2006), which suggests
the composite galaxies are in the transition state. By calculating a series of photoionization models,
Stasinska et al. (2006) suggested that the right wing in the BPT
diagram could be explained by the variation of the balance between massive stars and AGN powers along the wing.
In the BPT diagram, a galaxy moves toward right-up direction from the occupation of
the composite AGNs as the fraction of AGN's power increases.

On the other hand, although the multiply-classified galaxies mostly lie
below the theoretical demarcation curve in the [\ion{N}{2}]/H$\alpha$ vs. [\ion{O}{3}]/H$\beta$ diagram,
a fraction of them are located in the LINER region in the [\ion{O}{1}]/H$\alpha$ vs. [\ion{O}{3}]/H$\beta$ diagram.
The harder ionizing field and smaller AGN's power could be inferred in these multiply-classified galaxies,
because of the saturation of [\ion{N}{2}]/H$\alpha$ with increasing nebular metallicity in \ion{H}{2} region
(e.g., Kewley \& Dopita 2002; Pettini \& Pagel 2004). We suspect that the multiply-classified galaxies
associated with very old stellar populations are at the end of life of AGNs.

\section{Summary}
The goal of this paper is to study the evolutionary connection between AGN and its host galaxy in
partially obscured AGNs (i.e., Seyfert 1.8 and 1.9 galaxies). Not as the previous studies,
we can derive the properties of both AGN and its host together in these objects
by direct measurements rather than some tracers.
Taking advantage of the SDSS spectroscopic survey with large number of spectra of
emission-line galaxies,
we searched for Seyfert 1.8 and 1.9 galaxies (composite AGNs) from the emission-line
galaxies (MPA/JHU catalogs), which are located between the empirical and theoretical classification lines
in the [\ion{O}{3}]/H$\beta$ vs. [\ion{N}{2}]/H$\alpha$ diagram.
The analysis of these 85 composite AGNs allows us to draw results as follows:

\begin{enumerate}

\item The distribution of the composite AGNs on the [\ion{N}{2}]/H$\alpha$
vs. [\ion{O}{3}]/H$\beta$ diagram shows a
narrow vertical stripe with nearly constant line ratio [\ion{N}{2}]/H$\alpha$.
In contrast, the composite AGNs can be separated in the other two BPT diagrams, in which the
objects associated with the very old stellar populations are located in the occupation of LINERs.

\item The BPT diagrams allow us to further divide the 68 composite AGNs into two subgroups:
45 HII composite AGNs which are below the theoretical
separation lines in all the three
BPT diagrams; 23 multiply-classified galaxies which are classified differently
in the three BPT diagrams. All the HII composite AGNs
are associated with the young stellar populations ($<$1Gyrs),
while a very wide and possibly bimodal distribution of
the age of the stellar population is identified in
the multiply-classified galaxies. Totally, $\sim$80\% of the composite AGNs
are associated the young stellar populations within post-starburst phase ($D_n(4000)\leq1.5$). We suspect that
the composite AGNs associated with the very old stellar populations are at the end of life of AGNs.

\item The broad H$\alpha$ emission inferred $L/L_{\mathrm{Edd}}$ forms a smooth sequence
with the age of the stellar population of AGN's host as assessed by $D_n(4000)$,
which confirms the previous studies using other
tracers (e.g., $L_{\mathrm{[OIII]}}/\sigma_*^4$, infrared color $\alpha(60,25)$). The
sequence degrades when $L/L_{\mathrm{Edd}}$ is replaced by $L_{\mathrm{H\alpha}}$.
In addition, the $L/L_{\mathrm{Edd}}$ roughly depends on the line ratio
[\ion{O}{1}]/H$\alpha$ which is sensitive to the hardness of the ionizing field.
Both theoretical population synthesis models and
[\ion{O}{2}] inferred SFR $\sim1-10M_{\odot}\mathrm{yr^{-1}}$ show a suppressed
current star formation in the composite AGNs. These results strongly support the evolutionary
scenario suggested in Wang \& Wei (2006).

\item We first find a significant correlation between the line ratio [\ion{O}{1}]/H$\alpha$
and $D_n(4000)$ in the composite AGNs. If it is confirmed, the
line ratio could be used to trace the age of stellar population in type I AGNs.

\end{enumerate}

\acknowledgments

We would like to thank the anonymous referee for his/her very useful comments and suggestions.
We thank D. W. Xu, C. N. Hao, Y. C. Liang and C. Cao for valuable discussions. We are
grateful to Todd A. Boroson and Richard F. Green for providing us with the
\ion{Fe}{2} template. This work was supported by the National Science Foundation
of China (grants 10503005 and 10473013). The SDSS achieve data are created and
distributed by the Alfred P. Sloan Foundation.


\section{Introduction}
\label{intro} Neutron stars can be included in the class of compact
stars which are stars with high mass and small radius and their
density will be very high. So we need to take into account effects
from general relativity, like the curvature of space--time in the
study of neutron stars. Baade and Zwicky were the first to propose
the idea of neutron stars (Baade and Zwicky, 1934),
 pointing out that they have very high density and small
radius. They also made a very important suggestion that neutron
stars would be formed in supernova explosions. We know that there
are two forces acting on a star, one of them is gravitation and the
second one arises from the pressure. The opposing force of
gravitation will be thermal pressure for ordinary stars and
degeneracy pressure for white dwarfs and neutron stars. Thus we can
obtain two coupled differential equations for mass and pressure
which tell their change with the radius of the star.
Tolman-Oppenheimer-Volkoff (TOV) equation (Tolman, 1939; Volkoff and
Oppenheimer, 1939)
 derived from the Einstein's field equation gives the pressure gradient
 of neutron stars. Numerical model of a neutron star with a
relativistic mean field theory was given by Walecka (Walecka and
Chin, 1974; Walecka, 1975).

Our Universe can be represented by SdS metric, since the recent
astronomical observations show that the expansion of our universe is
accelerating (Guzzo et al., 2008; Perlmutter et al., 1999; Riess et
al., 1998; Weinberg, 1972).
Thus Schwarzschild-de Sitter metric described with a positive
cosmological constant and having two horizons is considered to be a
model of this universe. In such a case, the use of TOV equation,
derived from the Einstein's field equation in the case of
accelerating universe, will be more apt for the theoretical
constructions of stellar models. Thus we will get a modified TOV
equation for pressure, for calculating the maximum mass of a neutron
star. Anisotropic models for compact self gravitating objects
have been studied extensively (Corchero, 2001; Dev and Gleiser,
2003; Herrera and Santos, 1997; Ivanov, 2002; Mak and Harko, 2003).
Our aim is to establish a modified TOV equation which is the prime
structure equation for compact stars.

In \S2 we will show the derivation of the new equation to replace
the TOV equation, starting with the standard metric for a
spherically symmetric star and the consecutive use of Einstein's
field equation in the case of accelerating universe. Then we will
show the numerical calculation of maximum mass of Neutron stars by
plotting the mass-radius relation.
\section{Derivation of the modified TOV equation}
\label{sec:1} We have the standard metric (Weinberg, 1972)
for a spherically symmetric star,
\begin{equation}
ds^{2}=-U(r)dt^{2}+V(r)dr^{2}+r^{2}d\theta ^{2}+r^{2}\sin ^{2}\theta
d\phi ^{2},  \label{1}
\end{equation}
where,
\begin{equation}
V(r)=\frac{1}{U(r)}.  \label{2}
\end{equation}

Ricci tensor is given by,
\begin{equation}
R_{\mu k}=\frac{\partial \Gamma _{\mu \lambda }^{\lambda }}{\partial
x^{k}}-\frac{\partial \Gamma _{\mu k}^{\lambda }}{\partial
x^{\lambda }}+\Gamma _{\mu \lambda }^{\eta }\Gamma _{k\eta
}^{\lambda }-\Gamma _{\mu k}^{\eta }\Gamma _{\lambda \eta }^{\lambda
}.  \label{3}
\end{equation}
Substituting for the Christoffel symbols, we will get the equations
as,
\begin{equation}
R_{rr}=-\frac{U^{\prime \prime }}{2U}-\frac{U^{\prime }}{4U}\left(
\frac{V^{\prime }}{V}+\frac{U^{\prime }}{U}\right)
-\frac{1}{r}\frac{V^{\prime }}{V},  \label{4}
\end{equation}
\begin{equation}
R_{\theta \theta }=-1+\frac{r}{2V}\left( -\frac{V^{\prime
}}{V}+\frac{U^{\prime }}{U}\right) +\frac{1}{V},  \label{5}
\end{equation}
\begin{equation}
R_{\phi \phi }\approx R_{\theta \theta },  \label{6}
\end{equation}
\begin{equation}
R_{tt}=-\frac{U^{\prime \prime }}{2V}+\frac{U^{\prime }}{4V}\left(
\frac{V^{\prime }}{V}+\frac{U^{\prime }}{U}\right) -\frac{U^{\prime
}}{rV}. \label{7}
\end{equation}
The result $R_{\phi \phi }\approx R_{\theta \theta }$ is the
consequence of the rotational invariance of the metric.

Einstein's field equation in the case of accelerating universe
(Weinberg, 1972)
is given by,
\begin{equation}
G_{\mu \nu }=-8\pi GT_{\mu \nu }+\Lambda g_{\mu \nu },  \label{8}
\end{equation}
where
\begin{equation}
G_{\mu \nu }=R_{\mu \nu }-\frac{1}{2}g_{\mu \nu }R.  \label{9}
\end{equation}
Substituting Eq. (\ref{9}) in Eq. (\ref{8}), we get
\begin{equation}
R_{\mu \nu }=-8\pi G\left({T}_{\mu \nu}-\frac{1}{2}g_{\mu \nu
}{T}_{\lambda }^{\lambda }-\frac{\Lambda }{8\pi G}g_{\mu \nu
}\right).   \label{10}
\end{equation}
From hydrodynamics, we have the relation for energy-momentum tensor
in the case of models with cosmological constant,
\begin{equation}
\tilde{T}_{\mu \nu }=\tilde{p}g_{\mu \nu }+\left(
\tilde{p}+\tilde{\epsilon }\right) u_{\mu }u_{\nu }, \label{11}
\end{equation}
where $\tilde{p}$-the modified pressure, $\tilde{\epsilon } $-the
modified total energy density and $u_{\mu }$-the velocity four
vector. Here
\begin{equation}
\tilde{p}=p-\frac{\Lambda }{8\pi G} \ \ \ \ \ \ \ \ \tilde{\epsilon
}=\epsilon +\frac{\Lambda }{8\pi G},  \label{12}
\end{equation}

Also the velocity vector is defined so that,
\begin{equation}
g^{\mu \nu }u_{\mu }u_{\nu }=-1.  \label{13}
\end{equation}%
Assuming fluid at rest, we have
\begin{equation}
u_{r}=u_{\theta }=u_{\phi }=0,  \label{14}
\end{equation}
\begin{equation}
u_{t}=-(-g^{tt})^{-\frac{1}{2}}=-\sqrt{U(r)}.  \label{15}
\end{equation}

Our assumptions of time independence and spherical symmetry imply
that $p$ and $\epsilon $ are functions only of the radical
co-ordinate $r$.

Using Eq. (\ref{10}), we will get
\begin{eqnarray}
R_{rr}=-\frac{U^{\prime \prime }}{2U}-\frac{U^{\prime }}{4U}\left(
\frac{V^{\prime }}{V}+\frac{U^{\prime }}{U}\right)
-\frac{1}{r}\frac{V^{\prime }}{V} \nonumber \\ =-4\pi G(\epsilon
-p)V-\Lambda V, \label{16}
\end{eqnarray}
\begin{eqnarray}
R_{\theta \theta }=-1+\frac{r}{2V}\left( -\frac{V^{\prime
}}{V}+\frac{U^{\prime }}{U}\right) +\frac{1}{V} \nonumber \\ =-4\pi
G(\epsilon -p)r^{2}-\Lambda r^{2}, \label{17}
\end{eqnarray}
\begin{eqnarray}
R_{tt}=-\frac{U^{\prime \prime }}{2V}+\frac{U^{\prime }}{4V}\left(
\frac{V^{\prime }}{V}+\frac{U^{\prime }}{U}\right) -\frac{U^{\prime
}}{rV} \nonumber \\ =-4\pi G(\epsilon +3p)U+\Lambda U,  \label{18}
\end{eqnarray}
Here single prime denotes first derivative and double prime denotes
second derivative.

We have the equation for hydrostatic equilibrium,
\begin{equation}
\frac{U^{\prime }}{U}=\frac{-2p^{\prime }}{p+\epsilon }.  \label{19}
\end{equation}

Let us derive the equation for V(r) alone. Using Eq. (\ref{16}), Eq.
(\ref{17}) and Eq. (\ref{18}), we will get
\begin{equation}
\frac{-V^{\prime }}{rV^{2}}-\frac{1}{r^{2}}+\frac{1}{Vr^{2}}=-8\pi
G\epsilon -\Lambda,   \label{20}
\end{equation}
which can be written as
\begin{equation}
\left( \frac{r}{V}\right) ^{\prime }=1-8\pi G\epsilon r^{2}-\Lambda
r^{2}. \label{21}
\end{equation}

Integrating,
\begin{equation}
V(r)=\left[ 1-\frac{2GM(r)}{r}-\frac{\Lambda }{3}r^{2}\right] ^{-1},
\label{22}
\end{equation}
where
\begin{equation}
M(r)=\int\limits_{0}^{r}4\pi r^{\shortmid 2}\epsilon (r^{\shortmid
})dr^{\shortmid }.  \label{23}
\end{equation}
Eq. (\ref{22}) tells that outside the star, the space-time looks
like SdS.

 Now, we can use Eq. (\ref{19}), Eq. (\ref{22}) and Eq.
(\ref{23}) to eliminate U(r) and V(r) from Eq. (\ref{17}), which
becomes

\begin{equation}
-r^{2}p^{\prime }=GM\left( r\right) \epsilon \left( r\right)
\frac{\left[ 1+\frac{p\left( r\right) }{\epsilon \left( r\right)
}\right] \left[ 1+\frac{4\pi r^{3}p\left( r\right) }{M\left(
r\right) }+\frac{\Lambda r^{3}}{GM\left( r\right) }\right] }{\left[
1-\frac{2GM\left( r\right) }{r}-\frac{\Lambda r^{2}}{3}\right]}.
\label{26}
\end{equation}

Thus the modified TOV equation is,
\begin{equation}
\frac{dp}{dr}=-\frac{GM\left( r\right) \epsilon \left( r\right)
}{r^{2}}\frac{\left[ 1+\frac{p\left( r\right) }{\epsilon \left(
r\right) }\right]\left[ 1+\frac{4\pi r^{3}p\left( r\right) }{M\left(
r\right) }+\frac{\Lambda r^{3}}{GM\left( r\right) }\right] }{\left[
1-\frac{2GM\left( r\right) }{r}-\frac{\Lambda r^{2}}{3}\right] }.
\label{27}
\end{equation}

This equation is different from TOV equation and this becomes TOV
equation for $\Lambda =0$.

Now we can use this new equation to find the maximum mass of neutron
stars. In this work, neutron stars are considered to be composed of
a Fermi gas of degenerate neutrons. The mass-radius and
density-radius relations are plotted and maximum mass limit is
checked.

We have the standard structure equations for a spherically symmetric
star which describe how the mass and the pressure of a star change
with radius:

\begin{equation}
\frac{dm}{dr}=\frac{4\pi r^{2}\epsilon (r)}{c^{2}},  \label{28}
\end{equation}

\begin{equation}
\frac{dp}{dr}=-\frac{G\epsilon (r)m(r)}{c^{2}r^{2}}, \label{29}
\end{equation}

where $G$\ is Newton's Gravitational constant, $\epsilon $\ is the
mass density, $m$\ is the mass up to radius $r$\ and $c$\ is the
velocity of light.

The modified TOV equation is similar to the differential equation
for the pressure, Eq. (\ref{29}), but here has three correction
factors. All three correction factors are larger than one, i.e. they
strengthen the term from Newtonian gravity. Here one has to solve
(Chandrasekhar, 1931)
the coupled differential Equations (\ref{28}) and (\ref{27}). We
must notice that while $\frac{dm}{dr}$\ is positive,
$\frac{dp}{dr}$\ must always be negative. Starting with certain
positive values for $m$\ and $p$\ in a small central region of the
star the mass will increase while the pressure decreases eventually
reaching zero. We can set $m(r=0)=0$\ and , in addition we specify
the central pressure $p(r=0)=p_{0}$\ in order to solve for Equations
(\ref{28}) and (\ref{27}). The behaviour of $m$\ and $p$\ as a
function of radius will become important for numerical calculations.

Both the coupled differential equations depend on energy density. So
we need a relation between pressure $p$ and energy density
$\epsilon$ which is known as the Equation of State(EoS) (Shapiro and
Teukolsky, 1983)
for the neutron stars. For mass calculation, we will use the well
known EoS:
\begin{equation}
p=K\epsilon^\gamma,  \label{32}
\end{equation}
where $K$ is a constant and $\gamma$ is the polytropic exponent.
$\gamma$=$\frac{4}{3}$ in general relativistic case (Silbar and
Reddy, 2004).

We will give the value of cosmological constant (Carmeli and
Kuzmenko, 2001; Carroll, 2000)
as $\Lambda$$\approx$$10^{-46}$$km^{-2}$.
\section{Result and Discussion}
\label{res}
 The mass-radius and density-radius relations obtained is
shown in Fig. \ref{1} and Fig. \ref{2}. For calculational purpose,
we put $G$=1 and $c$=1 and then pick the units of all dimensionful
quantities to be powers of kilometers.

Converting the results back from natural units, we will get the
maximum mass of neutron star $M$=0.74$M_{\odot}$ and radius
$R$=6.98$km$.

Comparing with the result that obtained when calculated using the
normal TOV equation, we found that the obtained values are lesser.
So we can conclude that there is a slight influence of cosmological
constant in the stellar structure. The influence of cosmological
constant will be higher if the value is higher than the recently
accepted one. Future advancements in dark energy research may
strengthen the importance of cosmological constant. The modified TOV
equation could also be used in the calculations of surface tension
of compact stars (Bagchi, et al.,
2005; Dey, et al., 1998)

\begin{figure}
  \includegraphics{fig1.eps}
\caption{Mass-Radius plot for neutron stars for a central density
$1\times 10^{16}gm/cm^{3}$.}
\label{fig:1}
\end{figure}

\begin{figure}
  \includegraphics{fig2.eps}
\caption{Density-Radius plot for neutron stars for a central density
$1\times 10^{16}gm/cm^{3}$.}
\label{fig:2}
\end{figure}