there is much reliable evidence for the existence of dark matter ( dm ) which is mainly associated with the missing mass problem .
astronomers and particle physicists seek to probe dm properties by direct laboratory experiments or indirect observations @xcite .
indirect techniques include efforts to detect gamma rays from dm annihilation using telescopes such as the fermi large area telescope ( fermi - lat , @xcite ) , the major atmospheric gamma - ray imaging cherenkov ( magic ) telescope @xcite , the high energy stereoscopic system ( h.e.s.s . ,
@xcite ) and the very energetic radiation imaging telescope array system ( veritas , @xcite ) ( see , e.g. , @xcite for a review ) .
it was first suggested by gondolo and silk @xcite that adiabatic growth of a bh at the center of a dm halo whose density had a singular power - law cusp @xmath4 with @xmath5 led to a high density dm region around the central bh , @xmath6 with @xmath7 .
this region is called a dm spike . inside the spike ,
dm annihilations are enhanced and produce the strong gamma - ray photon flux which could be detectable to the telescopes mentioned above .    however , subsequent studies pointed out that this spike could be weakened by dynamical processes such as mergers of host - galaxies , sub - halo accretion and passing of molecular clouds @xcite .
these processes transfer energy to the dm particles and destroy the structure of the dm spike .
then the annihilation rate in the spike is smaller than predicted in @xcite because it depends on the line - of - sight integral of the squared density of the spike .
if supermassive black holes ( smbhs ) have experienced mergers , they are unlikely to have surviving spike structure .
even this however is controversial because of the uncertainty in whether the final parsec problem for smbh mergers has been resolved phenomenologically @xcite or even theoretically @xcite . on the other hand , formation scenarios of intermediate - mass black holes ( imbh ) which allow
dm mini - spikes have been proposed @xcite .
if the imbh have never experienced mergers in the past , the dm mini - spike around the central imbh is likely to survive .
imbhs may exist in our universe @xcite , and even several hundreds would reside in the halo of the present - day milky way galaxy @xcite .
those imbhs in globular clusters are recognized as promising sources for the evolved laser interferometer space antenna ( elisa ) @xcite / the new gravitational wave observatory ( ngo ) @xcite and deci - hertz interferometer gravitational wave observatory ( decigo)@xcite .    in our previous work @xcite
, we demonstrated that a very tiny effect such as the gravitational pull of a dm mini - spike around an imbh indeed affects detectability of gw by elisa and thereby we could infer presence or absence of a dm mini - spike around an imbh using gw . specifically , when a stellar mass object inspirals into the central imbh , it is affected by the gravitational force of not only the central imbh but also the mini - spike . therefore the inspiral gw is modified by the mini - spike around the central imbh .
we found that the very tiny effect from gravitational pull of a dm mini - spike could have a large impact on detectability of the gw , thanks to the huge number of orbital cycles which the binary experienced in the elisa detection frequency band .
we also found that gw detectability strongly depends on the density profile of the dm mini - spike .    in this paper
, we extend our previous work in the following way .
we again consider gws emitted from a binary system consisting of a stellar mass object and an imbh harbored in a dm mini - spike , and calculate the gw waveform including the effect of both the gravitational potential and the dynamical friction on the falling stellar mass object in the dm mini - spike .
furthermore , we investigate how accurately the dm parameters are determined by the gw observations .
we find that the dm information contained in the waveform can be extracted with very good accuracy by gw observations if the central imbh has a steep density mini - spike .
we also investigate how the detection accuracy of the dm parameters changes depending on the masses of the binary components and the density profile of the dm mini - spike such as the power index and overall normalization .
recently @xcite made clear the importance of the dynamical friction on the gw waveform in a quite different context from ours , namely , a stellar mass object falling in a compact configuration of dm clouds @xcite .
also , @xcite has given a wide survey on astrophysical environmental effects on gw signals using order of magnitude estimates , concluding that astrophysical environmental effects such as accretion disks , magnetic fields , and dm halos do not obscure gravitational wave astrophysics , e.g. , precision measurements of binary masses and tests of general relativity @xcite ( see also @xcite ) .
to indicate one exception , our paper shall clearly show that , in the recently advocated dm mini - spike scenario , environmental effects do affect gw detectability @xcite and we can measure dm properties quite accurately from elisa gw detection , which will be shown through a detailed study using a matched filtering technique and fisher matrix analysis .
we stress that while gamma - ray observations from dm annihilations can only work if the dm is a weakly interacting massive particle ( wimp ) , the gw observations we proposed should be widely applicable to any type of dm particles . furthermore
, matter is almost completely transparent to gws unlike electromagnetic waves because of the smallness of their gravitational cross - section .
hence , the gws carry pure information on the dm from the mini - spike to the detector .
future gw experiments will probe the structure of the dm mini - spike and will even offer a hint on the nature of the dm particle .
the rest of the paper is organized as follows .
[ sec : mini - halo_model ] presents the dm mini - spike model and candidates for the stellar mass object . in sec .
[ sec : gw_waveform ] we derive the gw waveform from the system which we consider and the observational errors of the waveform parameters are calculated in sec . [ sec : parameter_resolution ] . finally our conclusions are given in sec .
[ sec : conclusion ] .
we assume that the initial dm mini - halo profile which leads to the dm mini - spike after the adiabatic growth of the imbh is approximately described by the navarro - frenk - white ( nfw ) profile @xcite @xmath8 where @xmath9 is the radius , @xmath10 is the mass density and the subscript `` s '' stands for the scaling .
navarro , frenk and white obtained this profile via cosmological n - body numerical simulation and numerical fitting of obtained dm halo profiles around clusters of galaxies .
surprisingly , their simulations showed that every dm halo around a cluster of galaxies follows the nfw profile when normalized properly by @xmath11 and @xmath12 .
later work , however , shows that the inner slope may be slightly steeper than the nfw one ( @xmath13 ) and may not even be universal @xcite . in any case
, we refer to the nfw profile and the derived parameters listed in the table [ tbl : dmparameters ] below as our reference model in the following for simplicity .
we study how accurately the dark matter parameters can be measured in our reference model in detail , then extend our analysis to different sets of values of the dm parameters to take into account ambiguities in the dm distribution around an imbh .
the nfw parameters @xmath11 and @xmath12 are related to the cluster mass and concentration parameters by    @xmath14    where @xmath15 and @xmath16 is the virial radius and @xmath17 is the virial mass of the cluster , @xmath18 is the formation redshift of the cluster , @xmath19 is the matter density parameter , @xmath20 is the critical matter density of the universe and the function @xmath21 is the volume integral of the nfw profile @xmath22 ( see , e.g. , @xcite ) .
we used the fitting formula given by @xcite for the parameter @xmath23 : @xmath24 where @xmath25 @xcite .
the mass - concentration relation is taken from @xcite which fits the profiles of the clusters of galaxies obtained in their n - body simulations .
this result for clusters of galaxies may or may not apply for the mini - halo . in any case ,
concentration parameters between @xmath26 lead to qualitatively similar results and are given by the following relation .
@xmath27 where we assume @xmath28 from the result of @xcite .
the parameters @xmath29 and so on may be appropriately used when the overdensity @xmath23 equals @xmath30 .
however , we here assumed @xmath31 and so on for simplicity .
as will be shown later , the gw waveform depends on the dm mini - spike slope @xmath0 and some combination of the radius at which the mini - spike is established , @xmath32 , and the dm density there , @xmath33 . under the assumption of adiabatic growth , while the final power - law index of the dm mini - spike , @xmath0 , depends on the power - law index of the initial inner dm profile , the latter two depend on @xmath11 , @xmath12 and @xmath0 . for concreteness
, we adopt @xmath34 , @xmath35 in eqs .
( [ eq : mvir ] ) , ( [ eq : rhos ] ) , and ( [ eq : concentration ] ) @xcite and find @xmath36 , @xmath37 , and @xmath38 .
we proceed to discuss the dm profile of the mini - spike .
if the dm mini - halo initially has a cuspy profile @xmath4 with @xmath39 , then the adiabatic growth of the central imbh produces the dm mini - spike .
hence the dark matter profile becomes @xcite @xmath40 with    @xmath41    where @xmath33 is the normalization constant and @xmath32 is empirically defined by @xmath42 .
the radius @xmath43 is the distance of the gravitational influence of the central imbh with the mass @xmath44 and is approximately obtained by @xmath45 . instead of @xmath46 to estimate @xmath33 and @xmath32 .
] the slope of the dm mini - spike takes the value @xmath47 for @xmath39 . in the case of an initially nfw profile , @xmath48
, this gives rise to @xmath49 .
if the initial profile of the mini - halo is a uniform distribution , then the final profile after the adiabatic growth of the imbh would become a more gentle @xmath50 @xcite .
it is important to note that the final profile of the dm mini - halo depends on the formation history of the central imbh .
if the imbh has experienced disruptive processes such as mergers in the past , the mini - spike would be weakened or disappear .
for this reason , we do not specify the value of the power - law index @xmath0 of the dm mini - spike and treat it as a free parameter within the range @xmath51 . in the following , even if @xmath52 , we will still call the dm distribution close to the central imbh described by eq .
( [ eq : minispikedensity ] ) `` a dm mini - spike '' for the sake of simplicity .
indeed we will see that the `` dm mini - spike '' leaves its signature in the gw waveform when @xmath53 , but certainly does not when @xmath54 .
we will also assume different values of @xmath33 to study how the ambiguities in @xmath11 and @xmath12 mentioned above affect our results . finally , we take @xmath55 to be the innermost stable circular orbit of the central imbh , @xmath56
. it may be more precise to use @xmath57 @xcite , but such a change of @xmath55 does not alter at all the measurement accuracy of the dm parameters shown below .
the parameters of the dm density profile are summarized in the table [ tbl : dmparameters ] below .    .
our reference model parameters of the imbh , the dm mini - halo and the dm mini - spike .
@xmath58 : the total mass of the mini - halo , @xmath59 : the mass of the central intermediate mass black hole , @xmath60 : the formation redshift of the mini - halo , @xmath61 : the concentration of the mini - halo , @xmath16 : the virial radius of the mini - halo , @xmath12 : the nfw @xmath12 parameter of the mini - halo , @xmath11 : the nfw @xmath11 parameter of the mini - halo , @xmath43 : the radius at which @xmath62 , @xmath32 : the radius where the spike forms ( estimated by @xmath63 ) , and @xmath33 : the mini - halo mass density at @xmath32 . [ cols="^,^,^,^,^",options="header " , ]      before moving onto the calculation of the gw waveform , we discuss what can be a candidate for `` a stellar mass object '' .
let us consider a stellar mass object with mass @xmath64 denoted by @xmath65 , orbiting around an intermediate mass black hole @xmath66 with mass @xmath59 .
we consider the inspiral up to the the innermost stable circular orbit @xmath67 @xmath68 hence , the object @xmath65 should have a radius smaller than at most @xmath69 km . at the same time
, the tidal radius of @xmath65 orbiting @xmath66 at the orbital radius of @xmath67 is @xmath70 hence , this object must be either a black hole or a neutron star . alternatively ,
if we assume @xmath65 to be a white dwarf of radius @xmath71 km or a sun - like object of radius @xmath72 , the innermost orbital radius should be replaced by the radius below which the object @xmath65 is tidally destroyed @xmath73 as will be stated , we will consider the orbital radius of order @xmath74 or less , so we can not assume our stellar mass object to be a normal star with radius @xmath75 .
a white dwarf may be an interesting candidate since an electromagnetic counterpart may be expected when it is tidally disrupted ( e.g. , @xcite ) . yet
, here in this paper we assume a neutron star or a black hole when we refer to a stellar mass object .
let us consider a binary system which involves a small compact object with a mass of @xmath76 and an imbh with a mass of @xmath77 .
the mass of the stellar mass object @xmath64 is much smaller than the mass of ibmh @xmath78 .
so the reduced mass is approximately equal to @xmath64 and the barycenter position is approximately equal to the position of the imbh . by adopting a reference frame attached to the barycenter , the equation of motion of the radial relative separation between the stellar mass object and the imbh describes the motion of the former and
is given by @xmath79 where @xmath80 is the angular momentum of the stellar mass object per its mass , and @xmath81 and @xmath82 are defined by    @xmath83    the mass @xmath84 denotes the dm mass contained within the isco and is defined as @xmath85 . the first term on the right - hand side of eq .
( [ eq : eom ] ) describes the gravitational potential of the effective mass of the central imbh which is modified by the dm due to the absence of the dm within the isco , the second term accounts for the dm effect , and the third term represents a centrifugal force . here the dynamical friction force and the gw back reaction force are neglected because these effects are much smaller than the gravitational potential of the imbh .
we will introduce these effects to include an adiabatic evolution of the orbital radius in the next subsection .
we assume that the stellar mass object orbits in a circular manner for simplicity .
the orbital radius @xmath86 is obtained by solving @xmath87 in eq .
( [ eq : eom ] ) .
the orbital frequency @xmath88 is related to the angular momentum @xmath80 by @xmath89 , so we get @xmath90^{1/2}. \label{eq : orbital_frequency}\end{aligned}\ ] ] when a dm mini - spike is not present around the imbh , @xmath91 and @xmath92 , so eq . ( [ eq : orbital_frequency ] ) leads to the kepler s law @xmath93 .      in this subsection
, we introduce the gw back - reaction and the dynamical friction into the stellar mass object s orbit by taking the energy balance equation into account . when the stellar mass object orbits around the imbh , a part of its energy @xmath94 is converted into gw emission loss @xmath95 and dynamical friction loss @xmath96 .
thus the following energy balance equation is satisfied : @xmath97 as we will see in this subsection , this energy balance equation gives the time evolution of the orbital radius .
the resulting orbit can be regarded as a quasi - circular orbit because of the smallness of these dissipative effects .
the orbital energy @xmath94 is the sum of the kinetic energy and the gravitational potential of the stellar mass object , so we can calculate @xmath94 using eq .
( [ eq : orbital_frequency ] ) , @xmath98 where @xmath99 is the orbital velocity .
when we consider the evolution of the radius @xmath86 , @xmath100 does not vanish .
so the time derivative of eq .
( [ eq : orbital_energy ] ) gives the following equation , @xmath101 to the lowest order in the post newtonian expansion , the gravitational radiation energy is given by the quadrupole formula .
we apply the formula to the circular newtonian binary and obtain @xmath102 when the stellar mass object moves through the cloud of dm , it gravitationally interacts with dm particles .
this effect is called dynamical friction , sometimes referred to as gravitational drag which was first discussed by chandrasekar @xcite .
because of dynamical friction , the stellar mass object running through the dm halo is decelerated in the direction of its motion and loses its kinetic energy as well as its angular momentum .
the dynamical friction force is given by @xmath103 where @xmath99 is the velocity of the stellar mass object @xcite .
the coulomb logarithm @xmath104 is defined by @xmath105 where @xmath106 is the maximum impact parameter and @xmath107 is the typical velocity of the stellar mass object .
we take @xmath108 . from the expression of the dynamical friction force
, we obtain the friction loss , @xmath109    to find the numerical solution of the energy balance equation ( [ eq : energy_balance_eq ] ) easily , we introduce a dimensionless radius parameter @xmath110 defined by @xmath111 with @xmath112 using the above definition of @xmath110 , the energy balance equation ( [ eq : energy_balance_eq ] ) can be rewritten in the form of the differential equation of @xmath110 with respect to time @xmath113 as @xmath114 }        - c_{\text{df}}\dfrac{1}{\left(1+x^{3-\alpha}\right)^{1/2}\left[1+\left(4-\alpha \right)x^{3-\alpha}\right]x^{- 5/2 + \alpha } } , \label{eq : dxdt}\end{aligned}\ ] ] where the coefficients are defined by    @xmath115}. \label{eq : c_df}\end{aligned}\ ] ]    the coefficient @xmath116 is related to the gravitational radiation energy and the coefficient @xmath117 is related to the dynamical friction . in the case of the initially nfw profile , @xmath118 , the coefficients @xmath116 and @xmath119 are @xmath120 , c_{\text{df } } = 2.1 \times 10^{-8 } \ \left[1/\text{year}\right ] $ ] .
note that the dynamical friction coefficient @xmath117 is much larger than the gravitational radiation coefficient @xmath116 .
the gw waveform from the binary composed of two compact objects with masses @xmath64 and @xmath59 is given by    @xmath121    where @xmath122 is the distance to the source ( luminosity distance for a cosmological source ) , @xmath86 is the orbital radius , @xmath123 is the inclination angle , which is the angle between the line - of - sight and the rotational axis of the orbits , and @xmath124 is the gw frequency which is given by @xmath125@xcite .    the waveforms eqs .
( [ eq : gwwaveform_plus_without_backreaction ] ) and ( [ eq : gwwaveform_cross_without_backreaction ] ) are derived on the assumption that the motion of the source is described by a circular newtonian orbit .
but in fact , the radius @xmath86 and the frequency @xmath88 are not constant because the orbital energy @xmath94 decreases gradually due to both dynamical friction and the gw back - reaction . including these effects
, the orbit shrinks adiabatically and becomes a quasi - circular orbit .
so the radius @xmath86 and the frequency @xmath88 should be replaced by @xmath126 , @xmath127 and the phase @xmath128 should also be replaced by @xmath129 as defined by eq .
( [ eq : phase_t_def ] ) below .
thus , the gw waveform is expressed by    @xmath130 , \label{eq : gwwaveform_plus_including_backreaction } \\   h_{\times}\left(t\right )    & = \dfrac{1}{\sourcedistance } \dfrac{4g\mu \omega_s(t)^2 r(t)^2}{c^4 }   \cos \iota \sin\left [ \phi(t)\right ] , \label{eq : gwwaveform_cross_including_backreaction } \\
\phi \left(t\right )    & \equiv \int^t \omega_{\text{gw } } \left(t'\right ) \ dt ' .
\label{eq : phase_t_def}\end{aligned}\ ] ]    in order to discuss detectability and parameter accuracy in gw observations , it is convenient to work in the frequency domain
. the fourier transformation of the gw waveform is given by @xmath131 where @xmath132 is the gw frequency .
for simplicity , we consider a gw coming in the detector from the optimal direction for @xmath133 mode .
in such a situation , detector pattern function are @xmath134 and @xmath135 .
so the response of the detector to the gw is @xmath136 . using eq .
( [ eq : gwwaveform_plus_including_backreaction ] ) , we rewrite the gw waveform as    @xmath137    where @xmath138 is the time - dependent amplitude and @xmath139 is the time - dependent gw phase . in the above equations ,
we have introduced the retarded time @xmath140 . in the range of frequency
we are concerned with , the time - dependent amplitude @xmath141 varies slowly , while the time - dependent phase @xmath139 varies rapidly .
so , the fourier transformation of the gw waveform can be calculated approximately using the stationary phase method . in this method ,
the rapidly oscillating term is neglected and only the slowly oscillating term survives .
then the gw waveform in the fourier domain becomes    @xmath142^{1/2 } , \label{eq : waveformf}\\   & \psi(t )   =   2\pi f \dfrac{d}{c } + \tilde{\phi}(t ) - \dfrac{\pi}{4 } , \label{eq : psit_def}\\   & \tilde{\phi}(t ) \equiv   2\pi ft - \phi(t ) , \label{eq : tilde_phi_def}\end{aligned}\ ] ]    where the time @xmath113 is related to frequency by @xmath143 .
as we will discuss in appendix [ app : rewriting_waveform ] , the gw waveform eqs .
( [ eq : waveformf ] ) , ( [ eq : psit_def ] ) and ( [ eq : tilde_phi_def ] ) can be rewritten explicitly in the frequency domain as follows :    @xmath144^{-1/2 } , \label{eq : gwwaveform01}\\    & \mathcal{a}= \left(\dfrac{5}{24 } \right)^{1/2 } \dfrac{1}{\pi^{2/3 } }   \dfrac{c}{d } \left ( \dfrac{gm_c}{c^3}\right)^{5/6 }   \dfrac{1+\cos^2\iota}{2 } , \label{eq : gwwaveform02}\\    & \psi\left(f \right ) = 2\pi f\tilde{t}_c - \phi_c -\dfrac{\pi}{4 } - \tilde{\phi}\left(f \right ) , \label{eq : gwwaveform03}\\    & \tilde{\phi } \left(f\right ) = \dfrac{10}{3 }   \left ( \dfrac{8\pi gm_c}{c^3}\right)^{- 5/3 } \left [     - f \int_\infty^{f } df ' \
\dfrac{\chi^{11/2}}{f'{}^{11/3 } k \left ( 1 + \tilde{c } j \right )   }      + \int_\infty^{f } df ' \
\dfrac{\chi^{11/2 } } { f'{}^{8/3}k \left ( 1 + \tilde{c } j \right ) } \right ] , \label{eq : gwwaveform04}\\    & j\left(x\right ) = \dfrac{4x^{11/2-\alpha}}{\left ( 1 + x^{3-\alpha}\right)^{7/2 } } , \label{eq : gwwaveform05}\\    & k\left(x\right ) = \dfrac{\left ( 1 + x^{3-\alpha } \right)^{5/2 } \left(1+\alpha x^{3-\alpha}/3 \right)}{1 + \left(4-\alpha\right)x^{3-\alpha } } , \label{eq : gwwaveform06}\\    & \chi = \left ( \delta \varepsilon \right)^{1/(\alpha-3 ) } x , \label{eq : gwwaveform07}\\    & \delta = \left(\dfrac{gm_{\text{eff}}}{\pi^2 f^2 } \right)^{(3-\alpha)/3 } , \label{eq : gwwaveform09}\end{aligned}\ ] ]    where @xmath145 is the overall amplitude , @xmath146 is the chirp mass defined by @xmath147 , @xmath148 is the sum of the binary coalescence time @xmath149 and @xmath150 , @xmath151 is the phase at coalescence , @xmath0 is the power - law index of the dm mini - spike , @xmath152 is defined by @xmath153 , @xmath154 is a new frequency variable , @xmath110 is defined in eq .
( [ eq : def_x ] ) , and @xmath155 is defined by eq .
( [ eq : def_epsilon ] ) .
the dm information is encoded in the waveform eq .
( [ eq : gwwaveform11 ] ) through @xmath156 , @xmath157 , @xmath158 , @xmath152 and @xmath159 .
so if we take @xmath160 , @xmath161 , @xmath162 , @xmath163 , then eq . ( [ eq : gwwaveform11 ] )
becomes the waveform without the dm shown in eqs .
( [ eq : gw_without_dm_01])-([eq : gw_without_dm_04 ] ) .
as we will discuss in the next section , we consider a five - year observation by elisa which corresponds to @xmath165 . in this setup , @xmath166 is satisfied .
for example , we get @xmath167 for @xmath168 , @xmath169 , @xmath170 and the parameters @xmath33 , @xmath32 , and @xmath59 listed in the table [ tbl : dmparameters ] .
so @xmath171 can be treated as a small expansion parameter .
since the measurement errors of the physical parameters contained in gw are much more sensitive to the gw phase rather than its amplitude , we expand the gw waveform @xmath172 up to the first order in @xmath173 in the phase and up to the zero - th order in the amplitude . using an expansion of @xmath158 in @xmath173 , @xmath174 the gw waveform given by eqs .
( [ eq : gwwaveform01])-([eq : gwwaveform09 ] ) becomes    @xmath175 , \label{eq : gwwaveform13}\\    & l\left(f\right ) = 1 + 4c_\varepsilon \tilde{\delta}^ { ( 11 - 2\alpha ) /   \left[2\left(3-\alpha\right)\right ] } , \label{eq : gwwaveform14 } \\    & \tilde { \delta }   = \left(\dfrac{g}{\pi^2 f^2 } \right)^{(3-\alpha)/3 } , \label{eq : gwwaveform15}\\    & c_\varepsilon   = \meff^{\left(11 - 2\alpha\right)/6 } \tilde{c } \varepsilon^{(11 - 2\alpha)/\left[2\left(3-\alpha\right)\right ] } ,   \label{eq : gwwaveform16}\end{aligned}\ ] ]    where the overall amplitude @xmath145 is defined by eq . ( [ eq : gwwaveform02 ] ) and @xmath176 is defined by eq .
( [ eq : gwwaveform03 ] ) .
the upper bound of the integration in eq .
( [ eq : gwwaveform13 ] ) , @xmath177 in the elisa frequency band , is the gw frequency when the stellar mass object enters the innermost stable circular orbit .
hence , @xmath178 is in essence the gw cycles from the frequency @xmath132 to the coalescence .
the post - newtonian ( pn ) effects which are neglected in the above equations must be taken into account in real data analysis .
however , the frequency - dependence of the pn effect in the gw phase @xmath179 differs from that of the dm effect which depends on the power - law index @xmath0 .
so , the measurement accuracies of the dm parameters as we will discuss later would not be affected seriously by higher order terms in the pn expansion .    until the previous sections
, we have included both the dynamical friction and the gravitational pull of the dm mini - spike .
it is easily shown that the dynamical effect has much more impact on the measurement accuracy of the dm parameters than the dm mini - spike does @xcite , and the above expression indeed includes the dynamical friction but not the gravitational pull of the dm mini - spike .
in fact , within the approximation in this subsection and the following , the gravitational potential of the dm mini - spike shows its signature only in the imbh mass redefinition ( @xmath180 in @xmath146 in the above equations ) .
we however note that even such a tiny effect as the gravitational pull of the dm mini - spike do affect the detectability of gw thanks to the large number of the gw cycles in the elisa detection band @xcite .
it is important to note that the dm parameters appear only in @xmath0 and @xmath181 and that they are contained in the gw phase @xmath182 .
we make use of the above equations to calculate measurement errors of the waveform parameters in the next section .
we also define the phase difference @xmath183 by @xmath184 where @xmath182 defined by eq .
( [ eq : gwwaveform13 ] ) is the phase including the dm effect and @xmath185 defined by eq .
( [ eq : gw_without_dm_04 ] ) is the phase without the dm effect .
@xmath186 is shown in fig .
[ fig : phase_diff ] which indicates that the phase difference becomes significant for large @xmath0 and for the large gw frequency @xmath132 .
this is because in this case , the dm density near the central bh increases and the effect of the dm on the motion of the stellar mass object is significant . as we discussed in our previous paper @xcite
, the phase difference causes the mismatch between the waveform including the dm effect and the waveform without the dm effect .
the phase difference @xmath187 typically above 1 indicates the necessity to use the waveform including the dm effect as a template . as can be seen in fig .
[ fig : phase_diff ] , if the template without the dm effect is applied to the gw signal including the effect induced by the dm with @xmath188 , the resulting @xmath189 would degrade significantly .     against the power - law index @xmath0 , defined by eq .
( [ eq : phase_difference ] ) .
in essence , this is the difference between the accumulated phase from gw frequency @xmath132 and the binary coalescence with and without the dm mini - spike .
three different curves show @xmath190 for three different values of @xmath0 .
for instance , if detecting a binary gw from @xmath191hz to its coalescence , we would observe by a factor of @xmath192 more gw cycles in the case with a @xmath49 dm mini - spike than without any .
for this plot , we take @xmath170 and @xmath33 , @xmath32 , and @xmath59 are as listed in the table [ tbl : dmparameters ] . , width=302 ]
in this subsection , we give a brief review of parameter estimation ( see @xcite for more details ) .
let us consider detecting gws with a single detector .
the detector output @xmath193 can be written by the sum of the gw signal @xmath194 and detector noise @xmath195 : @xmath196 assuming that the detector noise is stationary , the correlation between different fourier components of the noise is expressed as @xmath197 where the angled brackets @xmath198 denote an ensemble average , the asterisk is complex conjugation and @xmath199 is a one - sided power spectral density of the detector noise . in this paper , we consider the gw observation using elisa which has the best sensitivity at around @xmath200 $ ] .
the noise spectral density of elisa is given by    @xmath201 ^ 2 , \label{eq : elisa_sn}\end{aligned}\ ] ]    where @xmath202 $ ] is the acceleration noise spectral density , @xmath203 $ ] is the shot noise spectral density , @xmath204 $ ] is the other measurement noise spectral density and @xmath205 $ ] is the separation between the spacecraft which is the length of its arms of the laser interferometer ( see @xcite for details ) .
it is convenient to introduce a noise - weighted inner product between two signals @xmath206 and @xmath207 by @xmath208 where @xmath209 denotes the real part and @xmath210 is the initial frequency . assuming that the detector noise is gaussian and stationary , the probability density of the detector noise is described by @xmath211 .
we can rewrite this expression in the form of detector signal @xmath193 and gws signal @xmath194 using eq .
( [ eq : detector_output ] ) as @xmath212 .    in the above case , @xmath194
is known , while in actual gw experiments , @xmath194 should be replaced with a template @xmath213 , where @xmath214 is a collection of unknown parameters . to determine the waveform parameters @xmath215 , it is necessary to search for the parameters which minimize the logarithm of the maximum likelihood ratio , @xmath216 . as a result of this process , we can infer the values of @xmath215 .
however , the expected values have statistical errors because the detector noise is a random process . these measurement errors @xmath217 of the waveform parameters are approximately described by the gaussian probability distribution for large @xmath189 , @xmath218 where @xmath219 is the normalization factor and @xmath220 is called the fisher information matrix defined by @xmath221 the inverse of the fisher matrix gives the root - mean - square ( rms ) errors of the waveform parameters @xmath222 : @xmath223 where @xmath224 denotes the diagonal elements of the inverse fisher matrix .
the inspiral gw waveform from the imbh surrounded by the dm mini - spike is described by six parameters which appear in eqs .
( [ eq : gwwaveform11 ] ) - ( [ eq : gwwaveform16 ] ) : the overall amplitude , @xmath145 ; the time constant , @xmath225 , which is the sum of the traveling time @xmath150 and the coalescence time @xmath149 ; the coalescence phase , @xmath151 ; the chirp mass , @xmath146 ; the two dm parameters , @xmath0 and @xmath181 .
note that the beam pattern function of elisa is neglected here because we are concerned with how the dm parameters are determined by gw observations but not with the angular resolution of elisa ( see @xcite for discussion of angular resolution ) .
the inner product between the derivatives of the waveform with respect to the parameters @xmath226 yields the values of the fisher matrix elements .
the derivatives with respect to @xmath227 and @xmath228 are calculated straightforwardly as follows :    @xmath229    the derivatives with respect to the dm parameters @xmath230 are obtained by applying the chain rule to the following equations :    @xmath231    where @xmath232 is defined in eq .
( [ eq : gwwaveform14 ] ) .
however , since the explicit expressions are complicated , we take the derivatives numerically .    next
, we derive the initial frequency at which the gw observation starts . in the presence of the dm mini - spike
, the stellar mass object orbiting the central imbh loses its angular momentum gradually due both to the dynamical friction and gw radiation reaction .
so the coalescence arises earlier than the case without the dm .
time evolution of the frequency is described by eq .
( [ eq : dfdt ] ) , @xmath233 ,   \label{eq : dfdtau}\end{aligned}\ ] ] where @xmath234 and @xmath235 is the time to the coalescence .
( note that @xmath236 . )
the lower bound @xmath237 of the integral in eq .
( [ eq : inner_product_def ] ) is required for calculating the inner product in the fisher matrix .
given that the gw is observed by elisa for 5 years prior to the coalescence , this bound is obtained by @xmath238\right).\end{aligned}\ ] ] by numerically solving eq .
( [ eq : dfdtau ] ) , we show the dependence of @xmath0 on @xmath210 in fig .
[ fig : fini ] .
this figure indicates that the dm mini - spike affects more strongly the motion of the stellar mass object for larger @xmath0 .
the initial frequency for a 5 year observation is almost constant for small @xmath0 due to the smallness of the effect of the dm .
conversely , the initial frequency drops sharply for large @xmath0 due to the dynamical friction from the dm .    .
we assume that the gw is detected by an elisa 5 year observation . for small @xmath0
, @xmath210 is almost constant . on the other hand , for large @xmath0
, @xmath210 drops sharply due mainly to the dynamical friction . for this plot ,
we take @xmath170 and @xmath33 , @xmath32 , and @xmath59 are as listed in the table [ tbl : dmparameters].,width=302 ]      in this section , we consider the most likely case that the dm mini - spike has an initially nfw profile , @xmath49 .
derivatives of the waveform with respect to the parameters given by eqs .
( [ eq : waveform_derivative_01])-([eq : waveform_derivative_06 ] ) are calculated numerically .
substitution of these results into eqs .
( [ eq : fisher_matrix_def ] ) and ( [ eq : rms_error ] ) gives rise to the rms errors @xmath239 , @xmath240 , @xmath241 , @xmath242 , @xmath243 and @xmath244 as follows .    @xmath245 \left(\dfrac{10}{s / n}\right ) ,   \\
\delta \phi_c                            & = 1.3 \left[\text{rad}\right ] \left(\dfrac{10}{s / n}\right ) ,   \\   \dfrac{\delta m_c}{m_c }                  & = 3.1 \times 10^{-7 } \left(\dfrac{10}{s / n}\right ) , \label{eq : chirpmass_error}\\   \dfrac{\delta \alpha}{\alpha }            & = 1.2 \times 10^{-6 } \left(\dfrac{10}{s / n}\right ) , \label{eq : alpha_error}\\   \dfrac{\delta c_\varepsilon}{c_\varepsilon } & = 5.9 \times 10^{-5 } \left(\dfrac{10}{s / n}\right ) .
\end{aligned}\ ] ]    here we take @xmath33 , @xmath32 , and @xmath59 from table [ tbl : dmparameters ] and @xmath170 .
these measurement errors are inversely proportional to @xmath189 .
so the waveform parameters are measurable with better accuracy for larger gw signals .
a notable feature of the above results is that the chirp mass @xmath146 and the two dm parameters @xmath0 and @xmath181 are determined much more accurately than the overall amplitude @xmath145 , the coalescence time @xmath148 and the coalescence phase @xmath151 .
this fact reflects that @xmath146 , @xmath0 and @xmath181 appear in the phase of the waveform @xmath182 . from eqs .
( [ eq : number_of_circles ] ) , ( [ eq : phase_t_def ] ) and ( [ eq : tilde_phi_def ] ) , the gw phase is proportional to the number of gw cycles which amplify the sensitivity to the parameters which appear in the phase @xmath182 by a factor @xmath246 .
thus , the fractional error of the chirp mass which is proportional to the phase is order of @xmath247 and the two dm parameters are also determined very accurately .
in fact , fig . [
fig : ncycle ] indicates the value of @xmath247 is about @xmath248 , which is consistent with the value of @xmath249 in eq .
( [ eq : chirpmass_error ] ) .
we also investigate the correlation between the parameters which appear in the phase , @xmath146 , @xmath250 and @xmath181 . figure .
[ fig : fisher_ellipse ] illustrates the fisher ellipses for @xmath251 and @xmath181 in @xmath252 . from the figures , we observe that @xmath251 and @xmath181 are strongly correlated with each other because all of them are contained in the phase .
however , they are not completely degenerate and are determined independently .
this fact can be traced to the difference of the frequency - dependence between @xmath251 and @xmath181 .     in the case where the initial dm halo has an nfw profile and
the final profile has the radial power - law index of @xmath49 through an adiabatic growth .
we assume @xmath33 , @xmath32 , and @xmath59 from the table [ tbl : dmparameters ] and @xmath170 . , width=566 ]    in the above discussion , the mass of the central imbh @xmath78 and that of the stellar mass object @xmath64 are fixed .
next we analyze the measurement errors for various values of @xmath64 and @xmath78 .
the results are shown in fig .
[ fig : error_virsus_mbh_mu ] .
the figure indicates the errors of the parameters in the phase @xmath182 increase linearly with the stellar mass object mass @xmath64 .
this behavior comes from the fact that the number of cycles @xmath253 decreases in proportion to the stellar mass object mass @xmath64 .
similarly , the larger is the mass of the imbh , the smaller the number of the orbital cycles the stellar mass object experienced in the five years prior to the coalescence within the elisa band .
for this reason , the measurement errors in @xmath146 , @xmath0 , and @xmath254 increase for a larger imbh mass as can be seen in fig .
[ fig : error_virsus_mbh_mu ] .     versus ( a ) the central bh mass @xmath78 and ( b ) the stellar mass object mass @xmath64 for @xmath255 and @xmath49 . for this plot , @xmath33 and @xmath32
are taken from the table [ tbl : dmparameters ] .
the other parameter is fixed to be @xmath256 in the left and @xmath257 in the right , respectively .
note that the both axes are in the logarithmic scales .
the solid line , the dashed line , the dashed - dotted line correspond to @xmath258 respectively . , width=680 ]      we now extend the analysis in the previous section where we considered the case of the initially nfw profile .
we next consider the general case without specifying the value of @xmath0 with @xmath59 , @xmath33 , and @xmath32 fixed to the values quoted in the table [ tbl : dmparameters ] .
the rms errors depend on the dm power - law index @xmath0 .
we show @xmath242 , @xmath243 and @xmath259 in fig .
[ fig : error_versus_alpha ] .    as shown in fig .
[ fig : error_versus_alpha ] , the accuracy of the dm parameters @xmath243 and @xmath260 are better for the larger @xmath0 .
this is because steeper density distributions contain more dm mass within the orbital radius ( see fig . 2 in @xcite ) .
in other words , the steeper density distribution has more impact on the motion of the stellar mass object and the gw waveform is modified more strongly by the dm mini - spike .
so the dm information can be extracted from the gw waveform if the dm mini - halo near the bh has a steep profile .
on the other hand , the measurement accuracy of the parameters which appear in the phase @xmath261 become worse in @xmath262 .
this feature can be explained by the number of gw cycles @xmath253 which will be discussed in the appendix [ app : ncycle ] .
there we show that @xmath253 falls sharply at @xmath263 ( see fig .
[ fig : ncycle ] ) .
the sensitivity to the parameters which appear in the phase @xmath261 is amplified by the number of circles @xmath253 in the frequency bandwidth of elisa .
for this reason , the measurement errors of @xmath251 and @xmath264 increase suddenly at @xmath263 , as is shown in fig . [
fig : error_versus_alpha ] .
we also note that this figure shows that we can measure the power - law index @xmath0 at 10 % level even for a moderately flat radial distribution with @xmath265 .
in fact , when considering the gravitational pull due to the dm potential only , it affects detectability of gw signals only for @xmath266 .
it is the dynamical friction that enables us to explore a flatter dm distribution than `` a dm mini - spike '' referred in the literature that has @xmath267 .
figure [ fig : error_versus_alpha_varrho ] shows the relative errors of the dm parameters , @xmath0 and @xmath181 for various values of @xmath268 as a function of @xmath0 . as can be seen in fig .
[ fig : error_versus_alpha_varrho ] , the relative errors for the fixed @xmath0 become smaller approximately linearly as the dm density increases .
this behaviour can be traced to the amount of the dm within the orbital radius of the stellar mass object .
it should be noted that the value of @xmath268 we adopt in this paper is derived under the assumption that the initial dm mini - halo profile is the nfw profile as discussed in the section [ sec : mini - halo_model ] .
even if the dm density is an order of magnitude more sparse than that indicated by the nfw profile , the power - law index @xmath0 is measurable with an accuracy of @xmath269 for @xmath270 .     versus the power - law index of the dm profile for @xmath252 in the case where the dm mini - spike harboring the imbh has a radially power - law profile .
the solid line , the dashed line , and the dashed - dotted line corresponds to @xmath271 respectively . for this plot , @xmath170 and the values of the parameters @xmath59 , @xmath33 , and @xmath32
are assumed as in the table [ tbl : dmparameters ] .
, width=302 ]     and ( b ) @xmath181 versus the power - law index of the dm profile for @xmath272 in the case where the dm mini - spike harboring the imbh has a radially power - law profile . the solid line , the dashed line and the dashed - dotted line corresponds to @xmath268 , @xmath273 , and @xmath274 respectively .
the value of @xmath268 is taken from the table [ tbl : dmparameters ] .
, width=680 ]
in this paper , we have investigated the measurement accuracy of dark matter ( dm ) parameters by gravitational ( gw ) observations .
we consider a binary system composed of an intermediate mass black hole ( imbh ) surrounded by a dm mini - spike and a stellar mass compact object .
the compact object falling into the central imbh is affected by the gravitational interaction of both the imbh and the dm mini - spike , namely the gravitational potential of both the imbh and the dm mini - spike , gravitational wave back - reaction and dynamical friction .
then the resulting inspiral gw is modified by the dm mini - spike in comparison with the case where the imbh has no dm mini - spike around it .
such a gw will be detected by future space - crafted detectors such as elisa / ngo .
we find that thanks to the dm parameters contained in the gw phase , the measurement errors of the dm parameters are very small for large power - law index of the mini - spike profile . to put it another way
, we can extract the dm parameters very accurately from the gw waveform using matched filtering if the dm mini - spike has a steep profile .
indeed , in our reference case as originally advocated by @xcite , we could determine the power - law index of the dm mini - spike radial profile with the 1 @xmath1 relative error of @xmath275 for a gw signal with signal - to - noise - ratio 10 and assuming 5 years observation with elisa , as shown in fig .
[ fig : fisher_ellipse ] and indicated by eq .
( [ eq : alpha_error ] ) .
we also investigated how accurately the dm parameters can be determined for various dm parameters and the masses of the imbh - stellar mass object binary surrounded by a dm mini - spike .
we have found that smaller the mass of the stellar mass object , that of the imbh , or the larger the power - law index of the dm mini - spike , we can measure dm parameters to better accuracy as shown in figs . [ fig : error_virsus_mbh_mu ] and [ fig : error_versus_alpha ] .
even a moderately flatter mini - spike with the radial distribution proportional to @xmath276 would still allow us to determine the power - law index to 10 % accuracy .
indirect dark matter searches in the gamma - ray band and through gw observation proposed in our previous @xcite and current papers are complementary to each other .
the gw observations we propose should be applicable to both very weakly annihilating and non - annihilating dm particles . even if the dm particles do not weakly interact with each other
, they affect the motion of the stellar mass object gravitationally and the resulting gw is modified by them .
gw are insensitive to absorption and scattering in the interstellar medium during the propagation unlike electromagnetic waves .
therefore gw observations offer information on the dm mini - spike directly . on the other hand ,
if the dm particles self - annihilate , an annihilation plateau may develop within a hubble time @xcite and the power - law index @xmath0 of the dm radial profile becomes effectively zero within a radius @xmath277 . for the case of the values of the parameters listed in the table [ tbl : dmparameters ] , the dm mass @xmath30gev , and its cross section @xmath278 @xcite , we find @xmath279pc which is much larger than the initial orbital radius at which the gw frequency from the binary enters the elisa detection frequency band .
hence in this case , gamma - ray searches are a better way to explore a dm mini - halo surrounding the imbh , as indicated by fig .
[ fig : error_versus_alpha ] .    in summary ,
the combination of gamma - ray observations with future gw observations will enable us to probe the structure of the dm mini - spike and even to offer hints that may clarify the nature of dm particles .
moreover , because the dm profile strongly depends on the formation history of the central imbh , both types of observation may shed light on how the imbh evolved with cosmic history .
the authors thank junichi yokoyama and enrico barausse for useful comments .
this work is supported by the grant - in - aid for jsps fellows
26@xmath2808636 ( ke ) , the grant - in - aid for young scientists no .
25800126 ( yi ) and the mext grant - in - aid for scientific research on innovative areas `` new developments in astrophysics through multi - messenger observations of gravitational wave sources '' ( grant number 24103005 ) ( yi ) .
our goal in this appendix is to rewrite the waveform eq .
( [ eq : waveformf ] ) in the form of an explicit function of gw frequency @xmath132 . according to eqs .
( [ eq : waveformf ] ) , ( [ eq : psit_def ] ) and ( [ eq : tilde_phi_def ] ) , the gw waveform from the binary system composed of the stellar mass object and the imbh surrounded by dm mini - halo is expressed by    @xmath281^{1/2 } , \label{eq : waveformf_appendix}\\   & \psi \left(t\right ) =   2\pi f   \dfrac{\sourcedistance}{c } + \tilde{\phi } \left(t
\right ) - \dfrac{\pi}{4 } , \label{eq : psit_def_appendix}\\   & \tilde{\phi } \left(t\right ) \equiv 2\pi ft - \phi \left(t\right ) , \label{eq : tilde_phi_def_appendix}\end{aligned}\ ] ]    where @xmath141 is the time - dependent amplitude defined by eq .
( [ eq : amplitude_t_def ] ) , @xmath139 is the time - dependent phase defined by eq .
( [ eq : phase_t_def ] ) , and @xmath122 is the distance to the source . we proceed as follows .
we start with the amplitude @xmath282 .
the frequency @xmath132 can be expressed in terms of the orbital radius @xmath86 which is related to the time @xmath113 by eq .
( [ eq : dxdt ] ) .
so the amplitude can be expressed as a function of the frequency @xmath132 through the relation between @xmath113 and @xmath132 .
next , we tackle the phase @xmath283 . finally ,
combining these results , we find the explicit expression of the gw waveform in the fourier domain .
the gw frequency @xmath284 which is defined by eq .
( [ eq : orbital_frequency ] ) is expanded in a taylor series in a power of @xmath86 : @xmath285^{1/2 } \nonumber\\    & = \dfrac{\sqrt{g m_{\text{eff}}}}{\pi } r^{-3/2}\left [ 1 +      \dfrac{1}{2 } r^{3-\alpha } \varepsilon -\dfrac{1}{8 }      r^{2(3-\alpha ) } \varepsilon^2 + \cdots \right].\end{aligned}\ ] ] inverting this equation , we obtain @xmath86 as a function of gw frequency and expanded in @xmath155 :    @xmath286 , \\ & \delta \equiv \left(\dfrac{gm_{\text{eff}}}{\pi^2 f^2 } \right)^{(3-\alpha)/3 } , \label{eq : delta_def } \end{aligned}\ ] ]    where we introduce a new frequency variable @xmath154 defined by eq .
( [ eq : delta_def ] ) for convenience . using the definition of @xmath110 given by eq .
( [ eq : def_x ] ) , the dimensionless radius parameter @xmath110 can be expanded in a power of @xmath155 :    @xmath287    where we introduce @xmath158 for convenience .
note that the function @xmath158 is equal to one when a dm mini - spike is not present around an imbh .    for later convenience , we rewrite @xmath288 which is defined by eq .
@xmath289 as follows .
@xmath290 } - c_{\text{df } } \dfrac{1}{\left(1+x^{3-\alpha } \right)^{1/2}\left [ 1+\left(4-\alpha \right)x^{3-\alpha } \right]x^{-5/2 + \alpha } } \nonumber\\    & = - c_{\text{gw } } f_{\text{gw}}\left ( x\right ) -   c_{\text{df } } f_{\text{df}}\left ( x\right ) \nonumber\\    & = - c_{\text{gw } }   f_{\text{gw}}\left(x\right ) \left [ 1 + \dfrac { c_{\text{df } } } { c_{\text{gw } } } \dfrac { f_{\text{df}}\left(x\right ) } { f_{\text{gw}}\left(x\right ) } \right ] \nonumber\\    & = - c_{\text{gw } }   f_{\text{gw}}\left(x\right ) \left [ 1 + \tilde{c } j\left(x\right ) \right ] , \label{eq : dxdt_x}\end{aligned}\ ] ] where functions @xmath291 and @xmath157 and a coefficient @xmath152 are defined by    @xmath292},\\    & f_{\text{df}}\left(x\right )      \equiv \dfrac{1}{\left(1+x^{3-\alpha } \right)^{1/2}\left [ 1+\left(4-\alpha \right)x^{3-\alpha } \right]x^{-5/2 + \alpha } } , \\    & j\left(x\right ) \equiv      \dfrac { f_{\text{df}}\left(x\right ) } { f_{\text{gw}}\left(x\right ) } = \dfrac{4x^{11/2-\alpha}}{\left ( 1 + x^{3-\alpha}\right)^{7/2 } } \label{eq : j_def},\\    & \tilde{c } \equiv
\dfrac { c_{\text{df } } } { c_{\text{gw}}}.\end{aligned}\ ] ]    the coefficient @xmath152 is the ratio of the dynamical friction coefficient to the gravitational radiation coefficient .
so @xmath152 includes the dm information .
next , we rewrite the second time derivative of @xmath293 , @xmath294 , as a function of @xmath110 . from eq .
( [ eq : phase_t_def ] ) , @xmath294 is expressed by @xmath295 } \dfrac{3+\alpha x^{3-\alpha}}{x^{5/2 } \left(1+x^{3-\alpha}\right)^{1/2}}\dfrac{dx\left(t\right)}{dt}.\label{eq : domegadt}\end{aligned}\ ] ] to move from the first line to the second , we have made use of eqs .
( [ eq : orbital_frequency ] ) and ( [ eq : def_x ] ) . the time derivative of @xmath110 displayed in the right - hand side of eq .
( [ eq : domegadt ] ) can be rewritten as a function of @xmath110 by eq .
( [ eq : dxdt_x ] ) .
so we can express @xmath294 as a function of @xmath110 : @xmath296 } c_{\text{gw } } \left [ 1 + \tilde{c } j\left(x\right ) \right ] \times f_{\text{gw}}\left(x\right ) \dfrac{3+\alpha x^{3-\alpha}}{x^{5/2 } \left(1+x^{3-\alpha}\right)^{1/2 } } \nonumber\\   & = \left(gm_{\text{eff } } \right)^{1/2}\varepsilon^{3/\left[2(3-\alpha)\right ] } c_{\text{gw } } \left [ 1 + \tilde{c } j\left(x\right ) \right ] \times \dfrac{3}{4}x^{-11/2 }    \dfrac{\left ( 1 + x^{3-\alpha } \right)^{5/2 } \left(1+\alpha x^{3-\alpha}/3 \right)}{1 + \left(4-\alpha\right)x^{3-\alpha } } \nonumber\\   & = \left(gm_{\text{eff } } \right)^{1/2 } \varepsilon^{3/\left[2(3-\alpha)\right ] } c_{\text{gw } } \left [ 1 + \tilde{c } j\left(x\right ) \right ] \times \dfrac{3}{4}x^{-11/2 } k\left(x\right ) , \label{eq : ddphi_rewriting}\end{aligned}\ ] ] where the function @xmath156 is defined by @xmath297 note that @xmath298 is equal to one when a dm mini - spike is not present around an imbh .
substitution of eqs .
( [ eq : orbital_frequency ] ) and ( [ eq : def_x ] ) into e. ( [ eq : amplitude_t_def ] ) gives @xmath299 after some algebra and simplification .
combining eqs .
( [ eq : a_rewriting ] ) and ( [ eq : ddphi_rewriting ] ) , we finally arrive at the final expression for the amplitude , @xmath300 }   c_{\text{gw}}^{-1/2 }         \left [ 1 + \tilde{c } j\left(x\right ) \right]^{-1/2 } x^{11/4}k\left(x\right)^{-1/2 } \nonumber\\    & = \sqrt{\dfrac{32\pi^5}{3 } } \dfrac{g\mu}{\sourcedistance c^4 } \left ( gm_{\text{eff}}\right)^{-1/4 } c_{\text{gw}}^{-1/2 } \varepsilon^{5/\left[4\left ( 3-\alpha\right)\right ] }            f^{2}x^{19/4 } \left[k\left(x\right ) \left ( 1 + \tilde{c } j\left(x\right ) \right ) \right]^{-1/2 }   \nonumber \\    & = \sqrt{\dfrac{5}{24 } } \dfrac{1}{\pi^{2/3 } }   \dfrac{c}{\sourcedistance } \left ( \dfrac{gm_c}{c^3}\right)^{5/6 } f^{-7/6 }       \chi^{19/4 }   \left[k\left(x\right ) \left ( 1 + \tilde{c } j\left(x\right ) \right ) \right]^{-1/2 } , \label{eq : amplitude_final}\end{aligned}\ ] ] where @xmath146 is defined by @xmath147 and is called the chirp mass . from the second line to the third line , we have used eqs .
( [ eq : def_epsilon ] ) , ( [ eq : c_gw ] ) and ( [ eq : x_delta_epsilon ] ) .
our next task is to express the phase @xmath283 given by eq .
( [ eq : psit_def_appendix ] ) as a function of @xmath301 . from eq .
( [ eq : ddphi_rewriting ] ) , the time derivative of frequency @xmath302 is expressed by @xmath303 } c_{\text{gw } } x^{-11/2 }   \left[k\left ( 1 + \tilde{c } j\right ) \right ] \nonumber \\    & = \dfrac{3}{5}\pi \left ( \dfrac{8\pi gm_c}{c^3}\right)^{5/3 } f^{11/3 } \chi^{-11/2 } \left[k \left ( 1 + \tilde{c } j \right ) \right ] .
\label{eq : dfdt}\end{aligned}\ ] ] we used eqs .
( [ eq : c_gw ] ) and ( [ eq : x_delta_epsilon ] ) to go from the second line to the third line .
( [ eq : dfdt ] ) , we get    @xmath304    where the constant of integration is determined by the initial condition of the gw phase .
collecting the above results , eqs .
( [ eq : amplitude_final ] ) , ( [ eq : phi_f ] ) and ( [ eq:2pift_f ] ) , we finally obtain the gw waveform in the frequency domain :    @xmath305^ { -1/2 } ,    \\    & \mathcal{a}= \left(\dfrac{5}{24 } \right)^{1/2 } \dfrac{1}{\pi^{2/3 } }   \dfrac{c}{\sourcedistance } \left ( \dfrac{gm_c}{c^3}\right)^{5/6 }   \dfrac{1+\cos^2\iota}{2 } , \\    & \psi\left(f \right ) = 2\pi f\left(t_c + \dfrac{\sourcedistance}{c}\right ) - \phi_c -\dfrac{\pi}{4 } - \tilde{\phi}\left(f \right ) , \\    & \tilde{\phi } \left(f\right ) = \dfrac{10}{3 }   \left ( \dfrac{8\pi gm_c}{c^3}\right)^{- 5/3 } \left [     - f \int_\infty^{f } df ' \
\dfrac{\chi^{11/2}}{f'{}^{11/3 } k \left ( 1 + \tilde{c } j \right )   }      + \int_\infty^{f } df ' \
\dfrac{\chi^{11/2}}{f'{}^{8/3 } k \left ( 1 + \tilde{c } j \right ) } \right],\end{aligned}\ ] ]    where @xmath145 is the overall amplitude , @xmath149 is the coalescence time , @xmath151 is the coalescence phase .
note that when a dm mini - spike is not present around an imbh , @xmath161 , @xmath306 , @xmath307 , so the waveform becomes    @xmath308    this is consistent with the waveform from the binary composed of two point - like compact object with mass @xmath64 and @xmath59 @xcite .
the detector sensitivity to the inspiral gws is closely related to the number of gw cycles .
that is because @xmath253 which is defined by eq .
( [ eq : number_of_circles ] ) is proportional to the gw phase which is defined by eq .
( [ eq : phase_t_def ] ) .
therefore snr strongly depends on the number of cycles @xmath253 .
the number of gw cycles in the frequency range @xmath309 $ ] is defined by @xmath310 where an overdot denotes the time derivative and @xmath302 can be calculated by eq .
( [ eq : dfdtau ] ) ( note that @xmath236 . ) .
we define the frequency bandwidth of elisa as the frequency range @xmath311 $ ] within which the square root of the noise spectral density is below half its minimum value : @xmath312 \right),\end{aligned}\ ] ] where @xmath313 is the frequency at which the elisa is most sensitive to the gws .
because we assume 5 year observation , depending on the binary configuration , @xmath314 may be smaller or larger than the initial frequency @xmath315 from which the inspiral gw frequency sweeps to the frequency at the innermost stable circular orbit , @xmath316 .
taking the initial frequency @xmath210 into account , the frequency bandwidth @xmath317 $ ] in which the inspiral gw sweeps is expressed by      and we obtain fig .
[ fig : f_min ] from eq .
( [ eq : number_of_circles ] ) . for @xmath319 ,
the initial frequency @xmath210 at which the gw start to be observed is within the full width at half minimum of @xmath320 .
so the minimum frequency @xmath321 which the inspiral gw spend in the detector bandwidth is equal to the initial frequency @xmath210 . on the other hand , for @xmath262 ,
the initial frequency @xmath210 is out of the full width at half minimum .
so @xmath321 is equal to the lower bound of the detector bandwidth of elisa @xmath314 . since @xmath322 is smaller than @xmath323 for all values of @xmath0 and in the cases we studied , @xmath324    using eqs .
( [ eq : dfdtau ] ) , ( [ eq : number_of_circles ] ) , ( [ eq : fmin ] ) and ( [ eq : fmax ] ) , the number of cycles @xmath253 is obtained in fig . [
fig : ncycle ] .
the figure shows that @xmath253 is almost constant for small @xmath0 but drops sharply for large @xmath0 .
this behavior of @xmath253 is explained by the fact that the dm has more influence on the motion of the stellar mass object for larger @xmath0 . for large @xmath0
, @xmath302 increases sharply as @xmath325 due to the dm effect and the gw frequency of the stellar mass object goes up rapidly through the frequency bandwidth of elisa .
it follows from this that larger @xmath0 leads to wider frequency band but to the less number of gw cycles near the best sensitivity of the detector .
the sensitivity to the gws is determined by the competition between these two effects .


