the physical mechanisms responsible for the generation of primordial perturbations seeding present - day large - scale structure , may leave their imprint in the form of small deviations from a gaussian distribution of the primordial perturbations .
searches for primordial non - gaussianity can thereby provide key information on the origin and evolution of cosmological structures ( e.g. , ref . @xcite and references therein ) .
although the standard single - field , slow - roll , canonical kinetic energy and adiabatic vacuum state inflation generates very small non - gaussianity , any inflationary model that deviates from this may entail a larger level of it ( refs .
@xcite and references therein ) .
deviations from gaussian initial conditions are often taken to be of the so - called local type and are parameterized by the constant dimensionless parameter @xmath7 @xcite : @xmath8 where @xmath9 denotes bardeen s gauge - invariant potential ( evaluated deep in the matter era ) and @xmath10 is a gaussian random field . is linearly extrapolated to @xmath11 , while in the cosmic microwave background ( cmb ) convention it is evaluated deep in the matter era .
thus , @xmath12 f^{\rm cmb}_{\rm nl}\sim1.3f^{\rm cmb}_{\rm nl}$ ] , where @xmath13 denotes the @xmath14-induced linear growth suppression factor
. in this paper we will use the cmb convention . ]
recently , a method for constraining non - gaussianity from large - scale structure surveys has been proposed @xcite , which exploits the fact that the clustering of dark matter halos where galaxies form on large scales is modified in a scale - dependent way by the presence of even small amount of non - gaussianity .
in particular , a non - gaussianity described by eq .
( [ eq : fnl ] ) , introduces a scale - dependent boost ( for @xmath15 and a suppression for @xmath16 ) of the halo power spectrum proportional to @xmath17 on large scales ( @xmath18mpc ) , which evolves roughly as @xmath19 .
large - scale structure ( lss ) surveys covering large volumes are needed to access the scales where the signal arises ( e.g. , ref . @xcite and references therein ) .
the possibility of using high redshift data to constrain non - gaussianity has also been addressed by means of hydrodynamic simulations of the intergalactic medium @xcite . among the many currently - available tracers of the lss , the radio sources from the nrao vla sky survey ( nvss ) @xcite and the quasar catalogue of the sloan digital sky survey release six ( sdss dr6 qsos ) @xcite are particularly interesting since they span large volumes extending out to substantial redshifts @xcite .
indeed these source samples were shown to provide tight constraints on primordial non - gaussianity by the cross - correlation measurements between cmb temperature fluctuations and the lss number density which could be used to detect the late integrated sachs - wolfe ( isw ) @xcite effect @xcite .    in recent papers
@xcite it was shown that the observed nvss and sdss dr6 qsos auto - correlation function ( acf ) and cross - correlation function ( ccf ) hint at a positive value of @xmath20 at more than @xmath21 .
this is because the acf is found to be still positive on angular scales @xmath22 , which , for the median source redshift ( @xmath23 ) , correspond to linear scales where the correlation function should be negative if the density fluctuation field is gaussian . a positive @xmath7 adds power on large angular scales , accounting for the observed acf .
in this paper , we will revisit the constraints on the primordial non - gaussianity from these two lss tracers using the power spectrum @xmath0 , instead of the correlation function @xmath24 .
while it is true that the angular correlation function and the angular power spectrum are spherical harmonic pairs and thus in principle they should enclose the same physical information , each of the two measurements is affected by systematic errors in different ways .
for example the correlation function is less affected by the survey mask than the power spectrum but error estimation for the power spectrum is much easier then for the correlation function .
a direct comparison of parameter constraints obtained via the correlation function and the power spectrum analyses offer therefore a way to quantify possible systematic effects .
we will also use the recently updated sdss ii lrg photometric redshift survey , megaz - lrg ( dr7 ) @xcite to constrain the primordial non - gaussianity . ref .
@xcite used the lrg angular power spectrum for a sample close to the sdss dr3 to constrain @xmath7 . here
we revisit the analysis using an improved bias model , improved lrg catalog , updated knowledge of the sources redshift distribution , up - to - date cmb maps and complementary analysis methods .
furthermore , we explore the constraints on physically motivated primordial non - gaussianity shapes that are different from the local case .
we consider non - gaussianities which effects on the halo clustering is well described by the equilateral and enfolded templates .
we consider the angular power spectra of these three lss tracers , as well as their cross - correlation power spectra with the cmb temperature fluctuation . finally , we consider the cubic correction on the halo bias , which can be motivated in scenarios like the curvaton model , in which a large cubic correction can be produced while simultaneously keeping the @xmath7 correction small @xcite .
the structure of the paper is as follows . in section [ theory ]
we review the effects of primordial non - gaussianity on the power spectrum @xmath0 .
section [ data ] contains the analysis of power spectrum @xmath0 for nvss radio sources , sdss dr6 qsos and megaz - lrg . in section [ method ]
we present the method used to derive constraints on @xmath7 .
section  [ result ] contains our main results .
we conclude with a discussion in section [ summary ] .
the effects of non - gaussianity on the source clustering properties arise because a non - zero @xmath7 affects the halo mass function and enhances the halo clustering on large scales .
the second effect is the dominant one on large scales .    in the presence of primordial non - gaussianity
, the mass function @xmath25 can be written in terms of the gaussian one @xmath26 , for which a good fit to the results of simulations is provided by e.g. , the sheth - tormen formula @xcite , multiplied by a non - gaussian correction factor @xcite : @xmath27~,\end{aligned}\ ] ] where the normalized skewness of the density field , @xmath28 , is @xmath29 , @xmath30 denotes the rms of the dark matter density field linearly extrapolated to @xmath11 and smoothed on scale @xmath31 corresponding to the lagrangian radius of a halo of mass @xmath32 , and @xmath33 is the critical density for ellipsoidal collapse , calibrated on n - body simulations @xcite . for high peaks ( @xmath34 ) and small @xmath7 ,
@xmath33 is slightly smaller than the critical density for spherical collapse , @xmath35 where @xmath36 is the linear growth factor , and @xmath37 evolves very weakly with redshift .
the large - scale halo bias is also modified by the presence of non - gaussianity @xcite : @xmath38 where the factor @xmath39 encloses the scale and halo mass dependence . here
, we consider three types of non - gaussianity .    in practice , we find that , on large scales , @xmath40 , is independent of the halo mass for the local type .
the factor @xmath39 for the enfolded template is proportional to @xmath41 , while the equilateral template gives an almost scale - independent @xmath39 .
it is important to note here that the equilateral and enfolded templates were created to have factorizable expressions that gave a good description of the actual bispectra , on average , over all configurations .
for the the halo bias effect , nearly squeezed configurations dominate . as discussed e.g. , in ref .
@xcite the templates may not necessarily offer a good approximation of the actual bispectra generated by higher - order derivatives or vacuum - state modifications .
for all inflationary single - field models the bispectrum scales either as @xmath41 or as @xmath42 in the squeezed limit ( i.e. either like the equilateral template or like the local case ) @xcite .
the orthogonal template instead scales like @xmath17 in the squeezed limit .
nevertheless it is interesting to consider such an intermediate case : inflationary models that go beyond single field can show a range of scaling in the squeezed limit ( e.g. , ref .
@xcite ) .
we assume that the large - scale linear halo bias for the gaussian case is @xcite : @xmath43+\frac{2p}{\delta_{\rm c}(z_{\rm f})d(z_{\rm o})}\left\{1+\left[\frac{q\delta^2_{\rm c}(z_{\rm f})}{\sigma^2_{\rm m}}\right]^p\right\}^{-1}~,\end{aligned}\ ] ] where @xmath44 is the halo formation redshift , and @xmath45 is the halo observation redshift . as we are interested in massive haloes , we expect that @xmath46 . here ,
@xmath47 and @xmath48 account for the non - spherical collapse and are a fit to numerical simulations ( see also refs .
@xcite ) . finally , the weighted effective halo bias is given by @xmath49 @xmath50 being the minimum halo mass hosting a source of the kind we are considering .
, of the different types of objects considered for our analysis .
nvss radio sources ( red dashed line ) , sdss dr6 qsos ( black solid line ) , and sdss dr6 lrgs ( blue dash - dotted line ) . ]
in this subsection we describe three lss probes , nvss radio sources , high redshift sdss dr6 qsos and megaz dr7 lrgs sample . in figure [ fig : redshift ] we show their redshift distribution , @xmath51 .
all distributions are normalized to unit integral .
the nvss @xcite offers the most extensive sky coverage ( @xmath52 of the sky to a completeness limit of about 3 mjy at 1.4 ghz ) .
we start by confining our analysis to nvss sources brighter than 10 mjy , since the surface density distribution of fainter sources suffers from declination - dependent fluctuations @xcite .
density gradients in the nvss catalog become increasingly unimportant as the source brightness threshold is increased .
figure 6 of ref .
@xcite , figure 1 of ref .
@xcite and figure 3 of ref .
@xcite showed variations of nvss source surface density as a function of declination for sources with different flux thresholds .
dim sources are strongly affected by density gradients , however this does not appreciably happen for the brighter sources .
furthermore , in our previous work we used the nvss sources with different flux thresholds @xmath53 mjy and @xmath54 mjy to compute the acf of nvss source and its ccf with cmb temperature fluctuations and found that the obtained acf and ccf results are stable @xcite .
also we mask the stripe @xmath55 , where the catalog may be substantially affected by galactic emissions . here
, we have explicitly checked that the choices of masked stripes with @xmath56 and @xmath57 have negligible impact on our final results . in order not to excessively reduce the sample , in our analysis we consider the nvss sources brighter than 10 mjy and the masked strip with @xmath55 to reduce the contamination from galactic emissions .
the nvss source surface density at this threshold is @xmath58 and the redshift distribution at this flux limit has been recently determined by ref .
their sample , complete to a flux density of 7.2 mjy , comprises 110 sources with @xmath59mjy , of which 78 ( 71% ) have spectroscopic redshifts , 23 have redshift estimates via the @xmath60@xmath61 relation for radio sources , and 9 were not detected in the @xmath60 band and therefore have a lower limit to @xmath61 .
we adopt here the smooth description of this redshift distribution given by ref .
@xcite , which is shown in the red dashed line of figure [ fig : redshift ] : @xmath62 the mean redshift of this sample is @xmath63 and the redshift range is @xmath64 .      the sdss dr6 qso catalog released by ref .
@xcite contains about @xmath65 objects with photometric redshifts ranging from @xmath66 to @xmath67 over a total area of @xmath68 ( @xmath69 of the whole sky ) .
we refer the reader to ref .
@xcite for a detailed description of the object selection with the non - parametric bayesian classification kernel density estimator ( nbc - kde ) algorithm .
we use the electronically - published table that contains only objects with the `` good '' flag with values in the range @xmath70 $ ] . the higher the value , the more probable for the object to be a real qso ( see  4.2 of ref .
@xcite for details ) . furthermore we restrict ourselves to the `` uvxts=1 '' , i.e. to qsos clearly showing a uv excess which should be a signature of a qso spectrum ( in this case we have @xmath71 qsos ) . as we know
, several systematics effects , including the galactic extinction by dust , sky brightness , number of point sources and poor seeing , could potentially affect both the observed auto-/cross - correlation power spectra . in our previous work we carefully checked for their contribution in calculations and in particular we considered extinction and point sources contamination .
we found that the dominant systematic effect is the extinction @xcite .
therefore , in order to minimize the effect of galactic extinction on the observed qso distribution , we use the extinction mask @xmath72 only in our following analysis @xcite .
these masks will remove about @xmath7320% of the considered area .
we fit the redshift distribution @xmath74 of the dr6 qso sample with a function of the form @xcite ( black solid line of figure [ fig : redshift ] ) : @xmath75~.\label{reddis}\ ] ] the best - fit values of the parameters are @xmath76 , @xmath77 , @xmath78 ; the mean redshift of the sample is @xmath79 .
we use the updated megaz lrg dr7 sample @xcite , which contains @xmath80 galaxies from the sdss dr7 in the redshift range @xmath81 with limiting magnitude @xmath82 . to reduce stellar contamination
, there is a variable of the megaz neural network estimator @xmath83 , defined such that @xmath84 if the object is a galaxy , and @xmath85 if it is a star . for a conservative analysis
, we choose a cut @xmath86 , which is reported to reduce stellar contamination below 2% while keeping 99.9% of the galaxies .
in addition to the sdss dr7 geometry mask , we also add two foreground masks to account for seeing ( removing pixels with median seeing in the red band larger than 1.4 arcsec ) and reddening ( removing pixels with median extinction in the red band @xmath87 ) @xcite .
the redshift distribution function ( blue dash - dotted line of figure [ fig : redshift ] ) in this case is found directly from the photometric redshifts that are given in the catalogue @xcite , which is similar to a gaussian function : @xmath88}~,\ ] ] where the mean redshift @xmath89 and the deviation @xmath90 .
a two dimensional density field @xmath91 defined over the full sky can be decomposed in a series of spherical harmonics @xmath92 and their corresponding coefficients @xmath93 : @xmath94 with @xmath95 on the full sky spherical harmonics are a complete and orthogonal basis set . however , in practice , a masked ( incomplete ) sky is observed .
thus the observed @xmath96 coefficients on the partially cut sky become @xcite : @xmath97 where @xmath98 is the position dependent , weight function imposed by the mask and , optionally , the adopted weighting scheme . from these coefficients
we can obtain the observed pseudo power spectrum @xmath99 , the pseudo-@xmath0 estimator @xcite : @xmath100 the pseudo power spectrum @xmath99 , given by the weighted spherical harmonic transform of a masked map , is clearly different from the full sky true power spectrum @xmath101 .
the expectation value of @xmath99 is related to the true @xmath101 by a convolution : @xmath102 where the coupling matrix @xmath103 , describing the mode coupling resulting from the weight function @xmath98 @xcite , can be expressed in terms of @xmath104 symbols as : @xmath105 @xmath106 is the angular power spectrum of the weight function : @xmath107 where the weight coefficients @xmath108 are : @xmath109 for the lss catalogues presented above , in practice , we construct the pixelized maps using the healpix software package @xcite with high resolution @xmath110 , yielding pixel areas of @xmath111 and estimate their pseudo power spectrum @xmath112 using the healpix function `` anafast '' .
we also subtract the shot noise contribution @xmath113 from the observed pseudo power spectrum , where @xmath114 is the surveyed area and @xmath115 is the observed number of sources . the pseudo power spectrum is measured up to @xmath116 , at these multipoles the magnitude of the estimated error bars is dominated by shot - noise . in order to avoid the effects of gauge corrections on the power spectrum on very large scales @xcite , we set the minimal multipole @xmath117 .
the effects of non - gaussianity are large on very large scales . by ignoring @xmath118 we are neglecting the scales where the signal is maximal .
however the theoretical predictions for the measured correlations is uncertain on very large scales as discussed in ref .
@xcite . in ref .
@xcite we estimated that these uncertainties could introduce a systematic error on @xmath7 of @xmath119 . here
we take a more conservative approach and simply ignore the largest scales , even if this implies a reduction of the signal - to - noise .      in the signal dominated limit
, the inverse of the power spectrum covariance matrix @xmath120 is @xcite : @xmath121 where @xmath0 is the theoretical angular power spectrum , @xmath122 denotes the shot noise term . @xmath0 and @xmath123 are the necessary error contributions from both cosmic variance and shot noise , respectively . for the full sky survey , @xmath124 and
the covariance matrix becomes diagonal : @xmath125      besides the angular power spectrum of each lss tracer , we also consider the cross - correlation power spectrum among them . in figure
[ fig : redshift ] , we show that the redshift distributions of these three tracers overlap .
we can thus expect some cross - correlation signal among them .
we also use the healpix software to estimate the cross - correlation power spectrum @xmath126 between tracer @xmath127 and @xmath128 : @xmath129 where @xmath130 are the observed coefficients of tracer @xmath131 and they are computed on the common sky area only . consequently , the inverse of the covariance matrix becomes : @xmath132 where the shot noise term @xmath133 is scale independent , @xmath114 is the common area of two surveys , @xmath134 is the observed number of sources in this common area .
similarly , we also estimate the cross - correlation power spectrum @xmath135 between the lss number density map and the cmb map .
while the planck satellite is operating @xcite , here , we use the the 7-years internal linear combination ( ilc ) map of the cmb provided by the wilkinson microwave anisotropy probe ( wmap ) team @xcite .
since an estimate of the uncertainty associated to this map is not provided , we have checked that using other cmb maps ( e.g. the 5 years ilc map by the wmap team ) does not change the results in any significant way . for the cross - check , we also use the v and w band all - sky 7-years cmb maps provided by wmap @xcite and smooth these maps to the resolution of ilc map ( 60 arcmins ) . using these maps
, we perform all the calculations again and find that our results are undistinguishable .
we adopt the wmap kq75 mask , excluding about 30% of the sky at low galactic latitude , to avoid most of the residual galactic contamination .
we still use the estimator of cross - correlation power spectrum @xmath136 : @xmath137 where @xmath138 are the observed coefficients of the cmb map @xmath139 , also computed on the common sky area only .
the inverse of the covariance matrix is : @xmath140 where @xmath141 is the theoretical angular power spectrum of cmb map , @xmath142 is measurement error of the cmb temperature power spectrum , in which the main source is the pixel noise which accounts for the imperfections in measured temperature on the sky .
pixel noise for wmap is reported as @xmath143 , where @xmath144 is the noise per observation and @xmath145 is the number of observations per given pixel @xcite . after getting the theoretical @xmath146 and @xmath147
, we include the correction of gaussian beam @xmath148 ( using the healpix gaussbeam routine ) : @xmath149 , @xmath150 .
we neglect the pixel window function since its contribution at @xmath151 is negligible at the resolution of the wmap ilc map which is @xmath152 arcmins fwhm .
we use the publicly available package camb@xmath153sources @xcite to calculate the theoretical power spectrum . in our analysis
, we use the halofit @xcite built - in routine for non - linear correction to obtain the fully - evolved , nonlinear matter power spectrum @xmath154 at any epoch .
we perform a global fitting of cosmological parameters , including @xmath7 , for the data of section [ data ] including also the datasets described below , using the cosmomc package @xcite , a markov chain monte carlo ( mcmc ) code , modified to calculate the theoretical acf and ccf .
we assume purely adiabatic initial conditions and a flat universe , with no tensor contribution to primordial fluctuations .
the following six cosmological parameters are allowed to vary with top - hat priors : the dark matter energy density parameter @xmath155 $ ] , the baryon energy density parameter @xmath156 $ ] , the primordial spectral index @xmath157 $ ] , the primordial amplitude @xmath158 \in [ 2.7,4.0]$ ] , the ratio ( multiplied by 100 ) of the sound horizon at decoupling to the angular diameter distance to the last scattering surface @xmath159 $ ] , and the optical depth to reionization @xmath160 $ ] . the pivot scale is set at @xmath161mpc@xmath162 and do not consider massive neutrinos and dynamical dark energy . besides these six basic cosmological parameters , we have two more parameters related to the power spectrum @xmath0 data : the non - gaussianity parameter , the minimal halo mass @xmath50 .
when we combine all the different probes together we have four parameters : three minimal halo masses for three lss surveys and the non - gaussianity parameter .
the model power spectrum @xmath163 is compared with the observed values @xmath99 , respectively , through the gaussian likelihood function @xcite : @xmath164 where the fisher matrix @xmath120 is the inverse of the power spectrum covariance matrix . while it is true that at very low @xmath165 a gaussian likelihood is not a good approximation for cmb data , at high @xmath165 the gaussian approximation works well even for cmb , for lss it is customary to adopt a gaussian likelihood .
the following cosmological data are also included in the fit : @xmath166 power spectra of cmb temperature and polarization anisotropies ; @xmath167 baryonic acoustic oscillations ( baos ) in the galaxy power spectra ; @xmath168 snia distance moduli .    to deal with the 7-years wmap ( wmap7 )
cmb temperature and polarization power spectra we use the routines for computing the likelihood supplied by the wmap team @xcite .
the wmap7 data are used only to improve the constraints on the six basic cosmological parameters , not to constrain @xmath7 .
the baos @xcite can , in principle , measure not only the angular diameter distance , @xmath169 , but also the expansion rate of the universe , @xmath170 .
however , the limited accuracy of current data only allows us to determine the ratio between the distance scale defined by ref .
@xcite : @xmath171^{1/3},\ ] ] and the comoving sound horizon at the baryon - drag epoch , @xmath172 ( see ref .
accurate determinations of the distance ratio @xmath173 , @xmath172 have been obtained by ref .
@xcite : @xmath174 we adopt these values as gaussian priors .
the snia distance moduli provide the luminosity distance as a function of redshift @xmath175 which provides strong constraints on the dark energy evolution . in this paper
we will use the latest snia data sets from the supernova cosmology project , `` union2 compilation '' which consists of 557 samples and spans the redshift range @xmath176 @xcite . in this data
set , they improved the data analysis method by using and refining the approach of their previous work @xcite . when comparing with the previous `` union compilation '' , they extended the sample with the supernovae from refs .
the authors also provide the covariance matrix of data with and without systematic errors and , in order to be conservative , we include systematic errors in our calculations . in the calculation of the likelihood from snia we marginalize over the nuisance parameter as done in refs .
@xcite .    furthermore , we add a prior on the hubble constant , @xmath177 km / s / mpc given by ref .
@xcite . finally , in the analyses we set the minimal halo mass at @xmath178 consistent with the results of ref . @xcite .
in this section , we present constraints on the primordial non - gaussianity parameter @xmath7 from the observed lss power spectra , in combination with the external data sets introduced above : the wmap7 , bao and snia .
these extra datasets are only used to constrain the underlying cosmology . in table
[ tab : i ] and figure [ fig : fnl ] we show the constraints on the primordial non - gaussianity parameter from different data combinations , after marginalizing over all the other parameters .
we also plot the observed power spectra data @xmath179 in figure [ fig : data1 ] and figure [ fig : data2 ] , together with the corresponding theoretical best fit models for comparison . for illustration purposes , here we bin the observed power spectra with the bin size @xmath180 .
we begin by considering the constraint on the local type of non - gaussianity @xmath7 from the angular power spectrum only . in the first panel of figure [ fig : data1 ]
, we plot the observed nvss angular power spectrum data , which is consistent with the previous work @xcite .
when using the nvss @xmath99 data only , we obtain the marginalized constraint on the non - gaussianity parameter : @xmath181 ( @xmath182 ) and @xmath183 ( @xmath21 ) , which is consistent with zero .
the minimal halo mass , @xmath184 ( @xmath182 ) , turns out to be remarkably close to our previous work @xcite .
we show the best fit model in this panel ( black dashed line ) . the amplitude of the power spectrum on large scales ( small @xmath165 ) is enhanced to fit the data points , by the positive value of @xmath7 . in our previous works @xcite
, we used the auto - correlation function ( acf ) of nvss radio sources to constrain @xmath7 and found that the nvss acf yields a positive @xmath7 at more than @xmath21 .
the reason for this difference is that we only use the power spectrum data at @xmath185 , and the power spectrum approach does not include the information on the total number of nvss radio sources @xmath186 ( which was folded in the acf analysis ) .
if we add a prior on the total number of sources @xmath187 : @xmath188 , we obtain a limit of @xmath189 at @xmath182 .
this result is perfectly consistent with previous works @xcite .
next , we use the estimated angular power spectrum of the sdss dr6 qso sample , which is shown in the second panel of figure [ fig : data1 ]
. there is an obvious excess power at large scales , which favors a positive @xmath7 : @xmath190 ( @xmath182 ) and @xmath191 ( @xmath21 ) .
similarly to our previous work @xcite , we also find evidence for positive @xmath7 at more than @xmath21 . however , ref .
@xcite , using sdss qso data , obtained non - gaussianity constraints consistent with zero at 95% confidence level . as we discussed before , ref .
@xcite used an extension of the dr3 qso sample that include sources that subsequently were released with the dr6 sample : we use the more complete and better calibrated final official sdss dr6 qso catalog release @xcite .
@xcite found that their qso sample at @xmath192 ( qso0 ) seem to suffer from contamination which they attribute to systematic calibration errors , discarded this sample from their analysis and only used the high redshift sample ( qso1 ) .
we have checked for this effect and find no evidence in the sample we use that the @xmath192 qsos have different contamination than the @xmath193 @xcite .
in addition , we remove three narrow stripes in the southern hemisphere of the galaxy and set the minimal multipole @xmath117 , which could reduce the effect of star contamination further .
finally , we consider the megaz dr7 lrg sample .
@xcite split this sample in four redshift bins according to each source photometric redshift estimate and found a large excess of power over the lowest multipoles ( @xmath194 ) in the angular power spectrum , which is also found in their previous analysis using megaz dr4 lrg sample @xcite . in our analysis , although we do not split the lrg sample in redshift bins , we still find this excess power on large scales , which is shown in the third panel of figure [ fig : data1 ] .
there are several explanations for this excess power that range from systematic errors as due to contaminants to new physics @xcite .
@xcite checked some possible systematic errors of this megaz lrg sample in details and found that systematics errors do not seem to be responsible for the excess power at large scales .
we also perform a similar check and find no significant changes in the power spectrum at large scales .
one of the possibilities is that this excess power could be induced by the primordial non - gaussianity . here
, we use the observed @xmath99 data of megaz dr7 lrg sample to constrain the primordial non - gaussianity for the first time .
we calculate the angular power spectrum on largest scales and account for the redshift distortion power as described in ref .
@xcite , although we set the minimal multipole @xmath117 .
after marginalizing over other free parameters , we obtain the constraint : @xmath195 which is consistent with zero at @xmath21 .
we also obtain a limit of @xmath196 at 68% confidence level . using eq .
( [ eq : nghalobias ] ) , we could obtain the limit of effective halo bias at the mean redshift of the survey : @xmath197 ( @xmath182 ) , which is consistent with other works @xcite , although in our analysis we directly use the full sample .
when we combine these three observed angular power spectra data together to constrain non - gaussianity , we obtain a limit of @xmath198 at 68% confidence level and @xmath199 at 95% confidence level , which implies the current angular power spectra of lss tracers favor @xmath15 at the @xmath200 level , since non - gaussianity adds power on large angular scales yielding a good fit to the observed data points .
as can be seen from figure [ fig : data1 ] and table [ tab : i ] , there is not a specific measurement that is driving the signal : all the different probes yield very consistent results and when combined the error bar on @xmath7 is reduced to yield a @xmath201 result .
comparing with the angular power spectrum , the constraining power of cross - correlation power spectrum on @xmath7 is much weaker @xcite . in the panels of figure [ fig : data2 ] , we plot the observed cross - correlation power spectra @xmath126 among three lss tracers and @xmath136 between the cmb map and three lss tracers .
these cross - correlation data points are consistent with zero within @xmath182 error bar but the best fit value @xmath202 is an excellent fit to the data points , as shown in figure [ fig : data2 ] . however , cross - correlation data give a larger error bar for @xmath7 than that from the angular power spectra data and are thus much less constraining .
first , we consider the single tracer case in combination with its cross - correlation with the cmb map .
when we combine the angular power spectrum and its cross - correlation power spectrum @xmath136 with cmb map to constrain the primordial non - gaussianity , we obtain @xmath203 , @xmath204 and @xmath205 at 68% confidence level for nvss radio sources , sdss dr6 qsos and megaz dr7 lrgs , respectively .
the results are consistent with those obtained from the angular power spectrum only , but those cross correlations does not seem to add much information .    if we consider all the angular and cross - correlation power spectra data sets of three lss tracers together , . ]
we obtain the constraint on the local type primordial non - gaussianity : @xmath206 this result is compatible with previous estimates @xcite , and with the wmap7 limits @xcite .
the lss data still implies a positive @xmath7 at more than @xmath21 . in figure
[ fig : data1 ] and figure [ fig : data2 ] , we plot the theoretical best fit model @xmath202 of our non - gaussianity calculations ( red solid lines ) .
these curves match the observed power spectra very well , especially on large scales .
generally , the correction to the standard halo bias due to the presence of primordial non - gaussianity is @xmath207 , where the function @xmath208 is given by @xcite : @xmath209 where @xmath210 is the halo mass which is related to the scale @xmath31 , @xmath211 denotes the expression for the primordial bispectrum of the bardeen potential @xmath9 , @xmath212 its power spectrum and @xmath213 .
in the local non - gaussian case , the expression for the correction to the standard gaussian halo bias simplify [ eq . ( [ eq : nghalobias ] ) ] . in this subsection
, we explore other two different types of primordial non - gaussianity given by the equilateral and enfolded templates .
the equilateral type of non - gaussianity @xcite , which can be generated from the inflationary models with higher - derivative operators of the inflaton , can be well described by the following template @xcite : @xmath214 , where @xmath215^{2/3}+(p_\phi(k_1)^{1/3}p_\phi(k_2)^{2/3}p_\phi(k_3)+5cyc)~.\ ] ] then we find that the correction @xmath208 in the equilateral non - gaussian case is almost scale - independent ( see figure 1 of ref .
@xcite ) . in principle
there is a mass dependence ( see discussion in ref .
@xcite ) , but the template and the actual physical bispectrum yield different effective @xmath7 ( denoted by @xmath216 ) for the halo bias . here , therefore , we consider a phenomenological scale - independent scaling and neglect the mass dependence .
the interpretation of the resulting constraint should be interpreted as a constraint on an effective @xmath7 , its connection to the physical bispectrum depends on the detailed model of inflation under consideration .
using all datasets , we obtain the constraint on the equilateral type of primordial non - gaussianity : @xmath217 which is consistent with the wmap7 results @xcite .
the template of another type of non - gaussianity `` enfolded '' , which arises in models with non - bunch - davies initial state @xcite or in effective field theories of inflation with higher - derivative interactions @xcite ( see also ref .
@xcite ) , is given by ref .
@xcite : @xmath218 , where @xmath219^{2/3}-(p_\phi(k_1)^{1/3}p_\phi(k_2)^{2/3}p_\phi(k_3)+5cyc)~.\ ] ] in figure 1 of ref .
@xcite , it can be seen that the correction in this type is proportion to @xmath220 on large scales .
no single field inflationary model has this scaling in the squeezed limit : all inflationary single field models considered so far have only the local or the template equilateral scaling in this limit @xcite .
this is nevertheless an interesting intermediate case to consider .
in fact non - single field models can yield all the intermediate scalings in the squeezed limit , see e.g. , ref .
we therefore consider a phenomenological scaling like @xmath220 which could cover this enfolded template case and neglect the mass dependence . in this case
, the whole data sets give the constraint on the effective non - gaussianity parameter for the enfolded template : @xmath221 in the wmap7 paper @xcite , a constraint on the orthogonal type of non - gaussianity @xcite is given .
the correction to the halo bias of this template is also proportional to @xmath220 on large scales , but have the opposite sign and a normalization factor of @xmath222 with a weak mass dependence which we neglect here .
therefore , based on the result of enfolded template , in our simple analysis , we obtain the constraint on the non - gaussianity parameter for the orthogonal template : @xmath223 ( @xmath182 ) and @xmath224 ( @xmath21 ) .
the orthogonal type of non - gaussianity in the squeezed limit however scales like the equilateral one and , as shown in ref .
@xcite its effect on the halo bias is almost degenerate with the equilateral non - gaussianity .
the enfolded type of non - gaussianity arising for example from modified initial state , in the squeezed limit scales like the local template and as shown in ref .
@xcite its effect on the halo bias is almost degenerate with the local non - gaussianity but with a different normalization @xmath225 .
this re - scaling can be used to re - interpret the local constraints into constraints for the modified initial state non - gaussianity .
finally , we consider the effect of the cubic correction on the halo bias model . in this case , the bardeen potential @xmath9 can be conveniently parameterized up to third order by : @xmath227 where @xmath226 is dimensionless , phenomenological parameter .
@xcite explored the effect of a cubic correction @xmath5 on the mass function and halo bias model in detail and gave the expression of the scale - dependent bias correction when @xmath228 and @xmath229 : @xmath230 where the normalized skewness of the density field @xmath231 and the parameter @xmath232 obtained by the fitting in ref . @xcite .
then we obtain the constraint on @xmath226 from all the data sets : @xmath233 the result is compatible with that obtained in refs .
@xcite .
c|c|c    datasets & +    wmap7+bao+sn &  @xmath182  &  @xmath21   +     + @xmath234(nvss)&@xmath235&@xmath236 $ ] + @xmath237(nvss)&@xmath238&@xmath239 $ ] +    @xmath234(qso)&@xmath240&@xmath241 $ ] + @xmath237(qso)&@xmath242&@xmath243 $ ] +    @xmath234(lrg)&@xmath244&@xmath245 $ ] + @xmath237(lrg)&@xmath246&@xmath247 $ ] +    @xmath234(all)&@xmath248&@xmath249 $ ] + @xmath250(all)&@xmath251&@xmath252 $ ] +     + @xmath250(all)&@xmath253&@xmath254 $ ] +     + @xmath250(all)&@xmath255&@xmath256 $ ] +     + @xmath250(all)&@xmath257&@xmath258 $ ] +     and @xmath21 limits on the local type of primordial non - gaussianity @xmath7 obtained from different data combinations .
the vertical line denotes the gaussian case @xmath228.[fig : fnl ] ]
searches for primordial non - gaussianity can provide key information on the physical mechanism that generate the primordial perturbations and can constrain specific inflationary models . a relatively new and promising method to detect primordial non - gaussianity exploits the fact that the clustering of dark matter halos ( which host massive galaxies and qsos ) is modified on large - scales in a scale - dependent way .
this effect is called non - gaussian halo bias . to access these large scales ,
large - scale structure surveys covering large volumes , high redshift and tracers of massive halos are most suited . among the currently available surveys , the nvss and the sdss qso ( in particular the dr6 catalog ) are the largest .
previous work had detected excess large - scale power in their auto - correlation function yielding a detection of a positive value of the non - gaussianity parameter @xmath7 at @xmath259 level .
here we have revisited the issue using the angular power spectrum rather than the correlation function . while the two statistics should enclose the same physical information they are affected by systematics effects in different ways : this complementary analysis offers thus a consistency check and a test for possible systematics .
we have also used the recently updated megaz - lrg survey : a previous release of the same survey had also been found to have excess large - scale power , but mostly at @xmath118 . in an effort to be as conservative as possible in our analysis , we excuded all multipoles at @xmath118 : these multipoles are affected by theoretical uncertainties and systematics such as calibration , sample contamination etc .
these also probe the largest scales , thus these are the multipoles where most of the non - gaussian signal is likely to be concentrated .
we combine the survey auto power spectra with the cross - power spectra among the surveys and of each of the survey with the cmb temperature map . in our analysis
we marginalize over the lcdm cosmological parameters with a prior imposed by wmap7 angular power spectra , bao , supernovae and hubble constant measurement .
we have investigated the constraints on the parameter @xmath7 , characterizing primordial non - gaussianity not only for the popular local shape but also for the equilateral , enfolded and orthogonal templates .
we have also considered constraints on a cubic non - gaussianity correction ( parameterized by the @xmath226 parameter ) .
we find no evidence for non - gaussianity for any `` shape '' which halo bias effect on large scale scales less steeply than @xmath17 and no evidence for non - zero @xmath226 .
however for the local type of non - gaussianity we find a @xmath260 signal for positive @xmath7 when all data - sets are combined .
the results are summarized in figure [ fig : fnl ] : no single data set drives the signal and the central @xmath7 values recovered are fully in agreement with previous analyses . due to the more conservative approach taken here ( and thus larger error - bars ) , the only single survey that yield more than @xmath261 detection is the sdss dr7 qso survey . the nvss survey central recovered @xmath7 value is also in agreement with both the qso measurement and previous analyses but the error - bar is larger .
the combination of the data sets yields : @xmath262 and @xmath263 at @xmath264 and @xmath21 .
thus the tantalizing hint for a positive local @xmath7 found from the auto - correlation function remains despite the different analysis .
we have investigated several possible sources of systematic errors such as the sources number density and the contamination by galactic emissions for nvss sources ; different choices for the cmb temperature fluctuations templates ; possible contamination of stars for the sdss and lrg samples .
our main results are stable when these systematic errors are considered .
such a result would have profound implications for inflationary mechanisms .
the improved statistical power of forthcoming large - scale structure surveys is needed to fully resolve the issue .
moreover , it will be also important to use the cmb data which will be provided by planck .
we acknowledge the use of the legacy archive for microwave background data analysis and the healpix package .
numerical analysis has been performed at the university of cambridge high performance computing service ( http://www.hpc.cam.ac.uk/ ) .
this research has been partially supported by the asi contract no .
i/016/07/0 cofis , the asi / inaf agreement i/072/09/0 for the planck lfi activity of phase e2 and a prin miur .
mv is supported by asi / aae , a prin miur , a prin inaf 2009 and the erc - stg `` cosmoigm '' .
support was given by the italian space agency through the asi contracts euclid - ic ( i/031/10/0 ) .
lv acknowledges support from grant fp7 erc - ideas phys.lss 240117 and micinn grant aya2008 - 03531 .
lv and sm thank the centro de ciencias de benasque pedro pascual where stimulating discussions about this work took place .
funding for the sdss and sdss - ii has been provided by the alfred p. sloan foundation , the participating institutions , the national science foundation , the u.s .
department of energy , the national aeronautics and space administration , the japanese monbukagakusho , the max planck society , and the higher education funding council for england .
the sdss web site is http://www.sdss.org/. the sdss is managed by the astrophysical research consortium for the participating institutions .
the participating institutions are the american museum of natural history , astrophysical institute potsdam , university of basel , university of cambridge , case western reserve university , university of chicago , drexel university , fermilab , the institute for advanced study , the japan participation group , johns hopkins university , the joint institute for nuclear astrophysics , the kavli institute for particle astrophysics and cosmology , the korean scientist group , the chinese academy of sciences ( lamost ) , los alamos national laboratory , the max - planck - institute for astronomy ( mpia ) , the max- planck - institute for astrophysics ( mpa ) , new mexico state university , ohio state university , university of pittsburgh , university of portsmouth , princeton university , the united states naval observatory , and the university of washington .


entanglement , nonlocality and quantum discord  ( qd ) are all different but somehow related aspects of the quantumness of the correlated composite quantum systems .
entanglement is a kind of quantum correlation and determines whether the given state is separable or not  @xcite . on the other hand
, the quantum systems may contain quantum correlations other than entanglement .
quantum discord , which is defined as the difference between the quantum versions of two analogous classically equivalent definitions of mutual information , captures all types of nonclassical correlations including entanglement  @xcite .
although entanglement and qd are equivalent for pure states  @xcite , qd is found to be nonzero for some separable states  @xcite .
such states were shown to be useful in the deterministic quantum computations with one pure qubit context  @xcite .
therefore qd is believed a new resource for quantum computation .
the nonlocality as measured by the violation of bell s inequalities is also a signature of the inseparability of a quantum state and signifies entangled states whose correlations can not be reproduced by a classical local hidden - variable model  @xcite .
however it must be emphasized here that the presence of entanglement does not always imply violation of bell s inequalities , while the violation of bell inequalities implies entanglement .
such entangled states violating bell s inequalities play a central role in some applications in quantum information science  @xcite , such as to guarantee the safety of device - independent key distrubition protocols in quantum cryptography  @xcite .
the relations among different measures of quantum correlations and their dynamical behaviour have been an active area of research with an aim to understand the fundamental question of what quantumness is as well as to characterize these correlations as useful resources for various quantum tasks  @xcite .
one of the most important obstacle to realizing quantum operations is the interaction between the qubits and the environment which tends to wipe out the quantumness of the system  @xcite .
characterization of the effect of environment on the system such as decoherence and dissipation is important .
the dynamics of quantum correlations such as entanglement , quantum discord and bell nonlocality for a two - qubit system subject to various environments have been studied by many groups recently  @xcite .
entanglement sudden death  ( esd ) is one of the new concepts that arise within the context of these studies  @xcite .
esd refers to death of entanglement between two qubits in a finite time while the coherence of the single qubit decays only exponentially . in
contrary to entanglement , qd is more robust in a decohering environment , actually it was shown that almost all states have non - zero quantum discord  @xcite .
qd was shown to decay exponentially in markovian environments , even at finite temperatures  @xcite and can become zero momentarily under non - markovian dynamics  @xcite .
qd was also experimentally investigated in an all - optical set up recently  @xcite .
furthermore , qd and entanglement were analyzed for a heisenberg spin chain with quantum phase transitions and qd was found to detect the critical points of the transition at finite temperatures while entanglement could not  @xcite .
very recently , sun _
_ investigated the induced quantum correlations between two qubits coupled to a common environment modelled by ising spin chain in a transverse field and during the time evolution , non - zero quantum discord is found to be created for specific initial states whereas entanglement can not be induced  @xcite . in ref .
@xcite quantum discord is found to be totally unaffected by non - dissipative noise for a long time interval under certain conditions , while entanglement decays to zero quickly .
these remarkable properties of quantum discord makes it more practical than the entanglement to characterize quantum correlations .
as mentioned , the nonlocality is curicial for some applications in quantum information .
therefore , many theoretical studies have been devoted to the comparision of entanglement and the violations of bell inequalities  @xcite .
it was demonstrated by several authors that the survival time for entanglement is much larger than that of the bell - inequality violation  @xcite . in ref .
@xcite , the nonlocal entanglement identified by the violations of a bell inequality is found to be preserved during the time evolution of a system consisting of two qubits which are embedded independently in zero - temperature bosonic reservoirs and are initially in an entangled mixed state .
_ showed that bell inequality might not be violated for a state with high entanglement for a two - qubit system subject to amplitude damping  @xcite .
moreover , the nonlocal entanglement between two qubits subject to classical dephasing environment modelled as ornstein - uhlenbeck noise has been investigated in refs .
@xcite , very recently , and it was demonstrated that the strong non - markovian effect can protect the nonlocal entanglement  @xcite . additionally , the nonlocal entanglement induced by the qubit - qubit interaction is found to be protected against the sudden death despite the strong influence of the environmental noise for a certain class of initial states regardless of whether the environment is markovian or non - markovian  @xcite .
although the dynamical relation between the entanglement and quantum discord as well as entanglement and bell nonlocality for a bipartite system subject to different environments have been investigated by a large number of groups as mentioned above , further investigations would help in clarifying the relation of these three correlation measures .
along these lines , in the present paper , we have investigated the dynamics of quantum correlations such as quantum discord , entanglement and bell nonlocality between two qubits in a lossy cavity .
very recently , entanglement , bell nonlocality and quantum discord have been studied in cavity qed  @xcite .
most of them assume identical atoms , or the atom - cavity resonant case , or only one excitation . here
, we have investigated the effect of atom - cavity detuning and cavity decay rate on the corresponding correlation measures for a system which is initially prepared in separable or maximally entangled state .
the most important results of the present study are : depending on the initial preparation of the system and the type of the detunings between the frequency of the cavity mode and the transition frequency of the qubits , the induced quantum discord is found to increase for a long time under cavity decay loss mechanism for certain conditions , while there is no increase in entanglement and bell nonlocality . in some situations ,
the quantum discord is induced even while the atoms remain separable at all times .
also , the preservation of induced entanglement and quantum discord was observed under strong dissipation for certain type of initial states .
the paper is organized as follows . in sec .
[ model ] , the model is introduced and the exact differential equations governing the dynamics of two atoms and single cavity mode are derived . the correlation measures ; concurrence , chsh - bell inequality and quantum discord are also introduced , briefly . in secs .
[ cavity1 ] and  [ cavity2 ] , the effect of cavity damping on the corresponding correlation measures are revealed for different detunings for the atoms initially prepared in separable or maximally entangled state .
the summary of the important results is given in sec .
[ conclusion ] .
here we consider two two - level atoms with excited states @xmath0 and ground states @xmath1 interacting with a single - mode cavity field of frequency @xmath2 .
the hamiltonian of this system is given as  @xcite : @xmath3 where @xmath4 and @xmath5 are the transition frequencies for the atoms @xmath6 and @xmath7 , @xmath8 and @xmath9 are the annihilation and creation operators for the field and @xmath10 and @xmath11 are the coupling constants to the cavity mode for the atoms @xmath6 and @xmath7 , respectively .
when maximum two - photon is considered to be present in the cavity  ( i.e. , the total excitation is equal to two ) and if there is no dissipation in the cavity modes , the closed form of the total state vector at any time may be written as  @xcite : @xmath12 where @xmath13  @xmath14 are the time - dependent probability amplitudes .
for example , @xmath15 represents that the atoms are in their ground states while there are two photons inside the cavity . for this case , the dynamics of the atoms and single mode cavity field can be solved by using time - dependent schrdinger equation and recently many studies have been reported for this case  @xcite . on the other hand , when cavity decay is taken into account , since the cavity decay only changes the number of photons inside the cavity , the atoms and the cavity mode can be found in any of the states  @xcite : @xmath16 .
then the density matrix of the system of two atoms interacting with a single mode cavity field will evolve according to the master equation  @xcite : @xmath17-\frac{\kappa}{2}\left(a^{\dagger}a\rho-2a\rho a^{\dagger}+\rho a^{\dagger}a\right),\end{aligned}\ ] ] where @xmath18 and @xmath19  @xmath20 are the spin inversion operators . here the first term in eq .
( [ mastereqn ] ) is generated by the hamiltonian given in eq .
( [ hamiltonian ] ) and represents the coherent evolution , while the other term denotes the incoherent dynamics ; @xmath21 is the rate of loss of photons from the cavity mode . by setting @xmath22 , @xmath23 , @xmath24 , @xmath25 , @xmath26 , @xmath27 , @xmath28 , @xmath29 as a basis and considering the eqs .
( [ hamiltonian ] ) and  ( [ mastereqn ] ) , the differential equations governing the dynamics of two atoms and the single mode cavity field can be obtained easily . in our analysis
we will restrict ourselves for the initial states : @xmath30 and @xmath31 , where @xmath32 .
one should note for these initial states the diagonal elements and some of the off - diagonal elements form a closed evolution .
they can be listed as : @xmath33 where @xmath34 and @xmath35 are called the detunings and describes the frequency difference between the energy levels of atoms and the cavity mode , @xmath36 and @xmath37 .
we will consider two particular detuning types : for @xmath38 the atoms will be called identical since the transition frequencies of the atoms are the same and the atoms will be called unidentical for @xmath39 case for which the individual atom transition frequencies are different  @xcite . in practice , the detunings might be slightly different due to broadening , in this work we will consider relatively large values of @xmath40
. one should note that it is quite difficult to obtain an analytic solution for this set of coupled differential equations , therefore we will solve these equations numerically for the considered initial states .
the density matrix of the two atoms is obtained by taking a partial trace over the cavity degrees of freedom : @xmath41 and used to calculate the quantum correlations as measured by chsh - bell inequality , concurrence and quantum discord between two atoms interacting with a single mode cavity field . for the considered initial states whose dynamics
is given by eq .
( [ solmaster ] ) , the reduced density matrix of the atoms in the two - qubit standart basis @xmath42 can be calculated as : @xmath43{cccc } a & 0 & 0 & 0 \\ 0 & b & e & 0 \\ o & e^ * & c & 0 \\ 0 & 0 & 0 & d \end{array}\right)\end{aligned}\ ] ] where @xmath44 and @xmath45 .
it should be noted that the density matrix in eq .
( [ reducedden ] ) has an x - form and the dynamics preserve its form .
the considered correlation measures for this type of states can be calculated easily  @xcite .
the concurrence  @xcite as an entanglement measure is given by  @xcite : @xmath46 the bell nonlocality measured by the chsh - inequality  @xcite can be found as  @xcite : @xmath47 where @xmath48 and @xmath49 .
qd is defined as  @xcite : @xmath50 where @xmath51 measures the total correlation between the atoms ; @xmath52 is the von - neumann entropy and @xmath53 is the reduced density matrix obtained by tracing @xmath54 over the subsystem @xmath55 .
the other quantity @xmath56 is the amount of classical correlations between the atoms defined as the maximum information one can get about the atom @xmath6 by performing measurements on the atom @xmath7 , or vice versa .
its calculation requires maximization over one of the subsystem and generally obtaining an analytic result is not an easy task  @xcite .
recently , for the density matrix having x - form the explicit expressions of qd and classical correlation are reported  @xcite .
we have used the expressions given in ref .
the classical correlation @xmath56 for the density matrix eq .
( [ reducedden ] ) is given as : @xmath57 where @xmath58 . and
qd is given as : @xmath59 where @xmath60 , with @xmath61 being the eigenvalues of @xmath54 , @xmath62 , @xmath63 , @xmath64 and @xmath65 .
we have also calculated qd numerically by considering projective measurements  @xcite and observed that the results obtained by using eq .
( [ qd2 ] ) and the numerical results agree perfectly .    as is well known , concurrence is equal to 1  ( 0 ) for a maximally entangled  ( separable ) state . although qd is equal to 1 for a maximally entangled state , it may or may not be equal to zero for a separable state because it was shown that even some separable states can carry non - zero quantum discord  @xcite . for pure states
qd is found to be equal to the entanglement of formation  @xcite . on the other hand
when @xmath66 the chsh - bell inequality is violated which signifies that the correlations can not be accessible by any classical local model  @xcite .
in this section , we will examine the effect of cavity decay rate , @xmath21 , on the quantum correlations as measured by concurrence , chsh - bell inequality and quantum discord . to do this
we will solve the differential equations in eq .
( [ solmaster ] ) numerically for the initial states : @xmath30 and @xmath31 , where @xmath32 and we will set @xmath67  ( symmetric coupling ) , @xmath68 for identical atoms and @xmath69 for unidentical atoms .
one should note that the detunings are chosen somewhat large because at higher detunings the qubit - cavity energy exchange is low  @xcite and a high - degree of cavity - induced quantum correlations can be created  @xcite .
quantum discord versus @xmath70 for unidentical atoms  ( @xmath69 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions entanglement and bell nonlocality are not induced , thus these are not plotted here and the inset in  ( c ) represents the total  @xmath77 and classical  @xmath56 correlations versus @xmath70 for @xmath74.,title="fig:",width=283 ]   quantum discord versus @xmath70 for unidentical atoms  ( @xmath69 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions entanglement and bell nonlocality are not induced , thus these are not plotted here and the inset in  ( c ) represents the total  @xmath77 and classical  @xmath56 correlations versus @xmath70 for @xmath74.,title="fig:",width=283 ]     quantum discord versus @xmath70 for unidentical atoms  ( @xmath69 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions entanglement and bell nonlocality are not induced , thus these are not plotted here and the inset in  ( c ) represents the total  @xmath77 and classical  @xmath56 correlations versus @xmath70 for @xmath74.,title="fig:",width=283 ]   quantum discord versus @xmath70 for unidentical atoms  ( @xmath69 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions entanglement and bell nonlocality are not induced , thus these are not plotted here and the inset in  ( c ) represents the total  @xmath77 and classical  @xmath56 correlations versus @xmath70 for @xmath74.,title="fig:",width=283 ]     quantum discord versus @xmath70 for unidentical atoms  ( @xmath69 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions entanglement and bell nonlocality are not induced , thus these are not plotted here and the inset in  ( c ) represents the total  @xmath77 and classical  @xmath56 correlations versus @xmath70 for @xmath74.,width=283 ]    we first consider the dynamics of concurrence and quantum discord for initially separable state @xmath71 and display @xmath78 and @xmath79 versus dimensionless time , @xmath70 in figs .  1 and  2 for unidentical and identical atoms , respectively .
for this initial state neither non - zero concurrence nor violation in bell inequality ever occurs for unidentical qubits , so their dynamics are not plotted in fig .  1 .
however , this does not indicate the loss of quantum correlations ; qd induced by the atom - field interaction oscillates nearly between 0 and 0.18 for @xmath72  ( fig .
a similar result is found by sun _
et al . , _ in ref .
@xcite where the authors studied the entanglement and quantum discord between two qubits which are independently coupled to a common environment modelled as an ising spin chain .
it can be seen from the analysis of figs .
1(b)-1(e ) that qd is created independently of the cavity decay rate , but its maximum value depends on the decay rate @xmath21 inversely and the effect of cavity decay on the dynamics of qd is to create a damped oscillatory behaviour ; damping in qd is proportional to the cavity decay rate .
also , in order to understand the nature of non - zero quantum discord , we have plotted the total correlation , @xmath77 , and classical correlation , @xmath56 , for @xmath74 in the inset of fig .
it should be noted that for this initial state , the total as well as classical correlations are also induced and @xmath77 is always greater than @xmath56 as a consequence non - zero quantum discord is present .
this inset also demonstrates that @xmath77 and @xmath56 are also damped oscillatorily as @xmath80 .
correlations in the system are gradually lost to the reservoir because of the cavity decay as found in the system studied in ref .
@xcite .
quantum discord and concurrence  ( insets in  ( a ) ,  ( b ) and  ( c ) ) versus @xmath70 for identical atoms  ( @xmath68 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions bell nonlocality is not induced , and for @xmath81 and @xmath76 , @xmath82 has no dynamics as well , thus they are not plotted here .
the insets are plotted under the same conditions as @xmath80.,title="fig:",width=283 ] quantum discord and concurrence  ( insets in  ( a ) ,  ( b ) and  ( c ) ) versus @xmath70 for identical atoms  ( @xmath68 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions bell nonlocality is not induced , and for @xmath81 and @xmath76 , @xmath82 has no dynamics as well , thus they are not plotted here .
the insets are plotted under the same conditions as @xmath80.,title="fig:",width=283 ]    quantum discord and concurrence  ( insets in  ( a ) ,  ( b ) and  ( c ) ) versus @xmath70 for identical atoms  ( @xmath68 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions bell nonlocality is not induced , and for @xmath81 and @xmath76 , @xmath82 has no dynamics as well , thus they are not plotted here .
the insets are plotted under the same conditions as @xmath80.,title="fig:",width=283 ] quantum discord and concurrence  ( insets in  ( a ) ,  ( b ) and  ( c ) ) versus @xmath70 for identical atoms  ( @xmath68 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions bell nonlocality is not induced , and for @xmath81 and @xmath76 , @xmath82 has no dynamics as well , thus they are not plotted here .
the insets are plotted under the same conditions as @xmath80.,title="fig:",width=283 ]    quantum discord and concurrence  ( insets in  ( a ) ,  ( b ) and  ( c ) ) versus @xmath70 for identical atoms  ( @xmath68 ) and @xmath71 initial state for the decay rates : ( a )  @xmath72 , ( b )  @xmath73 , ( c )  @xmath74 , ( d )  @xmath75 and ( e )  @xmath76 .
note that under these conditions bell nonlocality is not induced , and for @xmath81 and @xmath76 , @xmath82 has no dynamics as well , thus they are not plotted here .
the insets are plotted under the same conditions as @xmath80.,width=283 ]    we have calculated the dynamics of chsh inequality , concurrence and quantum discord for the same initial state @xmath71 for the identical atoms  ( @xmath68 ) case to understand the effect of the type of detuning for the same leaky cavity and plot the resulting @xmath79 and @xmath78 as a function of @xmath70 in figs .
2(a)-2(e ) and insets of figs .
2(a)-2(c ) , respectively .
as can be seen from the insets of figs .
2(a)-2(c ) , non - zero entanglement , although small , is created which decays to zero as cavity decay rate is increased and the induced entanglement has no dynamics for higher values of @xmath21  ( @xmath81 and @xmath76 ) .
bell inequality is found to be not violated for any value of @xmath21 , which is expected because entanglement is low and generally bell inequality violation is shown to be a hallmark of high @xmath82  @xcite , so @xmath83 is not displayed in fig .  2 .
one should note the difference in the dynamics of qd in the leaky cavity for unidentical and identical atom cases  ( figs .  1 and  2 ) .
while qd has an oscillatory decay type dynamics for the unidentical atoms  ( fig .  1 ) , it increases until the time @xmath84 and then decays for the identical atoms .
also , for the same decay parameters , qd is found to have a much longer lifetime for the identical atoms .
the higher leakage in the cavity decreases @xmath85 and the lifetime of qd after the time @xmath84 . on the other hand
, @xmath21 has no appreciable effect on the maximum value of @xmath80 .
also note that the fast oscillations in qd are gradually smoothed for larger values of @xmath21 and longer times as can be seen from figs .
2(b)-2(e ) .     chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) . here
the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]     chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]     chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) . here
the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) . here
the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]     chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) . here
the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]     chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) .
here the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) . here
the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]   chsh inequality  ( left plots ) , concurrence  ( middle plots ) and quantum discord  ( right plots ) versus @xmath70 for @xmath86 initial state and @xmath72  ( ( a ) ,  ( b ) ,  ( c ) ) , @xmath73  ( ( d ) ,  ( e ) ,  ( f ) ) , @xmath74  ( ( g ) ,  ( h ) ,  ( k ) ) , @xmath81  ( ( l ) ,  ( m ) ,  ( n ) ) and @xmath76  ( ( o ) ,  ( p ) ,  ( r ) ) . here
the main subfigures are for identical atoms with @xmath68 , while the insets are for unidentical atoms with @xmath69 under the same decay rate parameters .
note that for unidentical atoms @xmath83 is not violated , so they are not reported here.,title="fig:",width=192 ]    now we consider another initial separable state of the form @xmath87 which has one excitation in the field and the other excitation in one of the atoms to see the effect of initial partition of the excitation on the dynamics of the quantum correlations .
the time dependent of chsh - inequality , concurrence and quantum discord are displayed in the first , second and third columns of fig .  3 , respectively .
in these figures , insets show the dynamics for unidentical atoms  ( @xmath69 ) while main subfigures are for the identical atoms  ( @xmath68 ) . inspection of subfigures demonstrates the pronounced effect of qubits being identical or unidentical . for identical atoms ,
the atom - field interaction induces high values of @xmath82 and @xmath80 and bell inequality is also violated at the high values of entanglement for a high - quality cavity  ( fig .
3(a ) ,  3(b ) and  3(c ) ) . in the case of unidentical atoms , @xmath82 and @xmath80
are also created , but their maximum are small ; 0.018 and 0.13 , respectively , and @xmath83 is not violated . when unidentical atoms interact with a low - quality cavity , entanglement can cease to exist very quickly and qd shows a fast oscillatory decaying behaviour  ( see the insets in figs .
3(e)-3(r ) ) which is very similar to what we have found for @xmath71 initial state .
on the other hand , this situation is completely and significantly different in the case of identical atoms .
although chsh - bell inequality goes below 2 very quickly for @xmath88 and no longer is violated for @xmath89 , the dynamics of entanglement and quantum discord are very similar that both approach to constant values in an oscillatory manner ; the oscillations die down as the cavity decay rate is increased . in effect ,
cavity decay entangles initially separable state and the resulting entanglement as well as quantum discord are robust . in many studies , the dissipation of cavity fields were found to play a constructive role in the generation of quantum correlations as measured by concurrence and quantum discord  @xcite .
this situation can be understood by looking at the steady - state of the dynamics given in eq .
( [ mastereqn ] ) . for @xmath90 , @xmath91 and @xmath92 ,
the @xmath93 equation has a solution independent of @xmath21 and @xmath40 and the atom - atom reduced density matrix is given as : @xmath94 .
the concurrence for this state is 0.50 and quantum discord of it can be calculated from eq .
( [ qd2 ] ) to be 0.412 which are seen clearly from figs .
3(p ) and  3(r ) .
also , we want to emphasize that the dynamics of @xmath95 , @xmath78 and @xmath79 between identical atoms for the initial state @xmath96 is the same as that for @xmath97 , since we assume symmetric coupling  ( i.e. , @xmath98 ) between the atoms and the single mode cavity field .
the populations @xmath99
@xmath100 and the absolute value of the coherences @xmath101 given by eq .
( [ solmaster ] ) versus @xmath70 for the initial state @xmath86 and @xmath76 for identical ( @xmath68 )  ( a ) and unidentical ( @xmath69 )  ( b ) atoms .
( c ) shows the atomic purity  ( @xmath102 ) for the density matrix eq .
( [ reducedden ] ) versus @xmath70 for identical and unidentical atoms . here
the insets are for small time regions.,title="fig:",width=283 ]   the populations @xmath99
@xmath100 and the absolute value of the coherences @xmath101 given by eq .
( [ solmaster ] ) versus @xmath70 for the initial state @xmath86 and @xmath76 for identical ( @xmath68 )  ( a ) and unidentical ( @xmath69 )  ( b ) atoms .
( c ) shows the atomic purity  ( @xmath102 ) for the density matrix eq .
( [ reducedden ] ) versus @xmath70 for identical and unidentical atoms . here
the insets are for small time regions.,title="fig:",width=283 ]     the populations @xmath99
@xmath100 and the absolute value of the coherences @xmath101 given by eq .
( [ solmaster ] ) versus @xmath70 for the initial state @xmath86 and @xmath76 for identical ( @xmath68 )  ( a ) and unidentical ( @xmath69 )  ( b ) atoms .
( c ) shows the atomic purity  ( @xmath102 ) for the density matrix eq .
( [ reducedden ] ) versus @xmath70 for identical and unidentical atoms . here
the insets are for small time regions.,width=283 ]    to further elucidate the asymptotic constant entanglement and quantum discord for the initial state @xmath86 , we display the time dependence of the density matrix elements of the atom - field system given by eq .
( [ solmaster ] ) and the purity of the atomic subsystem for the density matrix eq .
( [ reducedden ] ) for identical and unidentical detuning cases in figs .
4(a)-4(c ) for @xmath76 . as can be seen from fig .
4(a ) , all populations , except @xmath103 and @xmath104 , become zero for a short time while @xmath104 goes to @xmath105 and @xmath106 go to @xmath107 asymptotically .
the key mechanism for the creation of the entanglement and qd in this case seems to be the coherence @xmath108 generated by the interaction of the atom with the dissipative cavity .
the absolute value of @xmath108 also reaches to 1/4 asymptotically as can be seen from the fig .
the evolved state is a trapping state which has the property that the probabilities of the different decay channels interfere destructively and the state is robust against the cavity loss mechanism  @xcite . generation or protecting entanglement by population trapping has been discussed by a number of groups before  @xcite .
the magnitude of @xmath21 has no effect on the long time limit of @xmath103 and @xmath104 , it only affects how long other population components can be nonzero .
the dynamics of atomic purity  ( which is equal to 1 for a pure state and 1/4 for a maximally mixed state in four dimensional hilbert space ) displayed in fig .
4(c ) as the solid line for @xmath68 indicates that the state approaches to a partially mixed state which has high entanglement and quantum discord as indicated before in figs .
3(p ) and  3(r ) . for unidentical atoms ,
the short time dynamics is similar to the identical detuning case that all states are populated for a short time , but the asymptotic state is now @xmath109 as can be seen from fig .
4(b ) . the atomic purity ( dotted line in fig .
4(c ) ) indicates that the state of the atoms firstly goes to a mixed state in a short time and then approaches to the pure state @xmath110 which is the reason of the loss of correlations between unidentical atoms in the asymptotic limit .      in this subsection
, we will investigate the effect of initial correlations in the atom - atom part of the state on the dynamics of quantum correlations .
chsh inequality  ( a ) , concurrence  ( b ) and quantum discord  ( c ) versus @xmath70 for unidentical atoms , @xmath69 and @xmath31 where @xmath32 . here
all subfigures include plots for @xmath72 , @xmath74 , @xmath81 and @xmath76.,title="fig:",width=283 ]   chsh inequality  ( a ) , concurrence  ( b ) and quantum discord  ( c ) versus @xmath70 for unidentical atoms , @xmath69 and @xmath31 where @xmath32 . here
all subfigures include plots for @xmath72 , @xmath74 , @xmath81 and @xmath76.,title="fig:",width=283 ]     chsh inequality  ( a ) , concurrence  ( b ) and quantum discord  ( c ) versus @xmath70 for unidentical atoms , @xmath69 and @xmath31 where @xmath32 . here
all subfigures include plots for @xmath72 , @xmath74 , @xmath81 and @xmath76.,width=283 ]    we take @xmath111 as the initial state which is the maximally entangled bell state for the atomic system .
qd and bell inequality violation are also maximal for this initial state .
we display @xmath95 , @xmath78 and @xmath79 as a function of @xmath70 for several cavity decay rates for unidentical  ( fig .  5 ) and identical  ( fig .
6 ) atoms . at @xmath72 ,
all correlations show oscillatory behaviour ; for the unidentical atoms @xmath83 , @xmath82 and @xmath80 oscillate nearly in the ranges [ @xmath112 , 2.2 ] , [ 1.0 , 0.63 ] and [ 1.0 , 0.58 ] , respectively  ( fig .  5 ) . on the other hand , @xmath83 and @xmath82 for @xmath72 have periodic sudden deaths and births for the identical detunings of the atoms  ( figs .
6(a ) and  6(b ) ) , while @xmath80 oscillates nearly between 0.3 and 1  ( fig .
6(c ) ) . for small values of @xmath21  @xmath113
, the effect of cavity decay on the dynamics of correlations can be described as underdamped oscillations for both identical and unidentical detuning cases . as @xmath21
is increased to @xmath114 , the oscillation becomes overdamped and finally at large @xmath21 s  ( @xmath76 ) damping becomes critical , i.e. @xmath82 and @xmath80 show exponential decay for the unidentical atom cases .
for the identical atoms at @xmath81 and @xmath76 concurrence has sudden death , while qd decays exponentially .
one should note that bell nonlocality is the most fragile quantum correlation considered here ; it suffers sudden death in a short time as long as @xmath91 .
chsh inequality  ( a ) and  ( d ) , concurrence  ( b ) and  ( e ) and quantum discord  ( c ) and  ( f ) versus @xmath70 for identical atoms , @xmath68 and @xmath31 where @xmath32 . here ( a ) , ( b ) and ( c ) are for @xmath72 , while ( d ) , ( e ) and ( f ) include plots for @xmath74 , @xmath81 and @xmath76.,title="fig:",width=192 ]   chsh inequality  ( a ) and  ( d ) , concurrence  ( b ) and  ( e ) and quantum discord  ( c ) and  ( f ) versus @xmath70 for identical atoms , @xmath68 and @xmath31 where @xmath32 . here ( a ) , ( b ) and ( c ) are for @xmath72 , while ( d ) , ( e ) and ( f ) include plots for @xmath74 , @xmath81 and @xmath76.,title="fig:",width=192 ]   chsh inequality  ( a ) and  ( d ) , concurrence  ( b ) and  ( e ) and quantum discord  ( c ) and  ( f ) versus @xmath70 for identical atoms , @xmath68 and @xmath31 where @xmath32 . here
( a ) , ( b ) and ( c ) are for @xmath72 , while ( d ) , ( e ) and ( f ) include plots for @xmath74 , @xmath81 and @xmath76.,title="fig:",width=192 ]     chsh inequality  ( a ) and  ( d ) , concurrence  ( b ) and  ( e ) and quantum discord  ( c ) and  ( f ) versus @xmath70 for identical atoms , @xmath68 and @xmath31 where @xmath32 . here ( a ) , ( b ) and ( c ) are for @xmath72 , while ( d ) , ( e ) and ( f ) include plots for @xmath74 , @xmath81 and @xmath76.,title="fig:",width=192 ]   chsh inequality  ( a ) and  ( d ) , concurrence  ( b ) and  ( e ) and quantum discord  ( c ) and  ( f ) versus @xmath70 for identical atoms , @xmath68 and @xmath31 where @xmath32 . here ( a ) , ( b ) and ( c ) are for @xmath72 , while ( d ) , ( e ) and ( f ) include plots for @xmath74 , @xmath81 and @xmath76.,title="fig:",width=192 ]   chsh inequality  ( a ) and  ( d ) , concurrence  ( b ) and  ( e ) and quantum discord  ( c ) and  ( f ) versus @xmath70 for identical atoms , @xmath68 and @xmath31 where @xmath32 . here ( a ) , ( b ) and ( c ) are for @xmath72 , while ( d ) , ( e ) and ( f ) include plots for @xmath74 , @xmath81 and @xmath76.,title="fig:",width=192 ]
we have investigated the dynamics of quantum correlations as quantified by entanglement , bell nonlocality and quantum discord for two qubits embedded in a common leaky cavity for three different initial states with double excitations .
the qubits were considered as identical or unidentical based on the detuning of their transition frequency with respect to the frequency of the field mode .
the main findings of the study can be summarized as : qubits being identical or not have profound influence on the dynamics of the considered correlations ; starting from the same initial state , the correlations show very different behaviour depending on whether the detunings of the individual qubits satisfy @xmath92 or @xmath115 . for unidentical qubits initially in their excited states and interacting with a high or low quality cavity , quantum discord is found to be induced for all considered cavity decay rates and it decays oscillatorily for a leaky cavity although the qubits remain unentangled for all times . for the same initial state and the leaky cavity ,
qd is found to increase in the long time for identical qubits , while entanglement is either `` never created '' or decays to zero in a short time .
one of the most interesting finding is the constant high entanglement and quantum discord state which is obtained as the long - time limit of the dynamics of @xmath87 initial state for the identical qubit case .
here the final state is independent of the cavity decay rate  ( if @xmath91 ) and the magnitude of the detuning and has high concurrence  ( @xmath116 ) and quantum discord  ( @xmath117 ) . from an analysis of the time dependence of the density matrix elements
it is seen that the considered initial state leads to a trapping state which happen to have high degree of quantum correlations . on the other hand , for the same initial state
, the correlations between the unidentical atoms are found to be destroyed quickly even for a small cavity decay rate .    for the atomic bell - type initial state
we have found esd for the identical qubits interacting with a low quality cavity , while concurrence decays exponentially for the unidentical qubits .
quantum discord is found to have asymptotical decay behaviour independent of the type of atoms for this initial state and leaky cavity . on the other hand ,
chsh - bell inequality violations survive only for a finite time for @xmath91 for both identical and unidentical qubits .
it is interesting to note that steady - state nonzero correlation state of @xmath118  ( or @xmath97 ) is not observed for the @xmath111 .
this shows that the stationary states depend strongly on the initial conditions  @xcite .
we have also demonstrated that bell nonlocality as manifested by chsh inequality is the most fragile correlation measure under dissipative environment and is violated only at the high values of entanglement .