compact binary ( neutron star - neutron star , ns@xmath1 , or black hole - neutron star , bh - ns ) mergers are prime sources of gravitational radiation .
the gw detectors ligo @xcite , virgo @xcite and geo600 @xcite are designed to optimally detect merger signals .
these detectors have been operational intermittently during the last few years reaching their nominal design sensitivity @xcite with detection horizons of a few dozen mpc for ns@xmath1 and almost a hundred mpc for bh - ns mergers ( the ligo - virgo collaboration adopts an optimal canonical distance of 33/70mpc ; @xcite ) .
both ligo and virgo are being upgraded now and by the end of 2015 are expected to be operational at sensitivities @xmath2 times greater than the initial ligo @xcite , reaching a few hundred mpc detection horizon for ns@xmath1 mergers and a gpc for bh - ns mergers ( 445/927 mpc are adopted by the ligo - virgo collaboration as canonical values ; @xcite ) .    understanding the observable em signature of compact binary mergers has several observational implications . first ,
once the detectors are operational it is likely that the first detection of a gw signal will be around or even below threshold .
detection of an accompanying em signal will confirm the discovery , thereby increasing significantly the sensitivity of gw detectors @xcite .
second , the physics that can be learned from observations of a merger event through different glasses is much greater than what we can learn through em or gw observations alone .
finally , even before the detectors are operational , detection of em signature will enable us to determine the expected rates , a question of outmost importance for the design and the operation policy of the advanced detectors .
the current constraints on the rates are rather loose . the last ligo and virgo runs provided only weak upper limits on the merger rates : 8700 myr@xmath3(@xmath4 corresponding to @xmath5 yr@xmath3 gpc@xmath6 for ns@xmath1 and @xmath7 myr@xmath3 ( @xmath4 ( @xmath8 ) for bh - ns @xcite .
estimates based on the observed binary pulsars in the galaxy are highly uncertain , with values ranging from @xmath9 @xcite .
it has been suggested @xcite that short gamma - ray bursts ( grbs ) arise from neutron star merger events .
the estimated rate of short grbs are indeed comparable to binary pulsar estimates @xcite . however , while appealing , the association is not proven yet @xcite .
if correct , the observed rate of short grbs , @xmath10 , provides a lower limit to the merger rate .
the true rate depends on a poorly constraint beaming angle , resulting in an uncertainty of almost two orders of magnitude .
there are no direct estimates of bh - ns mergers , as no such system has ever been observed , and here one has to relay only on a rather model dependent population synthesis ( e.g. , * ? ? ?
* ; * ? ? ?
possible em signals from coalescence events were discussed by several authors .
there are several suggestions @xcite of a prompt ( coinciding with the gw signal ) short lived em signals , mostly in low radio frequencies , whose amplitudes are highly uncertain .
@xcite suggested that the radioactive decay of ejected debris from the merger will drive a short lived supernova like event .
@xcite calculated the radioactive heating during this process self - consistently .
they find that if @xmath11 is ejected then the optical emission from a merger at @xmath12 mpc peaks after @xmath13 day at @xmath14 .
if the mass ejection is lower then the optical emission will be even fainter .
finding , and especially identifying the origin of , such rare and faint events in the crowded variable optical sky is an extremely challenging task , even for current and future optical searches like ptf @xcite , panstarr , and lsst @xcite .
an intriguing possibility is that mergers produce short grbs @xcite .
however , short grbs are expected to be beamed , and only rarely this em signal will point towards us . a beamed grb that is observed off - axis will inevitably produce a long lasting radio  orphan " afterglow @xcite . a key point in estimating
the detectability of grb orphan afterglows is that the well constrained observables are the _ isotropic _ equivalent energy of the flow and the rate of bursts that point towards earth .
however , the detectability of the orphan afterglows depends only on the _ total _ energy and _ true _ rate , namely on the poorly constraint jet beaming angle . @xcite
have shown that while large beaming increases the true rate it reduces the total energy , and altogether reduces the detectability of radio orphan afterglows .
this counterintuitive result makes the detectability of late emission from a decelerating jet , which produced a grb when it was still ultra - relativistic , less promising .
however , regardless of amount of ulrtra - relativistic outflow that is launched by compact binary mergers , and of whether they produce short grbs or not , mergers are most likely do launch an energetic sub - relativistic and mildly - relativistic outflows .
the interaction of these outflows with the surrounding matter will inevitably produce blast waves and possibly stronger radio counterparts than that of ultra - relativistic outflows .
below we first discuss ( in  2 ) the current estimates of mass and energy ejection from compact binary mergers . in
3 we calculate the radio emission resulting from the interaction of this ejecta ( sub - relativistic , mildly relativistic and off - axis relativistic ) with the surrounding interstellar matter ( ism ) .
the calculations follow to a large extent models of radio ib / c supernovae @xcite and long grb radio afterglows @xcite .
the success of radio supernova ( sn ) modeling , where the observations are superb , indicates that the microphysics is well constrained ant that equipartition parameters describe well the physical conditions .
hence the main uncertainty in the predicted radio signal is in the amount of matter ejected from the merger event and its velocity .
luckily the estimates of this important quantity can be significantly improved even using existing numerical models .
we discuss the observational implications for detectability of merger remnants in  4 .
we estimate the expected rates of detection of different outflows in  4.1 .
we devote in  4.2 a special attention to short grb orphan afterglows that are a special case of our model , in which the outflow is launched relativistically , but the radio emission peaks only during its mildly relativistic phase .
the estimates of orphan afterglows detectability is independent of whether they are the products of binary mergers or not . in
4.3 we examine possible other radio transients that may hinder the identification of merger remnants .
finally in  4.4 we examine blind transient searched done in the past and we identify rt 19870422 as a possible and even likely merger remnant .
we conclude in  5 .
numerical simulations of compact binary mergers have been carried out by various groups with two different approaches .
some ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) use newtonian dynamics ( modified to allow for gravitational radiation back - reaction ) with detailed microphysics .
others ( e.g. , * ? ? ?
* ; * ? ? ?
* ) use full general relativistic dynamics with different levels of approximate microphysics , with or without mhd . in almost all ns@xmath1 simulations
one finds an accretion disk surrounding a rapidly rotating massive object that eventually collapses to a black hole .
an exception is the recent general relativistic simulations of @xcite who find no disks in some configurations .
the system lifetime is at least a few dozen milliseconds and possibly longer .
the fate of an accretion disk in a bh - ns merger is expected to depend on the mass ratio , the bh spin and the ns compactness . in some cases
the disk is very small while in others it is substantial ( e.g. , * ? ? ?
* ; * ? ? ?
all simulations find some form of relativistic or sub - relativistic mass ejection .
first , matter is ejected as tidal tails during the first stages of the merger . in bh - ns mergers
the ejected energy can be very high and its velocity is mildly relativistic .
for example @xcite find @xmath15 c ejecta with @xmath16 erg , where c is the light speed . in a ns@xmath1 mergers a lower , but yet significant
, amount of energy can be ejected ( e.g. , @xcite find @xmath17 erg ) at a lower velocities of @xmath18 c. this mass ejection is expected also if no significant disk is formed .
disk formation leads to several additional outflow sources .
first , neutrino heating drives a wind from the disk surface ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the energy in this wind is substantial with predictions ranging between @xmath19 erg and @xmath20 erg for 0.01 - 0.1 m@xmath21 disk .
the outflow velocity is @xmath22 c from the outskirts of the disk and it is increasing , possibly up to relativistic velocities , for wind that is ejected from close to the central object . the mass ejection becomes even stronger when neutrino heating shuts - off and the wind is driven by viscous heating and by he - synthesis @xcite , leading to an ejection of @xmath23 of the initial disk mass at @xmath18 c. additional energy source is neutrino - antinuetrino annihilation above the disk , which can deposit up to @xmath24 erg , into an amount of mass that is not well constrained , leading possibly to a relativistic outflow .
finally , more speculative , but yet very plausible , source of outflow are em processes that tap the rotational energy of the central object , such as the blandford - znajek mechanism @xcite .
these are likely to produce relativistic outflows with an energy that can be as high as @xmath25 erg , and are the most probable engines of short grbs , if those are produced by compact binary mergers .
the conclusion is that a significant mass and energy ejection is a prediction of almost all compact binary merger modelings . in ns@xmath1 mergers
an ejection of @xmath26 erg at @xmath18 c is a fairly robust prediction .
faster ejecta ( relativistic or mildly relativistic ) with energy @xmath27 is also quite likely from inner parts of the 0.01 - 0.1 m@xmath21 disk that is typically found in simulations .
the outflow from bh - ns mergers was explored only by a few authors , but it is also seems to be significant and potentially even more energetic and at faster velocities than the outflow from ns@xmath1 mergers .
consider a spherical outflow with an energy @xmath28 and an initial lorentz factor @xmath29 , with a corresponding velocity @xmath30 , that propagates into a constant density , @xmath31 , medium .
if the outflow is not ultra relativistic , i.e. , @xmath32 it propagates at a constant velocity until , at @xmath33 , it reaches radius @xmath34 , where it collects a comparable mass to its own : @xmath35 and @xmath36 where we approximate @xmath37 and ignore relativistic effects . here and in the following , unless stated otherwise
, @xmath38 denotes the value of @xmath39 in c.g.s . units . at a radius @xmath40
the flow decelerates assuming the sedov - taylor self - similar solution , so the outflow velocity can be approximated as : @xmath41    if the outflow is collimated , highly relativistic and points away from a generic observer , as will typically happen if the mergers produce short grbs , the emission during the relativistic phase will be suppressed by relativistic beaming .
observable emission is produced only once the external shock decelerates to mildly relativistic velocities and the blast - wave becomes quasi spherical .
this takes place when @xmath42 namely at @xmath43 . from this radius
the hydrodynamics and the radiation become comparable to that of a spherical outflow with an initial lorentz factor @xmath44 .
this behavior is the source of the late radio grb orphan afterglows @xcite .
our calculations are therefore applicable for the detectability of mildly and non - relativistic outflows as well as for radio orphan grb afterglows .
emission from newtonian and mildly relativistic shocks is observed in radio sne and late phases of grb afterglows .
these observations are well explained by a theoretical model involving synchrotron emission of shock accelerated electrons in an amplified magnetic field .
the success of this model in explaining the detailed observations of radio ib / c sne ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) allows us to employ the same microphysics here .
energy considerations show that both the electrons and the magnetic field carry significant fractions of the total internal energy of the shocked gas , @xmath45 .
these values are consistent with those inferred from late radio afterglows of long grbs ( e.g. , * ? ? ?
* ; * ? ? ?
the observed spectra indicate that the distribution of the accelerated electrons lorentz factor , @xmath46 , is a power - law @xmath47 at @xmath48 where @xmath49 in mildly relativistic shocks ( e.g. , the radio emission from grb associated sne and late grb afterglows ) and @xmath50 in newtonian shocks ( as seen in typical radio sne ; ( * ? ? ? * and references therein ) ) . the value of @xmath51 is not observed directly but it can be calculated based on the total energy of the accelerated electrons , @xmath52    the radio spectrum generated by the shock is determined by two characteristic frequencies .
one is @xmath53 the typical synchrotron frequency of electrons with the typical ( also minimal ) lorentz factor @xmath51 .
the other is @xmath54 , the synchrotron self - absorption frequency .
we show below that since we are interested at the maximal flux at a given observed frequency , @xmath54 may play a role only if it is larger than @xmath55 .
its value in that case is @xmath56    , and frequencies , @xmath54 and @xmath55 .
the arrows show the temporal evolution of the characteristic values .
the temporal dependence before @xmath33 is noted below / to the left of the arrows while the temporal dependence after @xmath33 is noted above / to the right of the arrows .
note that the evolution of @xmath55 and @xmath57 , marked only in the left spectrum , is relevant for both spectra .
the evolution of @xmath54 , marked only in the right spectrum , is correct only when @xmath58 and is therefore relevant only in that spectrum.,title="fig:",width=491 ] +    figure [ fig : spec ] illustrates the two possible spectra , depending on the order of @xmath54 and @xmath55 .
the flux at any frequency can be found using these spectra and the unabsorbed synchrotron flux at @xmath55 : @xmath59 where @xmath60 is the distance to the source ( we neglect any cosmological effects ) .
note that this is the real flux at @xmath55 only if @xmath61 ( see figure [ fig : spec ] ) .
as long as the shock is moving with a constant velocity i.e. , at @xmath62 , the flux across the whole spectrum increases ( see equations [ eq num]-[eq fm ] ) .
the flux evolution at later times depends on the spectrum at @xmath33 , namely on @xmath63 and if @xmath64 then possibly on @xmath65 the flux at that time can be found using the unabsorbed synchrotron flux at @xmath66 : @xmath67    consider now a given observed frequency @xmath68 .
we are interested at the light curve near the peak flux at this frequency .
there are three possible types of light curves near the peak corresponding to : ( i ) @xmath69 , ( ii ) @xmath70 and ( iii ) @xmath71 . where we define @xmath72 as the frequency at which this equality will never take place . in that case
@xmath73 is the frequency at which this equality would have happened if @xmath29 would have been large enough ( see figure [ fig : cases ] ) ] @xmath74 . in figure
[ fig : cases ] we show a schematic sketch of the time evolution of @xmath54 and @xmath55 and the corresponding ranges of @xmath68 in which each of the cases is observed .     and @xmath55 in two cases , @xmath75 ( top ) and @xmath64 ( bottom ) . also marked is the value of @xmath73 .
the vertical dashed line marks @xmath76 .
the ranges of @xmath68 at which each of the cases is observed is separated by horizontal dashed lines and marked on the right .
note that in the bottom panel @xmath55 and @xmath54 are not crossing each other at @xmath77 and only two types of light curves , cases ( i ) and ( iii ) , can be observed.,title="fig:",width=453 ] +    to estimate the time and value of the peak flux we recall , that at all frequencies the flux increases until @xmath33 . in case ( i ) ,
@xmath69 , the deceleration time , @xmath76 , is also the time of the peak .
the reason is that while @xmath57 increases , @xmath55 decreases fast enough so that @xmath78 decreases after @xmath33 .
note that in that case @xmath54 plays no role since it decreases after deceleration .
overall , in this case the flux peaks at @xmath33 and @xmath79 .
in the two other cases , ( ii ) and ( iii ) , @xmath80 and/or @xmath81 and the flux keeps rising at @xmath77 until @xmath82 or @xmath83 , whichever comes last . to find out which one of the two frequencies
is it , we compare @xmath68 with @xmath73 . at @xmath84 , @xmath55 decreases faster than @xmath54 .
therefore in case ( ii ) where @xmath85 , the last frequency to cross @xmath68 is @xmath55 and the peak flux is observed when @xmath86 . in case ( iii )
where @xmath87 , and the last frequency to cross @xmath68 is @xmath54 and the peak flux is observed when @xmath88 .
now , it is straight forward to calculate the peak flux , @xmath89 and the time that it is observed , @xmath90 , for different frequencies .
it is also straight forward to calculate the flux temporal evolution prior and after @xmath33 using equations [ eq num]-[eq fm ] and the relation @xmath91 which holds at @xmath62 and @xmath92 at @xmath84 .
the peak fluxes , the times of the peak and the temporal evolution of the three different cases are summarized in table [ table1 ] .
the overall different light curves are depicted in fig .
[ fig : lightcurves ]        .[table1 ] the observed flux before and after @xmath90 in the three different regimes . @xmath93 the temporal evolution only during the last power - law segment before @xmath90 . at earlier times
the temporal evolution may be different . [ cols="<,^,^,^,^ " , ]     the most sensitive radio facilities are at frequencies of @xmath94 ghz and higher .
equations [ eq numdec ] and [ eq nuadec ] imply that in this frequency range , for most realistic scenarios , it is a case ( i ) light curve , i.e. , @xmath95 . therefore , a newtonian and mildly relativistic outflows as well as relativistic grb orphan afterglows peak at @xmath33 with : @xmath96    the regime of @xmath97 at lower radio frequencies ( @xmath98 ghz ) depends on the various parameters . if the outflow is newtonian or the density is low or the energy is low then @xmath99 mhz and equation [ eq fpeak reg1 ] is applicable .
otherwise low radio frequencies are in regime ( iii ) , i.e. , @xmath100 .
the flux peaks in this case at @xmath101 with @xmath102 in the last two equations we used @xmath103 ( other @xmath104 values in the range 2.1 - 3 yield slightly different numerical factors and power laws ) .    to date , the best observed signal from a mildly relativistic blast waves is the radio emission that follows grb associated sne .
the main difference is that in these cases the circum burst medium is typically a wind ( i.e. , @xmath105 ) and therefore the density at early times is much larger then in the ism and self absorption plays the main role in determining the light curve .
a good example for comparison of equation [ eq fpeak reg1 ] with observations is the light curve of sn 1998bw .
this light curve is observed at several frequencies at many epochs , enabling a detailed modeling that results in tight constraints of the blast wave and microphysical parameters .
@xcite find that at the time of the peak at @xmath94 ghz , about 40 days after the sn , taking @xmath106 , the energy in the blast wave is @xmath107 erg , its lorentz factor is @xmath108 and the external density at the shock radius is @xmath109 .
the peak is observed when @xmath110 and it depends only on these parameters ( it is only weakly sensitive to the density profile ) . therefore , equation [ eq fpeak reg1 ] is applicable in that case .
indeed , plugging these numbers into equation [ eq fpeak reg1 ] we obtain a flux of 20 mjy at the distance of sn 1998bw ( 40 mpc ) , compared to the observed flux of 30 mjy .
this is not surprising given that the model we use is based on that of radio sne .
in the discussion above we considered an outflow with a characteristic energy and velocity . however a compact binary merger may produce an outflow with an energy dependent velocity , e.g. , @xmath111 . if we consider a case ( i ) signal ( @xmath69 ; e.g. , above @xmath112 ghz ) then we find @xmath113 .
therefore if @xmath114 the flux is dominated by the mildly relativistic ejecta , assuming that the relativistic part of the outflow is not pointing towards the observer .
otherwise the flux is dominated by the slowest ejecta .
for canonical parameters the strongest signal is expected around @xmath94 ghz , conveniently where the sensitivity of radio telescopes is high . if the signal peaks at lower frequencies it decreases from the peak as @xmath115 . since 1.4 ghz receivers are ten times more sensitive than lower frequency ones they are still more likely to detect a signal even in this case .
therefore , @xmath94 ghz is the optimal frequency to look for radio remnants of compact binary mergers and we therefore consider here the delectability at a @xmath94 ghz survey .
the number of events in a single whole sky snapshot is @xmath116 , where @xmath117 is the detectable volume at the survey flux limit @xmath118 , @xmath119 is the time that the flux is above the detection limit and @xmath120 is the event rate .
since at @xmath94 ghz the relevant light curve is case ( i ) ( see  [ sec theory ] ) we use equations [ eq : tdec ] and [ eq fpeak reg1 ] , and the approximation @xmath121 , to find that the number of radio coalescence remnants in a single @xmath94 ghz whole sky snapshot is @xmath122 where @xmath123 mjy and @xmath124 is the merger rate in units of @xmath125 .
since the microphysical parameters are reasonably constrained by radio sne , the main uncertainty in the signal detectability is the blast - wave properties and the circum - merger density .
the latter is expected to vary significantly between different merger events , dropping down to @xmath126 for mergers that take place outside of their host galaxies .
however , the observed galactic double - ns population reveals that a significant fraction of compact binary mergers ( at least ns@xmath1 ) take place in the disks of milky - way - like galaxies where @xmath127 .    as discussed in
 2 there is a good chance that ns@xmath1 mergers eject mildly relativistic outflows with @xmath128 erg . in that case the em counterparts of these gw sources are detectable by current facilities . a deep radio survey that covers a significant fraction of the sky will most likely detect their radio remnants to a distance of @xmath129 mpc .
