n} }$ is complex, because it [@FB90; @SC10; @CS10] depends on the properties of each individual pulsar, the instrumentation, and radio frequency interference. We therefore model the noise phenomenologically using observational knowledge. In this paper, the noise part of a pulsar timing residual is modeled as a superposition of a white noise component and a red noise component, where the white noise is designated to the measurement uncertainty on the TOA due to radiometer noise and pulse jitter noise [@LKL12]. Timing residuals of several millisecond pulsars show clear evidence of temporally correlated noise, although its origin is not yet clear [@VBC09]. The red noise components are used to empirically model such effects. We further assume that the noise components are not correlated between any two different pulsars[^3].

For each pulsar, three parameters are used to characterize the noise spectrum. These parameters are the RMS level of the white noise ${ {\,^{i}\negthinspace}}\sigma_{\rm 
w}$, the RMS level for red noise ${ {\,^{i}\negthinspace}}\sigma_{\rm r}$, and the spectral index for the red noise ${ {\,^{i}\negthinspace}}\beta$. The white noise spectrum is $${ {\,^{i}\negthinspace}}S_{\rm
 w}(f)=\frac{{ {\,^{i}\negthinspace}}\sigma_{\rm w}}{\widetilde{\sigma_{\rm w}}}\,,
 \label{eq:whitespec}$$ and the red noise spectrum is $${ {\,^{i}\negthinspace}}S_{\rm
 r}(f)=\frac{{ {\,^{i}\negthinspace}}\sigma_{\rm r} f^{{ {\,^{i}\negthinspace}}\beta}}{ \widetilde{\sigma_{\rm
                r}}}\,,
 \label{eq:redspec}$$ where the $\widetilde{\sigma_{\rm w}}$ and $\widetilde{\sigma_{\rm
 r}}$ are for normalization. From the spectrum, one can derive the correlation between noise components $$\langle { {\,^{i}\negthinspace}}n_{k} { {\,^{j}\negthinspace}}n_{k'}\rangle=\left({ {\,^{i}\negthinspace}}\sigma_{\rm
 w}^2\delta_{kk'} +{ {\,^{i}\negthinspace}}\sigma_{\rm r}^2 g_{kk'}\
n } } $ is complex, because it [ @FB90; @SC10; @CS10 ] depends on the properties of each individual pulsar, the instrumentation, and radio frequency noise. We consequently model the noise phenomenologically using experimental knowledge. In this newspaper, the noise part of a pulsar timing remainder is model as a superposition of a blank noise component and a crimson noise component, where the white randomness is designated to the measurement uncertainty on the TOA due to radiometer noise and pulse jitter noise [ @LKL12 ]. Timing residuals of respective millisecond pulsars show clear evidence of temporally correlate noise, although its origin is not so far clear [ @VBC09 ]. The crimson noise components are used to empirically model such effects. We further wear that the noise components are not correlated between any two different pulsars[^3 ]. 

 For each pulsar, three parameters are used to characterize the noise spectrum. These parameters are the RMS level of the white noise $ { { \,^{i}\negthinspace}}\sigma_{\rm 
 w}$, the RMS level for red noise $ { { \,^{i}\negthinspace}}\sigma_{\rm r}$, and the spectral index for the red noise $ { { \,^{i}\negthinspace}}\beta$. The white noise spectrum is $ $ { { \,^{i}\negthinspace}}S_{\rm 
  w}(f)=\frac { { { \,^{i}\negthinspace}}\sigma_{\rm w}}{\widetilde{\sigma_{\rm w}}}\, , 
  \label{eq: whitespec}$$ and the red noise spectrum is $ $ { { \,^{i}\negthinspace}}S_{\rm 
  r}(f)=\frac { { { \,^{i}\negthinspace}}\sigma_{\rm roentgen } f^ { { { \,^{i}\negthinspace}}\beta } } { \widetilde{\sigma_{\rm 
                 r}}}\, , 
  \label{eq: redspec}$$ where the $ \widetilde{\sigma_{\rm w}}$ and $ \widetilde{\sigma_{\rm 
  r}}$ are for normalization. From the spectrum, one can derive the correlation coefficient between randomness components $ $ \langle { { \,^{i}\negthinspace}}n_{k } { { \,^{j}\negthinspace}}n_{k'}\rangle=\left ({ { \,^{i}\negthinspace}}\sigma_{\rm 
  w}^2\delta_{kk' } + { { \,^{i}\negthinspace}}\sigma_{\rm r}^2 g_{kk'}\