electronic transport in mesoscopic systems has been at the focus of experimental and theoretical interest during the last decade due to recent development in fabrication of small electronic devices and their interesting equilibrium and non - equilibrium properties .
especially interesting are the transport properties of a quantum dot ( qd ) under the influence of external time - dependent fields .
the high - frequency signals may be applied to a qd and the time - dependent fields will modify the tunneling current .
new effects have been observed and theoretically described , e.g. photon - assisted tunneling through small quantum dots with well - resolved discrete energy states [ 1 , 2 , 3 ] , photon - electron pumps [ 4 , 5 , 6 ] and others .
one can investigate the current flowing through a qd under periodic modulation of the qd electronic structure [ 7 ] or periodic ( non - periodic ) modulation of the tunneling barriers [ 6 ] and electron energy levels in both ( left and right ) electron reservoirs [ 8 ] ( see also [ 9 , 10 ] ) .
the progress of nanomaterials science has enabled the experimental study of the phase coherence of the charge carriers in many mesoscopic systems .
the asymmetric fano line shapes [ 11 ] are observed whenever resonant and nonresonant scattering paths interfere . in some nanostructures ,
e.g. in single - electron transistors , the fano resonances in the conductance were observed [ 12 ] , which imply that there are two paths for transfer of electrons between a source and a drain .
especially the recent experimental and theoretical study with a low - temperature scanning tunneling microscope ( stm ) of the single magnetic atom deposited on a metallic surface showed the asymmetric fano resonances in the tunneling spectra [ 1215 ] .
the stm measurements indicate that in tunneling of electrons between stm tip and a surface with a single impurity atom two different paths are present .
the electrons can tunnel between the tip and the adsorbate state and directly between the tip and the metal surface .
the electronic transport through a qd coupled to the electron reservoirs within a model with two electron tunneling channels was considered in ref .
[ 16 ] and it was shown , that transport of electrons through both channels leads to an asymmetric shape of the zero bias voltage conductance curves , which is typical behaviour for a fano resonance resulting from constructive and destructive interference processes for electrons transmitted through both channels .    in all papers
mentioned above and relating to the electron transport through a qd with the additional ( bridge ) transmission channel the external fields were not applied and the considered systems were driven out of equilibrium only by means of a dc voltage bias . in this paper we address the issue of a qd with a bridge channel between a source and a drain driven out of equilibrium by means of a dc voltage bias and additional time - dependent external fields . in this manner ,
our paper can be seen as generalization of ref . [ 9 ] to the case of a qd with the additional bridge channel in the presence of external microwave fields which are applied to the dot and two leads , respectively . in literature
, different theoretical approaches have been developed to treat the time - dependent , nonequilibrium electron transport processes in the mesoscopic systems .
it seems , that the most popular is the nonequilibrium green s function method .
however , these green s functions depend on the time arguments and for non trivial quantum models it is a rather difficult task to calculate them . in our treatment of the time - dependent tunneling through a mesoscopic system we use the evolution operator technique ( e.g. [ 17 , 18 ] )
. the final expressions for the qd charge and the current flowing in the system can be described in terms of the corresponding matrix elements of the evolution operator .
in our earlier work [ 19 ] we have considered the similar problem solving numerically the corresponding sets of the differential equations satisfied by the matrix elements of the evolution operator . due to the complexity and a large number of these equations we considered only a very limited number of intereting cases although we were able to take into consideration the electronic structure of the lead energy bands and the specific time - dependence of the qd - lead barriers . here
we give the analytical expressions for the qd charge and current assuming so called the wide - band limit approximation and the time - independent strength of the qd - lead barriers . as a result , due to these analytical forms , we are able to analyze the required characteristics of the considered system for the very broad class of parameters . additionally , due to the final forms given for some matrix elements of the evolution operator it is possible to build up the expressions for the qd charge or current in the form of the perturbation series .    in the next section we present the model and formalism and give the resulting expressions for equations for the corresponding matrix elements of the evolution operator . in section 3
we obtain the approximate solutions for all required matrix elements and give the final forms for the qd charge and the current flowing in the system .
the results are presented in section 4 which includes also the summary and a brief discussion .
we model the qd coupled to the left and right electron reservoirs with the additional bridge tunneling channel between them by the usually used hamiltonian @xmath0 , where & & h_1 = _ k _ _ k_(t ) a^+_k _ a_k _ + _ d(t ) a^+_d a_d , + & & v = _ k _ v_k_d(t ) a^+_k_a_d + h.c . + _
k_l , k_r v_k_lk_r(t ) a^+_k_l a_k_r + h.c .
the operators @xmath1 , @xmath2 are the annihilation ( creation ) operators of the electron in the lead @xmath3 ( @xmath4 ) and in the qd , respectively .
the couplings between qd and lead states and between both lead states are denoted by @xmath5 and @xmath6 , respectively . for simplicity
, the dot is characterized only by a single level @xmath7 and we have neglected the intradot electron - electron coulomb interaction .
we assume the case in which there exist microwave fields applied to the leads and qd .
in the adiabatic approximation our time - dependent driven system is described by @xmath8 , @xmath9 , i.e. the energy levels of the leads and qd are driven by the ac field with the frequency @xmath10 and the amplitudes @xmath11 and @xmath12 , respectively .    we describe the dynamical evolution of the charge localised on the qd and the current flowing in the system in terms of the time evolution operator @xmath13 ( in the interaction representation ) which satisfies the equation it u(t , t_0 ) = v(t)u(t , t_0 ) , where & & v(t ) = u_0(t , t_0 ) v(t ) u^+_0(t , t_0 ) , + & & u_0(t ) = t(i^t_t_0 dt h_1(t ) ) . here
we assume that the interaction between qd and leads and between both leads is switched on in the distant past @xmath14 , i.e. @xmath15 and @xmath16 equal to zero for @xmath17 and takes constant values for @xmath18 .
the qd charge and currents flowing in the system can be obtained from the knowledge of the appropriate matrix elements of the evolution operator @xmath13 .
the qd charge is given as follows ( cf .
[ 18 ] ) : n_d(t ) = n_d(t_0)|u_dd(t , t_0)|^2 + _ k _ n_k _ ( t_0)|u_dk _ ( t , t_0)|^2 , where @xmath19 and @xmath20 denote the matrix elements of @xmath13 calculated within the basis functions containing the electron single - particle functions of the leads and qd , @xmath21 , @xmath22 and @xmath23 , respectively . @xmath24 and @xmath25 represent the initial filling of the corresponding single - particle states
.    the tunneling current flowing , e.g. from the left lead into the qd and the right lead , @xmath26 , can be obtained from the time derivative of the total number of electrons in the left lead , @xmath27 ( cf . [ 9 ] ) , where n_l(t ) & = & _ k_l n_k_l(t ) = _ k_l(n_d(t_0)|u_k_l d(t , t_0)|^2 + & + & _
q_l n_q_l(t_0)|u_k_l , q_l(t , t_0)|^2 + _
k_r n_k_r(t_0)|u_k_l , k_r(t , t_0)|^2 ) .
let us begin with the calculations of the qd charge @xmath28 .
then we have to calculate the matrix elements @xmath29 and @xmath30 . using the identity operator @xmath31 the following set of coupled equations
can be obtained from eq .
( 3 ) : & & t u_dd(t , t_0 ) = -i_k_d v^*_k_d(t ) u_k_d(t , t_0 ) , + & & t u_k_l d(t , t_0 ) = -iv_k_l d(t ) u_dd(t , t_0 ) - i_k_r v_k_lk_r(t ) u_k_r d(t , t_0 ) , + & & t u_k_r d(t , t_0 ) = -iv_k_r d(t ) u_dd(t , t_0 ) - i_k_l v_k_rk_l(t ) u_k_l d(t , t_0 ) , where v_k_d(t ) & & k_|v(t)|d= v_k_d ( i(_- _ d)(t - t_0 ) + & & - i(_d - _ ) ( t - t_0)/ ) , + v_k_lk_r(t ) & & k_l|v(t)|k_r= v_k_lk_r ( i(_k_l - _ k_r)(t - t_0 ) - + & & - i(_l - _ r)(t - t_0)/ ) .
it is easy to show that the equation for @xmath29 can be written as follows : = -^t_t_0dt(k(t , t)u_dd(t,t_0 ) + _ k _ l_k_(t ,
t ) u_k_d(t,t_0 ) ) , where ( t , t ) = _ k _ v^*_k_d(t ) v_k_d(t ) , + l_k_l(t , t ) = _
k_r v^*_k_r(t ) v_k_rk_l(t ) , and the similar equation can be written for @xmath32 .    the formal solution of eq . ( 9 ) written in the form u_k_l d(t , t_0 ) = -i^t_t_0 dt v_k_l d(t ) u_dd(t,t_0 ) - i_k_r ^t_t_0 dt v_k_lk_r(t ) u_k_r d(t ) , ( and the similar equation for @xmath33 ) can be iterated giving     + & + & ^_j=2 ( -i)^j _ r_1,l_2,r_3,l_4,
,_j-1 ^t_t_0 dt_1
 ^t_j-1_t_0 dt_j v_lr_1(t_1 ) v_r_1 l_2(t_2 )