these will appear as transients , varying on a several weeks time scale , with an optically thin spectrum at @xmath130 ghz at all time .
the radio remnants will be identified in random places within their host galaxies , which should be easily detectable at that distance .
it will not be accompanied by any optical counterpart with similar variability time scales . given a gw trigger with localization of @xmath131 , a deep search of the error box will detect a remnant in cases that the surrounding density is not very low .
ns@xmath1 mergers are also expected to eject energetic , @xmath132 erg , outflow at lower velocities , @xmath133 .
the detectability is very sensitive to the velocity , @xmath134 for @xmath103 , and to the energy , @xmath135 .
therefore , even if there is no mildly relativistic component to the outflow , the number of detectable remnants that are dominated by slow moving ejecta is expected to be significant .
for example , assuming that an energy of @xmath136 erg is ejected at @xmath137 we expect @xmath138 .
the variability time scale of these transients is expected to be @xmath139 yr .
if the velocity is instead @xmath140c the number of events drops by an order of magnitude and the variability timescale increases to @xmath141 yr .
these time scales increase the difficulty in the identification of the remnants as transients , and require a long term survey
. the spectrum of these transients will be optically thin also at low frequencies ( @xmath142 mhz ) .
other characteristics of these transients ( e.g. , location within the host ) are expected to be similar to that of a mildly relativistic outflow .
a gw triggered search increases of course the probability to find a remnant and a @xmath143 erg outflow can be detected up to @xmath12 mpc even if its velocity is @xmath144c .
the detectability of bh - ns mergers is much harder to predict .
first due to their virtually unconstrained rate and second since the properties of the outflow are less certain .
the latter can be significantly improved by current and future merger simulations that put focus on the ejected mass . in any case
the potential of these mergers to throw out a considerable amount of energy , @xmath17 erg , at mildly relativistic velocities can make them detectable to their gw detection horizon , which is much further than the ns@xmath1 horizon .    finally , if compact binary mergers launch also collimated ultra - relativistic outflows , and produce short grb , then orphan short grb afterglows are also part of the post merger em signal .
the detectability of radio orphan afterglows can be estimated based on observations of short grbs and is independent of whether they are produced by mergers or not .
we discuss their detectability in the next subsection .
all together , the range of the current predictions is rather large , but with most parameters we expect detectable radio signals .
some of the new generation radio telescopes have large fields of view ( e.g. , askap with @xmath145 and apertif with @xmath146 ; @xcite ) and improved sensitivities , making them ideal for large scale sub - mjy blind survey .
the evla , which has a smaller field of view but a remarkable sensitivity , is the best facility for gw triggered search .
it is also currently the fastest radio - survey instrument and it can carry - out a sub - mjy blind survey .
all these observatories have a very good chance to detect compact binary merger remnants in dedicated blind searches .
in fact with very reasonable parameters ( e.g. , @xmath147 erg of mildly relativistic ejecta ) a sub - mjy whole sky survey can detect thousands of binary - merger radio remnants .
short grb outflow begins highly relativistic and probably highly beamed .
eventually it slows down ( see  [ sec theory ] ) and become detectable from all directions .
therefore , the rate estimate equation [ eq rate ] is also applicable for radio orphan afterglows when @xmath148 .
however some of the parameters in equation [ eq rate ] are not directly observable .
the observed quantities are isotropic equivalent @xmath46-ray energy , @xmath149 , and the rate of bursts that point to the observer @xmath150 , while equation [ eq rate ] depends on @xmath151 and @xmath152 , where @xmath153 is the fraction of the @xmath154 steradian covered by the jet and @xmath155 is the isotropic equivalent energy in the afterglow blast wave .
x - ray observations indicate that @xmath46-ray emission in short grbs is very efficient and that in general @xmath156 @xcite .
we assume that this is the case in the following discussion .
@xmath149 of short grbs ranges at least over four orders of magnitude ( @xmath157 erg ) .
the rate of observed short grbs is dominated by @xmath19 erg bursts , and the luminosity function can be well approximated by a power - law , at least in the range @xmath158 erg , such that @xmath159 where @xmath160 @xcite . plugging these into equation [ eq rate ] we obtain @xmath161
this equation is similar to equation 9 of @xcite , with the observed luminosity function already folded in .
narrower beamed bursts ( with lower @xmath162 ) are more numerous and they produce less total energy per burst .
the positive dependence of equation [ eq orphan_rate ] on @xmath162 implies that overall the lower energy is  winning " over the increased rate , and the detectability of narrower bursts is lower .
using , equation [ eq orphan_rate ] we can put a robust upper - limit on the orphans rate since all the parameters are rather well constraint by observations , with the exception of @xmath162 which is @xmath163 by definition .
therefore , assuming that short grbs are beamed , the detection of the common @xmath24 erg bursts in a blind survey , even with next generation radio facilities , is unlikely @xcite .
however , brighter events should be detectable .
if the beaming is energy independent , detectability increases with the burst energy .
the luminosity function possibly breaks around @xmath20 erg , in which case the orphans number is dominated by @xmath20 erg bursts . for @xmath164 we expect , from these bursts @xmath141 orphan afterglows at a 0.1 mjy in a single @xmath94 ghz whole sky snapshot .
so far we discussed detectability in a blind survey .
a followup dedicated search would be , of course , more sensitive .
if compact binary mergers produce short grbs than the energy of most gw detected bursts will be faint with @xmath165 erg .
the chance to detect their orphan afterglows again depended on their total energy and thus on @xmath162 .
equation [ eq fpeak reg1 ] shows that if @xmath164 then detection should be difficult but possible in a dedicated search mode . note
that since the energy of the burst is low , the radio emission will evolve quickly , reaching a peak and decaying on a week time scale , so a prompt and rather deep search will be needed .
@xcite carried out a 5 ghz survey looking for transients on timescales of a week to a year .
the survey sensitivity for transients with variability scale @xmath166 day is 0.37 mjy with an effective are of @xmath167 .
events with variability time scale of two months where surveyed at sensitivity of 0.2 mjy with an effective area of @xmath168 .
@xcite report the detection of 10 transients .
the most interesting of those , in our context , is rt 19870422 , which has a variability time scale of two months .
it is found within a star forming galaxy at a distance of 1.05 gpc , but at a significant offset from the host nucleus .
its luminosity and time scale are those expected from a @xmath169 erg mildly relativistic outflow that propagates in the ism .
it is , therefore , a prime candidate for a compact binary merger radio remnant .
based on this single event @xcite infer a best estimate rate of @xmath170 for rt 19870422-like events .
taking a @xmath171 poisson error , the best estimate translates to a range of @xmath172 , fully consistent with the estimates of compact binary mergers .
if this is indeed a merger remnant then , since for optically thin spectrum the fluxes at 1.4 ghz and 5 ghz are not very different , a sub - mjy 1.4 ghz whole sky survey would detect hundreds to thousands of radio remnants .
@xcite suggested that this transient is a radio sne similar to sn 1998bw , but brighter .
this is certainly a viable possibility , however , if true then this radio sn is brighter by an order of magnitude than the brightest radio sn ever observed before .
unfortunately , lacking optical search for a sn or a multi - wavelength measurement that determines the transient spectrum , it is impossible to rule out any of the two possibilities .    an additional interesting candidate is rt 19840613 .
it is variable on less than 7 days and it has a host galaxy at a distance of 140 mpc . even assuming a variability time scale of 7 days it is marginal as a merger remnant candidate . assuming 7 days variability @xcite find a rate of @xmath173 , which is again only marginally consistent with current estimates of compact binary merger rates .
therefore , while this may be a merger remnant it is not a very promising candidate .
@xcite suggest that this is also a sn 1998bw - like event . while this possibility can not be ruled out , the inferred rate is at least an order of magnitude larger than that of 98bw - like events ( note that radio sne as bright as sn 1998bw are very rare compared to typical radio sne ) .
the additional 8 events detected by @xcite have no clear host galaxies and are therefore probably not merger remnants .
a key issue with the detection of compact binary merger remnants in blind surveys is their identification . @xcite
present a census of the transient radio sky .
luckily the transient radio sky at @xmath94 ghz are relatively quite .
the main contamination source are radio active galactic nuclei ( agns ) , however their persistent emission is typically detectable in other wavelength and/or deeper radio observations .
moreover , the signal from a compact binary merger is expected to be located within its host galaxy ( otherwise the density is too low ) , but away from its center .
the host and the burst location within it , should be easily detectable at the relevant distances .
the only known , and guaranteed , transient 1.4 ghz source with similar properties are radio sne . among these typical radio sne
are the most abundant .
transient search over 1/17 of the sky with @xmath174 mjy at 1.4ghz @xcite finds one radio sn .
this rate translates to @xmath175 sne in a whole sky @xmath176 mjy survey .
these contaminators can be filtered in three ways .
first , by detection of the sn optical light
. however , the optical signal may be missed if it is heavily extincted , and given the large number of radio sne , misidentifying even a small fraction of them may render the survey useless for our purpose .
the second filter is the optically thick spectrum at high radio frequency ( @xmath141 ghz ) at early times , which is a result of the blast wave propagation in a wind .
thus , a multi - wavelength radio survey can identify radio sne .
the last filter is the luminosity - time scale relation of typical radio sne that is induced by the outflow velocity ( e.g. , figure 2 in * ? ? ?
type ii sn outflows are slow , @xmath177c , and therefore their radio emission is longer / fainter than that expected for merger remnants . the common type of ib / c radio sne is produced by @xmath178c blast waves but with much less energy than what we expect from a binary merger outflow , and therefore their radio emission is much fainter
. the combination of any two of these filters will hopefully be enough to identify all the typical radio sne .
slightly different contaminators are grb associated sne .
their outflows is as fast and as energetic as those that we expect from a binary merger and therefore their radio signature is similar in time scales and luminosities .
sn1998bw - like events are detectable by a 0.1 mjy survey at 1.4 ghz up to a distance of several hundred mpc for 40 days and their rate is @xmath179 @xcite , implying at least several sources at any whole sky snapshot . here
only the first filter ( sn optical ligth ) and possibly the second ( optically thick spectrum ) can be applied .
however , given the high optical luminosity of grb associated sne and their relatively low number this should be enough in order to filter them out .
these contaminators highlight the importance of a multi - wave length strategy where an optical survey accompany the radio survey to best utilize both surveys detections .
the results of @xcite implies that thousands of sources with properties similar to rt 19870422 are expected in a 5 ghz sub - mjy all sky survey
. if these events are optically thick during their whole evolution than these are not merger remnants and they could be easily filtered out . moreover , in that case their rate in a 1.4 ghz survey should be lower by two orders of magnitude .
if these events show a synchrotron optically thin spectrum and no optical counterparts , then they should be abundant also in a 1.4 ghz survey , and their origin is most likely compact binary mergers .    finally , radio is the place to look for blast waves in tenuous mediums , regardless of their origin .
any source of such explosion , being a binary merger , a grb or a sn , produces a radio signature .
therefore , all the strong explosions may be detectable is a deep radio survey , this include for example long grb on - axis and off - axis afterglows and giant flares from extra galactic soft gamma - repeaters .
the difference between the radio signatures of the different sources ( amplitude , spectrum and time evolution ) depends on the blast wave energy and velocity and on the external medium properties .
we thus will be able to identify the characteristics of binary mergers outflows . if , however , there is a different source of @xmath169 erg of mildly relativistic outflow that explodes in the ism it will be indistinguishable from binary mergers ( at least not in the radio ) .
currently we are not aware of any such source , with the exception of long grbs at the low end of the luminosity function , but these are too rare to contaminate a survey .
any other source of such outflows , if exist , will probably be a part of the family of collapsing / coalesing compact objects .
compact binaries are expected to eject sub - relativistic , mildly relativistic and possibly ultra - relativistic outflows as part of their merger process .
we have shown that these outflows will inevitably produce a long lived radio remnant .
these are the most robust predictions of an em counterpart of the merger gw signal .
the radio remnant appears weeks to years after the merger and remains bright for a similar time .
therefore , a trigger following a detection of gw signal can wait for a week after the event and no online triggering is needed . in addition
the long lasting remnants enable a detection in a blind survey .
for mildly relativistic outflows with @xmath19 erg that propagate in the ism we expect a few weeks radio transients with a 1.4 ghz flux of @xmath180 0.3 mjy from sources at @xmath129 mpc , the advanced ligo - virgo horizon for ns@xmath1 mergers .
the bh - ns gw horizon is farther , but current numerical simulations suggest they involve higher energy outflow resulting in a comparable flux .
follow up observations of gw candidate events , at a level of @xmath181jy are feasible and are very likely to show a radio transient for either ns@xmath1 merger of bh - ns merger .
we find that the optimal frequency to carry out a search for merger remnants is @xmath94 ghz .
assuming a mildly relativistic outflow with @xmath19 erg the canonical ns@xmath1 merger rate of 300 gpc@xmath6 yr@xmath3(and a range of 20 - @xmath182 gpc@xmath6 yr@xmath3 ) implies a detection of @xmath183 ( 1 - 1200 correspondingly ) radio ns@xmath1 remnants in a @xmath140 mjy all sky survey .
this rate depends quadratically on the outflow energy , so a very plausible ejected energy of @xmath136 erg increases the rate by two orders of magnitude , making them detectable even in a survey that covers only a small fraction of the sky or that is at a mjy sensitivity .
therefore carrying out a large field - of - view and sensitive ghz survey by currently available facilities has a great potential to constrain the rate of binary mergers , a piece of information that is of great importance for the design and operation of gw detectors .
even if mergers do not launch a significant mildly relativistic ejecta they are still expected to produce an energetic ( @xmath184 erg ) sub - relativistic ( 0.1 - 0.2)c outflows .
these outflows will also produce radio remnants .
these remnants will be fainter , detectable only to a distance of @xmath142 mpc at 0.1 mjy , and will evolve more slowly , on time scales of 3 - 10 yr .
these transients are also detectable at a rate of @xmath141 over the whole sky at 0.1 mjy , although identifying them as transients is harder and it requires a long term survey .
we estimate the detectability of short grb orphan afterglows , which may also be produced by compact binary mergers if they are launching also ultra - relativistic outflows .
these estimates are based on short grb observations and are therefore indifferent to whether short grbs are binary mergers or not .
the main uncertainty in the rate estimates is the grb beaming factor .
we find that assuming @xmath164 there are expected to be about 10 orphan afterglows at a 0.1 mjy in a single @xmath94 ghz whole sky snapshot .
the duration of these afterglows is several weeks
. if binary mergers are short grbs than a gw triggered event will most likely be of a low energy grb , @xmath165 erg , and a true energy , after beaming correction , that is even lower .
the radio orphan afterglow will probably still be detectable in a deep search .
however its variability time scale is short , about a week , so the search should be done promptly .
remarkably , the observed 5 ghz transient rt 19870422 , detected by @xcite fits very well our estimates of the typical expected properties of a compact binary merger radio remnant . at a distance of 1 gpc and a duration of two months
this transient is what expected from a mildly relativistic outflow with @xmath169 erg .
the rate inferred from this single event is also fully consistent with that of ns@xmath1 mergers .
this transient is an excellent candidate to be the first observed radio remnant of a merger .
unfortunately , one can not rule out the possibility , suggested by @xcite , that this is an especially bright radio sn .
note , however , that this interpretation requires a sn brighter by an order of magnitude than any radio sne previously observed .
simultaneous optical observations or multi - wavelength radio observations could have easily distinguished between the two possibilities .
the first could have determined if there was a sn or not .
the second could have distinguished between an optically thick radio spectrum expected in radio sne vs the optically thin spectrum expected in merger remnants at these frequencies at all times .
unfortunately , no such observations were available .
however , the rate implied by this even is very high and similar events should be detected in a sub - mjy survey of even a small fraction of the sky .
therefore , the nature of this type of events can be easily probed with current facilities .
our results show the great potential of @xmath94 ghz radio transient observations at the sub - mjy level for the detection of ns@xmath1 mergers . on the observational side
these predictions provide an excellent motivation for carrying out a whole sky sub - mjy survey using the evla or other upcoming radio telescopes .
the main source of contamination in such surveys would be radio supernova and those could be distinguished from compact binary mergers by their optical signal , spectrum and other characteristic properties .    while it is clear that compact binary mergers produce sub - relativistic to relativistic outflows , details of those outflows are not well determined at present .
this is to large extent because of lack of interest rather than because of specific difficulties in analyzing their properties .
our analysis elucidate the importance of a detailed quantitative estimates concerning these outflows , a task that is within the scope of current simulations .
we thank dale frail , shri kulkarni , andrew macfadyen , eran ofek and stephan rosswog for helpful discussions .
this research was supported by an erc advanced research grant , by the israeli center for excellent for high energy astrophysics , by the israel science foundation ( grant no .
174/08 ) and by an irg grant .


within the shubnikov phase of type ii superconductors the applied magnetic field enters the sample in the form of flux lines .
the standard mean - field type calculation@xcite shows that in an isotropic material two straight vortices repel at all distance scales , with an interaction strength @xmath1 , @xmath2 being the basic energy scale in the vortex matter , @xmath3 is the zero - order modified bessel function , @xmath0 is the inter - vortex distance , and @xmath4 the magnetic penetration length ( @xmath5 denotes the flux quantum ) .
however , it has recently been shown @xcite that in layered and strongly anisotropic superconductors the thermal fluctuations of the flux lines give rise to a long range attraction @xmath6 of the van der waals type between the vortices , where @xmath7 denotes the temperature and @xmath8 is the interlayer separation .
the strongly fluctuating and layered high temperature superconductors are particularly well suited to exhibit this attractive component in the vortex - vortex interaction .
alternatively , the attraction is induced through static vortex distortions due to an underlying pinning landscape , an effect recently studied by mukherji and nattermann@xcite and by volmer _
et al._@xcite    following a suggestion of nelson@xcite , the statistical mechanics of vortices can be mapped to the imaginary time quantum mechanics of two - dimensional ( 2d ) bosons .
the particular interaction between the flux lines renders the bosons charged ( with a charge screened on the scale of the london penetration depth @xmath4 ) .
this type of long range interaction can be formulated in terms of a massive gauge field theory@xcite . within the resulting 2d massive electrodynamics ,
the vortex matter acts as a dielectric medium and we can define a casimir problem , see fig .  1
: under the vortex @xmath9 boson mapping two half - spaces of vortices separated by a gap of width @xmath0 act as two dielectric planar media which attract each other via a casimir force .
= 8.5 cm    in the present work , we determine the casimir force between two dielectric half planes of charged bosons . for dilute media ,
the macroscopic casimir force can be related to the microscopic van der waals force between the media s constituents via pairwise summation . here
, we present a derivation of the van der waals force between vortex lines via this alternative route , calculating first the casimir force between two dielectric planar media in the boson picture and then reconstructing the van der waals force in the reductionist way .    since the discovery of the casimir effect@xcite in 1948 ,
several hundred papers have dealt with this phenomenon , disseminating its fascination into many branches of physics @xcite .
casimir forces between macroscopic bodies are a quantum effect caused by a shift in the zero - point energy of gauge field fluctuations such as the electromagnetic one .
the casimir effect is bound to a number of system properties , such as topology and dielectric permittivity , and the reduction to an analogous van der waals interaction is not always possible . however , the interpretation in terms of a van der waals attraction is possible for the case of rarefied media and appropriate geometries , such as the parallel plate setup@xcite . in our derivation of the van der waals force from the casimir effect
we will make use of these special conditions .    in the following
, we discuss the relationship between the casimir effect and the van der waals interaction within a path integral formulation ( sec .
ii ) before deriving the appropriate action for the 2d charged bosons from the london functional describing the vortices ( sec .
iii ) . in section iv , we briefly review the derivation of the van der waals force in the original vortex language and then proceed with the calculation of the 2d casimir effect in section v , the main section of the paper containing the new results .
we consider two parallel material slabs made from fluctuating dipoles and separated by a vacuum gap .
summing pairwise over all microscopic van der waals interactions between the dipoles provides the macroscopic casimir interaction between the slabs . on the other hand
, we can determine the dielectric properties of the individual slabs as produced by the fluctuating dipoles .
the specific boundary conditions due to the dielectric properties of the slabs influence the spectrum of the electromagnetic field confined in between . the change in the spectrum as a function of the separation of the slabs produces the casimir force .
this analogy between the casimir and the van der waals force is transparently brought out within a path integral formulation , see fig .  2 : assume the system under consideration can be described by an action @xmath10 $ ] depending on the gauge field @xmath11 and a particle current @xmath12 . carrying out the partial integration in the partition function @xmath13{\cal d}[{\bf j}]\ , e^{-{\cal s } [ { \bf j},{\bf a}]/\hbar}$ ] over the matter field @xmath12 or over the gauge field @xmath11
, we obtain an effective action @xmath14 describing the conjugate field alone : from the effective action @xmath15 $ ] describing the gauge field we can derive the casimir effect , while the interaction of the particle currents as described by @xmath16 $ ] will give us the van der waals attraction .
the casimir- and van der waals forces then can be related to one another via the pairwise summation of the interparticle forces@xcite : consider two @xmath8-dimensional homogeneous macroscopic dielectric bodies of density @xmath17 with parallel interfaces separated by a distance @xmath0 .
the two bodies attract one another due to a microscopic particle - particle interaction @xmath18 , @xmath19 , of the van der waals type .
the interaction energy can be written in the form @xmath20 where @xmath21 is a hypercube of size @xmath22 parallel to the interface and @xmath23 , with @xmath24 the @xmath25 dimensional in - plane coordinate , while @xmath26 is the coordinate along the direction perpendicular to the plane .
we then find for the casimir force density@xcite @xmath27 the result @xmath28 } { \gamma(\alpha/2 ) }                  \frac{1}{\alpha - d } \frac{1}{r^{\alpha - d}}. \label{f_of_r}\end{aligned}\ ] ] the result ( [ f_of_r ] ) then allows to infer the parameters @xmath29 and @xmath30 , characterizing the van der waals interaction @xmath31 , from the macroscopic casimir force density @xmath32 .
( 8.5,5 ) ( 0.1,0.1)(8.3,4.7 ) ( 2.25,2)(0,-1)0.75 ( 6.25,2)(0,-1)0.75 ( 4.25,0.75)(-1,0)1 ( 4.25,0.75)(1,0)1 ( 3.5,4)(-1,-1)1 ( 5,4)(1,-1)1 ( 3.8,4.25)@xmath10 $ ] ( 1.8,2.4)@xmath33 $ ] ( 5.9,2.4)@xmath16 $ ] ( 0.75,0.7)casimir effect ( 0.75,0.35)in dielectrics ( 5.75,0.7)van der waals ( 5.75,0.35)attraction ( 3.4,0.35)summation ( 3.6,1)pairwise ( 0.9,3.75)integrate ( 0.9,3.4)over @xmath12 in @xmath34 ( 6,3.75)integrate ( 6,3.4)over @xmath11 in @xmath34    after mapping the thermally fluctuating vortex matter to a system of 2d quantum charged bosons , we will arrive at an action @xmath35 $ ] with two gauge fields @xmath11 and @xmath36 , see eq .
( [ edr ] ) .
the casimir effect we are interested in here is the one produced by the fake gauge field @xmath11 : the integration over the matter field @xmath12 will produce the dielectric properties of the media , while the integration over the physical gauge field @xmath37 renders the fake field massive and hence will always exponentially confine the casimir force to finite distances .
we start from the london free energy in an isotropic superconductor@xcite @xmath38 = \frac{1}{8\pi}\int { \rm d}^3r \ , [ { \bf b}^2 + \lambda^2(\nabla \times { \bf b})^2 ] , \label{london}\ ] ] with @xmath39 the magnetic field and @xmath4 denoting the london penetration depth . in order to acquaint for the vortices ,
we add the current term @xmath40 with @xmath41 where @xmath42 , @xmath43 , the coordinates @xmath44 denote the position , the vectors @xmath45 the direction of the vortex lines , and @xmath5 is the unit of flux . ignoring screening , the interaction between the vortex lines is long ranged and thus can conveniently be expressed through a mediating gauge field @xmath11 . as usual , we introduce the gauge field @xmath11 as an auxiliary field such that @xmath46 ) = \exp(-\beta { \cal f})$ ] , with @xmath47 the inverse temperature@xcite ( we set the boltzmann constant @xmath48 to unity and fix the gauge through the condition @xmath49 ) , @xmath50 & = & \int { \rm d}^3{x } \ , \left[i{\bf a}\cdot\left({\bf j}-\frac{1}{\phi_0}(\nabla\times{\bf a})\right ) \right .
\label{fprime } \\ & & \quad + \left .
\frac{1}{8\pi}(\nabla \times { \bf a})^2 + \frac{1}{2g^2}(\nabla \times { \bf a})^2 \right ] , \nonumber\end{aligned}\ ] ] where @xmath51 and the line energy @xmath52 is the basic energy scale in the problem . in ( [ fprime ] ) , we have accounted for screening by introducing back into the model the real gauge field @xmath37 .
the anisotropy of uniaxial layered materials is most conveniently introduced through a rescaling of the scalar fake magnetic field @xmath53 @xcite , @xmath54 where the anisotropy factor @xmath55 is determined through the effective masses perpendicular ( @xmath56 ) and parallel ( @xmath57 ) to the @xmath58-plane ( the fake electric field @xmath59 remains unchanged ; the subscript ` @xmath60 ' identifies the planar component of a vector , @xmath61 ) .
we map the statistical mechanics of the vortex system to the imaginary time quantum mechanics of 2d bosons through the replacements@xcite @xmath62 ( imaginary boson time ) , @xmath63 ( the boson s planck constant ) , @xmath64 ( the boson action ) , and @xmath65 ( the light velocity in the boson system ) .
the boson partition function is given by @xmath66 { \cal d}[{\bf a}]{\cal d}[{\bf a } ] e^{-{\cal s}[\{{\bf r}_\mu\},{\bf a},{\bf a}]/\hbar^{\scriptscriptstyle b } } , \label{zzz}\ ] ] with @xmath67 and @xmath68 = \int \!\ !
{ \rm d}\tau \sum_\mu \left[\frac{m}{2 } \left(\partial_\tau { \bf r}_\mu ( \tau)\right)^2 - \mu^{\scriptscriptstyle b } \right],\label{plm}\\ & & { \cal s}_{\rm int } [ \{{\bf r}_\mu\},{\bf a},{\bf a } ] = \int \!\ ! { \rm d } \tau { \rm d}^2 r \left[i{\bf a}\cdot \left({\bf j}-\frac{1}{\phi_0 } ( \nabla \times { \bf a})\right ) \right . \label{edr } \\ & & \quad + \left
. \frac{1}{8\pi } ( \nabla \times { \bf a})^2 + \frac{1}{2g^2}\left((\nabla\times{\bf a})_{xy}^2 +   \frac{1}{\varepsilon^2 } ( \nabla\times { \bf a})_\tau^2\right ) \right ] .
\nonumber\end{aligned}\ ] ] here , @xmath69 is the chemical potential of the bosons , @xmath70 is the boson density , and @xmath71 is the boson mass ( we have introduced a term @xmath72 to account for the external magnetic field @xmath73 producing the vortices in the superconductor ; @xmath74 denotes the planar coherence length ) . in the thermodynamic limit , we have @xmath75 , corresponding to @xmath76 , i.e. , we are interested in the ground state physics of the boson system .
note that the free boson action @xmath77 contains the bare mass @xmath78 due to the vortex core energy .
the retarded self - interaction of the bosons via their gauge fields then produces the mass renormalization @xmath79 , where @xmath80 is the dispersive line tension@xcite of the vortex lines .
to set the stage , we briefly review the derivation of the van der waals interaction as presented in ref .  ; this will allow us to fix some flaws in the previous derivation and will provide us with a check on the results for the casimir force derived later .    for a simple qualitative analysis we consider two vortices in a layered superconductor ( with layers separated by the distance @xmath8 ) .
ignoring the coupling between the layers , the fluctuating pancake vortices interact via a logarithmic potential @xmath81 .
second order perturbation theory then provides us with a van der waals interaction @xmath82 , the energy scale @xmath83 being set by the driving action of the thermal fluctuations . at long distances @xmath84 a finite interlayer coupling changes this result as the josephson interaction reduces the cutoff @xmath85 on the @xmath86 modes to the new value @xmath87 ; the van der waals interaction crosses over to @xmath88 .
the two results have their analogue in the van der waals attraction between neutral atoms , where the interaction potential exhibits a crossover from a @xmath89 at short to a @xmath90 behavior at large distances@xcite .
we proceed with the derivation of the long range van der waals interaction between vortex lines / bosons . following the scheme in fig .   2
, we integrate over the gauge fields @xmath37 and @xmath11 in the partition function ( [ zzz ] ) and obtain the effective current - current interaction @xmath91 =   \frac{\varepsilon_0}{2}\int { \rm d}^3r \,{\rm d}^3 r'\ ,   { j}_\alpha({\bf r})v_{\alpha\beta}^{\rm int}({\bf r}-{\bf r } ' ) { j}_\beta({\bf r } ) \label{aa1}\ ] ] ( we return to the more natural statistical mechanics notation in this section , @xmath92 ) . inserting the expression ( [ current ] ) for the currents , we can cast eq .
( [ aa1 ] ) into the standard form@xcite @xmath93 = \frac{\varepsilon_0}{2 } \sum_{\mu,\nu } \int { \rm d}r_{\mu\alpha } { \rm d}r'_{\nu\beta }   v_{\alpha\beta}^{\rm int } ( { \bf r}_\nu - { \bf r}'_\mu ) , \label{aa2}\ ] ] with the vortex positions @xmath94 and the interaction potential @xmath95 , conveniently expressed within a fourier representation , @xmath96 with @xmath97 .
\label{vab}\end{aligned}\ ] ] here @xmath98 , @xmath99 and @xmath100 , @xmath101 .
the free energy ( [ aa2 ] ) can be split into the self energy part @xmath102 ( @xmath103 ) and the interaction part @xmath104 ( @xmath105 ) . restricting ourselves to two vortices a distance @xmath0 apart
, we obtain @xmath106}e^{ik_z(z_1-z_2 ) } \nonumber \\ & \times & [ v_{zz}^{\rm int}({\bf k } ) + t_{1\alpha}(z_1)t_{2\beta}(z_2 ) v_{\alpha \beta}^{\rm int}({\bf k } ) ] , \end{aligned}\ ] ] where we have split the vortex positions @xmath107 into a mean field part @xmath108 and a fluctuating part @xmath109 .
up to a constant , the free energy is given by @xmath110\rangle_0.\ ] ] the average @xmath111 has to be taken with respect to the self - energy @xmath102 of the free vortices . performing a cumulant expansion
, we obtain the effective vortex - vortex interaction in the form @xmath112/2t$ ] , where @xmath113 denotes the sample thickness .
splitting into longitudinal ( to the induction , @xmath114 ) and transverse ( @xmath115 ) parts , the longitudinal interaction produces the standard repulsive vortex - vortex interaction@xcite , to lowest order in @xmath116 , @xmath117 while higher orders in @xmath116 merely renormalize the prefactor .
the transverse part produces the van der waals interaction@xcite @xmath118 using the decomposition @xmath119 we can reduce the average in ( [ vdw_def ] ) to the simpler form @xmath120 ^ 2 \langle t^2(k_z ) \rangle_0 ^ 2,\ ] ] with the partial fourier transform @xmath121 in strongly anisotropic and layered material the single vortex mean squared amplitude of fluctuations @xmath122 is limited by the electromagnetic interaction through the dispersive elasticity @xmath80 , @xmath123 with @xmath124 the evaluation of the partial fourier transform ( [ vpartial ] ) is carried out in appendix a and making use of the results ( [ app_a_x ] ) and ( [ app_a_y ] ) in the limit @xmath125 , we find the van der waals interaction in the decoupled limit @xmath126 the interlayer distance @xmath8 provides the large @xmath86 cutoff on the thermal fluctuations of the individual lines .    in the continuous anisotropic case the single vortex
elasticity is still dominated by its electromagnetic contribution . at intermediate distances
@xmath127 we recover again the result ( [ v4_vdw ] ) as the josephson coupling is not effective yet . at large distances @xmath128 , however , the interlayer coupling becomes important through cutting off the @xmath86 modes at @xmath87 and we find the result @xmath129 these results coincide with the ones found by blatter and geshkenbein@xcite up to the numerical prefactor as they missed a factor @xmath130 in eq .
( [ t^2 ] ) , as well as a term in the partial fourier transform which contributes to the result ( [ v5_vdw ] ) to order @xmath131 .
the same mistakes show up in their results for the attraction of the vortex to the sample surface , eqs .
( 19 ) and ( 20 ) of ref .  .
we give the correct expressions here : @xmath132 for the decoupled case , and @xmath133 in the continuous anisotropic situation .
in this section , we derive the casimir force between two parallel cubes of vortex matter and then show , via the method of pairwise summation , that the results coincide with those found by means of the van der waals approach .
again , we begin with a simple analysis based on dimensional estimates , instructing us what to expect .
consider two @xmath8-dimensional hypercubes of size @xmath134 separated by a distance @xmath135 , see fig .  1 .
we wish to determine the casimir force @xmath32 per unit area acting between the interfaces .
being due to the zero point fluctuations of the gauge field , the only relevant dimensional quantities entering the expression for the force are planck s constant @xmath136 , the velocity of light @xmath137 , and the distance @xmath0 between the slabs .
the combination @xmath138 then provides us with a dimensionally correct expression for the force ; in two dimensions we then find @xmath139 for the retarded casimir effect . in the non - retarded case , @xmath140 and we have to replace the combination @xmath141 by the frequency scale @xmath142 where the matter becomes transparent ( for the vortex matter in a layered superconductor this ` frequency ' is given by the layer separation @xmath8 , @xmath143 ) .
dimensional estimates then give us for the force per unit length in two dimensions @xmath144 . in order to derive the corresponding microscopic inter - particle potential from these estimates ,
we apply the method of pairwise summation and obtain the result @xmath145 in the retarded case and @xmath146 in the non - retarded limit .
these simple estimates produce the correct power laws for the van der waals potential ; in order to find the correct sign and the complete prefactors , we have to go through the detailed analysis below .
the casimir energy is given through the difference in the sum over all cavity modes minus the free field contribution@xcite .
while the calculation is rather straightforward for the simple case of a metallic cavity , an appropriate formalism has to be set up to treat more general configurations involving dielectric media .
we start from a quadratic lagrangian density @xmath147 in @xmath8 dimensions , with @xmath148 a vector boson field with @xmath149 components and @xmath150 the green function matrix .
the partition function @xmath34 is expressed through the usual imaginary time path integral formalism @xmath151\ ,          \exp\left[-\frac{1}{\hbar }           \int { \rm d}^d x \int_0^{\hbar / t }           \!\!\!\!\!\!\!\!d\tau \ , { \cal l}[{\bf a}]\right ]           = [ \det { \bf g}]^{-1/2},\ ] ] and the free energy @xmath152 is given by @xmath153 we consider the classic parallel - plate geometry with isotropic media separated by a gap of width @xmath0 along the @xmath154-direction . for each polarization
@xmath155 we determine the eigenstates obeying the boundary conditions at the @xmath156-dimensional hypersurfaces placed at @xmath157 and @xmath0 .
the individual modes are characterized by their polarization @xmath158 , the ( matsubara ) frequency @xmath159 , @xmath160 , the transverse wave vector @xmath161 , and the longitudinal wave vector @xmath162 , @xmath163 .
given the set @xmath164 , the discrete wave vectors @xmath162 are obtained as the solutions of the boundary condition @xmath165 , where @xmath166 is a short hand for the other indices @xmath167 ( the function @xmath168 derives from the determinant associated with the set of boundary conditions ; for the classic casimir effect , @xmath169 with @xmath0 the distance between the metal plates ) . expressing the trace through all the indices we have @xmath170 where the prime on the sum indicates that we count the @xmath171 term with a weight @xmath172 . in ( [ f_ln ] )
we have made explicit the boundary condition at the hypersurfaces separating the media from the gap .
= 5.5 cm    we make use of cauchy s theorem to perform the sum over the longitudinal momenta @xmath162 : we first rewrite the sum in the form @xmath173          = \int { \rm d}y \sum_{d_{\{\alpha\}}(k)=0 }            \frac{1}{g_{\{\alpha\}}^{-1 } + y}.\ ] ] next , we let the boundary condition @xmath165 produce the desired sum in a cauchy loop integral in the complex @xmath174 plane , see fig .  3 .
with @xmath175 we have @xmath176 where it is understood that both @xmath177 and @xmath178 further depend on the other indices @xmath166 as well as on the parameter @xmath179 in the case of @xmath178 .
we define the function @xmath180 as the zeroes of the expression @xmath181   + y = 0 $ ] . taking the derivative of the last equation with respect to @xmath179 , @xmath182 , and using the result in ( [ c_loop ] ) , we obtain the desired result @xmath183 , \label{sum_k}\ ] ] where we have used that the green function exhibits the proper asymptotics @xmath184 and the symmetry @xmath185 . inserting the sum on @xmath174 , eq .
( [ sum_k ] ) with @xmath186 , back into the free energy expression ( [ f_ln ] ) we obtain @xmath187 , \label{f_lnd}\ ] ] where @xmath188 \equiv l^{d-1}\sum_{\nu = 1}^r          { \sum^{\infty}_{s=0 } } \!\!\phantom{|}^\prime \int \frac{{\rm          d}^{d-1}q}{(2\pi)^{d-1 } } \ , [ \dots]\ ] ] denotes the sum / integral over all modes @xmath189 .
the above formal manipulations require a regularization such that @xmath190 \rightarrow 1 $ ] .
we carry out the appropriate renormalization by the subtraction of the free energy of the free field without boundary conditions , corresponding to replacing the partition function through the ratio @xmath191 , with @xmath192 the partition function of the free field . for the latter ,
we define the ` boundary condition ' @xmath193 in such a way as to produce the fraction @xmath194 of the free field energy in the volume @xmath134 via the expression ( [ f_lnd ] ) ; for the classic casimir problem , @xmath195\}^{r / l}$ ] . with @xmath196 we finally obtain the regularized free energy of the casimir problem@xcite @xmath197 . \label{fren}\ ] ] we illustrate the use of this formalism with a brief derivation of the classic casimir result , the attraction between two parallel metallic plates of size @xmath198 and a distance @xmath0 apart .
the dispersion relation for the free electromagnetic field is @xmath199 , such that the function @xmath200 takes the form @xmath201 .
the boundary condition requiring the fields to vanish at @xmath157 and @xmath0 can be cast into the form @xmath202 , producing the modes with a longitudinal @xmath174 vector @xmath203 .
the ` boundary condition ' describing the free field takes the form @xmath204 , such that @xmath205 , where @xmath206 . combining the results for @xmath207 and @xmath208
, we have to carry out the mode summation over the logarithm of @xmath209 =   \left[1-\exp\left(-2r \sqrt{\frac{\xi_s^2}{c^2}+q^2}\right)\right]^2,\ ] ] where taking the square accounts for the two polarization modes of the electromagnetic field ( the single polarization mode at @xmath210 has been properly taken into account here ) .
the force density @xmath32 between the plates is given through the derivative of the energy ( [ fren ] ) with respect to @xmath0 , @xmath211 where we have replaced the matsubara sum through an integral over complex frequencies @xmath74 .
the remaining integrals are easily evaluated by changing to polar coordinates @xmath212 , @xmath213 and we obtain the classic result due to casimir@xcite @xmath214      in the following , we apply the above formalism to the 2d charged bosons described by the action @xmath215 as given in eqs .
( [ plm ] ) and ( [ edr ] ) . following the scheme outlined in fig .  2
, we first have to integrate over the currents @xmath216 .
we split the boson positions @xmath217 into the mean field part @xmath218 and a fluctuating part @xmath219 , @xmath220 . in carrying out the integration over the boson positions @xmath221 we can ignore the fluctuating part @xmath219 in the gauge field @xmath11 as the latter is smooth on the scale of the amplitude @xmath219 , @xmath222   \exp\bigg\{-\frac{1}{\hbar^{\scriptscriptstyle b } } \int{\rm d}\tau\sum_\mu \big[\frac{m}{2}(\partial_\tau { \bf u}_\mu)^2   \nonumber \\ & & \qquad\qquad\qquad\qquad\qquad + i(\partial_\tau { \bf u}_\mu,1)\cdot   { \bf a}({\bf r}_\mu^0,\tau)\big]\bigg\ } \nonumber \\ & & \propto \exp\left[-\frac{1}{\hbar^{\scriptscriptstyle b } } \!\ ! \int \ ! { \rm d}^2 r{\rm d}\tau   \left(\frac{1}{2 } \pi({\bf r } ) a^2_{xy}({\bf r } ) +   i \rho({\bf r } ) a_\tau({\bf r } )
\right)\right ] , \nonumber\end{aligned}\ ] ] where @xmath223 denote the polarizability and the density of the 2d bosons , respectively .
second , we have to integrate over the true gauge field @xmath37 ( we have fixed the gauge to @xmath49 ) .
the resulting term @xmath224 in the action then renders the fake gauge field @xmath11 massive ( physically , this mass term expresses the finite range @xmath4 of the interaction between the charged bosons / vortices ) .
third , we introduce the free green function matrix @xmath225 for the gauge field , in fourier representation , @xmath226 , \nonumber \\   \noalign{\vskip 5 pt } { \rm with } \quad ( { \bf g}_{xy})_{\alpha \beta}^{-1 }   & = & k^2 \delta_{\alpha\beta } + ( c^2 - 1 ) k^2 p_{\alpha \beta}({\bf k } ) , \nonumber \\ g_\tau^{-1 } & = & k^2 , \nonumber\end{aligned}\ ] ] and the transverse projector @xmath227 . combining the results of the above three steps , we arrive at the desired effective action for the fake gauge field @xmath11 , @xmath228 & = & \int\!\!{\rm d}^2 r { \rm d}\tau \left [ { \bf a}_{xy } \frac{1}{2g^2}\left({\bf g}_{xy}^{-1 } +   \frac{1}{\lambda^2 } + g^2\pi \right ) { \bf a}_{xy }   \right .
\nonumber \\ & & \left .
+ a_\tau \frac{1}{2g^2}\left(g_\tau^{-1 } +   \frac{1}{\lambda^2}\right)a_\tau + i\rho a_\tau \right ] .
\label{thing}\end{aligned}\ ] ] the euler - lagrange equations of this action determine the ( imaginary time ) field equations for the transverse and longitudinal parts of @xmath11 ( we remind that @xmath229 ) , @xmath230 and @xmath231 with @xmath232 .
in the last equation of ( [ trans_eom ] ) we have made use of the gauge condition @xmath233 and have ignored the time dependence in the longitudinal field component @xmath234 , see below .
the longitudinal field @xmath234 is generated by the source @xmath235 and has no dynamics on its own .
placing two bosons at a distance @xmath0 , a simple integration gives the vortex - vortex interaction ( c.f . , eq .
( [ rep_eq ] ) ) @xmath236 here , we ignore the fluctuations of the vortex position in the calculation of @xmath234 as they merely produce a small renormalization of the prefactor in the repulsive potential .
similarly , we can neglect such corrections in the calculation of the transverse modes below and set the r.  h.  s.  of ( [ trans_eom ] ) equal to zero .
the dynamical transverse field @xmath237 generates the casimir force . here
, we have in mind a geometry as sketched in fig .  1 , with two parallel cubes of vortex matter of density @xmath17 , separated by a vortex free region of thickness @xmath0 .
we choose the @xmath26-axis to lie perpendicular to the ` plates ' , the @xmath238-axis is directed along the vortices . within the boson language we deal with two parallel planes of 2d bosons lying in the @xmath60-plane , the direction along the vortices now transforming to the imaginary time coordinate ( below , we consider the limit @xmath239 , implying vortex lines of infinite length , @xmath240 ) .
note that in the present 2d case we have to consider only one polarization mode for the transverse gauge field @xmath237 , a consequence of the gauge condition @xmath49 .
going over to fourier space we can cast ( [ trans_eom ] ) into the form ( we remind the reader that @xmath241 ; @xmath242 and @xmath174 are the wave vectors along and perpendicular to the ` plates ' ) @xmath243 { \bf a}_{xy } = 0 , \label{trans_eom_k}\ ] ] from which we obtain the green function @xmath244 and the function @xmath208 , @xmath245 here , we have defined the ` dielectric ' constant @xmath246 where we go over to vortex parameters in the last equation ( @xmath247 ) and make use of the expression ( [ elasticity ] ) for the dispersive vortex elasticity @xmath248 .
next , we have to formulate the boundary conditions in terms of the zeroes of the function @xmath249 .
the translational invariance in the @xmath179-direction allows for a plane - wave ansatz @xmath250 , while the sequence of ` dielectric ' and ` vacuum ' regions along the @xmath26-axis , requires to match the plane waves @xmath251 , with ( see ( [ trans_eom_k ] ) ) @xmath252 the gauge condition @xmath253 is satisfied through the ansatz @xmath254 with piecewise constant amplitudes @xmath255 and @xmath256 .
the six coefficients are determined by the boundary conditions at @xmath257 and @xmath258 , requiring the parallel electric field @xmath259 and the magnetic field @xmath260 to be continuous , see appendix b. two further conditions force the fields to vanish at @xmath261 , @xmath262 . requiring the determinant of the resulting @xmath263 matrix problem to vanish ,
we find @xmath264 & = & e^{-ik_l l } e^{-ik_0 r } e^{-ik_r l } \label{dvdw } \\ & & \left[1 - \frac {    ( 1 + ( \lambda\xi)^{-2})k_l-(\epsilon_{\scriptscriptstyle v}+     ( \lambda\xi)^{-2})k                  } {    ( 1 + ( \lambda\xi)^{-2})k_l+(\epsilon_{\scriptscriptstyle v}+      ( \lambda\xi)^{-2})k                 } \right .     \nonumber \\ & & \times   \left .
\frac {    ( 1 + ( \lambda\xi)^{-2})k_r-(\epsilon_{\scriptscriptstyle v}+    ( \lambda\xi)^{-2})k                 } {    ( 1 + ( \lambda\xi)^{-2})k_r+(\epsilon_{\scriptscriptstyle v}+    ( \lambda\xi)^{-2})k                  } e^{2 i k r}\right ] .
\nonumber\end{aligned}\ ] ] for the free field , we find @xmath265 and combining the above results for @xmath208 , @xmath177 , and @xmath193 , eqs .
( [ fvdw ] ) , ( [ dvdw ] ) , and ( [ dvdw_0 ] ) , we can construct the function @xmath266 $ ] and obtain the formal expression for the free energy ( [ fren ] ) of the vortex casimir problem .
we note that the cut introduced through the definitions of @xmath267 and @xmath268 in eq .
( [ ki ] ) disappears from the boundary condition @xmath177 and hence does not contribute to the loop integral ( [ c_loop ] ) . as a result ,
the free energy expression ( [ fren ] ) for the casimir free energy remains valid .      in the present two dimensional situation the expression for the casimir energy eq .
( [ fren ] ) per length reads @xmath269,\ ] ] and produces the casimir force per length @xmath270 @xmath271^{-1 }                     } .
\nonumber\end{aligned}\ ] ] we introduce the new variables @xmath272 and arrive at the simplified expression @xmath273^{-1 }                     } .
\nonumber\end{aligned}\ ] ] eq .
( [ force_gen ] ) gives the two dimensional version for massive photons of the casimir force between dielectric media , first derived by lifshitz@xcite for the conventional 3d situation .
in the following we will be interested in the situation where the ` dielectric ' media are dilute ( dilute vortex matter at small magnetic induction @xmath274 ) . in this case , we expand ( [ force_gen ] ) in the small correction @xmath275 giving the deviation of the dielectric constant from its vacuum value .
furthermore , if the integral is dominated by large frequencies @xmath74 we can ignore the mass term everywhere . under these conditions , we approximate @xmath276 and @xmath277 , such that @xmath278 , while @xmath279 , and obtain the simplified force expression ( the weak logarithmic disperion in @xmath280 is approximated through an appropriate constant ) @xmath281 in ( [ force_gen ] ) , ( [ force_dil ] ) we have to account for three relevant frequency- or length scales . the first is introduced by the discrete nature of our superconductor : the layered structure limits the @xmath74-integration through the frequency cutoff @xmath282 ( this corresponds to the medium becoming transparent at high frequencies @xmath283 ) .
the second frequency scale is introduced through the dispersion in the dielectric constant : at low frequencies @xmath284 , the dielectric constant takes the form @xmath285 and the dielectric response changes over to a mass renormalization .
the third frequency scale @xmath286 appears through the exponential @xmath287 in the integrand of ( [ force_gen ] ) ; as @xmath74 goes beyond @xmath288 the exponent cuts off the integrand in ( [ force_gen ] ) .
thus , depending on the frequency @xmath74 , the media are either transparent , give a dielectric response , or produce a mass renormalization . and depending on the position of @xmath288 with respect to these response regimes , we will find a different behavior of the casimir force .
below , we will analyze the various regimes and use the pairwise summation technique to show that the results are in agreement with those obtained for the van der waals interaction potential in sec .
iv above .      at intermediate distances
@xmath127 , we have @xmath289 such that the cutoff on the @xmath74 integral is given by @xmath290 . as the integral
is dominated by large values of @xmath74 we can ignore the mass terms @xmath291 everywhere .
the exponential @xmath292 restricts the @xmath293-integral to values @xmath294 , admitting large values of @xmath293 such that we can drop the corrections to @xmath295 in ( [ force_dil ] ) , e.g. , @xmath296 . transforming @xmath293 to the new variable @xmath297 , we find the casimir force on intermediate scales @xmath298 using the result ( [ f_of_r ] ) of the pairwise summation and inserting the expression ( [ epsilon_v ] ) for the dielectric constant , we obtain the van der waals interaction describing the decoupled limit , @xmath299 in agreement with the result ( [ v4_vdw ] ) ( note that the density @xmath17 cancels in the final result for the van der waals interaction ) .      at large distances ,
the @xmath74-integral is cutoff on the scale @xmath288 , i.e. , the casimir force is limited by the distance @xmath0 through the exponential @xmath292 rather than by the transparency of the material .
we still can ignore the mass terms in ( [ force_gen ] ) as we assume that @xmath301 , providing us with a large regime for the @xmath74-integration where @xmath302 .
transforming the energy variable @xmath74 to @xmath303 , we obtain the following expression for the casimir force in the low density limit , @xmath304 following the usual scheme , this result produces the retarded van der waals interaction @xmath305 in agreement with eq .
( [ v5_vdw ] ) .
the situation at very large distances @xmath307 is most conveniently analyzed in the original formulation ( [ force_gen_0 ] ) .
the mode summation is limited by the exponential @xmath308 , implying the following restrictions on the integration variables @xmath74 and @xmath242 : @xmath309 and @xmath310 .
the integral @xmath311 $ ] then is given by the area @xmath312 times the @xmath313 limit of the integrand .
in the limit @xmath314 the mass terms are relevant and the dispersion in the dielectric constant ( [ epsilon_v ] ) produces the additional mass renormalization @xmath315 the casimir force then decays exponentially following the law @xmath316 for small densities , we can expand @xmath317 and obtain for the screened van der waals interaction in the regime @xmath306 the expression ( note that the pairwise summation formula ( [ f_of_r ] ) has to be modified for the exponential factor in the interaction @xmath318 ) @xmath319
we have derived the casimir force between two bodies of vortex matter and have inferred from the result the strength of the van der waals interaction between vortex lines via the method of pairwise summation .
while the calculation of the casimir force is carried out in the boson formulation , the van der waals attraction is determined in the vortex picture .
the agreement between the results is once more an illustration of the equivalence of the two formalisms @xcite .
the physical implications of these results lead to interesting modifications of the @xmath320-@xmath7 phase diagram of layered / anisotropic type ii superconductors@xcite : the attraction between the lines produces a generic vortex solid of density @xmath321 at low temperatures .
an interesting problem appearing in this context is the accurate determination of the entropic repulsion between the vortex lines , the latter giving an important contribution to the gibbs free energy .
various approaches to this problem have been discussed by blatter and geshkenbein@xcite , by volmer , mukherji , and nattermann@xcite , and by volmer and schwartz@xcite .
the transition from the meissner state to the van der waals vortex solid takes place via a sharp first - order transition at the lower critical field @xmath322 . the concurrent phase separation
then can be observed in a decoration experiment , though very clean samples are required .
furthermore , the phase separation leads to the concentration of the vortex lines into dense regimes separated by vortex free regions .
the interaction between the vortex domains then is given by the casimir force calculated here , thus providing a natural realization of the physics described in this paper .
a further remark is in place concerning the physical character of the van der waals interaction discussed here .
the usual van der waals interaction between neutral atoms arises from dipolar fluctuations in the charge distribution of the atoms , resulting in a force which is mediated through the scalar ( longitudinal ) potential , at least within the non - retarded regime at small distances . comparing to the van der waals force between the 2d charged bosons / vortex lines discussed here , we note that the elementary objects do not have an internal structure producing a fluctuating dipole .
the van der waals interaction then arises from _ current _ fluctuations and thus involves the ( transverse ) vector field . in principle
, such a ` transverse ' van der waals force is also present in an electron gas : while the longitudinal repulsive interaction is screened on the thomas - fermi length , the transverse attractive interaction of the van der waals type due to fluctuating currents survives at long distances ( and in principle induces a superconducting instability ) .
however , the transverse van der waals force in an electronic system involves the small parameter @xmath323 ( @xmath324 is the fermi velocity ) , rendering the effect small . in the 2d boson system
discussed here , the light velocity is given by the anisotropy of the superconductor , @xmath241 , and thus is a tunable parameter .
we thank vadim geshkenbein , a.  van otterlo , and christoph bruder for fruitful discussions .
in this section we sketch , by way of example , the calculation of @xmath325 .
we start from eq .
( [ vab ] ) and set @xmath326 ( here , we ignore the exponential cutoff function due to the vortex core ; the cases @xmath327 and @xmath328 , @xmath329 follow trivially ) @xmath330 performing a contour - integral over @xmath331 , we obtain @xmath332 with @xmath333 and @xmath334 . carrying out the integral@xcite over @xmath335 in @xmath336 we find @xmath337 the second integral @xmath338 can be calculated in the same manner , @xmath339 here , @xmath3 denotes the @xmath157-th order modified bessel functions of the second kind and @xmath340 .
collecting results and carrying out the derivatives , we obtain the desired result . here , we quote the final expressions for the case @xmath341 : @xmath342     \right .
\nonumber   \\ & &     \left .
\qquad\qquad\qquad   -\frac{a}{\lambda_c r } k_1\left(\frac{ar}{\lambda_c}\right )    \right\ } , \label{app_a_x } \\
v_{yy}^{\rm int}(r , k_z ) & = &    \frac{1}{2 \pi a^2 }    \left \{\frac{a^2}{\lambda_c^2}\left[k_0\left(\frac{ar}{\lambda_c}\right )            + \frac{\lambda_c}{ar}k_1\left(\frac{ar}{\lambda_c}\right )                               \right ]    \right .
\nonumber   \\ & &    \left .
\qquad\qquad\qquad -\frac{a}{\lambda r } k_1\left(\frac{ar}{\lambda}\right )    \right\}. \label{app_a_y}\end{aligned}\ ] ] the third component @xmath343 vanishes for @xmath344 . the important terms in ( [ app_a_x ] ) and ( [ app_a_y ] ) are those involving the large screening length @xmath345 in the argument of the bessel functions .
the first term @xmath346 in ( [ app_a_y ] ) has been missed in ref .  .
we sketch the main features of the 2d real time massive electrodynamics underlying the present work ( we denote the mass by @xmath347 and restrict ourselves to the isotropic case ) .
we introduce the metric @xmath348 and define the gauge field @xmath349 .
in addition , we define the derivatives @xmath350 , where @xmath351 .
the field tensor @xmath352 is given by @xmath353 and its dual takes the form @xmath354 ( @xmath355 is the antisymmetric tensor ) .
the lagrangian of the massive field coupled to an external current @xmath356 is given by @xmath357 with @xmath358 the coupling constant .
this lagrangian produces an imaginary time action @xmath359 =   \int \!\ ! { \rm d } \tau { \rm d}^2 r \left [   \frac{1}{2g^2}(\nabla\times{\bf a})^2 +   \frac{1}{2 g^2 \lambda^2}{\bf a}^2 +    { \bf a } { \bf j } \right].\ ] ] the transformation from real time to imaginary time is obtained via the formal rules @xcite @xmath360 , @xmath361 , @xmath362 .
the magnetic field remains unchanged , while the electric field transforms to @xmath363_{\perp } = i { \bf e}$ ] .
the free field real time lagrangian @xmath364 then goes over into the simple form @xmath365 within the imaginary time formalism .
the functional derivative of eq .
( [ l_aj ] ) provides us with the inhomogeneous maxwell equations ( the homogeneous one , @xmath366 , follows from the antisymmetric structure of @xmath352 ) , @xmath367 making use of gauss and stokes theorems for a loop encircling the boundary between two media , we obtain the boundary conditions @xmath368 continuous and @xmath369 continuous .
the quantization of the theory can be developed via the gupta - bleuler formalism within the lorentz gauge @xmath370 . in the case of a massive theory
this method leads to two polarization modes , while in a massless theory the gauge invariance reduces the number of polarization modes by one . within the present theory ,
the coupling of @xmath11 to the real @xmath37 field in the form @xmath371 produces a gauge invariant mass term @xmath372 .
as the gauge invariance is conserved in our theory we end up with a single polarization mode .