n} }$ ls complex, because it [@FB90; @SC10; @CS10] depends on the 'ropertjes of exch individual pulsar, the indteumenuction, and radio freqjency intvrference. We uherefore model tis noise phenojcnoloyirally using obsgrvational ktowledge. In thhs pcper, the noise part of a pulsar timigg resicuwl is modeled ws a fupediowition of a white noise compohent anv a red noise cpmponent, where the white nlise is designated to hhe measurenent yncertainty un the TOA due to radikmeter noise and pulse jitter nuise [@KKL12]. Timing rwsifoals of sevecal mijlisecond pulsars shmw cleat evidence of bempocallt correlated noise, alvhough its origin is not yet wlzar [@VBC09]. The red noise cimpongnts dre jwed to ekpjricalpy jodel such effects. We further assume thau trv noise compohents wrq not correlated between any two differtnt phlsars[^3].

For each pulsar, tyree parameters are uded to chwracterize the noise spectrum. These parameters ara the FMS lcvel if the white noise ${ {\,^{i}\negthinspace}}\sigma_{\rm 
w}$, the WJS lvvel for red noisc ${ {\,^{i}\negthinspace}}\sibmw_{\rk r}$, and the spgctral nhdsx for the red noide ${ {\,^{i}\nedthinwpace}}\beta$. The white noise spectrum is $${ {\,^{i}\begthinspace}}F_{\em
 w}(f)=\frac{{ {\,^{i}\negthindpace}}\sigma_{\ro w}}{\wodetikde{\sigma_{\rm w}}}\,,
 \label{eq:whicespec}$$ and the ref noise sldctrum is $${ {\,^{i}\negtfinxpdce}}S_{\rm
 r}(f)=\frac{{ {\,^{i}\negthinspacq}}\sigma_{\rm c} f^{{ {\,^{i}\uegthinsoace}}\neta}}{ \wydetilde{\sihma_{\rm
                r}}}\,,
 \lwbel{ee:radspec}$$ whege the $\widetilde{\sigma_{\rm w}}$ and $\wmvetilde{\sigma_{\rk
 s}}$ age for noxmalizstion. From thq spectrum, one can dexive tfe correlanion betwxen noise coiponents $$\langng { {\,^{i}\negthinspece}}n_{k} { {\,^{j}\nqgthunspqce}}n_{k'}\ravele=\left({ {\,^{i}\negthonspace}}\siyia_{\rm
 w}^2\delta_{kk'} +{ {\,^{i}\negthikspacg}}\sjgma_{\rm r}^2 g_{kk'}\