v__j-1 ( t_j)u_dd(t_j , t_0 ) .
+ & & here we have introduced the abbreviated form for the vector @xmath34 and replaced it by @xmath3 and @xmath35 . the equation for @xmath36
can be obtained from eq .
( 17 ) by interchanging @xmath37 . inserting these expressions for @xmath38 and @xmath39 into eq . ( 13 ) one obtains the closed integro - differential equation for @xmath29 + & - & ^_j=2 ( -i)^j-1 _ l_1,r_2,l_3,  ,_j ^t_t_0 dt_1  ^t_j-1_t_0 dt_j v^*_l_1(t ) v_l_1 r_2(t_1 )
v__j ( t_j ) u_dd(t_j , t_0 ) + & & + ( the  second  term  with  the  change  lr ) . equation ( 18 ) together with the expressions for @xmath40 and @xmath41 , eqs .
( 11 , 12 ) , and for @xmath42 written as follows ( t , t_1 ) = _
= l , r |v_|^2 d_(t - t_1 ) ( i_d(t - t_1 ) + i(_d - _ ) ( t - t_1)/ ) where @xmath43 a fourier transform of the @xmath3-th lead density of states , gives exact , closed equation for @xmath29 .
here we have assumed that @xmath44 does not depend on the wave vector @xmath34 and then the similar assumption will be made for @xmath45 .    under the wide - band limit ( wbl ) approximation ( e.g. [ 9 ] ) this equation
can be analytically solved and such solutions will be considered later .
formally , solving eq . ( 18 ) and inserting its solution to eq .
( 17 ) , the solutions for @xmath38 and @xmath39 can be obtained .
these functions are needed in the calculations of the first term of @xmath26 ( see eq .
7 ) .    in order to calculate @xmath28 we still need @xmath46 . writing down the set of closed equations for @xmath47 , @xmath48 and @xmath49 ( obtained on the basis of eq .
( 3 ) ) and performing similar calculations to those described above , one obtains for @xmath47 ( and similar equation for @xmath50 by interchanging @xmath37 ) : + & - & iv^*_r(t ) + ( -i)^2 ^t_t_0dt_1 _ l v^*_l(t ) v_lr(t_1 ) + & + & ^_j=2,4,6 ,  ( -i)^j+1 _ r_1,l_2,r_3,  ,l_j ^t_t_0dt_1  ^t_j-1_t_0 dt_j v^*_r_1(t ) v_r_1l_2(t_1 )
v_l_jr(t_j ) + & + & ^_j=3,5,7 ,
( -i)^j+1 _ l_1,r_2,l_3,
 ,l_j ^t_t_0dt_1  ^t_j-1_t_0 dt_j v^*_l_1(t ) v_l_1r_2(t_1 )
v_l_jr(t_j ) + & + & ( -i)^2 ^t_t_0 dt_1 _ l v^*_l(t ) v_l(t_1 ) u_dr(t_1,t_0 ) + & + & ^_j=2 ( -i)^j+1 _ l_1,r_2,l_3,
,_j ^t_t_0dt_1  ^t_j-1_t_0 dt_j v^*_l_1(t ) v_l_1r_2(t_1 )
v__j(t_j ) u_dr(t_j , t_0 ) , and for @xmath49 ( needed for calculation of the last term of @xmath51 , eq .
( 7 ) ) : + & + & ^_j=3,5,7  ( -i)^j _ r_1,l_2,r_3,  ,l_j-1 ^t_t_0dt_1^t_1_t_0 dt_2  ^t_j-1_t_0 dt_j v_lr(t_1)v_r_1l_2(t_2 )
v_l_j-1r(t_j ) + & - & i^t_t_0 dt_1 v_l(t_1 ) u_dr(t_1,t_0 ) + & + & ^_j=2 ( -i)^j _ r_1,l_2,  ,_j-1 ^t_t_0dt_1 ^t_1_t_0dt_2  ^t_j-1_t_0 dt_j v_lr_1(t_1 ) v_r_1l_2(t_2 )
v__j(t_j ) u_dr ( t_j , t_0 ) .
the analytical solutions of these equations under the wbl approximation will be discussed in the next section .    for calculation of @xmath52 one
needs the functions @xmath53 , @xmath49 and @xmath54 .
the first two functions are given in eqs .
( 17 ) and ( 21 ) and @xmath54 should be calculated from the set of coupled equations for @xmath50 , @xmath54 and @xmath55 obtained from eq .
the result is as follows : + & & - i^t_t_0 dt_1 v^*_l_1(t_1 ) u_dl(t_1,t_0 ) + & + & ^_j=2 ( -i)^j _ r_1,l_2,r_3,  ,_j-1 ^t_t_0 dt_1 ^t_1_t_0 dt_2  ^t_j-1_t_0 dt_j
v_l_1r_1(t_1 ) v_r_1l_2(t_2 )  v__j-1l(t_j)u_dl(t_j , t_0 ) , where @xmath50 is given by solving eq .
( 20 ) with the replacement @xmath37 .
in order to calculate the qd charge @xmath28 or current @xmath26 one has to solve , in the first step , the integro - differential equations satisfied by @xmath29 , eq . ( 18 ) , and @xmath56 , eq .
the other needed functions @xmath53 , @xmath54 and @xmath49 can be obtained from eqs .
( 17 , 21 , 22 ) inserting into them @xmath29 , @xmath57 and performing multiple time and @xmath58-vector integrations .
unfortunately , it is a very difficult task to solve the integro - differential equations and to perform these integrals in a general case when the leads are characterized by some density of state curves . here , we use the wbl approximation , under which all multiple time and @xmath58-vectors integrals can be performed without difficulties .
the wbl approximation has been widely used in calculating of many properties of mesoscopic systems ( e.g. [ 2 , 3 , 5 , 9 , 10 ] ) .
it is justified under the conditions that the qd energy level linewidth is much smaller than the bandwidth of the leads and that the density of states and hopping matrix elements vary slowly with energy .
furthermore , as we are not going to consider the case of the qd energy level lying close to the edges of the leads energy band or lying close to some singular structure present in the leads density of states , then application of the wbl approximation should be fully justified .
the conditions under which we perform our calculations are satisfied in most experimental constructions of mezoscopic systems . as a check
, we have performed the direct but time consuming numerical integration of eqs .
( 810 ) ( and similar equations for other functions ) for the rectangular leads density of states and did not find any differences in the results for the time - averaged qd charge or the currents flowing in the considered system .
therefore , in this section we consider the electron transport through the qd with the additional bridge tunnel over the dot within the wbl approximation . in this approach
the solutions of the integro - differential eqs .
( 19 ) and ( 20 ) in the analytical form can be obtained and the infinite sums of all terms in eqs .
( 17 ) and ( 21 ) can be calculated .
let us consider the equation for the matrix element @xmath29 , eq .
the function @xmath59 , eq .
( 19 ) , is approximated in wbl as follows + & & + i(_d - _ ) ( t - t_1)/ ) + & & _ 2(t - t_1 ) , where the leads density of states @xmath60 was replaced by the rectangular density of states with the  effective  bandwidth @xmath61    using similar approximation in calculations of the multiple integrals present in eq . ( 18 ) one obtains u_dd(t , t_0 ) = ( -2 - 2v^2d ^_j=1 ( -iv_rl / d)^j ) u_dd(t , t_0 ) , where @xmath62 , @xmath63 , @xmath64 , @xmath65 .    assuming @xmath66 the series can be summed up and finally the equation for @xmath29 reads u_dd(t , t_0 ) = -c_1 u_dd(t , t_0 ) , where @xmath67
.    similarly , eq .
( 20 ) becomes u_d(t , t_0 ) = -iv_dd(t)/(1 + iv_rl /d ) - c_1 u_d(t , t_0 ) .
the solutions of eqs .
( 25 ) and ( 26 ) are as follows : & & u_dd(t , t_0 ) = ( -c_1(t - t_0 ) ) , + & & u_d(t , t_0 ) = -i/(1 + iv_rl / d ) ^t_t_0
dt_1 v_d(t_1)(-c_1(t - t_1 ) ) and the qd charge @xmath28 can be easily obtained , eq .
it can be verified , that the first term of eq .
( 6 ) tends to zero as @xmath68 and finally for the qd charge we have n_d(t ) = 1 ( 1 + ( v_rl / d)^2 ) _
d f_()|a_(,t)|^2 , where a_(,t ) = -^t_t_0dt_1 ( i(_d - ) ( t - t_1 ) - i(_d - _ ) ( t - t_1)/- c_1 ( t - t_1 ) )    in the limit of vanishing bridge over the qd , eq .
( 29 ) reproduces the result of ref . [ 9 ] .
the current @xmath26 flowing from the left lead into the qd and the right lead is calculated from the evolution of the total number of electrons in the left lead ( see eq .
7 ) and one can read : j_l(t ) & = & 2 re\{n_d(t_0 ) _ k_l u^*_k_ld(t , t_0 ) ddt u_k_ld(t , t_0 ) .
+ & + & _ k_l
q_l n_q_l(t_0 ) u^*_k_l , q_l ( t , t_0 ) ddt u_k_l , q_l(t , t_0 ) + & + & . _
k_l , k_r n_k_r(t_0 ) u^*_k_l , k_r(t , t_0 ) ddt u_k_l , k_r(t , t_0 ) } .
the functions @xmath69 , @xmath70 and @xmath71 are calculated according to eqs .
( 17 ) , ( 22 ) and ( 21 ) , respectively , and after summing up of the corresponding multiple integrals one obtains : u_k_ld ( t , t_0 ) & = & -i1 + ix ^t_t_0dt_1 v_k_ld(t_1 ) u_dd(t_1,t_0 ) , + u_k_l , q_l(t , t_0 ) & = & _ k_l , q_l - v_rl x 1 + x^2 ^t_t_0 dt_1 \{i(_k_l - _ q_l ) ( t_1-t_0 ) } + & - & i1 + ix ^t_t_0 dt_1 v_k_ld(t_1 ) u_dq_l(t_1,t_0 ) , + u_k_l , k_r(t , t_0 ) & = & -i1 + ix^t_t_0 dt_1v_k_ld(t_1 ) u_dk_r(t_1,t_0 ) + & - & i1 + x^2 ^t_t_0 dt_1 v_k_lk_r(t_1 ) , where @xmath29 , @xmath72 and @xmath73 are given in eqs .
( 27 , 28 ) and @xmath74 .
one can verify , that the first term in eq .
( 31 ) tends to zero for @xmath68 as we have for this term n_d(t_0 ) |1/(1 + ix)|^2 /2
( -(t - t_0 ) re c_2 ) , where @xmath75 .
the second and third terms of the expression for @xmath26 , eq . ( 30 ) , together with eqs .
( 27 , 28 , 31 - 33 ) give finally the time averaged current : + & & + i m \{1 - x^2(1 + x^2)(1 + ix)^2 _
l df_l ( ) a_l(,t ) } + & & + re \{2x(1+x^2)(1 + ix ) _
l df_r ( ) a_r(,t)}. in the vanishing bridge channel case eq .
( 35 ) coincides with the results of ref . [ 9 ] : j^v_rl=0_l(t)= 2 n_d(t)+ _
l df_l ( ) i m a_l ( , t ) .
the current @xmath76 , eq .
( 36 ) , flowing from the left lead to the central region and to the right lead ( through the bridge channel ) is the superposition of four terms .
the first term corresponds to the current between two leads and this term is not disturbed by the qd .
the form of the second and third terms is the same as for @xmath77 , eq .
( 37 ) , except for the renormalization constants due to the additional tunneling channel .
note , that some additional renormalizations also occur due to @xmath78 which enters into the expression for @xmath28 , eq .
( 29 ) , and for @xmath79 , eq .
the last term of eq .
( 36 ) is the interference term due to the simultaneous tunneling through two channels .
here we show the numerical results of the time - averaged qd charge @xmath80 and the current @xmath76 and its derivatives with respect to the qd energy level position and the chemical potential @xmath81 ( or equivalently , with respect to the gate and source - drain voltages ) for different sets of parameters which characterize our system .
we assume the temperature @xmath82  k and the time - average of time - dependent quantities @xmath83 is defined by f(t)= _ 2 12 ^_- dt f(t ) , and because @xmath83 is a periodic function of time , we average it over the period @xmath84 .
we take the chemical potential of the right lead @xmath85 as the energy measurement point , @xmath86 . as the potential drop between
the left and right leads is given by @xmath87 and @xmath88 is the measured voltage between a source and a drain , then the derivatives of the current @xmath76 with respect to @xmath81 will correspond to the derivatives @xmath89 usually measured in experiments . in experiments the gate voltage controls the position of the qd energy level @xmath7 ( regardless
how complicated the relation between the gate voltage and @xmath7 is ) and for that reason to mimic measurements of the qd charge or current vs. the gate voltage we have calculated them vs. the position of the qd energy level @xmath7 .
the values of the hybridization matrix elements @xmath5 present in the hamiltonian do not enter the final expressions for the current or qd charge obtained within the wbl approximation .
usually the effective linewidth @xmath90 , responsible for the additional tunneling channel for which we should take some values in order to perform numerical calculations .
we have taken the values comparable with @xmath5 and estimated @xmath5 ( assuming its @xmath58-independence , @xmath91 ) using the relation @xmath92 , where @xmath61 is the @xmath3-lead s bandwidth and @xmath93 ( @xmath94 , @xmath64 was assumed ) .
we assumed the amplitude of the qd energy levels oscillation @xmath12 to be one half of @xmath95 and @xmath96 , if otherwise stated . in our calculations
the values @xmath78 were taken from the range ( 010 ) , in @xmath97 units .    in the first three figures , figs .
13 , we present the overall shape of the average current @xmath76 and the derivatves of the average current with respect to the qd energy level @xmath7 and the left chemical potential @xmath81 , respectively , against @xmath7 and @xmath81 .
the upper panels correspond to the @xmath98 case and the lower ones present the results obtained for the non - vanishing over - dot tunneling channel , @xmath99 .
there are quite visible differences between the case of a qd with and without additional bridge tunneling channel ( cf .
upper and lower panels in figs 13 ) . for better visualization of the peculiarities of the presented functions and for simpler discussion
let us consider the specific cuts of the surfaces given in figs . 1 - 3 .    at first ,
let us consider the dependence of @xmath76 vs. @xmath7 for given values of the left lead chemical potential @xmath81 . in fig .
4 we present such curves for different values of @xmath81  the subsequent curves beginning from the lower one correspond to @xmath100 up to the upper curve ( with the step @xmath101 ) obtained for @xmath102 .
the left ( right ) panel corresponds to @xmath98 ( @xmath103 ) . in the case of vanishing over - dot tunneling channel ( the left panel )
the current has a simple structure  a single peak localized in the middle between @xmath81 and @xmath85 for smaller values of @xmath81 .
the width of this peak increases with increasing @xmath104 ( @xmath86 ) and for greater values of @xmath104 the current is almost independent of @xmath7 localized inside the energy region between @xmath85 and @xmath81 . for non - vanishing over - dot tunneling ( fig .
4 , the right panel ) the curves @xmath76 become asymmetric . with increasing source - drain bias ,
the current possesses greater values in comparison with the @xmath98 case due to the direct tunneling between both leads .
note , however , that due to the interference effects the resulting @xmath76 curves are asymmetric .
the interference effects are most visible approximately for @xmath7 lying in the region ( @xmath105 ) .    in fig .
5 we show the average current @xmath76 vs. the left lead chemical potential @xmath81 for several values of @xmath7 . for vanishing @xmath78 the corresponding curves are nearly asymmetric with the asymmetry point @xmath106 . with increasing @xmath81 at fixed @xmath7 the current achieves a constant value depending on the position of the qd energy level @xmath107 with respect to the @xmath86 .
it means , that electrons which occupy the lead energy levels not too distant from @xmath7 take part in the tunneling process . for greater @xmath81
most of the lead energy levels lying far away from the @xmath7 are inactive in the tunneling between leads through the qd energy level .
however , at non - vanishing @xmath78 ( see the right panel of fig .
5 ) these lead energy levels are active and the current @xmath76 vs. @xmath81 is of much richer structure .
the current is nearly linearly growing with the increasing @xmath81 ( for larger @xmath81 ) because the tunneling through the qd can be neglected compared with the direct tunneling between both leads .
the clearly visible interference effects appear only for @xmath81 not too distant from @xmath7 .    in fig .
6 we show the derivatives of the average current vs. the qd energy level @xmath7 obtained for some values of @xmath81 .
these are the results of the intersection of the surfaces given in fig .
2 with the planes at constant values of @xmath81 or the results of the differentiation of curves shown in fig .
, the most visible differences between the results obtained for @xmath98 and @xmath99 are present for the qd energy level @xmath7 localized approximately between chemical potentials of both leads ( compare , for example , the curves calculated for @xmath102 ) .
7 presents the comparison of the @xmath108 vs. @xmath81 curves calculated for vanishing @xmath78 ( left panels ) and for @xmath103 ( right panels ) for two different values of the amplitudes @xmath95 ( @xmath109 , @xmath96 ) . at the vanishing value of @xmath78 , the shape of the curves is symmetrical in relation to the values @xmath106 although for greater @xmath95 some shoulders appear on both sides of the corresponding peaks in the distance @xmath110 from the curve centres . for non - vanishing @xmath78 , the corresponding curves are approximately asymmetric and for large @xmath81 they tend to constant , non - zero values corresponding to linear increasing of the current at large @xmath81 .
it is interesting that with the increasing amplitudes @xmath95 and @xmath12 very clear structures appear on both sides of the corresponding curves .
note , that all these curves can be obtained , for example , from the one calculated for @xmath111 and moved along the @xmath81-axis by the corresponding value ( equal to @xmath7 ) .    fig .
8 shows the dependence of the average current @xmath76 on the qd energy level @xmath7 and the direct coupling @xmath78 between both leads .
the upper and bottom panels correspond to different values of the amplitude @xmath95 ( @xmath109 ) .
a very distinct transition from the symmetric to nearly antisymmetric behaviour of the current vs. @xmath7 curves is observed with the increasing value of the over - dot additional coupling between the leads for both values of @xmath95 .
the larger amplitudes of the left lead and qd levels oscillations result only in some broadening of the characteristic features of the average current vs. @xmath78 and @xmath7 surface and do not introduce any additional structures on these surfaces ( for the range of parameters where @xmath112 and @xmath113 is not very small ) .    fig .
9 presents the average current @xmath76 obtained for the case in which only the qd energy level @xmath7 oscillates with some frequency ( @xmath114 ) and for the small source - drain voltage @xmath115 .
the subsequent panels a , b , c and d correspond to the increasing value of the amplitude @xmath12 and the different curves ( broken , thin , thick and very thick ) describe the qd without the over dot channel , @xmath98 , and with this channel at @xmath116 and 10 , respectively . for vanishing @xmath78 ( broken curves )
we observe for small amplitude @xmath12 only a central resonant peak ( fig .
9a ) . with increasing @xmath12 ,
the subsequent peaks appear and the distance between them and the central peak is an integer multiple of the frequency @xmath10 ( sidebands ) .
the location of peaks is independent of the amplitude @xmath12 but their relative intensity values change and with increasing @xmath12 the height of the central peak is reduced . for the relatively large amplitude @xmath12 , the heights of the two neighbouring peaks ( at @xmath117 )
are approximately equal to the height of the central peak .
if we take the additional over - dot tunneling channel into consideration , especially for small @xmath12 , the asymmetric shape of the current curve is observed and this asymmetry increases with the increasing strength of the over - dot coupling between both leads . with the increasing amplitude @xmath12 this asymmetry
is reduced largely due to the extra , photon - assisted tunneling peaks whose strength increases with the increasing @xmath12 . for sufficiently large values of @xmath12 and @xmath78 the functional dependence of the average current on the qd energy level @xmath7 is nearly the same for vanishing and non - vanishing over - dot tunneling channels .
there is only one difference  for large @xmath78 the corresponding curve is shifted to the higher values due to the direct channel between both leads .
the last two figures 10 , 11 are devoted to the analysis of the average current @xmath76 dependence on the oscillation period @xmath84 of the external fields . in fig .
10 we show the overall dependence of @xmath76 on @xmath84 and the qd energy level @xmath7 .
the upper ( lower ) part of the figure presents the results for the vanishing ( non - vanishing ) over dot channel between both leads .
the most visible differences between averaged currents calculated for @xmath118 and @xmath99 can be observed for @xmath119 , especially for small values of @xmath84 .
more detailed analysis of the @xmath76 dependence on the oscillation period @xmath84 is shown for some chosen values of the parameters @xmath7 and @xmath78 in fig .
the thin ( thick ) curves correspond to @xmath120 ( @xmath121 ) and the solid ( broken ) curves correspond to @xmath122 ( @xmath98 ) .
additionally , we give the resutls for two values of the amplitude @xmath95 ( @xmath109 , @xmath96 ) , i.e. for @xmath123 and @xmath124 , the left and the right parts of fig .
11 respectively .
we observe the characteristic average current oscillations damped with increasing @xmath84 for @xmath7 lying in the central part between the left and right chemical potentials ( see also [ 9 ] ) .
these oscillations are present for both @xmath98 and @xmath99 and more visible for greater amplitudes @xmath95 and @xmath12 but the maxima and minima of the oscillating average current are localized at the same value of the oscillation period .
note , that the existence of the additional over - dot tunneling channel results approximately in shifting the corresponding curves to higher values without any additional serious modifications . for the qd energy level lying away from the middle point between the left and right chemical potentials ,
the average current is still an oscillating function of @xmath84 although these oscillations are less transparent and their oscillation period is much longer .    to conclude , we have provided a detailed investigation of a qd connected to two leads with an additional over - dot tunneling channel .
the harmonic external microwave fields were considered as applied to the qd and two leads which result in time - dependence of the corresponding qd and leads energy levels . the qd charge and
the average current flowing in this system were calculated within the evolution operator technique .
the corresponding matrix elements of the evolution operator required for the calculation of the qd charge and current were presented in the form of the infinite series of multiple time integrals of the functions containing the information about the coupling between the qd and leads or in the form of the integro - differential equation . applying the wbl aproximation we were able to obtain all required evolution operator matrix elements in closed forms and
give the final analytical expressions for the qd charge and current .
we have performed the extended numerical calculations for the qd charge , the average current and the derivatives of the current with respect to the gate and source - drain voltages ( fig .
the most spectacular influence of the additional bridge tunneling channel is visible in the @xmath76 dependence vs. the position of the qd energy level at the constant source - drain voltage .
going from the vanishing values of the over - dot tunneling channel strength to the non - vanishing one due to the interference effects , the corresponding curve transforms from the symmetric to nearly antisymmetric shape
. similar influence of the non - vanishing @xmath78 is visible in the dependence of @xmath108 vs. @xmath81 ( fig .
the characteristic behaviour of the average current vs. the qd energy level position at the small source - drain voltage is observed for the case when the external oscillating field is applied only to the qd ( fig .
9 ) . for the vanishing over - dot tunneling channel at the small amplitude @xmath12 ,
the main resonant peak is only observed and with the increasing amplitude @xmath12 the next peaks localized at @xmath7 equal to the multiplicity of @xmath10 appear corresponding to the photon - assisted tunneling . in that case
the current is a symmetric curve centered around the main resonant peak . for the nonvanishing over - dot tunneling channel and small amplitudes @xmath12 , the character of the dependence of the average current on the qd energy level position transforms with increasing @xmath78 from the symmetrical to nearly antisymmetrical behaviour .
this tendency to the antisymmetrical behaviour with increasing @xmath78 at small amplitude @xmath12 is reduced with increasing @xmath12 . for sufficiently large amplitude
@xmath12 the overall behaviour of the average current vs. @xmath7 is very similar for different values of @xmath78 and does not manifest the tendency for the antisymmetry with increasing @xmath78 .
1 .   t. h. oosterkamp , l.p .
kouwenhoven , a.e.a .
koden , n.c .
van der vaart , c.j.p.m .
harmans , phys .
lett . * 78 * , 1536 ( 1997 ) 2 .
qing - feng sun , tsung - han lin , phys .
* b56 * 3591 ( 1997 ) 3 .
qing - feng sun , jian wang , tsung - han lin , phys . rev . * b58 * , 13007 ( 1998 ) 4 .   c.a .
stafford , n.s .
wingreen , phys .
lett . * 76 * , 1916 ( 1996 ) 5 .
qing - feng sun , tsung - han lin , j. phys . : condens .
matter * 9 * , 3043 ( 1997 ) 6 .
kouwenhoven , a.t .
johnson , n.c .
van der vaart , a. van der enden , c.j.p.m .
harmans , c.t .
foxon , z. phys .
b , condens . matter * 85 * , 381 ( 1991 ) 7 .
zhao , j. wang , eur .
j. * b9 * , 513 ( 1999 ) 8 .
y. goldin , y. avishai , phys .
rev . * b61 * , 16750 ( 2000 ) 9 .
jauho , n.s .
wingreen , y. meir , phys . rev . *
b50 * , 5528 ( 1994 ) 10 .
qing - feng sun , tsung - han lin , j. phys . : condens .
matter * 9 * , 4875 ( 1998 ) 11 .
u. fano , phys . rev . *
124 * , 1866 ( 1969 ) 12 .
j. gres , d. goldhaber - gordon , s. heemeyer , m.a .
kastner , h. shtrikman , d. mahalu , u. meirav , phys . rev . *
b62 * , 2188 ( 2000 ) 13 .
m. plihal , j.w .
gadzuk , phys . rev . *
b63 * , 085404 ( 2001 ) 14 . v. madhavan , w. chen , t. jamneala , m.f .
crommie , n.s .
wingreen , science * 280 * , 567 ( 1998 ) 15 .
o. ujsaghy , j. kroha , l. szunyogh , a. zawadowski , phys .
lett . * 85 * , 2557 ( 2000 ) 16 .
buka , p. stefaski , phys .
* 86 * , 5128 ( 2001 ) 17 .
m. tsukada , n. shima , in : `` dynamical processes and ordering on solid surfaces '' , eds .
a. yoshimori and m. tsukada ( springer , berlin 1995 ) , p. 34
grimley , v.c.j .
bhasu , k.l .
sebastian , surf .
* 121*,305 ( 1983 ) 19 .
t. kwapiski , r. taranko , physica e , in press     against the left lead chemical potential @xmath81 and qd energy level @xmath7 for @xmath98 ( the qd without the over - dot channel , the upper panel ) and for @xmath103 ( lower panel ) .
@xmath86 , @xmath125 , @xmath126 , @xmath127 , @xmath96 , @xmath128 and all energies are given in @xmath97 units.,scaledwidth=90.0% ]         against @xmath7 for given values of @xmath81 ( beginning from @xmath100 up to @xmath102 ) .
the left and right panels correspond to @xmath98 and @xmath103 , respectively , and the other parameters as in fig .
1.,scaledwidth=100.0% ]     against @xmath81 for given values @xmath7 ( beginning from @xmath130 up to @xmath131 ) .
the broken curves correspond to @xmath111 .
the left ( right ) panel corresponds to @xmath98 ( @xmath132 .
the other parameters as in fig .
1.,scaledwidth=100.0% ]     with respect to the qd energy level @xmath7 for given values of @xmath81 ( beginning from @xmath100 up to @xmath102 ) .
the broken curves correspond to @xmath100 .
the left ( right ) panel corresponds to @xmath98 ( @xmath103 ) and the other parameters as in fig .
1.,scaledwidth=100.0% ]    , @xmath108 , for given values of @xmath7 ( beginning from @xmath130 up to @xmath131 ) .
the left ( right ) panels correspond to @xmath98 ( @xmath103 ) and upper ( lower ) panels correspond to @xmath126 ( @xmath133 ) .
the other parameters as in fig .
1.,scaledwidth=90.0% ]       against @xmath7 for the oscillating qd energy level at @xmath78 = 0 , 4 , 7 , 10  broken , thin , thick and very thick curves , respectively .
the panels a , b , c and d correspond to @xmath136 and 7 , respectively .
@xmath137 , @xmath138 , @xmath125 , @xmath139 , @xmath140.,scaledwidth=90.0% ]     against @xmath7 and the period of the time oscillation of the external field for @xmath98  the upper panel and for @xmath103  the lower panel .
@xmath125 , @xmath141 , @xmath124 , @xmath142 , @xmath96.,scaledwidth=90.0% ]     against the period of the time oscillation of the external field for @xmath98  the broken curves , and for @xmath122  the solid curves .
the thin ( thick ) curves correspond to @xmath121 ( @xmath120 ) .
the left ( right ) panel corresponds to @xmath123 , @xmath143 ( @xmath124 , @xmath142 ) .
@xmath96 , @xmath125 , @xmath141.,scaledwidth=100.0% ]