n} }$ is complex, because it [@FB90; depends the properties each individual pulsar, interference. therefore model the phenomenologically using observational In this paper, the noise part a pulsar timing residual is modeled as a superposition of a white noise and a red noise component, where the white noise is designated to the uncertainty the due radiometer noise and pulse jitter noise [@LKL12]. Timing residuals of several millisecond pulsars show clear evidence temporally correlated noise, although its origin is not clear [@VBC09]. The red components are used to empirically such We further that noise are not correlated any two different pulsars[^3]. For each pulsar, three parameters are used to characterize the noise spectrum. These are the of the noise {\,^{i}\negthinspace}}\sigma_{\rm the RMS level noise ${ {\,^{i}\negthinspace}}\sigma_{\rm r}$, and the the red noise ${ {\,^{i}\negthinspace}}\beta$. The white noise is $${ w}(f)=\frac{{ {\,^{i}\negthinspace}}\sigma_{\rm w}}{\widetilde{\sigma_{\rm w}}}\,, \label{eq:whitespec}$$ and red noise spectrum is $${ {\,^{i}\negthinspace}}S_{\rm r}(f)=\frac{{ {\,^{i}\negthinspace}}\sigma_{\rm f^{{ {\,^{i}\negthinspace}}\beta}}{ \widetilde{\sigma_{\rm r}}}\,, \label{eq:redspec}$$ where the $\widetilde{\sigma_{\rm w}}$ and $\widetilde{\sigma_{\rm r}}$ are for normalization. spectrum, one can derive correlation between noise $$\langle {\,^{i}\negthinspace}}n_{k} {\,^{j}\negthinspace}}n_{k'}\rangle=\left({ w}^2\delta_{kk'} +{ r}^2 g_{kk'}\
n} }$ is complex, because it [@FB90; @SC10; @CS10] Depends on tHe proPerTieS oF eacH indIvidual pulsar, tHE insTrumentation, and radio frEquenCy INterFErEnce. WE therefORe MODel ThE nOisE pHEnOmenoLogIcally uSing observAtiOnAl knowledge. IN ThIs paper, the NoiSe part of a pulSar Timing ReSidUAl is mOdeLed as A superPOsitioN of a white NoISe compONent and A REd NoisE component, where thE WhITe noise is desigNated tO tHE mEASurEmeNt uncertaiNtY on thE tOA due tO RaDIOMetER noise and pulsE jitter noisE [@lKL12]. timing ReSidUAls of sEveraL mILliSecond pulsaRs shOw clear evIdence OF temporALly corrElated NoiSe, aLthoUGh ItS orIgIN is NOt Yet CLeaR [@VBC09]. The rEd NoIse coMponENTS Are uSed To emPiricAlly model such EffEcts. wE fuRther AssumE thaT tHe noiSe compOnentS aRe not correlated BetwEen any two DifFeRenT pUlsarS[^3].

for eacH puLsaR, three pArameteRS arE uSED To Characterize the noiSe SPEcTrum. ThesE paramETeRs ARe the RMS LeVel Of thE WHite nOise ${ {\,^{I}\NeGthinspaCe}}\sigmA_{\Rm 
W}$, tHe RMS leVeL for reD nOisE ${ {\,^{i}\nEgthiNSpacE}}\sigma_{\Rm r}$, and thE specTRal index for the REd noise ${ {\,^{i}\negthINsPACe}}\BEta$. THe wHite noise spEctrUM is $${ {\,^{i}\NegtHInSpaCE}}S_{\rm
 w}(F)=\frac{{ {\,^{I}\nEGtHInspace}}\sigma_{\rm w}}{\wideTiLde{\sigMa_{\rm w}}}\,,
 \Label{eq:whitesPec}$$ and the rED NOise specTrum IS $${ {\,^{i}\NEgthinspace}}S_{\rm
 R}(f)=\fraC{{ {\,^{i}\negthinsPAce}}\sigma_{\Rm r} f^{{ {\,^{i}\NegthinsPace}}\beta}}{ \wIDEtilde{\siGma_{\Rm
                r}}}\,,
 \LabEl{eQ:REdSpec}$$ where the $\wIDEtilDe{\Sigma_{\rm W}}$ anD $\widetiLde{\SigMa_{\rM
 r}}$ aRe For normalIzation. FRoM tHe SpEctRum, onE Can derivE tHe cOrRelAtion BEtween Noise CompOnEnTS $$\laNgle { {\,^{i}\neGThINSpacE}}n_{K} { {\,^{j}\NegtHinSpAce}}n_{k'}\RangLE=\leFt({ {\,^{i}\negtHinspace}}\sIgmA_{\Rm
 w}^2\dElTa_{Kk'} +{ {\,^{i}\negtHinspace}}\sigma_{\Rm R}^2 g_{kk'}\
n} }$ is complex, becauseit [@FB90;  @SC1 0;@CS 10 ] de pend s on the prope r ties  of each individual pu lsar,  t h e in s tr ument ation,a nd r adi ofr equ en c yinter fer ence. W e therefor e m od el the noise ph enomenolog ica lly using ob ser vation al  kn o wledg e.In th is pap e r, the  noise pa rt of a p u lsar ti m i ng  res idual is modeleda sa  superposition  of awh i te n ois e c omponent a nd  a re d  noisec om p o n ent ,  where the wh ite noise i s  de signat ed  to the me asure me n t u ncertaintyon t he TOA du e to r a diomete r  noiseand pu lse  ji tter no is e [ @L K L12 ] .Tim i ngresidual sof  seve ralm i l l isec ond  pul sarsshow clear ev ide nceo f t empor allycorr el atednoise,  alth ou gh its origin i s no t yet cle ar[@ VBC 09 ]. Th e  red n ois e c omponen ts areu sed  t o e mp irically model suc he f fe cts. Wefurthe r  a ss u me thatth e n oise c ompon ents ar e not co rrelat e dbe tween a ny  two d if fer ent  puls a rs[^ 3].

F or eachpulsa r , three parame t ers are usedt oc h ar a cter ize  the noisespec t rum.  The s epar a meter s are  t h eR MS level of the whi te  noise  ${ { \,^{i}\negthi nspace}}\s i g m a_{\rm w}$, th e  RMS level for  rednoise ${ { \ ,^{i}\ne gthin space}}\ sigma_{\r m r}$, and  th e s pec tra l in dex for the r e d  noi se  ${ {\, ^{i }\negth ins pac e}} \be ta $. The wh ite nois esp ec tr umis $$ {  {\,^{i} \n egt hi nsp ace}} S _{\rm w}(f )=\f ra c{ {  {\ ,^{i}\n e gt h i nspa ce }} \sig ma_ {\ rm w} }{\w i det ilde{\s igma_{\rm  w} } }\,, 
\l abel{eq :whitespec}$$  a nd the red  n ois e spec t r um is $$ { {\,^{i}\negthinspace} } S_{\rm r} (f)=\ frac {{ {\,^{i }\n egthin spa c e}}\si gma_{\ rm r}  f ^{{ { \,^{i } \ ne gth in space}}\be t a }}{  \wid et ilde {\sigma _{\rm
             r}}}\,,
 \lab el{ eq:r e d sp ec} $ $w her et he$ \ widetilde{\sigm a_{\rm w}} $a nd  $\widetil d e{\ si gma_{\r m
 r}}$  aref or norm alization . From th espec t r um,  one can d erive th e correla t ion b e tw een n ois e comp on ent s $$\ langle { { \,^{i }\negt hi nspace }}n_{ k}  { {\,^{ j}\negthinspace}}n_{k'} \rangl e=\le ft( { {\,^{i} \ne g thi nspace}}\ sigm a_{\rm
 w} ^2\ del ta_{k k'} +{ {\ ,^{i } \n egt h inspa ce}} \ sigma_{\r m  r }^2 g _{ kk'}\
n} }$_is complex,_because it [@FB90; @SC10;_@CS10] depends_on_the properties_of_each individual pulsar,_the instrumentation, and_radio frequency interference. We_therefore model the_noise_phenomenologically using observational knowledge. In this paper, the noise part of a pulsar timing_residual_is modeled_as_a_superposition of a white noise_component and a red noise_component, where_the white noise is designated to the measurement_uncertainty_on the TOA_due to radiometer noise and pulse jitter noise [@LKL12]._Timing residuals of several millisecond pulsars_show clear evidence_of_temporally_correlated noise, although its_origin is not yet clear [@VBC09]._The red noise components are used_to empirically model such effects. We further_assume that the noise components are_not correlated between any two_different pulsars[^3].

For_each pulsar, three parameters are_used to characterize_the noise_spectrum. These parameters_are the RMS level of the_white noise ${_{\,^{i}\negthinspace}}\sigma_{\rm 
w}$, the RMS level for_red_noise ${ {\,^{i}\negthinspace}}\sigma_{\rm_r}$,_and_the spectral_index for the_red_noise ${_{\,^{i}\negthinspace}}\beta$._The white noise spectrum is $${_{\,^{i}\negthinspace}}S_{\rm
_w}(f)=\frac{{ {\,^{i}\negthinspace}}\sigma_{\rm w}}{\widetilde{\sigma_{\rm w}}}\,,
 \label{eq:whitespec}$$ and the_red noise spectrum is_$${_{\,^{i}\negthinspace}}S_{\rm
 r}(f)=\frac{{ {\,^{i}\negthinspace}}\sigma_{\rm r}_f^{{ {\,^{i}\negthinspace}}\beta}}{ \widetilde{\sigma_{\rm
  _     _ _ _   r}}}\,,
 \label{eq:redspec}$$ where the $\widetilde{\sigma_{\rm w}}$ and $\widetilde{\sigma_{\rm
_r}}$ are for normalization. From the_spectrum, one can derive_the correlation_between_noise components $$\langle_{_{\,^{i}\negthinspace}}n_{k} {_{\,^{j}\negthinspace}}n_{k'}\rangle=\left({ {\,^{i}\negthinspace}}\sigma_{\rm
 w}^2\delta_{kk'} +{ {\,^{i}\negthinspace}}\sigma_{\rm r}^2 g_{kk'}\